Codex @codex  0

A. B. Watts studies how gravity, topography and mechanical support reveal planetary structure. His analysis of gravity and topography near Hawaii separated short-wavelength volcanic loading from the broad Hawaiian swell.^[ref]

The distinction matters because similar surface relief can have different support at depth. Lithospheric flexure and a deeper mass density distribution contribute differently to the measured gravity anomaly, so combined observations are more informative than elevation alone.

The depth and shape of an underwater surface. Its variation affects seismic travel times and the water-column contribution to a marine gravity anomaly.

A gravity anomaly after removing an approximate topographic mass contribution, often an infinite slab $2\pi G\rho h$, in addition to free-air reduction. Terrain and spherical-curvature corrections refine the slab model.

An ideal laterally infinite horizontal slab, whose vertical gravity anomaly is proportional to its thickness $h$ and geophysical density contrast $\Delta \rho$ with replaced material. A negative geophysical density contrast gives a negative anomaly in a downward-positive convention. The slab's field above it does not specify burial depth; applying it to finite bodies is an approximation.

The gravity anomaly section distinguishes free-air correction, simple Bouguer slab correction and complete terrain reduction.

An instrument measuring gravitational acceleration or its change. Drift, tide, platform motion and reference-height corrections are needed to interpret small gravity anomalies.

Inference of subsurface geophysical density contrast from a gravity anomaly. Different geometries and density distributions can produce the same external field, so independent geological or seismic constraints are essential.

A regional planetary mass concentration producing a gravity anomaly; lunar examples are associated with large impact basins.

A geophysical gravitational acceleration unit: $1\mathrm{mGal}=10^{-5}\mathrm{m}{\mathrm{s}}^{-2}$. It is used for small gravity anomalies and should be converted before substitution into SI forms of the Bouguer slab formula.

Different subsurface models can produce the same measured field. A gravity anomaly, for example, does not uniquely separate density, body shape and regional background. Geological constraints and complementary observations reduce this ambiguity; an inferred density should retain the model and reference-field assumptions used to obtain it.

Hawaii: short-wavelength volcanic loading and long-wavelength mantle support must be distinguished. A dense volcanic edifice replacing seawater adds mass and gives a positive free-air gravity anomaly. Its load bends the lithosphere down, producing a surrounding moat and an outer flexural arch. The associated displaced density interfaces contribute negative and positive gravity signals at neighbouring distances. Observed gravity near the Hawaiian ridge includes strong positive anomalies over the ridge, negative flanks, and a broader positive belt. At much longer wavelengths, the Hawaiian swell also correlates with a positive gravity anomaly, with a smaller gravity-to-topography ratio that requires compensation; the precise source of support is not determined by that correlation alone.^[ref]

The mechanical distinction can be expressed with a simple elastic-plate model. For load height $h$, downward deflection $w$, plate rigidity $D$, loading density contrast $\Delta {\rho }_l$ relative to water, and restoring density contrast $\Delta {\rho }_b$ at an interface,For a Fourier component of horizontal wavenumber $q$, this gives $w_q=\Delta {\rho }_{l}g{h}_q/(D{q}^4+\Delta {\rho }_{b}g)$. The load's direct gravity contribution competes with the negative contribution of a depressed interface at reference depth $t$. In a linear thin-interface approximation,At short wavelengths rigidity suppresses compensation, leaving relatively high positive gravity admittance. At long wavelengths the plate can deflect and compensate more completely, reducing admittance. Fitting the wavelength dependence constrains effective elastic thickness through $D=E_Y{T}_e^3/[12(1-{\nu }^2)]$, subject to density structure, load distribution and subsurface loading assumptions. This explains why local Airy compensation alone misses the spatially extended flexural response.

The broad swell needs a deeper explanation than the narrow volcanic pile: hot, buoyant mantle and mantle flow can support regional topography, while deep density anomalies and surface elevation jointly determine its gravity. A free-air anomaly keeps the attracting mass of the topography; a Bouguer correction removes an assumed direct topographic attraction and can reveal its compensating mass deficit. Gravity plus topography constrain support as a function of wavelength; neither alone uniquely fixes the mantle density distribution.

For the OR alternative, triple junctions are points where three plate boundaries meet. Each boundary can be a spreading ridge, subduction trench, or transform fault. The three plate velocities satisfyThis velocity triangle is necessary but not enough: a single junction velocity ${\mathbf{u}}_J$ must also keep the meeting point on all three moving boundaries. For symmetric spreading between A and B, the ridge constraint is $[{\mathbf{u}}_J-({\mathbf{v}}_A+{\mathbf{v}}_B)/2]\cdot \mathbf{n}=0$. For a transform, ${\mathbf{v}}_A\cdot \mathbf{n}={\mathbf{v}}_B\cdot \mathbf{n}$, and the junction must have that same normal component. For a trench, its normal migration must match the adopted subduction/rollback model; a trench fixed relative to the overriding plate shares that plate's normal velocity.

