Under the surface: Pressure-induced planetary-scale waves,volcanic lightning,and gaseous clouds caused by the submarine eruption of Hunga Tonga-Hunga Ha'apai volcano

We present a narrative of the eruptive events culminating in the cataclysmic January 15, 2022 eruption of Hunga Tonga-Hunga Ha'apai Volcano by synthesizing diverse preliminary seismic, volcanological, sound wave, and lightning data available within the first few weeks after the eruption occurred. The first hour of eruptive activity produced fast-propagating tsunami waves, long-period seismic waves, loud audible sound waves, infrasonic waves, exceptionally intense volcanic lightning and an unsteady volcanic plume that transiently reached—at 58 ?km—the Earth's mesosphere. Energetic seismic signals were recorded worldwide and the globally stacked seismogram showed episodic seismic events within the most intense periods of phreatoplinian activity, and they correlated well with the infrasound pressure waveform recorded in Fiji. Gravity wave signals were strong enough to be observed over the entire planet in just the first few hours, with some circling the Earth multiple times subsequently. These large-amplitude, long-wavelength atmospheric disturbances come from the Earth's atmosphere being forced by the magmatic mixture of tephra, melt and gasses emitted by the unsteady but quasi-continuous eruption from 0402±1–1800 UTC on January 15, 2022. Atmospheric forcing lasted much longer than rupturing from large earthquakes recorded on modern instruments, producing a type of shock wave that originated from the interaction between compressed air and ambient (wavy) sea surface. This scenario differs from conventional ideas of earthquake slip, landslides, or caldera collapse-generated tsunami waves because of the enormous (～1000x) volumetric change due to the supercritical nature of volatiles associated with the hot, volatile-rich phreatoplinian plume. The time series of plume altitude can be translated to volumetric discharge and mass flow rate. For an eruption duration of ～12 ?h, the eruptive volume and mass are estimated at 1.9 ?km³ and ～2 900 ?Tg, respectively, corresponding to a VEI of 5–6 for this event. The high frequency and intensity of lightning was enhanced by the production of fine ash due to magma—seawater interaction with concomitant high charge per unit mass and the high pre-eruptive concentration of dissolved volatiles. Analysis of lightning flash frequencies provides a rapid metric for plume activity and eruption magnitude. Many aspects of this eruption await further investigation by multidisciplinary teams. It represents a unique opportunity for fundamental research regarding the complex, non-linear behavior of high energetic volcanic eruptions and attendant phenomena, with critical implications for hazard mitigation, volcano forecasting, and first-response efforts in future disasters.     

桂西南那坡盆地河口组沉凝灰岩的时代与来源及其大地构造意义SCIEI北大核心CSCD

桂西南那坡盆地位于八布-Song Hien构造带,发育早-中三叠世岩浆岩和巨厚海相碎屑岩系,是研究华南地块南缘古特提斯造山过程的关键地区。本文对那坡盆地中三叠统河口组中-上部的两套沉凝灰岩开展了详细的岩石学、地球化学和锆石U-Pb-Hf同位素分析研究。两件沉凝灰岩样品LA-ICP-MS锆石U-Pb定年结果分别为240.6±1.0Ma和236.3±0.8Ma,指示河口组碎屑岩沉积时代为中三叠世拉丁期,部分为晚三叠世卡尼期早期。全岩地球化学结果显示,该沉凝灰岩属于亚碱性英安质,富集轻稀土(LREEs)和大离子亲石元素(LILEs)、Rb、Th和U,亏损高场强元素(HFSEs),具有显著的Nb、Ta和Ti负异常。沉凝灰岩锆石具有较低的ε_(Hf)(t)值(-18.6~-6.3)和较老的地壳模式年龄(t_(DM)^(C)=1664~2434Ma),结合锆石微量元素特征认为其形成于俯冲环境。整体上,那坡盆地河口组沉凝灰岩表现出俯冲作用形成的岛弧火山岩地球化学特征,推断其来源于桂西南凭祥地区的中-晚三叠世酸性火山岩。结合已有区域地质调查成果,认为华南地块南缘八布-Cao Bang古特提斯分支洋盆于晚二叠世至晚三叠世早期向北俯冲,随后自西向东呈“剪刀”式逐渐闭合。华南地块与北越地块完全碰撞拼合的时间发生在晚三叠世中晚期。     

