Monitoring the 2020 Magna, Utah, Earthquake Sequence with Nodal Seismometers and Machine Learning

Author(s):  
Ben Baker ◽  
Monique M. Holt ◽  
Kristine L. Pankow ◽  
Keith D. Koper ◽  
Jamie Farrell

Abstract Immediately following the 18 March 2020 Mww 5.7 Magna, Utah, earthquake, work began on installing a network of three-component, 5 Hz geophones throughout the Salt Lake Valley. After six days, 180 geophones had been sited within 35 km of the epicenter. Each geophone recorded 250 samples per second data onsite for ∼40 days. Here, we integrate the geophone data with data from the permanent regional seismic network operated by the University of Utah Seismograph Stations (UUSS). We use machine learning (ML) methods to create a new catalog of arrival time picks, earthquake locations, and P-wave polarities for 18 March 2020–30 April 2020. We train two deep-learning U-Net models to detect P waves and S waves, assigning arrival times to maximal posterior probabilities, followed by a two-step association process that combines deep learning with a grid-based interferometric approach. Our automated workflow results in 142,000 P picks, 188,000 S picks, and over 5000 earthquake locations. We recovered 95% of the events in the UUSS authoritative catalog and more than doubled the total number of events (5000 vs. 2300). The P and S arrival times generated by our ML models have near-zero biases and standard deviations of 0.05 s and 0.09 s, respectively, relative to corresponding analyst times picked at backbone stations. We also use a deep-learning architecture to automatically determine 70,000 P-wave first motions, which agree with 93% of 5876 hand-picked up or down first motions from both the backbone and nodal stations. Overall, the use of ML led to large increases in the number of arrival times, especially S times, that will be useful for future tomographic studies, as well as the discovery of thousands more earthquakes than exist in the UUSS catalog.

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D283-D291 ◽  
Author(s):  
Peng Liu ◽  
Wenxiao Qiao ◽  
Xiaohua Che ◽  
Xiaodong Ju ◽  
Junqiang Lu ◽  
...  

We have developed a new 3D acoustic logging tool (3DAC). To examine the azimuthal resolution of 3DAC, we have evaluated a 3D finite-difference time-domain model to simulate a case in which the borehole penetrated a rock formation boundary when the tool worked at the azimuthal-transmitting-azimuthal-receiving mode. The results indicated that there were two types of P-waves with different slowness in waveforms: the P-wave of the harder rock (P1) and the P-wave of the softer rock (P2). The P1-wave can be observed in each azimuthal receiver, but the P2-wave appears only in the azimuthal receivers toward the softer rock. When these two types of rock are both fast formations, two types of S-waves also exist, and they have better azimuthal sensitivity compared with P-waves. The S-wave of the harder rock (S1) appears only in receivers toward the harder rock, and the S-wave of the softer rock (S2) appears only in receivers toward the softer rock. A model was simulated in which the boundary between shale and sand penetrated the borehole but not the borehole axis. The P-wave of shale and the S-wave of sand are azimuthally sensitive to the azimuth angle variation of two formations. In addition, waveforms obtained from 3DAC working at the monopole-transmitting-azimuthal-receiving mode indicate that the corresponding P-waves and S-waves are azimuthally sensitive, too. Finally, we have developed a field example of 3DAC to support our simulation results: The azimuthal variation of the P-wave slowness was observed and can thus be used to reflect the azimuthal heterogeneity of formations.


Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. T221-T231 ◽  
Author(s):  
Christine E. Krohn ◽  
Thomas J. Murray

The top 6 m of the near surface has a surprisingly large effect on the behavior of P- and S-waves. For unconsolidated sediments, the P-wave velocity gradient and attenuation can be quite large. Computer modeling should include these properties to accurately reproduce seismic effects of the near surface. We have used reverse VSP data and computer simulations to demonstrate the following effects for upgoing P-waves. Near the surface, we have observed a large time delay, indicating low velocity ([Formula: see text]), and considerable pulse broadening, indicating high attenuation ([Formula: see text]). Consequently, shallowly buried geophones have greater high-frequency bandwidth compared with surface geophones. In addition, there is a large velocity gradient in the shallow near surface (factor of 10 in 5 m), resulting in the rotation of P-waves to the vertical with progressively smaller amplitudes recorded on horizontal phones. Finally, we have found little indication of a reflection or ghost from the surface, although downgoing reflections have been observed from interfaces within the near surface. In comparison, the following have been observed for upgoing S-waves: There is a small increase in the time delay or pulse broadening near the surface, indicating a smaller velocity gradient and less change in attenuation. In addition, the surface reflection coefficient is nearly one with a prominent surface ghost.


