Geology topics

David R. Shelly

Publications and source records attributed to David R. Shelly.

At least 55 records · Page 3Linked to original sources

USGS2020

Leveraging deep learning in global 24/7 real-time earthquake monitoring at the National Earthquake Information Center

Machine‐learning algorithms continue to show promise in their application to seismic processing. The U.S. Geological Survey National Earthquake Information Center (NEIC) is exploring the adoption of these tools to aid in simultaneous local, regional, and global real‐time earthquake monitoring. As a first step, we describe a simple framework to incorporate deep‐learning tools into NEIC operations. Automatic seismic arrival detections made from standard picking methods (e.g., short‐term average/long‐term average [STA/LTA]) are fed to trained neural network models to improve automatic seismic‐arrival (pick) timing and estimate seismic‐arrival phase type and source‐station distances. These additional data are used to improve the capabilities of the NEIC associator. We compile a dataset of 1.3 million seismic‐phase arrivals that represent a globally distributed set of source‐station paths covering a range of phase types, magnitudes, and source distances. We train three separate convolutional neural network models to predict arrival time onset, phase type, and distance. We validate the performance of the trained networks on a subset of our existing dataset and further extend validation by exploring the model performance when applied to NEIC automatic pick data feeds. We show that the information provided by these models can be useful in downstream event processing, specifically in seismic‐phase association, resulting in reduced false associations and improved location estimates.

Seismological Research Letters↗

USGS2020

Geodetic measurements of slow slip events southeast of Parkﬁeld, CA

Tremor and low-frequency earthquakes are presumed to be indicative of surrounding slow, aseismic slip that is often below geodetic detection thresholds. This study uses data from borehole seismometers and long-baseline laser strainmeters to observe both the seismic and geodetic signatures of episodic tremor and slip on the Parkfield region of the San Andreas Fault near Cholame, CA. The observed occurrence rates of both the tremors and co-located families of low-frequency earthquakes are not steady but instead exhibit quasiperiodic bursts of increased activity. We show that these periods of elevated seismic activity correlate with statistically significant stacked strain signals consisting of 44 slow-slip events. Modeled individual slow-slip events and their total summed moment, which are constrained by seismic signals and stacked strain, respectively, indicate that the individual moment magnitudes of these events range from 4.6–5.2. We find that the measured geodetic signal likely precedes the seismic signal by several hours, consistent with the aseismic slip preceding and driving the observed seismic tremor activity. We confirm that strike-slip faults, in addition to subduction zones, are capable of producing episodic tremor and slip.

California↗

USGS2020

A high-resolution seismic catalog for the initial 2019 Ridgecrest Earthquake sequence: Foreshocks, aftershocks, and faulting complexity

I use template matching and precise relative relocation techniques to develop a high-resolution earthquake catalog for the initial portion of the 2019 Ridgecrest earthquake sequence, from July 4-16, encompassing the foreshock sequence and the first 10+ days of aftershocks following the M w 7.1 mainshock. Using 13,525 routinely cataloged events as waveform templates, I detect and precisely locate a total of 34,091 events. Precisely located earthquakes reveal numerous cross-cutting fault structures with dominantly perpendicular SW- and NW-strikes. Foreshocks of the M w 6.4 event appear to align on a NW-striking fault. Aftershocks of the M w 6.4 event suggest that it further ruptured this NW-striking fault, as well as the SW-striking fault where surface rupture was observed. Finally, aftershocks of the M w 7.1 show a highly complex distribution, illuminating a primary NW-striking fault zone consistent with surface rupture, but also numerous cross-cutting SW-striking faults. Aftershock relocations suggest that the M w 7.1 event ruptured adjacent to the previous NW-striking rupture of the M w 6.4, perhaps activating a subparallel structure southwest of the earlier rupture. Both the NW and SE rupture termini of the M w 7.1 rupture exhibited multiple-fault branching, with particularly high rates of aftershocks and multiple fault orientations in the dilatational quadrant NE of the NW rupture terminus.

California↗

USGS2020

Post-1978 tumescence at Long Valley Caldera, California: A geophysical perspective

