Geology topics

Michael Manga

Publications and source records attributed to Michael Manga.

At least 19 recordsLinked to original sources

USGS2026

Water volumes, heat flow, and solute discharge from Old Faithful Geyser eruptions, Yellowstone National Park, USA

The iconic Old Faithful Geyser in Yellowstone National Park, USA, has attracted a significant amount of research because of the relative regularity and impressive size of its eruptions. Numerous studies have included observations, measurements, and analyses that informed models of geyser eruptions. However, fundamental quantities, including the associated mass and heat discharged, remain poorly constrained. In April 2025 we measured the volume of water from 45 Old Faithful Geyser eruptions using a portable flume in an outflow channel and specific conductance measurements in the Firehole River. We used high-speed video to perform velocimetry, measured changes in water chemistry to calculate the volume of water evaporated along the outflow channels, and used thermodynamic calculations to estimate the volume of water erupted as steam and to quantify the geyser's heat output. The calculated average volume of water discharged by Old Faithful Geyser in each eruption is 27.9 ± 9.4 m 3 , with no relation between eruption volume and the length of the preceding eruption interval. Video analysis of the eruption's liquid-dominated phase yields similar volumes of 21–30 m 3 . The calculated heat flow from the geyser is 2.2–2.4 MW and the average annual discharge of chloride, fluoride, and arsenic are 63 tons, 3.9 tons, and 241 kg, respectively. Average annual silica deposition rate on the geyser cone and along the outflow channels is 7 tons. This study provides a methodology for future studies at geysers worldwide and a baseline for monitoring future activity changes at Old Faithful.

Wyoming↗

USGS2025

Memory and jamming in fault zone sediments

Many subsurface processes involve transitions in granular material states, from arrested to creeping to flowing. Experiments and frameworks for idealized systems reveal that granular fabrics develop during shearing, co-evolve with applied stress, and govern such transitions. We use microtomography to test whether fabrics at two San Andreas fault sites reflect slip history and whether idealized frameworks extend to nature. Near-surface sediments within the fault zone transition between deformation patterns over the seismic cycle, including bulk/localized grain re-arrangements, individual grain fracturing, and localized zones of fracturing. Aseismic and co-seismic shearing produce distinct preferred grain orientations. Co-seismic fabrics can be preserved after centuries of aseismic strain, aseismic fabrics may be overprinted, and grain size and coordination number influence the fabrics. Idealized frameworks, namely anisotropic critical state theory, frictional jamming, and material memory, can explain our observations, and fault zone sediments likely undergo cycles of memory creation and erasure that influence rigidity spatiotemporally.

California↗

USGS2024

The evolution of a young ocean within Mimas

The fractured, young surfaces on confirmed ocean worlds such as Europa and Enceladus suggest that ocean-bearing moons with relatively thin overlying ice shells should be easy to identify. Hence, the discovery that Mimas’ rotation state is best explained by an internal ocean seems challenging to reconcile with its heavily cratered surface. Previous studies have shown that an internal ocean is compatible with Mimas’ geology as long as the ice shell has been thinning throughout the surface age. This scenario has yet to be placed into context within Mimas’ thermal-orbital history, in particular, the link between tidal dissipation and orbit circularization. Here, we model Mimas’ coupled thermal-orbital evolution, and the implications of an emerging ocean on Mimas’ geology, to determine whether a thinning ice shell is feasible. We find that the ice shell can thin – and the ocean grow – even as the eccentricity decays due to tidal dissipation as long as the ocean formed within the past 10 – 15 million years and the onset of melting occurred when Mimas’ eccentricity was between 2.5 and 3 times the present-day value. As Mimas’ eccentricity continues to decay, the ocean will likely enter an epoch of freezing, leading to fractures within the ice shell and eruptions of ocean material at the surface, before freezing completely. These results suggest that Mimas, Enceladus, and Tethys may represent three points along a continuum of ocean initiation, growth, and loss.

Earth and Planetary Science Letters↗

USGS2023

Pre-eruptive outgassing and pressurization, and post-fragmentation bubble nucleation, recorded by vesicles in breadcrust bombs from vulcanian activity at Guagua Pichincha Volcano, Ecuador

