ORIGINAL RESEARCH published: 19 November 2020 doi: 10.3389/feart.2020.574032
A System Dynamics Approach to Understanding the deep Magma Plumbing System Beneath Dominica (Lesser Antilles)
Clara Solaro 1, Hélène Balcone-Boissard 2*, Daniel Joseph Morgan 3, Georges Boudon 1, Caroline Martel 4 and Léa Ostorero 1
1Université de Paris, Institut de Physique du Globe de Paris, CNRS, Paris, France, 2Institut des Sciences de la Terre de Paris (ISTeP), UMR 7193, CNRS-Sorbonne Université, Paris, France, 3Institute of Geophysics and Tectonics, School of Earth & Environment, University of Leeds, Leeds, United Kingdom, 4Institut des Sciences de la Terre d’Orléans (ISTO), UMR 7327, Université d’Orléans-CNRS/INSU-BRGM, Orléans, France
To understand the dynamics of magmatic systems, one must first seek to characterize the time-dependent behavior of magma storage and ascent. Herein, we do this through a combination of the Crystal System Approach and careful study of Fe-Mg interdiffusion in orthopyroxene. This allows us to trace the pre-eruptive dynamics of magma plumbing Edited by: systems, both in space and time. We apply this novel approach on two large silicic Jacob B. Lowenstern, eruptions (about 3–5km3 DRE/eruption) that occurred in the central part of Dominica Volcano Disaster Assistance Program ∼ (USGS), United States Island (Lesser Antilles Arc): the eruptions of Layou ( 51 ka) from Morne Diablotins, and Reviewed by: Roseau (∼33 ka) from Morne Trois Pitons-Micotrin. For the Roseau eruption, two Dawn Ruth, magmatic environments (MEs) are identified on the basis of orthopyroxene United States Geological Survey, composition, with a dominant reverse-zoning pattern from 50 to 54 to 54–59 mol% United States Geoff Kilgour, enstatite (En), indicating interaction with hotter magma. For the Layou eruption, three GNS Science, New Zealand MEs are observed as represented by three populations of pyroxenes: En47-51, En51-53 *Correspondence: and En53-58. The normal-zoning pathway from En51-53 to En47-51 is significantly Hélène Balcone-Boissard helene.balcone_boissard@ registered by crystals, interpreted as convective mixing in a zoned reservoir. The reverse- sorbonne-universite.fr zoning pathway from En47-51 to En51-53 and also En53-58 is also significantly present, Specialty section: supporting the mixing within the zoned reservoir but also suggesting mixing with a hotter This article was submitted to magma, possibly stored in another part of a sub-volcanic mush. The crystal and glass Volcanology, a section of the journal compositions (melt inclusion and matrix glass) from both studied eruptions suggest Frontiers in Earth Science heating and mixing between different magma pockets located within the mush that were Received: 18 June 2020 the dominant process for mobilizing eruptible magma. In parallel, we constrain the Accepted: 14 October 2020 associated pre-eruptive timescales by modeling the diffusive relaxation of Fe-Mg Published: 19 November 2020 Citation: chemical gradients that originated within the zonation of the same orthopyroxene Solaro C, Balcone-Boissard Hè, crystals. Diffusion modeling was considered along the b-axis of 66 zoned Morgan DJ, Boudon G, Martel C and orthopyroxene crystals for these two eruptions, at a magmatic temperature of 850 ± Ostorero L (2020) A System Dynamics ° Approach to Understanding the deep 25 C. In light of these results, we propose that the Layou and the Roseau magma Magma Plumbing System Beneath reservoirs were rejuvenated and heated by ∼25–50°C about 10 years prior to eruption by Dominica (Lesser Antilles). the injection of an underplating, hotter magma, creating the observed dominant reverse- Front. Earth Sci. 8:574032. doi: 10.3389/feart.2020.574032 zoning patterns of the erupted orthopyroxenes. We thus have evidence that silicic mush
Frontiers in Earth Science | www.frontiersin.org 1 November 2020 | Volume 8 | Article 574032 Solaro et al. Dynamics Magma Plumbing System Dominica
can be re-mobilized over timescales of decades prior to eruption, as previously suggested for Santorini and Taupo volcanoes.
Keywords: orthopyroxenes, Fe-Mg interdiffusion, timescale, plumbing system, crystal system analysis, Dominica (Lesser Antilles), ignimbritic eruptions
INTRODUCTION (2004) and Chakraborty (2008). Because of the various diffusion ratesofelementsincrystals(Costa et al., 2020), this technique allows How magma is stored at depth and how magma bodies assemble investigation of timescales spanning minutes (Charlier et al., 2012) prior to eruption are two questions that are still widely debated. to thousands of years (Allan et al., 2017) and, being applicable to The architecture and associated dynamics of magma plumbing single crystals, it can yield timescales of specific processes. Several systems for silicic volcanoes have been in the last decade envisaged studies, covering various volcanic settings and involving different as transcrustal magma systems (Cashman et al., 2017)following compositions of magmas and eruptive dynamics have already been the pioneering work of Bachmann and Bergantz (2003). conducted with different element/crystal pairs, including olivine In the process of sequential accretion of sills by which large (Costa and Chakraborty, 2004; Costa and Morgan, 2010; Kahl et al., crustal magma reservoirs form (Annen, 2009; Bachmann and 2011; Kahl et al., 2015), plagioclase (Zellmer et al., 2003; Zellmer Huber, 2016; Blundy and Annen, 2016), the resident magma can et al., 2003; Ginibre et al., 2004; Ruprecht and Wörner, 2007; cool and partially crystallize in between successive magma pulses Saunders et al., 2010; Shcherbakov et al., 2011; Druitt et al., to become a crystal mush (Marsh, 2006). Crystallinity stays in the 2012), pyroxenes (Morgan et al., 2004; Allan et al., 2013; Barker 50–60% range and the solidus is not reached quickly because of et al., 2016), magnetite (Boudon et al., 2015)andquartz(Matthews the release of latent heat during crystallization (Marsh, 1989). et al., 2012; Chamberlain et al., 2014). This type of highly-crystallized crustal reservoir raises questions In the Lesser Antilles Arc, the volume of magma involved in regarding magma mobility and the timescale over which magma most of the eruptions is relatively low along the volcanic arc can be remobilized through the rejuvenation of the mush to (generally <0.5 km3) as evidenced by many historical and become eruptible melt. Several mechanisms have been proposed prehistorical eruptions (Lindsay et al., 2005). Dominica is an to enable mush rejuvenation and possible melt extraction: self- exception: it is the sole island of the arc with several active mixing (Couch et al., 2001), gas sparging (Bachmann and volcanoes. Together with the neighboring islands of Guadeloupe Bergantz, 2003; Bachmann and Huber, 2016; Huber et al., and Martinique, it concentrates half of the active volcanoes of the 2010), or convective stirring (Huber et al., 2009). In any case, arc within its central part, only accounting for 1/5 of its total length underplating mafic magmas bringing heat and/or volatiles are (Figure 1A). The pioneer works on the volcanic activity of considered as major players in the reactivation of mush. Heat and Dominica recognized a voluminous silicic eruption, dated volatile transfers change the rheological properties of the mush between 50 and 30 ka and called “the Roseau Tuff” (Sigurdsson, (buoyancy, crystallinity), favoring its mobility (formation of melt 1972; Carey and Sigurdsson, 1980). Approximatively 58 km3 Dense connections, volatile foams). The efficiency of the enthalpy Rock Equivalent (DRE) of magma has been first estimated for this (±mass) exchange with newly injected recharge magma of eruption, making it as the most voluminous eruption of the arc. higher enthalpy necessary to unlock a crystal mush highly Though recent works now recognized several eruptions, in depends not only on the size but also on the composition and particular in the central part of Dominica (Lindsay et al., 2005; temperature of the intruding magma body. Estimated timescales Smith et al., 2013; Boudon et al., 2017), with necessarily lower suggest that the reactivation timescale is much longer than the emitted volumes, these volumes still contrast with those recognized following homogenization timescale (Huber et al., 2012). for the other active volcanoes. For instance, the ongoing eruption of Detailed petrological investigations of the compositions of Soufrière Hill’s of Montserrat only emitted 1 km3 of magma crystals, their melt inclusions and the residual melt provide (Wadge et al., 2014) and the famous devastating eruption of strong constraints on the dynamics of the magmatic plumbing Montagne Pelée (Martinique) in 1902–1905 was responsible of system in terms of architecture and timescales. Although such the emission of only 0.2–0.3 km3 (Boudon, 1983). Dominica studies are major undertakings, they are highly informative (Ruth volcanic activity thus offers the unique opportunity to et al., 2018; Albert et al., 2019). Using the method developed by Palm investigate the magma dynamics in the plumbing system prior (2005), the pre-eruptive dynamics of recent eruptions of Mt. Etna to repeated large silicic eruptions, only recognized on this island. has been explored using the sequential variations in olivine crystal Here, the aim is to identify and characterize the pre-eruptive compositions, represented schematically in a box model and called disrupting events leading to two successive large silicic eruptions Crystal System Analysis (CSA) (Kahl et al., 2011; Kahl et al., 2015). (∼3–5km3 DRE; Boudon et al., 2017) in the central part of Coupling CSA with diffusion timescales can be used to show variable Dominica, Layou (∼51 ka) and Roseau (∼33 ka), and their evolution of magmatic interactions prior to eruption and showed associated timescales. For that purpose, we present an original variable magmatic pathway systematics between eruptions. In recent study using orthopyroxenes that combine a variant of the Crystal decades, intracrystalline diffusion modeling has been increasingly System Analysis method and intracrystalline diffusion modeling. used to study timescales of pre-eruptive processes, built on the More precisely, we performed a detailed petrological pioneering works of Costa et al. (2003), Costa and Chakraborty investigation of pumices ejected during the first eruptive
Frontiers in Earth Science | www.frontiersin.org 2 November 2020 | Volume 8 | Article 574032 Solaro et al. Dynamics Magma Plumbing System Dominica
FIGURE 1 | Dominica and studied ignimbritic eruptions. (A) Lesser Antilles arc, with Dominica in the central part of the arc. (B) Dominica and location of the sampled outcrops (stars) for Layou and Roseau ignimbrites, after Boudon et al. (2017). Red: deposits from Morne Trois Pitons–Micotrin eruptive center, with dark red deposits corresponding to Roseau eruption. Blue: deposits from Morne Diablotins eruptive center, with dark blue deposits corresponding to Layou eruption. (C) Layou Plinian fallout deposit. (D) Roseau fallout deposit. (E) Layou pyroclastic density current deposit. (F) Roseau pyroclastic density current deposit.
Plinian phase of these two major eruptions. We focus our analysis GEOLOGICAL CONTEXT AND PREVIOUS on orthopyroxenes, which crystallize over a wide range of WORK temperature, pressure and oxygen fugacity. Thus, orthopyroxene crystals may be good tracers of evolving crustal Geological Context storage conditions. This work builds upon a revised The Lesser Antilles Arc results from the westward subduction of chronostratigraphy of the eruptions (Boudon et al., 2017), the North-South American plates beneath the Caribbean plate, at experimental investigations of the pre-eruptive magma storage an average rate of 2 cm/yr (Wadge, 1984). The volcanic arc is conditions (Solaro et al., 2019), and a detailed study of major, 850 km long, and extends from the island of Saba in the north to trace, and volatile element compositions of melt inclusions Grenada in the south, with an eastward convexity. Dominica is (Balcone-Boissard et al., 2018). The results are discussed for located in the central part of the arc (Figure 1). their implications in terms of the more general framework of North of Dominica, the arc splits into two branches: the Outer the magma plumbing system architecture and dynamics beneath or Older arc (Eocene to mid-Oligocene) to the north-east and an Dominica. Furthermore, we assess the rates of silicic magma Inner or Recent arc (Miocene to Recent) trending north-west. rejuvenation prior to voluminous explosive volcanic eruptions. This jump in the northern sector has been ascribed to the ancient
Frontiers in Earth Science | www.frontiersin.org 3 November 2020 | Volume 8 | Article 574032 Solaro et al. Dynamics Magma Plumbing System Dominica subduction of an aseismic ridge in the center of the arc causing a 2019), with some high An spikes. Amphiboles in Layou are ferri- rupture in the subducting plate and a displacement of volcanic magnesio-hornblende with Mg# [Mg/(Mg + Fetot)] ranging from activity to the west in the early-Miocene (∼20 Ma; Bouysse and 0.52 to 0.60. Clinopyroxenes are augites with Mg# ranging from Westercamp, 1990). South of Dominica, the two branches merge, 0.48 to 0.71. Fe-Ti oxides are titanomagnetites (Mag71–75;Mag creating a single volcanic arc active since the Miocene. The active magnetite mol%) and ilmenites (Ilm87;Ilm ilmenite, mol%). The parts of the arc (Saba to Dominica in the North, and Dominica to orthopyroxene compositions, the focus of this study, are presented Grenada in the South) have twelve, currently-active volcanic hereafter. centers, the best-known being Soufrière (Guadeloupe), Montagne Pelée (Martinique), Soufrière Hills (Montserrat) and Melt Inclusions Soufrière (St.-Vincent). Conversely to the other islands, Only glassy Melt Inclusions (MIs) with no post-entrapment Dominica is characterized by four, potentially-active volcanic modifications have been investigated, as detailed in Balcone- centers, from North to South: Morne aux Diables, Morne Boissard et al. (2018). MIs are ubiquitous in orthopyroxenes Diablotins, Morne Trois Pitons-Micotrin and Morne Plat Pays (sizes up to ∼50 µm). In plagioclase, scarce MIs are smaller than Volcanic Complex (Lindsay et al., 2005; Figure 1B). 40 µm. The MIs are all rhyolitic regardless of their host crystal, Five ignimbritic eruptions are currently recognized on with some variations in Na2O+K2O and SiO2 contents. MIs are Dominica: Grand Bay, Roseau, Grand Fond, Layou and among the most water-rich yet recorded (up to 8 wt% H2O). Grande Savane ignimbrites (Lindsay et al., 2005; Smith et al., Carbon dioxide contents are generally low (<650 ppm), although 2013; Howe et al., 2014; Boudon et al., 2017; Figure 1B). These MIs trapped at high pressures contain higher CO2 contents (up to deposits have been recently re-investigated by Boudon et al. 3,500–4,000 ppm). Some low-pressure (<100 MPa) inclusions 14 – (2017) who derive new C ages for the Layou, Roseau and have elevated CO2 contents (up to 1,100 1,150 ppm), ∼ ∼ ∼ fl fl Grand Fond ignimbrites respectively at 51, 33 and 24 ka; suggestive of uxing of shallow magmas with CO2-rich uids. – 3 each eruption involved magma volumes of the order of 3 5km CO2-trace element systematics indicate that melts were volatile- DRE. On the basis of geological arguments and different saturated at the time of entrapment (Balcone-Boissard et al., geochemical signatures observed in Roseau magmas compared 2018). to Layou ones, the Roseau ignimbrite originates from Morne Trois Pitons–Micotrin. The Layou ignimbrite likely originated Magma Storage Conditions from the Morne Diablotins eruptive center. The deposit of each Phase equilibrium experiments starting from glasses of Roseau eruption is subdivided into a basal Plinian fallout (<10% of total and Layou bulk compositions have been performed in an erupted volumes) (Figures 1C,D), overlain by thick (up to 100 m) internally-heated pressure vessel, spanning a temperature pumiceous pyroclastic density current deposits that vary in grain range of ∼800–900°C, a pressure range of 200–400 MPa (under size, welding and thickness distribution (Figures 1E,F). water saturation and under-saturation), and fixed oxygen fugacity of ≈ NNO +0.6 (NNO Nickel- Nickel oxide buffer) (Solaro et al., 2019). By experimentally reproducing the natural melt and Previous Works: Petrology of Pumice Clasts phenocryst rim compositions, phenocryst assemblages, phase and Pre-eruptive Magma Storage Conditions proportions, and water contents in MIs, the magma storage Mineralogy conditions have been determined to be close to a pressure of Whole rock, melt inclusions, and residual glass compositions ∼ ± ° – 400 MPa, a temperature of 850 25 C with 7 8 wt% H2O (major and trace elements) have been studied together with the dissolved in melt. These parameters will be used for the crystal contents of the natural samples of the Roseau and Layou diffusion modeling below. eruptions (Boudon et al., 2017; Balcone-Boissard et al., 2018). Magma of the Roseau eruption is intermediate between andesite – ∼ – and dacite (62 63 wt% SiO2, 4.5 5.0 wt% Na2O+K2O) whereas MATERIAL AND METHODS > ∼ magma from the Layou eruption is dacitic ( 64 wt% SiO2, 5wt% Na2O+K2O). For both eruption deposits, the residual glass is Sample Types and Preparation – – rhyolitic (76.5 78.0 wt% SiO2,5.57.0 wt% Na2O+K2O). The Studied pumices and crystals come from the field campaigns samples are porphyritic with ∼30 ± 5 vol% phenocrysts (Boudon performed previously for chronostratigraphic, melt inclusions et al., 2017) consisting of plagioclase, orthopyroxene, amphibole, and experimental petrological studies (Boudon et al., 2017; clinopyroxene, Ti-magnetite and ilmenite embedded in a Balcone-Boissard et al., 2018; Solaro et al., 2019). The samples microlite-free residual glass (Solaro et al., 2019). Modal consist of white, homogeneous pumice clasts collected from the composition of the magma of each eruption shows that basal Plinian fallout deposits of Layou and Roseau (Figures plagioclase is the most abundant phenocryst (70%), followed by 1C,D). Whole rock samples were gently crushed to liberate orthopyroxene (15–17%), amphibole (up to 17% for Layou, though crystals and sieved. Orthopyroxene (751 for Layou, 701 for <1% for Roseau) and, in lesser proportions, clinopyroxene and Fe- Roseau) crystals were separated by hand-picking under Ti oxides. Rare quartz phenocrysts are also recognized in the Layou binocular microscope and mounted in epoxy resin. deposit. Plagioclase crystals show oscillatory zoning patterns and Orthopyroxenes were oriented during mounting to obtain b-c resorption textures, with anorthite content (An, expressed in mol sections needed for diffusion modeling along the b-axis of %) spanning a broad range between An42 to An90 (Solaro et al., crystals.
Frontiers in Earth Science | www.frontiersin.org 4 November 2020 | Volume 8 | Article 574032 Solaro et al. Dynamics Magma Plumbing System Dominica
FIGURE 2 | From crystal to greyscale profile. (A) High resolution backscattering image [Back-Scattered Electron (BSE) image] of an orthopyroxene oriented with the c-axis in N-S direction. The core and the rim parts of the crystal are identified; (B) Detailed BSE image of the main core-rim boundary; (C) Core-to-rim greyscale profile: averaged greyscale values over the selected area of the boundary (mean width: 20–50 µm) in (B) using the free software ImageJ; (D) Microprobe Mg# (Mg# Mg/ [Mg + Fetot)] core-to-rim traverse over the selected area of the crystal. Spatial resolution is of 2 µm. Uncertainty: 0.2 mol% on Mg; (E) Final calibrated profile and modeled profile resulting from the intercalibration of Back-Scattered Electron greyscale values on microprobe analyses, assuming a step function as the initial condition (dashed lines). Data are presented in Supplementary Tables S2, S3 in supplementary material.
Analytical Techniques Diffusion Modeling Textural Investigations The method (Figures 2A–E) is based on the modeling as a All crystals were investigated using Back-Scattered Electron diffusion profile of the core-to-rim traverse (Figure 2E) (BSE) imaging via Scanning Electron Microscopy (SEM; obtained by coupling a high analytical resolution profile by ZEISS-Supra 55, Sorbonne Université, Paris) with an Electron Probe Microanalyses (Figure 2D) and a high spatial accelerating voltage of 20 kV, a beam current of 8 nA and a resolution profile by SEM (Figure 2C), following the counting time in order to reduce noise impact on each image and methodology developed by Allan et al. (2013) and thus better identify the crystal zonation (scan speed: 7; number of Couperthwaite et al. (2020). passages per scan: 8; Figure 2). The SEM images were processed From each BSE image, a greyscale profile was extracted using the free software ImageJ (http://rsb.info.nih.gov/ij/). (Figure 2C), averaging over a selected area. In effect, this Selected zones on orthopyroxene crystals have been imaged averages together many individual pixel traverses across the under high resolution to obtain a high spatial resolution image, necessary to reduce thermal noise from the backscatter profile (Figures 2B C). detector (Figure 2B). Profiles have been measured across the c- axis, in the b-crystal direction and away from the terminal faces of Electron Probe Microanalyses crystals, which are more sensitive to recrystallization (Figure 2A). Major elements (Si, Ti, Al, Mn Fe, Mg, Ca, Na, K, P) in The modeling was accomplished using 1-D modeling, an orthopyroxene crystals were analyzed by electron microprobe acceptable procedure as the profiles are very short (Cameca SX-Five and SX-100, Service Camparis, Paris) with an (∼20–30 µm) compared to the scale of the crystal (∼300 µm), accelerating voltage of 15 kV, a beam current of 10 nA and a and diffusion is effectively normal to the long face with no focused beam. Counting time on peak and background was set at additional components. As the diffusion coefficient of Ganguly 10 s for all elements, except for Fe and Mg with 60 s for and Tazzoli (1994) is composition-dependent, modeled profiles orthopyroxenes. Zoned orthopyroxene crystals have been develop asymmetry that is more pronounced as the analyzed along 2 µm step rim-to-rim or rim-to-core traverses compositional contrast becomes greater. The modeling of such (in regions previously picked out via high-resolution BSE profiles can be eased by exploiting the property that 1-D diffusion imaging; Figure 2D). profiles proportionally share symmetry, so that if the
Frontiers in Earth Science | www.frontiersin.org 5 November 2020 | Volume 8 | Article 574032 Solaro et al. Dynamics Magma Plumbing System Dominica