Island Volcano Research Articles

Dynamics and Detection of Pulsed Tremor at Whakaari (White Island), Aotearoa New Zealand

AbstractVolcanic tremor is a crucial indicator for assessing the state and hazard potential of volcanic systems. At Whakaari (White Island volcano, Aotearoa New Zealand), a pulsed tremor signal emerged after a hydrothermal explosion in August 2012. The tremor accompanied the extrusion of a lava dome, before gradually disappearing prior to the onset of renewed hydrothermal activity in January 2013. We interpret this seismic signal to represent discrete gas transfers from a magmatic intrusion toward a permeable cap—possibly a hydrothermal seal—in the upper layers of Whakaari's hydrothermal system. Such tremor may thus be associated with heightened potential for hazardous explosive activity but is difficult to detect using conventional seismic monitoring parameters. To highlight the emergence of subtle periodic signals, we experiment with Lomb‐Scargle periodograms (LS). LS detect the tremor 5 days before it becomes visible in seismograms, thus facilitating the recognition of such elusive seismic patterns.

Understanding Multi-Hazard Interactions and Impacts on Small-Island Communities: Insights from the Active Volcano Island of Ternate, Indonesia

Drawing on a case study from Ternate Island, a densely populated volcanic island in Eastern Indonesia, this research illustrates how multi-hazards and extreme weather events are likely to compound and cascade, with serious consequences for sustainable development in small island context. At the heart of Ternate Island sits the active Gamalama volcano, posing a constant eruption threat. Its location within the Ring of Fire further exposes the island to the risks of tsunamis and earthquakes. Additionally, the island’s physical features make it highly susceptible to flooding, landslides, and windstorms. Rapid urbanization has led to significant coastal alterations, increasing exposure to hazards. Ternate’s small-island characteristics include limited resources, few evacuation options, vulnerable infrastructure, and inadequate resilience planning. Combining GIS multi-hazard mapping with a structured survey in 60 villages in Ternate, this case study investigates the multi-hazard exposure faced by the local population and land coverage. The findings suggest significant gaps between village chiefs’ perceptions of the types of hazards and the multi-hazard assessment in each village. Out of 60 village chiefs surveyed, 42 (70%) are aware of earthquake risks, 17 (28%) recognize tsunami threats, and 39 see volcanoes as a danger. GIS assessments show that earthquakes could impact all villages, tsunamis could affect 46 villages (77%), and volcanoes could threaten 39 villages. The hazard map indicates that 32 villages are at risk of flash floods and 37 are at risk of landslides, and extreme weather could affect all villages. Additionally, 42 coastal villages on Ternate Island face potential extreme wave and abrasion disasters, but only 18 chiefs acknowledge extreme weather as a threat. The paper argues that addressing the cognitive biases reflected in the perceptions of community leaders requires transdisciplinary dialogue and engagement.

Myth busting: Was Pulau Tiga really first created by a mud volcano eruption in 1897?

Pulau Tiga is Malaysia’s largest mud volcano, and the entire island is famously reported to have suddenly appeared in September 1897. The story of Pulau Tiga’s “birth” following two large earthquakes in the Philippines is widespread, being found in scientific papers, media stories, Wikipedia, and thousands of websites. This study “fact checks” whether Pulau Tiga really first appeared in 1897. The event is described in detail in one 1986 scientific study, and the idea of Pulau Tiga’s sudden appearance is supported by instances of islands being created by large mud volcano eruptions, such as Pulau Batu Hairan offshore Malaysia in 1988. However, maps of Borneo published prior to 1897 indicate that Pulau Tiga is charted on maps dating back to 1554 and specifically named on more than 40 maps. These maps conclusively prove that Pulau Tiga does not suddenly appear in 1897 and also indicate that Pulau Tiga has had approximately the same shape for centuries, demonstrating that the island was also not partially formed in 1897. A review of the newspapers and scientific reports from 1897 to 1904 reveals that two new islands did emerge offshore northern Borneo on 21 September 1897, but neither of them is Pulau Tiga. New mud volcano islands temporarily appeared off the western Klias Peninsula near Bukit Tomboh and off the northern tip of Sabah at the approximate position of the Pulau Batu Hairan. The origin of Pulau Tiga’s 1897 birth “myth” is proposed to be a 1986 study by McManus and Tate, who appear to have confused reports of an island forming off the Klias Peninsula in 1897 with an eruption of Pulau Tiga. This study definitively “busts” the myth of Pulau Tiga’s 1897 birth and unearths details of Malaysia’s history of large mud volcano eruptions that have been largely forgotten by the geoscience community.

