An official website of the United States government
Abstract We develop an automated processing procedure to derive a new catalog of earthquake locations, magnitudes, and potencies and analyze 9 years of data between 2008 and 2016 in the San Jacinto fault‐zone region. Our procedure accounts for detailed 3‐D velocity structure using a probabilistic global‐search location inversion and obtains high‐precision relative event locations using differential travel times measured by cross‐correlating waveforms. The obtained catalog illuminates spatiotemporal seismicity patterns in the fault zone with observations for 108,800 earthquakes in the magnitude range −1.8 to 5.4. Inside a focus region consisting of an 80‐km by 50‐km rectangle oriented parallel to the main fault trace, we estimate a 99% detection rate of earthquakes with magnitude 0.6 and greater and detect and locate about 60% more events than those present in the Southern California Seismic Network catalog. The results provide the most complete catalog available for the focused study region during the analyzed period and include both deeper events and very shallow patches of seismicity not present in the regional catalog. The seismicity exhibits a variety of complex patterns that contain important information on deformation processes in the region. The fraction of event pairs with waveforms having cross‐correlation coefficients ≥0.95 is only about 3%, indicating diverse processes operating in the fault zone.
more »
« less
Discrimination between icequakes and earthquakes in southern Alaska: an exploration of waveform features using Random Forest algorithm
SUMMARY This study examines the feature space of seismic waveforms often used in machine learning applications for seismic event detection and classification problems. Our investigation centres on the southern Alaska region, where the seismic record captures diverse seismic activity, notably from the calving of marine-terminating glaciers and tectonic earthquakes along active plate boundaries. While the automated discrimination of earthquakes and glacier quakes is our nominal goal, this data set provides an outstanding opportunity to explore the general feature space of regional seismic phases. That objective has applicability beyond ice quakes and our geographic region of study. We make a noteworthy discovery that features rooted in the spectral content of seismic waveforms consistently outperform statistical and temporal features. Spectral features demonstrate robust performance, exhibiting resilience to class imbalance while being minimally impacted by factors such as epicentral distance and signal-to-noise ratio. We also conduct experiments on the transferability of the model and find that transferability primarily depends on the appearance of the waveforms. Finally, we analyse misclassified events and find examples that are identified incorrectly in the original regional catalogue.
more »
« less
- PAR ID:
- 10654846
- Publisher / Repository:
- Oxford Academic, Geophysical Journal International
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 237
- Issue:
- 2
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 1189 to 1207
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY We present a computational technique to model hydroacoustic waveforms from teleseismic earthquakes recorded by mid-column Mermaid floats deployed in the Pacific, taking into consideration bathymetric effects that modify seismo-acoustic conversions at the ocean bottom and acoustic wave propagation in the ocean layer, including reverberations. Our approach couples axisymmetric spectral-element simulations performed for moment-tensor earthquakes in a 1-D solid Earth to a 2-D Cartesian fluid–solid coupled spectral-element simulation that captures the conversion from displacement to acoustic pressure at an ocean-bottom interface with accurate bathymetry. We applied our workflow to 1129 seismograms for 682 earthquakes from 16 Mermaids (short for Mobile Earthquake Recording in Marine Areas by Independent Divers) owned by Princeton University that were deployed in the Southern Pacific as part of the South Pacific Plume Imaging and Modeling (SPPIM) project. We compare the modelled synthetic waveforms to the observed records in individually selected frequency bands aimed at reducing local noise levels while maximizing earthquake-generated signal content. The modelled waveforms match the observations very well, with a median correlation coefficient of 0.72, and some as high as 0.95. We compare our correlation-based traveltime measurements to measurements made on the same data set determined by automated arrival-time picking and ray- traced traveltime predictions, with the aim of opening up the use of Mermaid records for global seismic tomography via full-waveform inversion.more » « less
-
Abstract The San Fernando Valley (SFV), part of the Los Angeles metropolitan area, is a seismically active urban environment. Large-magnitude earthquakes, such as the 1994 Mw 6.7 Northridge event that occurred on a blind fault beneath the valley, caused significant infrastructure damage in the region, underscoring the need for enhanced seismic monitoring to improve the identification of buried faults and hazard evaluation. Currently, the Southern California Earthquake Data Center operates four broadband instruments within the valley; however, the network’s ability to capture small earthquakes beneath the region may be limited. To demonstrate how this data gap can be filled, we use recordings from the SFV array, comprised of 140 nodal instruments with interstation distances ranging from 0.3 to 2.5 km that recorded for one month. High-anthropogenic noise levels in urbanized areas tend to conceal earthquake signals; therefore, we applied a previously developed machine learning model fine-tuned on similar waveforms to detect events and pick seismic phases. In a two-step event association workflow, isolated phase picks were first culled, which eliminated false positive detections and reduced computational runtime. We located 62 events within a 209 km radius of our array with magnitudes ranging from ML 0.13 to 4, including 36 new events that were undetected by the regional network. One event cluster reveals a previously unidentified (5.3 km × 4 km) blind fault zone located ∼5 km beneath the southern part of the valley. Seismicity from this zone is rare in the regional catalog (<3 events per year), despite producing a Mb 4.4 event in 2014. Our results highlight the benefits of detecting small-magnitude seismicity for hazard estimation. Temporary nodal arrays can identify critical gaps in regional monitoring and guide site selection for permanent stations. In addition, our workflow can be applied to complement seismic monitoring in other urban settings.more » « less
-
Abstract The spectra of earthquake waveforms can provide important insight into rupture processes, but the analysis and interpretation of these spectra is rarely straightforward. Here we develop a Bayesian framework that embraces the inherent data and modeling uncertainties of spectral analysis to infer key source properties. The method uses a spectral ratio approach to correct the observedS‐wave spectra of nearby earthquakes for path and site attenuation. The objective then is to solve for a joint posterior probability distribution of three source parameters—seismic moment, corner frequency, and high‐frequency falloff rate—for each earthquake in the sequence, as well as a measure of rupture directivity for select target events with good azimuthal station coverage. While computationally intensive, this technique provides a quantitative understanding of parameter tradeoffs and uncertainties and allows one to impose physical constraints through prior distributions on all source parameters, which guide the inversion when data is limited. We demonstrate the method by analyzing in detail the source properties of 14 different target events of magnitude M5 in southern California that span a wide range of tectonic regimes and fault systems. These prominent earthquakes, while comparable in size, exhibit marked diversity in their source properties and directivity, with clear spatial patterns, depth‐dependent trends, and a preference for unilateral directivity. These coherent spatial variations source properties suggest that regional differences in tectonic setting, hypocentral depth or fault zone characteristics may drive variability in rupture processes, with important implications for our understanding of earthquake physics and its relation to hazard.more » « less
-
SUMMARY Seismic tomography is a principal method for studying mantle structure, but imaging of Earth’s wave speed anomalies is conditioned by seismic wave sampling. Global models use misfit criteria that may strive for balance between portions of the data set but can leave important regional domains underserved. We evaluate two full-waveform global tomography wave speed models, GLAD-M25 (Global adjoint tomography model) and SEMUCB-WM1 (whole-mantle tomography model derived from fully numerical spectral element method forward modelling), in the mantle below the Pacific Ocean. The region of the South Pacific Superswell contains multiple hotspots which may be fed by plumes anchored in the Large Low Shear-Velocity Province at the base of the mantle. The uneven distribution of seismic receivers worldwide leaves several candidate plumes beneath various hotspots poorly resolved. We assess the regional quality of GLAD-M25 relative to its global performance using a partition of the seismic waveform data used in its construction. We evaluate synthetic waveforms computed using the spectral-element method to determine how well they fit the data according to a variety of criteria measured across multiple seismic phases and frequency bands. The distributions of traveltime anomalies that remain in GLAD-M25 are wider for trans-Pacific paths than globally, suggesting comparatively insufficiently resolved seismic velocity structure in the region of interest. Hence, Pacific-centred regional inversions, based on (augmented) subsets of the global data set have the potential to enhance the resolution of velocity structure. We compare GLAD-M25 and SEMUCB-WM1 by cross-validation with a new, independent, data set. Our results reveal that short- and long-wavelength structure is captured differently by the two models. Our findings lead us to recommend focusing future model iteration on and around the Pacific Superswell and adding data that sample new corridors, especially using ocean sensors, to better constrain seismic velocity structure in this area of significant geodynamic complexity.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/gji/ggae106
