TY - JOUR A1 - Negi, Sanjay S. A1 - Paul, Ajay A1 - Cesca, Simone A1 - Kamal, A1 - Kriegerowski, Marius A1 - Mahesh, P. A1 - Gupta, Sandeep T1 - Crustal velocity structure and earthquake processes of Garhwal-Kumaun Himalaya: Constraints from regional waveform inversion and array beam modeling JF - Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth N2 - In order to understand present day earthquake kinematics at the Indian plate boundary, we analyse seismic broadband data recorded between 2007 and 2015 by the regional network in the Garhwal-Kumaun region, northwest Himalaya. We first estimate a local 1-D velocity model for the computation of reliable Green's functions, based on 2837 P-wave and 2680 S-wave arrivals from 251 well located earthquakes. The resulting 1-D crustal structure yields a 4-layer velocity model down to the depths of 20 km. A fifth homogeneous layer extends down to 46 km, constraining the Moho using travel-time distance curve method. We then employ a multistep moment tensor (MT) inversion algorithm to infer seismic moment tensors of 11 moderate earthquakes with Mw magnitude in the range 4.0–5.0. The method provides a fast MT inversion for future monitoring of local seismicity, since Green's functions database has been prepared. To further support the moment tensor solutions, we additionally model P phase beams at seismic arrays at teleseismic distances. The MT inversion result reveals the presence of dominant thrust fault kinematics persisting along the Himalayan belt. Shallow low and high angle thrust faulting is the dominating mechanism in the Garhwal-Kumaun Himalaya. The centroid depths for these moderate earthquakes are shallow between 1 and 12 km. The beam modeling result confirm hypocentral depth estimates between 1 and 7 km. The updated seismicity, constrained source mechanism and depth results indicate typical setting of duplexes above the mid crustal ramp where slip is confirmed along out-of-sequence thrusting. The involvement of Tons thrust sheet in out-of-sequence thrusting indicate Tons thrust to be the principal active thrust at shallow depth in the Himalayan region. Our results thus support the critical taper wedge theory, where we infer the microseismicity cluster as a result of intense activity within the Lesser Himalayan Duplex (LHD) system. KW - Critical taper wedge KW - Lesser Himalayan Duplex KW - Out-of-sequence thrust Y1 - 2017 U6 - https://doi.org/10.1016/j.tecto.2017.05.007 SN - 0040-1951 SN - 1879-3266 VL - 712 SP - 45 EP - 63 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Cesca, Simone A1 - Heimann, Sebastian A1 - Kriegerowski, Marius A1 - Saul, Joachim A1 - Dahm, Torsten T1 - Moment tensor inversion for nuclear explosions BT - what can we learn from the 6 January and 9 September 2016 Nuclear Tests, North Korea? JF - Seismological research letters N2 - Two nuclear explosions were carried out by the Democratic People’s Republic of North Korea in January and September 2016. Epicenters were located close to those of the 2006, 2009, and 2013 previous explosions. We perform a seismological analysis of the 2016 events combining the analysis of full waveforms at regional distances and seismic array beams at teleseismic distances. We estimate the most relevant source parameters, such as source depth, moment release, and full moment tensor (MT). The best MT solution can be decomposed into an isotropic source, directly related with the explosion and an additional deviatoric term, likely due to near‐source interactions with topographic and/or underground facilities features. We additionally perform an accurate resolution test to assess source parameters uncertainties and trade‐offs. This analysis sheds light on source parameters inconsistencies among studies on previous shallow explosive sources. The resolution of the true MT is hindered by strong source parameters trade‐offs, so that a broad range of well‐fitting MT solutions can be found, spanning from a dominant positive isotropic term to a dominant negative vertical compensated linear vector dipole. The true mechanism can be discriminated by additionally modeling first‐motion polarities at seismic arrays at teleseismic distances. A comparative assessment of the 2016 explosion with earlier nuclear tests documents similar vertical waveforms but a significant increase of amplitude for the 2016 explosions, which proves that the 9 September 2016 was the largest nuclear explosion ever performed in North Korea with a magnitude Mw 4.9 and a shallow depth of less than 2 km, although there are no proofs of a fusion explosion. Modeling transversal component waveforms suggests variable size and orientation of the double‐couple components of the 2009, 2013, and 2016 sources. Y1 - 2017 U6 - https://doi.org/10.1785/0220160139 SN - 0895-0695 SN - 1938-2057 VL - 88 IS - 2A SP - 300 EP - 310 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Petersen, Gesa Maria A1 - Cesca, Simone A1 - Kriegerowski, Marius T1 - Automated quality control for large seismic networks BT - implementation and application to the AlpArray seismic network JF - Seismological research letters N2 - As a consequence of the rapid growing worldwide seismic data set, a huge variety of automatized data-processing methods have been developed. To perform automatized waveform-based seismological studies aiming for magnitudes or source process inversion, it is crucial to identify network stations with erroneous transfer functions, gain factors, or component orientations. We developed a new tool dedicated to automated station quality control of dense seismic networks and arrays. The python-based AutoStatsQ toolbox uses the pyrocko seismic data-processing environment. The toolbox automatically downloads data and metadata for selected teleseismic events and performs different tests. As a result, relative gain factors, sensor orientation corrections, and reliable frequency bands are computed for all stations in a chosen time period. Relative gain factors are calculated for all stations and events in a time domain based on maximum P-phase amplitudes. A Rayleigh-wave polarization analysis is used to identify deviating sensor orientations. The power spectra of all stations in a given frequency range are compared with synthetic ones, accessing Global Centroid Moment Tensor (CMT) solutions. Frequency ranges of coinciding synthetic and recorded power spectral densities (PSDs) may serve as guidelines for choosing band-pass filters for moment tensor (MT) inversion and help confirm the corner frequency of the instrument. The toolbox was applied to the permanent and temporary AlpArray networks as well as to the denser SWATH-D network, a total of over 750 stations. Stations with significantly deviating gain factors were identified, as well as stations with inverse polarity and misorientations of the horizontal components. The tool can be used to quickly access network quality and to omit or correct stations before MT inversion. Electronic Supplement: List of teleseismic events and tables of median, mean, and standard deviation of relative gain factors, and figures of relative gain factors of all event-station pairs, waveform example showing inverse polarity of horizontal components on ZS.D125, histograms of median, mean, and standard deviation of the correction angles, examples of synthetic and recorded frequency spectra of ZS.D046 and NI.VINO. Y1 - 2019 U6 - https://doi.org/10.1785/0220180342 SN - 0895-0695 SN - 1938-2057 VL - 90 IS - 3 SP - 1177 EP - 1190 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Karamzadeh, Nasim Toularoud A1 - Kühn, Daniela A1 - Kriegerowski, Marius A1 - López-Comino, José Ángel A1 - Cesca, Simone A1 - Dahm, Torsten T1 - Small-aperture array as a tool to monitor fluid injection- and extraction-induced microseismicity BT - applications and recommendations JF - Acta Geophysica N2 - The monitoring of microseismicity during temporary human activities such as fluid injections for hydrofracturing, hydrothermal stimulations or wastewater disposal is a difficult task. The seismic stations often cannot be installed on hard rock, and at quiet places, noise is strongly increased during the operation itself and the installation of sensors in deep wells is costly and often not feasible. The combination of small-aperture seismic arrays with shallow borehole sensors offers a solution. We tested this monitoring approach at two different sites: (1) accompanying a fracking experiment in sedimentary shale at 4km depth and (2) above a gas field under depletion. The small-aperture arrays were planned according to theoretical wavenumber studies combined with simulations considering the local noise conditions. We compared array recordings with recordings available from shallow borehole sensors and give examples of detection and location performance. Although the high-frequency noise on the 50-m-deep borehole sensors was smaller compared to the surface noise before the injection experiment, the signals were highly contaminated during injection by the pumping activities. Therefore, a set of three small-aperture arrays at different azimuths was more suited to detect small events, since noise recorded on these arrays is uncorrelated with each other. Further, we developed recommendations for the adaptation of the monitoring concept to other sites experiencing induced seismicity. KW - Microseismic monitoring KW - Induced seismicity KW - Array seismology KW - Shallow borehole sensors Y1 - 2019 U6 - https://doi.org/10.1007/s11600-018-0231-1 SN - 1895-6572 SN - 1895-7455 VL - 67 IS - 1 SP - 311 EP - 326 PB - Springer CY - Cham ER - TY - JOUR A1 - Maghsoudi, Samira A1 - Cesca, Simone A1 - Hainzl, Sebastian A1 - Kaiser, Diethelm A1 - Becker, Dirk A1 - Dahm, Torsten T1 - Improving the estimation of detection probability and magnitude of completeness in strongly heterogeneous media, an application to acoustic emission (AE) JF - Geophysical journal international N2 - Reliable estimations of magnitude of completeness (M-c) are essential for a correct interpretation of seismic catalogues. The spatial distribution of M-c may be strongly variable and difficult to assess in mining environments, owing to the presence of galleries, cavities, fractured regions, porous media and different mineralogical bodies, as well as in consequence of inhomogeneous spatial distribution of the seismicity. We apply a 3-D modification of the probabilistic magnitude of completeness (PMC) method, which relies on the analysis of network detection capabilities. In our approach, the probability to detect an event depends on its magnitude, source receiver Euclidian distance and source receiver direction. The suggested method is proposed for study of the spatial distribution of the magnitude of completeness in a mining environment and here is applied to a 2-months acoustic emission (AE) data set recorded at the Morsleben salt mine, Germany. The dense seismic network and the large data set, which includes more than one million events, enable a detailed testing of the method. This method is proposed specifically for strongly heterogeneous media. Besides, it can also be used for specific network installations, with sensors with a sensitivity, dependent on the direction of the incoming wave (e.g. some piezoelectric sensors). In absence of strong heterogeneities, the standards PMC approach should be used. We show that the PMC estimations in mines strongly depend on the source receiver direction, and cannot be correctly accounted using a standard PMC approach. However, results can be improved, when adopting the proposed 3-D modification of the PMC method. Our analysis of one central horizontal and vertical section yields a magnitude of completeness of about M-c approximate to 1 (AE magnitude) at the centre of the network, which increases up to M-c approximate to 4 at further distances outside the network; the best detection performance is estimated for a NNE-SSE elongated region, which corresponds to the strike direction of the low-attenuating salt body. Our approach provides us with small-scale details about the capability of sensors to detect an earthquake, which can be linked to the presence of heterogeneities in specific directions. Reduced detection performance in presence of strong structural heterogeneities (cavities) is confirmed by synthetic waveform modelling in heterogeneous media. KW - Seismic attenuation KW - Statistical seismology Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggt049 SN - 0956-540X VL - 193 IS - 3 SP - 1556 EP - 1569 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Cesca, Simone A1 - Braun, Thomas A1 - Maccaferri, Francesco A1 - Passarelli, Luigi A1 - Rivalta, Eleonora A1 - Dahm, Torsten T1 - Source modelling of the M5-6 Emilia-Romagna, Italy, earthquakes (2012 May 20-29) JF - Geophysical journal international N2 - On 2012 May 20 and 29, two damaging earthquakes with magnitudes M-w 6.1 and 5.9, respectively, struck the Emilia-Romagna region in the sedimentary Po Plain, Northern Italy, causing 26 fatalities, significant damage to historical buildings and substantial impact to the economy of the region. The earthquake sequence included four more aftershocks with M-w, >= 5.0, all at shallow depths (about 7-9 km), with similar WNW-ESE striking reverse mechanism. The timeline of the sequence suggests significant static stress interaction between the largest events. We perform here a detailed source inversion, first adopting a point source approximation and considering pure double couple and full moment tensor source models. We compare different extended source inversion approaches for the two largest events, and find that the rupture occurred in both cases along a subhorizontal plane, dipping towards SSW Directivity is well detected for the May 20 main shock, indicating that the rupture propagated unilaterally towards SE. Based on the focal mechanism solution, we further estimate the co-seismic static stress change induced by the May 20 event. By using the rate-and-state model and a Poissonian earthquake occurrence, we infer that the second largest event of May 29 was induced with a probability in the range 0.2-0.4. This suggests that the segment of fault was already prone to rupture. Finally, we estimate peak ground accelerations for the two main events as occurred separately or simultaneously. For the scenario involving hypothetical rupture areas of both main events, we estimate M-w = 6.3 and an increase of ground acceleration by 50 per cent. The approach we propose may help to quantify rapidly which regions are invested by a significant increase of the hazard, bearing the potential for large aftershocks or even a second main shock. KW - Earthquake dynamics KW - Earthquake source observations Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggt069 SN - 0956-540X VL - 193 IS - 3 SP - 1658 EP - 1672 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Sen, Ali Tolga A1 - Cesca, Simone A1 - Bischoff, Monika A1 - Meier, Thomas A1 - Dahm, Torsten T1 - Automated full moment tensor inversion of coal mining-induced seismicity JF - Geophysical journal international N2 - Seismicity induced by coal mining in the Ruhr region, Germany, has been monitored continuously over the last 25 yr. In 2006, a dense temporary network (HAMNET) was deployed to locally monitor seismicity induced by longwall mining close to the town of Hamm. Between 2006 July and 2007 July, more than 7000 events with magnitudes M-L from -1.7 to 2.0 were detected. The spatiotemporal distribution of seismicity shows high correlation with the mining activity. In order to monitor rupture processes, we set up an automated source inversion routine and successfully perform double couple and full moment tensor (MT) inversions for more than 1000 events with magnitudes above M-L -0.5. The source inversion is based on a full waveform approach, both in the frequency and in the time domain, providing information about the centroid location, focal mechanism, scalar moment and full MT. Inversion results indicate a strong dominance of normal faulting focal mechanisms, with a steeper plane and a subhorizontal one. Fault planes are oriented parallel to the mining stopes. We classify the focal mechanisms based on their orientation and observe different frequency-magnitude distributions for families of events with different focal mechanisms; the overall frequency-magnitude distribution is not fitting the Gutenberg-Richter relation. Full MTs indicate that non-negligible opening tensile components accompanied normal faulting source mechanisms. Finally, extended source models are investigated for largest events. Results suggest that the rupture processes mostly occurred along the subvertical planes. KW - Geomechanics KW - Fracture and flow KW - Earthquake source observations KW - Seismicity and tectonics Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggt300 SN - 0956-540X SN - 1365-246X VL - 195 IS - 2 SP - 1267 EP - 1281 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Heimann, Sebastian A1 - Gonzalez, Alvaro A1 - Wang, Rongjiang A1 - Cesca, Simone A1 - Dahm, Torsten T1 - Seismic characterization of the Chelyabinsk Meteor's terminal explosion JF - Seismological research letters Y1 - 2013 U6 - https://doi.org/10.1785/0220130042 SN - 0895-0695 VL - 84 IS - 6 SP - 1021 EP - 1025 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Zhao, Peng A1 - Kuhn, Daniela A1 - Oye, Volker A1 - Cesca, Simone T1 - Evidence for tensile faulting deduced from full waveform moment tensor inversion during the stimulation of the Basel enhanced geothermal system JF - Geothermics : an international journal of geothermal research and its applications N2 - Our study presents the results of a moment tensor inversion of 19 microseismic events with M-L between 2.0 and 3.4, associated with the stimulation operation of an enhanced geothermal reservoir in Basel, Switzerland, in 2006. We adopt a three-step procedure to retrieve point source solution parameters based on full waveform inversion. The inversion is performed by fitting displacement amplitude spectra and displacement seismograms in the first and second step, respectively, assuming a double couple source model and thus obtaining focal solutions for all 19 events. Our results are in agreement with focal mechanisms from a previous study, which employed P wave first-motion polarities from more than 40 stations, whereas our solutions are achieved using full waveform data recorded by less than 10 surface stations. In the last step, a full moment tensor inversion is performed. The results from the moment tensor inversion show an improvement on the waveform fitting compared to the double couple models, which is verified by an F-test. We investigate the stability of the moment tensor solutions by employing different velocity models. The isotropic components of the moment tensor solutions of some events are not negligible, suggesting source volume changes due to fluid injection. Events with significant isotropic components occurred mainly during the stimulation phase and close to the injection well. On the other hand, events that occurred in the post-stimulation phase are predominantly pure shear failure and located further away from the well bore. These spatio-temporal patterns can be explained by the influence of pore pressure variations during and after the hydraulic stimulation at the geothermal site. (C) 2014 Elsevier Ltd. All rights reserved. KW - Induced seismicity KW - Hydraulic fracture KW - Crack opening and closure KW - Source mechanism KW - F-test Y1 - 2014 U6 - https://doi.org/10.1016/j.geothermics.2014.01.003 SN - 0375-6505 SN - 1879-3576 VL - 52 SP - 74 EP - 83 PB - Elsevier CY - Oxford ER - TY - JOUR A1 - Passarelli, Luigi A1 - Rivalta, Eleonora A1 - Cesca, Simone A1 - Aoki, Yosuke T1 - Stress changes, focal mechanisms, and earthquake scaling laws for the 2000 dike at Miyakejima (Japan) JF - Journal of geophysical research : Solid earth N2 - Faulting processes in volcanic areas result from a complex interaction of pressurized fluid-filled cracks and conduits with the host rock and local and regional tectonic setting. Often, volcanic seismicity is difficult to decipher in terms of the physical processes involved, and there is a need for models relating the mechanics of volcanic sources to observations. Here we use focal mechanism data of the energetic swarm induced by the 2000 dike intrusion at Miyakejima (Izu Archipelago, Japan), to study the relation between the 3-D dike-induced stresses and the characteristics of the seismicity. We perform a clustering analysis on the focal mechanism (FM) solutions and relate them to the dike stress field and to the scaling relationships of the earthquakes. We find that the strike and rake angles of the FMs are strongly correlated and cluster on bands in a strike-rake plot. We suggest that this is consistent with optimally oriented faults according to the expected pattern of Coulomb stress changes. We calculate the frequency-size distribution of the clustered sets finding that focal mechanisms with a large strike-slip component are consistent with the Gutenberg-Richter relation with a b value of about 1. Conversely, events with large normal faulting components deviate from the Gutenberg-Richter distribution with a marked roll-off on its right-hand tail, suggesting a lack of large-magnitude events (M-w>5.5). This may result from the interplay of the limited thickness and lower rock strength of the layer of rock above the dike, where normal faulting is expected, and lower stress levels linked to the faulting style and low confining pressure. KW - dike intrusion KW - dike-induced stresses KW - dike-induced seismicity KW - Miyakejima intrusion KW - Gutenberg-Richter relationship Y1 - 2015 U6 - https://doi.org/10.1002/2014JB011504 SN - 2169-9313 SN - 2169-9356 VL - 120 IS - 6 SP - 4130 EP - 4145 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Maghsoudi, Samira A1 - Hainzl, Sebastian A1 - Cesca, Simone A1 - Dahm, Torsten A1 - Kaiser, Diethelm T1 - Identification and characterization of growing large-scale en-echelon fractures in a salt mine JF - Geophysical journal international N2 - The spatiotemporal seismicity of acoustic emission (AE) events recorded in the Morsleben salt mine is investigated. Almost a year after backfilling of the cavities from 2003, microevents are distributed with distinctive stripe shapes above cavities at different depth levels. The physical forces driving the creation of these stripes are still unknown. This study aims to find the active stripes and track fracture developments over time by combining two different temporal and spatial clustering techniques into a single methodological approach. Anomalous seismicity parameters values like sharp b-value changes for two active stripes are good indicators to explain possible stress accumulation at the stripe tips. We identify the formation of two new seismicity stripes and show that the AE activities in active clusters are migrated mostly unidirectional to eastward and upward. This indicates that the growth of underlying macrofractures is controlled by the gradient of extensional stress. Studying size distribution characteristic in terms of frequency-magnitude distribution and b-value in active phase and phase with constant seismicity rate show that deviations from the Gutenberg-Richter power law can be explained by the inclusion of different activity phases: (1) the inactive period before the formation of macrofractures, which is characterized by a deficit of larger events (higher b-values) and (2) the period of fracture growth characterized by the occurrence of larger events (smaller b-values). KW - Earthquake source observations KW - Statistical seismology Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggt443 SN - 0956-540X SN - 1365-246X VL - 196 IS - 2 SP - 1092 EP - 1105 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Passarelli, Luigi A1 - Hainzl, Sebastian A1 - Cesca, Simone A1 - Maccaferri, Francesco A1 - Mucciarelli, Marco A1 - Rößler, Dirk A1 - Corbi, Fabio A1 - Dahm, Torsten A1 - Rivalta, Eleonora T1 - Aseismic transient driving the swarm-like seismic sequence in the Pollino range, Southern Italy JF - Geophysical journal international N2 - Tectonic earthquake swarms challenge our understanding of earthquake processes since it is difficult to link observations to the underlying physical mechanisms and to assess the hazard they pose. Transient forcing is thought to initiate and drive the spatio-temporal release of energy during swarms. The nature of the transient forcing may vary across sequences and range from aseismic creeping or transient slip to diffusion of pore pressure pulses to fluid redistribution and migration within the seismogenic crust. Distinguishing between such forcing mechanisms may be critical to reduce epistemic uncertainties in the assessment of hazard due to seismic swarms, because it can provide information on the frequency-magnitude distribution of the earthquakes (often deviating from the assumed Gutenberg-Richter relation) and on the expected source parameters influencing the ground motion (for example the stress drop). Here we study the ongoing Pollino range (Southern Italy) seismic swarm, a long-lasting seismic sequence with more than five thousand events recorded and located since October 2010. The two largest shocks (magnitude M-w = 4.2 and M-w = 5.1) are among the largest earthquakes ever recorded in an area which represents a seismic gap in the Italian historical earthquake catalogue. We investigate the geometrical, mechanical and statistical characteristics of the largest earthquakes and of the entire swarm. We calculate the focal mechanisms of the M-l > 3 events in the sequence and the transfer of Coulomb stress on nearby known faults and analyse the statistics of the earthquake catalogue. We find that only 25 per cent of the earthquakes in the sequence can be explained as aftershocks, and the remaining 75 per cent may be attributed to a transient forcing. The b-values change in time throughout the sequence, with low b-values correlated with the period of highest rate of activity and with the occurrence of the largest shock. In the light of recent studies on the palaeoseismic and historical activity in the Pollino area, we identify two scenarios consistent with the observations and our analysis: This and past seismic swarms may have been 'passive' features, with small fault patches failing on largely locked faults, or may have been accompanied by an 'active', largely aseismic, release of a large portion of the accumulated tectonic strain. Those scenarios have very different implications for the seismic hazard of the area. KW - Seismicity and tectonics KW - Statistical seismology KW - Dynamics: seismotectonics Y1 - 2015 U6 - https://doi.org/10.1093/gji/ggv111 SN - 0956-540X SN - 1365-246X VL - 201 IS - 3 SP - 1553 EP - 1567 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Maghsoudi, Samira A1 - Cesca, Simone A1 - Hainzl, Sebastian A1 - Dahm, Torsten A1 - Zöller, Gert A1 - Kaiser, Diethelm T1 - Maximum Magnitude of Completeness in a Salt Mine JF - Bulletin of the Seismological Society of America N2 - In this study, we analyze acoustic emission (AE) data recorded at the Morsleben salt mine, Germany, to assess the catalog completeness, which plays an important role in any seismicity analysis. We introduce the new concept of a magnitude completeness interval consisting of a maximum magnitude of completeness (M-c(max)) in addition to the well-known minimum magnitude of completeness. This is required to describe the completeness of the catalog, both for the smallest events (for which the detection performance may be low) and for the largest ones (which may be missed because of sensors saturation). We suggest a method to compute the maximum magnitude of completeness and calculate it for a spatial grid based on (1) the prior estimation of saturation magnitude at each sensor, (2) the correction of the detection probability function at each sensor, including a drop in the detection performance when it saturates, and (3) the combination of detection probabilities of all sensors to obtain the network detection performance. The method is tested using about 130,000 AE events recorded in a period of five weeks, with sources confined within a small depth interval, and an example of the spatial distribution of M-c(max) is derived. The comparison between the spatial distribution of M-c(max) and of the maximum possible magnitude (M-max), which is here derived using a recently introduced Bayesian approach, indicates that M-max exceeds M-c(max) in some parts of the mine. This suggests that some large and important events may be missed in the catalog, which could lead to a bias in the hazard evaluation. Y1 - 2015 U6 - https://doi.org/10.1785/0120140039 SN - 0037-1106 SN - 1943-3573 VL - 105 IS - 3 SP - 1491 EP - 1501 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Pro, C. A1 - Buforn, Elisa A1 - Cesca, Simone A1 - Sanz de Galdeano, C. A1 - Udias, A. T1 - Rupture process of the Lorca (southeast Spain) 11 May 2011 (M (w)=5.1) earthquake JF - Journal of seismology N2 - On 11 May 2011, a M (w) = 5.1 earthquake shook the town of Lorca (SE Spain) causing a disproportionately large damage for its magnitude. In order to contribute to knowledge of the behavior of the active faults present in the region and define the parameters which control their motion, we made a detailed study of the rupture process of this earthquake from inversion of body waves at regional and teleseismic distances. Ground motion displacements obtained in this way are in agreement with near-field strong motion data and GPS observations recorded in Lorca. We have obtained a partly bilateral rupture propagating to WSW (238A degrees, 54A degrees, 59A degrees) with 27 cm of maximum slip and shallow focus (4 km). The fault plane orientation corresponds to that of the Cejo de los Enamorados Fault located NE of the Lorca town and parallel to the Alhama de Murcia Fault. The distribution of slip on the fault plane can explain the lack of any observed surface rupture as we found that the rupture started at 4-km depth along a plane dipping at 54A degrees, with motion propagating upward to stop at 1.5 km below the surface. The strong motion and GPS data recorded near the epicenter are in agreement with the maximum slip on the fault. Directivity effects and the extreme shallowness of the rupture could explain the considerable damage that the earthquake caused in the town of Lorca. KW - Seismic slip distribution KW - Directivity KW - Ground motion and GPS Y1 - 2014 U6 - https://doi.org/10.1007/s10950-014-9421-8 SN - 1383-4649 SN - 1573-157X VL - 18 IS - 3 SP - 481 EP - 495 PB - Springer CY - Dordrecht ER - TY - JOUR A1 - Cesca, Simone A1 - Sen, Ali Tolga A1 - Dahm, Torsten T1 - Seismicity monitoring by cluster analysis of moment tensors JF - Geophysical journal international N2 - We suggest a new clustering approach to classify focal mechanisms from large moment tensor catalogues, with the purpose of automatically identify families of earthquakes with similar source geometry, recognize the orientation of most active faults, and detect temporal variations of the rupture processes. The approach differs in comparison to waveform similarity methods since clusters are detected even if they occur in large spatial distances. This approach is particularly helpful to analyse large moment tensor catalogues, as in microseismicity applications, where a manual analysis and classification is not feasible. A flexible algorithm is here proposed: it can handle different metrics, norms, and focal mechanism representations. In particular, the method can handle full moment tensor or constrained source model catalogues, for which different metrics are suggested. The method can account for variable uncertainties of different moment tensor components. We verify the method with synthetic catalogues. An application to real data from mining induced seismicity illustrates possible applications of the method and demonstrate the cluster detection and event classification performance with different moment tensor catalogues. Results proof that main earthquake source types occur on spatially separated faults, and that temporal changes in the number and characterization of focal mechanism clusters are detected. We suggest that moment tensor clustering can help assessing time dependent hazard in mines. KW - Persistence KW - memory KW - correlations KW - clustering KW - Earthquake source observations Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggt492 SN - 0956-540X SN - 1365-246X VL - 196 IS - 3 SP - 1813 EP - 1826 PB - Oxford Univ. Press CY - Oxford ER - TY - INPR A1 - Cesca, Simone A1 - Dost, Bernard A1 - Oth, Adrien T1 - Preface to the special issue "Triggered and induced seismicity: probabilities and discrimination" T2 - Journal of seismology Y1 - 2013 U6 - https://doi.org/10.1007/s10950-012-9338-z SN - 1383-4649 VL - 17 IS - 1 SP - 1 EP - 4 PB - Springer CY - Dordrecht ER - TY - GEN A1 - Dahm, Torsten A1 - Becker, Dirk A1 - Bischoff, Monika A1 - Cesca, Simone A1 - Dost, B. A1 - Fritschen, R. A1 - Hainzl, Sebastian A1 - Klose, C. D. A1 - Kuhn, D. A1 - Lasocki, S. A1 - Meier, Thomas A1 - Ohrnberger, Matthias A1 - Rivalta, Eleonora A1 - Wegler, Ulrich A1 - Husen, Stephan T1 - Recommendation for the discrimination of human-related and natural seismicity T2 - Journal of seismology N2 - Various techniques are utilized by the seismological community, extractive industries, energy and geoengineering companies to identify earthquake nucleation processes in close proximity to engineering operation points. These operations may comprise fluid extraction or injections, artificial water reservoir impoundments, open pit and deep mining, deep geothermal power generations or carbon sequestration. In this letter to the editor, we outline several lines of investigation that we suggest to follow to address the discrimination problem between natural seismicity and seismic events induced or triggered by geoengineering activities. These suggestions have been developed by a group of experts during several meetings and workshops, and we feel that their publication as a summary report is helpful for the geoscientific community. Specific investigation procedures and discrimination approaches, on which our recommendations are based, are also published in this Special Issue (SI) of Journal of Seismology. KW - Triggered seismicity KW - Induced seismicity Y1 - 2013 U6 - https://doi.org/10.1007/s10950-012-9295-6 SN - 1383-4649 VL - 17 IS - 1 SP - 197 EP - 202 PB - Springer CY - Dordrecht ER - TY - GEN A1 - Passarelli, Luigi A1 - Hainzl, Sebastian A1 - Cesca, Simone A1 - Maccaferri, Francesco A1 - Mucciarelli, Marco A1 - Roessler, Dirk A1 - Corbi, Fabio A1 - Dahm, Torsten A1 - Rivalta, Eleonora T1 - Aseismic transient driving the swarm-like seismic sequence in the Pollino range, Southern Italy (vol 201, pg 1553, 2015) T2 - Geophysical journal international KW - Seismicity and tectonics KW - Statistical seismology KW - Dynamics: seismotectonics Y1 - 2016 U6 - https://doi.org/10.1093/gji/ggv425 SN - 0956-540X SN - 1365-246X VL - 204 SP - 365 EP - 365 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - del Fresno, Carmen A1 - Dominguez Cerdena, Itahiza A1 - Cesca, Simone A1 - Buforn, Elisa T1 - The 8 October 2011 Earthquake at El Hierro (M-w 4.0): Focal Mechanisms of the Mainshock and Its Foreshocks JF - Bulletin of the Seismological Society of America N2 - We have studied the focal mechanism of an M-w 4.0 earthquake that occurred on 8 October 2011 in the southwest of El Hierro (Canary Islands), the largest shock of the swarm that preceded the submarine eruption of El Hierro 2011-2012. The joint focal mechanism solution of 34 foreshocks has also been obtained. The results for the mainshock are consistent with a pure double-couple mechanism of a strike-slip motion with a reverse component and a focal depth of 12-13 km. The stress pattern obtained from the focal mechanism indicates horizontal compression in a north-northwest-south-southeast direction, parallel to the southern ridge of the island, and a quasi-horizontal extension in an east-west direction. Similar stress pattern is derived from the joint solution of the foreshocks. The occurrence of this family of earthquakes at the moment of the maximum strain rate of the pre-eruptive swarm suggests that their rupture process is related to tectonic stress, which led to the eruption only two days later, 5 km away from the mainshock epicenter. Y1 - 2015 U6 - https://doi.org/10.1785/0120140151 SN - 0037-1106 SN - 1943-3573 VL - 105 IS - 1 SP - 330 EP - 340 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Grigoli, Francesco A1 - Cesca, Simone A1 - Rinaldi, Antonio Pio A1 - Manconi, Andrea A1 - Lopez-Comino, José Ángel A1 - Clinton, John F. A1 - Westaway, Rob A1 - Cauzzi, Carlo A1 - Dahm, Torsten A1 - Wiemer, Stefan T1 - The November 2017 M-w 5.5 Pohang earthquake BT - a possible case of induced seismicity in South Korea JF - Science N2 - The moment magnitude (M-w) 5.5 earthquake that struck South Korea in November 2017 was one of the largest and most damaging events in that country over the past century. Its proximity to an enhanced geothermal system site, where high-pressure hydraulic injection had been performed during the previous 2 years, raises the possibility that this earthquake was anthropogenic. We have combined seismological and geodetic analyses to characterize the mainshock and its largest aftershocks, constrain the geometry of this seismic sequence, and shed light on its causal factors. According to our analysis, it seems plausible that the occurrence of this earthquake was influenced by the aforementioned industrial activities. Finally, we found that the earthquake transferred static stress to larger nearby faults, potentially increasing the seismic hazard in the area. Y1 - 2018 U6 - https://doi.org/10.1126/science.aat2010 SN - 0036-8075 SN - 1095-9203 VL - 360 IS - 6392 SP - 1003 EP - 1006 PB - American Assoc. for the Advancement of Science CY - Washington ER - TY - JOUR A1 - Cesca, Simone A1 - Grigoli, Francesco A1 - Heimann, Sebastian A1 - Dahm, Torsten A1 - Kriegerowski, Marius A1 - Sobiesiak, M. A1 - Tassara, C. A1 - Olcay, M. T1 - The M-w 8.1 2014 Iquique, Chile, seismic sequence: a tale of foreshocks and aftershocks JF - Geophysical journal international N2 - The 2014 April 1, M-w 8.1 Iquique (Chile) earthquake struck in the Northern Chile seismic gap. With a rupture length of less than 200 km, it left unbroken large segments of the former gap. Early studies were able to model the main rupture features but results are ambiguous with respect to the role of aseismic slip and left open questions on the remaining hazard at the Northern Chile gap. A striking observation of the 2014 earthquake has been its extensive preparation phase, with more than 1300 events with magnitude above M-L 3, occurring during the 15 months preceding the main shock. Increasing seismicity rates and observed peak magnitudes accompanied the last three weeks before the main shock. Thanks to the large data sets of regional recordings, we assess the precursor activity, compare foreshocks and aftershocks and model rupture preparation and rupture effects. To tackle inversion challenges for moderate events with an asymmetric network geometry, we use full waveforms techniques to locate events, map the seismicity rate and derive source parameters, obtaining moment tensors for more than 300 events (magnitudes M-w 4.0-8.1) in the period 2013 January 1-2014 April 30. This unique data set of fore- and aftershocks is investigated to distinguish rupture process models and models of strain and stress rotation during an earthquake. Results indicate that the spatial distributions of foreshocks delineated the shallower part of the rupture areas of the main shock and its largest aftershock, well matching the spatial extension of the aftershocks cloud. Most moment tensors correspond to almost pure double couple thrust mechanisms, consistent with the slab orientation. Whereas no significant differences are observed among thrust mechanisms in different areas, nor among thrust foreshocks and aftershocks, the early aftershock sequence is characterized by the presence of normal fault mechanisms, striking parallel to the trench but dipping westward. These events likely occurred in the shallow wedge structure close to the slab interface and are consequence of the increased extensional stress in this region after the largest events. The overall stress inversion result suggests a minor stress rotation after the main shock, but a significant release of the deviatoric stress. The temporal change in the distribution of focal mechanisms can also be explained in terms of the spatial heterogeneity of the stress field: under such interpretation, the potential of a large megathrust earthquake breaking a larger segment offshore Northern Chile remains high. KW - Earthquake source observations KW - South America Y1 - 2016 U6 - https://doi.org/10.1093/gji/ggv544 SN - 0956-540X SN - 1365-246X VL - 204 SP - 1766 EP - 1780 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Cesca, Simone A1 - Grigoli, Francesco A1 - Heimann, Sebastian A1 - Gonzalez, Alvaro A1 - Buforn, Elisa A1 - Maghsoudi, Samira A1 - Blanch, Estefania A1 - Dahm, Torsten T1 - The 2013 September-October seismic sequence offshore Spain: a case of seismicity triggered by gas injection? JF - Geophysical journal international N2 - A spatially localized seismic sequence originated few tens of kilometres offshore the Mediterranean coast of Spain, close to the Ebro river delta, starting on 2013 September 5, and lasting at least until 2013 October. The sequence culminated in a maximal moment magnitude M-w 4.3 earthquake, on 2013 October 1. The most relevant seismogenic feature in the area is the Fosa de Amposta fault system, which includes different strands mapped at different distances to the coast, with a general NE-SW orientation, roughly parallel to the coastline. However, no significant known historical seismicity has involved this fault system in the past. The epicentral region is also located near the offshore platform of the Castor project, where gas is conducted through a pipeline from mainland and where it was recently injected in a depleted oil reservoir, at about 2 km depth. We analyse the temporal evolution of the seismic sequence and use full waveform techniques to derive absolute and relative locations, estimate depths and focal mechanisms for the largest events in the sequence (with magnitude mbLg larger than 3), and compare them to a previous event (2012 April 8, mbLg 3.3) taking place in the same region prior to the gas injection. Moment tensor inversion results show that the overall seismicity in this sequence is characterized by oblique mechanisms with a normal fault component, with a 30A degrees low-dip angle plane oriented NNE-SSW and a subvertical plane oriented NW-SE. The combined analysis of hypocentral location and focal mechanisms could indicate that the seismic sequence corresponds to rupture processes along shallow low-dip surfaces, which could have been triggered by the gas injection in the reservoir, and excludes the activation of the Amposta fault, as its known orientation is inconsistent with focal mechanism results. An alternative scenario includes the iterated triggering of a system of steep faults oriented NW-SE, which were identified by prior marine seismics investigations. KW - Earthquake dynamics KW - Earthquake source observations Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggu172 SN - 0956-540X SN - 1365-246X VL - 198 IS - 2 SP - 941 EP - 953 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Grigoli, Francesco A1 - Cesca, Simone A1 - Amoroso, Ortensia A1 - Emolo, Antonio A1 - Zollo, Aldo A1 - Dahm, Torsten T1 - Automated seismic event location by waveform coherence analysis JF - Geophysical journal international N2 - Automated location of seismic events is a very important task in microseismic monitoring operations as well for local and regional seismic monitoring. Since microseismic records are generally characterized by low signal-to-noise ratio, automated location methods are requested to be noise robust and sufficiently accurate. Most of the standard automated location routines are based on the automated picking, identification and association of the first arrivals of P and S waves and on the minimization of the residuals between theoretical and observed arrival times of the considered seismic phases. Although current methods can accurately pick P onsets, the automatic picking of the S onset is still problematic, especially when the P coda overlaps the S wave onset. In this paper, we propose a picking free earthquake location method based on the use of the short-term-average/long-term-average (STA/LTA) traces at different stations as observed data. For the P phases, we use the STA/LTA traces of the vertical energy function, whereas for the S phases, we use the STA/LTA traces of a second characteristic function, which is obtained using the principal component analysis technique. In order to locate the seismic event, we scan the space of possible hypocentral locations and origin times, and stack the STA/LTA traces along the theoretical arrival time surface for both P and S phases. Iterating this procedure on a 3-D grid, we retrieve a multidimensional matrix whose absolute maximum corresponds to the spatial coordinates of the seismic event. A pilot application was performed in the Campania-Lucania region (southern Italy) using a seismic network (Irpinia Seismic Network) with an aperture of about 150 km. We located 196 crustal earthquakes (depth < 20 km) with magnitude range 1.1 < M-L < 2.7. A subset of these locations were compared with accurate manual locations refined by using a double-difference technique. Our results indicate a good agreement with manual locations. Moreover, our method is noise robust and performs better than classical location methods based on the automatic picking of the P and S waves first arrivals. KW - Time-series analysis KW - Inverse theory KW - Earthquake source observations KW - Seismicity and tectonics KW - Early warning Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggt477 SN - 0956-540X SN - 1365-246X VL - 196 IS - 3 SP - 1742 EP - 1753 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Niemz, Peter A1 - Dahm, Torsten A1 - Milkereit, Claus A1 - Cesca, Simone A1 - Petersen, Gesa Maria A1 - Zang, Arno T1 - Insights into hydraulic fracture growth gained from a joint analysis of seismometer-derived tilt signals and scoustic emissions JF - Journal of geophysical research : Solid earth N2 - Hydraulic fracturing is performed to enhance rock permeability, for example, in the frame of geothermal energy production or shale gas exploitation, and can potentially trigger induced seismicity. The tracking of increased permeabilities and the fracturing extent is often based on the microseismic event distribution within the stimulated rock volume, but it is debated whether the microseismic activity adequately depicts the fracture formation. We are able to record tilt signals that appear as long-period transients (<180 s) on two broadband seismometers installed close (17-72 m) to newly formed, meter-scale hydraulic fractures. With this observation, we can overcome the limitations of the microseismic monitoring alone and verify the fracture mapping. Our analysis for the first time combines a catalog of previously analyzed acoustic emissions ([AEs] durations of 20 ms), indirectly mapping the fractures, with unique tilt signals, that provide independent, direct insights into the deformation of the rock. The analysis allows to identify different phases of the fracturing process including the (re)opening, growth, and aftergrowth of fractures. Further, it helps to differentiate between the formation of complex fracture networks and single macrofractures, and it validates the AE fracture mapping. Our findings contribute to a better understanding of the fracturing processes, which may help to reduce fluid-injection-induced seismicity and validate efficient fracture formation.
Plain Language Summary Hydraulic fracturing (HF) describes the opening of fractures in rocks by injecting fluids under high pressure. The new fractures not only can facilitate the extraction of shale gas but can also be used to heat up water in the subsurface in enhanced geothermal systems, a corner stone of renewable energy production. The fracture formation is inherently accompanied by small, nonfelt earthquakes (microseismic events). Occasionally, larger events felt by the population can be induced by the subsurface operations. Avoiding such events is important for the acceptance of HF operations and requires a detailed knowledge about the fracture formation. We jointly analyze two very different data sets recorded during mine-scale HF experiments: (a) the tilting of the ground caused by the opening of the fractures, as recorded by broadband seismometers-usually deployed for earthquake monitoring-installed close to the experiments and (b) a catalog of acoustic emissions, seismic signals of few milliseconds emitted by tiny cracks around the forming hydraulic fracture. The novel joint analysis allows to characterize the fracturing processes in greater detail, contributing to the understanding of the physical processes, which may help to understand fluid-injection-induced seismicity and validate the formation of hydraulic fractures. KW - hydraulic fracturing KW - fracture growth KW - tilt KW - acoustic emissions KW - injections KW - broadband seismometer Y1 - 2021 U6 - https://doi.org/10.1029/2021JB023057 SN - 2169-9313 SN - 2169-9356 VL - 126 IS - 12 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Grigoli, Francesco A1 - Cesca, Simone A1 - Vassallo, Maurizio A1 - Dahm, Torsten T1 - Automated seismic event location by travel-time stacking an application to mining induced seismicity JF - Seismological research letters Y1 - 2013 U6 - https://doi.org/10.1785/0220120191 SN - 0895-0695 VL - 84 IS - 4 SP - 666 EP - 677 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Dahm, Torsten A1 - Cesca, Simone A1 - Hainzl, Sebastian A1 - Braun, Thomas A1 - Krüger, Frank T1 - Discrimination between induced, triggered, and natural earthquakes close to hydrocarbon reservoirs: A probabilistic approach based on the modeling of depletion-induced stress changes and seismological source parameters JF - Journal of geophysical research : Solid earth N2 - Earthquakes occurring close to hydrocarbon fields under production are often under critical view of being induced or triggered. However, clear and testable rules to discriminate the different events have rarely been developed and tested. The unresolved scientific problem may lead to lengthy public disputes with unpredictable impact on the local acceptance of the exploitation and field operations. We propose a quantitative approach to discriminate induced, triggered, and natural earthquakes, which is based on testable input parameters. Maxima of occurrence probabilities are compared for the cases under question, and a single probability of being triggered or induced is reported. The uncertainties of earthquake location and other input parameters are considered in terms of the integration over probability density functions. The probability that events have been human triggered/induced is derived from the modeling of Coulomb stress changes and a rate and state-dependent seismicity model. In our case a 3-D boundary element method has been adapted for the nuclei of strain approach to estimate the stress changes outside the reservoir, which are related to pore pressure changes in the field formation. The predicted rate of natural earthquakes is either derived from the background seismicity or, in case of rare events, from an estimate of the tectonic stress rate. Instrumentally derived seismological information on the event location, source mechanism, and the size of the rupture plane is of advantage for the method. If the rupture plane has been estimated, the discrimination between induced or only triggered events is theoretically possible if probability functions are convolved with a rupture fault filter. We apply the approach to three recent main shock events: (1) the M-w 4.3 Ekofisk 2001, North Sea, earthquake close to the Ekofisk oil field; (2) the M-w 4.4 Rotenburg 2004, Northern Germany, earthquake in the vicinity of the Sohlingen gas field; and (3) the M-w 6.1 Emilia 2012, Northern Italy, earthquake in the vicinity of a hydrocarbon reservoir. The three test cases cover the complete range of possible causes: clearly human induced, not even human triggered, and a third case in between both extremes. KW - induced seismicity KW - probabilistic discrimination KW - hydrocarbon field KW - triggered earthquake KW - seismic hazard KW - earthquake Y1 - 2015 U6 - https://doi.org/10.1002/2014JB011778 SN - 2169-9313 SN - 2169-9356 VL - 120 IS - 4 SP - 2491 EP - 2509 PB - American Geophysical Union CY - Washington ER - TY - GEN A1 - Kriegerowski, Marius A1 - Cesca, Simone A1 - Ohrnberger, Matthias A1 - Dahm, Torsten A1 - Krüger, Frank T1 - Event couple spectral ratio Q method for earthquake clusters BT - application to northwest Bohemia T2 - Postprints der Universität Potsdam Mathematisch-Naturwissenschaftliche Reihe N2 - We develop an amplitude spectral ratio method for event couples from clustered earthquakes to estimate seismic wave attenuation (Q-1) in the source volume. The method allows to study attenuation within the source region of earthquake swarms or aftershocks at depth, independent of wave path and attenuation between source region and surface station. We exploit the high-frequency slope of phase spectra using multitaper spectral estimates. The method is tested using simulated full wave-field seismograms affected by recorded noise and finite source rupture. The synthetic tests verify the approach and show that solutions are independent of focal mechanisms but also show that seismic noise may broaden the scatter of results. We apply the event couple spectral ratio method to northwest Bohemia, Czech Republic, a region characterized by the persistent occurrence of earthquake swarms in a confined source region at mid-crustal depth. Our method indicates a strong anomaly of high attenuation in the source region of the swarm with an averaged attenuation factor of Qp < 100. The application to S phases fails due to scattered P-phase energy interfering with S phases. The Qp anomaly supports the common hypothesis of highly fractured and fluid saturated rocks in the source region of the swarms in northwest Bohemia. However, high temperatures in a small volume around the swarms cannot be excluded to explain our observations. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 683 KW - west bohemia KW - attenuation tomography KW - swarm earthquakes KW - focal zone KW - parameters KW - locations KW - fault Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-426029 IS - 683 ER - TY - JOUR A1 - Kriegerowski, Marius A1 - Cesca, Simone A1 - Ohrnberger, Matthias A1 - Dahm, Torsten A1 - Krüger, Frank T1 - Event couple spectral ratio Q method for earthquake clusters BT - application to northwest Bohemia JF - Solid Earth N2 - We develop an amplitude spectral ratio method for event couples from clustered earthquakes to estimate seismic wave attenuation (Q-1) in the source volume. The method allows to study attenuation within the source region of earthquake swarms or aftershocks at depth, independent of wave path and attenuation between source region and surface station. We exploit the high-frequency slope of phase spectra using multitaper spectral estimates. The method is tested using simulated full wave-field seismograms affected by recorded noise and finite source rupture. The synthetic tests verify the approach and show that solutions are independent of focal mechanisms but also show that seismic noise may broaden the scatter of results. We apply the event couple spectral ratio method to northwest Bohemia, Czech Republic, a region characterized by the persistent occurrence of earthquake swarms in a confined source region at mid-crustal depth. Our method indicates a strong anomaly of high attenuation in the source region of the swarm with an averaged attenuation factor of Qp < 100. The application to S phases fails due to scattered P-phase energy interfering with S phases. The Qp anomaly supports the common hypothesis of highly fractured and fluid saturated rocks in the source region of the swarms in northwest Bohemia. However, high temperatures in a small volume around the swarms cannot be excluded to explain our observations. KW - west bohemia KW - attenuation tomography KW - swarm earthquakes KW - focal zone KW - parameters KW - locations KW - fault Y1 - 2019 U6 - https://doi.org/10.5194/se-10-317-2019 SN - 1869-9529 IS - 10 SP - 317 EP - 328 PB - Copernicus Publications CY - Göttingen ER - TY - JOUR A1 - Niemz, Peter A1 - Cesca, Simone A1 - Heimann, Sebastian A1 - Grigoli, Francesco A1 - von Specht, Sebastian A1 - Hammer, Conny A1 - Zang, Arno A1 - Dahm, Torsten T1 - Full-waveform-based characterization of acoustic emission activity in a mine-scale experiment BT - a comparison of conventional and advanced hydraulic fracturing schemes JF - Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society N2 - Understanding fracturing processes and the hydromechanical relation to induced seismicity is a key question for enhanced geothermal systems (EGS). Commonly massive fluid injection, predominately causing hydroshearing, are used in large-scale EGS but also hydraulic fracturing approaches were discussed. To evaluate the applicability of hydraulic fracturing techniques in EGS, six in situ, multistage hydraulic fracturing experiments with three different injection schemes were performed under controlled conditions in crystalline rock at the Aspo Hard Rock Laboratory (Sweden). During the experiments the near-field ground motion was continuously recorded by 11 piezoelectric borehole sensors with a sampling rate of 1 MHz. The sensor network covered a volume of 30x30x30 m around a horizontal, 28-m-long injection borehole at a depth of 410 m. To extract and characterize massive, induced, high-frequency acoustic emission (AE) activity from continuous recordings, a semi-automated workflow was developed relying on full waveform based detection, classification and location procedures. The approach extended the AE catalogue from 196 triggered events in previous studies to more than 19600 located AEs. The enhanced catalogue, for the first time, allows a detailed analysis of induced seismicity during single hydraulic fracturing experiments, including the individual fracturing stages and the comparison between injection schemes. Beside the detailed study of the spatio-temporal patterns, event clusters and the growth of seismic clouds, we estimate relative magnitudes and b-values of AEs for conventional, cyclic progressive and dynamic pulse injection schemes, the latter two being fatigue hydraulic fracturing techniques. While the conventional fracturing leads to AE patterns clustered in planar regions, indicating the generation of a single main fracture plane, the cyclic progressive injection scheme results in a more diffuse, cloud-like AE distribution, indicating the activation of a more complex fracture network. For a given amount of hydraulic energy (pressure multiplied by injected volume) pumped into the system, the cyclic progressive scheme is characterized by a lower rate of seismicity, lower maximum magnitudes and significantly larger b-values, implying an increased number of small events relative to the large ones. To our knowledge, this is the first direct comparison of high resolution seismicity in a mine-scale experiment induced by different hydraulic fracturing schemes. KW - Fracture and flow KW - Spatial analysis KW - Statistical methods KW - Time-series analysis KW - Induced seismicity Y1 - 2020 U6 - https://doi.org/10.1093/gji/ggaa127 SN - 0955-419X SN - 1365-246X VL - 222 IS - 1 SP - 189 EP - 206 PB - Oxford Univ. Press CY - Oxford ER -