Search NASA⌕ Search

SEARCH · Search NASA

Results for “explosive volcanism”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Potential of Stratospheric Aerosol Injection to Reduce the Climatic Risks of Explosive Volcanic Eruptions

Sulfur-rich volcanic eruptions happen sporadically. If Stratospheric Aerosol Injection (SAI) were to be deployed, it is likely that explosive volcanic eruptions would happen during such a deployment. Here we use an ensemble of Earth System Model simulations to show how changing the injection strategy post-eruption could be used to reduce the climate risks of a large volcanic eruption; the risks are also modified even without any change to the strategy. For a medium-size eruption (10 Tg-SO 2 ) comparable to the SAI injection rate, the volcanic-induced cooling would be reduced if it occurs under SAI, especially if artificial sulfur dioxide injections were immediately suspended. Alternatively, suspending injection only in the eruption hemisphere and continuing injection in the opposite would reduce shifts in precipitation in the tropical belt and thus mitigate eruption-induced drought. Finally, we show that for eruptions much larger than the SAI deployment, changes in SAI strategy would have minimal effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mafic explosive volcanism at Llaima Volcano: 3D x-ray microtomography reconstruction of pyroclasts to constrain shallow conduit processes

Abstract Mafic volcanic activity is dominated by effusive to mildly explosive eruptions. Plinian and ignimbrite-forming mafic eruptions, while rare, are also possible; however, the conditions that promote such explosivity are still being explored. Eruption style is determined by the ability of gas to escape as magma ascends, which tends to be easier in low-viscosity, mafic magmas. If magma permeability is sufficiently high to reduce bubble overpressure during ascent, volatiles may escape from the magma, inhibiting violent explosive activity. In contrast, if the permeability is sufficiently low to retain the gas phase within the magma during ascent, bubble overpressure may drive magma fragmentation. Rapid ascent may induce disequilibrium crystallization, increasing viscosity and affecting the bubble network with consequences for permeability, and hence, explosivity. To explore the conditions that promote strongly explosive mafic volcanism, we combine microlite textural analyses with synchrotron x-ray computed microtomography of 10 pyroclasts from the 12.6 ka mafic Curacautín Ignimbrite (Llaima Volcano, Chile). We quantify microlite crystal size distributions (CSD), microlite number densities, porosity, bubble interconnectivity, bubble number density, and geometrical properties of the porous media to investigate the role of magma degassing processes at mafic explosive eruptions. We use an analytical technique to estimate permeability and tortuosity by combing the Kozeny-Carman relationship, tortuosity factor, and pyroclast vesicle textures. The groundmass of our samples is composed of up to 44% plagioclase microlites, > 85% of which are < 10 µm in length. In addition, we identify two populations of vesicles in our samples: (1) a convoluted interconnected vesicle network produced by extensive coalescence of smaller vesicles (> 99% of pore volume), and (2) a population of very small and completely isolated vesicles (< 1% of porosity). Computed permeability ranges from 3.0 × 10 −13 to 6.3 × 10 −12 m 2 , which are lower than the similarly explosive mafic eruptions of Tarawera (1886; New Zealand) and Etna (112 BC; Italy). The combination of our CSDs, microlite number densities, and 3D vesicle textures evidence rapid ascent that induced high disequilibrium conditions, promoting rapid syn-eruptive crystallization of microlites within the shallow conduit. We interpret that microlite crystallization increased viscosity while simultaneously forcing bubbles to deform as they grew together, resulting in the permeable by highly tortuous network of vesicles. Using the bubble number densities for the isolated vesicles (0.1-3 −3 × 10 4 bubbles per mm 3 ), we obtain a minimum average decompression rate of 1.4 MPa/s. Despite the textural evidence that the Curacautín magma reached the percolation threshold, we propose that rapid ascent suppressed outgassing and increased bubble overpressures, leading to explosive fragmentation. Further, using the porosity and permeability of our samples, we estimated that a bubble overpressure > 5 MPa could have been sufficient to fragment the Curacautín magma. Other mafic explosive eruptions report similar disequilibrium conditions induced by rapid ascent rate, implying that syn-eruptive disequilibrium conditions may control the explosivity of mafic eruptions more generally.

Valdivia, Pedro↗
DOE OSTI2022

Discriminating Types of Volcanic Electrical Activity: Toward an Eruption Detection Algorithm

