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: 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: 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.

Description: The upper Quaternary pyroclastic flow deposits of Laacher See volcano show compositional and structural facies variations on four different scales: (1) eruptive units of pyroclastic flows, composed of many flow units; (2) depositional cycles of as many as five flow units; flow units containing (3) regional intraflow-unit facies; and (4) local intraflow-unit subfacies. These facies can be explained by successively overlapping processes beginning in the magma column and ending with final deposition. The pyroclastic flow deposits thus reflect major aspects of the eruptive history of Laacher See volcano: (a) drastic changes in eruptive mechanism due to increasing access of water to the magma chamber and (b) change in chemical composition and crystal and gas content as evacuation of a compositionally zoned magma column progressed. The four scales of facies result from four successive sets of processes: (1) differentiation in the magma column and external factors governing the mechanism of eruption; (2) temporal variations of factors inducing eruption column collapse; (3) physical conditions in the eruption column and the way in which its collapse proceeds; and (4) interplay of flow-inherent and morphology-induced transport mechanics.

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: 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: 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: 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: [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.

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.