Data and early results from temporary seismic arrays for monitoring and investigating magmatic processes beneath Mt. Halla and Ulleung Island volcanoes, South Korea

Data and early results from temporary seismic arrays for monitoring and investigating magmatic processes beneath Mt. Halla and Ulleung Island volcanoes, South Korea

Iron valence systematics in clinopyroxene crystals from ocean island basalts

The valence state of Fe plays a vital role in setting and recording the oxidation state of magmas, commonly expressed in terms of oxygen fugacity (fO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{\ extrm{O}_{2}}$$\\end{document}). However, our knowledge of how and why fO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{\ extrm{O}_{2}}$$\\end{document} varies within and between magmatic systems remains patchy because of diverse challenges associated with estimating the valence state of Fe in glasses and minerals routinely. Here we investigate Fe valence systematics in magmatic clinopyroxene crystals from ocean island basalts (OIBs) erupted in Iceland and the Azores to explore whether they record information about magma Fe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document} contents and magmatic fO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{\ extrm{O}_{2}}$$\\end{document} conditions. Although many studies assume that all Fe in augitic clinopyroxene crystals from OIBs occurs as Fe2+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{2+}$$\\end{document}, we find that up to half of the total Fe present can occur as Fe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document}, with crystals from alkali systems typically containing more Fe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document} than those from tholeiitic systems. Thus, Fe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document} is a major if under-appreciated constituent of augitic clinopyroxene crystals erupted from ocean island volcanoes. Most Fe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document} in these crystals is hosted within esseneite component (CaFe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document}AlSiO6\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{6}$$\\end{document}), though some may be hosted in aegirine component (NaFe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document}Si2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document}O6\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{6}$$\\end{document}) in crystals from alkali systems. Observations from samples containing quenched matrix glasses suggest that the incorporation of Fe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document} is related to the abundance of tetrahedrally coordinated Al (IV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathrm {^{IV}}$$\\end{document}Al), implying some steric constraints over Fe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document} partitioning between clinopyroxene and liquid (i.e., DFe2O3cpx-liq\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$D\\mathrm {^{{cpx-liq}}_{{Fe_{2}O_{3}}}}$$\\end{document} values), though this may not be an equilibrium relationship. For example, IV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathrm {^{IV}}$$\\end{document}Al-rich {hk0}\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\{hk0\\}$$\\end{document} prism sectors of sector-zoned crystals contain more Fe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document} than IV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathrm {^{IV}}$$\\end{document}Al-poor {1¯11}\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\{\\bar{1}11\\}$$\\end{document} hourglass sectors. Moreover, IV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathrm {^{IV}}$$\\end{document}Al-rich compositions formed during disequilibrium crystallisation are enriched in Fe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document}. Apparent clinopyroxene-liquid Fe2+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{2+}$$\\end{document}–Mg exchange equilibria (i.e., KD,Fe2+-Mgcpx-liq\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$K\\mathrm{{_{D, {Fe^{2+}-Mg}}^{cpx-liq}}}$$\\end{document} values) are similarly affected by disequilibrium crystallisation in our samples. Nonetheless, it is possible to reconcile our observed clinopyroxene compositions with glass Fe valence systematics estimated from olivine-liquid equilibria if we assume that KD,Fe2+-Mgcpx-liq\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$K\\mathrm{{_{D, {Fe^{2+}-Mg}}^{cpx-liq}}}$$\\end{document} values lies closer to experimentally reported values of 0.24-\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$-$$\\end{document}0.26 than values of ∼\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim$$\\end{document}0.28 returned from a general model. In this case, olivine-liquid and clinopyroxene-liquid equilibria record equivalent narratives, with one of our glassy samples from Iceland recording evolution under fO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{\ extrm{O}_{2}}$$\\end{document} conditions about one log unit above fayalite-magnetite-quartz (FMQ) equilibrium (i.e., ∼\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim$$\\end{document}FMQ+1) and our glassy Azorean sample recording evolution under significantly more oxidising conditions (≥\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\ge$$\\end{document}FMQ+2.5) before experiencing syn-eruptive reduction, likely as a result of SO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document} degassing; our other glassy sample from Iceland was also affected by reductive SO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document} degassing. Overall, our findings demonstrate that the Fe valence systematics of clinopyroxene crystals can record information about the conditions under which OIBs evolve, but that further experimental work is required to properly disentangle the effects of magma composition, disequilibrium and fO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{\ extrm{O}_{2}}$$\\end{document} conditions on clinopyroxene-liquid equilibria involving Fe2+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{2+}$$\\end{document} and Fe3+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{3+}$$\\end{document}.