Abstract We present a method for classification of the two distinct types of electrical activity that occur during an explosive volcanic eruption: vent discharges and lightning. Vent discharges occur at the onset of an explosion and create a distinctive radio frequency signature called continual radio frequency. Seconds to minutes after the onset of the eruption, lightning begins to occur throughout the eruption column. We use logistic regression to classify a radio frequency impulse as being part of either a lightning flash or a period of continual radio frequency. The classifier uses the number of peaks in the amplitude envelope from 1 ms windows before and after the impulsive very high frequency waveform, with an average accuracy of 97.9%. We propose that this method could be used in an algorithm to determine when explosive eruptions occurred by identification of the distinctive signatures of vent discharges and lightning.

58 GEOSCIENCES↗
DOE OSTI2023

Volcanic Imprints in Last-Millennium Land Summer Temperatures in the Circum-North Atlantic Area

Summer cooling is one of the most direct consequences of explosive volcanic eruptions that can affect ecosystems and human societies. Recent studies revealed a multiyear cooling impact on hemispheric and global summer temperatures after tropical eruptions, yet the volcanic responses appear to vary on regional scales. Here, in this work, we revisit volcano-induced summer cooling in eastern Canada and northern and central Europe by applying superposed epoch analysis on CMIP6-PMIP4 simulations and millennial temperature reconstructions based on tree-ring density. We then examine potential causes modulating region-specific volcanic impact. While confirming that, on average, tropical eruptions over the last millennium have induced a longer cooling (>4 yr) than eruptions from extratropical Northern Hemisphere in all three North Atlantic regions, we show that the peak magnitude of cooling is stronger in eastern Canada. We also find that the detected volcanic temperature anomalies can be strongly affected by the selection and number of volcanic events and nonvolcanic signals embedded in the climate time series. This study highlights the risks of using highly noisy proxy records to investigate volcanic impacts, especially in regions with strong unforced climate variability. The CMIP6-PMIP4 simulations generally agree with the three reconstructions on the average response to tropical eruptions, but their performance is poorer regarding the production of significant cooling after extratropical eruptions. Our results further suggest that the particular sensitivity to tropical eruptions in eastern Canada is likely related to increased sea ice surrounding Quebec–Labrador associated with the positive Arctic Oscillation and North Atlantic Oscillation formed during the first posteruption winter.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Computational Prediction of Infrasound Arrival Times and Directions from Stationary and Moving Impulsive Sources

This report addresses the need to predict infrasound signal arrival times and back azimuths at monitoring stations, enabling more focused and efficient searches within recorded waveform data. The primary challenge is estimating expected signal arrival windows for stationary and moving acoustic sources, such as chemical explosions, volcanic eruptions, meteoroids, and spacecraft re-entry events. To address this challenge, a reproducible methodology is described that uses simplified propagation speeds for boundary layer, tropospheric, stratospheric, and thermospheric atmospheric waveguides. While the Python source code itself is not freely available, this document provides detailed, step-by-step instructions, and equations enabling users to replicate and adapt the method independently. The method reliably predicts signal arrival intervals and back azimuths, thereby supporting rapid detection and accurate interpretation of infrasound events. Results demonstrate that this method effectively identifies plausible signal arrival intervals and directions, facilitating faster event detection and more reliable interpretation. This methodology directly supports atmospheric monitoring, planetary defense, and forensic analysis of explosive atmospheric events.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Infrasound From a Buried Chemical Explosion Recorded on a Balloon in the Lower Stratosphere

Abstract Infrasound waveforms generated by natural and anthropogenic phenomena contain important clues about the size and nature of the event. We show that sensors on balloons in the lower stratosphere can record faithful representations of the near‐source acoustic wave field at unprecedented range. The acoustic signature of a buried chemical explosion recorded at a range of 56 km and an altitude of 21.8 km was nearly identical to that recorded on the ground 0.5 km from the epicenter, but absent on a ground sensor located 46 km away. Our results demonstrate that balloon‐borne infrasound techniques greatly increase the range at which well‐preserved acoustic representations of near‐source physics can be acquired, and that their propagation is simple to model. Our work has implications for monitoring remote regions of the earth for explosions, volcanic eruptions, and other phenomena. It also supports the prospect of balloon‐based infrasound seismology on Venus.

Bowman, Daniel C.↗
DOE OSTI2021

Synthetic Evaluation of Infrasonic Multipole Waveform Inversion