2009 ◽  
Vol 24 (16n17) ◽  
pp. 3191-3225 ◽  
Author(s):  
KOJI HARADA ◽  
HIROFUMI KUBO ◽  
ATSUSHI NINOMIYA

We extend our Wilsonian renormalization group (RG) analysis on the pionless nuclear effective field theory in the two-nucleon sector in two ways; on the one hand, (1) we enlarge the space of operators up to including those of [Formula: see text] in the S waves, and, on the other hand, (2) we consider the RG flows in higher partial waves (P and D waves). In the larger space calculations, we find, in addition to nontrivial fixed points, two "fixed lines" and a "fixed surface" which are related to marginal operators. In the higher partial wave calculations, we find similar phase structures to that of the S waves, but there are two relevant directions in the P waves at the nontrivial fixed points and three in the D waves. We explain the physical meaning of the P-wave phase structure by explicitly calculating the low-energy scattering amplitude. We also discuss the relation between the Legendre flow equation which we employ and the RG equation by Birse, McGovern and Richardson, and possible implementation of power divergence subtraction in higher partial waves.


Geophysics ◽  
1992 ◽  
Vol 57 (11) ◽  
pp. 1444-1452 ◽  
Author(s):  
Guy W. Purnell

High‐velocity layers (HVLs) often hinder seismic imaging of deeper reflectors using conventional techniques. A major factor is often the unusual energy partitioning of waves incident at an HVL boundary from lower‐velocity material. Using elastic physical modeling, I demonstrate that one effect of this factor is to limit the range of dips beneath an HVL that can be imaged using unconverted P‐wave arrivals. At the same time, however, partitioning may also result in P‐waves outside the HVL coupling efficiently with S‐waves inside. By exploiting some of the waves that convert upon transmission into and/or out of the physical‐model HVL, I am able to image a much broader range of underlying dips. This is accomplished by acoustic migration tailored (via the migration velocities used) for selected families of converted‐wave arrivals.


1995 ◽  
Vol 85 (1) ◽  
pp. 254-268 ◽  
Author(s):  
Jie Zhang ◽  
Charles A. Langston

Abstract Teleseismic broadband P and S waves recorded at the NARS station NE06 (Dourbes, Belgium) are shown to exhibit strong anomalous particle motion not attributable to instrument miscalibration or malfunction. Azimuthally varying radial and tangential components have been observed on 38 recordings after vector rotation of horizontal P waves into the ray direction. The tangenital P waves attain amplitudes comparable to the radial components from the east with negative polarity and west with positive polarity, but tend to be zero in the north and south, suggesting major discontinuities in the crust dipping southward. The SH wave from the east contains a large SPmP phase, an S-to-P conversion at the free surface and then reflected back to the surface from the Moho. The polarity of this SPmP phase presents further evidence for a southward-dipping Moho. We employ ray theory for three-dimensionally dipping interfaces to compute the P-wave response. Linear inverse theory with smoothness constraints is applied to the simultaneous inversions of P-wave receiver functions for four different backazimuths. Through the progressive change of interface strike and dip and the inversion of layer shear-wave velocities, a dipping crustal model that is consistent with both the observed waveforms and results of previous local geophysical surveys has been determined. The results suggest a large velocity contrast in the shallow structure near the surface, another major interface at a depth of 12 km with dip of 10°, and a seismically transparent unit below the interface. The interface at a depth of 12 km reportedly emerges at the Midi fault 50 km north of the station NE06.


Author(s):  
Joel D. Simon ◽  
Frederik J. Simons ◽  
Jessica C. E. Irving

Abstract Mobile Earthquake Recorder in Marine Areas by Independent Divers (MERMAID) is a passively drifting oceanic diving float that transmits acoustic pressure records from global earthquakes within hours or days of their rupture. The onboard algorithm used for the detection and identification of signals from the hydrophone prioritizes the recovery of ∼1 Hz teleseismic P waves, which are useful for seismic imaging of Earth’s mantle. Two years into a mission that launched 50 MERMAIDs to map 3D mantle wavespeed anomalies with high resolution under the Pacific in French Polynesia, it is clear that the data returned contain much information beyond the first-arriving seismic P phases. These include acoustic conversions from S waves, surface waves, T waves, and inner- and outer-core phases, generated by earthquakes heard across the globe—and sounds from otherwise unidentified events occurring in remote and uninstrumented parts of the world’s oceans. Our growing database of automatically accumulating ∼240 s long-triggered segments contains a treasure trove for geophysicists interested in seismology beyond P-wave tomography. Furthermore, equipped with two-way communication capabilities, MERMAID can entertain requests to deliver data from its 1 yr buffer. In this article, we highlight the data classes and categories in MERMAID’s “extended-utility” catalog.


