Geology topics

Josef Dufek

Publications and source records attributed to Josef Dufek.

10 recordsLinked to original sources

USGS2026

PDCD-DAT: A global database of pyroclastic density current deposit field data

Pyroclastic density currents represent one of the deadliest hazards posed by active volcanoes. Analysis of their deposits provides valuable insights into their internal dynamics and informs numerical simulations of pyroclastic density currents which underpin many volcanic hazard assessments. We present PDCD-DAT, a global database of pyroclastic density current deposit characteristics compiled from peer-reviewed literature. The database includes both quantitative datasets (e.g., grain size, density, bedform dimensions, thickness) and qualitative descriptors (e.g., sedimentary structures, lithofacies). PDCD-DAT includes data from 85 source publications, covering 97 eruptions or eruptive phases, and 214 individual depositional units from 55 globally distributed volcanoes. Eruptions recorded in the database range from VEI 1–8. We highlight examples of potential applications of the database, which include (i) comparison of single deposit case studies to global datasets, (ii) informing input parameters and conditions for numerical and analogue models of pyroclastic density currents, (iii) validation of numerical and analogue models against a wide variety of natural case study deposit architectures, and (iv) estimating hazard impact metrics of pyroclastic density currents from past eruptions. We show that the database represents a useful tool for improving our ability to model pyroclastic density currents, predict their associated hazards, and understand the relationships between the internal dynamics of pyroclastic density currents and the properties of their deposits. PDCD-DAT is integrated with the FlowDat Mass Flow Database, which provides a sustainable platform for the database. We aim for PDCD-DAT to be expanded in the future through addition of pyroclastic density current deposit datasets from new field studies conducted by the volcanology research community.

Journal of Applied Volcanology↗

USGS2026

The magmatic-hydrothermal system of the Three Sisters volcanic cluster, Oregon, imaged from field gravity measurements

From 2019 to 2024, gravity surveys were conducted at the Three Sisters volcanic cluster (TSVC), measuring 246 gravity sites using a spring relative gravimeter. We calculated the residual Bouguer anomaly and identified three main zones with negative anomalies, ranging from −4 to −8 mGal, located southwest and west of South Sister, within an area that has been uplifting for the past two decades. After inversion, we obtain a 3D density model of the subsurface and identify low-density bodies extending from the surface down to 3 km. We estimate a total of 15 km 3 of crustal bodies with density close to 2 g/cm 3 that could store up to ~5 km 3 of water, forming an extensive hydrothermal system beneath the TSVC. We explore the possible combinations of melt compositions and temperatures that could create a bulk density close to our reference crustal density (2.5 g/cm 3 ) using MELTS thermodynamic simulations. Our results indicate that a magmatic mush with as little as 15% partial melt of bulk rhyolitic composition or as much as 52%–57% partial melt of a bulk dacitic composition could be stored in a magmatic system under TSVC without generating a detectable gravity anomaly. Episodic magma injections at the base of the magmatic system, such as the 1998–2000 intrusion at ~6 km depth, would bring heat and gas to the hydrothermal system while maintaining a low melt fraction in the magmatic mush, as imaged at other Cascade volcanoes.

Oregon↗

USGS2025

The transition from melt accumulation to eruption initiation recorded by orthopyroxene Fe-Mg diffusion timescales in late Holocene rhyolites, South Sister volcano, Oregon Cascade Range

South Sister volcano, Oregon Cascade Range, USA, has repeatedly erupted rhyolite since ca. 40 ka. The youngest such eruptions are the ca. 2 ka Rock Mesa and Devils Chain rhyolites, erupted several hundred years apart from two multi-vent complexes separated by 3–6 km. Fe-Mg interdiffusion models of orthopyroxene rims from both rhyolites produce timescales up to several-thousand years, but dominantly decades-to-centuries. Notably, the timescales of step-normal zoned orthopyroxene rims (i.e., normally zoned with a steep chemical gradient) from the Rock Mesa rhyolite are longer than those of reversely zoned crystals, whereas the Devils Chain produced mostly decadal timescales for both zoning types. Despite the proximity and broadly similar products of these episodes, their respective timescales indicate distinct sequences of events leading up to each eruption. The Rock Mesa timescales record centuries of magma chamber growth followed by decades of predominantly magma rejuvenation, reorganization, and destabilization. In contrast, the Devils Chain episode was preceded by a single episode of coupled rhyolite extraction, rejuvenation, and hybridization. Rare, high-An plagioclase cores and evidence of reheating implicate cryptic emplacement of mafic magma at the base of the rhyolite reservoirs. However, the diffusion timescales do not unequivocally support a single magma recharge event that affected both. Fluid fluxing and the reorganization of melt into buoyant magma chambers likely provided the source of increasing pressurization that initiated each eruption after several decades. Geodetic models of ongoing deformation west of South Sister could consider these processes in addition to magma emplacement.

