Eruption Column Research Articles

10 keV electron irradiation of methane ices at ocean world surface temperatures

The surface environmental conditions of many of the icy ocean worlds in the Solar System are outside of the thermal stability field of methane ice. However, mechanisms may exist that entrain or trap methane within water ice that arrives at the surface via plume eruptions or resurfacing events. At the surface, the ice grains may experience charged particle irradiation that could lead to the production of new complex hydrocarbons. We investigated the irradiation products of 12CH4 versus 13CH4, the temperature dependence of CH4 radiation processing and product stability, and the temperature stability of CH4 when mixed with H2O in a cryogenic vacuum environment. Significantly, methane ice experimentally encased in water ice was irradiated at 100 K in an ultra-high vacuum environment and produced similar radiation-formed secondary phases as pure methane ice irradiated at 20 K. The methane radiation products produced outside of the methane ice thermal stability field, while encased in a water ice film, were most notably 13CO2 and a candidate hydrocarbon feature. Evidence for complex hydrocarbons, such as C2H6, was observed during warming and sublimation of the sample.

Risk predictive score and cord morphology classification for intraoperative neuromonitoring alerts in kyphosis surgery

BACKGROUNDIntraoperative neuromonitoring (IONM) alert is one of the worrying events of kyphosis corrective surgery, which can result in a postoperative neurological deficit. To our knowledge, there is no risk prediction score to predict such events in patients undergoing kyphosis surgery. PURPOSETo develop a new preoperative MRI-based cord morphology classification (CMC) and risk prediction score for predicting IONM alerts in patients with kyphotic deformity. STUDY DESIGNRetrospective analysis of prospectively collected data. PATIENT SAMPLEAbout 114 patients undergoing surgical correction for kyphotic deformity. OUTCOME MEASURESIntraoperative neuromonitoring alerts and postoperative neurological status using AIS grading. METHODSKyphotic deformity patients undergoing posterior spinal fusion were retrospectively reviewed. Based on the morphology of the spinal cord and surrounding CSF in MRI, there are 5 types of cord. Type 1 (normal cord): circular cord with surrounding visible CSF between the cord and the apex, Type 2 (flattened cord): cord with <50% distortion at the apex with obliteration of the anterior CSF; Type 3 (deformed cord): cord with >50% distortion at the apex with complete obliteration of the surrounding CSF; Type 4 (stretched cord): the cord is stretched and atrophied over the apex of the curve. Type 5 (translated cord): horizontal translation of the cord at the apex with buckling collapse of the vertebral column. Preoperative radiographs were used to measure the preoperative sagittal cobbs angle, sagittal deformity angular ratio (S-DAR), sagittal vertical axis (SVA), apex of the curve, and type of kyphosis. Clinical data like the duration of symptoms, clinical signs of myelopathy, neurological status (AIS grade), grade of myelopathy using the mJOA score, and type of osteotomy were documented. Multivariate logistic regression was used to determine the risk factors for IONM alerts and the risk prediction score was developed which was validated with new cohort of 30 patients. RESULTSA total of 114 patients met the inclusion criteria. IONM alerts were documented in 33 patients (28.9%), with full recovery of the signal in 25 patients and a postoperative deficit in 8 patients. Rate of IONM alerts was significantly higher in Type 5 (66%), followed by Type 4 (50%), Type 3 (21.1%), Type 2 (11.1%), and Type 1 (11.1%) (p-value<.001). Based on multiple logistic regression, 7 factors, namely preoperative neurological status, mJOA score≤6, presence of signs of myelopathy, apex of the curve above T5, preoperative sagittal cobbs, S-DAR, and MRI-based CMC, were identified as risk predictors. The value for the risk factors varies from 0 to 4, and the maximum total risk score was 13. The cut-off value of 6 had good sensitivity (84.9%) and specificity (77.8%) indicating a high risk for IONM alerts. The AUC of the predictive model was 0.92, indicating excellent discriminative ability. CONCLUSIONWe developed and validated a risk predictive score that identifies patients at risk of IONM alerts during kyphosis surgery. Identification of such high-risk patients (risk score≥6) helps in proper evaluation and preoperative counselling and helps in providing a proper evidence-based reference for treatment strategies.

