In this study, I demonstrate that distributed acoustic sensing (DAS) raw strain rate data can directly be used to estimate spectral source parameters through an Empirical Green's Function (EGF) deconvolution analysis. Previously, DAS had been widely used in passive seismology to image the subsurface and analyze ground motion variations by converting strain or strain rate to particle velocity or acceleration prior to analysis. In this study, spectral analysis is applied to the PoroTomo joint DAS and seismic Nodal array in the Brady Hot Springs geothermal field to obtain source parameters for two M4 earthquakes via EGF analysis, where nearly collocated smaller events are used as an EGF to remove path and site effects. The EGF workflow is applied to raw DAS strain rate data without conversion to particle velocities and raw Nodal seismic data. The DAS and Nodal results are very consistent with similar features of spectral ratios, corner frequencies and moment ratios for the same event pairs. The uncertainty due to stacked spectral measurement is much lower on the DAS array, suggesting better stability of spectral shape measurement, possibly due to the much denser spatial sampling. The uncertainty due to model fitting is similar between DAS and Nodal arrays with slightly lower uncertainty on the DAS array. These observations demonstrate potential for directly using the strain rate measurements from DAS arrays for earthquake source characterizations.
Instrument Response Removal and the 2020 MLg 3.1 Marlboro, New Jersey, Earthquake
Abstract To better understand earthquakes as a hazard and to better understand the interior structure of the Earth, we often want to measure the physical displacement, velocity, or acceleration at locations on the Earth’s surface. To this end, a routine step in an observational seismology workflow is the removal of the instrument response, required to convert the digital counts recorded by a seismometer to physical displacement, velocity, or acceleration. The conceptual framework, which we briefly review for students and researchers of seismology, is that of the seismometer as a linear time-invariant system, which records a convolution of ground motion via a transfer function that gain scales and phase shifts the incoming signal. In practice, numerous software packages are widely used to undo this convolution via deconvolution of the instrument’s transfer function. Here, to allow the reader to understand this process, we start by taking a step back to fully explore the choices made during this routine step and the reasons for making them. In addition, we introduce open-source routines in Python and MATLAB as part of our rflexa package, which identically reproduce the results of the Seismic Analysis Code, a ubiquitous and trusted reference. The entire workflow is illustrated on data recorded by several instruments on Princeton University campus in Princeton, New Jersey, of the 9 September 2020 magnitude 3.1 earthquake in Marlboro, New Jersey.
more »
« less
- Award ID(s):
- 1736046
- NSF-PAR ID:
- 10294368
- Date Published:
- Journal Name:
- Seismological Research Letters
- ISSN:
- 0895-0695
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
SUMMARY -
Abstract Geological processes at subduction zones and their associated geohazards (e.g., megathrust earthquakes, submarine landslides, tsunamis, and arc volcanism) are, to a large extent, controlled by the structure, physical properties and fluid content of the subducting plate, the accreted sediments, and the overriding plate. In these settings, modern seismic modeling and imaging techniques based on controlled-source, multicomponent ocean-bottom seismometer (OBS) data are some of the best tools available for determining the subseafloor elastic properties, which can be linked to the aforementioned properties. Here, we present CASIE21-OBS, a controlled-source marine wide-angle OBS data set recently collected across the Cascadia convergent margin as part of the larger CAscadia Seismic Imaging Experiment 2021 (CASIE21). The main component of CASIE21 is a long-offset multichannel seismic (MCS) survey of the Cascadia margin conducted in June–July 2021 onboard R/V M.G. Langseth (cruise MGL2104) aiming to characterize the incoming plate, the plate interface geometry and properties, and the overlying sediment stratigraphy and physical properties. CASIE21-OBS was conducted during R/V M.G. Langseth cruise MGL2103 (May 2021) and R/V Oceanus cruise OC2106A (June–July 2021). It consisted of 63 short-period four-component OBSs deployed at a total 120 stations along 10 across-trench profiles extending from ∼50 km seaward of the deformation front to the continental shelf, and from offshore northern Vancouver Island to offshore southern Oregon. The OBSs recorded the airgun signals of the CASIE21-MCS survey as well as natural seismicity occurring during the deployment period (24 May 2021 19:00 UTC–9 July 2021 09:00 UTC). The OBS data are archived and available at the Incorporated Research Institutions for Seismology Data Management Center under network code YR_2021 for continuous time series (miniSEED) and identifier 21-008 for assembled data set (SEG-Y).more » « less
-
We document the performance of new ATONA (‘aA to nA’) amplifiers installed on an Isotopx Phoenix thermal ionisation mass spectrometer (TIMS) at Princeton University and evaluate their suitability for high-precision analyses of Pb and U isotopes in pg- to ng-size samples characteristic for U–Pb geochronology. The new amplifiers are characterised by low and stable noise levels comparable to 10 12 to 10 13 ohm resistors, response time <0.5 s, exceptional gain stability <1 ppm and a vast dynamic range theoretically allowing to quantify signals from aA (10 −18 A) to nA (10 −9 A) level. We measured a set of Pb standards, synthetic U–Pb solutions and natural zircons at currents of 2 × 10 −16 to 2 × 10 −12 A (corresponding to intensities of 20 μV to 200 mV relative to a 10 11 ohm amplifier) to assess the utility of ATONA in replacing ion counting for the smallest samples. The results show a clear precision benefit of using ATONA-Faraday detection over Daly ion counting for ion currents of >10 −14 A (1 mV relative to a 10 11 ohm amplifier or ca. 60 kcps). As such currents are routinely achievable for major Pb peaks of interest ( 205–208 Pb) in natural samples containing more than ca. 10 pg Pb* (radiogenic Pb), we expect ATONA-Faraday detection to find broad applications in U–Pb geochronology. Its practical use for low-blank, radiogenic samples continues to require ion counting for 204 Pb, either with a fixed Faraday–ion counter gain or using a dynamic two-step ( e.g. FaraDaly) method. Routine adoption of ATONA-Faraday collection in place of ion counting for most major Pb and U isotopes has the potential to increase sample throughput and precision, both improving the accessibility of isotope dilution (ID)-TIMS geochronology and pushing this technique towards better reproducibility.more » « less
-
null (Ed.)Abstract Local seismic events recorded by the large-N Incorporated Research Institutions for Seismology Community Wavefield Experiment in Oklahoma are used to estimate Moho reflections near the array. For events within 50 km of the center of the array, normal moveout corrections and receiver stacking are applied to identify the PmP and SmS Moho reflections on the vertical and transverse components. Corrections for the reported focal depths are applied to a uniform event depth. To stack signals from multiple events, further static corrections of the envelopes of the Moho reflected arrivals from the individual event stacks are applied. The multiple-event stacks are then used to estimate the pre-critical PmP and SmS arrivals, and an average Poisson’s ratio of 1.77±0.02 was found for the crust near the array. Using a modified Oklahoma Geological Survey (OGS) velocity model with this Poisson’s ratio, the time-to-depth converted PmP and SmS arrivals resulted in a Moho depth of 41±0.6 km. The modeling of wide-angle Moho reflections for selected events at epicenter-to-station distances of 90–135 km provides additional constraints, and assuming the modified OGS model, a Moho depth of 40±1 km was inferred. The difference between the pre-critical and wide-angle Moho estimates could result from some lateral variability between the array and the wide-angle events. However, both estimates are slightly shallower than the original OGS model Moho depth of 42 km, and this could also result from a somewhat faster lower crust. This study shows that local seismic events, including induced events, can be utilized to estimate properties and structure of the crust, which, in turn, can be used to better understand the tectonics of a given region. The recording of local seismicity on large-N arrays provides increased lateral phase coherence for the better identification of precritical and wide-angle reflected arrivals.more » « less
-
Yortsos, Yannis (Ed.)more » « less
Abstract Transfer learning (TL), which enables neural networks (NNs) to generalize out-of-distribution via targeted re-training, is becoming a powerful tool in scientific machine learning (ML) applications such as weather/climate prediction and turbulence modeling. Effective TL requires knowing (1) how to re-train NNs? and (2) what physics are learned during TL? Here, we present novel analyses and a framework addressing (1)–(2) for a broad range of multi-scale, nonlinear, dynamical systems. Our approach combines spectral (e.g. Fourier) analyses of such systems with spectral analyses of convolutional NNs, revealing physical connections between the systems and what the NN learns (a combination of low-, high-, band-pass filters and Gabor filters). Integrating these analyses, we introduce a general framework that identifies the best re-training procedure for a given problem based on physics and NN theory. As test case, we explain the physics of TL in subgrid-scale modeling of several setups of 2D turbulence. Furthermore, these analyses show that in these cases, the shallowest convolution layers are the best to re-train, which is consistent with our physics-guided framework but is against the common wisdom guiding TL in the ML literature. Our work provides a new avenue for optimal and explainable TL, and a step toward fully explainable NNs, for wide-ranging applications in science and engineering, such as climate change modeling.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1785/0220210118
