The 2.6-2.3 ka explosive eruptive period of the Pululahua dome complex, Ecuador: insights from pyroclast analysis.

Pululahua is an active volcano located 15 km north of Quito, Ecuador, that comprises sixteen dacitic-andesitic lava domes and a 13 km2 sub-rectangular depression formed between ~ 2.6 and ~ 2.3 ka. We use a detailed study of 70 flow and fall deposits that make up the pyroclastic sequence to show that the depression, previously classified as a caldera, was formed by numerous Vulcanian to (sub-) Plinian eruptions that destroyed both earlier and co-eruptive lava domes. We support this interpretation with field work, analysis of grain size distributions, density and components of 24 key deposits, supplemented by textural and petrologic analyses of 16 juvenile pyroclasts from throughout the pyroclastic sequence. These data document an alternation of (sub-) Plinian and Vulcanian eruptions dominated by denser juvenile material that preserves microtextural variations indicating changes in shallow level magma storage accompanying Vulcanian explosions. An exploratory examination of phenocryst textures and plagioclase and amphibole rim compositions suggests that much of the eruptive activity was driven by repeated inputs of less evolved magma into the Pululahua magmatic system. The inferred sequence of events provides a new hypothesis for the formation of the current morphology of Pululahua, including multiple episodes of both effusive and explosive eruptions accompanied by vent migration. Our findings offer an important insight into Pululahua's potential future hazard scenarios, which could affect millions of people. Supplementary Information: The online version contains supplementary material available at 10.1007/s00445-022-01590-4.

El Pululahua es un volcán activo situado 15 km al norte de Quito, Ecuador, que comprende dieciséis domos de lava dacítica-andesítica y una depresión sub-rectangular de 13 km² formada entre ~2,6 y ~2,3 mil años atrás. En este artículo realizamos un estudio detallado de 70 depósitos de flujos y caídas que componen la secuencia piroclástica para mostrar que la depresión, previamente clasificada como caldera, se formó en realidad por numerosas erupciones Vulcanianas a (sub-) Plinianas que destruyeron domos de lava pre-existentes y co-eruptivos. Fundamentamos esta interpretación con trabajo de campo, el análisis de densidad, componentes y distribución del tamaño de grano de 24 depósitos clave, además del análisis textural y petrológico de 16 piroclastos juveniles muestreados a lo largo de toda la secuencia. Los datos documentan una alternancia entre erupciones (sub-) Plinianas y Vulcanianas. La predominancia de material juvenil denso que conserva variaciones microtexturales dentro de los depósitos Vulcanianos indica cambios en el almacenamiento de magma a poca profundidad antes de estas erupciones. Un análisis exploratorio de las texturas de los fenocristales y de las composiciones de los bordes de plagioclasa y anfíbol sugiere que gran parte de la actividad eruptiva fue desencadenada por repetidas inyecciones de magma menos diferenciado al sistema magmático del Pululahua. La secuencia inferida de eventos eruptivos entre ~2,6 y ~2,3 mil años proporciona una nueva hipótesis para la formación de la actual morfología del Pululahua, incluyendo múltiples episodios de erupciones efusivas y explosivas acompañadas de la migración de los ventos eruptivos. Nuestros resultados aportan una nueva interpretación de los posibles futuros escenarios de peligro asociados al volcán Pululahua, los cuales podrían afectar a millones de personas.

Community established best practice recommendations for tephra studies-from collection through analysis.

Tephra is a unique volcanic product with an unparalleled role in understanding past eruptions, long-term behavior of volcanoes, and the effects of volcanism on climate and the environment. Tephra deposits also provide spatially widespread, high-resolution time-stratigraphic markers across a range of sedimentary settings and thus are used in numerous disciplines (e.g., volcanology, climate science, archaeology). Nonetheless, the study of tephra deposits is challenged by a lack of standardization that inhibits data integration across geographic regions and disciplines. We present comprehensive recommendations for tephra data gathering and reporting that were developed by the tephra science community to guide future investigators and to ensure that sufficient data are gathered for interoperability. Recommendations include standardized field and laboratory data collection, reporting and correlation guidance. These are organized as tabulated lists of key metadata with their definition and purpose. They are system independent and usable for template, tool, and database development. This standardized framework promotes consistent documentation and archiving, fosters interdisciplinary communication, and improves effectiveness of data sharing among diverse communities of researchers.

Prince Rupert's Drops: An analysis of fragmentation by thermal stresses and quench granulation of glass and bubbly glass.