Evolution of a large-scale phreatoplinian eruption: Constraints from the 40 ka caldera-forming eruption of Kutcharo volcano, eastern Hokkaido, Japan

“Phreatoplinian” is an explosive phreatomagmatic eruption style that is defined by the fragmentation of magma and widespread dispersal of the resulting fine ash and accretionary lapilli. These eruptions pose significant future risks at caldera volcanoes that host lakes and abundant groundwater. There have been no direct observations of a phreatoplinian eruption, therefore, constraining the detailed mechanisms and sequences of such events relies on studying the deposits of previous eruptions. In order to advance our understanding of these hazardous phenomena we conducted a case study of the 40 ka caldera-forming eruption (Kp I) from Kutcharo volcano in eastern Hokkaido, Japan. We subdivided Kp I eruption deposits into 7 units (Units 1 to 7 in ascending order). Units 1 to 6 are air fall deposits consisting of alternating thin pumice and thick silty ash layers with abundant spherical accretionary lapilli. Stratigraphically higher ash fall units are thicker, finer in grain-size, and more widely distributed. The maximum eruption column height and mass-discharge rate were calculated to be 40 km and 1.4 × 109 kg/s, respectively. Unit 7 is a climactic ignimbrite (76 km3), which is distributed widely over the area north of Kutcharo caldera.Unit 6 is the largest air fall unit and can be considered to have been deposited by a phreatoplinian eruption, given its abundant accretionary lapilli, wide dispersion, and high degree of fragmentation. Unit 6 had the highest mass discharge rate (1.4 × 109 kg/s), suggesting the interaction between magma and external water was most intense, and it is thought that a large eruption column covered eastern Hokkaido. In addition, Kp I eruption deposits commonly contain glass shards derived from fragmentation via both magma degassing and Molten Fuel Coolant Interaction (MFCI). To account for this observation, we infer that the conduit penetrated a large aquifer, and the margin of the ascending magma came into contact with this external water source. Due to repeated caldera-forming eruptions, intra-caldera filled deposits (hosting a large aquifer) likely played a key role in supplying external caldera lake water to a level near the fragmentation depth of H2O-saturated felsic magma. The occurrence of these intra-caldera conduit and caldera-lake systems may provide the required conditions for phreatoplinian eruptions at continental arc caldera volcanoes in Japan and globally.

Transition of dome formation to sudden explosive eruptions at Popocatépetl, Mexico: magnetic indicators

Transitions from effusive to explosive activity can increase hazards making it crucial to define early indicators such as changes in the magnetic signals. After more than 80 cycles of crater-dome extrusion and destruction from 1996 on, Popocatépetl volcano (Mexico) experienced changes in its behavior from March 15 to 18 July 2019, when no lava domes were observed. Some of the domes behaved as contained lava flows within the crater floor (pancakes) while others were more irregular-shaped. Activity decreased considerably over this 2019 interval except for the unexpected explosions in March and June, that produced ash plumes reaching up to 14,000 m a.s.l. In order to investigate the causes of the transition from effusive to explosive behavior in March and June, we analyzed the time series from the magnetic monitoring network at Popocatépetl volcano between October 2018 and December 2019. The raw signals were analyzed by weighted differences (WD) based on the elimination of non-local changes from the total intensity values of the geomagnetic field and the discrete-time continuous wavelet transform was used to evaluate the local variations of energy within the time series. The high energy periods (linked to negative magnetic anomalies) are induced by magma ascent associated with movement within the conduit. They indicate that the sudden explosions were due to the ascent of several magma batches that were slowed during ascent and were not able to reach the surface. Changes in the rheology of the lava are linked to the influx of several batches of magma with different compositions as well as to compaction by gas loss when ascending andesitic magma pushed out overlying more viscous degassed magma clearing the conduit, which can explain why these sudden explosions were more energetic. Several geophysical data sets as well as tephra compositions were integrated to support this conclusion. The correlated multiparameters also confirm that geomagnetic volcano monitoring has been essential in understanding the processes that drive the observed changes in eruptive behavior. We present new evidence for the detection of transient events produced by magma ascent and changes in the feeding system of Popocatépetl volcano with wavelet analysis. Detailed vulcanomagnetic processing, especially when it is correlated with other monitoring parameters, provides information on ascending magma and several conduit processes that would otherwise be camouflaged. Ascending batches may precede an eruption but they can also ascend in several pulses indicating how dome growth occurs.