Acoustic source inversions estimate the mass flow rate of volcanic explosions or yield of chemical explosions and provide insight into potential source directionality. However, the limitations of applying these methods to complex sources and their ability to resolve a stable solution have not been investigated in detail. We perform synthetic infrasound waveform inversions that use 3-D Green’s functions for a variety of idealized and realistic deployment scenarios using both a flat plane and Yasur volcano, Vanuatu as examples. We investigate the ability of various scenarios to retrieve the input source functions and relative amplitudes for monopole and multipole (monopole and dipole) inversions. Infrasound waveform inversions appear to be a robust method to quantify mass flow rates from simple sources (monopole) using deployments of infrasound sensors placed around a source, but care should be taken when analyzing and interpreting results from more complex acoustic sources (multipole) that have significant directional components. In the examples we consider the solution is stable for monopole inversions with a signal-to-noise ratio greater than five and the dipole component is small. For most scenarios investigated, the vertical dipole component of the multipole explosion source is poorly constrained and can impact the ability to recover the other source term components. Because multipole inversions are ill-posed for many deployments, a low residual does not necessarily mean the proper source vector has been recovered. Synthetic studies are able to help investigate the limitations and place bounds on information that may be missing using monopole and multipole inversions for potentially directional sources.

58 GEOSCIENCES↗
DOE OSTI2022

Bayesian Seismoacoustic Source Location: Acoustic Approach

Seismic waves and infrasound are key technologies in the International Monitoring System (IMS) to monitor explosive events in the solid Earth and atmosphere. Energetic man-made or natural events (e.g., chemical/nuclear explosions, volcanic eruptions, and earthquakes) near the Earth’s surface produce both ground motion and atmospheric pressure disturbances which propagate as seismic waves and infrasound, respectively. Seismic waves have been generally used to detect and identify underground and near-surface events (Myers, et. al., 2007), and infrasound are sensitive to events near the surface or in the atmosphere (Modrak et. al., 2010). Due to their different sensitivities to events, they can complement to each other to improve the event detection and discrimination. The framework of joint seismoacoustic event location has recently reviewed by theoretical research (Koch and Arrowsmith, 2019). Although the early applications showed promising results to improve the accuracy of event location, their application were still limited to a small set of events selected to prove the concepts, and practical capability of the method for operational purpose is not fully evaluated with data. Our final goal is to apply the method of seismoacoustic event location to a larger set of events and evaluate its applicability for operational event location in practice. To that end, we focus on developing and verifying acoustic source location method in this study.

58 GEOSCIENCES↗
DOE OSTI2026

Stratospheric aerosol forcing for CMIP7 – Part 1: optical properties for pre-industrial, historical, and scenario simulations

Stratospheric aerosols, most of which originate from explosive volcanic sulfur emissions into the stratosphere, are a key natural driver of climate variability. They are thus a forcing provided by the Coupled Model Intercomparison Project (CMIP) Climate Forcings Task Team to climate modelling groups participating in phase 7 of CMIP. For the historical period, we provide two datasets covering 1750–2023: (i) a volcanic upper tropospheric-stratospheric sulfur emission dataset, documented in a companion paper; and (ii) a stratospheric sulfate aerosol optical property dataset, documented in this paper. For the satellite era (from 1979 onwards), stratospheric aerosol optical properties are derived from the Global Space-based Stratospheric Aerosol Climatology (GloSSAC) dataset. For the pre-satellite era (1750–1978), optical properties are derived from our volcanic SO 2 emission dataset using a new version of the reduced-complexity volcanic aerosol model Easy Volcanic Aerosol (Height) (EVA_H v2). A background, non-volcanic stratospheric aerosol climatology is derived from the 1998–2001 period with a trend over 1850–1978 accounting for increasing anthropogenic aerosols. A monthly stratospheric aerosol climatology is derived from the 1850–2021 average for both pre-industrial and scenario (future) simulations, with a 9-year ramp over 2022–2030 for scenario simulations to ensure a smooth transition from the historical period. CMIP7 simulations should use version 2.2.1 of our pre-industrial and historical datasets, and version 2.2.2 of our scenario dataset. Our methodology to produce historical aerosol optical properties significantly differs from CMIP6 for the pre-satellite era, and the resulting forcings in turn largely differ. In particular, the CMIP6 dataset was mostly based on the sparse and uncertain pyrheliometer record, which resulted in strongly underrepresented emissions from small-to-moderate magnitude eruptions. The resulting bias is addressed in CMIP7, which is entirely emission-derived in the pre-satellite era and uses more recent ice-core-based volcanic sulfur emission inventories than CMIP6. Our approach results in an overall larger volcanic aerosol forcing for CMIP7, with the 1850–2014 mean mid-visible global mean stratospheric aerosol optical depth (SAOD) in CMIP7 (0.0138) being 29 % higher than in CMIP6 (0.0107). The pre-industrial mean of the same variable is 26 % higher in CMIP7 (0.0135, derived from the historical 1850–2021) than CMIP6 (0.0107, derived from the historical 1850–2014 mean). Using a reduced-complexity climate model, we simulate a global mean surface temperature that is 0.07 °C colder for 1850–1900 when using the CMIP7 dataset instead of CMIP6, whereas 2000–2014 is 0.03 °C warmer in CMIP7. Our dataset also exhibits lower forcing for 1960–1980, resulting in temperatures 0.06 °C warmer when averaged over 1960–1990, a period for which CMIP6 climate models exhibit a cold bias. Given the large uncertainties characterizing the dataset, in particular for the pre-satellite era, we advise against treating the CMIP7 or CMIP6 dataset as uniquely superior for any specific year and highlight the need for further evaluation. We conclude the study by discussing sources of uncertainty for the dataset, future research avenues to improve it, as well as requirements to operationalize the production of the dataset, i.e. extend it and update it on an annual basis instead of every 5–7 years following CMIP cycles.