Correction to: The Origin and Differentiation of CO2-Rich Primary Melts in Ocean Island Volcanoes: Integrating 3D X-Ray Tomography with Chemical Microanalysis of Olivine-Hosted Melt Inclusions from Pico (Azores)

Correction to: The Origin and Differentiation of CO2-Rich Primary Melts in Ocean Island Volcanoes: Integrating 3D X-Ray Tomography with Chemical Microanalysis of Olivine-Hosted Melt Inclusions from Pico (Azores)

Density structure of Kīlauea volcano: implications for magma storage and transport

SUMMARY A Bayesian linear regression to determine the bias in the Nafe–Drake relationship between compressional velocity and density provides an improved model for the density structure of Kīlauea volcano, Hawaiʻi. In previous work, we combined the results of seismic tomography with the Nafe–Drake relationship between compressional velocity and density to explain the large values of gravity disturbances overlying the summits and rift zones of the island's volcanoes. These results were used to determine mechanisms for gravitational instability of the island flanks. Here, we use laboratory measurements of the relationship of velocity and density for a wide range of Hawaiʻi island rocks as a prior in a Bayesian regression, with seismic tomography, to refine the 3-D density structure for Kīlauea volcano. This refined structure shows dense bodies (3220 kg m−3) between 5 and 8 km below sea level that underly regions of magma storage, found from geodetic and geophysical studies, beneath the summit and East Rift Zone of Kīlauea volcano. Above these bodies, density isosurfaces surround and cradle sources of pressure change determined from geodetic models, both at the summit and along the East Rift Zone. Continued subsidence of the summit following the 2018 eruption is aligned with a bowl-shaped density structure, formed primarily by density isosurfaces between 2800 and 2900 kg m−3 at 4–6 km depth. These surfaces underly the ∼3 km depth at which dyke injection initiates, are largely aseismic, and from their density values are inferred to contain high concentrations of olivine. Taken together, these density structures are consistent with an olivine-rich mush with variable porosity that increases in density with depth and provides a mechanism to form olivine cumulates both at the summit and along the rift zones. This structural framework for Kīlauea volcano is consistent with melt and mush transport occurring over a large range of depths to accommodate the growth and spreading of the volcano.

Microstructural and microchemical analysis of zircon in a syenite lithic fragment from Ulleung Island volcano, South Korea

Abstract Background The intricate textural patterns commonly observed in metamorphosed and recrystallized zircon (ZrSiO4) underscore the crucial necessity of understanding the underlying mechanisms governing their formation to ensure accurate interpretation of the chemical and isotope data they contain. This study employed a combination of microanalytical techniques, including electron backscattered diffraction (EBSD) analysis, electron microprobe (EMP) mapping, and scanning electron microscope (SEM) imaging, to investigate the processes of formation and modification of zircon in a late Pleistocene (~ 35 ka) syenite enclosed within the Nari Tephra Formation on Ulleung Island in South Korea. Findings Under cathodoluminescence (CL), zircons within the syenite reveal dark, featureless, or oscillatory-zoned cores containing numerous inclusions of britholite. These cores are partially or entirely replaced by inward-penetrating bright-CL domains that exhibit minimal inclusion presence. Despite these changes, the external morphologies of the zircons remain largely unchanged, and the faded oscillatory zoning is preserved in the replaced regions. EMP mapping discloses amoebiform micro-domains with high Y, U, and Th concentrations within the dark-CL cores, while the bright-CL domains are relatively deficient in these trace elements. Microstructural analysis of the zircons using EBSD mapping indicates no significant misorientation between the dark-CL cores and the bright-CL rims. Deformation-related low-angle boundaries by lattice distortion are clearly observed in certain grains, cutting across the discrete SEM and EMP domains, and often aligned along submicron pore trails. Conclusions Microstructural and microchemical analyses carried out in this study establish that the zircons within the Ulleung syenite have undergone subsolidus recrystallization, a process likely influenced by the presence of fresh melts or fluids. This recrystallization process could be attributed to either coupled dissolution and reprecipitation or thermoactivated particle and defect volume diffusion due to inherent lattice strain. The subsequent deformation observed in the zircons might be a result of increased stress within the magma system after the recrystallization.

