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

On mapping seafloor mineral deposits with central loop transient electromagnetics

Electromagnetic methods are commonly employed in exploration for land-based mineral deposits. A suite of airborne, land, and borehole electromagnetic techniques consisting of different coil and dipole configurations have been developed over the last few decades for this purpose. In contrast, although the commercial value of marine mineral deposits has been recognized for decades, the development of suitable marine electromagnetic methods for mineral exploration at sea is still in its infancy. One particularly interesting electromagnetic method, which could be used to image a mineral deposit on the ocean floor, is the central loop configuration. Central loop systems consist of concentric transmitting and receiving loops of wire. While these types of systems are frequently used in land-based or airborne surveys, to our knowledge neither system has been used for marine mineral exploration. The advantages of using central loop systems at sea are twofold: (1) simplified navigation, because the transmitter and receiver are concentric, and (2) simplified operation because only one compact unit must be deployed. We produced layered seafloor type curves for two particular types of central loop methods: the in-loop and coincident loop configurations. In particular, we consider models inspired by real marine mineral exploration scenarios consisting of overburdens 0 to 5 m thick overlying a conductive ore body 5 to 30 m thick. Modeling and resolution analyses showed that, using a 50 m(2) transmitting loop with 20 A of current, these two configurations are useful tools to determine the overburden depth to a conductive ore deposit and its thickness. In the most extreme case, absolute voltage errors on the order of 10 nV are required to resolve the base of a 30 m thick ore deposit. Whether such noise floors can be achieved in real marine environments remains to be seen

Structure of the mantle beneath the Alboran basin from magnetotelluric soundings

We present results of marine MT acquisition in the Alboran sea that also incorporates previously acquired land MT from southern Spain into our analysis. The marine data show complex MT response functions with strong distortion due to seafloor topography and the coastline, but inclusion of high resolution topography and bathymetry and a seismically defined sediment unit into a 3D inversion model has allowed us to image the structure in the underlying mantle. The resulting resistivity model is broadly consistent with a geodynamic scenario that includes subduction of an eastward trending plate beneath Gibraltar, which plunges nearly vertically beneath the Alboran. Our model contains three primary features of interest: a resistive body beneath the central Alboran, which extends to a depth of ~150 km. At this depth, the mantle resistivity decreases to values of ~100 Ohm-m, slightly higher than those seen in typical asthenosphere at the same depth. This transition suggests a change in slab properties with depth, perhaps reflecting a change in the nature of the seafloor subducted in the past. Two conductive features in our model suggest the presence of fluids released by the subducting slab or a small amount of partial melt in the upper mantle (or both). Of these, the one in the center of the Alboran basin, in the uppermost-mantle (20-30km depth) beneath Neogene volcanics and west of the termination of the Nekkor Fault, is consistent with geochemical models, which infer highly thinned lithosphere and shallow melting in order to explain the petrology of seafloor volcanics

Structure of the mantle beneath the Alboran Basin from Magnetotelluric Soundings

We present results of marine MT acquisition in the Alboran sea that also incorporates previously acquired land MT from southern Spain into our analysis. The marine data show complex MT response functions with strong distortion due to seafloor topography and the coastline, but inclusion of high resolution topography and bathymetry and a seismically defined sediment unit into a 3‐D inversion model has allowed us to image the structure in the underlying mantle. The resulting resistivity model is broadly consistent with a geodynamic scenario that includes subduction of an eastward trending plate beneath Gibraltar, which plunges nearly vertically beneath the Alboran. Our model contains three primary features of interest: a resistive body beneath the central Alboran, which extends to a depth of ∼150 km. At this depth, the mantle resistivity decreases to values of ∼100 Ohm‐m, slightly higher than those seen in typical asthenosphere at the same depth. This transition suggests a change in slab properties with depth, perhaps reflecting a change in the nature of the seafloor subducted in the past. Two conductive features in our model suggest the presence of fluids released by the subducting slab or a small amount of partial melt in the upper mantle (or both). Of these, the one in the center of the Alboran basin, in the uppermost‐mantle (20-30 km depth) beneath Neogene volcanics and west of the termination of the Nekkor Fault, is consistent with geochemical models, which infer highly thinned lithosphere and shallow melting in order to explain the petrology of seafloor volcanics

Tectonic Controls on Gas Hydrate Distribution off SW Taiwan

The northern part of the South China Sea is characterized by widespread occurrence of bottom simulating reflectors (BSR) indicating the presence of marine gas hydrate. Because the area covers both a tectonically inactive passive margin and the termination of a subduction zone, the influence of tectonism on the dynamics of gas hydrate systems can be studied in this region. Geophysical data show that there are multiple thrust faults on the active margin while much fewer and smaller faults exist in the passive margin. This tectonic difference matches with a difference in the geophysical characteristics of the gas hydrate systems. High hydrate saturation derived from ocean bottom seismometer data and controlled source electromagnetic data and conspicuous high‐amplitude reflections in P‐Cable 3D seismic data above the BSR are found in the anticlinal ridges of the active margin. In contrast all geophysical evidence for the passive margin points to normal to low hydrate saturations. Geochemical analyses of gas samples collected at seep sites on the active margin show methane with heavy δ13C isotope composition, while gas collected at the passive margin shows light carbon isotope composition. Thus, we interpret the passive margin as a typical gas hydrate province fuelled by biogenic production of methane and the active margin gas hydrate system as a system that is fuelled not only by biogenic gas production but also by additional advection of thermogenic methane from the subduction system

Hydrothermal circulation at the Cleft-Vance overlapping spreading center : results of a magnetometric resistivity survey

We report on a magnetometric resistivity sounding carried out in the overlapping spreading center between the Cleft and Vance segments of the Juan de Fuca Ridge. The data collected reveal a strong three dimensionality in the crustal electrical resistivity structure on wavelengths of a few kilometers. Areas of reduced crustal electrical resistivities, with values approaching that of seawater, are seen beneath the neovolcanic zones of both active spreading centers. We interpret these reduced resistivities as evidence of active hydrothermal circulation within the uppermost 1 km of hot, young oceanic crust

Mental Prayer for Catholic Youth

Drovers a many have laid waste my vineyard, trampled down my lands; the land I loved so, turned into a lonely wilderness 1 Desolate have they made it, and desolate it mourns for me now; a very picture of desolation, and all for the want of men with heeding hearts. 1 How often the holy prophet Jeremiah has been quoted, for his words are timeless. If he walked through the streets of the modern world, would he change his statement in any way? Perhaps not, unless to add a further cry of dismay that the Son of God continues to be ignored, forgotten, denied, persecuted, and condemned. No wonder the world is desolate, Man is too busy to think about his final destination, in spite of the fact that he is in a terrific hurry to get there