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 ~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

Data report: marine tephra compositions in the deep drilling cores of the South China Sea, IODP Expeditions 349 and 367/368

We present geochemical major and trace element glass data for tephra samples from International Ocean Discovery Program (IODP) Expeditions 349 and 367/368 from four drilling sites in the South China Sea. Overall, we obtained data for 55 samples and identified 46 as tephra layers, with dominant volcanic glass shards in the component inventory (in the 63–125 µ fraction). In total, we performed 720 single glass shard analyses using an electron microprobe for major element compositions, as well as 130 single glass shard analyses using laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) for trace element compositions. The compositions of the samples range from basaltic, (trachy-) andesitic to trachytic, and rhyolitic and fall mainly into the calc-alkaline and K-rich calc-alkaline magmatic series. One sample falls into the shoshonitic series. Tephras from Expedition 349 Site U1431 span the whole compositional range, whereas tephras from the other sites are limited to rhyolitic composition. Tephra ages, calculated applying sedimentation rates, range to ~2 Ma at Site U1431, ~0.8 Ma at Expedition 367 Site U1499, ~0.6 Ma at Expedition 368 Site U1501, and ~0.9 Ma at Expedition 368 Site U1505

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

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)

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

Zircon and Melt Extraction From a Long‐Lived and Vertically Extensive Magma System Underneath Ilopango Caldera (El Salvador)

The Tierra Blanca (TB) eruptive suite comprises the last four major eruptions of Ilopango caldera in El Salvador (≤45 ka), including the youngest Tierra Blanca Joven eruption (TBJ; ∼106 km3): the most voluminous event during the Holocene in Central America. Despite the protracted and productive history of explosive silicic eruptions at Ilopango caldera, many aspects regarding the longevity and the prevailing physicochemical conditions of the underlying magmatic system remain unknown. Zircon 238U-230Th geochronology of the TB suite (TBJ, TB2, TB3, and TB4) reveals a continuous and overlapping crystallization history among individual eruptions, suggesting persistent melt presence in thermally and compositionally distinct magma reservoirs over the last ca. 80 kyr. The longevity of zircon is in contrast to previously determined crystallization timescales of <10 kyr for major mineral phases in TBJ. This dichotomy is explained by a process of rhyolitic melt segregation from a crystal-rich refractory residue that incorporates zircon, whereas a new generation of major mineral phases crystallized shortly before eruption. Ti-in-zircon temperatures and amphibole geothermobarometry suggest that rhyolitic melt was extracted from different storage zones of the magma reservoir as indicated by distinct but synchronous thermochemical zircon histories among the TB suite eruptions. Zircon from TBJ and TB2 suggests magma differentiation within deeper and hotter parts of the reservoir, whereas zircon from TB3 and TB4 instead hints at crystallization in comparatively shallower and cooler domains. The assembly of the voluminous TBJ magma reservoir was also likely enhanced by cannibalization of hydrothermally altered components as suggested by low-δ18O values in zircon (+4.5 ± 0.3‰)

Depositional setting, provenance and tectonic-volcanic setting of Eocene-Recent deep-sea sediments of the oceanic Izu-Bonin forearc, NW Pacific (IODP Expedition 352)

