How sea level paces faulting at fast-spreading mid-ocean ridges

TL;DR

This study proposes a mechanism linking Pleistocene sea-level variability to abyssal hill fault spacing at fast-spreading mid-ocean ridges, supported by extended elastic unbending theory and numerical modeling.

Contribution

It introduces a quantitative model showing how sea-level-driven plate-thickness perturbations influence fault spacing, explaining the spectral peak near 41 ky.

Findings

Fault spacing correlates with Pleistocene sea-level cycles.

Small perturbations (~0.1%) in plate thickness can phase-lock faulting.

Model predictions match observed abyssal-hill spectral peaks.

Abstract

Abyssal hills, arguably the most extensive coherent pattern in Earth's surface topography, record the spacing of normal faults formed at mid-ocean ridges. At fast-spreading ridges, high-resolution bathymetry shows a pronounced spectral peak near 41 ky, coincident with obliquity-paced Pleistocene sea-level variability. The origin of this apparent orbital imprint on seafloor structure remains unresolved. We hypothesise that glacial-interglacial sea-level variability influences fault spacing by modulating plate thickness and the flexural stresses produced during plate unbending. Sea-level change alters mantle melting rates and magma supply at ridge axes, generating variations in the properties of the accreting plate. As the plate moves off axis, it unbends from its ingrown curvature, producing tensile fibre stresses that drive normal faulting. We hypothesise that small perturbations in elastic plate thickness modulate these stresses and thereby influence fault spacing. To test this, we extend the elastic unbending theory of Buck (2001) to include spatially variable plate thickness and yield-weakening viscoplastic flexure, which localises deformation into discrete kinks interpreted as faults. Linearised analysis shows that plate-thickness perturbations generate proportional fibre-stress variations. Numerical solutions demonstrate that perturbations as small as approximately 0.1 percent can phase-lock faulting to the imposed forcing. When driven by plate-thickness perturbations derived from the Pleistocene oxygen-isotope record, the model predicts fault spacings concentrated near 41 ky in the early Pleistocene and near 100 ky in the late Pleistocene, consistent with observed abyssal-hill spectra. These results provide a quantitative mechanism by which glacial-interglacial sea-level variability can be transmitted into tectonic structure.