skip to main content

Title: Crystals reveal magma convection and melt transport in dyke-fed eruptions

The processes and ranges of intensive variables that control magma transport and dyke propagation through the crust are poorly understood. Here we show that textural and compositional data of olivine crystals (Mg/Fe, Ni and P) from the tephra of the first months of Paricutin volcano monogenetic eruption (Mexico, 1943–1952) record fast growth and large temperature and oxygen fugacity gradients. We interpret that these gradients are due to convective magma transport in a propagating dyke to the Earth’s surface in less than a few days. The shortest time we have obtained is 0.1 day, and more than 50% of the calculated timescales are < 2 days for the earliest erupted tephra, which implies magma ascent rates of about 0.1 and 1 m s−1. The olivine zoning patterns change with the eruptive stratigraphy, and record a transition towards a more steady magma flow before the transition from explosive to effusive dynamics. Our results can inform numerical and experimental analogue models of dyke propagation, and thus facilitate a better understanding of the seismicity and other precursors of dyke-fed eruptions.

; ; ; ;
Publication Date:
Journal Name:
Scientific Reports
Nature Publishing Group
Sponsoring Org:
National Science Foundation
More Like this
  1. Hotspot tracks (quasilinear chains of seamounts, ridges, and other volcanic structures) provide important records of plate motions, as well as mantle geodynamics, magma flux, and mantle source compositions. The Tristan-Gough-Walvis Ridge (TGW) hotspot track, extending from the active volcanic islands of Tristan da Cunha and Gough through a province of guyots and then along Walvis Ridge to the Etendeka flood basalt province, forms one of the most prominent and complex global hotspot tracks. The TGW hotspot track displays a tight linear age progression in which ages increase from the islands to the flood basalts (covering ~135 My). Unlike Pacific tracks, which are simple chains of seamounts that are often compared to chains of pearls, the TGW track is alternately a steep-sided narrow ridge, an oceanic plateau, subparallel linear ridges and chains of seamounts, and areas of what appear to be randomly dispersed seamounts. The track displays isotopic zonation over the last ~70 My. The zonation appears near the middle of the track just before it splits into two to three chains of ridge- and guyot-type seamounts. The older ridge is also overprinted with age-progressive late-stage volcanism, which was emplaced ~30–40 My after the initial eruptions and has a distinct isotopicmore »composition. The plan for Expedition 391 was to drill at six sites, three along Walvis Ridge and three in the seamount (guyot) province, to gather igneous rocks to better understand the formation of track edifices, the temporal and geochemical evolution of the hotspot, and the variation in paleolatitudes at which the volcanic edifices formed. After a delay of 18 days to address a shipboard outbreak of the coronavirus disease 2019 (COVID-19) virus, Expedition 391 proceeded to drill at four of the proposed sites: three sites on the eastern Walvis Ridge around Valdivia Bank, an ocean plateau within the ridge, and one site on the lower flank of a guyot in the Center track, a ridge located between the Tristan subtrack (which extends from the end of Walvis Ridge to the island of Tristan da Cunha) and the Gough subtrack (which extends from Walvis Ridge to the island of Gough). One hole was drilled at Site U1575, located on a low portion of the northeastern Walvis Ridge north of Valdivia Bank. At this location, 209.9 m of sediments and 122.4 m of igneous basement were cored. The latter comprised 10 submarine lava units consisting of pillow, lobate, sheet, and massive lava flows, the thickest of which was ~21 m. Most lavas are tholeiitic, but some alkalic basalts were recovered. A portion of the igneous succession consists of low-Ti basalts, which are unusual because they appear in the Etendeka flood basalts but have not been previously found on Walvis Ridge. Two holes were drilled at Site U1576 on the west flank of Valdivia Bank. The first hole was terminated because a bit jammed shortly after penetrating igneous basement. Hole U1576A recovered a remarkable ~380 m thick sedimentary section consisting mostly of chalk covering a nearly complete sequence from Paleocene to Late Cretaceous (Campanian). These sediments display short and long cyclic color changes that imply astronomically forced and longer term paleoenvironmental changes. The igneous basement yielded 11 submarine lava units ranging from pillows to massive flows, which have compositions varying from tholeiitic basalt to basaltic andesite, the first occurrence of this composition recovered from the TGW track. These units are separated by seven sedimentary chalk units that range in thickness from 0.1 to 11.6 m, implying a long-term interplay of sedimentation and lava eruptions. Coring at Site U1577, on the extreme eastern flank of Valdivia Bank, penetrated a 154 m thick sedimentary section, the bottom ~108 m of which is Maastrichtian–Campanian (possibly Santonian) chalk with vitric tephra layers. Igneous basement coring progressed only 39.1 m below the sediment-basalt contact, recovering three massive submarine tholeiite basalt lava flows that are 4.1, 15.5, and >19.1 m thick, respectively. Paleomagnetic data from Sites U1577 and U1576 indicate that their volcanic basements formed just before the end of the Cretaceous Normal Superchron and during Chron 33r, shortly afterward, respectively. Biostratigraphic and paleomagnetic data suggest an east–west age progression across Valdivia Bank, becoming younger westward. Site U1578, located on a Center track guyot, provided a long and varied igneous section. After coring through 184.3 m of pelagic carbonate sediments mainly consisting of Eocene and Paleocene chalk, Hole U1578A cored 302.1 m of igneous basement. Basement lavas are largely pillows but are interspersed with sheet and massive flows. Lava compositions are mostly alkalic basalts with some hawaiite. Several intervals contain abundant olivine, and some of the pillow stacks consist of basalt with remarkably high Ti content. The igneous sequence is interrupted by 10 sedimentary interbeds consisting of chalk and volcaniclastics and ranging in thickness from 0.46 to 10.19 m. Paleomagnetic data display a change in basement magnetic polarity ~100 m above the base of the hole. Combining magnetic stratigraphy with biostratigraphic data, the igneous section is inferred to span >1 My. Abundant glass from pillow lava margins was recovered at Sites U1575, U1576, and U1578. Although the igneous penetration was only two-thirds of the planned amount, drilling during Expedition 391 obtained samples that clearly will lead to a deeper understanding of the evolution of the Tristan-Gough hotspot and its track. Relatively fresh basalts with good recovery will provide ample samples for geochemical, geochronologic, and paleomagnetic studies. Good recovery of Late Cretaceous and early Cenozoic chalk successions provides samples for paleoenvironmental study.« less
  2. Abstract

    Mafic volcanic activity is dominated by effusive to mildly explosive eruptions. Plinian and ignimbrite-forming mafic eruptions, while rare, are also possible; however, the conditions that promote such explosivity are still being explored. Eruption style is determined by the ability of gas to escape as magma ascends, which tends to be easier in low-viscosity, mafic magmas. If magma permeability is sufficiently high to reduce bubble overpressure during ascent, volatiles may escape from the magma, inhibiting violent explosive activity. In contrast, if the permeability is sufficiently low to retain the gas phase within the magma during ascent, bubble overpressure may drive magma fragmentation. Rapid ascent may induce disequilibrium crystallization, increasing viscosity and affecting the bubble network with consequences for permeability, and hence, explosivity. To explore the conditions that promote strongly explosive mafic volcanism, we combine microlite textural analyses with synchrotron x-ray computed microtomography of 10 pyroclasts from the 12.6 ka mafic Curacautín Ignimbrite (Llaima Volcano, Chile). We quantify microlite crystal size distributions (CSD), microlite number densities, porosity, bubble interconnectivity, bubble number density, and geometrical properties of the porous media to investigate the role of magma degassing processes at mafic explosive eruptions. We use an analytical technique to estimate permeability and tortuosity bymore »combing the Kozeny-Carman relationship, tortuosity factor, and pyroclast vesicle textures. The groundmass of our samples is composed of up to 44% plagioclase microlites, > 85% of which are < 10 µm in length. In addition, we identify two populations of vesicles in our samples: (1) a convoluted interconnected vesicle network produced by extensive coalescence of smaller vesicles (> 99% of pore volume), and (2) a population of very small and completely isolated vesicles (< 1% of porosity). Computed permeability ranges from 3.0 × 10−13to 6.3 × 10−12m2, which are lower than the similarly explosive mafic eruptions of Tarawera (1886; New Zealand) and Etna (112 BC; Italy). The combination of our CSDs, microlite number densities, and 3D vesicle textures evidence rapid ascent that induced high disequilibrium conditions, promoting rapid syn-eruptive crystallization of microlites within the shallow conduit. We interpret that microlite crystallization increased viscosity while simultaneously forcing bubbles to deform as they grew together, resulting in the permeable by highly tortuous network of vesicles. Using the bubble number densities for the isolated vesicles (0.1-3−3 × 104 bubbles per mm3), we obtain a minimum average decompression rate of 1.4 MPa/s. Despite the textural evidence that the Curacautín magma reached the percolation threshold, we propose that rapid ascent suppressed outgassing and increased bubble overpressures, leading to explosive fragmentation. Further, using the porosity and permeability of our samples, we estimated that a bubble overpressure > 5 MPa could have been sufficient to fragment the Curacautín magma. Other mafic explosive eruptions report similar disequilibrium conditions induced by rapid ascent rate, implying that syn-eruptive disequilibrium conditions may control the explosivity of mafic eruptions more generally.

    « less
  3. New bathymetric and gravity mapping, refined volume calculations and petrologic analyses show that the Hawaiian volcano P¯uh¯ahonu is the largest and hottest shield volcano on Earth. This ∼12.5-14.1 Ma volcano in the northwest Hawaiian Ridge (NWHR) is twice the size of Mauna Loa volcano (148 ±29 vs. 74.0 ×103km3), which was assumed to be not only the largest Hawaiian volcano but also the largest known shield volcano. We considered four testable mechanisms to increase magma production, including 1) thinner lithosphere, 2) slower propagation rate, 3) more fertile source, and 4) hotter mantle. The first three of these have been ruled out. The lithosphere was old (∼88 Myrs) when P¯uh¯ahonu was formed, and thus, too thick and cold to allow for greater extents of partial melting. The propagation rate was relatively fast when it erupted (87 km/Myr), so this is another unlikely reason. Source fertility was Kea-like and no more fertile than for other much smaller NWHR volcanoes. A hotter mantle remains the best mechanism to produce the large magma volumes and is consistent with the high forsteritic olivine phenocryst compositions (up to 91.8%) and the calculated high percent of melting (24%). Thus, the gargantuan size of P¯uh¯ahonu reflects its highmore »melting temperature, the highest reported for any Cenozoic basalt. A solitary wave within the Hawaiian plume is the probable cause of P¯uh¯ahonu’s higher melting temperature and the resulting increased volume flux given the absence of a more fertile source for P¯uh¯ahonu basalts, as found for many basalts from the Hawaiian Islands.« less
  4. Frey, Bonnie A. ; Kelley, Shari A. ; Zeigler, Kate E. ; McLemore, Virginia T. ; Goff, Fraser ; Ulmer-Scholle, Dana S. (Ed.)
    Monogenetic small-volume basaltic volcanoes are the most abundant subaerial volcanic landforms on Earth but are some of the most poorly understood systems. Their short durations, small volumes, and lack of recurrence make monitoring and hazard assessment difficult. The Zuni-Bandera volcanic field in western New Mexico contains small-volume basaltic centers erupting tholeiitic to alkalic basalts. Evidence shows no correlation of magma composition with eruption age, location, or volumetric output, prompting questions about the influence of magma ascent rates, magma storage conditions, and mantle source characteristics on lava compositions. Here, we present olivine major and minor element mineral chemistry from the 3200-year-old McCartys Flow, the youngest tholeiite basalt in the volcanic field. Olivine displays four phenocryst types with unique textures and major and minor element compositions. Multiple olivine types co-exist at the thin section scale. Major and minor element diffusion at frozen melt–phenocryst interfaces was modeled, revealing magma residence times ranging from 3–9 months. Type 3 olivine phenocrysts require step function initial conditions and record diffusion re-equilibration followed by magma mixing. These profiles indicate the magma resided in the reservoir for 10–15 years and accumulated from multiple batches of mixed magmas less than 10 days before the eruption. Our results show thatmore »primitive magmas in small-volume monogenetic volcanoes have complex lithospheric magmatic histories and stored in magma bodies influenced by an open system to develop different local chemical environments.« less
  5. Abstract

    Super-Earths are extremely common among the numerous exoplanets that have been discovered. The high pressures and temperatures in their interiors are likely to lead to long-lived magma oceans. If their electrical conductivity is sufficiently high, the mantles of Super-Earth would generate their own magnetic fields. With ab initio simulations, we show that upon melting, the behavior of typical mantle silicates changes from semi-conducting to semi-metallic. The electrical conductivity increases and the optical properties are substantially modified. Melting could thus be detected with high-precision reflectivity measurements during the short time scales of shock experiments. We estimate the electrical conductivity of mantle silicates to be of the order of 100 Ω−1 cm−1, which implies that a magnetic dynamo process would develop in the magma oceans of Super-Earths if their convective velocities have typical values of 1 mm/s or higher. We predict exoplanets with rotation periods longer than 2 days to have multipolar magnetic fields.