Elsevier

Chemical Geology

Volume 532, 20 January 2020, 119301
Chemical Geology

Modelling and physico-chemical constraints to the 4.5 ka Agnano-Monte Spina Plinian eruption (Campi Flegrei, Italy)

https://doi.org/10.1016/j.chemgeo.2019.119301Get rights and content

Abstract

The 4.5 ka trachytic Plinian eruption of Agnano-Monte Spina is the largest magnitude event of the past 5 ka at Campi Flegrei caldera. The complete eruptive sequence consists of six members, three of which, named A, B and D, are characterized by the association of fallout and pyroclastic density current (PDC) deposits well preserved at proximal locations.

In this study, we analyze the textural characteristics of the pumice clasts of the three major fallout deposits (A1, B1, D1) and of their associated PDC deposits (A2, B2, D2), and link them to the physical properties of magma in order to investigate conduit fluid dynamics. A combination of data (field work, grain-size and density measurements, vesicle number densities and size distributions, crystal content, water content) is used to set up the source term conditions for numerical simulations.

Each fall/PDC transition is accompanied by distinctive changes in textural properties of the juveniles, recognized by a lowering in vesicle number densities of about one order of magnitude (from 108 to 107 cm−3), indicating a significant decrease in the magma ascent rate. Melt inclusions show a marked decrease in volatile content recurrent at each fall/PDC transition and indicate that the three main pulses of the eruption were fed by distinct and progressively deeper magma batches. Numerical simulations, taking into account magma properties derived from the textural analyses, and variations in initial water content, show decreases in the exit velocities and Mass Discharge Rate (MDR) consistent with such fall/PDC transitions. Different initial water contents together with changes in conduit diameters allow us to simulate the different column heights reconstructed for the three Plinian phases. The reconstructed scenario for the Agnano-Monte Spina eruption involves a stop-start behavior and a top-down trigger for the most voluminous and intense eruptive episode D.

Introduction

Plinian eruptions commonly show transitions from phases of sustained convective columns to phases characterized by partial or total collapse generating pyroclastic density currents (PDCs) (e.g., Orsi et al., 1992; Wohletz et al., 1995; Rosi et al., 2001; Cioni et al., 2015; Pensa et al., 2015). Furthermore, while many Plinian sequences show the “classic” transition from a lower (basal) fall deposit into Pyroclastic Density Current (PDC) deposits (Cas and Wright, 1987), many others have been described as related to more uncommon stop-start behavior (e.g. Christiansen and Peterson, 1981; Giordano, 1998; Adams et al., 2001; Yasui and Koyaguchi, 2004; Miyaji et al., 2011).

Factors that determine sustained column to PDC transitions include variations in mass discharge rate as a consequence of changing volatile contents, magma composition and/or crystallinity, which in turn affect the rheology of magma, changes in conduit geometry and magma-water interaction (e.g. Sparks, 1978; Wilson et al., 1980; Cashman and Mangan, 1994; Klug and Cashman, 1994; Lyakhovsky et al., 1996; Papale, 1999; Rust et al., 2004; Gurioli et al., 2005; Houghton et al., 2010; Shea et al., 2011, 2012; Vinkler et al., 2012; Pardo et al., 2014).

The link between eruption column and PDCs (i.e., if the PDC is related to a total or partial column collapse) is of great interest for modeling conduit dynamics aimed at volcanic hazard assessment and risk mitigation. Considerable advancements in such studies have been achieved in recent years (e.g. Esposti Ongaro et al., 2006; Dellino et al., 2007; de’ Michieli Vitturi et al., 2010d; Suzuki and Koyaguchi, 2012; Carazzo et al., 2015; Trolese et al., 2019). Despite that, the critical conditions that drive the collapse of a Plinian column still need to be thoroughly explored. Textural analysis of pyroclasts has been proved to provide reliable proxies for magma ascent dynamics (e.g. Toramaru, 2006; Gurioli et al., 2008; Piochi et al., 2008; Costantini et al., 2010; Shea et al., 2010), and, coupled with experiments or parameterization of magma rheology (e.g., Costa, 2005; Llewellin and Manga, 2005; Mader et al., 2013; Vona et al., 2016), may help in understanding and constraining the transitions in eruptive styles (Adams et al., 2006; Sable et al., 2006; Moitra et al., 2013; La Spina et al., 2016). The 4.5 ka trachytic Plinian eruption of Agnano-Monte Spina (A-MS; de Vita et al., 1999; Blockley et al., 2008; Smith et al., 2011) provides an excellent opportunity to study an uncommon explosive sequence related to at least three different eruptive pulses, each initiated by a buoyant phase later evolving into a collapsing phase. This eruption has been taken as reference eruption scenario for the next high scale eruption at Campi Flegrei (e.g. Orsi et al., 2004, 2009; Bevilacqua et al., 2015, 2017) by the Italian Civil Protection Agency and therefore it represents a very important event for the hazard studies of one of the most populated active volcanic areas on Earth.