Oregon↗

USGS2024

Zircon constraints on the eruptive sequence and magma evolution of rhyolites at South Sister volcano, Oregon

We present 230 Th- 238 U crystallization ages and trace element compositions for zircons spanning the late Pleistocene to Holocene rhyolite eruptive record at South Sister volcano in the central Oregon Cascade Range. Most zircon ages are between 100 and 20 ka, with very few in secular equilibrium (>350 ka). The weighted mean of zircon ages for the two oldest South Sister rhyolites, 31.5 ± 2.1 and 39.1 ± 2.4 ka, are significantly younger than the associated 40 Ar/ 39 Ar ages, 47.4 ± 9.7 and 51.4 ± 9.7 ka. We propose that these 40 Ar/ 39 Ar dates, performed on plagioclase separates, are compromised by a subtle amount of excess Ar and therefore the younger weighted mean zircon ages yield more reliable eruption ages. These results imply that the interval of rhyolite eruption at South Sister during the late Pleistocene was both shorter and more productive than previously thought and that eruption at South Sister initiated after Middle Sister. Compositionally, zircons from the Pleistocene rhyolites are broadly similar and show down-temperature zircon and plagioclase crystallization trends. However, we argue that destabilized amphibole and titanite in a common mush also exert leverage on the Pleistocene zircon trace element compositions. Divergence in the Eu/Eu* ratio between the Pleistocene and Holocene lavas implies chemically distinct magma reservoirs originating from the Pleistocene rhyolite eruptive sequence and the Holocene eruptive sequence. This work suggests a higher flux of rhyolite volcanism than previously thought and characterizes magmatic storage distinctions between the Pleistocene and Holocene rhyolites, aiding in the assessment of future eruptive hazards at South Sister volcano.

Oregon↗

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↗

USGS2023

Dynamics of the December 2020 ash-poor plume formed by lava-water interaction at the summit of Kilauea Volcano, Hawaii

On 20 December 2020, after more than 2 years of quiescence at Kīlauea Volcano, Hawaiʻi, renewed volcanic activity in the summit crater caused boiling of the water lake over a period of ∼90 min. The resulting water-rich, electrified plume rose to 11–13 km above sea level, which is among the highest plumes on record for Kīlauea. Although conventional models would infer a high mass flux from explosive magma-water interaction, the plume was not associated with an infrasound signal indicative of “explosive” activity, nor did it produce a measurable ash-fall deposit. We use multisensor data to characterize lava-water interaction and plume generation during this opening phase of the 2020–21 eruption. Satellite, weather radar, and eyewitness observations revealed that the plume was rich in water vapor and hydrometeors but transported less ash than expected from its maximum height. Volcanic lightning flashes detected by ground-based cameras were confined to freezing altitudes of the upper cloud, suggesting that the ice formation drove the electrification of this plume. The low acoustic energy from lava-water interaction points to a weakly explosive style of hydrovolcanism. Heat transfer calculations show that the lava to water heat flux was sufficient to boil the lake within 90 min. Limited mixing of lava and water inhibited major steam explosions and fine fragmentation. Results from one-dimensional plume modeling suggest that the models may underpredict plume height due to overestimation of crosswind air-entrainment. Our findings shed light on an unusual style of volcanism in which weakly explosive lava-water interaction generated an outsized plume.

Hawaii↗

USGS2022

The nature of active magma reservoirs and storage underneath Cascade volcanoes

Volcanic activity is dictated by crustal pathways and reservoirs through which magma ascends and collects. The Cascades are a natural laboratory to study the variable influence of the crust on erupted magmas and subterranean reservoirs. The interaction of ascending magmas with the variable subduction geometry and tectonics of the overriding North American plate have given rise to a diversity of magmatic storage conditions. The confluence of geochemical and geophysical investigations emphasizes that most magmatic systems in the Cascades have been built at multiple levels in the crust as determined by tectonics, pre-existing structure, and magmatic flux from the mantle.

British Columbia, California, Oregon, Washington↗

USGS2017

Field-trip guide for exploring pyroclastic density current deposits from the May 18, 1980, eruption of Mount St. Helens, Washington

