Rapid ascent and emplacement of basaltic lava during the 2005–06 eruption of the East Pacific Rise at ca. 9°51′N as inferred from CO2 contents

© The Author(s), 2016. This is the author's version of the work and is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Earth and Planetary Science Letters 453 (2016): 152-160, doi:10.1016/j.epsl.2016.08.007.Eruption rates at the mid–ocean ridges (MORs) are believed to strongly control the morphology and length of lava flows emplaced along the ridge axis, and thus the structure and porosity of the upper oceanic crust. Eruption rate also represents one of the few tools to gain insight into the driving pressures within sub-ridge magmatic systems. As eruption rate is inferred to vary systematically along the global mid-ocean ridge, understanding of how to assess eruption rate in submarine systems and how it maps to observable features of the ridge axis would provide a powerful tool to understand Earth's largest magmatic system. Eruption rates at MORs are poorly constrained, however, because of a lack of direct observations, preventing the duration of an eruption to be quantified. This study uses decompression experiments of MORB samples and numerical modeling of CO2 degassing to reconstruct the timescales for magma ascent and lava emplacement during the 2005–06 eruption of the East Pacific Rise at ca. 9°51’N. Samples collected from the lava flow are all supersaturated in dissolved CO2 contents, but CO2 decreases with distance from the vent, presumably as a consequence of progressive CO2 diffusion into growing bubbles. Samples collected at the vent contain ~105 vesicles per cm3. Pieces of these samples were experimentally heated to 1225°C at high pressure and then decompressed at controlled rates. Results, plus those from numerical modeling of diffusive bubble growth, indicate that magma rose from the axial magma chamber to the seafloor in ≤1 hour and at a rate of ≥2–3 km hr-1. Our modeling, as validated by experimental decompression of MORB samples with ~106 vesicles cm-3, also suggests that CO2 degassed from the melt within ~10–100 minutes as the vesicular lava traveled ~ 1.7 km along the seafloor, implying a volumetric flow rate on order of 103–4 m3 s-1. Given an ascent rate of ≥0.2 m s-1, the width of a rectangular dike feeding the lava would have been ~1–2 meters wide. MORB samples from the Pacific ridge are generally more supersaturated in dissolved CO2 than those from slower spreading Atlantic and Indian ridges. Our results suggest that Pacific MORBs ascend to the seafloor faster than Atlantic or Indian MORBsThis project was partially funded by a grant to J.E.G. from the U.S. National Science Foundation (OCE-1333882).2017-08-2

NanoSIMS results from olivine-hosted melt embayments: Magma ascent rate during explosive basaltic eruptions

The explosivity of volcanic eruptions is governed in part by the rate at which magma ascends and degasses. Because the time scales of eruptive processes can be exceptionally fast relative to standard geochronometers, magma ascent rate remains difficult to quantify. Here we use as a chronometer concentration gradients of volatile species along open melt embayments within olivine crystals. Continuous degassing of the external melt during magma ascent results in diffusion of volatile species from embayment interiors to the bubble located at their outlets. The novel aspect of this study is the measurement of concentration gradients in five volatile elements (CO2, H2O, S, Cl, F) at fine-scale (5–10 μm) using the NanoSIMS. The wide range in diffusivity and solubility of these different volatiles provides multiple constraints on ascent timescales over a range of depths. We focus on four 100–200 μm, olivine-hosted embayments erupted on October 17, 1974 during the sub-Plinian eruption of Volcán de Fuego. H2O, CO2, and S all decrease toward the embayment outlet bubble, while F and Cl increase or remain roughly constant. Compared to an extensive melt inclusion suite from the same day of the eruption, the embayments have lost both H2O and CO2 throughout the entire length of the embayment. We fit the profiles with a 1-D numerical diffusion model that allows varying diffusivities and external melt concentrations as a function of pressure. Assuming a constant decompression rate from the magma storage region at approximately 220 MPa to the surface, H2O, CO2 and S profiles for all embayments can be fit with a relatively narrow range in decompression rates of 0.3–0.5 MPa/s, equivalent to 11–17 m/s ascent velocity and an 8 to 12 minute duration of magma ascent from ~ 10 km depth. A two stage decompression model takes advantage of the different depth ranges over which CO2 and H2O degas, and produces good fits given an initial stage of slow decompression (0.05–0.3 MPa/s) at high pressure (< 145 MPa), with similar decompression rates to the single-stage model for the shallower stage. The magma ascent rates reported here are among the first for explosive basaltic eruptions and demonstrate the potential of the embayment method for quantifying magmatic timescales associated with eruptions of different vigor

Nickel and helium evidence for melt above the core–mantle boundary

High ^(3)He/^(4)He ratios in some basalts have generally been interpreted as originating in an incompletely degassed lower-mantle source. This helium source may have been isolated at the core–mantle boundary region since Earth’s accretion. Alternatively, it may have taken part in whole-mantle convection and crust production over the age of the Earth; if so, it is now either a primitive refugium at the core–mantle boundary or is distributed throughout the lower mantle. Here we constrain the problem using lavas from Baffin Island, West Greenland, the Ontong Java Plateau, Isla Gorgona and Fernandina (Galapagos). Olivine phenocryst compositions show that these lavas originated from a peridotite source that was about 20 per cent higher in nickel content than in the modern mid-ocean-ridge basalt source. Where data are available, these lavas also have high ^(3)He/^(4)He. We propose that a less-degassed nickel-rich source formed by core–mantle interaction during the crystallization of a melt-rich layer or basal magma ocean, and that this source continues to be sampled by mantle plumes. The spatial distribution of this source may be constrained by nickel partitioning experiments at the pressures of the core–mantle boundary

