Milankovitch frequencies in tephra records at volcanic arcs: The relation of kyr-scale cyclic variations in volcanism to global climate changes

Highlights The increase in volcanic activity after the last glacial maximum observed on Iceland has led to one of the most fascinating hypothesis in science in the last decades: that deglaciation may force volcanism. We: - Re-analyzed four longer tephra records with the same statistical method and demonstrated that all contain the ∼41 kyr and ∼100 kyr Millankovitch periodicities. - The frequency spectra of the tephra and δ18O records are significantly correlated supporting the hypothesis that orbital-driven global climate changes interact with the volcanic eruption frequency regionally and globally. - However, the simultaneous analysis of the four best-characterized tephra records shows that correlations and associated time lags suffer from a number of uncertainties including the nature and quality of tephra time series, a wide range in geographic latitudes and geological settings, as well as applied statistical methods Therefore more precise tephra time series (preservation and age optimized) from different regions (glaciated versus non-glaciated) and geological settings (island arcs, continental arcs, intraplate) are needed together with standardized statistical analysis to decipher the impact of these factors on a global perspective of how climate may control volcanism. Abstract The increase in volcanic activity after the last glacial maximum observed on Iceland has led to one of the most fascinating hypothesis in science in the last decades: that deglaciation may force volcanism. Consequently, tephrostratigraphic records of sufficient length that cover multiple glacial cycles have been used to test whether such relationships hold systematically through the Quaternary. Here we review such tephra records that have been linked with climate proxy records such as δ18O in marine sediments, which is a measure of sea-level change and which is thought to be orbitally forced, as it exhibits the characteristic Milankovitch periodicities of precession (∼23 kyr), obliquity (∼41 kyr) and eccentricity (∼100 kyr). Statistical analyses have identified these periodicities also in long tephra records from different latitudes and geotectonic settings, as well as in compiled semi-global records. These studies detect Milankovitch periods in their tephra record, and also a phase shift relative to the δ18O record in such that periods of increased eruption frequencies coincide with the deglaciation period at the glacial/interglacial transition when ice and water loads on the lithosphere change most rapidly. However, there are also disparities in results and interpretations, which may be attributable to the different methods of analysis applied by the studies. We have therefore re-analyzed the four best-characterized tephra records by the same methods. We distinguish between analysis in the frequency domain, a novel approach, and analysis in the time domain, which has been used in previous studies. Analysis in the frequency domain identifies harmonic frequencies that arise from the binary nature of the tephra records and complicate the identification of primary frequencies. However, we show that all four records show spectral density peaks near the main Milankovitch periodicities of 41 and 100 kyr, and that they produce meaningful and significant statistical correlations with each other and the global δ18O record but not with random time series. Although the time-domain correlations with δ18O roughly confirm phase shifts implying peak volcanism during deglaciation, correlation coefficients arising from very noisy records are generally too low for precise constraints on the relative timing. These deficiencies presently hamper the recognition of the physical mechanisms through which global climate changes affect volcanism at both, high-latitude glaciated regions and low-latitude non-glaciated regions

100- kyr cyclicity in volcanic ash emplacement: evidence from a 1.1 Myr tephra record from the NW Pacific

It is a longstanding observation that the frequency of volcanism periodically changes at times of global climate change. The existence of causal links between volcanism and Earth's climate remains highly controversial, partly because most related studies only cover one glacial cycle. Longer records are available from marine sediment profiles in which the distribution of tephras records frequency changes of explosive arc volcanism with high resolution and time precision. Here we show that tephras of IODP Hole U1437B (northwest Pacific) record a cyclicity of explosive volcanism within the last 1.1 Myr. A spectral analysis of the dataset yields a statistically significant spectral peak at the similar to 100 kyr period, which dominates the global climate cycles since the Middle Pleistocene. A time-domain analysis of the entire eruption and delta O-18 record of benthic foraminifera as climate/sea level proxy shows that volcanism peaks after the glacial maximum and similar to 13 +/- 2 kyr before the delta O-18 minimum right at the glacial/interglacial transition. The correlation is especially good for the last 0.7 Myr. For the period 0.7-1.1 Ma, during the Middle Pleistocene Transition (MPT), the correlation is weaker, since the 100 kyr periodicity in the delta O-18 record diminishes, while the tephra record maintains its strong 100 kyr periodicity

100- kyr cyclicity in volcanic ash emplacement: evidence from a 1.1 Myr tephra record from the NW Pacific

