Geology topics

Brandon Schmandt

Publications and source records attributed to Brandon Schmandt.

10 recordsLinked to original sources

USGS2026

Small earthquake moment magnitude and implications for frequency–magnitude scaling of injection induced earthquakes of the Raton Basin

Accurate estimation of earthquake source parameters—such as moment magnitudes, corner frequencies, and stress drops—is essential for improving seismic hazard assessments and understanding earthquake physics. In this study, moment magnitudes ( M W ) are calculated for 31,581 earthquakes associated with wastewater injection in the Raton Basin (located along the border between northern New Mexico and southern Colorado) between 2016 and 2024 using radiative transfer theory to fit coda decay envelopes. Our results show that it is feasible to estimate moment magnitudes down to M W ~1 with coda envelopes from a small local monitoring network. Significant differences were found between M W and local magnitudes ( M L ) for small earthquakes ( M < 3.0). A linear relationship was optimized to convert M L to M W : M W = 0.7 M L + 0.96 and M W = 0.73 M L + 0.99 (for the events reported by the U.S. Geological Survey), which can be applied in future studies of Raton Basin seismicity. We find that b -values calculated employing different methods and using M L are approximately 1.0, while those using M W range from 1.2 to 1.4. A larger estimate of the b -value could influence interpretations of the statistical behavior of earthquakes associated with injection and consequently seismic hazard assessments based on a magnitude–frequency distribution. The potential differences between local versus moment magnitude-based earthquake statistics should be considered in other seismically active regions.

Colorado, New Mexico↗

USGS2025

Evidence for fluid pressurization of fault zones and persistent sensitivity to injection rate beneath the Raton Basin

Subsurface wastewater injection has increased the seismicity rate within the Raton Basin over more than two decades, with the basin-wide injection rate peaked between 2009-2015. To understand the evolution of injection-induced earthquakes, we systematically analyzed 2016-2024 broadband recordings with a machine-learning-based phase picker and constructed a catalog with 95,993 earthquakes (-1≤ M L ≤4.3). We then inverted for full centroid moment tensors (CMT) for 90 M L ≥ 2 events, with a special interest in constraining the non-double-couple components via probabilistic metrics. Both relocations and CMT solutions support basement-rooted normal faults, including graben and half-graben structures. Furthermore, we observe the non-double-couple components that imply elevated pore pressure in the fault zones. An earthquake cluster emerged in the north-central basin in 2023, preceded by ~1-yr of increased injection volume from wells within 15km. Despite a basin-wide decrease in the injection volume, we highlights the persistence of seismicity that remains to sensitive to injection rates within the Raton Basin.

Colorado, New Mexico↗

USGS2024

(Re)discovering the seismicity of Antarctica: A new seismic catalog for the southernmost continent

We apply a machine learning (ML) earthquake detection technique on over 21 yr of seismic data from on‐continent temporary and long‐term networks to obtain the most complete catalog of seismicity in Antarctica to date. The new catalog contains 60,006 seismic events within the Antarctic continent for 1 January 2000–1 January 2021, with estimated moment magnitudes (⁠Mw ⁠) between −1.0 and 4.5. Most detected seismicity occurs near Ross Island, large ice shelves, ice streams, ice‐covered volcanoes, or in distinct and isolated areas within the continental interior. The event locations and waveform characteristics indicate volcanic, tectonic, and cryospheric sources. The catalog shows that Antarctica is more seismically active than prior catalogs would indicate, examples include new tectonic events in East Antarctica, seismic events near and around the vicinity of David Glacier, and many thousands of events in the Mount Erebus region. This catalog provides a resource for more specific studies using other detection and analysis methods such as template matching or transfer learning to further discriminate source types and investigate diverse seismogenic processes across the continent.

Seismological Research Letters↗

USGS2024

Background seismic noise levels among the Caribbean network and the role of station proximity to coastline

