An official website of the United States government
Seismic array processing is routinely used to infer detailed earthquakeproperties of intermediate and large events, however, the sourceproperties of microseismicity often remain elusive. In this study, weuse high signal-to-noise ratio seismograms of 204 earthquakes induced bythe 6 km deep 2018 Espoo/Helsinki geothermal stimulation to evaluate thecapabilities of beamforming and back-projection array methods. We showthat mini array beamforming is sensitive to medium heterogeneity andrequires calibration to mitigate systematic slowness biases.A combinedand wave back-projection approach significantly improves depthresolution, reducing offsets to catalogue locations from km to m.Supported by numerical experiments, we demonstrate that back-projectionswimming patterns can constrain focal mechanisms. Our results imply thatback-projection of data collected over a wide azimuthal range can beused to monitor and characterize local-scale microseismicity, whereasbeamforming calibration requires independently obtained referenceobservations.
more »
« less
Source Properties of the Induced ML 0.0–1.8 Earthquakes from Local Beamforming and Backprojection in the Helsinki Area, Southern Finland
Abstract Seismic arrays constrain local wave propagation that can be used to infer earthquake source characteristics. Array processing is routinely used to infer detailed earthquake properties of intermediate and large events. However, the source properties of microseismicity often remain elusive. In this study, we use high signal-to-noise ratio seismograms of 204 ML 0.0–1.8 earthquakes induced by the 6 km deep 2018 Espoo/Helsinki geothermal stimulation to evaluate the performance and capabilities of beamforming and backprojection array methods. Using accurate travel-time-based event locations as a reference, we first show that miniarray beamforming is sensitive to medium heterogeneities and requires calibration to mitigate local systematic slowness biases. A catalog-based calibration significantly improves our multiarray beam raytracing estimates of source locations. Second, the application of the backprojection technique using P-wave signals with sufficient azimuthal coverage yields hypocenter estimates with generally good horizontal but poor vertical resolution. The short local source–receiver distances result in incomplete separation of P- and S-wave arrivals during backprojection. Numerical tests show that the relatively large S-wave amplitudes can influence coherent P-wave stacks, resulting in large location errors. Our combined P- and S-wave backprojection approach mitigates the influence of the large S-wave amplitude and improves the depth resolution significantly. The average depth offset to the reference catalog locations reduces from ≥1.4 km to ∼91 m. Third, 3D numerical simulations demonstrate that backprojection swimming patterns are not merely processing or configuration artifacts. We show that the swimming patterns correlate with and can resolve the source focal mechanism when the azimuthal wavefield sampling is sufficiently complete. Our work demonstrates that the backprojection techniques can help to better constrain important properties of local-scale microseismicity.
more »
« less
- Award ID(s):
- 2121568
- PAR ID:
- 10633055
- Publisher / Repository:
- SSA
- Date Published:
- Journal Name:
- Seismological Research Letters
- Volume:
- 96
- Issue:
- 1
- ISSN:
- 0895-0695
- Page Range / eLocation ID:
- 111 to 129
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We introduce BPMF (backprojection and matched filtering)—a complete and fully automated workflow designed for earthquake detection and location, and distributed in a Python package. This workflow enables the creation of comprehensive earthquake catalogs with low magnitudes of completeness using no or little prior knowledge of the study region. BPMF uses the seismic wavefield backprojection method to construct an initial earthquake catalog that is then densified with matched filtering. BPMF integrates recent machine learning tools to complement physics-based techniques, and improve the detection and location of earthquakes. In particular, BPMF offers a flexible framework in which machine learning detectors and backprojection can be harmoniously combined, effectively transforming single-station detectors into multistation detectors. The modularity of BPMF grants users the ability to control the contribution of machine learning tools within the workflow. The computation-intensive tasks (backprojection and matched filtering) are executed with C and CUDA-C routines wrapped in Python code. This leveraging of low-level, fast programming languages and graphic processing unit acceleration enables BPMF to efficiently handle large datasets. Here, we first summarize the methodology and describe the application programming interface. We then illustrate BPMF’s capabilities to characterize microseismicity with a 10 yr long application in the Ridgecrest, California area. Finally, we discuss the workflow’s runtime scaling with numerical resources and its versatility across various tectonic environments and different problems.more » « less
-
Abstract Earthquake focal mechanisms provide critical in-situ insights about the subsurface faulting geometry and stress state. For frequent small earthquakes (magnitude< 3.5), their focal mechanisms are routinely determined using first-arrival polarities picked on the vertical component of seismometers. Nevertheless, their quality is usually limited by the azimuthal coverage of the local seismic network. The emerging distributed acoustic sensing (DAS) technology, which can convert pre-existing telecommunication cables into arrays of strain/strain-rate meters, can potentially fill the azimuthal gap and enhance constraints on the nodal plane orientation through its long sensing range and dense spatial sampling. However, determining first-arrival polarities on DAS is challenging due to its single-component sensing and low signal-to-noise ratio for direct body waves. Here, we present a data-driven method that measures P-wave polarities on a DAS array based on cross-correlations between earthquake pairs. We validate the inferred polarities using the regional network catalog on two DAS arrays, deployed in California and each comprising ~ 5000 channels. We demonstrate that a joint focal mechanism inversion combining conventional and DAS polarity picks improves the accuracy and reduces the uncertainty in the focal plane orientation. Our results highlight the significant potential of integrating DAS with conventional networks for investigating high-resolution earthquake source mechanisms.more » « less
-
Abstract Accurate and (near) real-time earthquake monitoring provides the spatial and temporal behaviors of earthquakes for understanding the nature of earthquakes, and also helps in regional seismic hazard assessments and mitigations. Because of the increase in both the quality and quantity of seismic data, an automated earthquake monitoring system is needed. Most of the traditional methods for detecting earthquake signals and picking phases are based on analyses of features in recordings of an individual earthquake and/or their differences from background noises. When seismicity is high, the seismograms are complicated, and, therefore, traditional analysis methods often fail. With the development of machine learning algorithms, earthquake signal detection and seismic phase picking can be more accurate using the features obtained from a large amount of earthquake recordings. We have developed an attention recurrent residual U-Net algorithm, and used data augmentation techniques to improve the accuracy of earthquake detection and seismic phase picking on complex seismograms that record multiple earthquakes. The use of probability functions of P and S arrivals and potential P and S arrival pairs of earthquakes can increase the computational efficiency and accuracy of backprojection for earthquake monitoring in large areas. We applied our workflow to monitor the earthquake activity in southern California during the 2019 Ridgecrest sequence. The distribution of earthquakes determined by our method is consistent with that in the Southern California Earthquake Data Center (SCEDC) catalog. In addition, the number of earthquakes in our catalog is more than three times that of the SCEDC catalog. Our method identifies additional earthquakes that are close in origin times and/or locations, and are not included in the SCEDC catalog. Our algorithm avoids misidentification of seismic phases for earthquake location. In general, our algorithm can provide reliable earthquake monitoring on a large area, even during a high seismicity period.more » « less
-
The 2019 Ridgecrest, CA earthquake sequence has provided a unique opportunity and a rich dataset to understand earthquake source properties and near-fault structure. Using the high-quality seismic data provided by the SCEC Stress Drop Validation group, we first estimate the corner frequency of M2.0-4.5 earthquakes by applying the spectral ratio method based on empirical Green’s function (Liu et al., 2020). We relate corner frequency estimates to stress drops assuming the Brune source model and circular cracks. Our preliminary results show increasing median stress drops with magnitude for both P and S waves, from 1 MPa for M2.0 events to 10 MPa for M4.0 events, though the limited frequency bandwidth may cause underestimation for small events. The estimated moment magnitude is proportional to the catalog magnitude by a factor of 0.72, which is close to 0.74 estimated by Trugman (2020) for the Ridgecrest earthquake sequence. In the second part of the study, we examine the impact of fault zone structure on the azimuthal variation of the source spectra. Using kinematic simulations and observations of the 2003 Big Bear earthquake sequence, Huang et al. (2016) showed that fault damage zones can act as an effective wave guide and cause high-frequency wave amplification along directions close to fault strike. We use clusters of M1.5-3 earthquakes in the Ridgecrest region to further examine the azimuthal variation of the stacked source spectra and investigate if the near-source structure can affect our corner frequency estimates. We aim to develop robust methods that utilize high-quality seismic data to illuminate earthquake source processes and fault zone properties.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1785/0220240122
