Geology topics

Peter M. Shearer

Publications and source records attributed to Peter M. Shearer.

10 recordsLinked to original sources

USGS2024

SCEC/USGS Community Stress Drop Validation Study: How spectral fitting approaches influence measured source parameters

Spectral source parameters used to estimate an earthquake’s stress drop (Δσ) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress-Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low-frequency displacement (Ω0) and corner frequency ( f c ) by fitting a smooth function to benchmark displacement spectra. Participants used linear- or log-sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate Ω0, f c , and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in Ω0 and f c can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade-off between Ω0 and f c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω0 and somewhat higher f c compared to those assuming a Brune model, although resulting Δσ estimates are similar. When compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ.

California↗

USGS2024

SKHASH: A python package for computing earthquake focal mechanisms

We introduce a Python package for computing focal mechanism solutions. This algorithm, which we refer to as SKHASH, is largely based on the HASH algorithm originally written in Fortran over 20 yr ago. HASH innovated the use of suites of solutions, spanning the expected errors in polarities and takeoff angles, to estimate focal mechanism uncertainty. SKHASH benefits from new features with flexible input formats and allows users to take advantage of recent advances in constraining focal mechanisms for small magnitude or poorly recorded earthquakes. The 3D locations of earthquakes and the velocity models used are varied when finding acceptable solutions. As a result, source–receiver azimuths are reflective of errors from the earthquake locations and velocity models, in addition to the takeoff angles. Users can consider weighted P ‐wave first‐motion polarities derived from traditional or machine‐learning picks, cross‐correlation consensus, and/or imputation techniques using SKHASH. Focal mechanism solutions can also be further constrained using traditional, machine learning, and/or cross‐correlation consensus S / P amplitude ratios. With improved reporting of individual and collective P polarity and S / P amplitude misfits, users can better evaluate the success of the solutions and the quality of the measurements. The reporting also makes it easier to identify potential issues with metadata, including incorrectly reported station polarity reversals. In addition, by leveraging vectorized operations, taking advantage of an efficient backend Python C Application Programming Interface, and the use of a parallel environment, the Python SKHASH routine may compute mechanisms quicker than the HASH routine.

Seismological Research Letters↗

USGS2021

Does earthquake stress drop increase with depth in the crust?

We combine earthquake spectra from multiple studies to investigate whether the increase in stress drop with depth often observed in the crust is real, or an artifact of decreasing attenuation (increasing Q ) with depth. In many studies, empirical path and attenuation corrections are assumed to be independent of the earthquake source depth. We test this assumption by investigating whether a realistic increase in Q with depth (as is widely observed) could remove some of the observed apparent increase in stress drop with depth. We combine event spectra, previously obtained using spectral decomposition methods, for over 50,000 earthquakes (M0 to M5) from 12 studies in California, Nevada, Kansas and Oklahoma. We find that the relative high-frequency content of the spectra systematically increases with increasing earthquake depth, at all magnitudes. By analyzing spectral ratios between large and small events as a function of source depth, we explore the relative importance of source and attenuation contributions to this observed depth dependence. Without any correction for depth-dependent attenuation, we find a systematic increase in stress drop, rupture velocity, or both, with depth, as previously observed. When we add an empirical, depth-dependent attenuation correction, the depth dependence of stress drop systematically decreases, often becoming negligible. The largest corrections are observed in regions with the largest seismic velocity increase with depth. We conclude that source parameter analyses, whether in the frequency or time domains, should not assume path terms are independent of source depth, and should more explicitly consider the effects of depth-dependent attenuation.

Journal of Geophysical Research↗

USGS2020

Abundant spontaneous and dynamically triggered submarine landslides in the Gulf of Mexico

Submarine landslides that occur offshore are common along the U.S. continental margins. These mass wasting events can trigger tsunamis and hence potentially devastate coastal communities and damage offshore infrastructure. However, the initiation and failure processes of submarine landslides are poorly understood. Here, we identify and locate 85 previously unknown submarine landslides in the Gulf of Mexico from 2008 to 2015. Ten of these landslides failed spontaneously while the remaining 75 were dynamically triggered by passing seismic surface waves from distant earthquakes with magnitudes as small as ∼5. Our observations demonstrate ongoing submarine landslide activity in the Gulf of Mexico where dense energy industry infrastructure is present and that the region is prone to secondary seismic hazard despite the low local seismicity rate. Our results should facilitate future investigations to identify unstable offshore slopes, to illuminate dynamic processes of landslides, and perhaps to apply remote detection technology in tsunami warning systems.