The Origin and Differentiation of CO2-Rich Primary Melts in Ocean Island Volcanoes: Integrating 3D X-Ray Tomography with Chemical Microanalysis of Olivine-Hosted Melt Inclusions from Pico (Azores).

ABSTRACT Constraining the initial differentiation of primary mantle melts is vital for understanding magmatic systems as a whole. Chemical compositions of olivine-hosted melt inclusions preserve unique information about the mantle sources, crystallisation behaviour and volatile budgets of such melts. Crucially, melt inclusion CO2 contents can be linked to mantle CO2 budgets and inform us on Earth's carbon fluxes and cycles. However, determining total inclusion CO2 contents is not straightforward, as they often need to be reconstructed from CO2 dissolved in melts and CO2 stored in a vapour bubble. Here, we improve upon existing reconstruction methods by combining 3D X-ray computed tomography (CT) with geochemical microanalyses of major, trace and volatile elements. We show that in comparison to CT data, traditional reconstruction methods using 2D photomicrographs can underestimate CO2 budgets by more than 40%. We applied our improved methods to basaltic olivine-hosted melt inclusions from Pico volcano (Azores) in order constrain the formation and differentiation of volatile-rich primary melts in the context of a mantle plume. Results for these inclusions yielded 1935 to 9275 μg/g reconstructed total CO2, some of the highest values reported for ocean island volcanoes to date. Using these CO2 concentrations, we calculate entrapment pressures of 105 to 754 MPa that indicate a magma reservoir comprising stacked sills straddling the crust–mantle boundary. In the magma reservoir, crystallisation of volatile saturated melts drives extensive degassing, leading to fractionated CO2/Ba ratios of 3.5 to 62.2 and a loss of over 79% of primary mantle-derived CO2. Variabilities in trace elements (La, Y) show that differentiation occurred by concurrent mixing and crystallisation of two endmember melts, respectively depleted and enriched in trace elements. Geochemical models show that enriched endmember melts constitute 33 wt % of all melts supplied to the crust at Pico and that primary melts underwent 60% crystallisation prior to eruption. Mantle melting models indicate that the enriched and depleted primary melt endmembers are low- and high-degree melts of carbon-poor lherzolite and carbon-rich pyroxenite, respectively. Moreover, since deep magmas at Pico island are dominantly pyroxenite derived, their CO2-enrichement is mainly controlled by mantle source carbon content. Overall, our study illustrates that by combining 3D imaging, geochemical microanalyses and numerical modelling, melt inclusions provide a unique record of differentiation and storage of deep magmas, as well as mantle melting.

Long-lived shallow crustal magma reservoirs beneath Ulleung Island volcano (South Korea) and their rejuvenation recorded in zircon

The microanalytical data obtained from zircons collected on Ulleung Island, an emergent peak of a Quaternary alkaline magma system in the southwestern East Sea (Sea of Japan) back-arc basin, provide new insights into the origin of magma, magmatic processes, and the trigger of explosive eruptions. The U-Th ages of zircons extracted from a pumice layer formed during the intense Plinian eruption at ∼11 ka, known as the Ulleung-Oki tephra, show a wide range spanning the pre-eruption equilibrium state (>350 ka) up to the eruption event itself. Zircon grains derived from syenite and gabbro fragments entrained within the tephra yield weighted mean U-Th and U-Pb ages of 34.7 ± 2.3 ka and 1378 ± 22 ka, respectively. Under cathodoluminescence, zircons from the syenite reveal dark, featureless, or oscillatory-zoned cores, containing numerous inclusions of britholite-like Y-rich phases and pyrochlore. These cores are partially or entirely replaced by inclusion-free, brighter domains, suggesting that in situ dissolution and precipitation processes significantly influenced the transformation of zircon. The consistently low δ18O values (5.4–2.6‰) of zircons from Middle Pleistocene trachyte lavas, the Ulleung-Oki tephra, and lithic fragments indicate the persistent presence of shallow crustal magma reservoirs beneath Ulleung Island. These reservoirs have undergone substantial assimilation of hydrothermally altered rocks, with assimilation rates typically around 10–30% and reaching a maximum of approximately 50%. Zircon core-to-rim variations in oxygen isotopic composition, along with textural observations, suggest the rejuvenation of magma reservoirs through the injection of new melt, which likely served as a trigger for the explosive eruptions at Ulleung Island. The similarity in the range of εHf values between the core and rim (+7.4 to −4.3) implies the Hf isotopic consistency of the recharge melt.