Breadcrust bombs formed during Vulcanian eruptions are assumed to originate from the shallow plug or dome. Their rim to core texture reflects the competition between cooling and degassing timescales, which results in a dense crust with isolated vesicles contrasting with a highly vesicular vesicle network in the interior. Due to relatively fast quenching, the crust can shed light on pre- and syn-eruptive conditions prior to or during fragmentation, whereas the interior allows us to explore post-fragmentation vesiculation. Investigation of pre- to post-fragmentation processes in breadcrust bombs from the 1999 Vulcanian activity at Guagua Pichincha, Ecuador, via 2D and 3D textural analysis reveals a complex vesiculation history, with multiple, spatially localized nucleation and growth events. Large vesicles (Type 1), present in low number density in the crust, are interpreted as pre-eruptive bubbles formed by outgassing and collapse of a permeable bubble network during ascent or stalling in the plug. Haloes of small, syn-fragmentation vesicles (Type 2), distributed about large vesicles, are formed by pressurization and enrichment of volatiles in these haloes. The nature of the pressurization process in the plug is discussed in light of seismicity and ground deformation signals, and previous textural and chemical studies. A third population (Type 3) of post-fragmentation small vesicles appears in the interior of the bomb, and growth and coalescence of Type 2 and 3 vesicles causes the transition from isolated to interconnected bubble network in the interior. We model the evolution of viscosity, bubble growth rate, diffusion timescales, bubble radius and porosity during fragmentation and cooling. These models reveal that thermal quenching dominates in the crust whereas the interior undergoes a viscosity quench caused by degassing, and that the transition from crust to interior corresponds to the onset of percolation and development of permeability in the bubble network.

Guagua Pichincha Volcano↗

USGS2023

The relation between decadal droughts and eruptions of Steamboat Geyser in Yellowstone National Park, USA

In the past century, most eruptions of Steamboat Geyser in Yellowstone National Park's Norris Geyser Basin were mainly clustered in three episodes: 1961–1969, 1982–1984, and ongoing since 2018. These eruptive episodes resulted in extensive disturbance to surrounding trees. To characterize tree response over time as an indicator of geyser activity adjustments to climate variability, aerial and ground images were analyzed to document changes in tree coverage around the geyser since 1954. Radiocarbon dating of silicified tree remnants from within 14 m of the geyser vent was used to examine geyser response to possible variations in decadal to centennial precipitation patterns. We searched for atypical or absent growth rings in cores from live trees in years associated with large geyser eruptions. Photographs indicate that active eruptive phases have adversely affected trees up to 30 m from the vent, primarily in the dominant downwind direction. Radiocarbon dates indicate that the geyser formed before 1878, in contrast to the birthdate reported in historical documents. Further, the ages of the silicified trees cluster within three episodes that are temporally correlated with periods of relative drought in the Yellowstone region during the 15th–17th centuries. The discontinuous growth of trees around the geyser suggests that changes in eruptive patterns occur in response to decadal to multidecadal droughts. This inference is supported by the lack of silicified specimens with more than 20 annual rings and by the existence of atypical or missing rings in live trees during periods of extended geyser activity.

Wyoming↗

USGS2023

Coordinating science during an eruption: Lessons from the 2020–2021 Kīlauea volcanic eruption

Data collected during well-observed eruptions can lead to dramatic increases in our understanding of volcanic processes. However, the necessary prioritization of public safety and hazard mitigation during a crisis means that scientific opportunities may be sacrificed. Thus, maximizing the scientific gains from eruptions requires improved planning and coordinating science activities among governmental organizations and academia before and during volcanic eruptions. One tool to facilitate this coordination is a Scientific Advisory Committee (SAC). In the USA, the Community Network for Volcanic Eruption Response (CONVERSE) has been developing and testing this concept during workshops and scenario-based activities. The December 2020 eruption of Kīlauea volcano, Hawaii, provided an opportunity to test and refine this model in real-time and in a real-world setting. We present here the working model of a SAC developed during this eruption. Successes of the Kīlauea SAC (K-SAC) included broadening the pool of scientists involved in eruption response and developing and codifying procedures that may form the basis of operation for future SACs. Challenges encountered by the K-SAC included a process of review and facilitation of research proposals that was too slow to include outside participation in the early parts of the eruption and a decision process that fell on a small number of individuals at the responding volcano observatory. Possible ways to address these challenges include (1) supporting community-building activities between eruptions that make connections among scientists within and outside formal observatories, (2) identifying key science questions and pre-planning science activities, which would facilitate more rapid implementation across a broader scientific group, and (3) continued dialog among observatory scientists, emergency responders, and non-observatory scientists about the role of SACs. The SAC model holds promise to become an integral part of future efforts, leading in the short and longer term to more effective hazard response and greater scientific discovery and understanding.

Hawaii↗

USGS2022

Violent groundwater eruption triggered by a distant earthquake

