ALBERT

Description: Drilling at Integrated Ocean Drilling Program Site U1381 on the Cocos Ridge offshore Costa Rica recovered 67 primary Miocene (ca. 8 Ma to ca. 16.5 Ma) marine fallout ash layers. Geochemical, volcanological, and geological criteria link these ashes to Plinian eruptions that carried ash to at least 50–450 km from the Galápagos hotspot. These ash layers are the first documentation of highly explosive Miocene Galápagos hotspot volcanism. This volcanism is bimodal with two-thirds of the tephra layers generated by basaltic magmas (glass compositions 〈57 wt% SiO 2 ) and one-third by rhyolitic magmas. The temporal distribution of the tephra layers, inferred from sediment accumulation rates calibrated by 40 Ar/ 39 Ar and biostratigraphic ages, reveals a distinct increase in eruption frequency and hence increased volcanic activity of the Galápagos hotspot after 14 Ma which we interpret in the context of dynamic interaction between the Galápagos plume and spreading ridge.

Description: Lahars, debris flows, and sediment-rich floods are frequent and deadly hazards at all mountain-forming volcanoes. Their hazard potential is traditionally assessed through mass-conserving closed-system models, where peak conversion rates of potential energy to mechanical energy and hence maximum destruction potential are predicted to occur on the steepest volcano flanks. This belies evidence of extremely high-energy and deadly catastrophes caused by such flows at large distances from volcanoes. Here we use the first high-resolution record of a moving lahar to develop a new model of the temporally and spatially variable mass-flow structure. We show that bulk flow energy can grow dramatically in such systems over tens to hundreds of kilometers via momentum transfers from the lahar into water and particles along its path. We also demonstrate that dynamic transformations of such flows and their ultimate runout are primarily controlled by the mass flow front.

Description: Large explosive volcanic eruptions inject gases, aerosols, and fine ashes into the stratosphere, potentially influencing climate. Emissions of chlorine (Cl) and bromine (Br) from such large eruptions play an important role for catalytic destruction of ozone in the stratosphere, but hitherto the global effects of simultaneous catastrophic release of volcanic Br and Cl into the stratosphere have not been investigated. The Br release from 14 large explosive eruptions throughout Nicaragua covering an entire subduction zone segment in the past 70 ka was determined with petrologic methods. Melt inclusions in volcanic phenocrysts were analyzed using a new optimized synchrotron–X-ray fluorescence microprobe set-up. Single eruptions produced Br outputs of 4–600 kt, giving an average Br emission of 27 kt per eruption. Using the assumption that 10% of the emitted halogens reach the stratosphere, the average Br and Cl loading to the stratosphere would be 3 ppt and 1500 ppt, respectively, which together would account for 185% of the preindustrial equivalent effective stratospheric Cl loading. We thus conclude that many large tropical volcanic eruptions had and have the potential to substantially deplete ozone on a global scale, eventually forming future ozone holes.

Description: A rigorous detection of Milankovitch periodicities in volcanic output across the Pleistocene-Holocene ice age has remained elusive. We report on a spectral analysis of a large number of well-preserved ash plume deposits recorded in marine sediments along the Pacific Ring of Fire. Our analysis yields a statistically significant detection of a spectral peak at the obliquity period. We propose that this variability in volcanic activity results from crustal stress changes associated with ice age mass redistribution. In particular, increased volcanism lags behind the highest rate of increasing eustatic sea level (decreasing global ice volume) by 4.0 ± 3.6 k.y. and correlates with numerical predictions of stress changes at volcanically active sites. These results support the presence of a causal link between variations in ice age climate, continental stress field, and volcanism.

