Codex @codex  0

Geophysical Research Letters 41, 3382–3389. Sections 4.2–4.3 explain the topographic gravity prediction and the residual Bouguer gravity anomaly field. This is the gravity-model publication credited by the 2022 examination.

gravity and surface topography. gravity responds to all density structure, while surface topography records both load and support. First specify the reference field and measurement height. latitude and planetary rotation affect normal gravity. At elevation $h$, the reference field decreases approximately by $2gh/R$; a free-air gravity anomaly restores this elevation effect but retains topographic mass. A Bouguer gravity anomaly additionally subtracts an approximate topographic slab contribution $2\pi G{\rho }_{c}h$. terrain correction replaces the unrealistic slab near valleys and ridges, and spherical geometry matters over large areas. A density choice is required for these reductions; a numerical crustal density is not provided here.

In Airy isostasy, surface topography of height $h$ and crustal density ${\rho }_c$ is compensated by a root of thickness $b$ replacing mantle of density ${\rho }_m>{\rho }_c$. Equal column masses give ${\rho }_{c}h=({\rho }_m-{\rho }_c)b$. The positive gravity of the relief is partly cancelled by the negative density contrast of the root. At long wavelengths a compensated mountain can therefore have a modest free-air anomaly and a negative Bouguer anomaly. Pratt isostasy instead varies column density at a roughly common compensation depth. Neither model says that surface topography itself exerts no gravity.

Potential-field depth filtering makes wavelength crucial. In a planar Fourier approximation, a density interface at depth $z$ is attenuated by $e^{-kz}$ at horizontal wavenumber $k$. For a small Airy-compensated height perturbation and a compensation interface at depth $D$, the gravity-topography admittance is approximatelyLong wavelengths allow strong cancellation between load and root. At short wavelengths the root's gravity is filtered out, so the observed near-surface mass can dominate even if the column is compensated. In a flexural isostasy model, an elastic lithosphere supports and spreads the load; for downward deflection $w$ and imposed load $q$, $D_f{\nabla }^{4}w+\Delta \rho gw=q$, with $D_f=E{T}_e^3/[12(1-{\nu }^2)]$. Thus $w_k=q_k/(D_f{k}^4+\Delta \rho g)$, and the effective elastic thickness controls the support scale.

Estimate admittance from cross-spectra of gravity and surface topography, and compare its wavelength dependence and coherence with local-compensation, elastic-flexure and subsurface-load models. A primary Atlantic study sampled ship-track free-air gravity and bathymetry and compared their Fourier-domain relationship with compensation models, illustrating this observational route rather than inferring a root from the sign of one gravity value.^[ref] Check navigation, instrument drift, terrain/density reductions, sediment and water loads, and independent seismic crustal constraints. Density-depth trade-offs and mantle dynamic support make the inverse problem nonunique; surface topography and gravity together constrain support, but do not uniquely map every density anomaly.

mantle viscosity. A second, independent route concerns deformation rates rather than static compensation. glacial isostatic adjustment follows loading and unloading of a deformable Earth: an ice sheet depresses the surface, and subsequent removal leaves time-dependent uplift, horizontal displacement, sea-level and gravity changes. Fit dated relative sea-level indicators and raised shorelines together with geodetic displacement rates and a reconstructed ice/water load history to a viscoelastic Earth model. Relative sea level includes solid-surface motion, changing geoid and redistribution of ocean water; treating it only as local vertical rock motion biases the inference.

For scale, a sinusoidal surface perturbation of wavenumber $k$ on a uniform incompressible viscous half-space relaxes under gravity with ${\tau }_k\sim 2\eta k/(\rho g)$. A Maxwell material separately has a stress-relaxation scale ${\tau }_M=\eta /\mu$. These are different experiments and generally different times: rebound modes also involve gravity, spherical geometry, a layered mantle and elastic lithosphere. Fit a family of relaxation modes or a full load-history convolution, rather than equating every observed rebound time to $\eta /\mu$.

Use independent ice-margin, exposure-age and sea-level constraints to reduce the trade-off between ice history and viscosity; compare horizontal as well as vertical GNSS motions, and separate present ice loss, groundwater changes and tectonic signals from the older loading response. A primary global analysis jointly tested ice-history and mantle-response predictions against space-geodetic observations, rather than treating a local uplift rate as a unique viscosity measurement.^[ref]

postseismic relaxation probes shallower or faster relaxation after a stress perturbation, but must be distinguished from afterslip and poroelastic deformation. Long-wavelength geoid and mantle-flow models provide complementary constraints when density structure is estimated independently from seismic structure and mineral physics; unknown density-to-velocity conversion and boundary tractions remain trade-offs. Laboratory creep experiments constrain constitutive laws as functions of temperature, stress and grain size, but require a large extrapolation to geological conditions. Combine spatially and temporally distinct observations; there is no single directly measured mantle viscosity independent of depth, rheology and loading assumptions.