The amplitude and frequency content of background seismic noise is highly variable with geographic location. Understanding the characteristics and behavior of background seismic noise as a function of location can inform approaches to improve network performance and in turn increase earthquake detection capabilities. Here, we calculate power spectral density estimates in one‐hour windows for over 15 yr of vertical‐component data from the nine‐station Caribbean network (CU) and look at background noise within the 0.05–300 s period range. We describe the most visually apparent features observed at the CU stations. One of the most prominent features occurs in the 0.75–3 s band for which power levels are systematically elevated and decay as a function of proximity to the coastline. Further examination of this band on 1679 contiguous USArray Transportable Array stations reveals the same relationship. Such a relationship with coastal distance is not observed in the 4–8 s range more typical of globally observed secondary microseisms. A simple surface‐wave amplitude decay model fits the observed decay well with geometric spreading as the most important factor for stations near the coast (<∼50 km). The model indicates that power levels are strongly influenced by proximity to coastline at 0.75–3 s. This may be because power from nearshore wave action at 0.75–3 s overwhelms more distant and spatially distributed secondary microseism generation. Application of this basic model indicates that a power reduction of ∼25 dB can be achieved by simply installing the seismometer 25 km away from the coastline. This finding may help to inform future site locations and array design thereby improving network performance and data quality, and subsequently earthquake detection capabilities.

Seismological Research Letters↗

USGS2024

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗

USGS2024

The magmatic origin of the Columbia River Gorge, USA

Along subduction zones, high-relief topography is associated with sustained volcanism parallel to the plate margin. However, the relationship between magmatism and mountain building in arcs is poorly understood. Here, we study patterns of surface deformation and correlated fluvial knickpoints in the Columbia River Gorge to link long-term magmatism to the uplift and ensuing topographic development of the Cascade Range. An upwarped paleochannel exposed in the walls of the Gorge constrains unsteady deep magma flux, the ratio of intrusive to extrusive magmatic contributions to topography, and the impact of magmatism on Columbia River incision since 3.5 million years ago. Geophysical data indicate that deep magma influx beneath the arc axis is ongoing and not aligned with the current locations of volcanic edifices, representing a broad regional influence on arc construction.

Oregon, Washington↗

USGS2022

Revisiting the depth distribution of seismicity before and after the 2004–2008 eruption of Mount St. Helens

Hypocenter estimation at active volcanoes improves our understanding of their magmatic systems and indicates changing conditions at depth for continuously monitored volcanoes. The most active volcano in the Cascades Range, Mount St. Helens, has a multi-decadal earthquake catalog and it shows an apparent change in the depth distribution of seismicity before and after the 2004–2008 dome-building eruption and unrest sequence. We use two new resources to evaluate the accuracy of hypocenters and consequently the change in depth distribution of seismicity before and after the 2004–2008 eruption. First, we deployed a dense array of 136 three-component nodal seismographs for one month in 2017, including sub-arrays on the newly extruded dome and crater floor. Second, for events recorded during this month, we located their hypocenters using a three-dimensional (3D) wavefront-tracking location solver and a recently developed tomography model derived from active and passive source data. The relocated hypocenters are generally shallower and more concentrated beneath the crater compared to their catalog locations. The mean hypocenter movement from catalog locations is 2.86 km, with an averaged depth shift of 2.53 km upward. Comparison between 2017 hypocenters located using all of the available phase picks and those located using only the catalog picks from the permanent network suggests the improved hypocenters mostly resulted from the location solver and the 3D velocity model, with smaller changes due to the dense three-component array. Applying the same hypocenter estimation method to phase picks for all events from 1997 to 2004 and 2008–2021, we found a concentration of seismicity between sea level and ~ 1 km above it before and after the 2004–2008 eruption. The new results suggest that the seismogenic structure in the shallow magmatic system quickly re-equilibrated to its earlier state after the dome-building eruption.

Washington↗

USGS2020

Local earthquake Vp and Vs tomography in the Mount St. Helens region with the iMUSH broadband array

