Airfall volume of the 15 January 2022 eruption of Hunga volcano estimated from ocean color changes.

On 15 January 2022, Hunga volcano erupted, creating an extensive and high-reaching umbrella cloud over the open ocean, hindering traditional isopach mapping and fallout volume estimation. In MODIS satellite imagery, ocean surface water was discolored around Hunga following the eruption, which we attribute to ash fallout from the umbrella cloud. By relating intensity of ocean discoloration to fall deposit thicknesses in the Kingdom of Tonga, we develop a methodology for estimating airfall volume over the open ocean. Ash thickness measurements from 41 locations are used to fit a linear relationship between ash thickness and ocean reflectance. This produces a minimum airfall volume estimate of 1.8-0.4+0.3 km3. The whole eruption produced > 6.3 km3 of uncompacted pyroclastic material on the seafloor and a caldera volume change of 6 km3 DRE. Our fall estimates are consistent with the interpretation that most of the seafloor deposits were emplaced by gravity currents rather than fall deposits. Our proposed method does not account for the largest grain sizes, so is thus a minimum estimate. However, this new ocean-discoloration method provides an airfall volume estimate consistent with other independent measures of the plume and is thus effective for rapidly estimating fallout volumes in future volcanic eruptions over oceans. Supplementary Information: The online version contains supplementary material available at 10.1007/s00445-024-01744-6.

Autotrophic biofilms sustained by deeply sourced groundwater host diverse bacteria implicated in sulfur and hydrogen metabolism.

BACKGROUND: Biofilms in sulfide-rich springs present intricate microbial communities that play pivotal roles in biogeochemical cycling. We studied chemoautotrophically based biofilms that host diverse CPR bacteria and grow in sulfide-rich springs to investigate microbial controls on biogeochemical cycling. RESULTS: Sulfide springs biofilms were investigated using bulk geochemical analysis, genome-resolved metagenomics, and scanning transmission X-ray microscopy (STXM) at room temperature and 87 K. Chemolithotrophic sulfur-oxidizing bacteria, including Thiothrix and Beggiatoa, dominate the biofilms, which also contain CPR Gracilibacteria, Absconditabacteria, Saccharibacteria, Peregrinibacteria, Berkelbacteria, Microgenomates, and Parcubacteria. STXM imaging revealed ultra-small cells near the surfaces of filamentous bacteria that may be CPR bacterial episymbionts. STXM and NEXAFS spectroscopy at carbon K and sulfur L2,3 edges show that filamentous bacteria contain protein-encapsulated spherical elemental sulfur granules, indicating that they are sulfur oxidizers, likely Thiothrix. Berkelbacteria and Moranbacteria in the same biofilm sample are predicted to have a novel electron bifurcating group 3b [NiFe]-hydrogenase, putatively a sulfhydrogenase, potentially linked to sulfur metabolism via redox cofactors. This complex could potentially contribute to symbioses, for example, with sulfur-oxidizing bacteria such as Thiothrix that is based on cryptic sulfur cycling. One Doudnabacteria genome encodes adjacent sulfur dioxygenase and rhodanese genes that may convert thiosulfate to sulfite. We find similar conserved genomic architecture associated with CPR bacteria from other sulfur-rich subsurface ecosystems. CONCLUSIONS: Our combined metagenomic, geochemical, spectromicroscopic, and structural bioinformatics analyses of biofilms growing in sulfide-rich springs revealed consortia that contain CPR bacteria and sulfur-oxidizing Proteobacteria, including Thiothrix, and bacteria from a new family within Beggiatoales. We infer roles for CPR bacteria in sulfur and hydrogen cycling. Video Abstract.

Coordinating science during an eruption: lessons from the 2020-2021 Kilauea 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 Kilauea 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 Kilauea 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.

The History of Water in Martian Magmas From Thorium Maps.

Water inventories in Martian magmas are poorly constrained. Meteorite-based estimates range widely, from 102 to >104 ppm H2O, and are likely variably influenced by degassing. Orbital measurements of H primarily reflect water cycled and stored in the regolith. Like water, Th behaves incompatibly during mantle melting, but unlike water Th is not prone to degassing and is relatively immobile during aqueous alteration at low temperature. We employ Th as a proxy for original, mantle-derived H2O in Martian magmas. We use regional maps of Th from Mars Odyssey to assess variations in magmatic water across major volcanic provinces and through time. We infer that Hesperian and Amazonian magmas had ∼100-3,000 ppm H2O, in the lower range of previous estimates. The implied cumulative outgassing since the Hesperian, equivalent to a global H2O layer ∼1-40 m deep, agrees with Mars' present-day surface and near-surface water inventory and estimates of sequestration and loss rates.

Microstructural differences between naturally-deposited and laboratory beach sands.

