2026-05-12 | | Total: 4
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.
Travel-time tomography forces a trade-off between mesh resolution and stability in which the regularizer choice dominates what can be recovered. We introduce MIMIR, a differentiable framework that represents the 2D velocity field as a Fourier-feature neural network, replacing the grid-based slowness vector with a continuous, infinitely differentiable function. Prior neural-field tomography has staircased smooth fields under total-variation (TV) priors or oscillated near interfaces under $L^2$ Laplacian smoothing. We adopt second-order total generalized variation (TGV$^2$) and parametrize its auxiliary vector field as a second neural network jointly optimized with the velocity field, eliminating the inner Chambolle-Pock primal-dual loop that classically dominates TGV computation. On three synthetic benchmarks (Gaussian, horizontally layered, curved-fault inspired by OpenFWI) using cross-well acquisition, 5% travel-time noise, and five seeds, MIMIR-TGV$^2$ ties a classical FMM-LSMR baseline with auto-tuned hyperparameters on the Gaussian ($p=0.134$, paired $t$-test) and significantly outperforms it on layered ($p<0.0001$, 44% RMSE reduction) and curved-fault ($p=0.0002$, 33% reduction). Replacing TGV$^2$ with TV degrades performance on Gaussian ($p=0.004$) and layered ($p=0.003$); curriculum-annealed TV improves Gaussian RMSE by only 5.4%, confirming that TV's staircase bias is intrinsic to the regularizer rather than a scheduling artifact. The results empirically validate the Bredies-Kunisch-Pock prediction that piecewise-affine priors are better suited to subsurface velocity recovery than piecewise-constant TV priors. We argue that the central design choice in physics-informed neural-field inversion is not the network architecture but the regularizer. The full pipeline reproduces in under one hour on consumer hardware.
Full waveform inversion (FWI) can produce accurate subsurface velocity models. However, the lack of sufficiently low-frequency content in field data often causes cycle skipping and traps the inversion in local minima. The Hilbert-transform envelope (HTE) provides a low-frequency representation that helps mitigate cycle skipping, but it may be insufficient when the initial velocity model is highly inaccurate. To further enhance low-frequency information and reduce dependence on the initial model, we compute an approximate envelope using a sequence of 2D max-pooling operations. Compared with HTE, the resulting max-pooling-based approximate envelope (MPBAE) contains richer low-frequency components and better mitigates cycle skipping. We further combine the MPBAE loss with a shot patching strategy and exploit the inherent normalization property of the Euclidean loss to formulate the MPBAEP loss, in which each shot gather is divided into localized patches for misfit evaluation. This introduces local adjoint-source energy balancing, as the adjoint source associated with the Euclidean loss exhibits a normalization effect within each local region, thereby improving gradient balance and accelerating convergence. Numerical experiments on synthetic and field data demonstrate that MPBAE-FWI significantly outperforms HTE-FWI when the initial model is poor, while MPBAEP-FWI further improves inversion accuracy.
Electromagnetic (EM) sounding can constrain the electrical structure of Enceladus and, in turn, the salinity of its ocean and the porosity, fluid content, and thermal state of its hydrothermally active core. Here, we assess the feasibility of EM sounding at Enceladus using both global (orbiter) and local (lander) EM induction transfer functions. We provide a physical framework for modeling EM induction for 1-D and 3-D subsurface conductivity models and discuss how transfer functions can be estimated from global or local measurements of the magnetic and electric fields. We simulate 3-D induction effects arising from variations in ice-shell thickness. The magnitude of these effects in the magnetic field correlates with the ice-shell thickness at the surface and is strongly dependent on the ocean's conductivity. These magnetic variations, if observed, would favor a moderately to highly conductive ocean, providing lower bounds on salinity and volatile content. The absence of these effects indicates a thicker, more homogeneous ice shell and/or a lower-conductivity ocean. Given plausible magnitudes, a polar-orbiting mission with low-altitude measurements will be required to detect these effects. In summary, an orbiter will constrain global ocean conductivity using long-period induction and possibly map the ice thickness variations. The detailed EM sounding of both the hydrosphere and the core can be achieved by a lander-based broadband EM sounding at periods $\approx 10^1-10^5$ s to probe ocean salinity and thickness, as well as core properties including porosity, fluid content, and temperature.