INVESTIGATIONS ON INSAR DATA PROCESSING STANDARD FOR VOLCANO ISLAND MONITORING IN INDONESIA

Abstract. In critical situations such as disaster and humanitarian relief, the reliable geospatial data will present an appropriate decision on a timely manner. As fundamental datasets, large-scale topographic maps are are mandatory in order to perform geospatial analysis for many societal challenges. Indeed, Large Scale Topographic Mapping in Indonesia shall be accelerated in an innovative and productive way to consider the production efficiency, especially for the geospatial data acquisition as its primary source. Interferometric Synthetic Aperture Radar (InSAR) data acquisition presents some advantages, especially in conjunction with the security clearance and global weather as some main constraints for the geospatial data acquisition. In terms of disaster preparedness and emergency response, radar interferometry techniques can play an important role by generating a Digital Elevation Model (DEM) as an input to Decision Support Systems (DSS). In addition, the techniques of differential Interferometric Synthetic Aperture Radar (D-InSAR) can provide the earth surface deformation from time series Radar Datasets. A set of standards related to the InSAR data acquisition are essential, especially with regards to the improved sensors, and processing methods implemented by the worldwide space agencies nowadays. Hence, the precise definition of a standard is fundamental in order to prepare the InSAR data as one potential source of geospatial data and information, especially for LSTM data production. Generated time series of Digital Elevation Model (DEM) can identify the potential cause of 2018’s Sunda Strait tsunami (Tampubolon, 2019). It was related with the increasing eruption activity of Anak Krakatau (Child of Krakatoa). From DEM analysis, the Volcanic mudflow material avalanche has been detected before the tsunami hit the coastal surrounding area. In this paper, more advanced techniques was investigated by the combination between Time Series DEM Analysis and Ground Displacement detection. The focus on the integration between reliable TanDEM-X and up-to-date Sentinel 1A/B was proposed on another volcano island of Gamalama in Ternate. As TanDEM-X uses X-band as its medium wave, it can produce more reliable DEM as a reference data. On the other hand, Sentinel 1 uses C-band that can contribute to detect the ground deformation over the aforementioned volcano island. TanDEM-X uses a more robust and consistent approach to scan the earth surface by single data acquisition scheme, which is generally more accurate and reliable for the DEM generation. However, many factors influence the selection of a platform that can be best suited for radar interferometry. Initially, the German TanDEM-X Coregistered Singlelook Slant-range Complex (CoSSC) using bi-static InSAR approach has been implemented in order to generate DEMs comply with the LSTM specification in Indonesia. The basic principle is performing a simultaneous measurement of the same scene and identical doppler spectrum by using 2 sensors, thereby avoiding temporal decorrelation (DLR, 2012). Subsequently, the differential InSAR techniques by using Sentinel 1A/B has been performed in order to detect the ground deformation. Finally, tt can be shown that the reliable TanDEM-X data has increased the ground deformation resolution and accuracy in order to provide potential solution for disaster preparedness.

Ash Deposition Triggers Phytoplankton Blooms at Nishinoshima Volcano, Japan

AbstractVolcanoes that deposit eruptive products into the ocean can trigger phytoplankton blooms near the deposition area. Phytoplankton blooms impact the global carbon cycle, but the specific conditions and mechanisms that facilitate volcanically triggered blooms are not well understood, especially in low nutrient ocean regions. We use satellite remote sensing to analyze the chlorophyll response to an 8‐month period of explosive and effusive activity from Nishinoshima volcano, Japan. Nishinoshima is an ocean island volcano in a low nutrient low chlorophyll region of the Northern Pacific Ocean. From June to August 2020, during explosive activity, satellite‐derived chlorophyll‐a was detectable with amplitudes significantly above the long‐term climatological value. After the explosive activity ceased in mid‐August 2020, these areas of heightened chlorophyll concentration decreased as well. In addition, we used aerial observations and satellite imagery to demonstrate a spatial correlation between blooms and ash plume direction. Using a sun‐induced chlorophyll‐a fluorescence satellite product, we confirmed that the observed chlorophyll blooms are phytoplankton blooms. Based on an understanding of the nutrients needed to supply blooms, we hypothesize that blooms of nitrogen‐fixing phytoplankton led to a 1010–1012 g drawdown of carbon. Thus, the bloom could have significantly mediated the output of carbon from the explosive phase of the eruption but is a small fraction of anthropogenic CO2 stored in the ocean or the global biological pump. Overall, we provide a case study of fertilization of a nutrient‐poor ocean with volcanic ash and demonstrate a scenario where multi‐month scale deposition triggers continuous phytoplankton blooms across 1,000s of km2.

