Description: Paleo-climate records and geodynamic modelling indicate the existence of complex interactions between glacial sea level changes, volcanic degassing, and atmospheric CO2, which may have modulated the climate system's descent into the last ice age. Between ~85-70 ka, during an interval of decreasing axial tilt, the orbital component in global temperature records gradually declined, while atmospheric CO2, instead of continuing is long-term correlation with Antarctic temperature, remained relatively stable. Based on novel global geodynamic models and the joint interpretation of paleo-proxy data as well as biogeochemical simulations, we show that a sea level fall in this interval caused enhanced pressure-release melting in the uppermost mantle, which may have induced a surge in magma and CO2 fluxes from mid-ocean ridges and oceanic hotspot volcanoes. Our results reveal a hitherto unrecognised negative feedback between glaciation and atmospheric CO2 predominantly controlled by marine volcanism on multi-millennial (suborbital) timescales of ~ 5,000-15,000 years.

Description: Bathymetric data from oceanic transform faults and their associated fracture zones were compiled, providing high-resolution gridded seafloor topography. Data used in this compilation were open and archived at US American National Oceanographic and Atmospheric Administration (https://maps.ngdc.noaa.gov/viewers/bathymetry), Japan Agency for Marine-Earth Science and Technology (http://www.godac.jamstec.go.jp/darwin/e), and the German Datacenter for bathymetric data (https://www.bsh.de/EN/DATA/Oceanographic_Data_Center/Surveying_data/surveying_data_node). Data were processed and gridded using Multibeam System (https://www.mbari.org/products/research-software/mb-system) and can be displayed using Generic Mapping Tools (https://gmt.soest.hawaii.edu). All grids are in netCDF format. The compilation includes transform faults and fracture zones from the Northern and Southern East Pacific Rise, the Cosos-Nazca spreading center, Chile Rise and the Pacific Antarctic Ridge, the Southwest Indian Ridge, Central Indian Ridge and Southeast Indian Ridge as well as the Mid-Atlantic Ridge.

Description: Paleo-climate records and geodynamic modelling indicate the existence of complex interactions between glacial sea level changes, volcanic degassing and atmospheric CO2, which may have modulated the climate system's descent into the last ice age. Between ∼85 and 70 kyr ago, during an interval of decreasing axial tilt, the orbital component in global temperature records gradually declined, while atmospheric CO2, instead of continuing its long-term correlation with Antarctic temperature, remained relatively stable. Here, based on novel global geodynamic models and the joint interpretation of paleo-proxy data as well as biogeochemical simulations, we show that a sea level fall in this interval caused enhanced pressure-release melting in the uppermost mantle, which may have induced a surge in magma and CO2 fluxes from mid-ocean ridges and oceanic hotspot volcanoes. Our results reveal a hitherto unrecognized negative feedback between glaciation and atmospheric CO2 predominantly controlled by marine volcanism on multi-millennial timescales of ∼5,000-15,000 years.

Description: The division of the earth’s surface into continents and oceans is a consequence of plate tectonics but a geological paradox exists at continent-ocean boundaries. Continental plate is thicker and lighter than oceanic plate, floating higher on the mantle asthenosphere, but it can rift apart by thinning and heating to form new oceans. In theory, continental plate subsides in proportion to the amount it is thinned and subsequently by the rate it cools down. However, seismic and borehole data from continental margins like the Atlantic show that the upper surface of many plates remains close to sea-level during rifting, inconsistent with its thickness, and subsides after breakup more rapidly than cooling predicts. Here we use numerical models to investigate the origin and nature of this puzzling behaviour with data from the Kwanza Basin, offshore Angola. We explore an idea where the continental plate is made increasingly buoyant during rifting by melt produced and trapped in the asthenosphere. Using finite element simulation, we demonstrate that partially molten asthenosphere combined with other mantle processes can counteract the subsidence effect of thinning plate, keeping it elevated by 2-3 km until breakup. Rapid subsidence occurs after breakup when melt is lost to the embryonic ocean ridge.

Description: Previous studies have estimated that mantle serpentinization reactions generate H2 at a rate of 1010–1012 mol/yr along the global mid-ocean ridge (MOR) system. Here we present results of 3-D geodynamic simulations that predict rates of additional mantle serpentinization and H2 production at oceanic transform faults (OTF). We find that the extent and rate of mantle serpentinization increases with OTF length and is maximum at intermediate slip rates of 5 to 10 cm/yr. The additional global OTF-related production of H2 is found to be between 6.1 and 10.7 × 1011 mol/yr, which is comparable to the predicted background MOR rate of 4.1–15.0 × 1011 mol H2/yr. This points to oceanic transform faults as potential sites of intense fluid-rock interaction, where chemosynthetic life could be sustained by serpentinization reactions.

Description: We reanalyze existing paleodata of global mean surface temperature ΔTg and radiative forcing ΔR of CO2 and land ice albedo for the last 800,000 years to show that a state‐dependency in paleoclimate sensitivity S, as previously suggested, is only found if ΔTg is based on reconstructions, and not when ΔTg is based on model simulations. Furthermore, during times of decreasing obliquity (periods of land‐ice sheet growth and sea level fall) the multi‐millennial component of reconstructed ΔTg diverges from CO2, while in simulations both variables vary more synchronously, suggesting that the differences during these times are due to relatively low rates of simulated land ice growth and associated cooling. To produce a reconstruction‐based extrapolation of S for the future we exclude intervals with strong ΔTg‐CO2 divergence and find that S is less state‐dependent, or even constant (state‐independent), yielding a mean equilibrium warming of 2–4 K for a doubling of CO2.

Description: The conditions permitting mantle serpentinization during continental rifting are explored within 2-D thermotectonostratigraphic basin models, which track the rheological evolution of the continental crust, account for sediment blanketing effects, and allow for kinetically controlled mantle serpentinization processes. The basic idea is that the entire extending continental crust has to be brittle for crustal scale faulting and mantle serpentinization to occur. The isostatic and latent heat effects of the reaction are fully coupled to the structural and thermal solutions. A systematic parameter study shows that a critical stretching factor exists for which complete crustal embrittlement and serpentinization occurs. Increased sedimentation rates shift this critical stretching factor to higher values as sediment blanketing effects result in higher crustal temperatures. Sediment supply has therefore, through the temperature-dependence of the viscous flow laws, strong control on crustal strength and mantle serpentinization reactions are only likely when sedimentation rates are low and stretching factors high. In a case study for the Norwegian margin, we test whether the inner lower crustal bodies (LCB) imaged beneath the Møre and Vøring margin could be serpentinized mantle. Multiple 2-D transects have been reconstructed through the 3-D data set by Scheck-Wenderoth and Maystrenko (2011). We find that serpentinization reactions are possible and likely during the Jurassic rift phase. Predicted thicknesses and locations of partially serpentinized mantle rocks fit to information on LCBs from seismic and gravity data. We conclude that some of the inner LCBs beneath the Norwegian margin may be partially serpentinized mantle.

Description: We present quantitative modeling results for the effects of surface relief on hydrothermal convection at ocean-spreading centers investigating how vent site locations and subsurface flow patterns are affected by bathymetry induced sub-seafloor pressure variations. The model is based on a 2-D FEM solver for fluid flow in porous media and is used to simulate hydrothermal convection systematically in 375 synthetic studies. The results of these studies show that bathymetric relief has a profound effect on hydrothermal flow: bathymetric highs induce subsurface pressure variations that can deviate upwelling zones and favor venting at structural highs. The deviation angle from vertical upwelling can be expressed by a single linear dependence relating deviation angle to bathymetric slope and depth of the heat source. These findings are confirmed in two case studies for the East Pacific Rise at 9°30′N and Lucky Strike hydrothermal fields. In both cases, it is possible to predict the observed vent field locations only if bathymetry is taken into account. Our results thereby show that bathymetric relief should be considered in simulations of submarine hydrothermal systems and plays a key role especially in focusing venting of across axis hydrothermal flow onto the ridge axis of fast spreading ridges.

Description: The subduction of partially serpentinized oceanic mantle may potentially be the key geologic process leading to the regassing of Earth’s mantle and also has important consequences for subduction zone processes such as element cycling, slab deformation, and intermediate-depth seismicity. Little is known about the quantity of water that is retained in the slab during mantle serpentinization. Recent studies using thermodynamical and/or experimental models of subduction zone processes have assumed that the mantle is uniformly serpentinized to a depth determined from the equilibrium stability of serpentine minerals in P-T space. This approach yields an incomplete picture of the pattern of serpentinization that may occur during bending-related faulting; an initial state that is essential for quantifying subsequent dehydration processes. In order to provide further constraints on the pattern of hydration and the amount of water trapped in the subducting mantle, we build a 2-D reactive-flow model incorporating the kinetic rate-dependence of serpentinization based on experimental results. After simulating hydration processes at the trench outer-rise, we find that the water content in serpentinized mantle strongly depends on the age of the subducting lithosphere and subduction rate, with values ranging between 1.8x105 and 4.0x106 kgm-2 reactive water uptake into the subducting mantle column. Serpentinization also results in a reduction in surface heat flux towards the trench caused by advective downflow of seawater into the reaction region. Observed heat flow reductions are larger than the reduction due to the minimum-water downflow needed for partial serpentinization, predicting that active hydrothermal vents and chemosynthetic communities should also be associated with bend-fault serpentinization. Model results agree with previous studies that the lower plane of double Benioff zones can be generated due to dehydration of serpentinized mantle at depth. The depth-dependent pattern of serpentinization including reaction kinetics predicts a separation between the two Benioff planes consistent with seismic observations.