The orientation of, and contacts between, grains of sand reflect the processes that deposit the sands. Grain orientation and contact geometry also influence mechanical properties. Quantifying and understanding sand microstructure thus provide an opportunity to understand depositional processes better and connect microstructure and macroscopic properties. Using x-ray computed microtomography, we compare the microstructure of naturally-deposited beach sands and laboratory sands created by air pluviation in which samples are formed by raining sand grains into a container. We find that naturally-deposited sands have a narrower distribution of coordination number (i.e., the number of grains in contact) and a broader distribution of grain orientations than pluviated sands. The naturally-deposited sand grains orient inclined to the horizontal, and the pluviated sand grains orient horizontally. We explain the microstructural differences between the two different depositional methods by flowing water at beaches that re-positions and reorients grains initially deposited in unstable grain configurations.

Widespread deep seismicity in the Delaware Basin, Texas, is mainly driven by shallow wastewater injection.

Industrial activity away from plate boundaries can induce earthquakes and has evolved into a global issue. Much of the induced seismicity in the United States' midcontinent is attributed to a direct pressure increase from deep wastewater disposal. This mechanism is not applicable where deep basement faults are hydraulically isolated from shallow injection aquifers, leading to a debate about the mechanisms for induced seismicity. Here, we compile industrial, seismic, geodetic, and geological data within the Delaware Basin, western Texas, and calculate stress and pressure changes at seismogenic depth using a coupled poroelastic model. We show that the widespread deep seismicity is mainly driven by shallow wastewater injection through the transmission of poroelastic stresses assuming that unfractured shales are hydraulic barriers over decadal time scales. A zone of seismic quiescence to the north, where injection-induced stress changes would promote seismicity, suggests a regional tectonic control on the occurrence of induced earthquakes. Comparing the poroelastic responses from injection and extraction operations, we find that the basement stress is most sensitive to shallow reservoir hydrogeological parameters, particularly hydraulic diffusivity. These results demonstrate that intraplate seismicity can be caused by shallow human activities that poroelastically perturb stresses at hydraulically isolated seismogenic depths, with impacts on seismicity that are preconditioned by regional tectonics.

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.

Stable isotopes show that earthquakes enhance permeability and release water from mountains.

Hydrogeological properties can change in response to large crustal earthquakes. In particular, permeability can increase leading to coseismic changes in groundwater level and flow. These processes, however, have not been well-characterized at regional scales because of the lack of datasets to describe water provenances before and after earthquakes. Here we use a large data set of water stable isotope ratios (n = 1150) to show that newly formed rupture systems crosscut surrounding mountain aquifers, leading to water release that causes groundwater levels to rise (~11 m) in down-gradient aquifers after the 2016 Mw 7.0 Kumamoto earthquake. Neither vertical infiltration of soil water nor the upwelling of deep fluids was the major cause of the observed water level rise. As the Kumamoto setting is representative of volcanic aquifer systems at convergent margins where seismotectonic activity is common, our observations and proposed model should apply more broadly.

Pore-pressure diffusion, enhanced by poroelastic stresses, controls induced seismicity in Oklahoma.

Induced seismicity linked to geothermal resource exploitation, hydraulic fracturing, and wastewater disposal is evolving into a global issue because of the increasing energy demand. Moderate to large induced earthquakes, causing widespread hazards, are often related to fluid injection into deep permeable formations that are hydraulically connected to the underlying crystalline basement. Using injection data combined with a physics-based linear poroelastic model and rate-and-state friction law, we compute the changes in crustal stress and seismicity rate in Oklahoma. This model can be used to assess earthquake potential on specific fault segments. The regional magnitude-time distribution of the observed magnitude (M) 3+ earthquakes during 2008-2017 is reproducible and is the same for the 2 optimal, conjugate fault orientations suggested for Oklahoma. At the regional scale, the timing of predicted seismicity rate, as opposed to its pattern and amplitude, is insensitive to hydrogeological and nucleation parameters in Oklahoma. Poroelastic stress changes alone have a small effect on the seismic hazard. However, their addition to pore-pressure changes can increase the seismicity rate by 6-fold and 2-fold for central and western Oklahoma, respectively. The injection-rate reduction in 2016 mitigates the exceedance probability of M5.0 by 22% in western Oklahoma, while that of central Oklahoma remains unchanged. A hypothetical injection shut-in in April 2017 causes the earthquake probability to approach its background level by ∼2025. We conclude that stress perturbation on prestressed faults due to pore-pressure diffusion, enhanced by poroelastic effects, is the primary driver of the induced earthquakes in Oklahoma.

Controls on explosive-effusive volcanic eruption styles.