New biostratigraphical, geochemical, and magnetic evidence is synthesized with IODP Expedition 352 shipboard results to understand the sedimentary and tectono-magmatic development of the Izu–Bonin outer forearc region. The oceanic basement of the Izu–Bonin forearc was created by supra-subduction zone seafloor spreading during early Eocene (c. 50–51 Ma). Seafloor spreading created an irregular seafloor topography on which talus locally accumulated. Oxide-rich sediments accumulated above the igneous basement by mixing of hydrothermal and pelagic sediment. Basaltic volcanism was followed by a hiatus of up to 15 million years as a result of topographic isolation or sediment bypassing. Variably tuffaceous deep-sea sediments were deposited during Oligocene to early Miocene and from mid-Miocene to Pleistocene. The sediments ponded into extensional fault-controlled basins, whereas condensed sediments accumulated on a local basement high. Oligocene nannofossil ooze accumulated together with felsic tuff that was mainly derived from the nearby Izu–Bonin arc. Accumulation of radiolarian-bearing mud, silty clay, and hydrogenous metal oxides beneath the carbonate compensation depth (CCD) characterized the early Miocene, followed by middle Miocene–Pleistocene increased carbonate preservation, deepened CCD and tephra input from both the oceanic Izu–Bonin arc and the continental margin Honshu arc. The Izu–Bonin forearc basement formed in a near-equatorial setting, with late Mesozoic arc remnants to the west. Subduction-initiation magmatism is likely to have taken place near a pre-existing continent–oceanic crust boundary. The Izu–Bonin arc migrated northward and clockwise to collide with Honshu by early Miocene, strongly influencing regional sedimentation

Izu-Bonin-Mariana Rear Arc: The Missing Half of the Subduction Factory

4GT) lies in the western part of the Izu fore-arc basin, ~60 km east of the arc-front volcano Aogashima, ~170 km west of the axis of the Izu-Bonin Trench, 1.5 km west of Ocean Drilling Program (ODP) Site 792, and at 1776 meters below sea level (mbsl). It was drilled as a 150 m deep geotechnical test hole for potential future deep drilling (5500 meters below seafloor [mbsf]) at proposed Site IBM-4 using the D/V Chikyu. Core from Site U1436 yielded a rich record of Late Pleistocene explosive volcanism, including distinctive black glassy mafic ash layers that may record large-volume eruptions on the Izu arc front. Because of the importance of this discovery, Site U1436 was drilled in three additional holes (U1436B, U1436C, and U1436D), as part of a contingency operation, in an attempt to get better recovery on the black glassy mafic ash layers and enclosing sediments and to better constrain the thickness of the mafic ash layers. IODP Site U1437 is located in the Izu rear arc, ~330 km west of the axis of the IzuBonin Trench and ~90 km west of the arc-front volcanoes Myojinsho and Myojin Knoll, at 2117 mbsl. The primary scientific objective for Site U1437 was to characterize “the missing half of the subduction factory”; this was because numerous ODP/Integrated Ocean Drilling Program sites had been drilled in the arc to fore-arc region (i.e., ODP Site 782A Leg 126), but this was the first site to be drilled in the rear part of the Izu arc. A complete view of the arc system is needed to understand the formation of oceanic arc crust and its evolution into continental crust. Site U1437 on the rear arc had excellent core recovery in Holes U1437B and U1437D, and we succeeded in hanging the longest casing ever in the history of R/V JOIDES Resolution scientific drilling (1085.6 m) in Hole U1437E and cored to 1806.5 mbsf

South China Sea Rifted Margin Testing hypotheses for lithosphere thinning during continental breakup: Drilling at the South China Sea rifted margin