Did steam boost the height and growth rate of the giant Hunga eruption plume?

The eruption of Hunga volcano on 15 January 2022 produced a higher plume and faster-growing umbrella cloud than has ever been previously recorded. The plume height exceeded 58 km, and the umbrella grew to 450 km in diameter within 50 min. Assuming an umbrella thickness of 10 km, this growth rate implied an average volume injection rate into the umbrella of 330–500 km3 s−1. Conventional relationships between plume height, umbrella-growth rate, and mass eruption rate suggest that this period of activity should have injected a few to several cubic kilometers of rock particles (tephra) into the plume. Yet tephra fall deposits on neighboring islands are only a few centimeters thick and can be reproduced using ash transport simulations with only 0.1–0.2 km3 erupted volume (dense-rock equivalent). How could such a powerful eruption contain so little tephra? Here, we propose that seawater mixing at the vent boosted the plume height and umbrella growth rate. Using the one-dimensional (1-D) steady plume model Plumeria, we find that a plume fed by ~90% water vapor at a temperature of 100 °C (referred to here as steam) could have exceeded 50 km height while keeping the injection rate of solids low enough to be consistent with Hunga’s modest tephra-fall deposit volume. Steam is envisaged to rise from intense phreatomagmatic jets or pyroclastic density currents entering the ocean. Overall, the height and expansion rate of Hunga’s giant plume is consistent with the total mass of fall deposits plus underwater density current deposits, even though most of the erupted mass decoupled from the high plume. This example represents a class of high (> 10 km), ash-poor, steam-driven plumes, that also includes Kīlauea (2020) and Fukutoku-oka-no-ba (2021). Their height is driven by heat flux following well-established relations; however, most of the heat is contained in steam rather than particles. As a result, the heights of these water-rich plumes do not follow well-known relations with the mass eruption rate of tephra.

Rapid Primary Sulfate Aerosol Generation Observed With OP‐FTIR in the Eruptive Plume of the Fagradalsfjall Basaltic Eruption, Iceland, 2021

AbstractOpen‐Path Fourier‐Transform Infrared (OP‐FTIR) absorption spectroscopy is a powerful method for remote characterization of volcanic plume composition from safe distances. Many studies have used it to examine the composition of volcanic gas emitted at the surface, which is influenced by initial volatile contents and magma ascent/storage processes, and help to reveal the dynamics controlling surface activity. However, to evaluate the health hazard threats associated with volcanic emissions and their potential impact on wider atmospheric conditions, near‐source particle measurements are also key. Here we present a forward model and fitting algorithm which allows quantification of particle size and abundance. This was successfully applied to radiometrically uncalibrated OP‐FTIR spectra collected with a highly dynamic radiation source during the Fagradalsfjall eruption, Iceland, on 11 August 2021. Quantification of plume temperatures ranging from 350 to 650 K was essential to characterize the emission‐absorption behavior of SO2, enabling retrievals of particulate matter in the thermal infrared spectral window (750–1250 cm−1) in each spectrum. For the first time, we observe the rapid formation of primary aerosols in young plumes (only a few seconds old) with OP‐FTIR. Temperature‐dependent SO2/SO42− molar ratios range from 100 to 250, consistent with a primary formation mechanism controlled by cooling and entrainment of atmospheric gases. This novel aerosol spectrum retrieval opens new frontiers in field‐based measurements of sulfur partitioning and volcanic plume evolution, with the potential to improve volcano monitoring and quantification of air quality hazard assessments.

Evolving Particles in the 2022 Hunga Tonga—Hunga Ha'apai Volcano Eruption Plume