Tectonic evolution of the Late Paleozoic Early Mesozoic orogenic belt in the eastern segment of the northern margin of the North China Block: Evidence from meta volcanic rocks of Jianshanzi,northern Liaoning ProvinceSCIEI北大核心CSCD

LA-ICP-MS U-Pb dating, whole -rock geochemistry and Hf isotopic analyses for the Jianshanzi volcanic rocks in the Kaiyuan area are presented to explore the Late Paleozoic -Early Mesozoic tectonic evolution history of the eastern segment of the northern margin of the North China block. The detailed research indicates that the protolith of Jianshanzi metamorphic volcanic rocks includes rhyolite, dacite and andesite, they are calc-alkaline series and metaluminous-weak peraluminous rocks. According to the zircon U-Pb dating, the metamorphic rhyolite and metamorphic andesite formed in the Late Permian (256. 1 +/- 1. 5Ma, 252.4 +/- 1. 7Ma) and the Late Triassic (222. 4 +/- 1. 8Ma, 227. 9 +/- 1. 2Ma), respectively. They are characteristized by enrichment in light rare earth elements and large ion lithophile elements, and depletion of heavy rare earth elements and high field strength elements. Metamorphic rhyolite has an obvious negative Eu anomaly with epsilon(Hf) ( t) values of 10. 83 + 8. 10 and the second -stage model are 698 similar to 1967Ma. They are mainly due to the partial melting of the juvenile basaltic lower crust and the contamination of a small amount of ancient crustal materials. Metamorphic andesite has no obvious Eu anomaly, and the epsilon(Hf)( t) value and the second -stage model are +4.54 similar to + 6. 79 and 825 similar to 968 Ma, which may be the product of partial melting of basic lower crust. There are two periods of magmatism in Jianshanzi volcanic rocks: Late Permian -Early Triassic (247 similar to 256Ma) and Late Triassic (similar to 225Ma). The material composition of the Jianshanzi volcanic rocks is characterized by different times and different tectonic backgrounds. It is comprehensively believed that the eastern segment of the northern margin of the North China block was in an active continental margin tectonic environment during the Permian, and the Palaeo-Asian Ocean subducted in both directions; Mantle -derived magma ascends and emplaces, forming a nascent crust and causing crustal thickening. The Paleo-Asian Ocean finally closed in the early Middle Triassic, and the impact of collisional orogeny may have continued until the early Late Triassic.     

CO2 output and δ13C(CO2) from Mount Etna as indicators of degassing of shallow asthenosphere

 An estimated average CO₂ output from Etna's summit craters in the range of 13±3 Mt/a has recently been determined from the measured SO₂ output and measured CO₂/SO₂ molar ratios. To this amount the CO₂ output emitted diffusely from the soil (≈ 1 Mt/a) and the amount of CO₂ dissolved in Etna's aquifers (≈ 0.25 Mt/a) must be added. Data on the solubility of CO₂ in Etnean magmas at high temperature and pressure allow the volume of magma involved in the release of such an amount of this gas to be estimated. This volume of magma (≈ 0.7 km³/a) is approximately 20 times greater than the volume of magma erupted annually during the period 1971–1995. On the basis of C-isotopic data of CO₂ collected in the Etna area and of new hypotheses on the source of Mediterranean magmas, significant contributions of CO₂ from non-magmatic sources to the total output from Etna are unlikely. Such large outputs of CO₂ and also of SO₂ from Etna could be due to an anomalously shallow asthenosphere beneath the volcano that allows a continuous escape of gases toward the surface, even without migration of magma. Received: 7 August 1996 / Accepted: 9 November 1996     