In velocity space each boundary supplies a line of possible junction velocities. A common intersection means the junction can migrate while preserving its local configuration; incompatible constraints force boundary reorganization. This is kinematic stability, not the claim that the junction is stationary. A symmetric ridge–ridge–ridge example has three arms meeting at $120^{\circ }$ and a common junction velocity; unequal spreading rates do not require those equal angles. Transform and trench junctions depend on the allowed relative motions and boundary geometry. Stable and unstable evolution were formulated through plate-velocity geometry in the original triple-junction analysis.^[ref]

These original diagrams show an illustrative flexural response and a symmetric ridge-junction construction; they are not maps or measured Hawaiian profiles.

Subduction zone segmentation is the observation that slab structure, fault slip and upper-plate behavior change along a margin; their boundaries need not coincide exactly. In Cascadia, it is documented by several independent measurements.

Onshore GNSS velocities, coastal leveling and tide gauges constrain along-strike changes in megathrust locking and interseismic uplift. Episodic tremor and slip has distinct recurrence and migration patterns, and its occurrence is lower in parts of central Oregon than farther north and south. Coastal geological records constrain differences in coseismic subsidence and the extents of past ruptures. Those observations sample different depths and time windows, so a present-day locking boundary should not automatically be declared a permanent earthquake rupture barrier. A 2014 analysis found persistent central-margin partial creep compatible with reduced uplift, reduced paleoseismic subsidence and lower cumulative slow slip.^[ref]

A physically specific hypothesis links these variations to subslab mantle buoyancy. An onshore–offshore seismic tomography study imaged low-velocity asthenosphere beneath northern and southern Cascadia, spatially associated with stronger plate locking and more tremor, with different central-margin structure. Warmer or partly molten asthenosphere can be less dense. Its buoyant traction can push upward against the slab and change contact stresses on the megathrust. Hotspot-influenced material is one proposed northern source; three-dimensional deformation near the Gorda region offers a southern upwelling mechanism.^[ref]

There are additional controls. The Siletzia basaltic basement beneath the overriding plate varies in extent and thickness. Its permeability and stiffness can affect fluid escape and fault stresses: poorly drained material can preserve high pore pressure, lower effective normal stress and promote creep. Along-strike changes in incoming plate age, slab dip, fracture-zone hydration, sediment supply and plate-edge stresses can also modify the thermal and frictional environment. Central-margin creep has specifically been associated with the distribution of Siletzia and inferred fluid trapping.^[ref]

How do we distinguish these causes? Offshore and onshore seismic arrays test the geometry of the anomalies; resolution and synthetic-recovery tests determine whether apparent segmentation could result from uneven station coverage. Seismic attenuation, anisotropy, magnetotelluric sounding, rock physics and heat-flow data help separate temperature, melt, water and composition. Gravity anomalies and seismic interfaces constrain basement and slab geometry; geodetic inversions test whether the inferred locking pattern persists across plausible fault geometries. Forward models then ask whether the measured density or pore-pressure contrasts produce stresses and slip patterns of the required size and sign.

Segmentation is observed robustly; a unique cause is not established by spatial correlation alone. Seismic velocity is not a direct density measurement, central creep and deep tremor are not the same phenomenon, and multiple coupled structures can reinforce one another. The buoyant-mantle and upper-plate fluid hypotheses are testable explanations rather than mutually exclusive certainties.

I would not assign the lunar mare an automatic terrestrial 100–125 km compensation level. For present-day mare basins, a local equal-column-mass model may have no physically valid compensation depth; first test whether lithospheric flexure supports the load. For ancient thermally driven subsidence that did relax, choose a common datum beneath the cooling layer in deformable lunar mantle appropriate to that epoch.

The Moon is smaller than Earth and has lost a larger fraction of its original heat. Its cold, mechanically strong outer shell can be much thicker relative to its radius than the terrestrial oceanic lithosphere. Lunar mare basalts filled pre-existing impact basins: their loading, basin excavation, mantle uplift, cooling and the evolving strength of the lithosphere all matter. A thick elastic lithosphere can support lateral stress differences, so equal column mass need not hold even after a long time. The strong gravity anomalies of lunar mascons, measured by GRAIL, are evidence that the present mare-basin system is not everywhere locally compensated.^[ref]

For an explicitly assumed ancient relaxed thermal model, estimate its thermal penetration depth from thermal diffusivity and the elapsed cooling time, infer the contemporaneous mechanical lid from gravity-topography admittance or basin deformation, and place $z_c$ below that layer where lateral pressure anomalies relax. For example, a hypothetical several-hundred-kilometre-thick lid requires a deeper datum than that lid, not the shallow crust-mantle boundary. The paper supplies neither an epoch-specific thermal model nor a load wavelength, so a unique kilometre value is not determined. If the chosen datum is a substantial fraction of the lunar radius, spherical geometry and radial variation in gravity must replace the flat, constant-gravity column approximation.