It is now well established that earthquakes cause various hydrogeological responses at distances thousands of kilometers from the epicenter. What remains unexplained is the large amplitude and intensity of some responses. Following the 2004 Mw 9.1 Sumatra earthquake, groundwater 3,200 km from the epicenter erupted violently from a well and formed a water fountain reaching a height exceeding 60 m. We model the relevant processes by combining tidal analysis of groundwater level with numerical simulations using a two-dimensional finite-element model. We suggest that the eruption resulted from a combination of factors, including a rapid increase of crustal permeability and runaway CO 2 exsolution and bubble nucleation induced by the passage of seismic waves. Our results may have implications for some engineering applications such as oil production and CO 2 sequestration, and the eruption of hydrothermal features such as geysers.

Geophysical Research Letters↗

USGS2022

Outgassing through magmatic fractures enables effusive eruption of silicic magma

Several mechanisms have been proposed to allow highly viscous silicic magma to outgas efficiently enough to erupt effusively. There is increasing evidence that challenges the classic foam-collapse model in which gas escapes through permeable bubble networks, and instead suggests that magmatic fracturing and/or accompanying localized fragmentation and welding within the conduit play an important role in outgassing. The 2011–2012 eruption at Cordón Caulle volcano, Chile, provides direct observations of the role of magmatic fractures. This eruption exhibited a months-long hybrid phase, in which rhyolitic lava extrusion was accompanied by vigorous gas-and-tephra venting through fractures in the lava dome surface. Some of these fractures were preserved as tuffisites (tephra-filled veins) in erupted lava and bombs. We integrate constraints from petrologic analyses of erupted products and video analyses of gas-and-tephra venting to construct a model for magma ascent in a conduit. The one-dimensional, two-phase, steady-state model considers outgassing through deforming permeable bubble networks, magmatic fractures, and adjacent wall rock. Simulations for a range of plausible magma ascent conditions indicate that the eruption of low-porosity lava observed at Cordón Caulle volcano occurs because of significant gas flux through fracture networks in the upper conduit. This modeling emphasizes the important role that outgassing through magmatic fractures plays in sustaining effusive or hybrid eruptions of silicic magma and in facilitating explosive-effusive transitions.

Puyehue Cordon Caulle volcanic complex↗

USGS2021

Brittle fragmentation by rapid gas separation in a Hawaiian fountain

Brittle fragmentation, generating small pyroclasts from magma, is a key process determining eruptive style. How low-viscosity magma fragments within a rising fountain in a brittle manner, however, is not well understood. Here we describe a fragmentation process in Hawaiian fountains on the basis of observations from the 2018 lower East Rift Zone eruption of Kīlauea Volcano, Hawai’i. The dominant fragmentation mechanism is inertia driven and produces a population of large fluidal pyroclasts. However, when sufficient volcanic gas is released in the fountain, a subpopulation of smaller and more vesicular pyroclasts is generated and entrained into the gas-dominant convective plume. The size distribution of these pyroclasts is similar to that of brittlely fragmented solid materials. The erupted high-vesicularity pyroclasts sometimes preserve a deformed shape. These observations suggest that late-stage rapid expansion lowers the gas temperature adiabatically and cools the outer surface of liquid pyroclasts below the glass transition temperature. The rigid crust fragments as the hot interior attempts to expand due to further volatile diffusion from the melt into bubbles. Adiabatic expansion of volcanic gas occurs in all eruptions. Brittle fragmentation induced by rapid adiabatic cooling may be a widespread process, although of varying importance, in explosive eruptions.

Hawaii↗

USGS2021

The 2018 reawakening and eruption dynamics of Steamboat Geyser, the world’s tallest active geyser

Steamboat Geyser in Yellowstone National Park’s Norris Geyser Basin began a prolific sequence of eruptions in March 2018 after 34 y of sporadic activity. We analyze a wide range of datasets to explore triggering mechanisms for Steamboat’s reactivation and controls on eruption intervals and height. Prior to Steamboat’s renewed activity, Norris Geyser Basin experienced uplift, a slight increase in radiant temperature, and increased regional seismicity, which may indicate that magmatic processes promoted reactivation. However, because the geothermal reservoir temperature did not change, no other dormant geysers became active, and previous periods with greater seismic moment release did not reawaken Steamboat, the reason for reactivation remains ambiguous. Eruption intervals since 2018 (3.16 to 35.45 d) modulate seasonally, with shorter intervals in the summer. Abnormally long intervals coincide with weakening of a shallow seismic source in the geyser basin’s hydrothermal system. We find no relation between interval and erupted volume, implying unsteady heat and mass discharge. Finally, using data from geysers worldwide, we find a correlation between eruption height and inferred depth to the shallow reservoir supplying water to eruptions. Steamboat is taller because water is stored deeper there than at other geysers, and, hence, more energy is available to power the eruptions.