AbstractThe Multi‐angle Imaging SpectroRadiometer (MISR) aboard NASA's Terra satellite observed the Hunga Tonga—Hunga Ha'apai (HTHH) 15 January eruption plume on eight occasions between 15 and 23 January 2022. From the MISR multi‐angle, multi‐spectral imagery we retrieve aerosol plume height geometrically, along with plume‐level motion vectors, and derive radiometrically constraints on particle effective size, shape, and light‐absorption properties. Parts of two downwind aerosol layers were observed in different places and times, one concentrated in the upper troposphere (11–18 km ASL), and a mid‐stratosphere layer ∼23–30+ km ASL. After the initial day (1/15), the retrievals identified only spherical, non‐light‐absorbing particles, typical of volcanic sulfate/water particles. The near‐tropopause plume particles show constant, medium‐small (several tenths of a micron) effective size over 4 days. The mid‐stratosphere particles were consistently smaller, but retrieved effective particle size increased between 1/17 and 1/23, though they might have decreased slightly on 1/22. As a vast amount of water was also injected into the stratosphere by this eruption, models predicted relatively rapid sulfate particle growth from the modest amounts of SO2 gas injected by the eruption to high altitudes along with the water (Zhu et al., 2022, https://doi.org/10.5194/acp‐22‐10267‐2022). MISR observations up to 10 days after the eruption are consistent with these model predictions. The possible decrease in stratospheric particle size after initial growth was likely caused by evaporation, as the plume mixed with drier, ambient air. Particles in the lower‐elevation plume observed on 1/15 were larger than all the downwind aerosols and contained significant non‐spherical (likely ash) particles.

Density Stratification and Buoyancy Evolution in Pyroclastic Density Currents

AbstractPyroclastic density currents (PDCs) are density‐stratified along their vertical axis, with the near‐bed portion being denser than the upper portion, resulting from particle settling and ambient air entrainment at current margins. Whereas vertical density stratification likely influences mixing, sedimentation, and buoyancy of PDCs, many depth‐averaged models of PDC dynamics assume currents are well‐mixed. We investigated this discrepancy by performing sub‐aqueous laboratory experiments and conducted complementary numerical simulations to interrogate current dynamics at finer scales. Currents with small temperature difference with the ambient fluid become density‐stratified during propagation. The dynamics of such currents resemble two‐phase flows, in which particles move freely and particle concentration becomes stratified, but fluid density remains constant. Currents with large temperature difference with the ambient fluid, however, do not develop density stratification during propagation, due to current dynamics becoming dominated by the fluid phase and the lessening importance of particles. Currents that develop density stratification do not lift off from the bed within the domain of the setup, whereas poorly stratified currents do lift off, forming a rising plume. Strong density stratification within currents inhibits turbulence production, preventing entrained ambient fluid on current edges from mixing into current interiors. Poorly stratified currents are highly turbulent, have vigorous internal mixing, thereby achieving lift‐off. The strongly stratified currents are analogous to PDCs that result from eruption column collapse, maintaining fast velocity, low internal mixing, and high temperature over long distances. The poorly stratified currents are analogous to dilute ash‐cloud surges that develop atop basal avalanches, having short runout distances.

The 23–24 March 2021 lava fountain at Mt Etna, Italy

In 2021, more than 50 paroxysmal episodes occurred at the South-East Crater (SEC) of Mt Etna, Italy. The 23–24 March lava fountain was one of the longest episodes and began with weak Strombolian explosions, gradually transitioning to lava fountaining. The eruption intensity then dropped more slowly than in previous episodes, resulting in pulsating Strombolian explosions dominated by ash emission. Thirty-four tephra samples were used to reconstruct the fallout dispersal and estimate the total erupted mass. Grain size, textural, petrological and geochemical analyses indicate different features and were compared with the gas phase (SO2 and HCl) in the volcanic plume. By applying stochastic global optimization to simulations of the temporal evolution of the eruption column height and tephra dispersal and deposition, the total erupted mass retrieved (6.76 × 108 kg) matches well the total erupted mass estimation by the ground-based deposit (8.03 ± 2.38 × 108 kg), reducing the column height throughout the episode from 6.44 to 4.5 km above sea level and resulting in a mass eruption rate ranging from 1.96 × 105 to 8.18 × 103 kg/s. The unusual duration of the March episode and the characteristics of the erupted products point to the change in explosive style and magma fragmentation from fountaining to ash emission phases, associated with a slower magma supply inducing a change in magma rheology and a final, prolonged ash generation. Furthermore, this study showed that using observational data and the variation in eruption source parameters for numerical simulations can improve the accuracy of predicting the dispersal plume, thus mitigating the potential impact of longer paroxysmal episodes.

Near-real-time multiparametric seismic and visual monitoring of explosive activity at Sabancaya volcano, Peru

This study presents the development of a multiparametric system that utilizes artificial intelligence techniques to identify and analyze volcanic explosions in near real-time. The study analyzed 1343 explosions recorded between 2019 and 2021, along with seismic, meteorological, and visible image data from the Sabancaya volcano. Deep learning algorithms like the U-Net convolutional neural network were used to segment and measure volcanic plumes in images, while boosting-based machine learning ensembles were used to classify seismic events related to ash plumes. The findings demonstrate that these approaches effectively handle large amounts of data generated during seismic and eruptive crises. The U-Net network achieved precise segmentation of volcanic plumes with over 98% accuracy and the ability to generalize to new data. The CatBoost classifier achieved an average accuracy of 94.5% in classifying seismic events. These approaches enable the real-time estimation of eruptive parameters without human intervention, contributing to the development of early warning systems for volcanic hazards. In conclusion, this study highlights the feasibility of using seismic signals and images to detect and characterize volcanic explosions in near real-time, making a significant contribution to the field of volcanic monitoring.

Large-Scale Laboratory Experiments On The Wave Generation Due To The Collapse Of Partially And Fully Submerged Granular Columns

Landslides that occur in coastal environments can drive cascading consequences such as wave forces, flooding, and infrastructure damage to coastal communities. It can be difficult to classify these slides as subaerial or submarine, and the mechanics of wave generation associated with partially submerged failures are not well understood. Limited physical modelling has been conducted that encompasses both the triggering of granular landslides and subsequent waves associated with partially and fully submerged mass movements. To date, laboratory work investigating tsunamis generated by submarine landslides has focused on the wave formed in the direction of the mass movement (seaward direction) for rigid block experiments (eg. Rzadkiewicz, 1997) and deformable slide masses (e.g. Grilli et al., 2017, Takabatake, 2020, Bullard et al., 2023). From these experimental data sets, predictive relationships connecting slide acceleration, mass, and initial submergence depth to the amplitude of the wave formed have been presented for the seaward wave. Such relationships have not been presented for the landward directed wave, which propagates in the opposite direction of the submarine landslide motion. Further, not all landslides are easily classified as either subaerial or submarine. Consider the 2018 Anak Krakatoa landslide in which the sliding surface was estimated to be 100 m below sea level (Pakoksung et al., 2020), resulting in one third of the total collapse being submerged. In comparison to the end-member conditions of subaerial and submarine failures, the mechanics of wave generation associated with partially submerged failures is much less clear. Granular column collapse experiments provide an idealized experimental framework to explore momentum transfer processes and the resulting waves generated in partially submerged and fully submerged conditions. Work by Cabrera et al., (2020) made use of granular collapse experiments of partially to fully submerged columns to derive a continuous momentum-based function to estimate the maximum seaward wave amplitude based on the initial column submergence ratio (Hw/Ho). However, these experiments were conducted at a small-scale (Ho = 0.15 m) with a width of one particle (2.4 mm diameter). To address this research gap, a series of 22 large-scale granular collapse experiments were conducted by releasing columns of river stone (0.75 m and 0.50 m high) into a laboratory flume reservoir with water depths ranging up to 1.10 m to explore the wave generation and runup processes in both seaward and landward directions. The columns were released by a rapid pneumatically actuated vertical rising gate designed to enable the near instantaneous loss of support of the source volumes resulting in granular collapse. The failure mechanics were captured with high-speed cameras (Figure 1a,b) and wave amplitudes were measured using wave capacitance gauges (Figure 1c,d). This work also provides the first experimental data set of the landward propagating wave and runup associated with submerged granular collapse experiments. Overall, the seaward wave amplitudes measured in these highly-instrumented, large-scale physical models agree with empirical relationships developed in a previous study using smaller-scale models.

Video-based measurements of the entrainment, speed and mass flux in a wind-blown eruption column

On May 4 2010 a wind-blown ash plume issued from Eyjafjallajökull volcano in Iceland. Analysis of a 17-minute-long video recording of the eruption suggests that, within 2–2.5 km of the vent, the flow was moving with the wind and rising under buoyancy, following a trajectory directly analogous with laboratory experiments of turbulent buoyant plumes in a cross-flow. The radius of the time-averaged ash cloud grew with height z at a rate r = 0.48 z, corresponding to an entrainment coefficient 0.4, again consistent with laboratory experiments. By analysing the frames in the video and comparing the shape of the plume to that predicted by the model, we estimate that during the 17 minutes recorded, the eruption rate gradually decreased by about 43% from an initial rate of 1.11 × 104 kg s−1 to 0.63 × 104 kg s−1. The analysis reported herein opens the way to assess eruption rates and eruption column processes from video recordings during explosive volcanic eruptions.

Coupled peridynamic modeling of hydraulic fracturing in rocks

This paper presents a fully coupled peridynamics method for simulating hydraulic fracturing in rocks. The proposed method involves a rigorous coupling between the classical total-Lagrangian formulation and a new semi-Lagrangian formulation within the framework of peridynamics. The total-Lagrangian formulation is employed to simulate a rock subjected to fluid-driven fracturing, whereas the semi-Lagrangian formulation is used to solve the Navier–Stokes equations for fluids in a nonlocal manner. A coupling algorithm is developed between the two formulations to model fluid–solid interactions without causing unphysical penetration. The proposed method was validated through a simulation of a water column collapse problem. It was further validated using a Kristianovich–Geertsma–De Klerk (KGD) fracturing problem, where a plane strain fracture propagates owing to fluid injection. The proposed method reasonably captures the fracturing pattern, fracture propagation speed, and injection pressure better than the analytical solution. It offers a unified framework within the peridynamics theory for modeling solid, fluid, and fluid-driven fractures.

3D SPH implementation to large deformation of granular mass using an advanced hypoplastic constitutive model

In this study, an advanced hypoplastic model incorporating the critical state is integrated into smoothed particle hydrodynamics (SPH) for the numerical investigation of geotechnical large-deformation problems. Moreover, a return mapping strategy is proposed to regulate the stress state to avoid calculation failure. Furthermore, a closed-form solution to predicting the failure surface implicitly incorporated in the hypoplastic model is derived. The proposed SPH method is examined with element tests such as the oedometer test and biaxial compression test to observe a good performance of the SPH model. Finally, the SPH model is applied to numerically reproducing sand column collapse to investigate the effect of the initial void ratio on the behavior of soil.

Plowing mechanism of rapid flow-like loess landslides: Insights from MPM modeling

Plowing refers to the phenomenon by which a landslide displaces and pushes the path material forward. Despite its importance in shaping the dynamics of many rapid flow-like loess landslides, the mechanism of plowing is not fully understood. Therefore, this paper presents an improved numerical model based on the material point method (MPM). This model takes into account the interactions between the landslide mass and the path material. The reliability of this model was validated by simulating two granular column collapse experiments involving interaction with an erodible mass. The model was then applied to analyze the plowing dynamics of the Ximiaodian loess landslide in Jingyang County, Shaanxi Province, China. A series of well-designed simulations were conducted to investigate this specific landslide event, contributing to a deeper understanding of the plowing phenomenon in loess landslides. The simulation results show that our model can satisfactorily simulate the plowing process of the Ximiaodian landslide. The plowing process of a rapid-flow loess landslide is proposed to consist of four distinct stages: 1) the collision stage, when the landslide strikes the terrace, resulting in the bending of its strata; 2) the shear failure stage, characterized by significant plastic shear deformation within the bent strata, resulting in the formation of multiple shear failure planes; 3) the steady plowing stage, during which the deformation structure in the affected terrace stabilizes and plowing continues in a steady pattern; and 4) the stopping stage, when the propagation of the disturbed terrace gradually stops. Furthermore, the disturbed terrace zone exhibits three distinct deformation structures. In the rear part, a stable zone develops where the terrace layers are compressed. However, due to the compressive forces exerted by the overlying loess, no significant shear deformation is observed. The middle part is the shear failure zone, characterized by the formation of multiple shear failure planes. The front part is the upheaval zone, where the terrace layers are slightly bent and tilted.