Description: The 14 Ma caldera-forming composite ignimbrite P1 on Gran Canaria (Canary Islands) represents the first voluminous eruption of highly differentiated magmas on top of the basaltic Miocene shield volcano. Compositional zonation of the ignimbrite is the result of vertically changing proportions of four component magmas, which were intensely mixed during eruption: (1) Crystal-poor to highly phyric rhyolite (∼10 km3), (2) sodic trachyandesite through mafic to evolved trachyte (∼6 km3), (3) Na-poor trachyandesite (〈1 km3), and (4) basalt zoned from 5.2 to 4.3 wt % MgO (∼26 km3). P1 basalt is composed of two compositionally zoned magma batches, B2 basalt and B3 basalt. B3 basalt is derived from a mantle source depleted in incompatible trace elements compared to the shield basalt source. Basaltic magmas were stored in a reservoir probably underplating the crust, in which zoned B2 basaltic magma formed by mixing of “enriched” (shield) and “depleted” (B3) mafic melts and subsequent crystal fractionation. Evolved magmas formed in a shallow crustal chamber, whereas intermediate magmas formed at both levels. Abundant pyroxenitic to gabbroid cumulates in P1 support crystal fractionation as the major differentiation process. On the basis of major and trace element modeling, we infer two contemporaneous fractional crystallization series: series I from “enriched” shield basalt through Na-poor trachyandesite to rhyolite, and series II from “depleted” P1 basalt through sodic trachyandesite to trachyte. Series II rocks were significantly modified by selective contamination involving feldspar (Na, K, Ba, Eu, Sr), zircon (Zr) and apatite (P, Y, rare earth elements) components; apatite contamination also affected series I Na-poor trachyandesite. Substantial sodium introduction into sodic trachyandesite is the main reason for the different major element evolution of the two series, whereas their different parentage is mainly reflected in the high field strength trace elements. Selective element contamination involved not only rapidly but also slowly diffusing elements as well as different saturation conditions. Contamination processes thus variably involved differential diffusion, partial dissolution of minerals, partial melt migration, and trace mineral incorporation. Magma mixing between trachyte and rhyolite during their simultaneous crystallization in the P1 magma chamber is documented by mutual mineral inclusions but had little effect on the compositional evolution of both magmas. Fe-Ti oxide thermometry yields magmatic temperatures of around 850°C for crystal-poor through crystal-rich rhyolite, ∼815°C for trachyte and ∼850°–900°C for the trachyandesitic magmas. High 1160°C for the basalt magma suggest its intrusion into the P1 magma chamber only shortly before eruption. The lower temperature for trachyte compared to rhyolite and the strong crustal contamination of trachyte and sodic trachyandesite support their residence along the walls of the vertically and laterally zoned P1 magma chamber. The complex magmatic evolution of P1 reflects the transient state of Gran Canaria's mantle source composition and magma plumbing system during the change from basaltic to silicic volcanism. Our results for P1 characterize processes operating during this important transition, which also occurs on other volcanic ocean islands.

Description: Multiphase flow in basaltic volcanic conduits is investigated using analog experiments and theoretical approaches. Depending on gas supply, large gas bubbles (gas slugs) may rise through basaltic magma in regimes of distinct fluid‐dynamical behavior: ascent of single slugs, supplied slugs fed from the gas source during ascent, and periodic slug flow. An annular flow regime commences at the highest gas supply rates. A first set of experiments demonstrates that the growth of gas slugs due to hydrostatic decompression does not affect their ascent velocity and that excess pressure in the slugs remain negligible. The applicability of theoretical formulae describing slug ascent velocity as a function of liquid and conduit properties is evaluated in a second set of experiments. A third set of experiments with continuous gas supply into a cylindrical conduit are scaled to basaltic conditions over Morton, Eotvös, Reynolds, and Froude numbers. Gas flow rate and liquid viscosity are varied over the whole range of flow regimes to observe flow dynamics and to measure gas and liquid eruption rates. Foam generation by slug bursting at the surface and partial slug disruption by wake turbulence can modify the bubble content and size distribution of the magma. At the transition from slug to annular flow, when the liquid bridges between the gas slugs disappear, pressure at the conduit entrance drops by ∼60% from the hydrostatic value to the dynamic‐flow resistance of the annular flow, which may trigger further degassing in a stored magma to maintain the annular flow regime until the gas supply is exhausted and the eruption ends abruptly. Magma discharge may also terminate when magma ascent is hindered by wall friction in long volcanic conduits and the annular gas flow erodes all magma from the conduit. Supplied slugs are found to reach much higher rise velocities than unsupplied slugs and to collapse to turbulent annular flow upon bursting at the surface. A fourth set of experiments uses a conduit partially blocked by built‐in obstacles providing traps for gas pockets. Once gas pockets are filled, rising gas slugs deform but remain intact as they move around obstacles without coalescence or significant velocity changes. Bursting of bubbles coalescing with trapped gas pockets causes pressure signals at least 3 orders of magnitude more powerful than gas pocket oscillation induced by passing liquid. Our experiments suggest a refined classification of Strombolian and Hawaiian eruptions according to time‐dependant behavior into sporadically pulsating lava fountains (driven by stochastic rise of single slugs), periodically pulsating lava fountains (resulting from slug flow), and quasi‐steady lava fountains (oscillating at the frequency of annular‐flow turbulence).