We present new 3-D P wave and S wave velocity models of the upper 20 km of the Mount St. Helens (MSH) region. These were obtained using local-source arrival time tomography from earthquakes and explosions recorded at 70 broadband stations deployed as part of the imaging Magma Under St. Helens (iMUSH) project and augmented by several data sets. Principal features of our models include (1) low P wave and S wave velocities along the St. Helens seismic zone to depths of at least 20 km corresponding to high conductivity imaged by iMUSH magnetotelluric studies. This delineates a zone of weakness that magma can exploit at the location of MSH; (2) a 5- to 7-km diameter, 6–15 km deep, 3–6% negative P wave and S wave velocity anomaly beneath MSH, consistent with previous estimates of the source region for recent eruptions. We interpret this as a magma storage region containing up to 15–20 km 3 of partial melt, which is about 5 times more than the largest documented eruption at MSH; (3) a broad region of low P wave velocity below 10-km depth extending between Mount Adams and Mount Rainier along and to the east of the main Cascade arc, which is likely due to high-temperature arc crust and possible presence of fluids or melt; (4) several anomalies associated with surface-mapped features, including high-velocity igneous units such as the Spud Mountain and Spirit Lake plutons and low velocities in the Chehalis sedimentary basin and the Indian Heaven volcanic field. Our results place further constraints on the geometry of these features at depth.

Washington↗

USGS2018

Rayleigh and S wave tomography constraints on subduction termination and lithospheric foundering in central California

The crust and upper mantle structure of central California have been modified by subduction termination, growth of the San Andreas plate boundary fault system, and small-scale upper mantle convection since the early Miocene. Here we investigate the contributions of these processes to the creation of the Isabella Anomaly, which is a high seismic velocity volume in the upper mantle. There are two types of hypotheses for its origin. One is that it is the foundered mafic lower crust and mantle lithosphere of the southern Sierra Nevada batholith. The alternative suggests that it is a fossil slab connected to the Monterey microplate. A dense broadband seismic transect was deployed from the coast to the western Sierra Nevada to fill in the least sampled areas above the Isabella Anomaly, and regional-scale Rayleigh and S wave tomography are used to evaluate the two hypotheses. New shear velocity (Vs) tomography images a high-velocity anomaly beneath coastal California that is sub-horizontal at depths of ∼40–80 km. East of the San Andreas Fault a continuous extension of the high-velocity anomaly dips east and is located beneath the Sierra Nevada at ∼150–200 km depth. The western position of the Isabella Anomaly in the uppermost mantle is inconsistent with earlier interpretations that the Isabella Anomaly is connected to actively foundering foothills lower crust. Based on the new Vs images, we interpret that the Isabella Anomaly is not the dense destabilized root of the Sierra Nevada, but rather a remnant of Miocene subduction termination that is translating north beneath the central San Andreas Fault. Our results support the occurrence of localized lithospheric foundering beneath the high elevation eastern Sierra Nevada, where we find a lower crustal low Vs layer consistent with a small amount of partial melt. The high elevations relative to crust thickness and lower crustal low Vs zone are consistent with geological inferences that lithospheric foundering drove uplift and a ∼3–4 Ma pulse of basaltic magmatism.

Earth and Planetary Science Letters↗

USGS2014

Steady incision of Grand Canyon at the million year timeframe: A case for mantle-driven differential uplift

The Grand Canyon region provides an excellent laboratory to examine the interplay between river incision, magmatism, and the geomorphic and tectonic processes that shape landscapes. Here we apply U-series, Ar–Ar, and cosmogenic burial dating of river terraces to examine spatial variations in incision rates along the 445 km length of the Colorado River through Grand Canyon. We also analyze strath terrace sequences that extend to heights of several hundred meters above the river, and integrate these with speleothem constrained maximum incision rates in several reaches to examine any temporal incision variations at the million-year time frame. This new high-resolution geochronology shows temporally steady long-term incision in any given reach of Grand Canyon but significant variations along its length from 160 m/Ma in the east to 101 m/Ma in the west. Spatial and temporal patterns of incision, and the long timescale of steady incision rule out models where geomorphic controls such as climate oscillations, bedrock strength, sediment load effects, or isostatic response to differential denudation are the first order drivers of canyon incision. The incision pattern is best explained by a model of Neogene and ongoing epeirogenic uplift due to an eastward propagating zone of increased upper mantle buoyancy that we infer from propagation of Neogene basaltic volcanism and a strong lateral gradient in modern upper mantle seismic structure.

Arizona↗