One of the biggest challenges in volcanic hazard assessment is to understand how and why eruptive style changes within the same eruptive period or even from one eruption to the next at a given volcano. This review evaluates the competing processes that lead to explosive and effusive eruptions of silicic magmas. Eruptive style depends on a set of feedback involving interrelated magmatic properties and processes. Foremost of these are magma viscosity, gas loss and external properties such as conduit geometry. Ultimately, these parameters control the speed at which magmas ascend, decompress and outgas en route to the surface, and thus determine eruptive style and evolution.

Timing of oceans on Mars from shoreline deformation.

Widespread evidence points to the existence of an ancient Martian ocean. Most compelling are the putative ancient shorelines in the northern plains. However, these shorelines fail to follow an equipotential surface, and this has been used to challenge the notion that they formed via an early ocean and hence to question the existence of such an ocean. The shorelines' deviation from a constant elevation can be explained by true polar wander occurring after the formation of Tharsis, a volcanic province that dominates the gravity and topography of Mars. However, surface loading from the oceans can drive polar wander only if Tharsis formed far from the equator, and most evidence indicates that Tharsis formed near the equator, meaning that there is no current explanation for the shorelines' deviation from an equipotential that is consistent with our geophysical understanding of Mars. Here we show that variations in shoreline topography can be explained by deformation caused by the emplacement of Tharsis. We find that the shorelines must have formed before and during the emplacement of Tharsis, instead of afterwards, as previously assumed. Our results imply that oceans on Mars formed early, concurrent with the valley networks, and point to a close relationship between the evolution of oceans on Mars and the initiation and decline of Tharsis volcanism, with broad implications for the geology, hydrological cycle and climate of early Mars.

The largest deep-ocean silicic volcanic eruption of the past century.

The 2012 submarine eruption of Havre volcano in the Kermadec arc, New Zealand, is the largest deep-ocean eruption in history and one of very few recorded submarine eruptions involving rhyolite magma. It was recognized from a gigantic 400-km2 pumice raft seen in satellite imagery, but the complexity of this event was concealed beneath the sea surface. Mapping, observations, and sampling by submersibles have provided an exceptionally high fidelity record of the seafloor products, which included lava sourced from 14 vents at water depths of 900 to 1220 m, and fragmental deposits including giant pumice clasts up to 9 m in diameter. Most (>75%) of the total erupted volume was partitioned into the pumice raft and transported far from the volcano. The geological record on submarine volcanic edifices in volcanic arcs does not faithfully archive eruption size or magma production.

Monitoring reservoir response to earthquakes and fluid extraction, Salton Sea geothermal field, California.

Continuous monitoring of in situ reservoir responses to stress transients provides insights into the evolution of geothermal reservoirs. By exploiting the stress dependence of seismic velocity changes, we investigate the temporal evolution of the reservoir stress state of the Salton Sea geothermal field (SSGF), California. We find that the SSGF experienced a number of sudden velocity reductions (~0.035 to 0.25%) that are most likely caused by openings of fractures due to dynamic stress transients (as small as 0.08 MPa and up to 0.45 MPa) from local and regional earthquakes. Depths of velocity changes are estimated to be about 0.5 to 1.5 km, similar to the depths of the injection and production wells. We derive an empirical in situ stress sensitivity of seismic velocity changes by relating velocity changes to dynamic stresses. We also observe systematic velocity reductions (0.04 to 0.05%) during earthquake swarms in mid-November 2009 and late-December 2010. On the basis of volumetric static and dynamic stress changes, the expected velocity reductions from the largest earthquakes with magnitude ranging from 3 to 4 in these swarms are less than 0.02%, which suggests that these earthquakes are likely not responsible for the velocity changes observed during the swarms. Instead, we argue that velocity reductions may have been induced by poroelastic opening of fractures due to aseismic deformation. We also observe a long-term velocity increase (~0.04%/year) that is most likely due to poroelastic contraction caused by the geothermal production. Our observations demonstrate that seismic interferometry provides insights into in situ reservoir response to stress changes.

Surface uplift and time-dependent seismic hazard due to fluid injection in eastern Texas.

Observations that unequivocally link seismicity and wastewater injection are scarce. Here we show that wastewater injection in eastern Texas causes uplift, detectable in radar interferometric data up to >8 kilometers from the wells. Using measurements of uplift, reported injection data, and a poroelastic model, we computed the crustal strain and pore pressure. We infer that an increase of >1 megapascal in pore pressure in rocks with low compressibility triggers earthquakes, including the 4.8-moment magnitude event that occurred on 17 May 2012, the largest earthquake recorded in eastern Texas. Seismic activity increased even while injection rates declined, owing to diffusion of pore pressure from earlier periods with higher injection rates. Induced seismicity potential is suppressed where tight confining formations prevent pore pressure from propagating into crystalline basement rocks.