Geophysical Research Letters↗

USGS2017

Source spectral properties of small-to-moderate earthquakes in southern Kansas

The source spectral properties of injection-induced earthquakes give insight into their nucleation, rupture processes, and influence on ground motion. Here we apply a spectral decomposition approach to analyze P -wave spectra and estimate Brune-type stress drop for more than 2000 M L 1.5–5.2 earthquakes occurring in southern Kansas from 2014 to 2016. We find that these earthquakes are characterized by low stress drop values (median ∼0.4MPa) compared to natural seismicity in California. We observe a significant increase in stress drop as a function of depth, but the shallow depth distribution of these events is not by itself sufficient to explain their lower stress drop. Stress drop increases with magnitude from M1.5–M3.5, but this scaling trend may weaken above M4 and also depends on the assumed source model. Although we observe a nonstationary, sequence-specific temporal evolution in stress drop, we find no clear systematic relation with the activity of nearby injection wells.

Kansas↗

USGS2015

Source mechanism of small long-period events at Mount St. Helens in July 2005 using template matching, phase-weighted stacking, and full-waveform inversion

Long-period (LP, 0.5-5 Hz) seismicity, observed at volcanoes worldwide, is a recognized signature of unrest and eruption. Cyclic LP “drumbeating” was the characteristic seismicity accompanying the sustained dome-building phase of the 2004–2008 eruption of Mount St. Helens (MSH), WA. However, together with the LP drumbeating was a near-continuous, randomly occurring series of tiny LP seismic events (LP “subevents”), which may hold important additional information on the mechanism of seismogenesis at restless volcanoes. We employ template matching, phase-weighted stacking, and full-waveform inversion to image the source mechanism of one multiplet of these LP subevents at MSH in July 2005. The signal-to-noise ratios of the individual events are too low to produce reliable waveform-inversion results, but the events are repetitive and can be stacked. We apply network-based template matching to 8 days of continuous velocity waveform data from 29 June to 7 July 2005 using a master event to detect 822 network triggers. We stack waveforms for 359 high-quality triggers at each station and component, using a combination of linear and phase-weighted stacking to produce clean stacks for use in waveform inversion. The derived source mechanism pointsto the volumetric oscillation (~10 m 3 ) of a subhorizontal crack located at shallow depth (~30 m) in an area to the south of Crater Glacier in the southern portion of the breached MSH crater. A possible excitation mechanism is the sudden condensation of metastable steam from a shallow pressurized hydrothermal system as it encounters cool meteoric water in the outer parts of the edifice, perhaps supplied from snow melt.

Washington↗

USGS2014

Three-dimensional seismic velocity structure of Mauna Loa and Kilauea volcanoes in Hawaii from local seismic tomography

We present a new three-dimensional seismic velocity model of the crustal and upper mantle structure for Mauna Loa and Kilauea volcanoes in Hawaii. Our model is derived from the first-arrival times of the compressional and shear waves from about 53,000 events on and near the Island of Hawaii between 1992 and 2009 recorded by the Hawaiian Volcano Observatory stations. The V p model generally agrees with previous studies, showing high-velocity anomalies near the calderas and rift zones and low-velocity anomalies in the fault systems. The most significant difference from previous models is in V p / V s structure. The high- V p and high- V p / V s anomalies below Mauna Loa caldera are interpreted as mafic magmatic cumulates. The observed low- V p and high- V p / V s bodies in the Kaoiki seismic zone between 5 and 15 km depth are attributed to the underlying volcaniclastic sediments. The high- V p and moderate- to low- V p / V s anomalies beneath Kilauea caldera can be explained by a combination of different mafic compositions, likely to be olivine-rich gabbro and dunite. The systematically low- V p and low- V p / V s bodies in the southeast flank of Kilauea may be caused by the presence of volatiles. Another difference between this study and previous ones is the improved V p model resolution in deeper layers, owing to the inclusion of events with large epicentral distances. The new velocity model is used to relocate the seismicity of Mauna Loa and Kilauea for improved absolute locations and ultimately to develop a high-precision earthquake catalog using waveform cross-correlation data.

Hawaii↗

USGS2014

High-precision relocation of long-period events beneath the summit region of Kı̄lauea Volcano, Hawai‘i, from 1986 to 2009