The eruptive history of the Tequila volcanic field, western Mexico: ages, volumes, and relative proportions of lava types

The eruptive history of the Tequila volcanic field (1600 km²) in the western Trans-Mexican Volcanic Belt is based on ⁴⁰Ar/³⁹Ar chronology and volume estimates for eruptive units younger than 1 Ma. Ages are reported for 49 volcanic units, including Volcán Tequila (an andesitic stratovolcano) and peripheral domes, flows, and scoria cones. Volumes of volcanic units 1 Ma were obtained with the aid of field mapping, ortho aerial photographs, digital elevation models (DEMs), and ArcGIS software. Between 1120 and 200 kyrs ago, a bimodal distribution of rhyolite (~35 km³) and high-Ti basalt (~39 km³) dominated the volcanic field. Between 685 and 225 kyrs ago, less than 3 km³ of andesite and dacite erupted from more than 15 isolated vents; these lavas are crystal-poor and show little evidence of storage in an upper crustal chamber. Approximately 200 kyr ago, ~31 km³ of andesite erupted to form the stratocone of Volcán Tequila. The phenocryst assemblage of these lavas suggests storage within a chamber at ~2–3 km depth. After a hiatus of ~110 kyrs, ~15 km³ of andesite erupted along the W and SE flanks of Volcán Tequila at ~90 ka, most likely from a second, discrete magma chamber located at ~5–6 km depth. The youngest volcanic feature (~60 ka) is the small andesitic volcano Cerro Tomasillo (~2 km³). Over the last 1 Myr, a total of 128±22 km³ of lava erupted in the Tequila volcanic field, leading to an average eruption rate of ~0.13 km³/kyr. This volume erupted over ~1600 km², leading to an average lava accumulation rate of ~8 cm/kyr. The relative proportions of lava types are ~22–43% basalt, ~0.4–1% basaltic andesite, ~29–54% andesite, ~2–3% dacite, and ~18–40% rhyolite. On the basis of eruptive sequence, proportions of lava types, phenocryst assemblages, textures, and chemical composition, the lavas do not reflect the differentiation of a single (or only a few) parental liquids in a long-lived magma chamber. The rhyolites are geochemically diverse and were likely formed by episodic partial melting of upper crustal rocks in response to emplacement of basalts. There are no examples of mingled rhyolitic and basaltic magmas. Whatever mechanism is invoked to explain the generation of andesite at the Tequila volcanic field, it must be consistent with a dominantly bimodal distribution of high-Ti basalt and rhyolite for an 800 kyr interval beginning ~1 Ma, which abruptly switched to punctuated bursts of predominantly andesitic volcanism over the last 200 kyrs.Electronic Supplementary Material Supplementary material is available in the online version of this article at Editorial responsility: J. Donnelly-NolanThis revised version was published online in January 2005 with corrections to Tables 1 and 3.An erratum to this article can be found at     

Volcanic risk ranking for Auckland, New Zealand. I: Methodology and hazard investigation

Volcanic eruptions typically produce a number of hazards, and many regions are at risk from more than one volcano or volcanic field. So that detailed risk assessments can be carried out, it is necessary to rank potential volcanic hazards and events in terms of risk. As it is often difficult to make accurate predictions regarding the characteristics of future eruptions, a method for ranking hazards and events has been developed that does not rely on precise values. Risk is calculated individually for each hazard from each source as the product of likelihood, extent and effect, based on the parameters order of magnitude. So that multiple events and outcomes can be considered, risk is further multiplied by the relative probability of the event occurring (probability_e) and the relative importance of the outcome (importance_o). By adding the values obtained, total risk is calculated and a ranking can be carried out.This method was used to rank volcanic hazards and events that may impact the Auckland Region, New Zealand. Auckland is at risk from the Auckland volcanic field, Okataina volcanic centre, Taupo volcano, Tuhua volcano, Tongariro volcanic centre, and Mt. Taranaki volcano. Relative probabilities were determined for each event, with the highest given to Mt. Taranaki. Hazards considered were, for local events: tephra fall, scoria fall and ballistic impacts, lava flow, base surge and associated shock waves, tsunami, volcanic gases and acid rain, earthquakes and ground deformation, mudflows and mudfills, lightning and flooding; and for distal events: tephra fall, pyroclastic flows, poisonous gases and acid rain, mudflows and mudfills, climate variations and earthquakes. Hazards from each source were assigned values for likelihood, with the largest for tephra fall from all sources, earthquakes and ground deformation, lava flows, scoria fall and base surge for an Auckland eruption on land, and earthquakes and ground deformation from an Auckland eruption in the ocean. The largest values for extent were for tephra fall and climate variation from each of the distal centres. However, these parameters do not give a true indication of risk. In a companion paper the effect of each hazard is fully investigated and the risk ranking completed.     

