A Detection of Milankovitch Frequencies in Global Volcanic Activity
Total Page:16
File Type:pdf, Size:1020Kb
A detection of Milankovitch frequencies in global volcanic activity Steffen Kutterolf1*, Marion Jegen1, Jerry X. Mitrovica2, Tom Kwasnitschka1, Armin Freundt1, and Peter J. Huybers2 1Collaborative Research Center (SFB) 574, GEOMAR, Wischhofstrasse 1-3, 24148 Kiel, Germany 2Department of Earth & Planetary Sciences, Harvard University, Cambridge, Massachusetts 02138, USA ABSTRACT 490 k.y. in the Central American Volcanic Arc A rigorous detection of Milankovitch periodicities in volcanic output across the Pleistocene- (CAVA; Fig. 1). We augment these data with Holocene ice age has remained elusive. We report on a spectral analysis of a large number of 42 tephra layers extending over ~1 m.y. found well-preserved ash plume deposits recorded in marine sediments along the Pacifi c Ring of in Deep Sea Drilling Project (DSDP) and Fire. Our analysis yields a statistically signifi cant detection of a spectral peak at the obliquity Ocean Drilling Program (ODP) Legs offshore period. We propose that this variability in volcanic activity results from crustal stress changes of Central America. The marine tephra records associated with ice age mass redistribution. In particular, increased volcanism lags behind in Central America are dated using estimated the highest rate of increasing eustatic sea level (decreasing global ice volume) by 4.0 ± 3.6 k.y. sedimentation rates and/or through correlation and correlates with numerical predictions of stress changes at volcanically active sites. These with radiometrically dated on-land deposits results support the presence of a causal link between variations in ice age climate, continental (e.g., Kutterolf et al., 2008; also see the GSA stress fi eld, and volcanism. Data Repository1). The second data subset includes DSDP, INTRODUCTION Our data set is comprised of two subsets. The ODP, and IODP (Integrated Ocean Drill- Volcanic activity varies over a wide range of fi rst involves records from Central America and ing Program) drill core records from other temporal scales, from cycles with less than one- consists, in part, of 49 eruptions preserved in sites along the ROF (Fig. 1; Table DR1 in the year periods in single volcanic systems, to inter- tephra layers from 56 cross-correlated gravity Data Repository). The records extend back vals extending out to plate tectonic time scales cores. These data provide an almost complete to ca. 1 Ma, but rarely include the shallowest (Cambray and Cadet, 1994; Kennett et al., 1977; record of large explosive eruptions (VEI >5 sediments with ages up to ca. 100 ka, and have Mason et al., 2004; Paterne et al., 1990). Con- [Volcanic Explosivity Index]) over the past thus been complemented by data from gravity nections between volcanism, the global carbon cycle and climate appear to be well established at the longest of these time scales (Walker et al., 1981). Moreover, on seasonal-to-decadal time Aleutian Basin Bering scales, volcanic eruptions are known to infl u- Sea of Sea ence climate (e.g., Hansen et al., 1992; Robock, Asia Okhotsk Gulf of North Alaska 2000), and perhaps vice versa (Novell et al., Alaska America 2006; Rampino et al., 1979). We focus on the Kamchatka Aleutian Arc intermediate time scales (103–105 yr) relevant to Sea of Japan Westerly variations during the ice age. There is evidence stratospheric winds that subaerial volcanism increased signifi cantly during the last deglaciation phase of the ice age (Huybers and Langmuir, 2009; Jull and 15°N Nankai Mc Kenzie, 1996), and a connection between Philippines CAVA climate and volcanism over broader ice age time Pacific scales has been established in regional studies Ocean Easterly stratospheric winds Ecuador (e.g., Jellinek et al., 2004; McGuire et al., 1997; South America Novell et al., 2006; Paterne et al., 1990). How- Peru ever, a general link between glacial cycles and 15°S Tonga global volcanic activity has remained elusive. Australia Westerly TEPHRA RECORD FROM THE “RING stratospheric winds OF FIRE” We analyze marine records of widely dis- New Zealand persed, and well-preserved, tephra layers associated with the Circum-Pacifi c volcanic chain referred to as the “Ring of Fire” (ROF). Figure 1. Bathymetric map of Pacifi c Ocean with red triangles indicating active arc volca- The ROF accounts for about half the global noes along the Pacifi c “Ring of Fire.” Red dots mark DSDP/ODP/IODP (Deep Sea Drilling length of active plate subduction. Ash plumes Project /Ocean Drilling Program / Integrated Ocean Drilling Program) coring sites and Central of highly explosive Plinian eruptions travel American volcanic arc (CAVA) gravity core sites. ODP/DSDP sites are as follows: MD87–2121 and 1124 (New Zealand); 836 and 837 (Tonga Arc); 767–770 (Philippines); 1173 and 1174 (Nan- through the atmosphere and deposit volcanic kai Trough); 192 (Kamchatka); 178, 179, 189, and 190 (Aleutian Islands and Alaska); 1239 and tephras downwind from the eruption sites; 1237 (Ecuador and Peru); Legs 67, 170 and 205 (Central America). Arrows and dashed lines sites included in this study lie 200–500 km mark the predominant global stratospheric wind directions. from their respective sources in the prevailing stratospheric wind directions (Fig. 1). 1GSA Data Repository item 2013055, data, additional fi gures explaining the methods, and an extended er- ror analysis of the data set and statistics, is available online at www.geosociety.org/pubs/ft2013.htm, or on re- *E-mail: [email protected]. quest from [email protected] or Documents Secretary, GSA, P.O. Box 9140, Boulder, CO 80301, USA. GEOLOGY, February 2013; v. 41; no. 2; p. 227–230; Data Repository item 2013055 | doi:10.1130/G33419.1 | Published online 30 November 2012 GEOLOGY© 2012 Geological | February Society 2013 of America.| www.gsapubs.org For permission to copy, contact Copyright Permissions, GSA, or [email protected]. 227 cores. These tephra layers are also dated using ash sequence, and the peak within 14% of the time series is shown in Figure 2B (and by the estimated sedimentation rates. obliquity frequency that rises furthest above the black line in Fig. DR3D). In addition to the peak The time series of ash deposits at a specifi c background continuum in a fractional sense is within the obliquity band, the spectrum exhibits location is based on the individual drill site recorded. To build up a null distribution, this pro- some additional concentration of energy near with the most continuous and complete tephra cedure is repeated 100,000 times, and indicates the 1 per 100 k.y. and 1 per 23 k.y. frequencies. record. Preference is given to drill sites on the that the 6.3 ratio identifi ed for the actual 1 per Although these concentrations contain no sig- incoming plate since these are generally less 44 k.y. peak would be realized less than 1% of nifi cant peaks, the presence of excess energy is affected by erosion than sites on the continental the time, according to our null model. Because in accord with other indicators of late Pleisto- slope. The record is complemented by core data no other process is expected to concentrate cene climate. from adjacent sites if possible. Tephra layers of energy in the band of frequencies that we con- the same age and depositional depth obtained sider, these results constitute a statistically sig- LINKING MILANKOVITCH CYCLES, from sites in close proximity (<50 km) are gen- nifi cant detection of variations in tephra related ICE AGE DYNAMICS AND VOLCANISM erally counted as one eruptive event. Tephras of to the quasi-periodic 1 per 41 k.y. changes in the The foregoing spectral results strongly sug- similar age recorded at sites at greater distance obliquity of Earth’s spin axis. gest that ice-age climate variations induce volca- from one another are counted as separate events. To further evaluate the spectral peak that we nism, and the physical link may be the changes In total, we have identifi ed 408 tephra layers at associate with obliquity, we compress all tephra in crustal stress associated with ice-ocean sites along the ROF. The estimated ages of the ages by a factor of 10%. Compression has the mass redistributions during the glacial cycles tephra layers have uncertainties ranging to 14% largest effect on the absolute ages of the oldest (Glazner et al., 1999; Huybers and Langmuir, of their age. These uncertainties are accounted samples, for which uncertainties are the largest. 2009; Jull and McKenzie, 1996; McGuire et al., for in subsequent analyses. For a detailed dis- The eruption frequency computed on the basis 1997; Sigvaldason et al., 1992). Volcanic output cussion and explanation of the dating and asso- of this compressed time series is given in Fig- would be expected to increase during periods in ciated errors see the Data Repository. ure 2A. The associated power spectrum of this which the local confi ning pressure decreases, FREQUENCY ANALYSIS We mapped the ROF data set into a binary 0.8 0.008 time series in which the occurrence of a tephra A 0.6 0.007 layer is given a value of unity. To display low- 0.4 0.006 frequency variability, we applied a stable-phase, 0.2 0.005 16 0 0.004 running average low-pass fi lter with a width of 14 −0.2 0.003 5 k.y. to the time series (red line in Fig. DR3A in −0.4 0.002 frequency 12 −0.6 0.001 the Data Repository). It is instructive to consider in eruption Variation 10 Normalized obliquity −0.8 0 the spectral estimates associated with a sequence 1 2 3 4 5 6 7 8 9 8 Time (kyr) of progressively older, overlapping 400 k.y. 6 segments of the tephra time series (Fig. DR4), 4 (eruptions/k.y.) which we obtain using a multitaper spectral 2 Eruption frequency analysis (Thomson, 1982; Percival and Walden, 0 1998).