Cascades Volcano Observatory

Abstract More than 3 m.y. of mafic volcanism near the Goat Rocks volcanic cluster in the southern Washington Cascade Range, USA, lends insight into the evolution of basalts and the subarc mantle at a long-lived, major arc volcanic locus. We contribute field observations, 40Ar/39Ar dates, paleomagnetic directions, and bulk rock and mineral compositions to characterize nine mafic units that erupted in association with the Goat Rocks volcanic cluster. The time frame of mafic volcanism, ca. 3.6 Ma to 60 ka, encompasses the lifespan of the central volcanic cluster (3.1 Ma to 115 ka), with a lull from ca. 2.7 Ma to 1.4 Ma. A climactic period of voluminous mafic activity and far-traveled lava flows, including construction of the Hogback Mountain shield volcano, coincided with voluminous andesite eruptions from the central volcanic cluster. The basaltic rocks in the Goat Rocks area are calc-alkaline to barely tholeiitic and have high field strength element depletion relative to large-ion lithophile elements characteristic of calc-alkaline basalts (CAB) of the Cascade volcanic arc. Unlike at neighboring andesitic volcanic centers (Mounts Adams, St. Helens, and Rainier), no other mafic end members such as high-aluminum olivine tholeiite (HAOT) or intraplate-type basalt (IPB) are present at or near the Goat Rocks volcanic cluster, although some of the calc-alkaline basalts in this study have IPB-like affinities. The Goat Rocks mafic units exhibit two main temporal trends in composition: (1) the most primitive basalts erupted earlier, compared to less primitive and more evolved compositions later, and (2) high field strength element concentrations are higher in the younger basalt units relative to the oldest two. In contrast to these temporal trends, the mafic units define two compositional groups that recur through time, a low-Sr and a high-Sr group, each with distinct trace element and Sr and Nd isotope ratios. Although radiogenic isotope ratios are generally aligned with High Cascades CAB and HAOT, some extend toward IPB of Mount Adams and Simcoe Mountains volcanic field. Olivine-dominated crystal fractionation at shallow pressure from a small range of parent magma compositions accounts for much of the variation among the basalts and basaltic andesites. A high-pressure fractionation model is plausible for only one of the youngest basalt units (basalt of Walupt Lake volcano). Mafic recharge and crustal assimilation accounts for the incompatible-element enriched composition of basaltic andesites erupted during construction of the largest andesitic centers, further supporting sustained basalt mass flux and thermal energy driving andesite genesis. We model the most primitive members of the Goat Rocks mafic units as partial melts of successively less depleted mantle in time. Variable degrees of fluxing with fluids and melts from subduction explain the distinction between high-Sr and low-Sr groups. We propose that mantle metasomatism by ancestral subduction and fluid-flux melting is heterogeneously distributed through the local subarc mantle and played a greater role in the genesis of the high-Sr basalt group. The limited range of primitive basalt types around the Goat Rocks volcanic cluster contrasts with the much greater diversity of basalts throughout the southern Washington to northern Oregon Cascade arc. On the other hand, the central volcanic cluster encompasses nearly the entire diversity observed at neighboring composite volcanoes. In the case of the Goat Rocks area at least, and perhaps attributable to the entire region, this means that the genesis of diverse intermediate magmas is independent from and does not require vastly different parental basalt compositions.

Abstract 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 and 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 ML 4.3). We then inverted for full centroid moment tensors (CMT) for 90 ML 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 15 km. Despite a basin‐wide decrease in the injection volume, we highlight the persistence of seismicity that remains sensitive to injection rates within the Raton Basin.

Abstract Detecting and cataloging seismic events are among the most fundamental tasks in seismology. Many standardized tools for these tasks exist, including the open-source package repeating earthquake detector in Python (REDPy). REDPy generates an organized catalog of seismic events from continuous waveform data, in which events are automatically separated into groups (“families”) by their waveform similarity through cross-correlation. REDPy also automatically generates various outputs that allow a user to visualize important trends in the catalog, which may be used in real time or in retrospective analyses to allow rapid identification of interesting features. The code was designed for near-real-time volcano monitoring but is applicable across a broad range of use cases in seismology and seismoacoustics. In this article, the utility and performance of REDPy are demonstrated on two highly seismogenic volcanic eruption sequences: the onset of the dome-building eruption of Mount St. Helens, Washington, from 2004 to 2005, and the entirety of the summit caldera collapse sequence of Kīlauea, Hawai‘i, in 2018. This article is meant to be a companion to the documentation of the code; in addition to detailing the basic required inputs, script functionality, and resulting outputs, the reasonings behind several important design decisions are also discussed.

Emplacement of the Siberian Traps large igneous province (LIP) around 252 Ma coincided with the most profound environmental disruption of the past 500 million years. The enormous volume of the Siberian Traps, its ability to generate greenhouse gases and other volatiles, and a temporal coincidence with extinction all suggest a causal link. Patterns of marine and terrestrial extinction/recovery are consistent with environmental stresses potentially triggered by the Siberian Traps. However, the nature of causal links between the LIP and mass extinction remains enigmatic. Understanding the origins, anatomy, and forcing potential of the Siberian Traps LIP and the spatiotemporal patterns of resulting stresses represents a critical counterpart to high-resolution fossil and proxy records of Permian–Triassic environmental and biotic shifts. This review provides a summary of recent advances and key questions regarding the Siberian Traps in an effort to illuminate what combination of factors made the Siberian Traps a uniquely deadly LIP. ▪ Large igneous provinces such as the Siberian Traps are capable of triggering global environmental destabilization. ▪ Greenhouse gases generated by Siberian Traps magmatism played a major role in driving the climate changes that triggered the end-Permian mass extinction. ▪ The end-Permian extinction fundamentally altered the evolutionary trajectory of Earth's biosphere, creating ecological space for many of the organisms seen today. ▪ Determining the timing and patterns of end-Permian marine and terrestrial mortality and recovery and the timing and character of Siberian Traps magmatism is key in understanding the causal link between magmatism and extinction. ▪ Understanding the cause(s) of past extinction events informs hypotheses about current and future environmental destabilization.

Summary Infrasound (acoustic waves below 20 Hz) can be used to detect, locate, and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height, and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local (<5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole, and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modeling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analyzing geophysical field-scale experimental data; however, especially at very local distances from the source (<5 km), three-dimensional simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications.