Pyroclastic Flows Research Articles

2D Geological Structure Identification in Mount Ciremai Geothermal Area using the GGMplus Data.

West Java’s Mount Ciremai holds substantial geothermal potential, marked by hot spring manifestations, which can be an alternative energy source. This research utilizes the GGMplus satellite gravity method, employing First Horizontal Derivative (FHD), Second Vertical Derivative (SVD), and 2D inversion to analyze the geothermal area’s geological structure and subsurface conditions. The results reveal fault structures, including reverse and normal fault types in Mount Ciremai. Specific details for both the Pejambon and Sangkanhurip lines and the Pejambon and Ciniru lines indicate distinct layers with varying densities (2.4 – 2.7 g/cm3) and depths. These layers comprise pyroclastic flow deposits, andesite-basaltic rock, caprock (clay rock), and a reservoir (sandstone), shedding light on the geothermal potential of this region.

A Hidden Eruption: The 21 May 2023 Paroxysm of the Etna Volcano (Italy)

On 21 May 2023, a hidden eruption occurred at the Southeast Crater (SEC) of Etna (Italy); indeed, bad weather prevented its direct and remote observation. Tephra fell toward the southwest, and two lava flows propagated along the SEC’s southern and eastern flanks. The monitoring system of the Istituto Nazionale di Geofisica e Vulcanologia testified to its occurrence. We analyzed the seismic and infrasound signals to constrain the temporal evolution of the fountain, which lasted about 5 h. We finally reached Etna’s summit two weeks later and found an unexpected pyroclastic density current (PDC) deposit covering the southern lava flow at its middle portion. We performed unoccupied aerial system and field surveys to reconstruct in 3D the SEC, lava flows, and PDC deposits and to collect some samples. The data allowed for detailed mapping, quantification, and characterization of the products. The resulting lava flows and PDC deposit volumes were (1.54 ± 0.47) × 106 m3 and (1.30 ± 0.26) × 105 m3, respectively. We also analyzed ground-radar and satellite data to evaluate that the plume height ranges between 10 and 15 km. This work is a comprehensive analysis of the fieldwork, UAS, volcanic tremor, infrasound, radar, and satellite data. Our results increase awareness of the volcanic activity and potential dangers for visitors to Etna’s summit area.

Magma intrusion process during pre-magmatic period (2010−2013) of Sinabung volcano as revealed by seismicity of volcano-tectonic and hybrid earthquakes

Sinabung volcano in Indonesia reactivated, with a series of phreatic eruptions, in August and September 2010. These eruptions were followed by various types of magmatic eruption at the summit, including lava dome growth, pyroclastic density currents, lava flow, and vulcanian eruptions, starting in December 2013. Prior to the magmatic eruptions, Sinabung exhibited two sequences of phreatic eruption and an increase in the seismicity of volcanic earthquakes. This study clarifies the progress of magma intrusion during the pre-magmatic period based on the hypocenter distribution, waveform similarity, seismic moment, rupture length, and stress drop for volcano-tectonic (VT) and hybrid earthquakes. It was found that the hypocenters of VT earthquakes are distributed north and northwest of the summit at a depth range of 1–10 km below sea level. Deep seismicity (4–10 km) alternated with shallow seismicity (2–4 km). The last shallow seismicity, from July to mid-December 2013, was different from previous seismicity in terms of higher intensity, as determined from the seismic moment, a temporary increase in the rupture length, and the migration of hypocenters towards the summit. The seismicity of VT earthquakes was altered by a swarm of hybrid earthquakes in mid-December. The hybrid earthquakes were smaller, in terms of seismic moment (mostly <3 × 1010 Nm), than the VT earthquakes. Their hypocenters were concentrated in the shallowest depth range (−0.5 to 1.5 km below sea level) directly below the summit, their source process was repeatable (they could be grouped into six earthquake families), their dominant peak was in a low-frequency range (2.5–4.5 Hz), and they had a relatively low stress drop (<0.15 MPa). This suggests that the swarm of hybrid earthquakes was induced by a frequent repeated fracture of fluid-filled cracks due to the intrusion of magma up to directly below the summit. The transition of the earthquake family and the change in its source parameters, namely an increase in the stress drop prior to the appearance of a lava dome and then a slight decrease, may reflect a gradual change in internal pressure in the hypocentral zone through the magma intrusion process.

Stratigraphy and structure of Chachahuén volcanic complex, southern Mendoza province, Argentina

The Chachahuén Volcanic Complex is an eroded Upper Miocene arc to back-arc volcanic system that was emplaced in the north-eastern edge of the Neuquén Basin, southern Mendoza province, Argentina. It was formed as the result of shallowing followed by steepening of the Nazca plate during Miocene-Pliocene, and is characterized by a volcanism which ranges from trachydacites to basalts. Thanks to erosion, the stratigraphy and the shallow plumbing systems of the volcano are well-exposed. A new evolution for the volcano stratigraphy is proposed, based on new field observations, satellite imagery, and incorporating the radiometric and petrological data of previous works. The stratigraphy is constrained around one main explosive event, recorded by a thick dacitic pyroclastic density current (PDC) deposit (the Corrales Ignimbrite) (>150m). The main units forming the Chachahuén Volcanic Complex are (1) the pre-dacite PDC deposits (trachydacitic to rhyolitic) Vizcachas Formation, (2) the thick Corrales Ignimbrite, (3) a post-dacite trachyandesitic PDC deposits (dominated by block and ash (BAF) deposits) and thick lava flows, and (4) thin mafic basaltic lava flows. The collapse of a large elliptical caldera occurred during the Trachyandesite phase. The main magma transport channels and feeders are sub-vertical dykes, which exhibit a radial distribution centred on the main caldera depression. Moreover, cryptodomes mainly preserved outside the rim of the caldera accommodated the transport and emplacement of the more silicic magma during the Trachyandesite phase.

SEATANI: hazards from seamounts in Southeast Asia, Taiwan, and Andaman and Nicobar Islands (eastern India)

Abstract. Submarine volcanism makes up approximately 85 % of volcanism that occurs on Earth, and its eruptions have the potential to cause several hazards including ash dispersal, pumice rafts, pyroclastic density currents, sector collapses, and tsunamis. Recent examples include the eruptions in Japan and in the Kingdom of Tonga in 2021 and 2022 respectively, but there has been little to no study of submarine volcanism in Southeast Asia and surroundings. Here we provide a compilation of 466 seamounts from the region, from different published sources, through the SEATANI dataset (Southeast Asia, Taiwan, and the Andaman and Nicobar Islands). We use this newly compiled dataset to assess on a regional level the seamount hazard potential and exposure potential as a springboard for future more quantitative hazard studies of the region. The hazard potential was assessed through seamount morphological and structural analyses, to determine the seamount evolution stage and grade of maturity. The exposure potential was evaluated with two different approaches: an areal analysis of the number of assets within a 100 km radius of each seamount and the development of a hazard-weighted seamount density map to highlight potential areas of interest for future more-in-depth studies. Our results show that there are several potentially hazardous seamounts in this region. Taiwan has the highest hazard and exposure potential, for all assets considered, while the Philippines, Indonesia, and Vietnam have relatively high exposure potential for submarine communication cables and ship traffic density. The results from this work serve as a first step towards Southeast Asia and neighbouring countries becoming more resilient against and prepared for submarine volcanic eruptions in the region.

Facies architecture, geochemistry and petrogenesis of Middle Triassic volcaniclastic deposits of Mt. Ivanščica (NW Croatia): evidence of bimodal volcanism in the Alpine-Dinaridic transitional zone

During the Middle Triassic, intensive volcanic activity took place along the eastern margin of Pangea, including the Greater Adria promontory, due to the Neotethyan oceanization. This resulted in the formation of various volcanic and volcaniclastic rock types. The region of NW Croatia, acting as a transition zone between the Southern Alps and the Dinarides, showcases the outcrops of these rocks. The present study investigates the facies of volcaniclastic rocks, the distribution of those facies, formation processes, as well as the genesis of the primary magma to gain a better understanding of the complex geodynamics of this region during the Middle Triassic. Six profiles across the Vudelja quarry front were surveyed using drone imaging and samples were collected for detailed petrographic and geochemical analyses. Two groups of volcaniclastic rocks were identified—mafic and intermediate/felsic. The former is represented by (I) autoclastic effusive facies and (II) resedimented autoclastic facies, while the latter is represented by (III) secondary pyroclastic facies. Mafic volcaniclastics were generated through basaltic effusions in marine environments, fragmentation in contact with seawater, mixing with shallow marine carbonate clasts, and subsequent redeposition in deeper marine areas. The secondary pyroclastic facies (III) consists of a regionally distributed felsic Pietra Verde tuff whose deposits may be related to pyroclastic density currents and syn-eruptive resedimentation by turbidite-like currents. Geochemical data indicate that parental magmas responsible for generating the mafic volcaniclastics had a calc-alkaline composition and originated in ensialic and mature arc settings of an active continental margin. The observed chemical composition is likely inherited from older, arc-related lithologies, associated with the subduction of the Paleotethys Ocean. Parental magmas are thought to have formed during continental rifting of the passive Middle Triassic margins of the Greater Adria through (i) partial melting of the heterogeneous lithospheric mantle, which had been metasomatized during an earlier Hercynian subduction, and (ii) subordinate processes related to the melting of the upper continental crust and subsequent fractionation. Ar/Ar dating on plagioclase separates yielded an age of 244.5 ± 2.8 Ma for mafic volcaniclastics. This aligns well with biostratigraphic ages of felsic tuffs which crop out on a broader regional scale of the Dinarides, the Southern Alps, and the Transdanubian Range. The overlapping ages obtained from radiometric dating of mafic volcaniclastics and biostratigraphic ages of the felsic Pietra Verde tuffs strongly suggest that the Greater Adria region experienced concurrent bimodal volcanism during the Middle Triassic.

Lahar events in the last 2000 years from Vesuvius eruptions – Part 1: Distribution and impact on densely inhabited territory estimated from field data analysis

Abstract. Lahars represent some of the most dangerous phenomena in volcanic areas for their destructive power, causing dramatic changes in the landscape with no premonitory signs and impacting the population and infrastructure. In this regard, the Campanian Plain turns out to be very prone to the development of these phenomena, since the slopes of the Somma–Vesuvius and Campi Flegrei volcanoes, along with the Apennine reliefs, are mantled by pyroclastic deposits that can be easily remobilized, especially after intense and/or prolonged rainfall. This study focuses on the analysis of pyroclastic fall and flow deposits and of the syn- and post-eruptive lahar deposits related to two sub-Plinian eruptions of Vesuvius in 472 CE (Pollena) and 1631. To begin with, historical and field data from the existing literature and from hundreds of outcrops were collected and organized into a database, which was integrated with several new pieces of data. In particular, stratigraphic, sedimentological (facies analysis and laboratory), and archeological analyses were carried out, in addition to rock magnetic investigations and impact parameter calculations. The new data are also referenced to the finding of ash beds in more distal areas, which were included in new isopach maps for the two sub-Plinian eruptions. The results show that for both eruptions the distribution of the primary deposits is wider than previously known. A consequence of these results is that a wider areal impact should be expected in terms of civil protection, as the sub-Plinian scenario is the reference one for a future large eruption of Vesuvius. Such a distribution of the pyroclastic deposits directly affects the one of the lahar deposits, also because a significant remobilization took place during and after the studied eruptions, which involved distal phreatomagmatic ash. From these integrated analyses, it was possible to constrain the timing of the deposition and the kind of deposits remobilized (pyroclastic fall vs. flow), and it was possible to calculate the velocities and dynamic pressures of the lahars and ultimately infer the lahar transport and emplacement mechanisms. The multidisciplinary approach adopted in this work shows how it is crucial to assess the impact of lahars in densely populated areas even at distances of several to tens of kilometers from active volcanoes. This especially applies to large parts of the densely populated areas around Somma–Vesuvius up to the nearby Apennine valleys.

Ash Fluidity and Regulation Mechanism of Iron-Rich Low-Rank Coal in High-Temperature Gasification Process

ABSTRACT To realize the high-efficiency utilization of iron-rich coal, it is necessary to conduct research on ash flow characteristics and regulatory mechanisms. The ash flow test was carried out, and we introduced additives such as Huajin red mud (HJ), pipe scale (GDG), diatomite (GZT), and Naomaohu coal (NMH). The results revealed that the elevated Fe2O3 content triggers the formation of metallic iron, leading to poor high-temperature fluidity. NMH and HJ additives have the ability to depolymerize SiO2 and Al2O3 structures. This generates higher contents of anorthite and sodium nepheline. The inclusion of GDG ash promoted the formation of cordierite and forsterite, contributing to improved slag flow. Caution is warranted when using a higher proportion of GDG ash additive, as excessive MgO can combine with Al2O3 to form a high-melting-point spinel mineral. GZT ash leads to an increase in quartz content. During cooling, mullite preferentially precipitates from the liquid phase, further reducing high-temperature fluidity.

First evidence of Pinaceae and Fagaceae in the fossil wood record of the České středohoří Mts. (Czech Republic): A comprehensive study of fossiliferous sites in pyroclastic rocks surrounding the late Oligocene Milá stratovolcano

A comprehensive anatomical and mineralogical study of fossil wood fragments from fields in the vicinity of Bečov and Břvany villages (NW Bohemia, Czech Republic) indicates that Taxodioxylon gypsaceum (Cupressaceae s.l.) predominates, but also identifies another coniferous wood: Pinuxylon parryoides (Pinaceae) and three angiosperms Quercoxylon böckhianum, Castanoxylon bavaricum and Lithocarpoxylon sp. (all Fagaceae). This paper therefore presents the first occurrence of Pinaceae and Fagaceae fossil wood in the volcanic rocks of the České Středohoří Mts. as well as its youngest palaeobotanical record in general, late Oligocene in age (26.56 ± 0.38 Ma). The samples were buried by alkaline pyroclastic deposits and were mineralized by carbonates. Two distinct depositional processes burying the fossil woods were identified. Closer to the vent, the woods occur in a near-vent pyroclastic fall deposits of the former pyroclastic cone, whereas more distant sites consist of pyroclastic flow deposits. Carbonate mineralization mostly consists of dolomite, but subordinate amounts of magnesite (likely the first time this is documented in fossil wood) as well as calcite and siderite are present. Only one sample collected in the same area, bearing clear signs of riverbed transport (Lithocarpoxylon sp.), was perfectly silicified, but its origin remains unclear.

A new experimental and numerical simulation method for measuring ash slag fluidity using VOF based on OpenFOAM