New streams and springs after the 2014 Mw6.0 South Napa earthquake.

Many streams and springs, which were dry or nearly dry before the 2014 Mw6.0 South Napa earthquake, started to flow after the earthquake. A United States Geological Survey stream gauge also registered a coseismic increase in discharge. Public interest was heightened by a state of extreme drought in California. Since the new flows were not contaminated by pre-existing surface water, their composition allowed unambiguous identification of their origin. Following the earthquake we repeatedly surveyed the new flows, collecting data to test hypotheses about their origin. We show that the new flows originated from groundwater in nearby mountains released by the earthquake. The estimated total amount of new water is ∼ 10(6) m(3), about 1/40 of the annual water use in the Napa-Sonoma area. Our model also makes a testable prediction of a post-seismic decrease of seismic velocity in the shallow crust of the affected region.

Formation of Box Canyon, Idaho, by megaflood: implications for seepage erosion on Earth and Mars.

Amphitheater-headed canyons have been used as diagnostic indicators of erosion by groundwater seepage, which has important implications for landscape evolution on Earth and astrobiology on Mars. Of perhaps any canyon studied, Box Canyon, Idaho, most strongly meets the proposed morphologic criteria for groundwater sapping because it is incised into a basaltic plain with no drainage network upstream, and approximately 10 cubic meters per second of seepage emanates from its vertical headwall. However, sediment transport constraints, 4He and 14C dates, plunge pools, and scoured rock indicate that a megaflood (greater than 220 cubic meters per second) carved the canyon about 45,000 years ago. These results add to a growing recognition of Quaternary catastrophic flooding in the American northwest, and may imply that similar features on Mars also formed by floods rather than seepage erosion.

Evidence for an ancient martian ocean in the topography of deformed shorelines.

A suite of observations suggests that the northern plains of Mars, which cover nearly one third of the planet's surface, may once have contained an ocean. Perhaps the most provocative evidence for an ancient ocean is a set of surface features that ring the plains for thousands of kilometres and that have been interpreted as a series of palaeoshorelines of different age. It has been shown, however, that topographic profiles along the putative shorelines contain long-wavelength trends with amplitudes of up to several kilometres, and these trends have been taken as an argument against the martian shoreline (and ocean) hypothesis. Here we show that the long-wavelength topography of the shorelines is consistent with deformation caused by true polar wander--a change in the orientation of a planet with respect to its rotation pole--and that the inferred pole path has the geometry expected for a true polar wander event that postdates the formation of the massive Tharsis volcanic rise.

Seismological constraints on a possible plume root at the core-mantle boundary.

Recent seismological discoveries have indicated that the Earth's core-mantle boundary is far more complex than a simple boundary between the molten outer core and the silicate mantle. Instead, its structural complexities probably rival those of the Earth's crust. Some regions of the lowermost mantle have been observed to have seismic wave speed reductions of at least 10 per cent, which appear not to be global in extent. Here we present robust evidence for an 8.5-km-thick and approximately 50-km-wide pocket of dense, partially molten material at the core-mantle boundary east of Australia. Array analyses of an anomalous precursor to the reflected seismic wave ScP reveal compressional and shear-wave velocity reductions of 8 and 25 per cent, respectively, and a 10 per cent increase in density of the partially molten aggregate. Seismological data are incompatible with a basal layer composed of pure melt, and thus require a mechanism to prevent downward percolation of dense melt within the layer. This may be possible by trapping of melt by cumulus crystal growth following melt drainage from an anomalously hot overlying region of the lowermost mantle. This magmatic evolution and the resulting cumulate structure seem to be associated with overlying thermal instabilities, and thus may mark a root zone of an upwelling plume.

Explosive volcanism may not be an inevitable consequence of magma fragmentation.

The fragmentation of magma, containing abundant gas bubbles, is thought to be the defining characteristic of explosive eruptions. When viscous stresses associated with the growth of bubbles and the flow of the ascending magma exceed the strength of the melt, the magma breaks into disconnected fragments suspended within an expanding gas phase. Although repeated effusive and explosive eruptions for individual volcanoes are common, the dynamics governing the transition between explosive and effusive eruptions remain unclear. Magmas for both types of eruptions originate from sources with similar volatile content, yet effusive lavas erupt considerably more degassed than their explosive counterparts. One mechanism for degassing during magma ascent, consistent with observations, is the generation of intermittent permeable fracture networks generated by non-explosive fragmentation near the conduit walls. Here we show that such fragmentation can occur by viscous shear in both effusive and explosive eruptions. Moreover, we suggest that such fragmentation may be important for magma degassing and the inhibition of explosive behaviour. This implies that, contrary to conventional views, explosive volcanism is not an inevitable consequence of magma fragmentation.