Zircon double-dating, trace element and O isotope analysis to decipher late Pleistocene explosive-effusive eruptions from a zoned ocean-island magma system, Ascension Island

In this first detailed study of zircon from Ascension Island, South Atlantic, we take a novel approach combining trace element and O isotope compositional data with double-dating (disequilibrium 238U–230Th and (U–Th)/He) to decipher timescales and dynamics of magmatic processes. The Echo Canyon (EC) sequence comprises small-volume explosive-effusive eruptions of trachyte that tapped a compositionally zoned magma system. Associated volcanic hazards may be constrained from the age of volcanism, duration of magma storage, and magma source and plumbing system character. Zircon U–Th–Pb dating of lithic lava clasts has revealed recurrent evolved volcanism at 1.34 and 0.6 Ma, and 95 ka. The (U–Th)/He zircon cooling ages indicate that most of the EC explosive-effusive sequence erupted in a brief episode at ca. 95 ka. Additionally, uniform 238U–230Th zircon crystallisation ages suggest moderately protracted magma storage with melt present at depth for at most 103–104 years before eruption. The enriched character of zircon trace element compositions, relative to MORB, in the absence of a continental crustal signature in the oxygen isotope values (δ18O range 2.67–5.63‰), suggests the presence of an enriched component in the EC magma source. Furthermore, low δ18O zircon compositions imply assimilation of high temperature hydrothermally altered country rock by the source magma. The mineral assemblage in crystal-poor pumices indicates equilibrium storage conditions: zircon saturation and Ti-in-zircon crystallisation temperatures are consistent with alkali feldspar-melt temperatures. Significantly, zircon crystals were preserved both as macrocryst inclusions and in the groundmass of EC explosive and effusive deposits. These rocks preserve evidence of magma evolution by fractional crystallisation. This process led to pre-eruptive compositional stratification, which is evidenced in the range of whole-rock major and trace element compositions and zircon Zr/Hf values. Notably, zircon crystallisation and cooling ages derived from pumice, lava, and accidental lithic lava clasts in highly explosive pyroclastic deposits, have revealed episodes of evolved magmatism that would otherwise have gone undetected. In addition, the zircon trace element and isotope compositions, in combination with the range of crystallisation ages, evidence progressively deeper tapping of less evolved magma stored in discrete lenses. Thus, a combined zircon geochronological-geochemical approach can place constraints on the ca. 0.6 Ma recurrence of past explosive-effusive pulses of millennial to decamillennial duration and their enriched magma sources. This information is relevant for assessing hazards and informing monitoring and forecasting efforts to assist in managing associated risks for small ocean island volcanoes with particularly vulnerable populations and infrastructure.

Progressive underplating of mafic material at mid-crustal depth beneath Ischia volcano, Italy

The destructive (Mw 3.9) earthquake of 21 August 2017 re-opened the question on where magma resides at the Ischia island volcano. The peculiar complexity of the seismic source initiated the debate on the involvement of fluid-related processes, but the magmatic origin of the event remains uncertain. Here we use ground displacement and broadband seismic data to investigate the magmatic system of Ischia volcano, where progressive underplating of mafic material at mid-crustal depth feeds the deep structures, which are characterized by exceptionally high seismic velocity and are connected with a shallow crystal mush. Although no direct evidence of large molten volumes was found, strong anisotropy suggests that the crystal mush is pervaded by magma-intruded dykes. We propose that the 2017 event was due to a negative tensile deformation caused by depressurization of supercritical fluids along a shallow southwest (SW-)-dipping fault defined by receiver functions (RFs) data, which acted as a valve regulating the overpressure of deep magmatic fluids.

Amplitude and inter-event time statistics for the island volcanoes Stromboli, Mount Etna, Yasur, and Whakaari

