Armin Freundt

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.

Emplacement of small-volume pyroclastic flows at Laacher See (East-Eifel, Germany)

We distinguish three eruptive units of pyroclastic flows (T1, T2, and T3; T for trass) within the late Quaternary Laacher See tephra sequence. These units differ in the chemical/mineralogical composition of the essential pyroclasts ranging from highly differentiated phonolite in T1 to mafic phonolite in T3. T1 and T2 flows were generated during Plinian phases, and T3 flows during a late Vulcanian phase. The volume of the pyroclastic flow deposits is about 0.6 km3. The lateral extent of the flows from the source vent decreases from > 10 km (T1) to < 4.5 km (T3). In the narrow valleys north of Laacher See, the total thickness of the deposits exceeds 60 m.At least 19 flow units in T1, 6 in T2, and 4 in T3 can be recognized at individual localities. Depositional cycles of 2 to 5 flow units are distinguished in the eruptive units. Thickness and internal structure of the flow units are strongly controlled by topography. Subfacies within flow units such as strongly enriched pumice and lithic concentration zones, dust layers, lapilli pipes, ground layers, and lithic breccias are all compositionally related to each other by enrichment or depletion of clasts depending on their size and density in a fluidized flow. While critical diameters of coarse-tail grading were found to mark the boundary between the coarse nonfluidized and the finer fluidized grain-size subpopulations, we document the second boundary between the fluidized and the very fine entrained subpopulations by histograms and Rosin-Rammler graphs. Grain-size distribution and composition of the fluidized middle-size subpopulations remained largely unchanged during transport.Rheological properties of the pyroclastic flows are deduced from the variations in flow-unit structure within the valleys. T1 flows are thought to have decelerated from 25 m/s at 4 km to < 15 m/s at 7 km from the vent; flow density was probably 600–900 kg/m3, and viscosity 5–50 P. The estimated yield strength of the flows of 200– > 1000 N/m2 is consistent with the divergence of lithic size/distance curves from purely Newtonian models; the transport of lithics must be treated as in a Bingham fluid. The flow temperature probably decreased from T1 (300°–500°C) to T3 (<200°C).A large-scale longitudinal variation in the flow units from proximal through medial to distal facies dominantly reflects temporal changes during the progressive collapse of an eruption column. Only a small amount of fallout tephra was generated in the T1 phase of eruption. The pyroclastic flows probably formed from relatively low ash fountains rather than from high Plinian eruption columns.

The entrainment of high-viscosity magma into low-viscosity magma in eruption conduits

We report experiments on the flow of two fluids of contrasting viscosity through a pipe in which low-viscosity fluid occupies the center of the pipe. The volume flux of the low-viscosity fluid in the pipe increased during an experiment but did not reach 100% in most cases. The transition from high- to low-viscosity-dominated outflow involved a drop in pressure gradient and an increase in flow rate due to reduced viscous resistance in the pipe. Initially, the central flow was thin and parallel-sided, but as its diameter increased the flow became unstable. A sequence of instabilities was observed during the course of each experiment, both in time and as a function of height in the pipe. In the most commonly observed instability the central flow adopted a helical geometry. The transition from parallel-sided to unstable flow first appeared at the top of the pipe and propagated downwards against the flow. Axisymmetric instabilities originating at the pipe entrance were also observed. All forms of instability exhibited entrainment of viscous fluid into the faster moving central flow. Entrainment was extensive early in the existence of the central flow, but later on the volume flux of lower-viscosity fluid in the central flow rose more rapidly than the rate of entrainment and the proportion of lower-viscosity fluid increased with time. These compositional changes determined the viscosity of the central flow which was found to control its diameter and velocity. In banded pumice deposits, silicic pumice without mafic component is commonly erupted alongside banded pumice blocks. We infer that banded pumice may correspond to the central flow in our experiments, i. e., that viscous magma has been incorporated into less viscous melt, and that pure acid pumice is derived from the outer flow. Changes in eruption style may be caused by variations in pressure gradient and flow rate due to changes in the viscosity of the melt in the conduit. Varied mafic/silicic proportions and degree of mixing in magmatic associations are controlled by the bulk volume erupted, discharge rate, initial temperature difference and aspect ratio of the conduit.

Experiments on conduit flow and eruption behavior of basaltic volcanic eruptions

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, Eotvos, 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).

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.

Lithic-enriched segregation bodies in pyroclastic flow deposits of Laacher See Volcano (East Eifel, Germany)

Small-volume (ca. 0.6 km3) pyroclastic flow deposits at Laacher See contain lithic breccias and two types of ground layers that differ significantly in their structure and composition from the main body of flow units. Lithic breccia bodies, up to 3.5 m thick, containing up to 85 weight% lithic blocks, occur locally at various distances from the vent. The deposition of these breccias was apparently governed by the strong influence of paleomorphology on the dynamics of the pyroclastic flows. The breccias were deposited at three main changes in bottom gradient along the path of the pyroclastic flows. The accumulation of large lithics is explained: (a) by compression of flows on the rising bottom close to the vent; (b) by thinning of flows accelerating over a steep incline; (c) by deceleration of the pre-concentrated lower part of flows in hydraulic jumps; and (d) possibly by a stationary vortex at the inner bend of a valley curvature. Poorly sorted lithic-rich ground layers, laterally highly variable in internal structure and composition, are restricted to marginal regions of the pyroclastic flow deposits within deep and narrow valleys. They are interpreted as having formed due to the extreme roughness of the valley walls, enforcing irregular turbulent flow and intense fluidization of the flow head, in which density-dominated segregation of lithics occurred. Wellsorted lapilli-rich ground layers of constant lateral thickness were probably generated by a more regularly moving, less intensely fluidized head of pyroclastic flows in which size-dominated segregation was effective but density-segregation was minor. A model of the temporal and longitudinal evolution of a flow head is proposed. Close to the vent, the head is exclusively erosive. With increasing distance, erosive power declines and erosion is paralleled by ground layer formation under strong fluidization. Further from the vent, the head ceases to erode while fluidization is still sufficient for ground layer formation. When fluidization declines to a level ineffective for segregation, ground layers terminate while the head advances and only terminates when plug-flow dominates.