Volcanic Skies: coupling explosive eruptions with atmospheric simulation to create consistent skyscapes

AbstractExplosive volcanic eruptions rank among the most terrifying natural phenomena, and are thus frequently depicted in films, games, and other media, usually with a bespoke once‐off solution. In this paper, we introduce the first general‐purpose model for bi‐directional interaction between the atmosphere and a volcano plume. In line with recent interactive volcano models, we approximate the plume dynamics with Lagrangian disks and spheres and the atmosphere with sparse layers of 2D Eulerian grids, enabling us to focus on the transfer of physical quantities such as temperature, ash, moisture, and wind velocity between these sub‐models. We subsequently generate volumetric animations by noise‐based procedural upsampling keyed to aspects of advection, convection, moisture, and ash content to generate a fully‐realized volcanic skyscape. Our model captures most of the visually salient features emerging from volcano‐sky interaction, such as windswept plumes, enmeshed cap, bell and skirt clouds, shockwave effects, ash rain, and sheathes of lightning visible in the dark.

The 15 January 2022 Hunga (Tonga) eruption: A gas-driven climactic explosion

An extraordinarily powerful, explosive eruption occurred from Hunga volcano in the Tonga island arc on 15 January 2022 and generated an eruption column 58 km high. The explosive eruption also generated atmospheric gravity waves, extreme runup tsunamis and quite unusual and destructive meteotsunamis. Together these place this VEI 6 eruption as, globally, one of the largest of the past 300 years.Based on the oceanic context of Hunga volcano, it has previously been assumed that the eruption was phreatomagmatic through a fuel-coolant Surtseyan-type interaction, but this is not supported by satellite imagery. Similarly, it has been suggested that a caldera-collapse was the eruption trigger, but this is not supported by bathymetric data or the seismicity recorded during the eruption. Here we develop a new model based on the observed energetics and time sequence of the eruption integrated with understanding of the internal structure of active volcanoes and their characteristic high flux discharges of volcanic gas.It has been shown elsewhere that magma-derived reactive gases (H2O, CO2, SO2, HCl, etc) aggressively alter the volcanic rocks in the core of a volcano leading to self-sealing of gas flow to the surface and consequent changes to deviatoric stress in the structure. Common minerals developed by these reactions include anhydrite (CaSO4), sulphides and silica (quartz), all of which have been recorded in volcanic ejecta including at Hunga.We here develop a first order numerical model that quantifies how the free discharge of such gas to the surface may progressively become choked by these sealing reactions leading to increased internal gas pressure. Hydraulic fracture of the seal occurs when the transmitted pressure of the compressed magmatic gas beneath the seal increases to a value greater than the lithostatic pressure plus the tensile strength of the sealed rock. This initiates the explosive release of compressed gas whose high-power discharge progressively develops and enlarges a crater. At the same time, the explosion feeds upon itself by generating larger pressure gradients in the pressurized gas within the fractured porous rock mass of the core of the volcano. Excavation of the crater may intersect high level intrusions and produce the pumice rafts that were observed after the eruption. The eruption itself diminished in intensity as the gas pressure in the reservoir declined.At Hunga, the eruption excavated an 850 m deep, 2-3 km diameter steep-walled crater. This volume may be assumed to approximate the volume of fractured porous rock (the control volume of the eruption) whose trapped gas was mined by the eruption until surrounding gas pressure was depleted. Our numerical model shows that the calculated potential energy of the trapped compressed gas matches the independent observations of the scale of the eruption. Sensor data have since shown that gas bubble flares continued for at least 6 months after the eruption indicating continued depletion of the gas reservoir of rocks surrounding the new crater. The systems-based, gas-driven model for the Hunga climactic eruption developed here also applies to Plinean-type eruptions on subaerial arc volcanoes such as at Pinatubo (Philippines) 1991.

