Atmosphere injection of sea salts during large explosive submarine volcanic eruptions.

The 15 January 2022 submarine eruption at Hunga volcano was the most explosive volcanic eruption in 140 years. It involved exceptional magma and seawater interaction throughout the entire submarine caldera collapse. The submarine volcanic jet breached the sea surface and formed a subaerial eruptive plume that transported volcanic ash, gas, sea salts and seawater up to ~ 57 km, reaching into the mesosphere. We document high concentrations of sea salts in tephra (volcanic ash) collected shortly after deposition. We also discuss the potential climatic consequences of large-scale injection of salts into the upper atmosphere during submarine eruptions. Sodium chloride in these volcanic plumes can reach extreme concentrations, and dehalogenation of chlorides and bromides poses the risk of long-term atmospheric and weather impact. Salt content in rapidly collected tephra samples may also be used as a proxy to estimate the water:magma ratio during eruption, with implications for quantification of fragmentation efficiency in submarine breaching events. The balance between salt loading into the atmosphere versus deposition in ash aggregates is a key factor in understanding the atmospheric and climatic consequences of submarine eruptions.

Heat flows in rock cracks naturally optimize salt compositions for ribozymes.

Catalytic nucleic acids, such as ribozymes, are central to a variety of origin-of-life scenarios. Typically, they require elevated magnesium concentrations for folding and activity, but their function can be inhibited by high concentrations of monovalent salts. Here we show that geologically plausible high-sodium, low-magnesium solutions derived from leaching basalt (rock and remelted glass) inhibit ribozyme catalysis, but that this activity can be rescued by selective magnesium up-concentration by heat flow across rock fissures. In contrast to up-concentration by dehydration or freezing, this system is so far from equilibrium that it can actively alter the Mg:Na salt ratio to an extent that enables key ribozyme activities, such as self-replication and RNA extension, in otherwise challenging solution conditions. The principle demonstrated here is applicable to a broad range of salt concentrations and compositions, and, as such, highly relevant to various origin-of-life scenarios.

Earthquakes indicated magma viscosity during Kilauea's 2018 eruption.

Magma viscosity strongly controls the style (for example, explosive versus effusive) of a volcanic eruption and thus its hazard potential, but can only be measured during or after an eruption. The identification of precursors indicative of magma viscosity would enable forecasting of the eruption style and the scale of associated hazards1. The unanticipated May 2018 rift intrusion and eruption of Kilauea Volcano, Hawai'i2 displayed exceptional chemical and thermal variability in erupted lavas, leading to unpredictable effusion rates and explosivity. Here, using an integrated analysis of seismicity and magma rheology, we show that the orientation of fault-plane solutions (which indicate a fault's orientation and sense of movement) for earthquakes preceding and accompanying the 2018 eruption indicate a 90-degree local stress-field rotation from background, a phenomenon previously observed only at high-viscosity eruptions3, and never before at Kilauea4-8. Experimentally obtained viscosities for 2018 products and earlier lavas from the Pu'u 'O'o vents tightly constrain the viscosity threshold required for local stress-field reorientation. We argue that rotated fault-plane solutions in earthquake swarms at Kilauea and other volcanoes worldwide provide an early indication that unrest involves magma of heightened viscosity, and thus real-time monitoring of the orientations of fault-plane solutions could provide critical information about the style of an impending eruption. Furthermore, our results provide insight into the fundamental nature of coupled failure and flow in complex multiphase systems.

In situ observation of the percolation threshold in multiphase magma analogues.

Magmas vesiculate during ascent, producing complex interconnected pore networks, which can act as outgassing pathways and then deflate or compact to volcanic plugs. Similarly, in-conduit fragmentation events during dome-forming eruptions create open systems transiently, before welding causes pore sealing. The percolation threshold is the first-order transition between closed- and open-system degassing dynamics. Here, we use time-resolved, synchrotron-source X-ray tomography to image synthetic magmas that go through cycles of opening and closing, to constrain the percolation threshold Φ C at a range of melt crystallinity, viscosity and overpressure pertinent to shallow magma ascent. During vesiculation, we observed different percolative regimes for the same initial bulk crystallinity depending on melt viscosity and gas overpressure. At high viscosity (> 106 Pa s) and high overpressure (~ 1-4 MPa), we found that a brittle-viscous regime dominates in which brittle rupture allows system-spanning coalescence at a low percolation threshold (Φ C ~0.17) via the formation of fracture-like bubble chains. Percolation was followed by outgassing and bubble collapse causing densification and isolation of the bubble network, resulting in a hysteresis in the evolution of connectivity with porosity. At low melt viscosity and overpressure, we observed a viscous regime with much higher percolation threshold (Φ C > 0.37) due to spherical bubble growth and lower degree of crystal connection. Finally, our results also show that sintering of crystal-free and crystal-bearing magma analogues is characterised by low percolation thresholds (Φ C = 0.04 - 0.10). We conclude that the presence of crystals lowers the percolation threshold during vesiculation and may promote outgassing in shallow, crystal-rich magma at initial stages of Vulcanian and Strombolian eruptions.

