Geology topics

Leif Karlstrom

Publications and source records attributed to Leif Karlstrom.

12 recordsLinked to original sources

USGS2025

State shifts in the deep Critical Zone drive landscape evolution in volcanic terrains

Understanding the near-surface environment where atmospheric and solid earth processes interact, often termed the “Critical Zone,” is important for assessing resources and building resilient societies. Here, we examine a volcanic landscape in the Oregon Cascade Range, an understudied Critical Zone setting that is host to major regional water resources, pervasive silicate weathering, and significant geohazards. We leverage a bedrock age chronosequence to show that the volcanic Critical Zone undergoes a structural shift, from depth extents of >1 km to meters, over timescales of ~1 My. We map an active groundwater volume comparable to major continental lakes, stored at the Cascade Range crest. This state shift makes volcanic landscape evolution a unique probe of deep coupling between Earth systems.

PNAS↗

USGS2024

Explosive 2018 eruptions at Kīlauea driven by a collapse-induced stomp-rocket mechanism

Explosive volcanic eruptions produce hazardous atmospheric plumes composed of tephra particles, hot gas and entrained air. Such eruptions are generally driven by magmatic fragmentation or steam expansion. However, an eruption mechanism outside this phreatic–magmatic spectrum was suggested by a sequence of 12 explosive eruptions in May 2018 at Kīlauea, Hawaii, that occurred during the early stages of caldera collapse and produced atmospheric plumes reaching 8 km above the vent. Here we use seismic inversions for reservoir pressure as a source condition for three-dimensional simulations of transient multiphase eruptive plume ascent through a conduit and stratified atmosphere. We compare the simulations with conduit ascent times inferred from seismic and infrasound data, and with plume heights from radar data. We find that the plumes are consistent with eruptions caused by a stomp-rocket mechanism involving the abrupt subsidence of reservoir roof rock that increased pressure in the underlying magma reservoir. In our model, the reservoir was overlain by a pocket of accumulated high-temperature magmatic gas and lithic debris, which were driven through a conduit approximately 600 m long to erupt particles at rates of around 3,000 m 3 s −1 . Our results reveal a distinct collapse-driven type of eruption and provide a framework for integrating diverse geophysical and atmospheric data with simulations to gain a better understanding of unsteady explosive eruptions.

Hawaii↗

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↗

USGS2023

Understanding the drivers of volcano deformation through geodetic model verification and validation

Volcano geodesy often involves the use of models to explain observed surface deformation. A variety of forward models are used, from analytical point sources to numerical simulations that consider complex magma system geometries, topography, and material properties. Various inversion methods can then be used to relate observed volcano data to models. Ideally, forward models should be verified through intercomparison, to check for implementation errors and quantify the error induced by any approximations used. Additionally, forward models and inversion methods should be validated through tests with synthetic and/or real data, to determine their ability to match data and estimate parameter values within uncertainty. However, to date, there have not been comprehensive verification and validation efforts in volcano geodesy. Here, we report on the first phase of the Drivers of Volcano Deformation (DVD) exercises, which were designed to build community involvement through web-based exercises involving calculations of static elastic displacement around pressurized magma reservoirs. The forward model exercises begin with a spherical reservoir in a homogeneous half space, then introduce topography, heterogeneous elastic properties, and spheroidal geometries. The inversion exercises provide synthetic noisy surface displacement data for a spherical reservoir in a homogeneous half space and assess consistency in estimates of reservoir location and volume/pressure change. There is variability in the results from both forward modeling and inversions, which highlights the strengths and limitations of different forward models, as well as the importance of inversion method choice and uncertainty quantification. This first phase of the DVD exercises serves as a community resource and will facilitate further efforts to develop standards of reproducibility.

Bulletin of Volcanology↗

USGS2022

Evolving magma temperature and volatile contents over the 2008–2018 summit eruption of Kīlauea Volcano