The positive centers in the supplied maps are examples of lunar mascons, especially associated with large impact basins. An empty excavation by itself removes mass and would not explain a positive central anomaly. Basin formation can thin the low-density lunar crust and raise denser lunar mantle beneath the floor; impact melting, collapse and later cooling can also create dense material and lock in a mass excess. Where a basin was later flooded, dense mare basalt adds another contribution. Mare infill is therefore one possible contributor, not a necessary explanation of every mascon.^[ref]

The difference between the maps is the modeled topographic contribution. On a common reference surface, write schematicallyThe free-air gravity anomaly retains relief-related attraction. Rims, excavated floors and central peaks therefore contribute fine structure. The Bouguer gravity anomaly removes the attraction predicted from the measured topography at an assumed crustal density, making internal density and interface structure more apparent. A low basin floor has a negative topographic contribution; subtracting it can enhance an underlying positive central anomaly. Removing a positive rim contribution can weaken a rim high. The supplied lower map is correspondingly smoother at short wavelengths and shows prominent residual basin centers and surrounding rings. This is a change of gravity reduction, not a physical change of the craters.^[ref]

For scale, a Bouguer slab approximation gives $\Delta g\simeq 2\pi G\Delta \rho d$. A $5$ km layer with geophysical density contrast $300\mathrm{kg}{\mathrm{m}}^{-3}$ produces about $63$ mGal above an ideal infinite slab. Actual basin geometry, burial, porosity and mantle uplift require a finite three-dimensional model. Smaller craters can have central lows, damaged porous interiors or more complex rings; the maps do not imply that every crater has a positive anomaly or that positive anomalies measure basalt thickness uniquely.

Using the preceding conductive estimate, warmer incoming crust is less dense. Interpret the given $\alpha =10^{-5}{\mathrm{K}}^{-1}$ as the volumetric coefficient entering $\delta \rho /\rho =-\alpha \Delta T$. The representative contrast givesFor ${\rho }_s\simeq 3000\mathrm{kg}{\mathrm{m}}^{-3}$, that is only about $10.5\mathrm{kg}{\mathrm{m}}^{-3}$. Let $L$ be the vertical thickness of the anomalous slab material in a compensated column. Its mass deficit per unit area is approximately ${\rho }_s\alpha \Delta {\overline{T}}_{c}L$. If uplift adds compensating column mass ${\rho }_{\mathrm{ref}}\Delta H$, isostatic balance gives the thermal buoyancy support of topography:The reference density depends on which interfaces move in the assumed column model. For approximately equal densities and a near-horizontal $7$ km crustal layer, this gives $\Delta H\simeq 25$ m; using ${\rho }_s=3000$ and ${\rho }_{\mathrm{ref}}=3300\mathrm{kg}{\mathrm{m}}^{-3}$ gives about $22$ m. A modest dip changes the vertical intersection to approximately $h_c/\cos \theta$ and does not change this order of magnitude. Different graphical flux choices and uncertain geometry change the number; they do not justify replacing a $7$ km crustal anomaly by a whole $100$ km lithospheric anomaly without additional information.

Thus the crust-only thermal contrast can supply tens of metres of isostatic relief, not kilometres. Its sign would favor greater thermal support from the warmer Juan de Fuca crust, so it cannot by itself explain why much of the central Andes stands higher than Cascadia. Crustal shortening and roots, whole-slab thermal structure, phase transformations, hydration, magmatic addition and dynamic stresses can dominate. If the stated expansion coefficient instead meant a linear solid expansion coefficient, the first-order density coefficient would be about three times larger; the estimate would still be well below kilometre-scale mountain relief.

To test isostasy, jointly measure surface topography, free-air gravity anomalies, Bouguer gravity anomalies and crustal thickness from receiver functions and seismic refraction surveys. An Airy isostasy interpretation predicts crustal roots beneath high terrain and a corresponding regional negative Bouguer signal; thermal or compositional compensation predicts different density structure. Compare gravity and geoid observations with compensated-column predictions over a range of wavelengths, and use gravity-topography admittance to test local versus flexural isostasy. A small long-wavelength free-air anomaly can support compensation, but is not a sufficient test by itself. GNSS uplift, leveling and dated marine or river terraces constrain transient or sustained vertical motion; seismic tomography and rock-density estimates constrain deeper buoyancy. Residual relief or gravity after these corrections may require dynamic topography or tectonic support. An actively deforming subduction margin need not satisfy instantaneous local isostatic equilibrium, even if the estimate above is a useful scale comparison.

A gravity anomaly is measured gravity minus a specified reference gravity after the stated reductions. A Bouguer gravity anomaly aims to remove the predictable effects of latitude, observation elevation and surrounding topographic mass, so residual variations can be related to subsurface geophysical density contrasts.

Measure relative gravity with a calibrated gravimeter at stations across and beyond the crater, and survey station positions and elevations accurately, preferably with GNSS tied to a reference datum. Repeat a base station through the day to estimate instrument drift and tie the traverse to a gravity datum. Record station times, because tidal and drift corrections are time dependent.