A compositional tipping point governing the mobilization and eruption style of rhyolitic magma.

The most viscous volcanic melts and the largest explosive eruptions on our planet consist of calcalkaline rhyolites. These eruptions have the potential to influence global climate. The eruptive products are commonly very crystal-poor and highly degassed, yet the magma is mostly stored as crystal mushes containing small amounts of interstitial melt with elevated water content. It is unclear how magma mushes are mobilized to create large batches of eruptible crystal-free magma. Further, rhyolitic eruptions can switch repeatedly between effusive and explosive eruption styles and this transition is difficult to attribute to the rheological effects of water content or crystallinity. Here we measure the viscosity of a series of melts spanning the compositional range of the Yellowstone volcanic system and find that in a narrow compositional zone, melt viscosity increases by up to two orders of magnitude. These viscosity variations are not predicted by current viscosity models and result from melt structure reorganization, as confirmed by Raman spectroscopy. We identify a critical compositional tipping point, independently documented in the global geochemical record of rhyolites, at which rhyolitic melts fluidize or stiffen and that clearly separates effusive from explosive deposits worldwide. This correlation between melt structure, viscosity and eruptive behaviour holds despite the variable water content and other parameters, such as temperature, that are inherent in natural eruptions. Thermodynamic modelling demonstrates how the observed subtle compositional changes that result in fluidization or stiffening of the melt can be induced by crystal growth from the melt or variation in oxygen fugacity. However, the rheological effects of water and crystal content alone cannot explain the correlation between composition and eruptive style. We conclude that the composition of calcalkaline rhyolites is decisive in determining the mobilization and eruption dynamics of Earth's largest volcanic systems, resulting in a better understanding of how the melt structure controls volcanic processes.

Crystal plasticity as an indicator of the viscous-brittle transition in magmas.

Understanding the flow of multi-phase (melt, crystals and bubbles) magmas is of great importance for interpreting eruption dynamics. Here we report the first observation of crystal plasticity, identified using electron backscatter diffraction, in plagioclase in andesite dome lavas from Volcán de Colima, Mexico. The same lavas, deformed experimentally at volcanic conduit temperature and load conditions, exhibit a further, systematic plastic response in the crystalline fraction, observable as a lattice misorientation. At higher stress, and higher crystal fraction, the amount of strain accommodated by crystal plasticity is larger. Crystal plastic distortion is highest in the intact segments of broken crystals, which have exceeded their plastic limit. We infer that crystal plasticity precludes failure and can punctuate the viscous-brittle transition in crystal-bearing magmas at certain shallow magmatic conditions. Since crystal plasticity varies systematically with imposed conditions, this raises the possibility that it may be used as a strain marker in well-constrained systems.

Volcanic ash supports a diverse bacterial community in a marine mesocosm.

Shallow-water coral reef ecosystems, particularly those already impaired by anthropogenic pressures, may be highly sensitive to disturbances from natural catastrophic events, such as volcanic eruptions. Explosive volcanic eruptions expel large quantities of silicate ash particles into the atmosphere, which can disperse across millions of square kilometres and deposit into coral reef ecosystems. Following heavy ash deposition, mass mortality of reef biota is expected, but little is known about the recovery of post-burial reef ecosystems. Reef regeneration depends partly upon the capacity of the ash deposit to be colonised by waterborne bacterial communities and may be influenced to an unknown extent by the physiochemical properties of the ash substrate itself. To determine the potential for volcanic ash to support pioneer bacterial colonisation, we exposed five well-characterised volcanic and coral reef substrates to a marine aquarium under low light conditions for 3 months: volcanic ash, synthetic volcanic glass, carbonate reef sand, calcite sand and quartz sand. Multivariate statistical analysis of Automated Ribosomal Intergenic Spacer Analysis (ARISA) fingerprinting data demonstrates clear segregation of volcanic substrates from the quartz and coral reef substrates over 3 months of bacterial colonisation. Overall bacterial diversity showed shared and substrate-specific bacterial communities; however, the volcanic ash substrate supported the most diverse bacterial community. These data suggest a significant influence of substrate properties (composition, granulometry and colour) on bacterial settlement. Our findings provide first insights into physicochemical controls on pioneer bacterial colonisation of volcanic ash and highlight the potential for volcanic ash deposits to support bacterial diversity in the aftermath of reef burial, on timescales that could permit cascading effects on larval settlement.

