Steffen Kutterolf

Pacific offshore record of plinian arc volcanism in Central America: 1. Along-arc correlations

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.

Pacific offshore record of plinian arc volcanism in Central America: 2. Tephra volumes and erupted masses

[1] 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.

A detection of Milankovitch frequencies in global volcanic activity

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.

Combined bromine and chlorine release from large explosive volcanic eruptions: A threat to stratospheric ozone?

Large explosive volcanic eruptions inject gases, aerosols, and fi ne ashes into the stratosphere, potentially infl uencing 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 fl uorescence 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.

Lifetime and cyclicity of fluid venting at forearc mound structures determined by tephrostratigraphy and radiometric dating of authigenic carbonates

At convergent margins, fluids rise through the forearc in response to consolidation of the upper plate and dewatering of the subducting plate, and produce various cold-seep–related features on the seafloor (mud diapirs, mud mounds). At the Central American forearc, authigenic carbonates precipitated from rising fluids within such structures during active venting while typical mixed-mud sediments were ejected onto the surrounding seafloor where they became intercalated with normal pelagic background sediments, indicating that mud mounds evolved unsteadily through alternating active and inactive phases. Intercalated regional ash layers from Plinian eruptions at the Central American volcanic arc provide time marks that constrain the ages of mud ejection activity. U/Th dating of drill core samples of authigenic carbonate caps of mud mounds yields ages agreeing well with those constrained by ash layers and showing that carbonate caps grow inward rather than outward during active venting. Both dating approaches show that offshore Nicaragua and Costa Rica (1) active and inactive phases can occur simultaneously at neighboring mounds, (2) mounds along the forearc have individual histories of activity, but there are distinct time intervals when nearly all mounds have been active or inactive, (3) lifetimes of mounds reach several hundred thousand years, and (4) highly active periods last 10–50 k.y. with intervening periods of >10 k.y. of relative quiescence.

Tephra layers: A controlling factor on submarine translational sliding?

Submarine slope failures occur at all continental margins, but the processes generating different mass wasting phenomena remain poorly understood. Multibeam bathymetry mapping of the Middle America Trench reveals numerous continental slope failures of different dimensions and origin. For example, large rotational slumps have been interpreted to be caused by slope collapse in the wake of subducting seamounts. In contrast, the mechanisms generating translational slides have not yet been described. Lithology, shear strength measurements, density, and pore water alkalinity from a sediment core across a slide plane indicate that a few centimeters thick intercalated volcanic tephra layer marks the detachment surface. The ash layer can be correlated to the San Antonio tephra, emplaced by the 6000 year old caldera-forming eruption from Masaya-Caldera, Nicaragua. The distal deposits of this eruption are widespread along the continental slope and ocean plate offshore Nicaragua. Grain size measurements permit us to estimate the reconstruction of the original ash layer thickness at the investigated slide. Direct shear test experiments on Middle American ashes show a high volume reduction during shearing. This indicates that marine tephra layers have the highest hydraulic conductivity of the different types of slope sediment, enabling significant volume reduction to take place under undrained conditions. This makes ash layers mechanically distinct within slope sediment sequences. Here we propose a mechanism by which ash layers may become weak planes that promote translational sliding. The mechanism implies that ground shaking by large earthquakes induces rearrangement of ash shards causing their compaction (volume reduction) and produces a rapid accumulation of water in the upper part of the layer that is capped by impermeable clay. The water-rich veneer abruptly reduces shear strength, creating a detachment plane for translational sliding. Tephra layers might act as slide detachment planes at convergent margins of subducting zones, at submarine slopes of volcanic islands, and at submerged volcano slopes in lakes.

Geochemical approaches to the quantification of dispersed volcanic ash in marine sediment

Volcanic ash has long been recognized in marine sediment, and given the prevalence of oceanic and continental arc volcanism around the globe in regard to widespread transport of ash, its presence is nearly ubiquitous. However, the presence/absence of very fine-grained ash material, and identification of its composition in particular, is challenging given its broad classification as an “aluminosilicate” component in sediment. Given this challenge, many studies of ash have focused on discrete layers (that is, layers of ash that are of millimeter-to-centimeter or greater thickness, and their respective glass shards) found in sequences at a variety of locations and timescales and how to link their presence with a number of Earth processes. The ash that has been mixed into the bulk sediment, known as dispersed ash, has been relatively unstudied, yet represents a large fraction of the total ash in a given sequence. The application of a combined geochemical and statistical technique has allowed identification of this dispersed ash as part of the original ash contribution to the sediment. In this paper, we summarize the development of these geochemical/statistical techniques and provide case studies from the quantification of dispersed ash in the Caribbean Sea, equatorial Pacific Ocean, and northwest Pacific Ocean. These geochemical studies (and their sedimentological precursors of smear slides) collectively demonstrate that local and regional arc-related ash can be an important component of sedimentary sequences throughout large regions of the ocean.