Magma rheology and volatile contents exert primary and highly nonlinear controls on volcanic activity. Subtle changes in these magma properties can modulate eruption style and hazards, making in situ inference of their temporal evolution vital for volcano monitoring. Here, we study thousands of impulsive magma oscillations within the shallow conduit and lava lake of Kīlauea Volcano, Hawai‘i, USA, over the 2008–2018 summit eruptive sequence, encoded by “very-long-period” seismic events and ground deformation. Inversion of these data with a petrologically informed model of magma dynamics reveals significant variation in temperature and highly disequilibrium volatile contents over days to years, within a transport network that evolved over the eruption. Our results suggest a framework for inferring subsurface magma dynamics associated with prolonged eruptions in near real time that synthesizes petrologic and geophysical volcano monitoring approaches.

Hawaii↗

USGS2021

Time-evolving surface and subsurface signatures of Quaternary volcanism in the Cascades Arc: Reply

No abstract available.

Geology↗

USGS2020

The prevalence and significance of offset magma reservoirs at arc volcanoes

Determining the spatial relations between volcanic edifices and their underlying magma storage zones is fundamental for characterizing long-term evolution and short-term unrest. We compile centroid locations of upper crustal magma reservoirs at 56 arc volcanoes inferred from seismic, magnetotelluric, and geodetic studies. We show that magma reservoirs are often horizontally offset from their associated volcanic edifices by multiple kilometers, and the degree of offset broadly scales with reservoir depth. Approximately 20% of inferred magma reservoir centroids occur outside of the overlying volcano’s mean radius. Furthermore, reservoir offset is inversely correlated with edifice size. Taking edifice volume as a proxy for long-term magmatic flux, we suggest that high flux or prolonged magmatism leads to more centralized magma storage beneath arc volcanoes by overprinting upper crustal heterogeneities that would otherwise affect magma ascent. Edifice volumes therefore reflect the spatial distribution of underlying magma storage, which could help guide monitoring strategies at volcanoes

Geophysical Research Letters↗

USGS2020

Time-evolving surface and subsurface signatures of Quaternary volcanism in the Cascades arc

Increased resolution of data constraining topography and crustal structures provides new quantitative ways to assess province-scale surface-subsurface connections beneath volcanoes. We used a database of mapped vents to extract edifices with known epoch ages from digital elevation models (DEMs) in the Cascades arc (western North America), deriving volumes that likely represent ∼50% of total Quaternary eruptive output. Edifice volumes and spatial vent density correlate with diverse geophysical data that fingerprint magmatic influence in the upper crust. Variations in subsurface structures consistent with volcanism are common beneath Quaternary vents throughout the arc, but they are more strongly associated with younger vents. Geophysical magmatic signatures increase in the central and southern Cascade Range (Cascades), where eruptive output is largest and vents are closely spaced. Vents and correlated crustal structures, as well as temporal transitions in the degree of spatially localized versus distributed eruptions, define centers with lateral extents of ∼100 km throughout the arc, suggesting a time-evolving spatial focusing of magma ascent.

California, Oregon, Washington↗

USGS2020

Spatiotemporal seismic structure variations associated with the 2018 Kīlauea eruption based on temporary dense geophone arrays

During the 2018 Kīlauea volcanic eruption, lava erupted from a series of new fissures in the lower East Rift Zone more than 30 km away from the summit through a dike intrusion. Between late May and early August, variations in the effusion rate at the persistent eruptive vent (Fissure 8) were observed following near‐daily summit caldera collapse events. Targeting the ongoing eruptive activity and the subsurface magma movement, we deployed a temporary dense seismic array. The observed time‐lapse changes in seismic velocity associated with the response of the summit collapse in three areas are presented in this study. The results show (1) clear spatially dependent co‐collapse velocity reductions across the newly‐intruded dike structure, (2) a gradual post‐collapse velocity increase near Fissure 8 correlated with the surge of magma supply, and (3) a gradual post‐collapse velocity increase on the summit likely associated with reservoir pressurization and crustal welding.

Hawaii↗

USGS2015

The effect of pressurized magma chamber growth on melt migration and pre-caldera vent locations through time at Mount Mazama, Crater Lake, Oregon