Surface area, porosity and water adsorption properties of fine volcanic ash particles

Our understanding on how ash particles in volcanic plumes react with coexisting gases and aerosols is still rudimentary, despite the importance of these reactions in influencing the chemistry and dynamics of a plume. In this study, six samples of fine ash (<100 m) from different volcanoes were measured for their specific surface area, a_s, porosity and water adsorption properties with the aim to provide insights into the capacity of silicate ash particles to react with gases, including water vapour. To do so, we performed high-resolution nitrogen and water vapour adsorption/desorption experiments at 77 K and 303 K, respectively. The nitrogen data indicated a_s values in the range 1.1–2.1 m²/g, except in one case where a a_s of 10 m²/g was measured. This high value is attributed to incorporation of hydrothermal phases, such as clay minerals, in the ash surface composition. The data also revealed that the ash samples are essentially non-porous, or have a porosity dominated by macropores with widths >500 Å. All the specimens had similar pore size distributions, with a small peak centered around 50 Å. These findings suggest that fine ash particles have relatively undifferentiated surface textures, irrespective of the chemical composition and eruption type. Adsorption isotherms for water vapour revealed that the capacity of the ash samples for water adsorption is systematically larger than predicted from the nitrogen adsorption a_s values. Enhanced reactivity of the ash surface towards water may result from (i) hydration of bulk ash constituents; (ii) hydration of surface compounds; and/or (iii) hydroxylation of the surface of the ash. The later mechanism may lead to irreversible retention of water. Based on these experiments, we predict that volcanic ash is covered by a complete monolayer of water under ambient atmospheric conditions. In addition, capillary condensation within ash pores should allow for deposition of condensed water on to ash particles before water reaches saturation in the plume. The total mass of water vapour retained by 1 g of fine ash at 0.95 relative water vapour pressure is calculated to be ~10^–2 g. Some volcanic implications of this study are discussed.Editorial responsibility: J. Gilbert     

The Tambora eruption in 1815 provides a test on possible global climatic and chemical perturbations in the past

A coupled one-dimensional radiative-convective-photochemical diffusion model, which takes into account the influence of ocean inertia on global radiative perturbations is used to investigate the possible climatic and other atmospheric effects of a major volcanic eruption, thought to be similar in magnitude to that of the Tambora eruption, Indonesia, which took place in 1815. A volcanic cloud was introduced in the model stratosphere between 20–25 km and the global average peak aerosol optical thickness was assumed to be 0.25. Both the aerosol optical thickness and aerosol composition, which determine the optical properties, were allowed to vary in the model atmosphere during the life cycle of the volcanic cloud. The results indicate that the global average surface temperature decreases steadily from the date of eruption (7–12 April 1815) with maximum cooling of 1° K occurring in the spring of 1816. The calculations also show significant warming of the stratosphere, with temperature increasing up to 15° K at 25 km in less than six months after the date of eruption. The important effects of the Tambora eruption on stratospheric ozone and UV-B radiation at the surface are also mentioned.     

Evolution of a complex isolated dome system,Cerro Pizarro,central México