Long Valley Caldera has been restless since at least 1978. Prominent symptoms of this unrest include earthquake swarms and tumescence (inflation) centered on the resurgent dome. Over the years, interpretations of physical processes underlying this unrest have varied considerably. Results from a collection of geophysical studies infer the presence and/or active intrusion of magma in the crust. Geologic evidence, however, does not support recent magmatic activity in the caldera, leading to an interpretation that the caldera volcanic system is moribund, and the current unrest is a result of second boiling (aqueous fluids released during crystallization of the rhyolitic magma that produced the Bishop Tuff). Here, we examine the collective constraints provided by geophysical studies over the past four decades. Although the current geophysical evidence does not conclusively discriminate between unrest driven by recent crustal magmatic intrusion versus second boiling, it does provide evidence that a large volume of partial melt persists within the mid- and lower crust. This implies that the Long Valley Caldera system as a whole is long-lived, and the magma reservoir remains at least partially molten. In the shallow crust, the possibility of small pockets of magma remains. In aggregate, the data suggest commonality with other large, long-lived silicic caldera systems, such as Yellowstone and Campi Flegrei. Although the possibility of eruption within Long Valley Caldera remains unlikely, the geophysical evidence argues that Long Valley has an active magmatic system at depth and we must retain the possibility of eruptive hazards.

California↗

USGS2019

Anatomy of a caldera collapse: Kīlauea 2018 summit seismicity sequence in high resolution

The 2018 Kīlauea eruption and caldera collapse generated intense cycles of seismicity tied to repeated large seismic (Mw ~5) collapse events associated with magma withdrawal from beneath the summit. To gain insight into the underlying dynamics and aid eruption response, we applied waveform-based earthquake detection and double-difference location as the eruption unfolded. Here, we augment these rapid results by grouping events based on patterns of correlation-derived phase polarities across the network. From April 29 to August 6, bracketing the eruption, we used ~2800 events cataloged by the Hawaiian Volcano Observatory to detect and precisely locate 44,000+ earthquakes. Resulting hypocentroids resolve complex, yet coherent structures, concentrated at shallow depths east of Halema‘uma‘u crater, beneath the eventual eastern perimeter of surface collapse. Based on a preponderance of dilatational P-wave first motions and similarities with previously inferred dike structures, we hypothesize that failure was dominated by coupled shear and crack closure.

Hawaii↗

USGS2019

Snowmelt-triggered earthquake swarms at the margin of Long Valley Caldera, California

Fluids are well known to influence earthquakes, yet rarely are earthquakes convincingly linked to precipitation. Weak modulation or limited data often leads to ambiguous interpretations. In contrast, here we find that shallow seismicity in the Sierra Nevada range near Long Valley Caldera is strongly modulated by snowmelt. Over 33 years, shallow seismicity rates were ~37 times higher during very wet periods versus very dry periods. Relative earthquake relocations from a swarm in 2017 reveal downward migration from ~1- to 3-km depth along a steeply inclined plane. Steeply dipping strata may provide high-permeability pathways and faulting plane. Here we combine the correlated seismicity and hydrologic time series with the propagation observed in the relatively relocated earthquakes. From this combined evidence, we infer that pressure diffusion from groundwater recharge dramatically accelerated shallow seismicity rates, causing seismic swarms unrelated to volcanic processes.

California↗

USGS2018

Seismic evidence for significant melt beneath the Long Valley Caldera, California, USA

A little more than 760 ka ago, a supervolcano on the eastern edge of California (United States) underwent one of North America's largest Quaternary explosive eruptions. Over this ~6-day-long eruption, pyroclastic flows blanketed the surrounding ~50 km with more than 1400 km 3 of the now-iconic Bishop Tuff, with ashfall reaching as far east as Nebraska. Collapse of the volcano's magma reservoir created the restless Long Valley Caldera. Although no rhyolitic eruptions have occurred in 100 k.y., beginning in 1978, ongoing uplift suggests new magma may have intruded into the reservoir. Alternatively, the reservoir could be approaching final crystallization, with present-day uplift related to the expulsion of fluid from the last vestiges of melt. Despite 40 years of diverse investigations, the presence of large volumes of melt in Long Valley's magma reservoir remain unresolved. Here we show, through full waveform seismic tomography, a mid-crustal zone of low shear-wave velocity. We estimate the reservoir contains considerable quantities of melt, >1000 km 3 , at melt fractions as high as ~27%. While supervolcanoes like Long Valley are rare, understanding the volume and concentration of melt in their magma reservoirs is critical for determining their potential hazard.

California↗

USGS2018

Imaging a crustal low-velocity layer using reflected seismic waves from the 2014 earthquake swarm at Long Valley Caldera, California: The magmatic system roof?

