ABSTRACT Earthquakes generate infrasound in multiple ways. Acoustic coupling at the surface from vertical seismic velocity, termed local infrasound, is often recorded by infrasound sensors but has seen relatively little study. Over 140 infrasound stations have recently been deployed in Alaska. Most of these stations have single sensors, rather than arrays, and were originally installed as part of the EarthScope Transportable Array. The single sensor nature, paucity of ground-truth signals, and remoteness makes evaluating their data quality and utility challenging. In addition, despite notable recent advances, infrasound calibration and frequency response evaluation remains challenging, particularly for large networks and retrospective analysis of sensors already installed. Here, we examine local seismoacoustic coupling on colocated seismic and infrasound stations in Alaska. Numerous large earthquakes across the region in recent years generated considerable vertical seismic velocity and local infrasound that were recorded on colocated sensors. We build on previous work and evaluate the full infrasound station frequency response using seismoacoustic coupled waves. By employing targeted signal processing techniques, we show that a single seismometer may be sufficient for characterizing the response of an entire nearby infrasound array. We find that good low frequency (<1 Hz) infrasound station response estimates can be derived from large (Mw>7) earthquakes out to at least 1500 km. High infrasound noise levels at some stations and seismic-wave energy focused at low frequencies limit our response estimates. The response of multiple stations in Alaska is found to differ considerably from their metadata and are related to improper installation and erroneous metadata. Our method provides a robust way to remotely examine infrasound station frequency response and examine seismoacoustic coupling, which is being increasingly used in airborne infrasound observations, earthquake magnitude estimation, and other applications.
more »
« less
Comparing Sensitivities of Geodetic Processing Methods for Rapid Earthquake Magnitude Estimation
Abstract Rapid earthquake magnitude estimation from real-time space-based geodetic observation streams provides an opportunity to mitigate the impact of large and potentially damaging earthquakes by issuing low-latency warnings prior to any significant and destructive shaking. Geodetic contributions to earthquake characterization and rapid magnitude estimation have evolved in the last 20 yr, from post-processed seismic waveforms to, more recently, improved capacity of regional geodetic networks enabled real-time Global Navigation Satellite System seismology using precise point positioning (PPP) displacement estimates. In addition, empirical scaling laws relating earthquake magnitude to peak ground displacement (PGD) at a given hypocentral distance have proven effective in rapid earthquake magnitude estimation, with an emphasis on performance in earthquakes larger than ∼Mw 6.5 in which near-field seismometers generally saturate. Although the primary geodetic contributions to date in earthquake early warning have focused on the use of 3D position estimates and displacements, concurrent efforts in time-differenced carrier phase (TDCP)-derived velocity estimates also have demonstrated that this methodology has utility, including similarly derived empirical scaling relationships. This study builds upon previous efforts in quantifying the ambient noise of three-component ground-displacement and ground-velocity estimates. We relate these noise thresholds to expected signals based on published scaling laws. Finally, we compare the performance of PPP-derived PGD to TDCP-derived peak ground velocity (PGV), given several rich event datasets. Our results indicate that TDCP-PGV is more likely than PPP-PGD to detect intermediate magnitude (∼Mw 5.0–6.0) earthquakes, albeit with greater magnitude estimate uncertainty and across smaller epicentral distances. We conclude that the computationally lightweight TDCP-derived PGV magnitude estimation is complementary to PPP-derived PGD magnitude estimates, which could be produced at the network edge at high rates and with increased sensitivity to ground motion than current PPP estimates.
more »
« less
- Award ID(s):
- 1835791
- NSF-PAR ID:
- 10388975
- Date Published:
- Journal Name:
- Seismological Research Letters
- Volume:
- 93
- Issue:
- 3
- ISSN:
- 0895-0695
- Page Range / eLocation ID:
- 1497 to 1509
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Earthquake early warning (EEW) systems aim to forecast the shaking intensity rapidly after an earthquake occurs and send warnings to affected areas before the onset of strong shaking. The system relies on rapid and accurate estimation of earthquake source parameters. However, it is known that source estimation for large ruptures in real‐time is challenging, and it often leads to magnitude underestimation. In a previous study, we showed that machine learning, HR‐GNSS, and realistic rupture synthetics can be used to reliably predict earthquake magnitude. This model, called Machine‐Learning Assessed Rapid Geodetic Earthquake model (M‐LARGE), can rapidly forecast large earthquake magnitudes with an accuracy of 99%. Here, we expand M‐LARGE to predict centroid location and fault size, enabling the construction of the fault rupture extent for forecasting shaking intensity using existing ground motion models. We test our model in the Chilean Subduction Zone with thousands of simulated and five real large earthquakes. The result achieves an average warning time of 40.5 s for shaking intensity MMI4+, surpassing the 34 s obtained by a similar GNSS EEW model. Our approach addresses a critical gap in existing EEW systems for large earthquakes by demonstrating real‐time fault tracking feasibility without saturation issues. This capability leads to timely and accurate ground motion forecasts and can support other methods, enhancing the overall effectiveness of EEW systems. Additionally, the ability to predict source parameters for real Chilean earthquakes implies that synthetic data, governed by our understanding of earthquake scaling, is consistent with the actual rupture processes.
-
more » « less
Abstract Most large tsunamis are generated by earthquakes on offshore plate boundary megathrusts. The primary factors influencing tsunami excitation are the seismic moment, faulting geometry, and depth of the faulting. Efforts to provide rapid tsunami warning have emphasized seismic and geodetic methods for quickly determining the event size and faulting geometry. It remains difficult to evaluate the updip extent of rupture, which has significant impact on tsunami excitation. Teleseismic
P waves can constrain this issue; slip under deep water generates strongpwP water reverberations that persist as ringingP codaafter the directP phases from the faulting have arrived. Event‐averagedP coda/P amplitude measures at large epicentral distances (>80°), tuned to the dominant periods of deep waterpwP (~12–15 s), correlate well with independent models of whether slip extends to near the trench or not. Data at closer ranges (30° to 80°) reduce the time lag needed for inferring the updip extent of rupture to <15 min. Arrival ofPP andPPP phases contaminates closer distanceP codameasures, but this can be suppressed by azimuthal or distance binning of the measures. Narrowband spectral ratio measures and differential magnitude measures ofP codaand directP (m B ) perform comparably to broader band root‐mean‐square (RMS) measures.P coda/P levels for large nonmegathrust events are also documented. Rapid measurement ofP coda/P metrics after a large earthquake can supplement quick moment tensor determinations to enhance tsunami warnings; observation of largeP codalevels indicates that shallow submarine rupture occurred and larger than typical tsunami (for givenM W ) can be expected. -
null (Ed.)Abstract We measure pseudospectral and peak ground motions from 44 intermediate‐depth Mw≥4.9 earthquakes in the Cook Inlet region of southern Alaska, including those from the 2018 Mw 7.1 earthquake near Anchorage, to identify regional amplification features (0.1–5 s period). Ground‐motion residuals are computed with respect to an empirical ground‐motion model for intraslab subduction earthquakes, and we compute bias, between‐, and within‐event terms through a linear mixed‐effects regression. Between‐event residuals are analyzed to assess the relative source characteristics of the Cook Inlet earthquakes and suggest a difference in the scaling of the source with depth, relative to global observations. The within‐event residuals are analyzed to investigate regional amplification, and various spatial patterns manifest, including correlations of amplification with depth of the Cook Inlet basin and varying amplifications east and west of the center of the basin. Three earthquake clusters are analyzed separately and indicate spatial amplification patterns that depend on source location and exhibit variations in the depth scaling of long‐period basin amplification. The observations inform future seismic hazard modeling efforts in the Cook Inlet region. More broadly, they suggest a greater complexity of basin and regional amplification than is currently used in seismic hazard analyses.more » « less
-
SUMMARY The Ecuadorian forearc is a complex region of accreted terranes with a history of large megathrust earthquakes. Most recently, a Mw 7.8 megathrust earthquake ruptured the plate boundary offshore of Pedernales, Ecuador on 16 April 2016. Following this event, an international collaboration arranged by the Instituto Geofisico at the Escuela Politécnica Nacional mobilized a rapid deployment of 65 seismic instruments along the Ecuadorian forearc. We combine this new seismic data set with 14 permanent stations from the Ecuadorian national network to better understand how variations in crustal structure relate to regional seismic hazards along the margin. Here, we present receiver function adaptive common conversion point stacks and a shear velocity model derived from the joint inversion of receiver functions and surface wave dispersion data obtained through ambient noise cross-correlations for the upper 50 km of the forearc. Beneath the forearc crust, we observe an eastward dipping slow velocity anomaly we interpret as subducting oceanic crust, which shallows near the projected centre of the subducting Carnegie Ridge. We also observe a strong shallow positive conversion in the Ecuadorian forearc near the Borbon Basin indicating a major discontinuity at a depth of ∼7 km. This conversion is not ubiquitous and may be the top of the accreted terranes. We also observe significant north–south changes in shear wave velocity. The velocity changes indicate variations in the accreted terranes and may indicate an increased amount of hydration beneath the Manabí Basin. This change in structure also correlates geographically with the southern rupture limit of multiple high magnitude megathrust earthquakes. The earthquake record along the Ecuadorian trench shows that no event with a Mw >7.4 has ruptured south of ∼0.5°S in southern Ecuador or northern Peru. Our observations, along with previous studies, suggest that variations in the forearc crustal structure and subducting oceanic crust may influance the occurrence and spatial distribution of high magnitude seismicity in the region.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1785/0220210265