Cerro Pizarro is an isolated rhyolitic dome in the intermontane Serdán-Oriental basin, located in the eastern Trans-Mexican Volcanic Belt. Cerro Pizarro erupted ~1.1 km³ of magma at about 220 ka. Activity of Cerro Pizarro started with vent-clearing explosions at some depth; the resultant deposits contain clasts of local basement rocks, including Cretaceous limestone, ~0.46-Ma welded tuff, and basaltic lava. Subsequent explosive eruptions during earliest dome growth produced an alternating sequence of surge and fallout layers from an inferred small dome. As the dome grew both vertically and laterally, it developed an external glassy carapace due to rapid chilling. Instability of the dome during emplacement caused the partial gravitational collapse of its flanks producing various block-and-ash-flow deposits. After a brief period of repose, re-injection of magma caused formation of a cryptodome with pronounced deformation of the vitrophyric dome and the underlying units to orientations as steep as near vertical. This stage began apparently as a gas-poor eruption and no explosive phases accompanied the emplacement of the cryptodome. Soon after emplacement of the cryptodome, however, the western flank of the edifice catastrophically collapsed, causing a debris avalanche. A hiatus in eruptive activity was marked by erosion of the cone and emplacement of ignimbrite derived from a caldera to the north of Cerro Pizarro. The final growth of the dome growth produced its present shape; this growth was accompanied by multiple eruptions producing surge and fallout deposits that mantle the topography around Cerro Pizarro. The evolution of the Cerro Pizarro dome holds aspects in common with classic dome models and with larger stratovolcano systems. We suggest that models that predict a simple evolution for domes fail to account for possibilities in evolutionary paths. Specifically, the formation of a cryptodome in the early stages of dome formation may be far more common than generally recognized. Likewise, sector collapse of a dome, although apparently rare, is a potential hazard that must be recognized and for which planning must be done.Editorial responsibility: J. Gilbert     

Geometry and growth of sill complexes: insights using 3D seismic from the North Rockall Trough

Doleritic sill complexes, which are an important component of volcanic continental margins, can be imaged using 3D seismic reflection data. This allows unprecedented access to the complete 3D geometry of the bodies and an opportunity to test classic sill emplacement models. The doleritic sills associated with basaltic volcanism in the North Rockall Trough occur in two forms. Radially symmetrical sill complexes consist of a saucer-like inner sill at the base with an arcuate inclined sheet connecting it to a gently inclined, commonly ragged, outer rim. Bilaterally symmetrical sill complexes are sourced by magma diverted from a magma conduit feeding an overlying volcano. With an elongate, concave upwards, trough-like geometry bilaterally symmetrical sills climb away from the magma source from which they originate. Both sill complex types can appear as isolated bodies but commonly occur in close proximity and consequently merge, producing hybrid sill complexes. Radial sill complexes consist of a series of radiating primary flow units. With dimensions up to 3 km, each primary flow unit rises from the inner saucer and is fed by primary magma tube. Primary flow units contain secondary flow units with dimensions up to 2 km, each being fed by a secondary magma tube branching from the primary magma tube. Secondary flow units in turn are composed of 100-m scale tertiary flow units. A similar branching hierarchy of flow units can also be seen in bilaterally symmetrical sill complexes, with their internal architecture resembling an enlarged version of a primary flow unit from a radial sill complex. This branching flow pattern, as well as the interaction between flow units of varying orders, provides new insights into the origin of the structures commonly seen within sill complexes and the hybrid sill bodies produced by their merger. The data demonstrate that each radially symmetrical sill complex is independently fed from a source located beneath the centre of the inner saucer, grows by climbing from the centre outwards and that peripheral dyking from the upper surface is a common feature. These features suggest a laccolith emplacement style involving peripheral fracturing and dyking during inner saucer growth and thickening. The branching hierarchy of flow units within bilaterally symmetrical sill complexes is broadly similar to that of primary flow units within a radially symmetrical sill complex, suggesting that the general features of the laccolith emplacement model also apply.Editorial responsibility: J. Stix