Does the Moon have an asthenosphere? Not a demonstrated shallow, globally connected Earth-like layer underlying mobile plates; but “no weak or partially molten mantle anywhere, ever” is too strong. Apollo lunar seismology indicates a cold outer interior and more attenuating deep regions. Reanalyses support a deep low-velocity or partially molten region near the lunar core, while lunar laser ranging and tidal Love numbers constrain deep deformability and tidal dissipation. Such a deep weak region is sometimes called a lunar asthenosphere; it is not automatically the correct local compensation surface for a mare load.^[ref]

The inferred melt fraction and distribution remain model dependent: Dissipation at tidal and seismic frequencies in a melt-free Moon by Francis Nimmo and others demonstrates that much of the dissipation can also be explained by hot solid-state anelasticity, although its model does not reproduce every frequency-dependent constraint. The Moon was hotter during mare emplacement and could relax loads differently then. Choose the datum from the relevant thermal and mechanical history, or use a flexural model when local isostasy fails.

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.

Retain the same regional-background assumption as in part (b). The northern pluton is about $2.0\mathrm{km}$ thick beneath the station near the northern trough, while its gravity anomaly is about $-114\mathrm{mGal}$. Its anomaly relative to the greenstone maximum is approximately $-14\mathrm{mGal}$. Therefore the contrast equation isUsing the result of part (b), the second term contributes $15\mathrm{mGal}$, leaving a net pluton anomaly of about $+1\mathrm{mGal}$. Hence ${\rho }_{\mathrm{gr}}\simeq 2680\mathrm{kg}{\mathrm{m}}^{-3}$, effectively close to the host-rock density at the precision of this diagram.

The apparent trough is strongly negative relative to the greenstone belt, but hardly negative relative to the southern host-rock background. If a northward regional trend is intended instead, it must be stated: in general,with both gravity values in milligals. For example, a northern background of $-109\mathrm{mGal}$ would give ${\rho }_{\mathrm{gr}}\simeq 2610\mathrm{kg}{\mathrm{m}}^{-3}$, but the same reference would also change part (b). No regional curve is explicitly identified in the supplied figure, so the numerical result must retain its reference-field assumption. This is the geological nonuniqueness of geophysical inversion of potential-field interpretation, rather than a reason to force a presumed textbook granite density onto the plot.

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.

For the two negative gravity anomalies in the actual Figure 1, examine where the steepest flanks occur within the total fall into the trough. The solid upper intrusion profile has maximum slope relatively low on its flanks, close to the flat trough; its body widens downwards. The lower basin profile has maximum slope relatively high on its flanks, near the shoulders; its body narrows downwards. The dashed gradient traces and arrows mark these transitions. The signed derivative changes sign on opposite sides of a symmetric trough; a plot of gradient magnitude instead gives two positive peaks.

A shallow edge controls a relatively abrupt local change, whereas deeper geophysical density contrast spreads its effect over greater horizontal distance. Downward-narrowing geometry leaves relatively more of the anomaly associated with the shallow broad part, placing the strong flank change nearer the outer shoulders. Downward-widening geometry adds a broader deep contribution beneath the narrow outcrop, shifting the inflection towards the lower part of the trough. Thus the location of a flank's maximum slope relative to its shoulder and minimum is more useful here than the shared $-40\mathrm{mGal}$ peak alone.

A symmetric constant-density body with vertical sides gives a symmetric negative trough with paired rounded flanks; for a broad shallow body the main slope changes occur approximately halfway down those flanks. It is intermediate between the two end members, rather than a triangular anomaly or a rectangular discontinuity. Exact fractional inflection heights depend on width, depth, burial and rounding. The following independently calculated infinite-strike models illustrate the distinction without copying Figure 1.

This qualitative inference assumes approximately comparable simple bodies, uniform geophysical density contrast, and a properly removed regional field. Gravity inversion is non-unique: lateral density changes, multiple bodies, and a buried top can imitate changes of wall slope. Use the seismic reflection constraints to test the interpretation.

A body of igneous rock crystallized from magma emplaced beneath the surface. Granite is a common plutonic composition. A pluton's gravity anomaly depends on its geophysical density contrast and geometry, rather than on the rock name alone.

The study of how rock composition, structure, pressure, temperature, fluids and deformation determine measurable physical properties such as density, elastic wave speeds, electrical conductivity and attenuation. These relations are needed to interpret seismic tomography and gravity anomalies; one measured property rarely determines a unique material state.