Pyroclastic density currents (PDCs) are one of the most dangerous phenomena associated with explosive volcanism. To help constrain damage potential, a combination of field studies, laboratory experiments, and numerical modeling are used to establish conditions that influence PDC dynamics and depositional processes, including runout distance. The objective of this field trip is to explore field relations that may constrain PDCs at the time of emplacement. The PDC deposits from the May 18, 1980, eruption of Mount St. Helens are well exposed along the steep flanks (10–30° slopes) and across the pumice plain (5–12° slopes) as far as 8 km north of the volcano. The pumice plain deposits represent deposition from a series of concentrated PDCs and are primarily thick (3–12 m), massive, and poorly sorted. In contrast, the steep east-flank deposits are stratified to cross-stratified, suggesting deposition from PDCs where turbulence strongly influenced transport and depositional processes. The PDCs that descended the west flank were largely nondepositional; they maintained a higher flow energy and carrying capacity than PDCs funneled through the main breach, as evidenced by the higher concentration of large blocks in their deposits. The PDC from the west flank collided with PDCs funneled through the breach at various points along the pumice plain. Evidence for flow collision will be explored and debated throughout the field trip. Evidence for substrate erosion and entrainment is found (1) along the steep eastern flank of the volcano, which has a higher degree of rough, irregular topography relative to the west flanks where PDCs were likely nonerosive, (2) where PDCs encountered debris-avalanche hummocks across the pumice plain, and (3) where PDCs eroded and entrained material deposited by PDCs produced during earlier phases of the eruption. Two features interpreted as large-scale (tens of meters wide) levees and a large (~200 m wide) channel scour-and-fill feature provide the first evidence of self-channelization within PDCs sustained for minutes to tens of minutes (total volume of deposits is ~0.12 km 3 ; area covered is ~15.5 km 2 ; Rowley and others, 1981). Our ability to interpret the deposits of PDCs is critical for understanding transport and depositional processes that control PDC dynamics. The results of extensive work on the May 18, 1980, PDC deposits show that slope and irregular topography strongly influence PDC flow path, dynamics, criticality (for example, supercritical versus subcritical), carrying capacity, and erosive capacity. However, the influence of these conditions on ultimate flow runout and damage potential warrants further exploration through the combination of field, experimental, and numerical approaches.

Washington↗

USGS2017

Thermal and petrologic constraints on lower crustal melt accumulation under the Salton Sea Geothermal Field

In the Salton Sea region of southern California (USA), concurrent magmatism, extension, subsidence, and sedimentation over the past 0.5 to 1.0 Ma have led to the creation of the Salton Sea Geothermal Field (SSGF)—the second largest and hottest geothermal system in the continental United States—and the small-volume rhyolite eruptions that created the Salton Buttes. In this study, we determine the flux of mantle-derived basaltic magma that would be required to produce the elevated average heat flow and sustain the magmatic roots of rhyolite volcanism observed at the surface of the Salton Sea region. We use a 2D thermal model to show that a lower-crustal, partially molten mush containing < 20 – 40 % "> <20–40% interstitial melt develops over a ∼ 10 5 "> ∼105 -yr timescale for basalt fluxes of 0.008 to 0.010 m 3 / m 2 / yr "> 0.010m3/m2/yr (∼0.0008 to ∼0.001 km / 3 yr "> /3yr injection rate) given extension rates at or below the current value of ∼0.01 m/yr ( Brothers et al., 2009 ). These regions of partial melt are a natural consequence of a thermal regime that scales with average surface heat flow in the Salton Trough, and are consistent with seismic observations. Our results indicate limited melting and assimilation of pre-existing rocks in the lower crust. Instead, we find that basalt fractionation in the lower crust produces derivative melts of andesitic to dacitic composition. Such melts are then expected to ascend and accumulate in the upper crust, where they further evolve to give rise to small-volume rhyolite eruptions (Salton Buttes) and fuel local spikes in surface heat flux as currently seen in the SSGF. Such upper crustal magma evolution, with limited assimilation of hydrothermally altered material, is required to explain the slight decrease in δ 18 O "> δ18O values of zircons (and melts) that have been measured in these rhyolites.

California↗

USGS2013

Prolonged magmatic activity on Mars inferred from the detection of felsic rocks

Rocks dominated by the silicate minerals quartz and feldspar are abundant in Earth’s upper continental crust. Yet felsic rocks have not been widely identified on Mars, a planet that seems to lack plate tectonics and the associated magmatic processes that can produce evolved siliceous melts on Earth. If Mars once had a feldspar-rich crust that crystallized from an early magma ocean such as that on the Moon, erosion, sedimentation and volcanism have erased any clear surface evidence for widespread felsic materials. Here we report near-infrared spectral evidence from the Compact Reconnaissance Imaging Spectrometer for Mars onboard the Mars Reconnaissance Orbiter for felsic rocks in three geographically disparate locations on Mars. Spectral characteristics resemble those of feldspar-rich lunar anorthosites, but are accompanied by secondary alteration products (clay minerals). Thermodynamic phase equilibrium calculations demonstrate that fractional crystallization of magma compositionally similar to volcanic flows near one of the detection sites can yield residual melts with compositions consistent with our observations. In addition to an origin by significant magma evolution, the presence of felsic materials could also be explained by feldspar enrichment by fluvial weathering processes. Our finding of felsic materials in several locations on Mars suggests that similar observations by the Curiosity rover in Gale crater may be more widely applicable across the planet.

Nature Geoscience↗