Wyoming↗

USGS2021

The structure and volume of large geysers in Yellowstone National Park, USA and the mineralogy and chemistry of their silica sinter deposits

Siliceous sinter is formed by biogenic and abiogenic opal deposition around hot springs and geysers. Using Structure-from-Motion photogrammetry we generated three-dimensional models of Giant and Castle Geysers from the Upper Geyser Basin of Yellowstone National Park. We use these models to calculate an approximate mass of sinter for each (~2 and ~ 5 kton, respectively) and estimate a range of plausible long-term deposition rates for Castle Geyser (470 to 940 kg·yr −1 ). We estimate ~2% of the silica discharged from Castle Geyser is deposited as sinter in the cone and proximal terraces. We collected 15 sinter samples following the stratigraphy of each geyser from an older terrace to a younger cone and examined them using a variety of analytical methods. We find that young opaline sinter with a water content of <12 wt% (from loss on ignition) contains higher concentrations of major and trace elements, notably As, Sb, Rb, Ga and Cs, relative to older dehydrated sinter. Rare earth element (REE) concentrations in sinter are 2–3 orders of magnitude higher than in the thermal water from which they are deposited. Sinter deposits are enriched in light REE, Gd and Yb when normalized to concentrations in thermal water and enriched in Eu, Tm, and Yb when normalized to the underlying rhyolite. Sinter samples with the highest REE concentrations are also enriched in organic material, implying either microbial uptake of REE, or that organic molecules are efficient ligands that form metal complexes.

Wyoming↗

USGS2021

Why study geysers?

Scientific research for more than two centuries has improved our understanding of Earth’s geysers. This knowledge provides insights into volcanic processes, the origin and environmental limits of life on Earth and potentially Mars, and on geysers on icy outer solar system satellites. Continued scientific research will help us understand and protect these natural wonders that attract millions of tourists annually.

Eos, American Geophysical Union↗

USGS2020

Origin and properties of hydrothermal tremor at Lone Star Geyser, Yellowstone National Park, USA

Geysers are rare geologic features that intermittently discharge liquid water and steam driven by heating and decompression boiling. The cause of variability in eruptive styles and the associated seismic signals are not well understood. Data collected from five broadband seismometers at Lone Star Geyser, Yellowstone National Park are used to determine the properties, location, and temporal patterns of hydrothermal tremor. The tremor is harmonic at some stages of the eruption cycle and is caused by near‐periodic repetition of discrete seismic events. Using the polarization of ground motion, we identify the location of tremor sources throughout several eruption cycles. During preplay episodes (smaller eruptions preceding the more vigorous major eruption), tremor occurs at depths of 7–10 m and is laterally offset from the geyser's cone by ~5 m. At the onset of the main eruption, tremor sources migrate laterally and become shallower. As the eruption progresses, tremor sources migrate along the same path but in the opposite direction, ending where preplay tremor originates. The upward and then downward migration of tremor sources during eruptions are consistent with warming of the conduit followed by evacuation of water during the main eruption. We identify systematic relations among the two types of preplays, discharge, and the main eruption. A point‐source moment tensor fit to low‐frequency waveforms of an individual tremor event using half‐space velocity models indicates average V S ≳ 0.8 km/s, source depths ~4–20 m, and moment tensors with primarily positive isotropic and compensated linear vector dipole moments.

Wyoming↗

USGS2020

Yellowstone's Old Faithful Geyser shut down by a severe 13th century drought

To characterize eruption activity of the iconic Old Faithful Geyser in Yellowstone National Park over past centuries, we obtained 41 new radiocarbon dates of mineralized wood preserved in the mound of silica that precipitated from erupted waters. Trees do not grow on active geyser mounds, implying that trees grew on the Old Faithful Geyser mound during a protracted period of eruption quiescence. Rooted stumps and root crowns located on higher parts of the mound are evidence that at the time of tree growth, the geyser mound closely resembled its current appearance. The range of calibrated radiocarbon dates (1233–1362 CE) is coincident with a series of severe multidecadal regional droughts toward the end of the Medieval Climate Anomaly, prior to the onset of the Little Ice Age. Climate models project increasingly severe droughts by mid‐21st century, suggesting that geyser eruptions could become less frequent or completely cease.