The waveforms generated by the 2014 Long Valley Caldera earthquake swarm recorded at station MLH show clear reflected waves that are often stronger than direct P and S waves. With waveform analyses, we discover that these waves are reflected at the top of a low-velocity body, which may be residual magma from the ∼767 ka caldera-forming eruption. The polarity of the reflection compared to direct P and S waves suggests that the reflection is S P waves ( S from hypocenters to reflector and then convert to P waves to the surface). Because the wavefields are coherent among different earthquakes and hold high signal-to-noise ratios, we apply them to a wavefield migration method for imaging reflectors. The depth of the imaged magmatic system roof is around 8.2 km below the surface. This is consistent with previous studies. Even though we use only one station and waveforms from one earthquake swarm, the dense cluster of accurately located earthquakes provides a high-resolution image of the roof.

California↗

USGS2018

Constraints on friction, dilatancy, diffusivity, and effective stress from low-frequency earthquake rates on the deep San Andreas Fault

Families of recurring low-frequency earthquakes (LFEs) within nonvolcanic tremor on the San Andreas Fault in central California are sensitive to tidal stresses. LFEs occur at all levels of the tides, are strongly correlated and in phase with the ~200 Pa shear stresses, and weakly and not systematically correlated with the ~2 kPa tidal normal stresses. We assume that LFEs are small sources that repeatedly fail during shear within a much larger scale, aseismically slipping fault zone and consider two different models of the fault slip: (1) modulation of the fault slip rate by the tidal stresses or (2) episodic slip, triggered by the tides. LFEs are strongly clustered with duration much shorter than the semidiurnal tide; they cannot be significantly modulated on that time scale. The recurrence times of clusters, however, are many times longer than the semidiurnal, leading to an appearance of tidal triggering. In this context we examine the predictions of laboratory-observed triggered frictional (dilatant) fault slip. The undrained end-member model produces no sensitivity to the tidal normal stress, and slip onsets are in phase with the tidal shear stress. The tidal correlation constrains the diffusivity to be less than ~1 × 10 −6 /s and the product of the friction and dilatancy coefficients to be at most 5 × 10 −7 , orders of magnitude smaller than observed at room temperature. In the absence of dilatancy the effective normal stress at failure would be about ~55 kPa. For this model the observations require intrinsic weakness, low dilatancy, and lithostatic pore fluid.

Journal of Geophysical Research↗

USGS2018

Kinematics of the 2015 San Ramon, California earthquake swarm: Implications for fault zone structure and driving mechanisms

Earthquake swarms represent a sudden increase in seismicity that may indicate a heterogeneous fault-zone, the involvement of crustal fluids and/or slow fault slip. Swarms sometimes precede major earthquake ruptures. An earthquake swarm occurred in October 2015 near San Ramon, California in an extensional right step-over region between the northern Calaveras Fault and the Concord–Mt. Diablo fault zone, which has hosted ten major swarms since 1970. The 2015 San Ramon swarm is examined here from 11 October through 18 November using template matching analysis. The relocated seismicity catalog contains ∼4000 events with magnitudes between − 0.2 < M d < 3.6 "> −0.2< M d <3.6 . The swarm illuminated three sub-parallel, southwest striking and northwest dipping fault segments of km-scale dimension and thickness of up to 200 m. The segments contain coexisting populations of different focal-mechanisms, suggesting a complex fault zone structure with several sets of en échelon fault orientations. The migration of events along the three planar structures indicates a complex fluid and faulting interaction processes. We searched for correlations between seismic activity and tidal stresses and found some suggestive features, but nothing that we can be confident is statistically significant.

California↗

USGS2018

Deep fluid pathways beneath Mammoth Mountain, California, illuminated by migrating earthquake swarms

Although most volcanic seismicity is shallow (within several kilometers of the surface), some volcanoes exhibit deeper seismicity (10 to 30+ km) that may reflect active processes such as magma resupply and volatile transfer. One such volcano is Mammoth Mountain, California, which has also recently exhibited high rates of CO 2 discharge at the surface. We perform high-resolution earthquake detection and relocation to reveal punctuated episodes of rapidly propagating seismicity at mid-crustal depths along a narrow fracture zone surrounding a body of partial melt. We infer that these earthquakes track dike intrusions or fluid pressure pulses associated with CO 2 exsolution, suggesting that the deep plumbing system of Mammoth Mountain is an active conduit for fluid transport from the base of the crust to the surface.

California↗

USGS2017

A 15-year catalog of more than 1 million low-frequency earthquakes: Tracking tremor and slip along the deep San Andreas Fault