Description: We studied the tephra inventory of 18 deep sea drill sites from six DSDP/ODP legs (Legs 84, 138, 170, 202, 205, 206) and two IODP legs (Legs 334 and 344) offshore the southern Central American Volcanic Arc (CAVA). Eight drill sites are located on the incoming Cocos plate and ten drill sites on the continental slope of the Caribbean plate. In total we examined ∼840 ash-bearing horizons and identified ∼650 of these as primary ash beds of which 430 originated from the CAVA. Correlations of ash beds were established between marine cores and with terrestrial tephra deposits, using major and trace element glass compositions with respect to relative stratigraphic order. As a prerequisite for marine-terrestrial correlations we present a new geochemical data set for significant Neogene and Quaternary Costa Rican tephras. Moreover, new Ar/Ar ages for marine tephras have been determined and marine ash beds are also dated using the pelagic sedimentation rates. The resulting correlations and provenance analyses build a tephrochronostratigraphic framework for Costa Rica and Nicaragua that covers the last 〉8 Myr. We define 39 correlations of marine ash beds to specific tephra formations in Costa Rica and Nicaragua; from the 4.15 Ma Lower Sandillal Ignimbrite to the 3.5 ka Rincón de la Vieja Tephra from Costa Rica, as well as another 32 widely distributed tephra layers for which their specific region of origin along Costa Rica and Nicaragua can be constrained.

Description: We collected 56 marine gravity cores from the Pacific seafloor offshore Central America which contain a total of 213 volcanic ash beds. Ash-layer correlations between cores and with their parental tephras on land use stratigraphic, lithologic, and compositional criteria. In particular, we make use of our newly built database of bulk-rock, mineral, and glass major and trace element compositions of plinian and similarly widespread tephras erupted since the Pleistocene along the Central American Volcanic Arc. We thus identify the distal ashes of 11 Nicaraguan, 8 El Salvadorian, 6 Guatemalan, and 1 Costa Rican eruptions. Relatively uniform pelagic sedimentation rates allow us to determine ages of 10 previously undated tephras by their relative position between ash layers of known age. Linking the marine and terrestrial records yields a tephrostratigraphic framework for the Central American volcanic arc from Costa Rica to Guatemala. This is a useful tool and prerequisite to understand the evolution of volcanism at a whole-arc scale.

Description: Sediment gravity cores collected from the Pacific seafloor offshore Central America contain numerous distal ash layers from plinian-type eruptions at the Central American Volcanic Arc dating back to more than 200 ka. In part 1 of this contribution we have correlated many of those ash layers between cores and with 26 tephras on land. The marine ash layers cover areas of up to 106 km2 in the Pacific Ocean and represent a major fraction (60–90%) of the erupted tephra volumes because the Pacific coast lies within a few tens of kilometers downwind from the volcanic arc. Combining our own mapping efforts on land and published mapping results with our marine data yields erupted volumes of all major tephras along the arc that range from ∼1 to 420 km3. Recalculated to erupted magma mass, the widespread tephras account for 65% of the total magma output at the arc. Complementing our tephra data with published volumes of the arc volcanic edifices and volcano ages, we calculate the long-term average magma eruption rates for each volcano. Moreover, we use incompatible element variations to calculate the cumulate masses that were fractionated during variable degrees of differentiation. This yields a minimum estimate of long-term average magma production rate at each volcano, because intrusives without surface expression and losses by erosion are not accounted for. Peak magma production rates increase from Costa Rica to Guatemala, but there is considerable scatter within each region and large differences even between neighboring volcanoes.

Description: Pacific drill sites offshore Central America provide the unique opportunity to study the evolution of large explosive volcanism and the geotectonic evolution of the continental margin back into the Neogene. The temporal distribution of tephra layers established by tephrochonostratigraphy in Part 1 indicates a nearly continuous highly explosive eruption record for the Costa Rican and the Nicaraguan volcanic arc within the last 8 M.y. The widely distributed marine tephra layers comprise the major fraction of the respective erupted tephra volumes and masses thus providing insights into regional and temporal variations of large-magnitude explosive eruptions along the southern Central American Volcanic Arc (CAVA). We observe three pulses of enhanced explosive magmatism between 0-1 Ma at the Cordillera Central, between 1-2 Ma at the Guanacaste and at 〉3 Ma at the Western Nicaragua segments. Averaged over the long-term the minimum erupted magma flux (per unit arc length) is ∼0.017 g/ms. Tephra ages, constrained by Ar-Ar dating and by correlation with dated terrestrial tephras, yield time-variable accumulation rates of the intercalated pelagic sediments with four prominent phases of peak sedimentation rates that relate to tectonic processes of subduction erosion. The peak rate at 〉2.3 Ma near Osa particularly relates to initial Cocos Ridge subduction which began at 2.91±0.23 Ma as inferred by the 1.5 M.y. delayed appearance of the OIB geochemical signal in tephras from Barva volcano at 1.42 Ma. Subsequent tectonic re-arrangements probably involved crustal extension on the Guanacaste segment that favored the 2-1 Ma period of unusually massive rhyolite production.