A numerical computational model of heated ash flow is proposed using the Volume of Fluid (VOF) method based on OpenFOAM. Established numerical model can accurately calculate the critical flow temperature (Tc) and critical flow velocity (Vc) of ash slag with the error of 1.5 % and 19 % from the experiment, respectively. Forty-five types of ash slag with widely varied properties in China were calculated based on the established model. Results show that high levels (>30 wt%) of Fe2O3 and CaO improve the flowability of the ash slag, and low levels of Al2O3 (10 wt%) also contribute to the flowability of the ash slag. Tc and Vc of the ash were distributed in the range of 1150 ∼ 1400 °C and 0.1 ∼ 1.0 mm/min, respectively. Tc and Vc change in opposite trends with the component content increases, indicating that a higher Tc corresponds to a lower Vc and fluidity reduction. Mathematical models for calculating the Tc and Vc of ash were obtained based on non-linear multiple regression and correlation coefficient, and prediction accuracy was within 10 % and 30 % for Tc and Vc, respectively. The proposed numerical simulation method and prediction model can provide an accurate reference for a wide range of ash fluidity. Also, it was demonstrated that elemental Na has a positive effect on ash flowability based on the model, which is important for measuring slag flowability in slag-tap furnaces.

Fluidization in pyroclastic flows

Pyroclastic flows are mixtures of gas and particles generated by gravitational collapse of lava domes or fall-back of eruption columns during explosive eruptions. They are major natural hazards because they can propagate far away from the source at speeds that largely exceed those of most natural or anthropogenic granular flows. Hazard mitigation requires improved understanding of their complex dynamics to predict flow path, velocity, and run-out distance.The relevance of fluidization of granular solids in pyroclastic flows to their enhanced mobility has received substantive support. Nonetheless, broad uncertainties still characterize fundamental aspects like the prevailing nature and source of the fluidizing gas, the role of particle-phase stresses, turbulence and shear, the effect of solids polydispersity and fragmentation. Moreover, the dynamics of the frontal zone of the flow deserves deeper investigation due to its relevance to air-cushioning/entrainment and the establishment of self-sustained motion-induced fluidization.

Stratigraphy and depositional history of the Aguja Formation (Upper Cretaceous, Campanian) of West Texas, southwestern USA

Abstract Although the Aguja Formation (West Texas, southwestern USA) and its fossil vertebrate fauna have been known for over a century, its basic stratigraphic requisites (type area and type section) have not been formally documented. The formation is herein subdivided into a series of formal members, and a lectostratotype section is proposed. Lithostratigraphic and biostratigraphic subdivisions are documented and integrated with geochronologic data to provide an age model for the formation. Four terrestrial vertebrate biozones are proposed. There are at least four major depositional intervals represented in the Aguja and intertonguing Pen Formations. An initial progradational deltaic succession is recorded by the La Basa Sandstone and lower part of the Abajo Shale Members of the Aguja Formation. A second phase of deposition resulted in a retrogradational shoreface succession that includes the upper part of the Abajo Shale, overlying Rattlesnake Mountain Sandstone Member, and lower part of the McKinney Springs Tongue of the Pen Formation, up to a skeletal phosphate bed interpreted to represent the maximum flooding surface. The third phase of deposition comprises a progradational deltaic succession that includes the upper part of the McKinney Springs Tongue, Terlingua Creek Sandstone Member of the Aguja Formation, and lower part of the Alto Shale Member of the Aguja Formation. This third succession records eastward migration of the strandline and withdrawal of the Western Interior Seaway from the Big Bend region. The fourth phase of deposition comprises a series of aggradational fluvial channel and floodplain successions that form the upper part of the Alto Shale Member and is coincident with redirection of stream flow to the southeast. This interval is much thicker in the central part of the Big Bend region, thins to the southwest and northeast, and likely records initial subsidence in the Laramide Tornillo Basin. The upper part of this succession was also contemporaneous with a series of basaltic pyroclastic eruptions, the westernmost expression of the Balcones igneous province. A dramatic constriction in the southern entrance to the Western Interior Seaway through the Gulf of Mexico occurred during this final phase in deposition of the Aguja Formation and corresponds to a shift of stream flow southeastward and to an outbreak of local pyroclastic eruptions. Regional uplift associated with this episode of magmatism is likely responsible for closing the southern aperture of the Western Interior Seaway.

Tsunamis generated by pyroclastic flows: experimental insights into the effect of the bulk flow density

Tsunamis generated by pyroclastic flows: experimental insights into the effect of the bulk flow density