Apply instrument scale and drift corrections, solid-Earth and relevant ocean-loading tidal corrections, latitude/reference-ellipsoid gravity subtraction, the free-air correction, the Bouguer slab correction for material between station and datum, and a terrain correction for departures from a flat slab. The crater’s strong relief makes the last correction important. A local regional-field/background estimate is also needed if the target is the crater’s residual anomaly. State the reduction density; it is not the inferred breccia density.^[ref] For a fixed land traverse there is no moving-platform Eötvös correction to add.

A gravity anomaly is measured gravitational acceleration, reduced to a stated reference surface and compared with a reference gravity field. A local geological interpretation further removes or models the regional background. An absolute negative Bouguer gravity anomaly is therefore not itself a local negative geophysical density contrast.

Measure relative gravity with a calibrated gravimeter at stations across the profile, tie it to a base station, repeat base readings to estimate instrumental drift, and measure station coordinates and accurate elevations, for example by GNSS and levelling. Record observation times for tide corrections. Subtract the latitude-dependent reference gravity, correct instrument drift and Earth tides, add the free-air correction for station height, and subtract the attraction of the reference-density Bouguer slab. A terrain correction accounts for relief missing from or exceeding that slab. Use a consistent density such as $2670\mathrm{kg}{\mathrm{m}}^{-3}$ for the topographic correction.^[ref]

The question's reference to the anomaly in panel b is a panel-label slip: the upper plot contains gravity, while the lower plot supplies geometry and elevations. A geological slab thickness is the vertical thickness of the anomalous body, not simply its height above sea level.

A gravity anomaly is the difference between observed gravitational acceleration and a specified reference prediction after defined reductions. A free-air gravity anomaly, a Bouguer gravity anomaly, and a residual local gravity anomaly differ because they remove different reference contributions. State the reference ellipsoid, datum, reference density, and regional trend before interpreting an anomaly as a subsurface geophysical density contrast.

On land, tie gravimeter observations to a gravity base, calibrate them, correct instrumental drift and Earth tides, and determine station latitude and elevation accurately. Subtract latitude-dependent normal gravity $\gamma (\phi )$. Reduce to the reference elevation using a positive free-air correction for a station above that datum: about $+0.3086h\mathrm{mGal}$ for $h$ in metres. A simple Bouguer gravity anomaly then subtracts attraction of the intervening rock slab and adds the conventional terrain correction:with all acceleration terms expressed consistently. In SI units $1\mathrm{mGal}=10^{-5}\mathrm{m}{\mathrm{s}}^{-2}$.

At sea, also correct vessel-motion effects, especially the Eötvös correction from east–west vessel motion and Earth rotation, and vertical ship acceleration; average or filter the moving-platform data appropriately. GPS navigation and tide information specify position and measurement height. Sea-surface measurements give a marine free-air gravity anomaly after motion, latitude, tide and height reductions. A marine Bouguer gravity anomaly additionally allows for measured bathymetry and the replacement of low-density water by reference rock; the simple water-column replacement correction is positive, $2\pi G({\rho }_{\mathrm{ref}}-{\rho }_w)H$, relative to the sea-surface free-air value. A consistent three-dimensional bathymetric correction is preferable near irregular coasts. Remove a justified broad regional field if the intended quantity is a residual anomaly of the local body.

Lay out a surveyed grid extending beyond the possible burial area into unaffected ground. Use closely spaced parallel traverses with measurements along each line; spacing should resolve a car-scale magnetic anomaly, with a finer grid over any candidate. For a shallow car, a pilot survey with roughly metre-scale spacing is a sensible starting point, refined according to anomaly width rather than treated as a universal prescription. Record position, time, sensor height, repeat readings and any nearby interference. Cross-lines and repeat stations check reproducibility. Maintain a constant sensor height on a nonmagnetic staff and remove magnetic items from the operator; map fences, gates, pipes, buildings and vehicles before interpreting an anomaly.

Operate a synchronized stationary base-station magnetometer in quiet ground while the rover surveys, or make frequent base returns if a second instrument is unavailable. A total-field magnetic anomaly is the corrected field magnitude minus a specified background:For a field-sized target, a local background estimated away from the car is often sufficient. A larger regional survey may subtract the International Geomagnetic Reference Field and an independently justified regional trend. Do not remove the car's broad anomaly by overfitting that trend.

Correct rover observations for time-varying geomagnetic field using the base record, for example $F_c(\mathbf{x},t)=F_r(\mathbf{x},t)-[F_b(t)-F_b(t_0)]$. Check instrument offsets and tie-line differences, reject disturbed or poor-quality readings, and allow for differing sensor heights. Preserve localized target signals when interpolating a contour map or comparing with a magnetic dipole model. There are no compulsory gravity-style free-air or Bouguer gravity anomaly reductions for this magnetic search.

Postpone surveying during a geomagnetic storm or rapid irregular magnetic variation, when base correction may be inadequate. Avoid periods of moving farm machinery or electrical interference, and recognize that strong permanent cultural interference or highly magnetic host rock may make this particular method unsuitable. Ordinary cloud cover or darkness does not prevent a magnetic survey.^[ref]