When volcanic eruptions involve interaction with external water (hydrovolcanism), the result is an ash-rich and energetic volcanic plume, as illustrated dramatically by the January 2022 Tonga eruption. The origin of the high explosive energy of these events remains an important question. We investigate this question by studying Prince Rupert's Drops (PRDs)-tadpole-shaped glass beads formed by dripping molten glass into water-which have long fascinated materials scientists because the great strength of the head contrasts with the explosivity of the metastable interior when the tail is broken. We show that the fragment size distribution (FSD) produced by explosive fragmentation changes systematically with PRD fragmentation in air, water, and syrup. Most FSDs are fractal over much of the size range, scaling that can be explained by the repeated fracture bifurcation observed in three-dimensional images from microcomputed tomography. The shapes of constituent fragments are determined by their position within the original PRD, with platey fragments formed from the outer (compressive) shell and blocky fragments formed by fractures perpendicular to interior voids. When molten drops fail to form PRDs, the glass disintegrates by quench granulation, a process that produces fractal FSDs but with a larger median size than explosively generated fragments. Critically, adding bubbles to the molten glass prevents PRD formation and promotes quench granulation, suggesting that granulation is modulated by heterogeneous stress fields formed around the bubbles during sudden cooling and contraction. Together, these observations provide insight into glass fragmentation and potentially, processes operating during hydrovolcanism.

Data on the 1902 Plinian eruption of Santa María volcano, Guatemala.

The study of historic volcanic eruptions is often complicated by the lack of recorded primary data and observations of such events. In the case of large-magnitude historic eruptions, these types of data are important to better understand not only the physical nature of these rare events but also the volcanic and social impacts that follow. In this paper, we compile contemporary data on the Santa María Plinian eruption of 1902, in Guatemala. The data supplement those presented in the original research article [1] but individually provide an interesting and useful compilation of eyewitness testimonies, scientific studies and newspaper reports. We identify key contemporary sources containing quantitative data as well as various qualitative reports that we convert to quantitative measurements through a simple classification scheme. We also compile wind reanalysis data from the time of the eruption to display wind direction and speed with height. Both the data and the description of the methods of data analysis can aid future studies of qualitative (historic, eyewitness) to quantitative data conversion, as well as studies investigating this important eruption.

Deep roots for mid-ocean-ridge volcanoes revealed by plagioclase-hosted melt inclusions.

The global mid-ocean ridge system is the most extensive magmatic system on our planet and is the site of 75 per cent of Earth's volcanism1. The vertical extent of mid-ocean-ridge magmatic systems has been considered to be restricted: even at the ultraslow-spreading Gakkel mid-ocean ridge under the Arctic Ocean, where the lithosphere is thickest, crystallization depths of magmas that feed eruptions are thought to be less than nine kilometres2. These depths were determined using the volatile-element contents of melt inclusions, which are small volumes of magma that become trapped within crystallizing minerals. In studies of basaltic magmatic systems, olivine is the mineral of choice for this approach2-6. However, pressures derived from olivine-hosted melt inclusions are at odds with pressures derived from basalt major-element barometers7 and geophysical measurements of lithospheric thickness8. Here we present a comparative study of olivine- and plagioclase-hosted melt inclusions from the Gakkel mid-ocean ridge. We show that the volatile contents of plagioclase-hosted melt inclusions correspond to much higher crystallization pressures (with a mean value of 270 megapascals) than olivine-hosted melt inclusions (with a mean value of 145 megapascals). The highest recorded pressure that we find equates to a depth 16.4 kilometres below the seafloor. Such higher depths are consistent with both the thickness of the Gakkel mid-ocean ridge lithosphere and with pressures reconstructed from glass compositions. In contrast to previous studies using olivine-hosted melt inclusions, our results demonstrate that mid-ocean-ridge volcanoes may have magmatic roots deep in the lithospheric mantle, at least at ultraslow-spreading ridges.

The significance of plagioclase textures in mid-ocean ridge basalt (Gakkel Ridge, Arctic Ocean).

Textures and compositions of minerals can be used to infer the physiochemical conditions present within magmatic systems. Given that plagioclase is an abundant phase in many magmatic systems, understanding the link between texture and process is vital. Here, we present a database of textural and compositional data for > 1800 plagioclase crystals in mid-ocean ridge basalt from the Gakkel Ridge (Arctic Ocean) to investigate the physiochemical conditions and processes that govern the formation of plagioclase textures and compositions. The Gakkel basalts have high modal crystal contents (up to 50%). The crystal cargo is complex, with both individual plagioclase and glomerocrysts showing large variations in crystal habit, zoning and resorption. The most common types of zoning are reverse and patchy; we attribute patchy zoning to infilling following either skeletal growth or resorption. Resorption is abundant, with multiple resorption events commonly present in a single crystal, and results from both magmatic recharge and decompression. Periods of strong undercooling, distinct to quench crystallisation, are indicated by matured skeletal crystals and thin normally zoned melt inclusion-rich bands following resorption. Individual samples often contain diverse textural and compositional plagioclase groups. Furthermore, most plagioclase is not in equilibrium with its host melt. Finally, the porous open structures of some glomerocrysts suggest that they represent pieces of entrained disaggregated mush. We interpret this to indicate that the crystal cargo is not generally phenocrystic in origin. Instead, plagioclase crystals that formed in different parts of a mush-dominated plumbing system were entrained into ascending melts. The textures of individual crystals are a function of their respective histories of (under)cooling, magma mixing and decompression. The morphologies of melt inclusion trapped in the plagioclase crystals are associated with specific host crystal textures, suggesting a link between plagioclase crystallisation processes and melt inclusion entrapment. The database of plagioclase presented herein may serve as a template for the interpretation of plagioclase textures in magmatic systems elsewhere.