Challenges in block-and-ash flow hazard assessment: The July 10–11, 2015 eruption of Volcán de Colima, Mexico

On July 10 and 11, 2015, Volcán de Colima, one of the most active Mexican volcanoes, underwent its largest eruption since the 1913 Plinian event. Highly mobile block-and-ash flows (BAF) originated from the dome collapse and subsequent extrusion of bubble gas-rich magma from the conduit. The volume and runout of flows were at least one order of magnitude higher (106 m3) than the events observed previously in the 2004–2005 eruptive episode (105 m3) and only comparable with the pyroclastic density currents emplaced during the Vulcanian eruptive phase that preceded the 1913 Plinian eruption. The TITAN2D code is here used to understand the factors controlling the mobility of the dense BAF undercurrents emplaced during the 2015 eruption. The selection of different input parameters was based on field data, on the seismic signal of the BAFs recorded by a broadband station located a few meters away from the channel, and by considering the different eruptive mechanisms of the July 10 eruptive phase (collapsing dome) with respect to the July 11 event (rapid extrusion of degassed magma directly from the conduit). A collapsing pile best fits the July 10 BAF, with an initial spread on the proximal reach and subsequent emplacement along the Montegrande ravine. Instead, for the July 11 BAF, an influx source best replicates the observed distribution, including the overspilling down the San Antonio ravine at the middle reach and the aggradation at the final depositional fan. Both events were better reproduced with the Voellmy-Salm rheology with respect to the Mohr-Coulomb model. Still, the mobility of the July 11 BAF was controlled mainly by the channel reduced capacity due to the 10 July BAF infilling. This was assessed by modeling the 11 July BAF on the real terrain and on a DEM whose altitude was modified based on the maximum flow depth obtained from the simulation of the former flow. Deep knowledge of the eruptive mechanism and flow behavior during the 2015 eruption allowed us to understand the main factors controlling concentrated pyroclastic density current mobility during sustained, multi-phase eruptions, especially for valley-confined BAF. The results presented here set the basis to improve the hazard zonation for pyroclastic density currents at Volcán de Colima, which should be based on a multifactor, probabilistic analysis.

Rainfall-induced Gravity Movement of the Unzen-Fugendake Volcanic Dome Analysis combining Ground-Radar Interferometry and XRAIN Rainfall radar system

Lava domes created by volcanoes often cause pyroclastic and debris flows, which have a significant impact on the surrounding infrastructure and population and have been the subject of much research. However, because volcanic domes tend not to survive the eruptions that form them, the instability of domes that survive eruptions such as Unzen Fugendake in Japan is both a poorly understood process and a danger. Therefore, the present contribution aims to (1) Displacement and precipitation from 2018 to 2020 for lava domes at Unzen volcano and their relationship to earthquakes, and (2) Haar wavelet analysis to understand the response of displacement to precipitation. The method is based on dome displacement from the Unzen Ground Based Synthetic Aperture Radar system and 48-hour rainfall from the MP radar rain gauge system. As a result, the authors confirmed the following: (1) precipitation of 150 mm or more in 48 hours tends to increase the vertical fluctuation of the dome, and even rainfall of less than 100 mm per 48 hours has a similar effect when it is repeated in an intensive manner; (2) After precipitation, major dome displacement can take days or weeks to occur, and is not instantaneous like the dome collapses in Soufriere and Merapi.

Volcano generated tsunami recorded in the near source

Volcano sector collapse and pyroclastic density currents are common phenomena on active volcanoes and potentially a fatal source of tsunami waves which constitute a serious hazard for local as well as distant coastal population. Several examples in recent history, warn us on the urgent need to improve our mitigation counter-actions when tsunamis have volcanic origin. However, instrumental record of tsunami generated by mass movement along a volcano flank are still rare and not well understood yet. Small tsunamis (≤1 m) induced by pyroclastic density currents associated to violent explosions of Stromboli volcano were recorded in near-source conditions (<1.6 km). We show how tsunami waveform remains unaltered regardless of the two orders of variability in the landslide volume and dynamics. This unprecedented record is also providing the lesson to develop unconventional warning strategies necessary when the tsunamigenic source is expected to be very close (<10 minutes) to densely populated coasts and with a limited time to issue an alert based on simulation of wave propagation and inundation.