Quantifying tidal versus non-tidal stresses in driving time-varying fluxes of Enceladus' plume eruptions

In this study we investigate the relative importance of endogenic vs. exogenic stresses in controlling the location and timing of active ice-shell deformation on Enceladus, which is expressed by cyclic plume eruptive flux along active fault zones (i.e., the tiger stripes). Although the variation of the eruption flux on Enceladus follows the periodicity of the diurnal tide, it remains unclear why there is a consistent phase delay of the observed peak eruption when compared to the predicted peak tidal stress. Here, we explore whether endogenic stresses in the ice shell are capable of explaining this observed phase delay. To achieve this goal, we performed geologic mapping along the tiger-stripe faults that host the erupting plumes. Using the fault kinematics established from our mapping, we determine the general stress state (i.e., the principal-stress directions) along the tiger-stripe faults. This knowledge in turn forms the basis for inferring the most likely plume-eruption mechanism. Our mapping shows that the tiger-stripe fractures are not tensile cracks but are instead left-slip fault zones locally displaying extensional fissures. This insight leads to a hypothesis that strike-slip faults and their local tensile cracks experience simultaneous shear and tensile failure, and that the tensional opening reaches maximum at the time of the peak plume flux. We quantified this hypothesis using a stress decomposition model that assesses (1) the relative importance in magnitude between the tectonic stress and tidal stress exerted on the tiger-stripe faults and (2) the role of ice shell properties such as the shear strengths, tensile strengths, and ice shell thickness in controlling the eruption phase delay. Using laboratory-determined ice strengths and the best estimate of the ice shell thickness at the South Polar Terrain of Enceladus, which hosts the tiger-stripe faults, our model results indicate that the endogenic tectonic stress is comparable in magnitude to the tidal stress. Although we cannot rule out warm-ice convection, true polar wander, and non-synchronous rotation as causes of endogenic stresses, the large variation in ice shell thickness makes the lateral gravitational-potential gradient the most plausible source of the endogenic stress required by our model results.

