Hydrothermal flow at oceanic spreading centres accounts for about ten per cent of all heat flux in the oceans and controls the thermal structure of young oceanic plates. It also influences ocean and crustal chemistry, provides a basis for chemosynthetic ecosystems, and has formed massive sulphide ore deposits throughout Earth's history. Despite this, how and under what conditions heat is extracted, in particular from the lower crust, remains largely unclear. Here we present high-resolution, whole-crust, two- and three-dimensional simulations of hydrothermal flow beneath fast-spreading ridges that predict the existence of two interacting flow components, controlled by different physical mechanisms, that merge above the melt lens to feed ridge-centred vent sites. Shallow on-axis flow structures develop owing to the thermodynamic properties of water, whereas deeper off-axis flow is strongly shaped by crustal permeability, particularly the brittle-ductile transition. About 60 per cent of the discharging fluid mass is replenished on-axis by warm (up to 300 degrees Celsius) recharge flow surrounding the hot thermal plumes, and the remaining 40 per cent or so occurs as colder and broader recharge up to several kilometres away from the axis that feeds hot (500-700 degrees Celsius) deep-rooted off-axis flow towards the ridge. Despite its lower contribution to the total mass flux, this deep off-axis flow carries about 70 per cent of the thermal energy released at the ridge axis. This combination of two flow components explains the seismically determined thermal structure of the crust and reconciles previously incompatible models favouring either shallower on-axis or deeper off-axis hydrothermal circulation.

[1] The lowermost section of the continental superdeep drill hole German Continental Deep Drilling Program (KTB) (south Germany) has been investigated for the first time by vertical seismic profiling (VSP). The new VSP samples the still accessible range of 6-8.5 km depth. Between 7 and 8.5 km depth, the drill hole intersects a major cataclastic fault zone which can be traced back to the Earth's surface where it forms a lineament of regional importance, the Franconian line. To determine the seismic properties of the crust in situ, in particular within and around this deep fault zone, was one of the major goals of the VSP. For the measurements a newly developed high-pressure/high-temperature borehole geophone was used that was capable of withstanding temperatures and pressures up to 260°C and 140 MPa, respectively. The velocity-depth profiles and reflection images resulting from the VSP are of high spatial resolution due to a small geophone spacing of 12.5 m and a broad seismic signal spectrum. Compared to the upper part of the borehole, we found more than 10% decrease of the P wave velocity in the deep, fractured metamorphic rock formations. P wave velocity is $5.5 km/s at 8.5 km depth compared to 6.0-6.5 km/s at more shallow levels above 7 km. In addition, seismic anisotropy was observed to increase significantly within the deep fracture zone showing more than 10% shear wave splitting and azimuthal variation of S wave polarization. In order to quantify the effect of fractures on the seismic velocity in situ we compared lithologically identical rock units at shallow and large depths: Combining seismic velocity and structural logs, we could determine the elastic tensors for three gneiss sections. The analysis of these tensors showed that we need fracture porosity in the percent range in order to explain seismic velocity and anisotropy observed within the fault zone. The opening of significant pore space around 8 km depth can only be maintained by differential tectonic stress combined with intense macroscopic fracturing. VSP reflection imaging based on PP and PS converted reflected waves showed that the major fault system at the KTB site is wider and more complex than previously known. The so-called SE1 reflection previously found in two-and three-dimensional surface seismic surveys corresponds to the top of an $1 km wide fault system. Its lower portion was not illuminated by surface seismic acquisition geometry. VSP imaging shows that the fault zone comprises two major and a number of smaller SE dipping fault planes and several conjugate fracture planes. The previously recognized upper fault plane is not associated with a strong velocity anomaly but indicates the depth below which the dramatic velocity decrease starts. Regarding the complexly faulted crustal section of the KTB site as a whole, we found that fluctuation spectra of rock composition and seismic velocity show similar patterns. We could verify that a significant amount of P wave energy is continuously converted into shear energy by forward scattering a...

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.

High-temperature (>300 ºC) off-axis hydrothermal systems found along the slow-spreading Mid-Atlantic Ridge are apparently consistently located at outcropping fault zones. While preferential flow of hot fluids along highly permeable, fractured rocks seems intuitive, such efficient flow inevitably leads to the entrainment of cold ambient seawater. The temperature drop this should cause is difficult to reconcile with the observed high-temperature black smoker activity and formation of associated massive sulfide ore deposits. Here we combine newly acquired seismological data from the high-temperature, off-axis Logatchev 1 hydrothermal field (LHF1) with numerical modeling of hydrothermal flow to solve this apparent contradiction. The data show intense off-axis seismicity with focal mechanisms suggesting a fault zone dipping from LHF1 toward the ridge axis. Our simulations predict high-temperature venting at LHF1 only for a limited range of fault widths and permeability contrasts, expressed as the fault's relative transmissibility (the product of the two parameters). The relative transmissibility must be sufficient to "capture" a rising hydrothermal plume and redirect it toward LHF1 but low enough to prevent extensive mixing with ambient cold fluids. Furthermore, the temperature drop associated with any high permeability zone in heterogeneous crust may explain why a significant part of hydrothermal discharge along slow-spreading ridges occurs at low temperatures.