Abstract:A key initial step in geophysical imaging is to devise an effective means of mapping the sensitivity of an observation to the model parameters, that is to compute its Fréchet derivatives or sensitivity kernel. In the absence of any simplifying assumptions and when faced with a large number of free parameters, the adjoint method can be an effective and efficient approach to calculating Fréchet derivatives and requires just two numerical simulations. In the Glacial Isostatic Adjustment problem, these consist of a forward simulation driven by changes in ice mass and an adjoint simulation driven by fictitious loads that are applied at the observation sites. The theoretical basis for this approach has seen considerable development over the last decade. Here, we present the final elements needed to image 3-D mantle viscosity using a dataset of palaeo sea-level observations. Developments include the calculation of viscosity Fréchet derivatives (i.e. sensitivity kernels) for relative sea-level observations, a modification to the numerical implementation of the forward and adjoint problem that permits application to 3-D viscosity structure, and a recalibration of initial sea level that ensures the forward simulation honours present-day topography. In the process of addressing these items, we build intuition concerning how absolute sea-level and relative sea-level observations sense Earth's viscosity structure and the physical processes involved. We discuss examples for potential observations located in the near field (Andenes, Norway), far field (Seychelles), and edge of the forebulge of the Laurentide ice sheet (Barbados). Examination of these kernels: (1) reveals why 1-D estimates of mantle viscosity from far-field relative sea-level observations can be biased; (2) hints at why an appropriate differential relative sea-level observation can provide a better constraint on local mantle viscosity and (3) demonstrates that sea-level observations have non-negligible 3-D sensitivity to deep mantle viscosity structure, which is counter to the intuition gained from 1-D radial viscosity Fréchet derivatives. Finally, we explore the influence of lateral variations in viscosity on relative sea-level observations in the Amundsen Sea Embayment and at Barbados. These predictions are based on a new global 3-D viscosity inference derived from the shear-wave speeds of GLAD-M25 and an inverse calibration scheme that ensures compatibility with certain fundamental geophysical observations. Use of the 3-D viscosity inference leads to: (1) generally greater complexity within the kernel; (2) an increase in sensitivity and presence of shorter length-scale features within lower viscosity regions; (3) a zeroing out of the sensitivity kernel within high-viscosity regions where elastic deformation dominates and (4) shifting of sensitivity at a given depth towards distal regions of weaker viscosity. The tools and intuition built here provide the necessary framework to explore inversions for 3-D mantle viscosity based on palaeo sea-level data.

Lateral and radial viscosity structure beneath Fennoscandia inferred from seismic and magnetotelluric observations

Thermal evolution of the North Atlantic lithosphere: New constraints from magnetic anomaly inversion with a fractal magnetization model

GIA imaging of 3-D mantle viscosity based on palaeo sea level observations - Part I: Sensitivity kernels for an Earth with laterally varying viscosity

A Crustal LVZ in Iceland Revealed by Ambient Noise Multimodal Surface Wave Tomography

Variation in bridgmanite grain size accounts for the mid-mantle viscosity jump

A Nearly Water-Saturated Mantle Transition Zone Inferred from Mineral Viscosity

Effect of Mantle Viscosity Structures on Simulations of Geoid Anomalies in the Ross Sea Area

Constraining mantle viscosity structure for a thermochemical mantle using the geoid observation

Vertical Displacements and Sea‐Level Changes in Eastern North America Driven by Glacial Isostatic Adjustment: An Ensemble Modeling Approach

Migration and accumulation of hydrous mantle incipient melt in the Earth's asthenosphere: Constraints from in-situ falling sphere viscometry measurements

Gravity, topography, and melt generation rates from simple 3D models of mantle convection

Composition Versus Temperature Induced Velocity Heterogeneities in a Pyrolitic Lower Mantle

A Common Diffusional Mechanism for Creep and Grain Growth in Polymineralic Rocks: Application to Lower Mantle Viscosity Estimates

Mantle Structure and Flow Beneath the Central-Western US: Constraints from Anisotropic Tomography

Timescales of glacial isostatic adjustment in Greenland: is transient rheology required?

Geodetic Observations of Time‐Variable Glacial Isostatic Adjustment in Southeast Alaska and Its Implications for Earth Rheology

Role of mush complex viscosity in modulating axial topography in mid-oceanic ridges

Effect of Rheology Structure on Flow Patterns Driven by Density Anomaly in the Mantle

Constraining Earth's nonlinear mantle viscosity using plate-boundary resolving global inversions

Present‐day Plate Motion Constraint on Mantle Rheology and Convection

Constraining models of glacial isostatic adjustment in eastern North America