Refine
Year of publication
Language
- English (65)
Is part of the Bibliography
- yes (65) (remove)
Keywords
- Earthquake source observations (11)
- Seismicity and tectonics (6)
- Induced seismicity (5)
- Time-series analysis (5)
- Statistical seismology (4)
- Earthquake dynamics (3)
- Fracture and flow (3)
- fault (3)
- induced seismicity (3)
- Body waves (2)
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.
Rupture directivity, implying a predominant earthquake rupture propagation direction, is typically inferred upon the identification of 2D azimuthal patterns of seismic observations for weak to large earthquakes using surface-monitoring networks. However, the recent increase of 3D monitoring networks deployed in the shallow subsurface and underground laboratories toward the monitoring of microseismicity allows to extend the directivity analysis to 3D modeling, beyond the usual range of magnitudes. The high-quality full waveforms recorded for the largest, decimeter-scale acoustic emission (AE) events during a meter-scale hydraulic fracturing experiment in granites at similar to 410 m depth allow us to resolve the apparent durations observed at each AE sensor to analyze 3D-directivity effects. Unilateral and (asymmetric) bilateral ruptures are then characterized by the introduction of a parameter kappa, representing the angle between the directivity vector and the station vector. While the cloud of AE activity indicates the planes of the hydrofractures, the resolved directivity vectors show off-plane orientations, indicating that rupture planes of microfractures on a scale of centimeters have different geometries. Our results reveal a general alignment of the rupture directivity with the orientation of the minimum horizontal stress, implying that not only the slip direction but also the fracture growth produced by the fluid injections is controlled by the local stress conditions.
A review of source models to further the understanding of the seismicity of the Groningen field
(2022)
The occurrence of felt earthquakes due to gas production in Groningen has initiated numerous studies and model attempts to understand and quantify induced seismicity in this region. The whole bandwidth of available models spans the range from fully deterministic models to purely empirical and stochastic models. In this article, we summarise the most important model approaches, describing their main achievements and limitations. In addition, we discuss remaining open questions and potential future directions of development.
Salt diapirs are common features of sedimentary basins. If close to the surface, they can bear a significant hazard due to possible dissolution sinkholes, karst formation and collapse dolines or their influence on ground water chemistry. We investigate the potential of ambient vibration techniques to map the 3-D roof morphology of shallow salt diapirs. Horizontal-to-vertical (H/V) spectral peaks are derived at more than 900 positions above a shallow diapir beneath the city area of Hamburg, Germany, and are used to infer the depth of the first strong impedance contrast. In addition, 15 small-scale array measurements are conducted at different positions in order to compute frequency-dependent phase velocities of Rayleigh waves between 0.5 and 25 Hz. The dispersion curves are inverted together with the H/V peak frequency to obtain shear-wave velocity profiles. Additionally, we compare the morphology derived from H/V and array measurements to borehole lithology and a gravity-based 3-D model of the salt diapir. Both methods give consistent results in agreement with major features indicated by the independent data. An important result is that H/V and array measurements are better suited to identify weathered gypsum caprocks or gypsum floaters, while gravity-derived models better sample the interface between sediments and homogeneous salt. We further investigate qualitatively the influence of the 3-D subsurface topography of the salt diapir on the validity of local 1-D inversion results from ambient vibration dispersion curve inversion.
A mid-aperture broad-band test array (OBS array DOCTAR) was deployed from June 2011 to April 2012 about 100 km north of the Gloria fault in the Eastern North Atlantic in about 5000 m water depth. In addition arrays were installed on Madeira Island and in western Portugal mainland. For the first time in the Eastern North Atlantic, we recorded a large number of high frequency Po and So waves from local and regional small and moderate earthquakes (M-L < 4). An incoherent beamforming method was adapted to scan continuous data for such Po and So arrivals applying a sliding window waveform migration and frequency-wavenumber technique. We identify about 320 Po and 1550 So arrivals and compare the phase onsets with the ISC catalogue (ISC 2015) for the same time span. Up to a distance of 6 degrees to the DOCTAR stations all events listed in the ISC catalogue could be associated to Po and So phases. Arrivals from events in more than 10 degrees distance could be identified only in some cases. Only few Po and/or So arrivals were detected for earthquakes from the European and African continental area, the continental shelf regions and for earthquakes within or northwest of the Azores plateau. Unexpectedly, earthquake clusters are detected within the oceanic plates north and south of the Gloria fault and far from plate boundaries, indicating active intraplate structures. We also observe and locate numerous small magnitude earthquakes on the segment of the Gloria fault directly south of DOCTAR, which likely coincides with the rupture of the 25 November 1941 event. Local small magnitude earthquakes located beneath DOCTAR show hypocentres up to 30 km depth and strike-slip focal mechanisms. A comparison with detections at temporary mid-aperture arrays on Madeira and in western Portugal shows that the deep ocean array performs much better than the island and the continental array regarding the detection threshold for events in the oceanic plates.
We conclude that sparsely distributed mid-aperture seismic arrays in the deep ocean could decrease the detection and location threshold for seismicity with M-L < 4 in the oceanic plate and might constitute a valuable tool to monitor oceanic plate seismicity.
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.
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.
The Dead Sea region has faced substantial environmental challenges in recent decades, including water resource scarcity, similar to 1 m annual decreases in the water level, sinkhole development, ascending-brine freshwater pollution, and seismic disturbance risks. Natural processes are significantly affected by human interference as well as by climate change and tectonic developments over the long term. To get a deep understanding of processes and their interactions, innovative scientific approaches that integrate disciplinary research and education are required. The research project DESERVE (Helmholtz Virtual Institute Dead Sea Research Venue) addresses these challenges in an interdisciplinary approach that includes geophysics, hydrology, and meteorology. The project is implemented by a consortium of scientific institutions in neighboring countries of the Dead Sea (Israel, Jordan, Palestine Territories) and participating German Helmholtz Centres (KIT, GFZ, UFZ). A new monitoring network of meteorological, hydrological, and seismic/geodynamic stations has been established, and extensive field research and numerical simulations have been undertaken. For the first time, innovative measurement and modeling techniques have been applied to the extreme conditions of the Dead Sea and its surroundings. The preliminary results show the potential of these methods. First time ever performed eddy covariance measurements give insight into the governing factors of Dead Sea evaporation. High-resolution bathymetric investigations reveal a strong correlation between submarine springs and neo-tectonic patterns. Based on detailed studies of stratigraphy and borehole information, the extension of the subsurface drainage basin of the Dead Sea is now reliably estimated. Originality has been achieved in monitoring flash floods in an arid basin at its outlet and simultaneously in tributaries, supplemented by spatio-temporal rainfall data. Low-altitude, high resolution photogrammetry, allied to satellite image analysis and to geophysical surveys (e.g. shear-wave reflections) has enabled a more detailed characterization of sinkhole morphology and temporal development and the possible subsurface controls thereon. All the above listed efforts and scientific results take place with the interdisciplinary education of young scientists. They are invited to attend joint thematic workshops and winter schools as well as to participate in field experiments. (C) 2015 The Authors. Published by Elsevier B.V.
Earthquake source arrays
(2020)
A collection of earthquake sources recorded at a single station, under specific conditions, are considered as a source array (SA), that is interpreted as if earthquake sources originate at the station location and are recorded at the source location. Then, array processing methods, that is array beamforming, are applicable to analyse the recorded signals. A possible application is to use source array multiple event techniques to locate and characterize near-source scatterers and structural interfaces. In this work the aim is to facilitate the use of earthquake source arrays by presenting an automatic search algorithm to configure the source array elements. We developed a procedure to search for an optimal source array element distribution given an earthquake catalogue including accurate origin time and hypocentre locations. The objective function of the optimization process can be flexibly defined for each application to ensure the prerequisites (criteria) of making a source array. We formulated four quantitative criteria as subfunctions and used the weighted sum technique to combine them in one single scalar function. The criteria are: (1) to control the accuracy of the slowness vector estimation using the time domain beamforming method, (2) to measure the waveform coherency of the array elements, (3) to select events with lower location error and (4) to select traces with high energy of specific phases, that is, sp- or ps-phases. The proposed procedure is verified using synthetic data as well as real examples for the Vogtland region in Northwest Bohemia. We discussed the possible application of the optimized source arrays to identify the location of scatterers in the velocity model by presenting a synthetic test and an example using real waveforms.
The design of an array configuration is an important task in array seismology during experiment planning. Often the array response function (ARF), which depends on the relative position of array stations and frequency content of the incoming signals, is used as the array design criterion. In practice, additional constraints and parameters have to be taken into account, for example, land ownership, site-specific noise levels or characteristics of the seismic sources under investigation. In this study, a flexible array design framework is introduced that implements a customizable scenario modelling and optimization scheme by making use of synthetic seismograms. Using synthetic seismograms to evaluate array performance makes it possible to consider additional constraints. We suggest to use synthetic array beamforming as an array design criterion instead of the ARF. The objective function of the optimization scheme is defined according to the monitoring goals, and may consist of a number of subfunctions. The array design framework is exemplified by designing a seven-station small-scale array to monitor earthquake swarm activity in Northwest Bohemia/Vogtland in central Europe. Two subfunctions are introduced to verify the accuracy of horizontal slowness estimation; one to suppress aliasing effects due to possible secondary lobes of synthetic array beamforming calculated in horizontal slowness space and the other to reduce the event’s mislocation caused by miscalculation of the horizontal slowness vector. Subsequently, a weighting technique is applied to combine the subfunctions into one single scalar objective function to use in the optimization process.
Small-aperture array as a tool to monitor fluid injection- and extraction-induced microseismicity
(2019)
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.
Despite its high-seismogenic potential, the details of the seismogenic processes of Zagros Simply Folded Belt (SFB) remains debated. Three large earthquakes (M-w 7.3, 5.9 and 6.3) struck in the Lurestan arc of the Zagros SFB in 2017 and 2018. The sequence was recorded by seismic stations at regional, and teleseismic distances. Coseismic surface displacements, measured by Sentinel-1A/B satellites, provide additional data and a unique opportunity to study these earthquakes in detail. Here, we complement previous studies of the coseismic slip distribution of the 12 November 2017 M-w 7.3 Ezgeleh earthquake by a detailed analysis of its aftershocks, and we analysed the rupture process of the two interrelated earthquakes (25 August 2018 M-w 5.9 Tazehabad and the 25 November 2018 M-w 6.3 Sarpol-e Zahab earthquakes). We model the surface displacements obtained from Interferometric Synthetic Aperture Radar (InSAR) measurements and seismic records. We conduct non-linear probabilistic optimizations based on joint InSAR and seismic data to obtain finite-fault rupture of these earthquakes. The Lurestan arc earthquakes were followed by a sustained aftershock activity, with 133 aftershocks exceeding M-n 4.0 until 30 December 2019. We rely on the permanent seismic networks of Iran and Iraq to relocate similar to 700 M-n 3 + events and estimate moment tensor solutions for 85 aftershocks down to M-w 4.0. The 2017 Ezgeleh earthquake has been considered to activate a low-angle (similar to 17 degrees) dextral-thrust fault at the depth of 10-20 km. However, most of its aftershocks have shallow centroid depths (8-12 km). The joint interpretation of finite source models, moment tensor and hypocentral location indicate that the 2018 Tazehabad and Sarpol-e Zahab earthquakes ruptured different strike-slip structures, providing evidence for the activation of the sinistral and dextral strike-slip faults, respectively. The deformation in the Lurestan arc is seismically accommodated by a complex fault system involving both thrust and strike-slip faults. Knowledge about the deformation characteristics is important for the understanding of crustal shortening, faulting and hazard and risk assessment in this region.
Data recorded by distributed acoustic sensing (DAS) along an optical fibre sample the spatial and temporal properties of seismic wavefields at high spatial density. Often leading to massive amount of data when collected for seismic monitoring along many kilometre long cables. The spatially coherent signals from weak seismic arrivals within the data are often obscured by incoherent noise. We present a flexible and computationally efficient filtering technique, which makes use of the dense spatial and temporal sampling of the data and that can handle the large amount of data. The presented adaptive frequency-wavenumber filter suppresses the incoherent seismic noise while amplifying the coherent wavefield. We analyse the response of the filter in time and spectral domain, and we demonstrate its performance on a noisy data set that was recorded in a vertical borehole observatory showing active and passive seismic phase arrivals. Lastly, we present a performant open-source software implementation enabling real-time filtering of large DAS data sets.
The occurrence of deep low-frequency (DLF) microearthquakes beneath volcanoes is commonly attributed to mass transport in the volcanic plumbing system and used to infer feeding channels from and into magma reservoirs. The key question is how magmas migrate from depth to the shallow crust and whether magma reservoirs are currently being recharged. For the first time since the improvement of the local seismic networks in the East Eifel region (Rhineland-Palatinate, Germany), we detect and locate recurrent DLF earthquakes in the lower crust and upper mantle beneath the Laacher See Volcano (LSV), using a joint data set of permanent sensors and a temporary deployment. So far, eight DLF earthquake sequences were observed in four distinct clusters between 10 and 40 km depth. These clusters of weak events (M-L< 2) align along an approximately 80. southeast dipping line south of the LSV. Moment tensor solutions of these events have large shear components, and the irregular dispersion and long coda of body waves indicate interaction processes between shear cracks and fluids. We find a rotation of P-axes orientation for shallow tectonic earthquakes compared to DLF events, indicating that the stress field in the depth interval of DLF events might favour a vertical migration of magma or magmatic fluids. The caldera of the LSV was formed by the last major eruption of the East Eifel Volcanic Field only 12.9 kyr ago, fed by a shallow magma chamber at 5-8 km depth and erupting a total magma volume of 6.7 km(3). The observed DLF earthquake activity and continuous volcanic gas emissions around the LSV indicate an active magmatic system, possibly connected with an upper mantle melt zone.
The computation of such synthetic GFs is computationally and operationally demanding. As a consequence, the onthe-fly recalculation of synthetic GFs in each iteration of an optimisation is time-consuming and impractical. Therefore, the pre-calculation and efficient storage of synthetic GFs on a dense grid of source to receiver combinations enables the efficient lookup and utilisation of GFs in time-critical scenarios. We present a Python-based framework and toolkit - Pyrocko-GF - that enables the pre-calculation of synthetic GF stores, which are independent of their numerical calculation method and GF transfer function. The framework aids in the creation of such GF stores by interfacing a suite of established numerical forward modelling codes in seismology (computational back ends). So far, interfaces to back ends for layered Earth model cases have been provided; however, the architecture of Pyrocko-GF is designed to cover back ends for other geometries (e.g. full 3-D heterogeneous media) and other physical quantities (e.g. gravity, pressure, tilt). Therefore, Pyrocko-GF defines an extensible GF storage format suitable for a wide range of GF types, especially handling elasticity and wave propagation problems. The framework assists with visualisations, quality control, and the exchange of GF stores, which is supported through an online platform that provides many pre-calculated GF stores for local, regional, and global studies. The Pyrocko-GF toolkit comes with a well-documented application programming interface (API) for the Python programming language to efficiently facilitate forward modelling of geophysical processes, e.g. synthetic waveforms or static displacements for a wide range of source models.
Oceanic lithospheric S-wave velocities from the analysis of P-wave polarization at the ocean floor
(2016)
Our knowledge of the absolute S-wave velocities of the oceanic lithosphere is mainly based on global surface wave tomography, local active seismic or compliance measurements using oceanic infragravity waves. The results of tomography give a rather smooth picture of the actual S-wave velocity structure and local measurements have limitations regarding the range of elastic parameters or the geometry of the measurement. Here, we use the P-wave polarization (apparent P-wave incidence angle) of teleseismic events to investigate the S-wave velocity structure of the oceanic crust and the upper tens of kilometres of the mantle beneath single stations. In this study, we present an up to our knowledge new relation of the apparent P-wave incidence angle at the ocean bottom dependent on the half-space S-wave velocity. We analyse the angle in different period ranges at ocean bottom stations (OBSs) to derive apparent S-wave velocity profiles. These profiles are dependent on the S-wave velocity as well as on the thickness of the layers in the subsurface. Consequently, their interpretation results in a set of equally valid models. We analyse the apparent P-wave incidence angles of an OBS data set which was collected in the Eastern Mid Atlantic. We are able to determine reasonable S-wave-velocity-depth models by a three-step quantitative modelling after a manual data quality control, although layer resonance sometimes influences the estimated apparent S-wave velocities. The apparent S-wave velocity profiles are well explained by an oceanic PREM model in which the upper part is replaced by four layers consisting of a water column, a sediment, a crust and a layer representing the uppermost mantle. The obtained sediment has a thickness between 0.3 and 0.9 km with S-wave velocities between 0.7 and 1.4 km s(-1). The estimated total crustal thickness varies between 4 and 10 km with S-wave velocities between 3.5 and 4.3 km s(-1). We find a slight increase of the total crustal thickness from similar to 5 to similar to 8 km towards the South in the direction of a major plate boundary, the Gloria Fault. The observed crustal thickening can be related with the known dominant compression in the vicinity of the fault. Furthermore, the resulting mantle S-wave velocities decrease from values around 5.5 to 4.5 km s(-1) towards the fault. This decrease is probably caused by serpentinization and indicates that the oceanic transform fault affects a broad region in the uppermost mantle. Conclusively, the presented method is useful for the estimation of the local S-wave velocity structure beneath ocean bottom seismic stations. It is easy to implement and consists of two main steps: (1) measurement of apparent P-wave incidence angles in different period ranges for real and synthetic data, and (2) comparison of the determined apparent S-wave velocities for real and synthetic data to estimate S-wave velocity-depth models.
Receiver functions (RF) have been used for several decades to study structures beneath seismic stations. Although most available stations are deployed on shore, the number of ocean bottom station (OBS) experiments has increased in recent years. Almost all OBSs have to deal with higher noise levels and a limited deployment time (approximate to 1year), resulting in a small number of usable records of teleseismic earthquakes. Here we use OBSs deployed as midaperture array in the deep ocean (4.5-5.5km water depth) of the eastern mid-Atlantic. We use evaluation criteria for OBS data and beamforming to enhance the quality of the RFs. Although some stations show reverberations caused by sedimentary cover, we are able to identify the Moho signal, indicating a normal thickness (5-8km) of oceanic crust. Observations at single stations with thin sediments (300-400m) indicate that a probable sharp lithosphere-asthenosphere boundary (LAB) might exist at a depth of approximate to 70-80km which is in line with LAB depth estimates for similar lithospheric ages in the Pacific. The mantle discontinuities at approximate to 410km and approximate to 660km are clearly identifiable. Their delay times are in agreement with PREM. Overall the usage of beam-formed earthquake recordings for OBS RF analysis is an excellent way to increase the signal quality and the number of usable events.
Marine seismology usually relies on temporary deployments of stand alone seismic ocean bottom stations (OBS), which are initialized and synchronized on ship before deployment and re-synchronized and stopped on ship after recovery several months later. In between, the recorder clocks may drift and float at unknown rates. If the clock drifts are large or not linear and cannot be corrected for, seismological applications will be limited to methods not requiring precise common timing. Therefore, for example, array seismological methods, which need very accurate timing between individual stations, would not be applicable for such deployments.
We use an OBS test-array of 12 stations and 75 km aperture, deployed for 10 months in the deep sea (4.5-5.5 km) of the mid-eastern Atlantic. The experiment was designed to analyse the potential of broad-band array seismology at the seafloor. After recovery, we identified some stations which either show unusual large clock drifts and/or static time offsets by having a large difference between the internal clock and the GPS-signal (skew).
We test the approach of ambient noise cross-correlation to synchronize clocks of a deep water OBS array with km-scale interstation distances. We show that small drift rates and static time offsets can be resolved on vertical components with a standard technique. Larger clock drifts (several seconds per day) can only be accurately recovered if time windows of one input trace are shifted according to the expected drift between a station pair before the cross-correlation. We validate that the drifts extracted from the seismometer data are linear to first order. The same is valid for most of the hydrophones. Moreover, we were able to determine the clock drift at a station where no skew could be measured. Furthermore, we find that instable apparent drift rates at some hydrophones, which are uncorrelated to the seismometer drift recorded at the same digitizer, indicate a malfunction of the hydrophone.