Registration of the atmospheric effect of the Hunga Tonga volcano eruption

The paper presents the results of recording of acoustic waves, caused by the Hunga Tonga volcano eruption in the South Pacific Ocean on January 15, 2022, in Eastern Siberia at a distance of about 11230 km from the eruption. The received acoustic signal is interpreted as a set of atmospheric waves in a wide range of oscillations. The structure of the signal is similar to signals from the previously known powerful sources: the thermonuclear explosion on Novaya Zemlya in 1961 and the explosion of the Tunguska meteorite in 1908. The acoustic signal was preceded by three trains of low-frequency damped oscillations. We assume that these three trains of oscillations are associated with three important stages in the Hunga Tonga volcano eruption: 1) destruction of Tonga island and formation of an underwater caldera; 2) release of hot magma from the caldera to the ocean surface and release of a large volume of superheated steam into the atmosphere 3) formation of a layered structure from a mixture of superheated steam, ash, and tephra on the ocean surface and formation of an eruptive convective column. Successive phases of the eruption might have contributed to the excitation of acoustic vibrations in a wide range of periods including Lamb waves, internal gravity waves (IGW), and infrasound. We compare the structure of the acoustic signal received in Siberia at a distance of more than 11000 km from the volcano and that of the acoustic signal recorded in Alaska at a distance of more than 9300 km. Using the solution of the linearized Korteweg — de Vries equation, we estimate the energy released during the volcanic eruption.

