Sensitivity of seafloor bathymetry to climate-driven fluctuations in mid-ocean ridge magma supply

Mid-ocean ridge unfaulting revealed by magmatic intrusions

Mid-ocean ridges (MORs) are quintessential sites of tectonic extension1-4, at which divergence between lithospheric plates shapes abyssal hills that cover about two-thirds of the Earth's surface5,6. Here we show that tectonic extension at the ridge axis can be partially undone by tectonic shortening across the ridge flanks. This process is evidenced by recent sequences of reverse-faulting earthquakes about 15 km off-axis at the Mid-Atlantic Ridge and Carlsberg Ridge. Using mechanical models, we show that shallow compression of the ridge flanks up to the brittle failure point is a natural consequence of lithosphere unbending away from the axial relief. Intrusion of magma-filled fractures, which manifests as migrating swarms of extensional seismicity along the ridge axis, can provide the small increment of compressive stress that triggers reverse-faulting earthquakes. Through bathymetric analyses, we further find that reverse reactivation of MOR normal faults is a widely occurring process that can reduce the amplitude of abyssal hills by as much as 50%, shortly after they form at the ridge axis. This 'unfaulting' mechanism exerts a first-order influence on the fabric of the global ocean floor and provides a physical explanation for reverse-faulting earthquakes in an extensional environment.

Magmatism controls global oceanic transform fault topography

Oceanic transform faults play an essential role in plate tectonics. Yet to date, there is no unifying explanation for the global trend in broad-scale transform fault topography, ranging from deep valleys to shallow topographic highs. Using three-dimensional numerical models, we find that spreading-rate dependent magmatism within the transform domain exerts a first-order control on the observed spectrum of transform fault depths. Low-rate magmatism results in deep transform valleys caused by transform-parallel tectonic stretching; intermediate-rate magmatism fully accommodates far-field stretching, but strike-slip motion induces across-transform tension, producing transform strength dependent shallow valleys; high-rate magmatism produces elevated transform zones due to local compression. Our models also address the observation that fracture zones are consistently shallower than their adjacent transform fault zones. These results suggest that plate motion change is not a necessary condition for reproducing oceanic transform topography and that oceanic transform faults are not simple conservative strike-slip plate boundaries.

Carbonate compensation depth drives abyssal biogeography in the northeast Pacific

Abyssal seafloor communities cover more than 60% of Earth's surface. Despite their great size, abyssal plains extend across modest environmental gradients compared to other marine ecosystems. However, little is known about the patterns and processes regulating biodiversity or potentially delimiting biogeographical boundaries at regional scales in the abyss. Improved macroecological understanding of remote abyssal environments is urgent as threats of widespread anthropogenic disturbance grow in the deep ocean. Here, we use a new, basin-scale dataset to show the existence of clear regional zonation in abyssal communities across the 5,000 km span of the Clarion-Clipperton Zone (northeast Pacific), an area targeted for deep-sea mining. We found two pronounced biogeographic provinces, deep and shallow-abyssal, separated by a transition zone between 4,300 and 4,800 m depth. Surprisingly, species richness was maintained across this boundary by phylum-level taxonomic replacements. These regional transitions are probably related to calcium carbonate saturation boundaries as taxa dependent on calcium carbonate structures, such as shelled molluscs, appear restricted to the shallower province. Our results suggest geochemical and climatic forcing on distributions of abyssal populations over large spatial scales and provide a potential paradigm for deep-sea macroecology, opening a new basis for regional-scale biodiversity research and conservation strategies in Earth's largest biome.

Influence of late Pleistocene sea-level variations on midocean ridge spacing in faulting simulations and a global analysis of bathymetry

We address the hypothesis that changes in Pleistocene sea level have consequences for abyssal-hill bathymetry. A model illustrates how reductions in magma production caused by rising sea level could trigger faults at 100-ky intervals at intermediate spreading-rate ridges (>2.3 cm/y) and at 41-ky intervals at faster spreading-rate ridges (>3.8 cm/y). Analysis of 17 different regional ridge systems gives characteristic length scales that closely align with the predictions from the faulting model. Furthermore, a robust spectral peak is found at the 41-ky obliquity period at fast-spreading ridges. Together, these results constitute strong evidence for a pervasive influence of Pleistocene variations in glaciation and sea level on the pattern of abyssal hills.