Low-frequency earthquakes (LFEs) are small, rapidly recurring slip events that occur on the deep extensions of some major faults. Their collective activation is often observed as a semi-continuous signal known as tectonic (or non-volcanic) tremor. This manuscript presents a catalog of more than 1 million LFEs detected along the central San Andreas Fault from 2001-2016. These events have been detected via a multi-channel matched-filter search, cross-correlating waveform templates representing 88 different LFE families with continuous seismic data. Together, these source locations span nearly 150 km along the central San Andreas Fault, ranging in depth from ~16-30 km. This accumulating catalog has been the source of numerous studies examining the behavior of these LFE sources and the inferred slip behavior of the deep fault. The relatively high temporal and spatial resolution of the catalog has provided new insights into properties such as tremor migration, recurrence, and triggering by static and dynamic stress perturbations. Collectively, these characteristics are inferred to reflect a very weak fault likely under near-lithostatic fluid pressure, yet the physical processes controlling the stuttering rupture observed as tremor and LFE signals remain poorly understood. This paper aims to document the LFE catalog assembly process and associated caveats, while also updating earlier observations and inferred physical constraints. The catalog itself accompanies this manuscript as part of the electronic supplement, with the goal of providing a useful resource for continued future investigations.

California↗

USGS2016

A new strategy for earthquake focal mechanisms using waveform-correlation-derived relative polarities and cluster analysis: Application to the 2014 Long Valley Caldera earthquake swarm

In microseismicity analyses, reliable focal mechanisms can typically be obtained for only a small subset of located events. We address this limitation here, presenting a framework for determining robust focal mechanisms for entire populations of very small events. To achieve this, we resolve relative P and S wave polarities between pairs of waveforms by using their signed correlation coefficients—a by-product of previously performed precise earthquake relocation. We then use cluster analysis to group events with similar patterns of polarities across the network. Finally, we apply a standard mechanism inversion to the grouped data, using either catalog or correlation-derived P wave polarity data sets. This approach has great potential for enhancing analyses of spatially concentrated microseismicity such as earthquake swarms, mainshock-aftershock sequences, and industrial reservoir stimulation or injection-induced seismic sequences. To demonstrate its utility, we apply this technique to the 2014 Long Valley Caldera earthquake swarm. In our analysis, 85% of the events (7212 out of 8494 located by Shelly et al . [2016]) fall within five well-constrained mechanism clusters, more than 12 times the number with network-determined mechanisms. Of the earthquakes we characterize, 3023 (42%) have magnitudes smaller than 0.0. We find that mechanism variations are strongly associated with corresponding hypocentral structure, yet mechanism heterogeneity also occurs where it cannot be resolved by hypocentral patterns, often confined to small-magnitude events. Small (5–20°) rotations between mechanism orientations and earthquake location trends persist when we apply 3-D velocity models and might reflect a geometry of en echelon, interlinked shear, and dilational faulting.

California↗

USGS2016

Fortnightly modulation of San Andreas tremor and low-frequency earthquakes

Earth tides modulate tremor and low-frequency earthquakes (LFEs) on faults in the vicinity of the brittle−ductile (seismic−aseismic) transition. The response to the tidal stress carries otherwise inaccessible information about fault strength and rheology. Here, we analyze the LFE response to the fortnightly tide, which modulates the amplitude of the daily tidal stress over a 14-d cycle. LFE rate is highest during the waxing fortnightly tide, with LFEs most strongly promoted when the daily stress exceeds the previous peak stress by the widest margin. This pattern implies a threshold failure process, with slip initiated when stress exceeds the local fault strength. Variations in sensitivity to the fortnightly modulation may reflect the degree of stress concentration on LFE-producing brittle asperities embedded within an otherwise aseismic fault.

California↗

USGS2016

3-D P- and S-wave velocity structure and low-frequency earthquake locations in the Parkfield, California region

To refine the 3-D seismic velocity model in the greater Parkfield, California region, a new data set including regular earthquakes, shots, quarry blasts and low-frequency earthquakes (LFEs) was assembled. Hundreds of traces of each LFE family at two temporary arrays were stacked with time–frequency domain phase weighted stacking method to improve signal-to-noise ratio. We extend our model resolution to lower crustal depth with LFE data. Our result images not only previously identified features but also low velocity zones (LVZs) in the area around the LFEs and the lower crust beneath the southern Rinconada Fault. The former LVZ is consistent with high fluid pressure that can account for several aspects of LFE behaviour. The latter LVZ is consistent with a high conductivity zone in magnetotelluric studies. A new Vs model was developed with S picks that were obtained with a new autopicker. At shallow depth, the low Vs areas underlie the strongest shaking areas in the 2004 Parkfield earthquake. We relocate LFE families and analyse the location uncertainties with the NonLinLoc and tomoDD codes. The two methods yield similar results.