Dynamic Soil-Foundation-Structure Interaction Analyses of Large Caissons

Large cellular reinforced concrete caissons exist as foundations of major long-span bridges across waterways in many parts of the country. This study was conducted to evaluate the important factors affecting the seismic response of large caissons. The paper presents the results of equivalent linear and non-linear analyses performed for a typical caisson idealized based on the cellular caisson at Pier W3 of the West San Francisco Bay Bridge subject to ground motion with a peak rock acceleration of 0.6 g. This caisson is 38.7 m (127 fi) long by 22.9 m (75 ft) wide submerged in about 32.6 m (107 ft) of water. It is embedded in 33.5 m (110 fi) of soil deposits and is founded on rock. Equivalent linear 3-D and 2-D analyses conducted in the direction of the short axis (longitudinal) were performed using a modified version of computer program SASSI. The results of these 3-D and 2-D analyses are similar. Non-linear analyses were performed for 2-D models using computer program FLAC. The results indicate that side gapping, base lifting, interface sliding, and soil yielding reduce the earth pressure, base bearing stress, caisson shear and bending moment, and caisson motions. However, the frequency characteristics of the responses appear to be relatively unaffected

Textural and geochemical constraints on andesitic plug emplacement prior to the 2004–2010 vulcanian explosions at Galeras volcano, Colombia

Hazardous sequences of vulcanian explosions are thought to result from the repeated emplacement and destruction of degassed, highly crystalline magma plugs in the shallow conduit of arc volcanoes. The processes governing the timing and magnitude of these explosions are thought to be related to magma ascent rate and efficiency of degassing and crystallisation. We study a rare suite of time-constrained ballistic bombs from the 2004–2010 period of activity of Galeras volcano to reconstruct magma plug architecture prior to six individual vulcanian explosions. We find that each plug was vertically stratified with respect to crystallinity, vesicularity and melt volatile content, melt composition and viscosity. We interpret this structure as resulting from multiple bubble nucleation events and degassing-driven crystallisation during multi-step ascent of the magma forming the plug, followed by spatially variable crystallisation within the plug under contrasting conditions of effective undercooling created by degassing. We propose that the shallow conduit evolved from more open degassing conditions during 2004–2008 to more closed conditions during 2009–2010. This resulted in explosions becoming smaller and less frequent over time during 2004–2008, then larger and more frequent over time during 2009–2010. This evolution was controlled by changing average ascent rates and is recorded by systematic changes in plagioclase microlite textures. Our results suggest that small volume vulcanian explosions (~ 105 m3) should generally be associated with longer repose times (hundreds of days) and produce ballistics characterised by small numbers of large, prismatic plagioclase microlites. Larger volume vulcanian explosions (~ 106 m3) should be associated with shorter repose times (tens of days) and produce ballistics characterised by high numbers of small, more tabular plagioclase microlites

Pre- and syn-eruptive degassing and crystallisation processes of the 2010 and 2006 eruptions of Merapi volcano, Indonesia

The 2010 eruption of Merapi (VEI 4) was the volcano’s largest since 1872. In contrast to the prolonged and effusive dome-forming eruptions typical of Merapi’s recent activity, the 2010 eruption began explosively, before a new dome was rapidly emplaced. This new dome was subsequently destroyed by explosions, generating pyroclastic density currents (PDCs), predominantly consisting of dark coloured, dense blocks of basaltic andesite dome lava. A shift towards open-vent conditions in the later stages of the eruption culminated in multiple explosions and the generation of PDCs with conspicuous grey scoria and white pumice clasts resulting from sub-plinian convective column collapse. This paper presents geochemical data for melt inclusions and their clinopyroxene hosts extracted from dense dome lava, grey scoria and white pumice generated during the peak of the 2010 eruption. These are compared with clinopyroxene-hosted melt inclusions from scoriaceous dome fragments from the prolonged dome-forming 2006 eruption, to elucidate any relationship between pre-eruptive degassing and crystallisation processes and eruptive style. Secondary ion mass spectrometry analysis of volatiles (H2O, CO2) and light lithophile elements (Li, B, Be) is augmented by electron microprobe analysis of major elements and volatiles (Cl, S, F) in melt inclusions and groundmass glass. Geobarometric analysis shows that the clinopyroxene phenocrysts crystallised at depths of up to 20 km, with the greatest calculated depths associated with phenocrysts from the white pumice. Based on their volatile contents, melt inclusions have re-equilibrated during shallower storage and/or ascent, at depths of ~0.6–9.7 km, where the Merapi magma system is interpreted to be highly interconnected and not formed of discrete magma reservoirs. Melt inclusions enriched in Li show uniform “buffered” Cl concentrations, indicating the presence of an exsolved brine phase. Boron-enriched inclusions also support the presence of a brine phase, which helped to stabilise B in the melt. Calculations based on S concentrations in melt inclusions and groundmass glass require a degassing melt volume of 0.36 km3 in order to produce the mass of SO2 emitted during the 2010 eruption. This volume is approximately an order of magnitude higher than the erupted magma (DRE) volume. The transition between the contrasting eruptive styles in 2010 and 2006 is linked to changes in magmatic flux and changes in degassing style, with the explosive activity in 2010 driven by an influx of deep magma, which overwhelmed the shallower magma system and ascended rapidly, accompanied by closed-system degassing