It is established that changes in sea level influence melt production at midocean ridges, but whether changes in melt production influence the pattern of bathymetry flanking midocean ridges has been debated on both theoretical and empirical grounds. To explore the dynamics that may give rise to a sea-level influence on bathymetry, we simulate abyssal hills using a faulting model with periodic variations in melt supply. For 100-ky melt-supply cycles, model results show that faults initiate during periods of amagmatic spreading at half-rates >2.3 cm/y and for 41-ky melt-supply cycles at half-rates >3.8 cm/y. Analysis of bathymetry across 17 midocean ridge regions shows characteristic wavelengths that closely align with the predictions from the faulting model. At intermediate-spreading ridges (half-rates >2.3 cm/y and ≤3.8 cm/y) abyssal hill spacing increases with spreading rate at 0.99 km/(cm/y) or 99 ky (n= 12; 95% CI, 87 to 110 ky), and at fast-spreading ridges (half-rates >3.8 cm/y) spacing increases at 38 ky (n= 5; 95% CI, 29 to 47 ky). Including previously published analyses of abyssal-hill spacing gives a more precise alignment with the primary periods of Pleistocene sea-level variability. Furthermore, analysis of bathymetry from fast-spreading ridges shows a highly statistically significant spectral peak (P < 0.01) at the 1/(41-ky) period of Earth’s variations in axial tilt. Faulting models and observations both support a linkage between glacially induced sea-level change and the fabric of the sea floor over the late Pleistocene.

Spatial and Temporal Variations in Crustal Structure in the Eurasian Basin

To understand the tectonic–magmatic history, crustal structure and crustal accretion mode of the Eurasian Basin in the Arctic, we calculated the crustal thickness, residual bathymetry (RB) and non-isostatic topography of the Eurasian Basin by using the latest bathymetry, free-air gravity anomaly, crustal age and sediment thickness data. The tectonic–magmatic process of the Eurasian Basin can be divided into two stages. During magnetic isochrons 24-13 (C24-C13), the crustal structure of the basin was characterized by regional variations, while the crustal structure of the basin was characterized by local variations during magnetic isochrons 13-0 (C13-0). On the whole, the western part of the basin had a thicker crustal thickness and higher RB than the eastern part of the basin during C24-C13, which should result from the northward movement of Greenland. During magnetic isochrons 24-20 (C24-C20), the crustal structure of the eastern part of the basin had abnormally strong asymmetry. We speculate that there may be mantle upwelling beneath the Kara Sea Shelf in the south of the Eastern Eurasian Basin, which provides a large amount of melt for the crustal accretion of the southern part of the Eastern Eurasian Basin. The melt focusing supply could generate abnormally thick crust (>7 km) during magma enhancement period. The Western Eurasian Basin had stronger spatial variability and more frequency asymmetric polarity reversal than the Eastern Eurasian Basin during magnetic isochrons 6-0 (C6-0). We attribute this to the inflow of the North Atlantic mantle.

Megafaunal variation in the abyssal landscape of the Clarion Clipperton Zone

The potential for imminent polymetallic nodule mining in the Clarion Clipperton Fracture Zone (CCZ) has attracted considerable scientific and public attention. This concern stems from both the extremely large seafloor areas that may be impacted by mining, and the very limited knowledge of the fauna and ecology of this region. The environmental factors regulating seafloor ecology are still very poorly understood. In this study, we focus on megafaunal ecology in the proposed conservation zone 'Area of Particular Environmental Interest 6' (study area centred 17°16'N, 122°55'W). We employ bathymetric data to objectively define three landscape types in the area (a level bottom Flat, an elevated Ridge, a depressed Trough; water depth 3950-4250 m) that are characteristic of the wider CCZ. We use direct seabed sampling to characterise the sedimentary environment in each landscape, detecting no statistically significant differences in particle size distributions or organic matter content. Additional seafloor characteristics and data on both the metazoan and xenophyophore components of the megafauna were derived by extensive photographic survey from an autonomous underwater vehicle. Image data revealed that there were statistically significant differences in seafloor cover by nodules and in the occurrence of other hard substrata habitat between landscapes. Statistically significant differences in megafauna standing stock, functional structuring, diversity, and faunal composition were detected between landscapes. The Flat and Ridge areas exhibited a significantly higher standing stock and a distinct assemblage composition compared to the Trough. Geomorphological variations, presumably regulating local bottom water flows and the occurrence of nodule and xenophyophore test substrata, between study areas may be the mechanism driving these assemblage differences. We also used these data to assess the influence of sampling unit size on the estimation of ecological parameters. We discuss these results in the contexts of regional benthic ecology and the appropriate management of potential mining activities in the CCZ and elsewhere in the deep ocean.

Fast Deposition of Aligning Edge-On Polymers for High-Mobility Ambipolar Transistors