Complex facies architecture of unusual pyroclastic density current deposits at the Rincón de la Vieja Volcano, Costa Rica

Pyroclastic Density Currents (PDCs) constitute a complex, multiphase, and multiscale volcanic process generated in a wide spectrum of eruptive scenarios and being among the most deadly and destructive volcanic hazards. Over the north flank of the Rincón de la Vieja Volcano, Costa Rica, previous works identified PDC deposits related to the plinian Río Blanco eruption, which is the largest Holocene explosive event recorded at this volcano. However, so far, there have been no in-depth studies aimed at understanding the eruptive style, and dynamics of the parental PDCs. This work presents, for the first time, a detailed facies analysis and architecture reconstruction (plus a geochemical characterization and novel C14 age determinations) of these outstanding deposits. The presence of cauliflower bombs, juvenile mega-blocks, scorias, and pumices, coupled with the scarcity of accidental and/or accessory lithics, is interpreted here as evidence of an unusual eruptive style. We propose that the Río Blanco eruption was a single complex eruption consisting of two successive phases, each involving the drainage of distinct magma batches. The PDCs were related to the sudden spill out of a lava-pool (Arenal-type PDC), during the collapse of an ultraproximal pyroclastic cone that was built in an initial plinian phase. The collapse of this ultraproximal cone produced a leading-edge pumice-rich PDC, which increased its mobility through fragmentation-induced fluidization. This leading edge PDC was followed by successive high concentrated, mafic, modified granular flow pulses derived from the lava-pool spill out. These modified granular flows were highly channelized in proximal zones, but when they reached the main break in slope, the basal underflow froze en masse while its upper, more diluted part started to decouple and move in a more unconfined manner. After the decoupling, the material feeding towards the basal zone of the surviving current produced a new underflow by rapid suspension sedimentation (“surge-derived pyroclastic flows”). Finally, our results indicate that the generation of Arenal-type PDCs is an underestimated scenario that should be considered in future short-term hazard and risk assessments at the Rincón de la Vieja Volcano, as well as in many other stratovolcanoes worldwide.

Modeling SO2 dispersion from future eruptions in the Auckland Volcanic Field, New Zealand

Auckland city (pop. 1.7 M) is Aotearoa New Zealand’s largest city and an important economic hub. The city is built upon the active intraplate basaltic Auckland Volcanic Field (AVF). An AVF eruption would cause considerable impacts. An important component of volcanic risk management is assessing the likely volcanic hazards to help inform emergency planning and other preparedness activities. Previous volcanic hazard assessments for the AVF, particularly those for emergency planning scenarios, have modeled multiple volcanic hazards including lava flows, pyroclastic density currents, ballistic projectiles and tephra fall. Despite volcanic gas being an important and impactful hazard from intraplate basaltic field eruptions, there has been limited consideration of volcanic gas in AVF hazard assessment to date. This project is one of the first to quantitatively assess potential volcanic gas hazards for an explosive eruption scenario. For basaltic volcanism, sulfur dioxide (SO2) gas is typically the most consequential volcanic gas emitted. The aim of this exploratory study was to model SO2 dispersion from a high impact eruption during weather conditions conducive to high ground level pollutant concentrations. Since ground level SO2 concentrations are influenced by complex wind patterns resulting from interactions of locally driven flow circulations and topographically influenced weather, we modeled SO2 dispersion using the HYSPLIT model, a state-of-the art hybrid Eulerian and Lagrangian dispersion model widely used for volcanic gases, using high-resolution meteorological forcing fields given by the Weather Research and Forecasting (WRF) model.Modeled air parcel trajectories and ground level SO2 concentrations illustrate the effect of the converging sea breeze winds on SO2 dispersion. Under worst-case dispersion conditions, extensive areas of up to hundreds of square kilometers to the north and northwest of the eruption location would exceed New Zealand short-term (24 h) air quality standards and guidelines for SO2, indicating heightened health risks to downwind communities. Using this numerical modeling approach, this work presents a methodology for future applications to other AVF eruption scenarios, with a wider range of meteorological conditions that can help in exploring consequences for health services such as anticipated emergency department respiratory admissions.