Wyoming↗

USGS2020

Radiocarbon dating of silica sinter and postglacial hydrothermal activity in the El Tatio geyser field

The El Tatio geothermal field in the Chilean Altiplano contains hydrothermal silica sinter deposits overlaying glacial and volcanic units, providing an opportunity to constrain the timing of deglaciation and volcanic activity in an area with sparse absolute chronologies. We obtained 51 new radiocarbon ages and δ 13 C values on the organic material trapped in these sinter deposits. Based on the δ 13 C values, we exclude 29 samples for possible contamination with bacterial mats that incorporate old carbon. We infer that hydrothermal activity initiated ~27 ka ago and has been nearly continuous ever since. The ages of the oldest sinter deposits coincide with ages of moraines that stabilized after the most recent deglaciation. Whereas late Pleistocene sinters are broadly distributed in the field, Holocene deposits are found around active hydrothermal features. Although recent volcanism is absent in the vicinity of El Tatio, persistent hydrothermal discharge implies a long‐lived magmatic heat source.

El Tatio Geothermal Field↗

USGS2020

Dating silica sinter (geyserite): A cautionary tale

We describe a new effort to date hydrothermal silica sinter deposits (geyserite) from the Upper Geyser Basin of Yellowstone National Park using 14C of co-deposited organic matter, U-series and cosmogenic 10Be methods. A majority of the samples were collected from stratigraphic sections, mainly at Riverside, Giant, and Castle Geysers. Ages obtained from 41 14C analyses range from modern to 12.1 cal ka BP. Nearly all the 14C ages show inconsistencies with their stratigraphic positions, and several replicate 14C analyses from the same sample result in significantly different ages. The δ13C values of the organic material in the sinter range from -26.6‰ to -12.7‰. The more enriched values are attributed to microbial fixation of dissolved inorganic carbon (DIC), which has heavier δ13C values and is 14C-depleted relative to atmospheric CO2, leading to apparent older ages. U-series analyses on 4 samples yielded ages between 2.2 and 7.4 ka. Large 230Th/U age uncertainties in the sinter, due to low uranium concentrations along with elevated 232Th and associated initial 230Th, make these ages imprecise for use on Holocene deposits. A single cosmogenic 10Be exposure age of 596±18 ka is considerably older than the age of underlying rhyolite and is thus unreliable. This apparent old age results from contamination by meteoric 10Be trapped in the opal that overprints the very small amount of cosmogenic 10Be. By presenting the problems we encountered and discussing their probable cause, this paper highlights the difficulty in obtaining reliable, high-precision geochronological data necessary to use sinter deposits as paleoenvironmental and paleo-hydrothermal archives.

Wyoming↗

USGS2018

Stronger peak ground motion, beyond the threshold to initiate a response, does not lead to larger stream discharge responses to earthquakes

The impressive number of stream gauges in Chile, combined with a suite of past and recent large earthquakes, makes Chile a unique natural laboratory to study several streams that recorded responses to multiple seismic events. We document changes in discharge in eight streams in Chile following two or more large earthquakes. In all cases, discharge increases. Changes in discharge occur for peak ground velocities greater than about 7–11 cm/s. Above that threshold, the magnitude of both the increase in discharge and the total excess water do not increase with increasing peak ground velocities. While these observations are consistent with previous work in California, they conflict with lab experiments that show that the magnitude of permeability changes increases with increasing amplitude of ground motion. Instead, our study suggests that streamflow responses are binary.

Geophysical Research Letters↗

USGS2018

Hydrothermal discharge from the El Tatio basin, Atacama, Chile

El Tatio in northern Chile is one of the best-studied geothermal fields in South America. However, there remain open questions about the mass and energy budgets, water recharge rates and residence time in the subsurface, origin of dissolved solutes, and processes affecting the phase and chemical composition of groundwater and surface water. We measured and sampled surface manifestations of the geothermal system (geysers perpetual spouters, mud pools/volcanoes, and non-eruptive hot springs) and meteoric water. From the isotopic composition we infer that the thermal water has a meteoric origin that is different from the composition of local meteoric water. The absence of detectable tritium in thermal waters indicates that most of the recharge occurred pre-1950. Boiling and steam separation from the deep reservoir appear to be the main subsurface processes affecting the thermal fluids. A large amount of heat is lost to the atmosphere by evaporation from surface water and by steam emitted from erupting geysers. Using the chloride inventory method, we estimate thermal water discharge to be 218 to 234 L/s, and the advective heat flow to be 120 to 170 MW.

El Tatio Basin↗