Detailed analyses of past major and minor seismo-volcanic events can help to understand the eruptive behavior of volcanoes and the underlying physical and chemical processes. Catalogs of these eruptions and, specifically, seismo-volcanic events may be generated using continuous seismic recordings at stations in the proximity of volcanoes. Here, we apply a recently-developed automated approach Adaptive-Window Volcanic Event Selection Analysis Module (AWESAM) to seismic data from Stromboli (Italy), Mount Etna (Italy), Yasur (Vanuatu) and Whakaari (New Zealand). We perform an inter-event time analysis to identify characteristic patterns in the events’ recurrence time and the volcanic activity. Using this identical approach for all volcanoes, we were able to discover that despite their differing types and activity, they exhibit similar statistical behaviors. For Whakaari, we noticed a bimodal inter-event time distribution for large events. Since this observation is based on single station data, further in-depth investigations are needed once more data is available in future. We also derive a new amplitude-frequency relationship from seismo-volcanic events. With this relation, we can confirm a change in slope for large events at Stromboli, which is based on 10 years of data. Additionally, we apply a classification model to events from Stromboli to differentiate between low-period (LP) events and high-frequency (HF) events and found an alternating behavior in the frequency of these events before and after the two paroxysms in 2019.

Two Distinct Magma Storage Regions at Ambrym Volcano Detected by Satellite Geodesy

AbstractThe flux of eruptible magma into a magmatic plumbing system influences eruption size and timing. If magma transfer is possible between two hydraulically‐connected magma lenses, system destabilization can tap a larger magma volume than stored in any one melt lens. This study identifies two distinct magma reservoirs beneath Ambrym, a basaltic island volcano in Vanuatu, during the time period February 2019 to January 2022. Using InSAR time series and a data assimilation approach, we estimate pressure changes within two reservoirs (located 5–7 and 4–6 km b.s.l.). Furthermore, a theoretical model demonstrates that the reservoirs may not currently be hydraulically connected, despite evidence of physical mixing of magma derived from each reservoir during the December 2018 eruption. These findings further our understanding of how magmatic plumbing systems at basaltic calderas may change after rift‐zone eruptions.

Multi-station volcano tectonic earthquake monitoring based on transfer learning

Introduction: Developing reliable seismic catalogs for volcanoes is essential for investigating underlying volcanic structures. However, owing to the complexity and heterogeneity of volcanic environments, seismic signals are strongly affected by seismic attenuation, which modifies the seismic waveforms and their spectral content observed at different seismic stations. As a consequence, the ability to properly discriminate incoming information is compromised. To address this issue, multi-station operational frameworks that allow unequivocal real-time management of large volumes of volcano seismic data are needed.Methods: In this study, we developed a multi-station volcano tectonic earthquake monitoring approach based on transfer learning techniques. We applied two machine learning systems—a recurrent neural network based on long short-term memory cells (RNN–LSTM) and a temporal convolutional network (TCN)—both trained with a master dataset and catalogue belonging to Deception Island volcano (Antarctica), as blind-recognizers to a new volcanic environment (Mount Bezymianny, Kamchatka; 6 months of data collected from June to December 2017, including periods of quiescence and eruption).Results and discussion: When the systems were re-trained under a multi correlation-based approach (i.e., only seismic traces detected at the same time at different seismic stations were selected), the performances of the systems improved substantially. We found that the RNN-based system offered the most reliable recognition by excluding low confidence detections for seismic traces (i.e., those that were only partially similar to those of the baseline). In contrast, the TCN-based network was capable of detecting a greater number of events; however, many of those events were only partially similar to the master events of the baseline. Together, these two approaches offer complementary tools for volcano monitoring. Moreover, we found that our approach had a number of advantages over the classical short time average over long time-average (STA/LTA) algorithm. In particular, the systems automatically detect VTs in a seismic trace without searching for optimal parameter settings, which makes it a portable, scalable, and economical tool with relatively low computational cost. Moreover, besides obtaining a preliminary seismic catalog, it offers information on the confidence of the detected events. Finally, our approach provides a useful tentative label for subsequent analysis carried out by a human operator. Ultimately, this study contributes a new framework for rapid and easy volcano monitoring based on temporal changes in monitored seismic signals.

Developing a South Pacific tephra framework: Initial results from a Samoan Holocene sequence