Long-period (0.5–5 Hz, LP) seismicity has been recorded for decades in the summit region of Kı̄lauea Volcano, Hawai‘i, and is postulated as linked with the magma transport and shallow hydrothermal systems. To better characterize its spatiotemporal occurrence, we perform a systematic analysis of 49,030 seismic events occurring in the Kı̄lauea summit region from January 1986 to March 2009 recorded by the &sim;50-station Hawaiian Volcano Observatory permanent network. We estimate 215,437 P wave spectra, considering all events on all stations, and use a station-averaged spectral metric to consistently classify LP and non-LP seismicity. We compute high-precision relative relocations for 5327 LP events (43% of all classified LP events) using waveform cross correlation and cluster analysis with 6.4 million event pairs, combined with the source-specific station term method. The majority of intermediate-depth (5–15 km) LPs collapse to a compact volume, with remarkable source location stability over 23 years indicating a source process controlled by geological or conduit structure.

Hawaii↗

USGS2011

Scattered P'P' waves observed at short distances

We detect previously unreported 1 Hz scattered waves at epicentral distances between 30° and 50° and at times between 2300 and 2450 s after the earthquake origin. These waves likely result from off-azimuth scattering of PKPbc to PKPbc in the upper mantle and crust and provide a new tool for mapping variations in fine-scale (10 km) mantle heterogeneity. Array beams from the Large Aperture Seismic Array (LASA) clearly image the scattered energy gradually emerging from the noise and reaching its peak amplitude about 80 s later, and returning to the noise level after 150 s. Stacks of transverse versus radial slowness (ρt, ρr) show two peaks at about (2, -2) and (-2,-2) s/°, indicating the waves arrive along the major arc path (180° to 360°) and significantly off azimuth. We propose a mantle and surface PKPbc to PKPbc scattering mechanism for these observations because (1) it agrees with the initiation time and distinctive slowness signature of the scattered waves and (2) it follows a scattering path analogous to previously observed deep-mantle PK•KP scattering (Chang and Cleary, 1981). The observed upper-mantle scattered waves and PK•KP waves fit into a broader set of scattered waves that we call P′•d•P′, which can scatter from any depth, d, in the mantle.

Bulletin of the Seismological Society of America↗

USGS2004

The global short-period wavefield modelled with a Monte Carlo seismic phonon method

At high frequencies (∼1 Hz), much of the seismic energy arriving at teleseismic distances is not found in the main phases (e.g. P , PP , S , etc.) but is contained in the extended coda that follows these arrivals. This coda results from scattering off small-scale velocity and density perturbations within the crust and mantle and contains valuable information regarding the depth dependence and strength of this heterogeneity as well as the relative importance of intrinsic versus scattering attenuation. Most analyses of seismic coda to date have concentrated on S -wave coda generated from lithospheric scattering for events recorded at local and regional distances. Here, we examine the globally averaged vertical-component, 1-Hz wavefield (>10° range) for earthquakes recorded in the IRIS FARM archive from 1990 to 1999. We apply an envelope-function stacking technique to image the average time–distance behavior of the wavefield for both shallow (≤50 km) and deep (≥500 km) earthquakes. Unlike regional records, our images are dominated by P and P coda owing to the large effect of attenuation on PP and S at high frequencies. Modelling our results is complicated by the need to include a variety of ray paths, the likely contributions of multiple scattering and the possible importance of P -to- S and S -to- P scattering. We adopt a stochastic, particle-based approach in which millions of seismic phonons are randomly sprayed from the source and tracked through the Earth. Each phonon represents an energy packet that travels along the appropriate ray path until it is affected by a discontinuity or a scatterer. Discontinuities are modelled by treating the energy normalized reflection and transmission coefficients as probabilities. Scattering probabilities and scattering angles are computed in a similar fashion, assuming random velocity and density perturbations characterized by an exponential autocorrelation function. Intrinsic attenuation is included by reducing the energy contained in each particle as an appropriate function of traveltime. We find that most scattering occurs in the lithosphere and upper mantle, as previous results have indicated, but that some lower-mantle scattering is likely also required. A model with 3 to 4 per cent rms velocity heterogeneity at 4-km scale length in the upper mantle and 0.5 per cent rms velocity heterogeneity at 8-km scale length in the lower mantle (with intrinsic attenuation of Q α = 450 above 200 km depth and Q α = 2500 below 200 km) provides a reasonable fit to both the shallow- and deep-earthquake observations, although many trade-offs exist between the scale length, depth extent and strength of the heterogeneity.

Geophysical Journal International↗