Aubry, Thomas J. [Univ. of Exeter, Devon (United K↗
DOE OSTI2023

Lightning Rings and Gravity Waves: Insights Into the Giant Eruption Plume From Tonga's Hunga Volcano on 15 January 2022

Abstract On 15 January 2022, Hunga Volcano in Tonga produced the most violent eruption in the modern satellite era, sending a water‐rich plume at least 58 km high. Using a combination of satellite‐ and ground‐based sensors, we investigate the astonishing rate of volcanic lightning (>2,600 flashes min −1 ) and what it reveals about the dynamics of the submarine eruption. In map view, lightning locations form radially expanding rings. We show that the initial lightning ring is co‐located with an internal gravity wave traveling >80 m s −1 in the stratospheric umbrella cloud. Buoyant oscillations of the plume's overshooting top generated the gravity waves, which enhanced turbulent particle interactions and triggered high‐current electrical discharges at unusually high altitudes. Our analysis attributes the intense lightning activity to an exceptional mass eruption rate (>5 × 10 9 kg s −1 ), rapidly expanding umbrella cloud, and entrainment of abundant seawater vaporized from magma‐water interaction at the submarine vent.

58 GEOSCIENCES↗
DOE OSTI2025

Infrasonic directivity of monopole, dipole and bipole ground-surface reflected sources

Infrasound (acoustic waves below 20 Hz) can be used to detect, locate and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local (⁠<5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modelling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analysing geophysical field-scale experimental data; however, especially at very local distances from the source (⁠<5 km), 3-D simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications.

Infrasound↗
DOE OSTI2025

The Hunga Volcanic Eruption Atmospheric Impacts Report

On 15 January 2022 a highly explosive eruption of the Hunga volcano occurred in the Kingdom of Tonga in the South Pacific Ocean (175°24’ W, 20°33’ S). The Volcanic Explosivity Index (VEI) 6 eruption originated from a shallow submarine vent, making it distinct from large subaerial eruptions of recent decades (e.g., 1982 El Chichón, 1991 Mt. Pinatubo). In particular, seawater enhanced explosivity and dampened sulfur dioxide (SO 2 ) emissions. The eruption was the culmination of ~1 month of precursory activity; however, the timing and size of the eruption were unexpected, partly due to the challenges of monitoring submarine volcanoes. The stratospheric hydration caused by the eruption was unprecedented in magnitude, altitude, and duration in the satellite record. This Executive Summary reflects the current assessment of the Hunga eruption and its impact on the climate system. We report key observations of the eruption and its aftermath, as well as simulations of its impact by global chemistry-climate models. The Hunga eruption had an unprecedented impact on the stratosphere and mesosphere due to the plume height and large water content, which increased the global stratospheric water vapour burden by 10%. Most of this water has remained in the atmosphere into 2025. However, Hunga’s net impact on surface climate was small compared to that of earlier large-magnitude volcanic eruptions, due to limited sulfate aerosol loading in the stratosphere and the high altitude of the water vapour injection.

58 GEOSCIENCES↗
DOE OSTI2022

Volcanic electrification: recent advances and future perspectives

Abstract The electrification of volcanic plumes has been described intermittently since at least the time of Pliny the Younger and the 79 AD eruption of Vesuvius. Although sometimes disregarded in the past as secondary effects, recent work suggests that the electrical properties of volcanic plumes reveal intrinsic and otherwise inaccessible parameters of explosive eruptions. An increasing number of volcanic lightning studies across the last decade have shown that electrification is ubiquitous in volcanic plumes. Technological advances in engineering and numerical modelling, paired with close observation of recent eruptions and dedicated laboratory studies (shock-tube and current impulse experiments), show that charge generation and electrical activity are related to the physical, chemical, and dynamic processes underpinning the eruption itself. Refining our understanding of volcanic plume electrification will continue advancing the fundamental understanding of eruptive processes to improve volcano monitoring. Realizing this goal, however, requires an interdisciplinary approach at the intersection of volcanology, atmospheric science, atmospheric electricity, and engineering. Our paper summarizes the rapid and steady progress achieved in recent volcanic lightning research and provides a vision for future developments in this growing field.

58 GEOSCIENCES↗
DOE OSTI2023

Distributed Acoustic Sensing as a Monitoring Tool at LANL [Slides]

DAS records a variety of signals, including (1) earthquakes and volcanic eruptions, (2) explosions (DAG, mining, etc.), (3) acoustic waves, (4) ocean (waves, whales), (5) vehicle, ship and pedestrian traffic, (6) glaciology and landslides, and (7) oil and gas. Main questions include: (1) what is the extent of the signal content? (2) how does DAS compare to traditional sensors? (3) how do interrogators, fiber types, and fiber depth burials impact measurements? (4) special applications offshore to overcome the lack of instrumentation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Multi-Event Explosive Seismic Source for the 2022 M w 6.3 Hunga Tonga Submarine Volcanic Eruption

The eruption of the Hunga Tonga–Hunga Ha’apai submarine volcano on 15 January 2022 produced a variety of geophysical responses, including a significant seismic signal. We study the seismic source process of this event by inverting for moment tensors (MTs) using regional surface waves (Rayleigh, Love). By comparing inversion results for the eruption with eight nearby earthquakes, we show that it is possible to discriminate MT source types. Our inversion yields a shallow explosive source for the eruption and reveals the importance of trade-offs among depth, magnitude, and source type. We illustrate these trade-offs by representing the misfit variations over the eigenvalue lune. Finally, we invert for the source-time function of the sequence of explosions that occurred in the first minutes of the eruption. The multi-event source-time function comprises four subevents spanning ~270 s, with a total magnitude estimate of M w 6.34 ± 0.10.

58 GEOSCIENCES↗
DOE OSTI2024

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Permeability scaling relationships of volcanic tuff from core to field scale measurements

A recent chemical explosive test in P-Tunnel at the Nevada National Security Site, Nevada, USA, was conducted to better understand how signals propagate from explosions in the subsurface. A primary signal of interest is the migration of gases that can be used to differentiate chemical from nuclear explosions. Gas migration is highly dependent on the rock permeability which is notoriously difficult to determine experimentally in the field due to a potentially large dependence on the scale over which measurements are made. Here, we present pre-explosion permeability estimates to characterize the geologic units surrounding the recent test. Permeability measurements were made at three scales of increasing size: core samples (≈2 cm), borehole packer system tests (≈1 m), and a pre-shot cavity pressurization test (> 10 m) across ten tuff units. Permeability estimates based on core measurements showed little difference from borehole packer tests. However, permeability in most rock units calibrated from cavity pressurization tests resulted in higher permeability estimates by up to two orders of magnitude. Here, we demonstrate that the scale of the measurement significantly impacts the characterization efforts of hydraulic properties in volcanic tuff, and that local-scale measurements (< 10 m scale) do not incorporate enough heterogeneity to accurately predict field-scale flow and mass transport.

Environmental sciences↗
DOE OSTI2023

The 2022 Hunga eruption exercises the nuclear explosion monitoring community

The 2022 Hunga eruption was a global event recorded by all three Comprehensive Test Ban Treaty Organization (CTBTO) International Monitoring System (IMS) geophysical waveform sensor technologies: seismic, infrasound, and hydroacoustic. Large volcanic eruptions are rare and additionally present an opportunity to test seismoacoustic detection, location, and identification methods for nuclear explosion monitoring.

58 GEOSCIENCES↗