ABSTRACTTephra preserved in sediments form useful isochronous marker layers, linking disparate geological, palaeoenvironmental and archaeological records. The application of tephrochronology is greatly enhanced through the detection of macroscopically invisible tephra (cryptotephra). Here, we identify two discrete cryptotephra in Samoan lake sediments, the first identification of cryptotephra in the region outside of New Zealand. Geochemical data suggest one ash layer is from a local Samoan source, providing the first data on an eruption of this age, adding to knowledge of the local volcanic record. The second has a distinctive rhyolitic glass composition, which matches either that of Raoul Island in the Kermadec Arc (1800 km south of Samoa), or two currently submarine volcanoes in the Tongan Arc, ‘Volcano F’ and Lateiki/Metis Shoal (550 and 700 km south of Samoa, respectively). In all possible source cases, this points to a regionally significant eruption of a Kermadec–Tongan volcano at ca. 10 000 a bp. The study marks the first step in the establishment of a South Pacific tephra framework that can be used to answer questions about the synchronicity of changes in hydroclimate, vegetation and early Polynesian migration patterns, as well as providing more information on the volcanic history of Pacific island volcanoes.

The upper Pleistocene (1.8–0.7 Ma) explosive eruptive history of Las Cañadas, ocean-island volcano, Tenerife

While most ocean-island volcanism is effusive, recent evidence has demonstrated that intraplate ocean island volcanoes can exhibit protracted explosive histories, with catastrophic eruption styles and hazardous behaviour more typically associated with volcanoes in continental and plate-margin settings. Tenerife is the largest explosive ocean-island volcano on Earth, with a prolonged (∼2 Ma) post-erosional history of caldera-forming, plinian and ignimbrite eruptions of evolved composition. The 0.7–1.8 Ma succession with 20 newly defined formations is reported for southern Tenerife, Canary Islands. In the last 2 Myr, the Las Cañadas volcano has produced >42 pumice-fall eruptions, 21 with extensive ignimbrites, and 12 inferred caldera collapse events. Pyroclastic density currents have repeatedly travelled more than15 km from source to the ocean, filling valleys and burying extensive interfluves. A robust whole-rock chemistry dataset, selected mineral chemistry, coupled with new 40Ar/39Ar ages of units through the pyroclastic stratigraphy, allow recognition of magmatic trends within the system on the order of 100 ky. The catastrophic explosive eruptions form three, 0.2–0.5 Myr-duration clusters (the Ucanca, Guajara and Diego-Hernandez) that do not appear to correspond simply with geochemical cycles, or to cycles of increasing eruption size or explosivity as has been previously proposed. During the clusters, large eruption frequencies averaged 1 every 20–45 kyrs. The eruption clusters were separated by hiatuses of ∼240–260 kyr, recorded by soils and unconformities, and may reflect marked changes in geographic dispersals following giant landslide breaches in Las Cañadas caldera wall. Two concurrent evolutionary magmatic trends are distinguished: one producing crystal-rich magmas, the other formed the cooler crystal-poor magmas: both spanned over a million years until 0.66 Ma, when the former ceased.

Identification of Seismo‐Volcanic Regimes at Whakaari/White Island (New Zealand) Via Systematic Tuning of an Unsupervised Classifier

AbstractWe present an algorithm based on Self‐Organizing Maps (SOM) and k‐means clustering to recognize patterns in a continuous 12.5‐year tremor time series recorded at Whakaari/White Island volcano, New Zealand (hereafter referred to as Whakaari). The approach is extendable to a variety of volcanic settings through systematic tuning of the classifier. Hyperparameters are evaluated by statistical means, yielding a combination of “ideal” SOM parameters for the given data set. Extending from this, we applied a Kernel Density Estimation approach to automatically detect changes within the observed seismicity. We categorize the Whakaari seismic time series into regimes representing distinct volcano‐seismic states during recent unrest episodes at Whakaari (2012/2013, 2016, and 2019). There is a clear separation in classification results between background regimes and those representing elevated levels of unrest. Onset of unrest is detected by the classifier 6 weeks before the August 2012 eruption, and ca. 3.5 months before the December 2019 eruption, respectively. Regime changes are corroborated by changes in commonly monitored tremor proxies as well as with reported volcanic activity. The regimes are hypothesized to represent diverse mechanisms including: system pressurization and depressurization, degassing, and elevated surface activity. Labeling these regimes improves visualization of the 2012/2013 and 2019 unrest and eruptive episodes. The pre‐eruptive 2016 unrest showed a contrasting shape and nature of seismic regimes, suggesting differing onset and driving processes. The 2016 episode is proposed to result from rapid destabilization of the shallow hydrothermal system, while rising magmatic gases from new injections of magma better explain the 2012/2013 and 2019 episodes.