In this study, we apply a multidisciplinary approach, using geological data (chemistry and texture of the emitted products) to set up the source term conditions for the modelling of the three major fallouts (A1, B1, D1) and of the related PDCs (A2, B2, D2) (de Vita et al., 1999; Dellino et al., 2001, 2004a). The aim of the study is to relate the textural and physical features of the deposit to intrinsic parameters of the magma during ascent along the conduit and to determine to what extent, if any, changes in magma properties may have affected the course of the A-MS eruption.

Section snippets

Campi Flegrei volcanic system and Agnano-Monte Spina eruption

The Campi Flegrei caldera is a nested structure which results from two main collapses that followed the Campanian Ignimbrite (39 ka; Fisher et al., 1993; Orsi et al., 1996; Rosi et al., 1996; Civetta et al., 1997; Fedele et al., 2008) and the Neapolitan Yellow Tuff eruptions (14 ka; Orsi et al., 1992, 1995; Wohletz et al., 1995; Deino et al., 2004). The activity during the last 15 ka was concentrated in three epochs (15-9.5, 8.6-8.2, 5.5-3.8 ka BP) separated by two quiescence periods (Di Vito

Sampling

We selected two key stratigraphic sections that well represent the volcanic stratigraphy of A-MS eruption: "Cavone degli Sbirri" (CS) and "Eremo dei Camaldoli" (EC) (Fig. 1A and Appendix A). These two sites are located at a distance of 2 and 5 km respectively from the eruptive center (Monte Spina). Both sections are located along the ENE-WSW axis of maximum dispersal of the main fallout deposits. The main fallout and PDC deposits have been sampled in the two sites. Bottom and top parts of the B1

Chemistry

The bulk rocks from the A-MS deposit ranges in composition from alkali-trachyte to trachyte (de Vita et al., 1999; Arienzo et al., 2010; Smith et al., 2011; Fig. 1c). Pumices are from near-aphyric to porphyritic (up to 13% of crystal in layer B1). Glass compositions of each A-MS layer and average analyses for all units are reported in Table 1 and Fig. 1c, respectively. The glass compositions, reported in Fig. 1c and Table 1, show a general decrease in silica content (from 62.3 to 59.8 wt% SiO2)

Numerical simulations of conduit dynamics during AM-S eruption

Based on the measured textural and chemical properties of the erupted products, we performed some 1-D numerical simulations using Confort 15 software (Campagnola et al., 2016a), in order to investigate the eruptive dynamics during magma rise and constrain the style and conditions of fragmentation, focusing in particular on the transitions from sustained to collapsing column. Confort 15 is an open source numerical model for steady-state magma flow through vertical and cylindrical conduits open

Conclusions

The 4.5 ka Agnano-Monte Spina eruption is the highest magnitude event of the past 5 ka at Campi Flegrei caldera. Despite the small volume of erupted magma (1 km3, Di Renzo et al., 2011), this event caused significant environmental impact over a very large area. According to different studies (e.g. Orsi et al., 2004; Bevilacqua et al., 2015), a future vent most likely will open in the area between the Agnano and San Vito plains and a future explosive eruption, in spite of the magnitude of the

Acknowledgements

This work was funded by the DPC-INGV Project V1-UR RM3 (resp. C. Romano). The Grant of Excellence Departments, MIUR-Italy (ARTICOLO 1, COMMI 314 - 337 LEGGE 232/2016) is also gratefully acknowledged. We thank Sergio Lo Mastro (LIME, Roma-tre, Italy) for the valuable technical support with SEM analyses. We also thank Stefania Sicola for valuable help with Raman data acquisition. Prof. Lucia Gurioli and an anonymous reviewer are strongly thanked for their thoughtful comments.

References (123)

  • G. De Natale et al.

    Seismicity and 3-D substructure at Somma–Vesuvius volcano: evidence for magma quenching

    Earth Planet. Sci. Lett.

    (2004)
  • S. de Vita et al.

    The Agnano – Monte Spina eruption (4100 years BP) in the restless Campi Flegrei caldera (Italy)

    J. Volcanol. Geotherm. Res.

    (1999)
  • M. de’ Michieli Vitturi et al.

    Transient effects of magma ascent dynamics along a geometrically variable dome-feeding conduit

    Earth Planet. Sci. Lett.

    (2010)
  • A. Deino et al.

    The age of the Neapolitan Yellow Tuff caldera-forming eruption (Campi Flegrei caldera-Italy) assessed by 40Ar/ 39Ar dating method

    J. Volcanol. Geotherm. Res.

    (2004)
  • C. Del Gaudio et al.

    Unrest episodes at Campi Flegrei: a reconstruction of vertical ground movements during 1905-2009

    J. Volcanol. Geotherm. Res.

    (2010)
  • P. Dellino et al.

    Interaction between particles transported by fallout and surge in the deposits of the Agnano–Monte Spina eruption (Campi Flegrei, Southern Italy)

    J. Volcanol. Geotherm. Res.

    (2004)
  • P. Dellino et al.

    Turbulent boundary layer shear flows as an approximation of base surges at Campi Flegrei (Southern Italy)

    J. Volcanol. Geotherm. Res.

    (2004)
  • D. Di Genova et al.

    Effect of iron and nanolites on Raman spectra of volcanic glasses: a reassessment of existing strategies to estimate the water content

    Chem. Geol.

    (2017)
  • V. Di Matteo et al.

    Water solubility in trachytic melts

    Chem. Geol.

    (2004)
  • A. Di Piazza et al.

    Unsteady magma discharge during the “El Retiro” subplinian eruption (Turrialba volcano, Costa Rica): insights from textural and petrological analyses

    J. Volcanol. Geotherm. Res.

    (2019)
  • V. Di Renzo et al.

    The magmatic feeding system of the Campi Flegrei caldera: architecture and temporal evolution

    Chem. Geol.

    (2011)
  • M.A. Di Vito et al.

    Volcanism and deformation since 12,000 years at the Campi Flegrei caldera (Italy)

    J. Volcanol. Geotherm. Res.

    (1999)
  • F. Dobran

    Nonequilibrium flow in volcanic conduits and application to the eruptions of Mt. St. Helens on May 18, 1980, and Vesuvius in AD 79

    J. Volcanol. Geotherm. Res.

    (1992)
  • R.V. Fisher et al.

    Mobility of a large-volume pyroclastic flow - emplacement of the Campanian ignimbrite

    Italy. J. Volcanol. Geotherm. Res.

    (1993)
  • B.F. Houghton et al.

    Diverse patterns of ascent, degassing, and eruption of rhyolite magma during the 1.8ka Taupo eruption, New Zealand: evidence from clast vesicularity

    J. Volcanol. Geotherm. Res.

    (2010)
  • E.W. Llewellin et al.

    Bubble suspension rheology and implications for conduit flow

    J. Volcanol. Geotherm. Res.

    (2005)
  • H.M. Mader et al.

    The rheology of two-phase magmas: a review and analysis

    J. Volcanol. Geotherm. Res.

    (2013)
  • L.G. Mastin et al.

    A multidisciplinary effort to assign realistic source parameters to models of volcanic ash-cloud transport and dispersion during eruptions

    J. Volcanol. Geotherm. Res.

    (2009)
  • D. Mele et al.

    Hazard of pyroclastic density currents at the Campi Flegrei Caldera (Southern Italy) as deduced from the combined use of facies architecture, physical modeling and statistics of the impact parameters

    J. Volcanol. Geotherm. Res.

    (2015)
  • N. Miyaji et al.

    High-resolution reconstruction of the Hoei eruption (AD 1707) of Fuji volcano, Japan

    J. Volcanol. Geotherm. Res.

    (2011)
  • D.J. Morgan et al.

    On estimating crystal shape for crystal size distribution analysis

    J. Volcanol. Geotherm. Res.

    (2006)
  • G. Orsi et al.

    The Neapolitan Yellow Tuff, a large-magnitude trachytic phreatoplinian eruption: eruptive dynamics, magma withdrawal and caldera collapse

    J. Volcanol. Geotherm. Res.

    (1992)
  • G. Orsi et al.

    Step-filling and development of a three-layer magma chamber: the Neapolitan Yellow Tuff case history

    J. Volcanol. Geotherm. Res.

    (1995)
  • G. Orsi et al.

    The restless, resurgent Campi Flegrei nested caldera (Italy): constraints on its evolution and configuration

    J. Volcanol. Geotherm. Res.

    (1996)
  • G. Orsi et al.

    Short-term ground deformations and seismicity in the resurgent Campi Flegrei caldera (Italy): an example of active block-resurgence in a densely populated area

    J. Volcanol. Geotherm. Res.

    (1999)
  • G. Orsi et al.

    Long-term forecast of eruption style and size at Campi Flegrei caldera (Italy)

    Earth Planet. Sci. Lett.

    (2009)
  • P. Papale et al.

    The role of magma composition and water content in explosive eruptions. 1. Conduit ascent dynamics

    J. Volcanol. Geotherm. Res.

    (1998)
  • C. Romano et al.

    The dry and hydrous viscosities of alkaline melts from Vesuvius and Phlegrean Fields

    Chem. Geol.

    (2003)
  • N.K. Adams et al.

    The physical volcanology of the 1600 eruption of Huaynaputina, southern Peru

    Bull. Volcanol.

    (2001)
  • N.K. Adams et al.

    Abrupt transitions during sustained explosive eruptions: examples from the 1912 eruption of Novarupta

    Alaska. Bull. Volcanol.

    (2006)
  • AGIP

    Geologia e Geofisica del sistema geotermico dei Campi Flegrei

    (1987)
  • F. Alfano et al.

    Insights into eruption dynamics from textural analysis: the case of the May, 2008, Chaitén eruption

    Bull. Volcanol.

    (2012)
  • A. Bevilacqua et al.

    Quantifying volcanic hazard at Campi Flegrei caldera (Italy) with uncertainty assessment: 1. Vent opening maps

    J. Geophys. Res. Solid Earth

    (2015)
  • A. Bevilacqua et al.

    The effects of vent location, event scale, and time forecasts on pyroclastic density current hazard maps at Campi Flegrei Caldera (Italy)

    Front. Earth Sci.

    (2017)
  • J. Blundy et al.

    Petrologic reconstruction of magmatic systems variables and processes

  • A. Burgisser et al.

    Experimental constraints on degassing and permeability in volcanic conduit flow

    Bull. Volcanol.

    (2005)
  • S. Campagnola

    Large-Scale Plinian Eruptions of the Colli Albani and the Campi Flegrei Volcanoes: Insights from Textural and Rheological Studies

    (2014)
  • S. Campagnola et al.

    Confort 15 model of conduit dynamics: applications to Pantelleria Green Tuff and Etna 122 BC eruptions

    Contrib. Mineral. Petrol.

    (2016)
  • S. Carey et al.

    Quantitative models of the fallout and dispersal of tephra from volcanic eruption columns

    Bull. Volcanol.

    (1986)
  • R.F. Cas et al.

    Volcanic Successions: Modern and Ancient

    (1987)
  • Cited by (13)

    • Micro-Raman water calibration in ultrapotassic silicate glasses: Application to phono-tephrites and K-foidites of Colli Albani Volcanic District (Central Italy)

      2022, Chemical Geology
      Citation Excerpt :

      The knowledge of water content in natural glasses is fundamental for understanding the numerous physical and chemical processes occurring before, during, and after the eruption (Bachmann et al., 2009; Berry et al., 2008; Blundy and Cashman, 2005; Di Genova et al., 2017; Dingwell, 2006; Hartley et al., 2014; Kennedy et al., 2005; Kotov et al., 2021; Métrich et al., 2010; Romano et al., 2020; Schiavi et al., 2018).

    • Microscale textural heterogeneity and tip-streaming instability in alkaline magmas: Evidence in tube pumices from Campi Flegrei, Italy

      2021, Journal of Volcanology and Geothermal Research
      Citation Excerpt :

      For magma quenched in AMS fallout pumices, pre-eruptive temperature, pressure and dissolved water content were ca. 900 °C, 100 MPa and 2–2.5 wt%, respectively [Piochi et al., 2008; Arienzo et al., 2010; Romano et al., 2020]. On the basis of temperature, water content and texturally determined crystal abundance (mostly 10% in both microvesicular and tube pumices and generally not more than 20% on average; see also Supplement), magma viscosity was estimated at ca. 104.5–105 Pa s [Piochi et al., 2008; Romano et al., 2020]. Melt viscosity referred to AMS products (0.7 wt% of water content) is at ca. 105.5 Pa s [Piochi et al., 2008].

    View all citing articles on Scopus
    View full text