Petrogenesis of rhyolite‐trachyte‐basalt composite ignimbrite P1, Gran Canada, Canary Islands

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

Eruption and emplacement of a basaltic welded ignimbrite during caldera formation on Gran Canaria

The 14.1 Ma old composite ignimbrite cooling unit P1 (45 km3) on Gran Canaria comprises a lower mixed rhyolite-trachyte tuff, a central rhyolite-basalt mixed tuff, and a slightly rhyolite-contaminated basaltic tuff at the top. The basaltic tuff is compositionally zoned with (a) an upward change in basalt composition to higher MgO content (4.3–5.2 wt.%), (b) variably admixed rhyolite or trachyte (commonly <5 wt.%), and (c) an upward increasing abundance of basaltic and plutonic lithic fragments and cognate cumulate fragments. The basaltic tuff is divided into three structural units: (I) the welded basaltic ignimbrite, which forms the thickest part (c. 95 vol.%) and is the main subject of the present paper; (II) poorly consolidated massive, bomb- and block-rich beds interpreted as phreatomagmatic pyroclastic flow deposits; and (III) various facies of reworked basaltic tuff. Tuff unit I is a basaltic ignimbrite rather than a lava flow because of the absence of top and bottom breccias, radial sheet-like distribution around the central Tejeda caldera, thickening in valleys but also covering higher ground, and local erosion of the underlying P1 ash. A gradual transition from dense rock in the interior to ash at the top of the basaltic ignimbrite reflects a decrease in welding; the shape of the welding profile is typical for emplacement temperatures well above the minimum welding temperature. A similar transition occurs at the base where the ignimbrite was emplaced on cold ground in distal sections. In proximal sections the base is dense where it was emplaced on hot felsic P1 tuff. The intensity of welding, especially at the base, and the presence of spherical particles and of mantled and composite particles formed by accretion and coalescence in a viscous state imply that the flow was a suspension of hot magma droplets. The flow most likely had to be density stratified and highly turbulent to prevent massive coalescence and collapse. Model calculations suggest eruption through low pyroclastic fountains (<1000 m high) with limited cooling during eruption and turbulent flow from an initial temperature of 1160°C. The large volume of 26 km3 of erupted basalt compared with only 16 km3 of the evolved P1 magmas, and the extremely high discharge rates inferred from model calculations are unusual for a basaltic eruption. It is suggested that the basaltic magma was erupted and emplaced in a fashion commonly only attributed to felsic magmas because it utilized the felsic P1 magma chamber and its ring-fissure conduits. Evolution of the entire P1 eruption was controlled by withdrawal dynamics involving magmas differing in viscosity by more than four orders of magnitude. The basaltic eruption phase was initially driven by buoyancy of the basaltic magma at chamber depth and continued degassing of felsic magma, but most of the large volume of basalt magma was driven out of the reservoir by subsidence of a c. 10 km diameter roof block, which followed a decrease in magma chamber pressure during low viscosity basaltic outflow.

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.

Mixing of rhyolite, trachyte and basalt magma erupted from a vertically and laterally zoned reservoir, composite flow P1, Gran Canaria

The 14.1 Ma composite welded ignimbrite P1 (45 km3 DRE) on Gran Canaria is compositionally zoned from a felsic lower part to a basaltic top. It is composed of four component magmas mixed in vertically varying proportions: (1) Na-rhyolite (10 km3) zoned from crystal-poor to highly phyric; (2) a continuously zoned, evolved trachyte to sodic trachyandesite magma group (6 km3); (3) a minor fraction of Na-poor trachyandesite (<1 km3); and (4) nearly aphyric basalt (26 km3) zoned from 4.3 to 5.2 wt% MgO. We distinguish three sites and phases of mixing: (a) Mutual mineral inclusions show that mixing between trachytic and rhyolitic magmas occurred during early stages of their intratelluric crystallization, providing evidence for long-term residence in a common reservoir prior to eruption. This first phase of mixing was retarded by increasing viscosity of the rhyolite magma upon massive anorthoclase precipitation and accumulation. (b) All component magmas probably erupted through a ring-fissure from a common upper-crustal reservoir into which the basalt intruded during eruption. The second phase of mixing occurred during simultaneous withdrawal of magmas from the chamber and ascent through the conduit. The overall withdrawal and mixing pattern evolved in response to pre-eruptive chamber zonation and density and viscosity relationships among the magmas. Minor sectorial variations around the caldera reflect both varying configurations at the conduit entrance and unsteady discharge. (c) During each eruptive pulse, fragmentation and particulate transport in the vent and as pyroclastic flows caused additional mixing by reducing the length scale of heterogeneities. Based on considerations of magma density changes during crystallization, magma temperature constraints, and the pattern of withdrawal during eruption, we propose that eruption tapped the P1 magma chamber during a transient state of concentric zonation, which had resulted from destruction of a formerly layered zonation in order to maintain gravitational equilibrium. Our model of magma chamber zonation at the time of eruption envisages a basal high-density Na-poor trachyandesite layer that was overlain by a central mass of highly phyric rhyolite magma mantled by a sheath of vertically zoned trachyte-trachyandesite magma along the chamber walls. A conventional model of vertically stacked horizontal layers cannot account for the deduced density relationships nor for the withdrawal pattern.