Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, October 15th until 12:00 AM ET on Friday, October 16th due to maintenance. We apologize for the inconvenience.


This content will become publicly available on March 16, 2027

Title: Volcanic Ash Affected by Lightning During the Hunga 2022 Eruption
Abstract The 15 January 2022 eruption of Hunga volcano (Kingdom of Tonga) produced the most lightning ever documented during an explosive eruption to date. This study estimates the mass of erupted tephra that may be structurally or electromagnetically affected by the lightning, based upon lightning peak current, channel length, and ash plume particle concentration. The lightning channels totaled 1.67 million m3in volume and contained 548 kg of volcanic ash at a calculated plume concentration of 0.328 g/m3. From this total, 54.8 kg of ash may display physical evidence in the form of lightning‐induced textures, such as lightning‐induced volcanic spherules, but this is an insignificant fraction of the total airfall deposit (10−8%). However, the total mass of ash exposed to magnetic flux densities exceeding Earth's ambient surface values is 2.24 × 1014 g, corresponding to just over a third (34.9%) of the calculated mass of total airfall (5.21 × 1014 g). This study reveals that even though physical evidence of volcanic lightning may be limited, ash particles will still be affected by the electromagnetic fields generated by the lightning discharge. The extent of these effects will be a function of lightning properties, ash properties, and location of the ash in relation to the discharge channel.  more » « less
Award ID(s):
2151806
PAR ID:
10673137
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
American Geophysical Union
Date Published:
Journal Name:
Journal of Geophysical Research: Atmospheres
Volume:
131
Issue:
5
ISSN:
2169-897X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract 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}$$ 1.8-0.4+0.3 km3. The whole eruption produced > 6.3 km3of uncompacted pyroclastic material on the seafloor and a caldera volume change of 6 km3DRE. 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. 
    more » « less
  2. The 21–22 June 2019 eruption of Raikoke volcano, Russia, provided an opportunity to explore how spatial trends in volcanic lightning locations provide insights into pulsatory eruption dynamics. Using satellite-derived plume heights, we examine the development of lightning detected by Vaisala’s Global Lightning Dataset (GLD360) from eleven, closely spaced eruptive pulses. Results from one-dimensional plume modeling show that the eruptive pulses with maximum heights 9–16.5 km above sea level were capable of producing ice in the upper troposphere, which contributed variably to electrification and volcanic lightning. A key finding is that lightning locations not only followed the main dispersal direction of these ash plumes, but also tracked a lower-level cloud derived from pyroclastic density currents. We show a positive relationship between umbrella cloud expansion and the area over which lightning occurs (the ‘lightning footprint’). These observations suggest useful metrics to characterize ongoing eruptive activity in near real-time. 
    more » « less
  3. Abstract Deposition of volcanic ash is thought to impact marine biogeochemical cycling by adding soluble iron (Fe) to the surface ocean. The magnitude of this input is a function of the amount of ash deposited, the total Fe content in the ash, and ash‐derived Fe's fractional solubility. However, the relative importance of chemical composition, acidic processing by the volcanic plume, and ash particle size in determining solubility is unclear. We paired an aerosol leach meant to provide an upper limit for fractional Fe solubility with chemical analyses of ash from the Cumbre Vieja (CV) and La Soufrière eruptions, which both impacted the North Atlantic in 2021. Fe in the ash samples is <6% soluble, but Fe fractional solubility in CV ash is approximately triple that of La Soufrière ash. Compared to La Soufrière, a larger proportion of the Fe in CV ash is in silicate rather than oxide minerals, which release more soluble Fe. Elevated levels of surficial fluorine (F) also suggest that CV ash was subjected to a more fluorine‐rich eruption plume and underwent more acidic processing. Particle size does not appear to be a primary control on Fe release. We estimate that the CV eruption had a much larger impact on dissolved Fe (DFe) concentration in the surface ocean than the La Soufrière eruption because of differences in soluble Fe content and particle deposition velocity. These differences may help explain why some eruptions elicit a biological response in the ocean while others do not. 
    more » « less
  4. During a volcanic eruption, fine ash particles collide with each other and form larger aggregates. Ash aggregation plays an important role in volcanic plume modeling, as larger aggregates fall closer to the vent exit, while smaller particles can travel far before falling out. Preferential concentration can promote ash aggregation by increasing the likelihood of collisions among the particles. This process is challenging to model due to the complexities associated with turbulence and particle tracking. To improve these models, we investigate the effects of several properties on the clustering of particles in a jet flow, including the Reynolds number, Stokes number, mass loading, and relative humidity. We vertically eject compressed air embedded with particles into a quiescent chamber at Reynolds numbers between 5 × 103 and 10 × 103. We study hollow glass and solid nickel particles, each ~ 13 μm in diameter, resulting in Stokes numbers between 1.0 and 9.4. Particle mass loadings range from one to three percent, while relative humidity varies from 40 to 75%. Particles are illuminated using a laser sheet, and images are captured with a high-resolution camera. Clusters are identified by examining the densities of particles in the series of instantaneous images. At Stokes numbers near one, the number of clusters is more than an order of magnitude higher, with nearly twice as many clusters downstream. In addition, the average area and perimeter are larger than at high Stokes numbers. However, when normalized, the distribution of clusters is largely independent of Reynolds number, Stokes number, mass loading, and relative humidity. 
    more » « less
  5. Abstract. Plume-SPH provides the first particle-based simulation ofvolcanic plumes. Smoothed particle hydrodynamics (SPH) has several advantagesover currently used mesh-based methods in modeling of multiphase freeboundary flows like volcanic plumes. This tool will provide more accurateeruption source terms to users of volcanic ash transport anddispersion models (VATDs), greatly improving volcanic ash forecasts. The accuracy ofthese terms is crucial for forecasts from VATDs, and the 3-D SPH modelpresented here will provide better numerical accuracy. As an initial effortto exploit the feasibility and advantages of SPH in volcanic plume modeling,we adopt a relatively simple physics model (3-D dusty-gas dynamic modelassuming well-mixed eruption material, dynamic equilibrium and thermodynamicequilibrium between erupted material and air that entrained into the plume,and minimal effect of winds) targeted at capturing the salient features of avolcanic plume. The documented open-source code is easily obtained andextended to incorporate other models of physics of interest to the largecommunity of researchers investigating multiphase free boundary flows ofvolcanic or other origins. The Plume-SPH code (https://doi.org/10.5281/zenodo.572819) also incorporates several newly developed techniques inSPH needed to address numerical challenges in simulating multiphasecompressible turbulent flow. The code should thus be also of general interestto the much larger community of researchers using and developing SPH-basedtools. In particular, the SPH−ε turbulence model is used to capturemixing at unresolved scales. Heat exchange due to turbulence is calculated bya Reynolds analogy, and a corrected SPH is used to handle tensile instabilityand deficiency of particle distribution near the boundaries. We alsodeveloped methodology to impose velocity inlet and pressure outlet boundaryconditions, both of which are scarce in traditional implementations of SPH. The core solver of our model is parallelized with the message passinginterface (MPI) obtaining good weak and strong scalability using novel techniquesfor data management using space-filling curves (SFCs), object creationtime-based indexing and hash-table-based storage schemes. These techniques areof interest to researchers engaged in developing particles in cell-typemethods. The code is first verified by 1-D shock tube tests, then bycomparing velocity and concentration distribution along the central axis andon the transverse cross with experimental results of JPUE (jet or plume thatis ejected from a nozzle into a uniform environment). Profiles of severalintegrated variables are compared with those calculated by existing 3-D plumemodels for an eruption with the same mass eruption rate (MER) estimated forthe Mt. Pinatubo eruption of 15 June 1991. Our results are consistent withexisting 3-D plume models. Analysis of the plume evolution processdemonstrates that this model is able to reproduce the physics of plumedevelopment. 
    more » « less