1974 ◽  
Vol 64 (5) ◽  
pp. 1501-1507 ◽  
Author(s):  
D. J. Sutton

Abstract A fall in P-wave velocity before the Gisborne earthquake of March 4, 1966 is indicated by arrival-time residuals of P waves from distant earthquakes recorded at the Gisborne seismograph station. Residuals were averaged over 6-month intervals from 1964 to 1968 and showed an increase of about 0.5 sec, implying later arrival times. The change began about 480 days before the earthquake. This precursory time interval is about that expected for an earthquake of this magnitude (ML = 6.2), but unlike most other reported instances, there was no obvious delay between the return of the velocity to normal and the occurrence of the earthquake. Similar analyses were carried out over the same period for two other New Zealand seismograph stations; at Karapiro there was no significant variation in mean residuals, and at Wellington the scatter was too large for the results to be meaningful. The Gisborne earthquake had a focus in the lower crust, about 25 km deep and was deeper than other events for which such precursory drops in P-wave velocity have been reported.


2013 ◽  
Vol 94 (1) ◽  
pp. 51-63 ◽  
Author(s):  
Neil P. Lareau ◽  
Erik Crosman ◽  
C. David Whiteman ◽  
John D. Horel ◽  
Sebastian W. Hoch ◽  
...  

The Persistent Cold-Air Pool Study (PCAPS) was conducted in Utah's Salt Lake valley from 1 December 2010 to 7 February 2011. The field campaign's primary goal was to improve understanding of the physical processes governing the evolution of multiday cold-air pools (CAPs) that are common in mountain basins during the winter. Meteorological instrumentation deployed throughout the Salt Lake valley provided observations of the processes contributing to the formation, maintenance, and destruction of 10 persistent CAP episodes. The close proximity of PCAPS field sites to residences and the University of Utah campus allowed many undergraduate and graduate students to participate in the study. Ongoing research, supported by the National Science Foundation, is using the PCAPS dataset to examine CAP evolution. Preliminary analyses reveal that variations in CAP thermodynamic structure are attributable to a multitude of physical processes affecting local static stability: for example, synoptic-scale processes impact changes in temperatures and cloudiness aloft while variations in boundary layer forcing modulate the lower levels of CAPs. During episodes of strong winds, complex interactions between the synoptic and mesoscale f lows, local thermodynamic structure, and terrain lead to both partial and complete removal of CAPs. In addition, the strength and duration of CAP events affect the local concentrations of pollutants such as PM2.5.


2018 ◽  
Vol 29 ◽  
pp. 00019
Author(s):  
Katarzyna Hubicka ◽  
Jakub Sokolowski

Seismic event consists of surface waves and body waves. Due to the fact that the body waves are faster (P-waves) and more energetic (S-waves) in literature the problem of their analysis is taken more often. The most universal information that is received from the recorded wave is its moment of arrival. When this information is obtained from at least four seismometers in different locations, the epicentre of the particular event can be estimated [1]. Since the recorded body waves may overlap in signal, the problem of wave onset moment is considered more often for faster P-wave than S-wave. This however does not mean that the issue of S-wave arrival time is not taken at all. As the process of manual picking is time-consuming, methods of automatic detection are recommended (these however may be less accurate). In this paper four recently developed methods estimating S-wave arrival are compared: the method operating on empirical mode decomposition and Teager-Kaiser operator [2], the modification of STA/LTA algorithm [3], the method using a nearest neighbour-based approach [4] and the algorithm operating on characteristic of signals’ second moments. The methods will be also compared to wellknown algorithm based on the autoregressive model [5]. The algorithms will be tested in terms of their S-wave arrival identification accuracy on real data originating from International Research Institutions for Seismology (IRIS) database.


1974 ◽  
Vol 64 (6) ◽  
pp. 1621-1627 ◽  
Author(s):  
J. C. Savage

abstract A comprehensive set of body-wave spectra has been calculated for the Haskell fault model generalized to a circular fault surface. These spectra are used to show that in practice the P-wave corner frequency (ƒp) may exceed the S-wave corner frequency (ƒs) when near-sonic or transonic rupture propagation obtains. The explanation appears to be that in such cases ƒs is so large that it is not identified within the recorded band, but rather a secondary corner is mistaken for ƒs. As a consequence of failing to detect the true asymptotic trend, the high-frequency falloff of the spectrum with frequency is substantially less for S waves than for P waves. This explanation appears to be consistent with the demonstration by Molnar, Tucker, and Brune (1973) that ƒp may exceed ƒs.


Sign in / Sign up

Export Citation Format

Share Document