It is a longstanding observation that the frequency of volcanism periodically changes at times of global climate change. The existence of causal links between volcanism and Earth’s climate remains highly controversial, partly because most related studies only cover one glacial cycle. Longer records are available from marine sediment profiles in which the distribution of tephras records frequency changes of explosive arc volcanism with high resolution and time precision. Here we show that tephras of IODP Hole U1437B (northwest Pacific) record a cyclicity of explosive volcanism within the last 1.1 Myr. A spectral analysis of the dataset yields a statistically significant spectral peak at the ~100 kyr period, which dominates the global climate cycles since the Middle Pleistocene. A time-domain analysis of the entire eruption and δ18O record of benthic foraminifera as climate/sea level proxy shows that volcanism peaks after the glacial maximum and ∼13 ± 2 kyr before the δ18O minimum right at the glacial/interglacial transition. The correlation is especially good for the last 0.7 Myr. For the period 0.7–1.1 Ma, during the Middle Pleistocene Transition (MPT), the correlation is weaker, since the 100 kyr periodicity in the δ18O record diminishes, while the tephra record maintains its strong 100 kyr periodicity

Miocene to Holocene marine tephrostratigraphy offshore northern Central America and southern Mexico: Pulsed activity of known volcanic complexes

We studied the tephra inventory of fourteen deep sea drill sites of three DSDP and ODP legs drilled offshore Guatemala and El Salvador (Legs 67, 84, 138), and one leg offshore Mexico (Leg 66). Marine tephra layers reach back from the Miocene to the Holocene. We identified 223 primary ash beds and correlated these between the drill sites, with regions along the volcanic arcs, and to specific eruptions known from land. In total, 24 correlations were established between marine tephra layers and to well‐known Quaternary eruptions from El Salvador and Guatemala. Additional 25 tephra layers were correlated between marine sites. Another 108 single ash layers have been assigned to source areas on land resulting in a total of 157 single eruptive events. Tephra layer correlations to independently dated terrestrial deposits provide new time markers and help to improve or confirm age models of the respective drill sites. Applying the respective sedimentation rates derived from the age models, we calculated ages for all marine ash beds. Hence, we also obtained new age estimates for eight known, but so far undated large terrestrial eruptions. Furthermore, this enables us to study the temporal evolution of explosive eruptions along the arc and we discovered five pulses of increased activity: 1) a pulse during the Quaternary, 2) a Pliocene pulse between 6 and 3 Ma, 3) a Late Miocene pulse between 10 and 7 Ma, 4) a Middle Miocene pulse between 17–11 Ma, and 5) an Early Miocene pulse (~>21 Ma)

Major changes in the post-glacial evolution of magmatic compositions and pre-eruptive conditions of Llaima Volcano, Andean Southern Volcanic Zone, Chile

Llaima is one of the most active volcanoes of the Chilean volcanic front with recent explosive eruptions in 2008 and 2009. Understanding how the volcano evolved to its present state is essential for predictions of its future behavior. The post-glacial succession of explosive volcanic eruptions of Llaima stratovolcano started with two caldera-forming eruptions at ∼16 and ∼15 ka, that emplaced two large-volume basaltic-andesitic ignimbrites (unit I). These are overlain by a series of fall deposits (unit II) changing from basaltic-andesitic to dacitic compositions with time. The prominent compositionally zoned, dacitic to andesitic Llaima pumice (unit III) was formed by a large Plinian eruption at ∼10 ka that produced andesitic surge deposits (unit IV) in its terminal phase. The following unit V represents a time interval of ∼8,000 years during which at least 30 basaltic to andesitic ash and lapilli fall deposits with intercalated volcaniclastic sediments and paleosols were emplaced. Bulk rock, mineral, and glass chemical data constrain stratigraphic changes in magma compositions and pre-eruptive conditions that we interpret in terms of four distinct evolutionary phases. Phase 1 (=unit I) magmas have lower large ion lithophile (LIL)/high field strength (HFS) element ratios compared to younger magmas and thus originated from a mantle source less affected by slab-derived fluids. They differentiated in a reservoir at mid-crustal level. During the post-caldera phase 2 (=units II–IV), relatively long residence times between eruptions allowed for increasingly differentiated magmas to form in a reservoir in the middle crust. Fractional crystallization led to volatile enrichment and oversaturation and is the driving force for the large Plinian eruption of the most evolved (unit III) dacite at Llaima, although replenishment by hot andesite probably triggered the eruption. During the subsequent phase 3 (=unit V >3 ka), frequent mafic replenishments at mid-crustal storage levels favored shorter residence times limiting erupted magma compositions to water-undersaturated basaltic andesites and andesites. At around 3 ka, the magma storage level for phase 4 (=unit V <3 ka to present) shifted to the uppermost crust where the hot magmas partly assimilated the granitic country rock. Although water contents of these basaltic andesites were low, the low-pressure storage facilitated water saturation before eruption. The change in magma storage level at 3 ka was responsible for the dramatic increase in eruption frequency compared to the older Llaima history. We suggest that the change from middle to upper crust magma storage is caused by a change in the stress regime below Llaima from transpression to tension