О регистрации атмосферного эффекта извержения вулкана Хунга-Тонга

The paper presents the results of recording of acoustic waves, caused by the Hunga Tonga volcano eruption in the South Pacific Ocean on January 15, 2022, in Eastern Siberia at a distance of about 11230 km from the eruption. The received acoustic signal is interpreted as a set of atmospheric waves in a wide range of oscillations. The structure of the signal is similar to signals from the previously known powerful sources: the thermonuclear explosion on Novaya Zemlya in 1961 and the explosion of the Tunguska meteorite in 1908. The acoustic signal was preceded by three trains of low-frequency damped oscillations. We assume that these three trains of oscillations are associated with three important stages in the Hunga Tonga volcano eruption: 1) destruction of Tonga island and formation of an underwater caldera; 2) release of hot magma from the caldera to the ocean surface and release of a large volume of superheated steam into the atmosphere 3) formation of a layered structure from a mixture of superheated steam, ash, and tephra on the ocean surface and formation of an eruptive convective column. Successive phases of the eruption might have contributed to the excitation of acoustic vibrations in a wide range of periods including Lamb waves, internal gravity waves (IGW), and infrasound. We compare the structure of the acoustic signal received in Siberia at a distance of more than 11000 km from the volcano and that of the acoustic signal recorded in Alaska at a distance of more than 9300 km. Using the solution of the linearized Korteweg — de Vries equation, we estimate the energy released during the volcanic eruption.