International Ocean Discovery Program Expedition 368 is the second of two consecutive cruises that form the South China Sea Rifted Margin program. Expeditions 367 and 368 share the common key objectives of testing scientific hypotheses of breakup of the northern South China Sea (SCS) margin and comparing its rifting style and history to other nonvolcanic or magma-poor rifted margins. Four primary sites were selected for the overall program: one in the outer margin high (OMH) and three seaward of the OMH on distinct, margin-parallel basement ridges. These three ridges are informally labeled A, B, and C. They are located within the continent-ocean transition (COT) zone ranging from the OMH to the interpreted steady-state oceanic crust (Ridge C) of the SCS. The main scientific objectives include 1. Determining the nature of the basement within crustal units across the COT of the SCS that are critical to constrain style of rifting, 2. Constraining the time interval from initial crustal extension and plate rupture to the initial generation of igneous ocean crust, 3. Constraining vertical crustal movements during breakup, and 4. Examining the nature of igneous activity from rifting to seafloor spreading. In addition, the sediment cores from the drill sites targeting primarily tectonic and basement objectives will provide information on the Cenozoic regional environmental development of the Southeast Asia margin. Expedition 368 was planned to drill at two primary sites (U1501 and U1503) at the OMH and Ridge C, respectively. However, based on drilling results from Expedition 367, Expedition 368 chose to insert an alternate site on Ridge A (Site U1502). In total, the expedition completed operations at four sites (U1501, U1502, U1504, and U1505). Site U1503, however, was not completed beyond casing to 990 m because of mechanical problems with the drilling equipment that limited the expedition from 25 May 2017 to the end of the expedition to operate with a drill string not longer than 3400 m. New alternate Site U1504 proposed during Expedition 367 met this condition. Site U1505 also met the operational constraints of the 3400 m drill string (total) and was an alternate site for the already drilled Site U1501. At Site U1501, we cored to 697.1 m in 9.4 days, with 78.5% recovery. We also drilled ahead for 433.5 m in Hole U1501D and then logged downhole data from 78.3 to 399.3 m. In 19.3 days at Site U1502, we penetrated 1679.0 m, set 723.7 m of casing and cored a total of 576.3 m with 53.5% recovery, and collected downhole log data from 785.3 to 875.3 m and seismic data through the 10¾ inch casing. At Site U1503, we penetrated 995.1 m, setting 991.5 m of 10¾ inch casing, but no cores were taken. At Site U1504, we took 40 rotary core barrel (RCB) cores over two holes. The cored interval between both holes was 277.3 m with 26.8% recovery. An 88.2 m interval was drilled in Hole U1504B. At Site U1505, we cored 668.0 m with 101.1% recovery. Logging data was collected from 80.1 to 341.2 m. Operations at this site covered 6.1 days. Except for Site U1505, we drilled to acoustic basement, which prior to the expedition, except for Site U1501, had been interpreted to be crystalline basement. A total of 6.65 days were lost due to mechanical breakdown or waiting on spare supplies for repair of drilling equipment. At Site U1501 on the OMH, coring ~45 m into the acoustic basement sampled highly lithified sandstone to conglomerate of presumed Mesozoic age overlain by siliciclastic Eocene pre- to synrift sediments of Oligocene age and topped by primarily carbonaceous postrift sediments of early Miocene to Pleistocene age. Site U1502 on Ridge A was cased to 723.7 m. At this site, we recovered 180 m of hydrothermally altered brecciated basalts comprising sheet and pillow lavas below deep-marine sediments of Oligocene to late Miocene age. Coring was not performed within the upper 380 m (~Pliocene-Pleistocene) at Site U1502. At Site U1503 on Ridge C, 991.5 m of casing was installed in preparation for the planned deep drilling to ~1800 m, but no coring was performed due to mechanical failures, and the site was abandoned without further activity. Coring at Site U1504 on the OMH ~45 km east of Site U1501 recovered metamorphic schist to gneiss (greenschist facies) below late Eocene (?) carbonate rocks (partly reef debris) and early Miocene to Pleistocene sediments. At Site U1505, we cored to 480.15 m through Pleistocene to late Oligocene mainly carbonaceous ooze followed at depth by early Oligocene to late Eocene siliciclastic sediments. Efforts were made at every drill site to correlate the core with the seismic data and seismic stratigraphic unconformities interpreted within the Eocene to Plio-Pleistocene sedimentary sequence prior to drilling. The predrilling interpretation of ages of these unconformities was in general confirmed by drilling results. As a result of the constraints on the length of drill string that could be deployed during the later part of Expedition 368, the secondary expedition objectives addressing the environmental history of the SCS and Southeast Asia received more focus than planned because these sites are located in shallower water depths and required less penetration depth. This forced change in emphasis, however, was without fatal consequences for the primary tectonic objectives. The two expeditions together provided solid evidence for a process of breakup that included vigorous synrift magmatism as opposed to the often-favored interpretation of the SCS margin as a magma-starved margin