Subduction initiation and ophiolite crust: new insights from IODP drilling

ABSTRACT International Ocean Discovery Program (IODP) Expedition 352 recovered a high-fidelity record of volcanism related to subduction initiation in the Bonin fore-arc. Two sites (U1440 and U1441) located in deep water nearer to the trench recovered basalts and related rocks; two sites (U1439 and U1442) located in shallower water further from the trench recovered boninites and related rocks. Drilling in both areas ended in dolerites inferred to be sheeted intrusive rocks. The basalts apparently erupted immediately after subduction initiation and have compositions similar to those of the most depleted basalts generated by rapid sea-floor spreading at mid-ocean ridges, with little or no slab input. Subsequent melting to generate boninites involved more depleted mantle and hotter and deeper subducted components as subduction progressed and volcanism migrated away from the trench. This volcanic sequence is akin to that recorded by many ophiolites, supporting a direct link between subduction initiation, fore-arc spreading, and ophiolite genesis.

Regional-scale input of dispersed and discrete volcanic ash to the Izu-Bonin and Mariana subduction zones

We have geochemically and statistically characterized bulk marine sediment and ash layers at Ocean Drilling Program Site 1149 (Izu-Bonin Arc) and Deep Sea Drilling Project Site 52 (Mariana Arc), and have quantified that multiple dispersed ash sources collectively comprise ~30-35% of the hemipelagic sediment mass entering the Izu-Bonin-Mariana subduction system. Multivariate statistical analyses indicate that the bulk sediment at Site 1149 is a mixture of Chinese Loess, a second compositionally distinct eolian source, a dispersed mafic ash, and a dispersed felsic ash. We interpret the source of these ashes as respectively being basalt from the Izu-Bonin Front Arc (IBFA) and rhyolite from the Honshu Arc. Sr-, Nd-, and Pb isotopic analyses of the bulk sediment are consistent with the chemical/statistical-based interpretations. Comparison of the mass accumulation rate of the dispersed ash component to discrete ash layer parameters (thickness, sedimentation rate, and number of layers) suggests that eruption frequency, rather than eruption size, drives the dispersed ash record. At Site 52, the geochemistry and statistical modeling indicates that Chinese Loess, IBFA, dispersed BNN (boninite from Izu-Bonin), and a dispersed felsic ash of unknown origin are the sources. At Site 1149 the ash layers and the dispersed ash are compositionally coupled, whereas at Site 52 they are decoupled in that there are no boninite layers, yet boninite is dispersed within the sediment. Changes in the volcanic and eolian inputs through time indicate strong arc- and climate-related controls.

Pacific offshore record of plinian arc volcanism in Central America: 3. Application to forearc geology

[1] Sediment gravity cores collected on the Pacific slope and incoming plate offshore Central America reach up to 400 ka back in time and contain numerous ash layers from plinian eruptions at the Central American Volcanic Arc. The compositionally distinct widespread ash layers form a framework of marker horizons that allow us to stratigraphically correlate the sediment successions along and across the Middle America Trench. Moreover, ash layers correlated with 26 known eruptions on land provide absolute time lines through these successions. Having demonstrated the correlations in part 1, we here investigate implications for submarine sedimentary processes. Average accumulation rates of pelagic sediment packages constrained by bracketing tephras of known age range from ∼1–6 cm/ka on the incoming plate to 30–40 cm/ka on the continental slope. There are time intervals in which the apparent pelagic sedimentation rates significantly vary laterally both on the forearc and on the incoming plate where steady conditions are usually expected. A period of unsteadiness at 17–25 ka on the forearc coincides with a period of intense erosion on land probably triggered by tectonic processes. Unsteady conditions on the incoming plate are attributed to bend faulting across the outer rise triggering erosion and resedimentation. Extremely low apparent sedimentation rates at time intervals >50–80 ka suggest stronger tectonic activity than during younger times and indicate bend faulting is unsteady on a longer timescale. Submarine landslides are often associated with ash layers forming structurally weak zones used for detachment. Ash beds constrain ages of >60 ka, ∼19 ka, and <6 ka for three landslides offshore Nicaragua. Phases of intense fluid venting at mud mounds produce typical sediments around the mound that become covered by normal pelagic sediment during phases of weak or no activity. Using intercalated ash layers, we determine for the first time the durations (several hundred to 9000 years) of highly active periods in the multistage growth history of mud mounds offshore Central America, which is essential to understand general mud-mound dynamics.