Submarine pyroclastic deposits from 7.3 ka caldera-forming Kikai-Akahoya eruption

Extensive submarine pyroclastic deposits were discovered around the Kikai caldera, Southwest Japan, which originated from the caldera-forming Kikai-Akahoya (K-Ah) eruption that occurred at 7.3 ka. This submarine pyroclastic unit originated from a pyroclastic density current (PDC) that entered the sea and was transformed into a subaqueous density current. Our high-resolution marine seismic reflection surveys reveal that the uppermost unit exhibiting prominent seismic facies thins radially away from the caldera and covers >4500 km2 of the seafloor. The estimated volume of the submarine deposits is >71 km3 based on the isopach map. Seismic sequence stratigraphic analyses and correlation of shard glass compositions from the piston core and remotely operated vehicle samples with those from subaerial outcrops identified the unit as pyroclastic deposits produced by the K-Ah eruption. The geomorphological features of the unit indicate a long runout distance (i.e., 40 km), exponential thinning with increasing distance from the source, and a depression-filling nature of the deposits suggesting that the unit was formed predominantly by a density current that entered the water. The quantitative comparisons of geometric shape and its thinning rate, including that for non-volcanic density currents, indicate that the subaerial PDC entering into the water was transported as a dilute density current in the water. We propose that the formation of these dilute density currents transformed by PDCs probably requires sufficient water depth, as inferred also from previous analog experiments. Adding on-land ignimbrite of 5 km3 and widespread tephra of 249–374 km3, our estimate of submarine pyroclastic deposits of >71 km3 leads to a total bulk volume of the K-Ah eruption of 332–457 km3 (133–183 km3 in DRE), making it probably the largest eruption of the Holocene.

Constraining Formation Hypotheses for Irregular Mare Patches on the Moon With Orbital Reflectance Spectra

AbstractIrregular mare patches (IMPs) are enigmatic volcanic features on the Moon's surface, whose lack of cratering and crisp appearance imply they formed &lt;100 Ma ago, ∼1 Ga after the expected turnoff of lunar volcanism. Multiple contrasting formation hypotheses have been put forth to explain their young appearance, including recent emplacement via eruptions of juvenile volcanic material or outgassing, versus ancient volcanic deposits that were emplaced billions of years ago but only appeared young due to highly porous material. If IMPs formed recently, this would require a reinterpretation of lunar thermal evolution. To help constrain formation hypotheses, we provide a comprehensive mineralogical analysis of IMPs using visible to near infrared hyperspectral data from the Moon Mineralogy Mapper. IMPs appear spectrally dominated by high‐calcium pyroxene and are spectrally similar to their host mare and fresh craters. IMPs do not show clear indications of significant glass, implying that pyroclastic eruption was not significant in IMP formation. Based on spectral comparisons to terrestrial magma foam, we find that this also contradicts the glassiness expected for a magma foam exposed at the surface. Thus, we find it is unlikely that IMPs are composed of recently erupted material and may instead be the result of recent or ongoing surface modification of materials similar in composition and likely contemporaneously emplaced with the mare. We favor previous hypotheses that collapse processes or drainage into subsurface voids or porous materials may have been the major drivers of IMP surface rejuvenation, supported by their proximity to collapse features.