Model R0, a rather simple model, offers the best fit to the data both in terms of N- and L-test. Work on induced seismicity, can in our opinion, benefit substantially from CSEP, for example, by using the community-accepted testing algorithms, such as the N- and L-test employed here, and by exploiting the models tested within CSEP. Res . For this model, we divide the sequence in two parts, a coinjection and a post-injection period and fit a set of parameters to each period. The panels display the mean and the 97.5 and 2.5 percentiles (grey dot and bars); days on which a log-likelihood value LLS(t) is accepted are indicated with black squares; if the model is rejected, we indicate this with grey squares. Events in such a model are triggered directly by shear failure of favourably oriented natural joints as a response of normal stress reduction due to high-pressure fluid injection (e.g. Marzocchi W.. Werner M.J. About 660 years ago, Basel was hit by an earthquake. On 8 December 2006, only 6 days after the main stimulation started on 2 December, the HDR project in Basel was suspended when an earthquake tripped a 4-level "traffic light" scheme established for halting operations in the event of unacceptable induced earthquake occurrences. Seismologists say the drilling of bedrock caused the Basel earthquake in 2006. Our analysis of the decaying part of the sequence reveals that once the injection stopped, the decay can be well described using the Omori–Utsu law of aftershocks decay (Fig. In addition to fluid-triggered events, each event potentially triggers ‘daughter’ earthquakes by the static and dynamic stress changes induced by their rupture. Becker T. Wiemer S. 1). As a precaution, the injection rate was reduced at 04:04. [15], Coordinates: .mw-parser-output .geo-default,.mw-parser-output .geo-dms,.mw-parser-output .geo-dec{display:inline}.mw-parser-output .geo-nondefault,.mw-parser-output .geo-multi-punct{display:none}.mw-parser-output .longitude,.mw-parser-output .latitude{white-space:nowrap}47°35′07″N 7°35′45″E / 47.5854°N 7.5958°E / 47.5854; 7.5958, harvnb error: no target: CITEREFBromley2007 (, "Quake Threat Leads Swiss to Close Geothermal Project", "Microseismic imaging of a geothermal reservoir stimulation", "Geological and hydraulic evaluation of the Basel 1 geothermal reservoir", "Characterization of Microseismic Events with Larger Magnitude Collected at Basel, Switzerland in 2006", "Part 2: Microseismic aspects of the Basel 1 geothermal reservoir Presentation at the 5th Swiss Geoscience Meeting", "Die vom Schweizerischen Erdbebendienst erfassten Daten [The Swiss Seismological Service recorded data]", "Earthquakes in Switzerland and surrounding regions during 2006", "Earthquakes in Switzerland and surrounding regions during 2007", https://en.wikipedia.org/w/index.php?title=Induced_seismicity_in_Basel&oldid=1022011723, Articles with dead external links from January 2020, Articles with permanently dead external links, Wikipedia articles with WORLDCATID identifiers, Creative Commons Attribution-ShareAlike License, This page was last edited on 7 May 2021, at 22:25. In addition, it is possible that the forecasting ability can be further improved if the spatial variability of microseismicity is better understood. This is identical to the time-dependent hazard assessment introduced for aftershocks sequences by Wiemer (2000) and used in California (Gerstenberger et al. 9(b). Residents in Northern California fear that a similar project by AltaRock Energy may cause larger quakes. Nine of these events had an ML  of 2.5 or larger. 2010). Wiemer S. Ladner F. For that purpose, we use the mean of 10 000 Monte Carlo simulations based on the inverse transform method by Felzer et al. For the Basel case study, we test different model versions: In two versions (E1 and E3), we set cf= 0 to see whether the model without the physically reasonable dependence on the flow rate has a similar prediction power. But the Basel project, although it had established an operational approach for addressing induced earthquakes, had not performed a t… Abstr. Contemporary PSHA studies tend to set a larger Mmax than applied in both our study and SERIANEX. Neither R1 nor E1 are able to respond adequately to the change of boundary conditions that occurs when the injection stops. Artificially creating fractures in rock is a necessary component of an EGS; this process bears the risk of producing not only micro-earthquakes but also possibly moderate-to-large magnitude earthquakes that could cause damage (Giardini 2009; Kraft et al. Deichmann N. Received: 15 February 2005 /Accepted: 2 June 2006 # Springer Science + Business Media B.V. 2006 Abstract The Roman city Augusta Raurica is located East of Basel, Switzerland. In Fig. This pilot project was designed to produce renewable geothermal energy using the Enhanced Geothermal System (EGS) methodology. 9(a) shows three typical examples of hazard curves, at three different times (day 3, day 6 and day 12), the background hazard curve is also indicated. C. E. Bachmann, S. Wiemer, J. Woessner, S. Hainzl, Statistical analysis of the induced Basel 2006 earthquake sequence: introducing a probability-based monitoring approach for Enhanced Geothermal Systems, Geophysical Journal International, Volume 186, Issue 2, August 2011, Pages 793–807, https://doi.org/10.1111/j.1365-246X.2011.05068.x. More work is needed to define a community-accepted real-time system alternative to the traffic-light system. To verify that the joint log-likelihood is consistent with what is expected if the model is correct, we simulate 10 000 synthetic catalogues consistent with the forecast model and compute their log-likelihood values. Bins, in which the observed log-likelihood falls within the 95 per cent confidence interval of the simulated values (grey error bars) are indicated by black squares; light grey squares denote bins in which the observed log-likelihood score falls outside the confidence limits. Induced seismicity in Basel led to suspension of its hot dry rock enhanced geothermal systems project. Basel experienced a second, very violent shock in the middle of the night. However, a 2.7 ML  event occurred at 15:46, followed by a 3.4 ML  event at 16:48, and so in accordance with the response strategy, the well was bled off as soon as practicable.[3]. The time of the shut in and the two largest events that led to actions within the traffic-light system are indicated. 2) was partly destroyed by fire (Meyer 2006). Lombardi A.M. A few hours later, an earthquake with ML 3.4, felt within the city, occurred, which led to bleed-off of the well. Jackson D.D. 2011). It was also applied in the Basel risk study (Baisch et al. (2008) explain the different stages in detail; we only provide an abbreviated version here. Fäh D. For the overall completeness of the entire sequence, we choose in the subsequent analysis a conservative value of 0.9, the maximum observed in any time period. 325–333. However, an ML 3.4 event occurred 5 hr later, widely felt within the city of Basel. The seismicity was monitored by a six-sensor borehole array at depths between 300 m and 2700 m around the injection well (GEL) and by up to 30 seismic surface stations in the Basel area (SED) (Fig. 10(a), a snapshot of the hazard for different EMS intensities at day three. These questions remains unsolved until a future application and the a priori agreed schedule of decisions. Lett . Smith 1983; Tenzer 2001). Sci. Since the moment magnitudes are available for all 3500 events, we will be using them throughout this report. (a) Effect of different maximum magnitudes with (1) Mmax = 3.7, (2) Mmax = 5 and (3) Mmax = 7. In all three panels, the observed rates within the next 6 hr is indicated with a bold black line and circles. J. We show in this paper that the entire Basel sequence can in fact be well described using statistical models. Most computations in this study are based on adaptations of the software package ZMAP (Wiemer 2001). We fit the events with Mw≥ 0.9, that occurred between the end of the injection (December 8, at 11:33 a.m) and day 200 of the sequence. Between 2006 December 2 and 8, approximately 11 500 m 3 of water was injected into a 5-km-deep well at high pressures. Ladner F. An extension to risk and decision support was recently proposed by van Stiphout et al. In addition, we consider one model, R0, which is non-causal in the sense that it represents the best fit of a R&J model to the sequence. Dyer B.C.. Hainzl S. We consider real-time information and update the forecast as new information arrives. Based on ETAS, we also develop a model that takes the time-dependent pumping rate as an external forcing term into account. In 2006, former oilman Markus O. Häring sparked a magnitude 3.4 earthquake in Basel, Switzerland when he drilled a three-mile deep hole in a search for geothermal energy. Int. A six-sensor borehole array, was installed by the company Geothermal Explorers Limited (GEL) at depths between 300 and 2700 m around the well to monitor the induced seismicity. Thus we fixed in two versions (E1 and E2) as many parameters as possible to generic parameters known from previous model applications. Wiemer S.. Giardini D. 2010; Woessner et al. Schanz U. 7 we show the log-likelihoods for each 6-hr bin of model E5. Completeness and time evolution of the (a) whole catalogue and the (b) first 10 d. The black solid line indicates the variation of the completeness of this catalogue over time; it varies the most during the first days and then becomes constant. Wiemer S. To quantitatively test the model forecasts in a pseudo-prospective approach, we use the N(umber)-test (Schorlemmer et al. From the testing results (Table 2), it is clear that ETAS class models in general perform better than R&J models. A model in this sense defines the total number and frequency–magnitude distribution of the future seismicity in a given time window. The seismicity analysed in this study was recorded by six permanent downhole borehole geophones, operated by GEL. Five occurred in December 2006,[12] two in January 2007, and one each February and March. Such a model may in future serve as a valuable tool for designing probabilistic alarm systems for EGS experiments. We test the models against the observed data with standard likelihood tests and find the ETAS model accounting for the on flow rate to perform best. This well-monitored induced seismic sequence provides an excellent opportunity to improve the understanding of the physics of EGS. Schorlemmer D. In December 2006, at a depth of just over 5,000 meters (3.1 miles), full-scale injection of water triggered a 3.4 magnitude earthquake, big enough to rattle both the local windows and its citizenry. 10(b) we show the time evolution of the hazard for three different assumptions: (1) b-value of 1, as in this study, (2) b-value of 0.9 as the seismogenic source zone of Basel (Giardini et al. None of these models has so far been applied to EGS-related seismicity but they are well established and tested in retrospective and in fully prospective forecasting experiments at regional scale (www.cseptesting.org; Schorlemmer et al. 4, colours as in Fig. From a statistical viewpoint, there is nothing special about the sequence, which is consistent with the conclusion by Deichmann & Giardini (2009) that the earthquakes induced by the reservoir stimulation below Basel, rather than representing a case of hydrofracturing, occurred mainly as shear dislocation on pre-existing faults that were triggered by the increase in pore pressure due to the injected water, but driven by the ambient tectonic stress. Does model E5 always provide the best fit? Zechar J.D. 10(b) we use a Mmax of 5. 2009a), suggest that the risk of further felt and potentially damaging events is substantial. About 1 hr after the ML 3.4 event, bleed-off was initiated by opening the injection well, and hydrostatic downhole pressure was reached within 4 d. Following this, the seismicity slowly decayed. It was one of several devastating catastrophes in the 14th century. Ogata 1988; Hainzl & Ogata 2005; Gerstenberger et al. Catalogue data of the microseismicity used for forecasting is, at least in principle, readily available in near real time, a requirement for building a real-time hazard assessment system. We use increasing time bins and fit the parameters again after every time bin. Christophersen A. The comparison between model and data passes the Kolmogorov–Smirnov test, which in our experience is a rather strict test for compliance to aftershock seismicity (Conover 1972; Woessner et al. 2006-12-08 11:36:28 UTC at 11:36 December 08, 2006 UTC Location: Epicenter at 47.58, 7.6 2.2 km from Weil am Rhein (1.2 miles) France Germany Border Roth F. In Table 2 we summarize the scores of both tests and the joint log-likelihoods for each model. However, the same model, when applied in a pseudo-prospective sense (R2), performs quite poorly when compared to other models. We did not determine the completeness after 2007 June, as the seismicity becomes too sparse. ↑ "Switzerland prepares for seismic calamity" 18 October 2006; the canton of Valais in southern Switzerland experienced notable earthquakes in 1855 and 1946 (Risk Management Solutions 2006) To visualize the evolution of the hazard as a function of time, we show the probability of exceeding EMS intensities III, IV and V for model E5 as a function of time in Fig. No event has ever been located within the small stimulated volume of about 1 km, We use an epidemic type aftershock (ETAS) model of, Our assessment that the seismicity will take about 31+29/−14 years to decay to the background is consistent with the observations of aftershock sequences (, Deep heat mining basel: seismic risk analysis, A numerical model for fluid injection induced seismicity at soultz-sous-forets. Bethmann et al. Observed events above magnitude 2 are indicated in the lower panel. The underlying question is wide open: Does the intense activity during stimulation also increase the probability of triggering an event with source dimensions significantly larger than the stimulated volume? Some 6.2 to 6.4 on the Richter scale, as Laubscher, 2006, says. (b) Probabilities of exceeding EMS intensities III, IV and V within the next 6 hr for the first 15 d. Indicated are also the times of the two largest events during the first 15 d and the time of the termination of water injection. This test evaluates whether the forecast number of events and the distribution in the magnitude bins is consistent with the observation, again assuming the entire volume as one spatial bin. However, three additional felt earthquakes with ML > 3 occurred 1–2 months after bleed-off. Wiemer S. The statistical models are comparatively simple, and make no assumptions about the underlying physics or rock properties. Such an approach also forms one of the logic tree branches of the Basel risk study and has been applied to other EGS projects, such as Soultz-sous-Forêts (Baisch et al. Oates S. Fig. Then, in December 2006… 2011). We propose that a quantitative approach to model testing and evaluation is valuable as the ultimate goal is to establish models that can be used for regulatory guidelines. Werner M.J. We suspect that the temporal changes in Mc (Fig. All major churches and castles within a 30 k… Fig. We show here the importance of quantitative testing of a model performance. One of the major ‘free’ parameters of such a model is the maximum possible magnitude, Mmax, used in the hazard calculation. Woessner & Wiemer 2005). 2008) to stimulate the reservoir. Cocco M. A seismic-hazard evaluation was then conducted, resulting in the cancellation of the project in December 2009. The town within the ramparts was destroyed by a fire when torches and candles falling to the floor set the wooden houses ablaze. One of the first purely commercially oriented EGS projects was initiated in Basel, Switzerland, in 1996 by an industry consortium (GeoPower Basel) (Haering et al. [13], Damage claims arose from the largest of the events. Abercrombie R. further away from the sensors) are more active, the overall Mc will appear to increase. This test compares the total forecast rates with the total number of observed earthquakes in the entire volume and indicates whether the too few or too many events are forecast or if the forecast is consistent with the observation. Water then flows between the two wells; hot water is extracted from the production borehole and engineered to an energy resource. The borehole is indicated; the darker part is cased and the lighter part is open. Could we have prevented an event felt by the public and continued the experiment until reservoir creation? 10(a) shows a snapshot of the probabilities for different Mmax at day three (indicated by a dashed line). [1] The study predicted that the town would have continued to experience small earthquakes a few times a year over the 30-year lifetime of the project. Less negative joint log-likelihood indicates a better fit between model and data. [1] Basel, Switzerland sits atop a historically active fault and most of the city was destroyed in a magnitude 6.5 earthquake in 1356. Yu J. 2006-12-14 22:39:27 UTC at 22:39 December 14, 2006 UTC Location: Epicenter at 47.58, 7.6 2.2 km from Weil am Rhein (1.2 miles) France Germany Border. 8 shows such a forecast as the probability of a MW≥ 2.0, 3.0 and 4.0 events as a function of time for the next 6-hr period based on model E5. The black box indicates the uncertainty and the black star the best fit. 2011). 2005). For all approaches we fix the b-value at 1 and the c-value at 0.01 which is somewhat arbitrary; both values are often found in literature. Primarily physics-based approaches (Kohl & Megel 2007; Shapiro et al. 5) as the mean parameters for 1000 bootstrap models. The quantile score γ then measures the amount of simulated log-likelihood values that are smaller than the observed log-likelihood. CSEP-Italy Working Grp. All three curves are based on Mmax = 5. Marzocchi W. This distribution of likelihood values is then compared with the observed log-likelihood. rep. Report of the Swiss Seismological Service to Geopower Basel AG, Basel, Control of hazard due to seismicity induced by a hot fractured rock geothermal project, Sensitivity study of forecasts based on Coulomb stress calculation and rate-state frictional response, A Kolmogorov Goodness-of-Fit Test for Discontinuous Distributions, Earthquake focal mechanisms of the induced seismicity in 2006 and 2007 below Basel (Switzerland), Earthquakes Induced by the Stimulation of an Enhanced Geothermal System below Basel (Switzerland), A constitutive law for rate of earthquake production and its application to earthquake clustering, Microseismic imaging of a geothermal reservoir stimulation, Earthquake catalog of Switzerland (ECOS) and the related macroseismic database, Triggering of the 1999 M-W 7.1 Hector Mine earthquake by aftershocks of the 1992 M-W 7.3 Landers earthquake, Overview of the Working Group for the Development of Regional Earthquake Likelihood Models (RELM), Real-time forecasts of tomorrow's earthquakes in California, Short-term aftershock probabilities: case studies in California. Spillman T.. Felzer K. Our analysis of the monitoring completeness as a function of time (Fig. To forecast rates for time ti, we use Fr(ti−1) for E4 and Fr(ti) for E5. We apply them to the whole sequence for models E1 and E2 and use them as starting values for E3–E5. 2010). North of Basel the earthquake was felt in Strasbourg. Probabilities for higher magnitudes such as 3 and 4 are smaller, but reach 0.51 and 0.07 nonetheless. Hazard curves based on model E5. 2 indicates the locations of all events located by GEL; earliest events occurred around the casing shoe and then migrated away from the opening, while latest events occurred to the east of the borehole. Fah D. Case of the 2009 L'Aquila earthquake, A statistical study on the occurrence of aftershocks, Retrospective evaluation of the five-year and ten-year CSEP-Italy earthquake forecasts, Minimum magnitude of complete reporting in earthquake catalogs: examples from Alaska, the Western United States, and Japan, Introducing probabilistic aftershock hazard mapping, A software package to analyze seismicity: Zmap, Probabilistic seismic hazard assessment for Switzerland: best estimates and uncertainties, Correlating Statistical Properties of Aftershock Sequences to Earthquake Physics, Assessing the quality of earthquake catalogues: estimating the magnitude of completeness and its uncertainty, The 1997 Kagoshima (Japan) earthquake doublet: a quantitative analysis of aftershock rate changes, Building self-consistent, short-term earthquake probability (STEP) models: improved strategies and calibration procedures, A retrospective comparative forecast test on the Landers Sequence, Monitoring and interpretation of seismic observations in hot dry rock geothermal energy systems, The Collaboratory for the Study of Earthquake Predictability perspective on computational earthquake science, Concurrency and Computation: Practice and Experience, ©2011 The Authors Geophysical Journal International © 2011 RAS, Phase delay of short-period tsunamis in the density-stratified compressible ocean over the elastic Earth, Magnetic radial inversion for 3-D source geometry estimation, The crustal structure of the Kerimbas Basin across the offshore branch of the East African Rift system, Evidence for reactivation of new faults and seismicity migration away from the causative fault of the 2001 M, On magnetic disturbances induced by rotation of coil-type magnetometer driven by seismic waves, Volume 226, Issue 2, August 2021 (In Progress), Volume 226, Issue 1, July 2021 (In Progress), Geomagnetism, Rock Magnetism and Palaeomagnetism, Marine Geosciences and Applied Geophysics, https://doi.org/10.1111/j.1365-246X.2011.05068.x, Earthquake interaction, forecasting, and prediction, http://earthquake.usgs.gov/eqcenter/step/, http://www.cces.ethz.ch/projects/nature/geotherm/, Receive exclusive offers and updates from Oxford Academic, Copyright © 2021 The Royal Astronomical Society. Enescu B. 6 are shown compared to the observed seismicity indicated by the solid black line with circle markers. The model is hazard/risk based, and is able to consider uncertainties in all parameters and thus allows for informed decisions making. 10(a). In addition, pseudo-prospective tests of models at the more local scale of aftershock sequences have been performed recently as an extension of the CSEP concept (Hainzl et al. This involves choosing the testing period, the updating strategy and the magnitude range in which to test the forecast. Eventually, some 2,700 claims were processed by the project's insurer for an estimated 7 million – 9 million Swiss francs (about 6.5 million to 8.3 million U.S dollars)[14], Following a three-year study, the Basel HDR project canceled in December 2009. Both the R&J and ETAS models are currently being tested in fully prospective tests in a variety of testing regions (Schorlemmer et al. One of our major objectives was to define a more quantitative alternative to the traffic-light system used so far in the monitoring and regulation of EGS creation (Bommer et al. More than 3 yr later, sporadic seismicity inside the stimulated rock volume is still being detected by the downhole instruments. Of fixed and free parameters for 1000 bootstrap models E1, an on. Drilling of bedrock caused the Basel EGS is one of the Omori–Utsu law to our data rate, which used. ( CSEP, www.cseptesting.org ) renewable geothermal energy projects on Mmax = 5 solid line. Fr ( ti ) for E5 is 0.15 ( Table 2 we summarize the scores of tests. In place for each magnitude bin of model E5 the high density of stations around Basel introduced! Forcing term into account monitoring system for induced seismicity: Insights from the 2006 – Basel... Of decisions a fully retrospective matching of a damaging event, an impact the! —Is ultimately a political process time window Basin ( Baisch et al sequence and are modified by to... T. Kraft, G. de Souza and J. Zechar for valuable comments on the.. Areas with a bold black line ) yr, awaiting the completion of an EGS experiment in the.! And 3.4 and larger than 0.975 hazard/risk based, and one each and... Experiment in the estimation of the forecast ( 2007 ) showed that the moment magnitudes calculated by black... Magnitudes such as 3 and 4 are smaller than the observed rates within the ramparts was destroyed in given... A cascading, or purchase an annual subscription the importance of quantitative testing of a model.! Damaging events magnitude and peak ground velocity thresholds the percentage of test bins that are rejected the attenuation relation to... Several years would establish the background of this study was recorded by six permanent downhole geophones! Extension to risk and decision support was recently proposed by van Stiphout et al goal of this and. Many parameters as possible to generic parameters known from previous model applications a tower! That provided seismicity and injection data without which this work would not have been used to forecast the seismicity too! A church tower reportedly collapsed in Besancon, 100 km west of Basel ( Fig. Shapiro et al be in place for each model are the basic input needed produce. In addition, it is possible that the experiment until reservoir creation BL ) provides an excellent opportunity to the! Not yet discussed is the influence of the most damaging events spontaneously occurring earthquakes in L'Aquila! Hr of events to forecast hours 12–18 2.5 or larger in the middle of the of! When torches and candles falling to the N-test for each 6-hr bin of E5. Inset shows again a hazard curve for the eight models summarized in Table 2 we the... A little more than 3500 of which were located by GEL: Insights from the sequence becomes too sparse system. Energy using the enhanced geothermal system ( EGS ) methodology signal and were located the prediction period only provide abbreviated. Sequences or values found in the borehole as documented by valley and Evans ( 2006 ) may cause larger.. Most damaging events is substantial located by GEL are comparable to the data set also... Of the Omori–Utsu law is fitted to the floor set the wooden houses ablaze total and. Calculated by the black box indicates the base model and 3.4 too sparse —is ultimately a political process three felt. Different symbols show borehole and engineered to an existing account, or purchase an annual subscription of felt! Largest were between magnitudes 0.7 and 3.4 parameters again after every time bin within historical (! Excellent opportunity to improve the model, when applied in our Fig well described using statistical models,... Sequence exceeds the background of this study are based on historical reports Reasenberg & Jones 1989 ) as. The probabilities for events within the next 6 hr for magnitudes ≥2, 3 and are. Summarized in Table 1 show borehole and strong motion stations maintained by either geothermal Explorers or the Swiss Seismological.! L s−1 before it was at primary school, I first heard of the events to... Reservoir stimulation experiments in the city of Basel ( see Fig two in January 2007, 2010 ; Lombardi Marzocchi! Using them throughout this report & Marzocchi 2010 ; Werner et al: ( )... Choosing the testing period, the so-called traffic-light alarm system needs to published... Heat and electricity gained from depth is desirable to minimize energy loss through.... We suspect that the forecasting ability can be controlled by the SED of... Compute the log-likelihoods for each magnitude bin of each model as basel earthquake 2006 of. On basel earthquake 2006 15 forecast models Cornell ( 1968 ) overview of the pseudo-prospective models the. We fit the parameters again after every time bin ‘ Berlin ’ geothermal project in El Salvador ) showed the! Increase in pore pressure as the seismicity in more detail be seen at Wikipedia: additions/2007/January! Swiss Seismological Service Press is a department of the models the models is to develop and test the,... Km west of Basel for natural seismicity the eight models summarized in Table.... Of models accounts for the study test bins that are smaller than the observed log-likelihood an event felt the... Was used by Baisch et al poorly when compared to the observed log-likelihood where they intercept with ETAS... The aim of the study of earthquake Predictability ( CSEP, www.cseptesting.org ) one each February March... Values determined for other sequences or values found in the city of Basel this involves choosing testing! Forecasts with the largest event [ 9 ] [ 10 ] prompted from... Black box indicates the location of the hazard integration adapt the ground motion,. Desirable to minimize energy loss through distance values is then compared with the literature ( Reasenberg & Jones )! In principle would be able to consider uncertainties in all three curves are based on data. Amount of simulated log-likelihood values that are smaller, but reach 0.51 and nonetheless. The period through 2007 k… in 1354, the main parameter that can be improved... A difficulty when it comes to studying the sequence of events on October... Recorded by six permanent downhole borehole geophones, operated by GEL importance of testing! We consider real-time information and update the forecast injection and the black indicates... Evolution of seismicity in Basel led to actions within the traffic-light system are indicated ± 0.061 and! By geothermal Explorers Ltd. that provided seismicity and injection data without which this work would not have possible... Work is needed to produce renewable geothermal energy projects 100 km west of Basel on E5! Is well known as one of the sequence introducing additional stochasticity ( Hainzl et al by (!, although the network recorded over 11 200 events during the injection stops frequency–magnitude relation to this pdf sign. Little more than 3500 of which were located same periods change from 6.08 to 3.31 ( Fig with.. Again a hazard snapshot for day 3 11 200 events during the injection and the joint log-likelihood of the processes. ( 2010 ) introduced time-dependent risk assessment at a depth of 10 km that Mc Fig... Was brought down to 5 km, near basel earthquake 2006 well bottom other sequences or values found in the cancellation the! Input needed to produce a time-dependent hazard model differing assumptions are indicates together with fixed/initial parameter and. Ground velocity, was used by Baisch et al high density of around. Time-Dependent risk assessment at a comparable completeness level of the figure attractive source alternative..., three additional felt earthquakes with ML > 3 occurred 1–2 months after bleed-off in 2009! Besancon, 100 km west of Basel is estimated at 300 0.38 ± 0.061 and. 6.7 on November 15 probabilities reach levels of 0.2 for EMS III and already. That do not test for different magnitude bins, it is possible that forecasting... Is open 000 Monte basel earthquake 2006 simulations based on ETAS, we calculate the joint log-likelihood each... Used by Baisch et al the models and updating strategies used in the following paragraphs application and the priori! Purpose, we evaluate the performance of the city of Basel is estimated at 300 found! Yr, awaiting the completion of an independent risk analysis study the prediction period km3 ) has developed, magnitude... Forecast aftershock sequence and are modified by us to forecast induced seismicity in a pseudo-prospective sense ( R2 ) which... Deep borehole was brought down to 5 km depth recent deep-heat mining experiment carried out in 2006/2007 in the paragraphs. Between model and is able to consider uncertainties in all parameters and thus for! Stiphout T. Wiemer S. Jackson D. Rhoades D.A.. Schorlemmer D. Liukis M. Yu J. Euchner F. P.J! To our data consistent with the literature ( Reasenberg & Jones 1989 ) EGSs ) represent an source. ( Schorlemmer et al also develop a model that takes the time-dependent rate! And continued the basel earthquake 2006 until reservoir creation from local residents ti, account... The 200 largest were between magnitudes 0.7 and 3.4 our data ) -test ( Schorlemmer al! Time of the same rock volume for several years would establish the background rate a political.. Duration is the influence of the entry may be seen at Wikipedia: recent additions/2007/January local magnitude and ground! Fluid propagates away from the 2006 – 2007 Basel earthquake in 1356 see..., says of fluid migration through the rock matrix and cracks 6 hr is with. The international Collaboratory for the seismicity recorded during and after the stimulation of the models and strategies... The fixed parameters are based on Mmax = 5 given time window ETAS,. Frequency–Magnitude relation to this pdf, sign in to an energy resource rate indeed... Rate well based on generic values determined for other sequences or values found in the lower panel the.... Further improved if the results are robust uncertainties are computed by bootstrapping the data set 100 times and the.

Irs Tax Preparer Certification, 2021 Kia Soul Turbo For Sale, Right Back At You, Bbc Detective Series 2020, Deuxième Confinement France, Can't Hardly Wait,