The pattern of eruptions at long-lived volcanic centers provides a window into the co-evolution of crustal magma transport, tectonic stresses, and unsteady magma generation at depth. Mount Mazama in the Oregon Cascades has seen variable activity over the last 400 ky, including the 50 km 3 climactic eruption at ca. 7.7 ka that produced Crater Lake caldera. The physical mechanisms responsible for the assembly of silicic magma reservoirs that are the precursors to caldera-forming eruptions are poorly understood. Here we argue that the spatial and temporal distribution of geographically clustered volcanic vents near Mazama reflects the development of a centralized magma chamber that fed the climactic eruption. Time-averaged eruption rates at Mount Mazama imply an order of magnitude increase in deep magma influx prior to the caldera-forming event, suggesting that unsteady mantle melting triggered a chamber growth episode that culminated in caldera formation. We model magma chamber–dike interactions over &sim;50 ky preceding the climactic eruption to fit the observed distribution of surface eruptive vents in space and time, as well as petrologically estimated deep influx rates. Best fitting models predict an expanding zone of dike capture caused by a growing, oblate spheroidal magma chamber with 10–30 MPa of overpressure. This growing zone of chamber influence causes closest approaching regional mafic vent locations as well as more compositionally evolved Mazama eruptions to migrate away from the climactic eruptive center, returning as observed to the center after the chamber drains during the caldera-forming eruption.

Oregon↗

USGS2014

Eruptions at Lone Star geyser, Yellowstone National Park, USA: 2. Constraints on subsurface dynamics

We use seismic, tilt, lidar, thermal, and gravity data from 32 consecutive eruption cycles of Lone Star geyser in Yellowstone National Park to identify key subsurface processes throughout the geyser's eruption cycle. Previously, we described measurements and analyses associated with the geyser's erupting jet dynamics. Here we show that seismicity is dominated by hydrothermal tremor (~5–40 Hz) attributed to the nucleation and/or collapse of vapor bubbles. Water discharge during eruption preplay triggers high-amplitude tremor pulses from a back azimuth aligned with the geyser cone, but during the rest of the eruption cycle it is shifted to the east-northeast. Moreover, ~4 min period ground surface displacements recur every 26 ± 8 min and are uncorrelated with the eruption cycle. Based on these observations, we conclude that (1) the dynamical behavior of the geyser is controlled by the thermo-mechanical coupling between the geyser conduit and a laterally offset reservoir periodically filled with a highly compressible two-phase mixture, (2) liquid and steam slugs periodically ascend into the shallow crust near the geyser system inducing detectable deformation, (3) eruptions occur when the pressure decrease associated with overflow from geyser conduit during preplay triggers an unstable feedback between vapor generation (cavitation) and mass discharge, and (4) flow choking at a constriction in the conduit arrests the runaway process and increases the saturated vapor pressure in the reservoir by a factor of ~10 during eruptions.

Wyoming↗

USGS2013

Eruptions at Lone Star Geyser, Yellowstone National Park, USA, part 1: energetics and eruption dynamics

Geysers provide a natural laboratory to study multiphase eruptive processes. We present results from a four–day experiment at Lone Star Geyser in Yellowstone National Park, USA. We simultaneously measured water discharge, acoustic emissions, infraredintensity, and visible and infrared video to quantify the energetics and dynamics of eruptions, occurring approximately every three hours. We define four phases in the eruption cycle: 1) a 28 ± 3 minute phase with liquid and steam fountaining, with maximum jet velocities of 16–28 m s − 1 , steam mass fraction of less than &sim; 0.01. Intermittently choked flow and flow oscillations with periods increasing from 20 to 40 s are coincident with a decrease in jet velocity and an increase of steam fraction; 2) a 26 ± 8 minute post–eruption relaxation phase with no discharge from the vent, infrared (IR) and acoustic power oscillations gliding between 30 and 40 s; 3) a 59 ± 13 minute recharge period during which the geyser is quiescent and progressively refills, and 4) a 69 ± 14 minute pre–play period characterized by a series of 5–10 minute–long pulses of steam, small volumes of liquid water discharge and 50–70 s flow oscillations. The erupted waters ascend froma 160 − 170° C reservoir and the volume discharged during the entire eruptive cycle is 20.8 ± 4.1 m 3 . Assuming isentropic expansion, we calculate a heat output from the geyser of 1.4–1.5 MW, which is < 0.1% of the total heat output from Yellowstone Caldera.

Wyoming↗