Fast deposition of aligning ambipolar polymers for high-performance organic field-effect transistors (OFETs) and inverter circuits are highly desired for both scientific studies and industry applications. Here, large-area and ordered polymer films are prepared by a bar-coating method at a rate of 120 mm s^-1 in air. Atomic force microscopy and grazing-incidence wide-angle X-ray scattering analysis indicate uniform edge-on poly(fluoroisoindigo-difluorobithiophene-fluoroisoindigo-bithiophene) (PFIBI-BT) in 11.7 ± 1 nm film (≈5 layers). The elongated, uniformly oriented grains can reduce the adverse effects of the grain boundaries and facilitate charge transport in polymers. Furthermore, OFETs based on parallel film show high hole/electron mobilities up to 5.5/4.5 cm² V^-1 s^-1 , which are approximately nine times of the devices prepared by spin-coating. The gain of the inverter is as high as 174, which is one of the highest values in polymer inventers currently. These results demonstrate that the excellent bipolar performance of few-layer PFIBI-BT can be ensured while achieving the compatibility of the experimental process with industrial preparation.

Anomalous K-Pg-aged seafloor attributed to impact-induced mid-ocean ridge magmatism

Eruptive phenomena at all scales, from hydrothermal geysers to flood basalts, can potentially be initiated or modulated by external mechanical perturbations. We present evidence for the triggering of magmatism on a global scale by the Chicxulub meteorite impact at the Cretaceous-Paleogene (K-Pg) boundary, recorded by transiently increased crustal production at mid-ocean ridges. Concentrated positive free-air gravity and coincident seafloor topographic anomalies, associated with seafloor created at fast-spreading rates, suggest volumes of excess magmatism in the range of ~10⁵ to 10⁶ km³. Widespread mobilization of existing mantle melt by post-impact seismic radiation can explain the volume and distribution of the anomalous crust. This massive but short-lived pulse of marine magmatism should be considered alongside the Chicxulub impact and Deccan Traps as a contributor to geochemical anomalies and environmental changes at K-Pg time.

Magmatic pulse driven by sea-level changes associated with the Messinian salinity crisis

Between 5 and 6 million years ago, during the so-called Messinian salinity crisis, the Mediterranean basin became a giant salt repository. The possibility of abrupt and kilometre-scale sea-level changes during this extreme event is debated. Messinian evaporites could signify either deep- or shallow-marine deposits, and ubiquitous erosional surfaces could indicate either subaerial or submarine features. Significant and fast reductions in sea level unload the lithosphere, which can increase the production and eruption of magma. Here we calculate variations in surface load associated with the Messinian salinity crisis and compile the available time constraints for pan-Mediterranean magmatism. We show that scenarios involving a kilometre-scale drawdown of sea level imply a phase of net overall lithospheric unloading at a time that appears synchronous with a magmatic pulse from the pan-Mediterranean igneous provinces. We verify the viability of a mechanistic link between unloading and magmatism using numerical modelling of decompression partial mantle melting and dike formation in response to surface load variations. We conclude that the Mediterranean magmatic record provides an independent validation of the controversial kilometre-scale evaporative drawdown and sheds new light on the sensitivity of magmatic systems to the surface forcing.

Sea level fall during glaciation stabilized atmospheric CO2 by enhanced volcanic degassing

Paleo-climate records and geodynamic modelling indicate the existence of complex interactions between glacial sea level changes, volcanic degassing and atmospheric CO₂, 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 CO₂, 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 CO₂ fluxes from mid-ocean ridges and oceanic hotspot volcanoes. Our results reveal a hitherto unrecognized negative feedback between glaciation and atmospheric CO₂ predominantly controlled by marine volcanism on multi-millennial timescales of ∼5,000-15,000 years.

Response to Comment on "Sensitivity of seafloor bathymetry to climate-driven fluctuations in mid-ocean ridge magma supply"

Tolstoy reports the existence of a characteristic 100 thousand year (ky) period in the bathymetry of fast-spreading seafloor but does not argue that sea level change is a first-order control on seafloor morphology worldwide. Upon evaluating the overlap between tectonic and Milankovitch periodicities across spreading rates, we reemphasize that fast-spreading ridges are the best potential recorders of a sea level signature in seafloor bathymetry.

Comment on "Sensitivity of seafloor bathymetry to climate-driven fluctuations in mid-ocean ridge magma supply"

Olive et al (Reports, 16 October 2015, p. 310) and Goff (Technical Comment, 4 September 2015, p. 1065) raise important concerns with respect to recent findings of Milankovitch cycles in seafloor bathymetry. However, their results inherently support that the Southern East Pacific Rise is the optimum place to look for such signals and, in fact, models match those observations quite closely.

Response to Comment on "Sensitivity of seafloor bathymetry to climate-driven fluctuations in mid-ocean ridge magma supply"

Huybers et al present new bathymetric spectra from an intermediate-spreading ridge as evidence for a primary contribution of sea level cycles to the morphology of the seafloor. Although we acknowledge the possibility that sea level-modulated magmatic constructions may be superimposed on a first-order tectonic fabric, we emphasize the difficulty of deciphering these different contributions in the frequency domain alone.