California↗

USGS2016

Aftershocks of the 2014 South Napa, California, Earthquake: Complex faulting on secondary faults

We investigate the aftershock sequence of the 2014 M W 6.0 South Napa, California, earthquake. Low-magnitude aftershocks missing from the network catalog are detected by applying a matched-filter approach to continuous seismic data, with the catalog earthquakes serving as the waveform templates. We measure precise differential arrival times between events, which we use for double-difference event relocation in a 3D seismic velocity model. Most aftershocks are deeper than the mainshock slip, and most occur west of the mapped surface rupture. While the mainshock coseismic and postseismic slip appears to have occurred on the near-vertical, strike-slip West Napa fault, many of the aftershocks occur in a complex zone of secondary faulting. Earthquake locations in the main aftershock zone, near the mainshock hypocenter, delineate multiple dipping secondary faults. Composite focal mechanisms indicate strike-slip and oblique-reverse faulting on the secondary features. The secondary faults were moved towards failure by Coulomb stress changes from the mainshock slip. Clusters of aftershocks north and south of the main aftershock zone exhibit vertical strike-slip faulting more consistent with the West Napa Fault. The northern aftershocks correspond to the area of largest mainshock coseismic slip, while the main aftershock zone is adjacent to the fault area that has primarily slipped postseismically. Unlike most creeping faults, the zone of postseismic slip does not appear to contain embedded stick-slip patches that would have produced on-fault aftershocks. The lack of stick-slip patches along this portion of the fault may contribute to the low productivity of the South Napa aftershock sequence.

California↗

USGS2016

The 2015 Fillmore earthquake swarm and possible crustal deformation mechanisms near the bottom of the eastern Ventura Basin, California

The 2015 Fillmore swarm occurred about 6 km west of the city of Fillmore in Ventura, California, and was located beneath the eastern part of the actively subsiding Ventura basin at depths from 11.8 to 13.8 km, similar to two previous swarms in the area. Template‐matching event detection showed that it started on 5 July 2015 at 2:21 UTC with an M ∼1.0 earthquake. The swarm exhibited unusual episodic spatial and temporal migrations and unusual diversity in the nodal planes of the focal mechanisms as compared to the simple hypocenter‐defined plane. It was also noteworthy because it consisted of >1400 events of M ≥0.0, with M 2.8 being the largest event. We suggest that fluids released by metamorphic dehydration processes, migration of fluids along a detachment zone, and cascading asperity failures caused this prolific earthquake swarm, but other mechanisms (such as simple mainshock–aftershock stress triggering or a regional aseismic creep event) are less likely. Dilatant strengthening may be a mechanism that causes the temporal decay of the swarm as pore‐pressure drop increased the effective normal stress, and counteracted the instability driving the swarm.

California↗

USGS2016

Fluid-faulting evolution in high definition: Connecting fault structure and frequency-magnitude variations during the 2014 Long Valley Caldera, California earthquake swarm

An extended earthquake swarm occurred beneath southeastern Long Valley Caldera between May and November 2014, culminating in three magnitude 3.5 earthquakes and 1145 cataloged events on 26 September alone. The swarm produced the most prolific seismicity in the caldera since a major unrest episode in 1997-1998. To gain insight into the physics controlling swarm evolution, we used large-scale cross-correlation between waveforms of cataloged earthquakes and continuous data, producing precise locations for 8494 events, more than 2.5 times the routine catalog. We also estimated magnitudes for 18,634 events (~5.5 times the routine catalog), using a principal component fit to measure waveform amplitudes relative to cataloged events. This expanded and relocated catalog reveals multiple episodes of pronounced hypocenter expansion and migration on a collection of neighboring faults. Given the rapid migration and alignment of hypocenters on narrow faults, we infer that activity was initiated and sustained by an evolving fluid pressure transient with a low-viscosity fluid, likely composed primarily of water and CO 2 exsolved from underlying magma. Although both updip and downdip migration were observed within the swarm, downdip activity ceased shortly after activation, while updip activity persisted for weeks at moderate levels. Strongly migrating, single-fault episodes within the larger swarm exhibited a higher proportion of larger earthquakes (lower Gutenberg-Richter b value), which may have been facilitated by fluid pressure confined in two dimensions within the fault zone. In contrast, the later swarm activity occurred on an increasingly diffuse collection of smaller faults, with a much higher b value.

California↗