Architecture and dynamics of magma reservoirs.

This introductory article provides a synopsis of our current understanding of the form and dynamics of magma reservoirs in the crust. This knowledge is based on a range of experimental, observational and theoretical approaches, some of which are multidisclipinary and pioneering. We introduce and provide a contextual background for the papers in this issue, which cover a wide range of topics, encompassing magma storage, transport, behaviour and rheology, as well as the timescales on which magma reservoirs operate. We summarize the key findings that emerged from the meeting and the challenges that remain. The study of magma reservoirs has wide implications not only for understanding geothermal and magmatic systems, but also for natural oil and gas reservoirs and for ore deposit formation. This article is part of the Theo Murphy meeting issue 'Magma reservoir architecture and dynamics'.

Mafic glass compositions: a record of magma storage conditions, mixing and ascent.

The trans-crustal magma system paradigm is forcing us to re-think processes responsible for magma evolution and eruption. A key concept in petrology is the liquid line of descent (LLD), which relates a series of liquids derived from a single parent, and therefore tracks the inverse of the crystallization path. It is common practice to attribute multiple magma compositions, and/or multiple melt compositions (from melt inclusions and matrix glass), to a single LLD. However, growing evidence for rapid, and often syn-eruptive, assembly of multiple magma components (crystals and melts) from different parts of a magmatic system suggests that erupted magma and melt compositions will not necessarily represent a single LLD, but instead may reflect the multiple paths in pressure-temperature space. Here, we use examples from mafic magmatic systems in both ocean island and arc settings to illustrate the range of melt compositions present in erupted samples, and to explore how they are generated, and how they interact. We highlight processes that may be deduced from mafic melt compositions, including the mixing of heterogeneous primitive liquids from the mantle, pre-eruptive magma storage at a range of crustal and sub-Moho depths, and syn-eruptive mixing of melts generated from these storage regions. The relative dominance of these signatures in the glasses depends largely on the water content of the melts. We conclude that preserved melt compositions provide information that is complementary to that recorded by the volatile contents of crystal-hosted melt inclusions and coexisting mineral compositions, which together can be used to address questions about both the pre- and syn-eruptive state of volcanic systems. This article is part of the Theo Murphy meeting issue 'Magma reservoir architecture and dynamics'.

Rapid crystallization during recycling of basaltic andesite tephra: timescales determined by reheating experiments.

Microcrystalline inclusions within microlite-poor matrix are surprisingly common in low intensity eruptions around the world, yet their origin is poorly understood. Inclusions are commonly interpreted as evidence of crystallization along conduit margins. Alternatively, these clasts may be recycled from low level eruptions where they recrystallize by heating within the vent. We conducted a series of experiments heating basaltic andesite lapilli from temperatures below the glass transition (~690 °C) to above inferred eruption temperatures (>1150 °C) for durations of 2 to >60 minutes. At 690 °C < T < 800 °C, crystallization is evident after heating for ~20 minutes; at T > 800 °C, crystallization occurs in <5 minutes. At T ≥ 900 °C, all samples recrystallize extensively in 2-10 minutes, with pyroxenes, Fe-oxides, and plagioclase. Experimental crystallization textures closely resemble those observed in natural microcrystalline inclusions. Comparison of inclusion textures in lapilli from the active submarine volcano NW Rota-1, Mariana arc and subaerial volcano Stromboli suggest that characteristic signatures of clast recycling are different in the two environments. Specifically, chlorine assimilation provides key evidence of recycling in submarine samples, while bands of oxides bordering microcrystalline inclusions are unique to subaerial environments. Correct identification of recycling at basaltic vents will improve (lower) estimates of mass eruption rate and help to refine interpretations of eruption dynamics.

Vertically extensive and unstable magmatic systems: A unified view of igneous processes.

Volcanoes are an expression of their underlying magmatic systems. Over the past three decades, the classical focus on upper crustal magma chambers has expanded to consider magmatic processes throughout the crust. A transcrustal perspective must balance slow (plate tectonic) rates of melt generation and segregation in the lower crust with new evidence for rapid melt accumulation in the upper crust before many volcanic eruptions. Reconciling these observations is engendering active debate about the physical state, spatial distribution, and longevity of melt in the crust. Here we review evidence for transcrustal magmatic systems and highlight physical processes that might affect the growth and stability of melt-rich layers, focusing particularly on conditions that cause them to destabilize, ascend, and accumulate in voluminous but ephemeral shallow magma chambers.