Spatial analysis of Mount St. Helens tephra leachate compositions: implications for future sampling strategies.

Tephra particles in physically and chemically evolving volcanic plumes and clouds carry soluble sulphate and halide salts to the Earth's surface, ultimately depositing volcanogenic compounds into terrestrial or aquatic environments. Upon leaching of tephra in water, these salts dissolve rapidly. Previous studies have investigated the spatial and temporal variability of tephra leachate compositions during an eruption in order to gain insight into the mechanisms of gas-tephra interaction which emplace those salts. However, the leachate datasets analysed are typically small and may poorly represent the natural variability and complexity of tephra deposits. Here, we have conducted a retrospective analysis of published leachate analyses from the 18 May 1980 eruption of Mount St. Helens, Washington, analysing the spatial structure of the concentrations and relative abundances of soluble Ca, Cl, Na and S across the deposits. We have identified two spatial features: (1) concentrated tephra leachate compositions in blast deposits to the north of the volcano and (2) low S/Cl and Na/Cl ratios around the Washington-Idaho border. By reference to the bulk chemistry and granulometry of the deposit and to current knowledge of gas-tephra interactions, we suggest that the proximal enrichments are the product of pre-eruptive gas uptake during cryptodome emplacement. We speculate that the low S/Cl and Na/Cl ratios reflect a combination of compositional dependences on high-temperature SO2 uptake and preferential HCl uptake by hydrometeor-tephra aggregates, manifested in terrestrial deposits by tephra sedimentation and fallout patterns. However, despite our interrogation of the most exhaustive tephra leachate dataset available, it has become clear in this effort that more detailed insights into gas-tephra interaction mechanisms are prevented by the prevalent poor temporal and spatial representativeness of the collated data and the limited characterisation of the tephra deposits. Future leachate studies should aim to extensively sample across tephra deposit limits whilst simultaneously characterising deposit stratigraphy and tephra chemistry, mineralogy and granulometry, taking steps to ensure the quality and comparability of collected leachate datasets.

Seismogenic lavas and explosive eruption forecasting.

Volcanic dome-building episodes commonly exhibit acceleration in both effusive discharge rate and seismicity before explosive eruptions. This should enable the application of material failure forecasting methods to eruption forecasting. To date, such methods have been based exclusively on the seismicity of the country rock. It is clear, however, that the rheology and deformation rate of the lava ultimately dictate eruption style. The highly crystalline lavas involved in these eruptions are pseudoplastic fluids that exhibit a strong component of shear thinning as their deformation accelerates across the ductile to brittle transition. Thus, understanding the nature of the ductile-brittle transition in dome lavas may well hold the key to an accurate description of dome growth and stability. Here we present the results of rheological experiments with continuous microseismic monitoring, which reveal that dome lavas are seismogenic and that the character of the seismicity changes markedly across the ductile-brittle transition until complete brittle failure occurs at high strain rates. We conclude that magma seismicity, combined with failure forecasting methods, could potentially be applied successfully to dome-building eruptions for volcanic forecasting.

High-load, high-temperature deformation apparatus for synthetic and natural silicate melts.

A unique high-load, high-temperature uniaxial press was developed to measure the rheology of silicate melts and magmatic suspensions at temperature up to 1050 degrees C. This new apparatus is designed to operate at constant stresses (up to 300 kN) or constant strain rates (approximately 10(-7) to 10(0) s(-1)) and further allows us to carry on experiments on samples with high viscosities (approximately 10(8) to 10(12) Pa s). The rheological instrument represents an advance in that it accommodates homogeneously heated samples (+/-2 degrees C) of voluminous sizes (up to 790 cm(3)) which permit the insertion of thermocouples to monitor temperature distribution evolutions during measurements. At last this setup allows for accurate measurements of viscosity of natural multiphase materials at strain rates and temperatures common to natural systems. The apparatus aspires to precisely (1) describe the onset of non-Newtonian behavior and its evolution with increasing strain rate until the point of rupture in the brittle regime, (2) constrain the effect of crystals and bubbles on the viscosity, and (3) record heating dissipated through viscous deformation. Here, we present a series of measurements on NIST standard material SRM 717a to calibrate the instrument. We couple the viscosity determined via Gent's equation with certified viscosity data of the standard material to calibrate this state-of-the-art apparatus. This work shows that we can resolve the viscosity of voluminous melt sample within 0.06 logarithmic unit and furthermore present the detection of minor viscous dissipation for a high-temperature, high strain rate experiment.