Codex @codex  0

Marcos A. Pimenta coauthored the 2015 Raman investigation of isotopic composition in graphene.^[ref]

The harmonic oscillator supplies a useful reference model. With the same effective bond stiffness $k$ and an effective vibrating mass $m$, the mode frequency is proportional to $\sqrt{k/m}$. Replacing carbon-12 by carbon-13 changes mass without changing the number of protons; comparing isotopes can therefore help separate mass effects from changes in the bonding model.

For an original unchanged-stiffness estimate, a mode at $1580\mathrm{c}{\mathrm{m}}^{-1}$ in ideal carbon-12 material would move toin the corresponding ideal carbon-13 model. The approximately $3.92\%$ frequency reduction is smaller than the $8.33\%$ mass increase because frequency depends on the inverse square root of mass.

If stiffness changes too, the ratio instead contains both effects:A frequency shift alone then cannot identify a stiffness change without composition information. This connects Raman spectroscopy to an inverse problem: a measured peak position constrains a combination of model parameters. The numerical frequencies above are an illustrative calculation, not a reproduction of Pimenta's experimental data.

Martin Benning coauthored an account of modern regularization of an inverse problem for inverse problems with Martin Burger. It explains how prior information and stabilizing penalties make ill-posed inversion tractable, including nonlinear regularization for linear forward maps. Existence depends on conditions such as compactness and lower semicontinuity, not merely on writing down a minimization problem.^[ref]

Martin Burger coauthored an exposition of regularization of an inverse problem with Martin Benning. Their treatment of inverse problems connects nonlinear penalties, prior information and stability of reconstruction. The framework distinguishes exact inversion from a controlled approximation that remains meaningful when measurements contain noise or the inverse map is discontinuous.^[ref]

Matthias Morzfeld is a researcher in computational and mathematical geophysics whose work combines models with observations. His Scripps Oceanography profile identifies research on data assimilation, inverse problems, Monte Carlo methods and uncertainty, and records a mechanical-engineering doctorate at UC Berkeley in2011. His author-maintained UC San Diego homepage describes a particular interest in how sampling methods behave in high-dimensional spaces.^[ref]

His earlier work with Fai Ma and Ali Imam studied decoupling of viscously damped linear systems in oscillatory free vibration. Their2009 paper constructs a time-varying coordinate transformation by synchronizing phases of damped modes, extending classical modal analysis beyond its usual classical-damping setting.^[ref] This addresses how a damping matrix can keep motion coupled even when undamped normal modes have been identified.

An original Bayesian inference example illustrates the model–observation combination underlying his later research interests. Suppose a scalar state has prior distribution $X\sim \mathcal{N}(0,4)$ and an independent measurement error has a Gaussian distribution $\epsilon \sim \mathcal{N}(0,1)$, with observed $Y=X+\epsilon =3$. The posterior distribution has density proportional toThus $X\mid Y=3\sim \mathcal{N}(2.4,0.8)$. The estimate moves toward the measurement while retaining the information and uncertainty from the prior. The posterior variance is lower than either the prior variance4 or the measurement-error variance1 in this stipulated independent Gaussian model.

These chosen numbers are not a result from Morzfeld's papers. They isolate a simple update mechanism; nonlinear dynamics, sparse observations and large state spaces can make the corresponding computation much harder. His work develops and studies numerical tools for such model–data problems rather than assuming that a small Gaussian calculation suffices for an Earth-science application.^[ref]

The author-maintained homepage provides a primary biographical reference and links a profile image.

Michael D. Fayer is a physical chemist studying molecular dynamics with ultrafast optical methods. His Stanford-hosted curriculum vitae explicitly gives this name; the D. is retained without further expansion. His research description connects two-dimensional infrared spectroscopy with molecular interactions in liquids, confined water and other structured materials.^[ref]

With Kyungwon Kwak, Sungnam Park and Ilya J. Finkelstein, Michael D. Fayer coauthored the 2007 center-line method for extracting dynamical information from two-dimensional infrared spectroscopy. This is a shared paper contribution, not evidence for separate individual experimental roles.^[ref]

The scientific question can be illustrated by an original frequency-fluctuation model. Let $\delta \omega (t)$ be the deviation of a molecular vibrational frequency from its mean, and suppose its vibrational frequency-frequency correlation function isThe normalized memory is $C(T)/C(0)=e^{-T/{\tau }_c}$. With stipulated ${\tau }_c=3\mathrm{ps}$, it falls to about 0.368 after 3 ps and 0.0498 after 9 ps. An ordinary mean frequency would not reveal that memory time. These are model calculations, not measurements of confined water.

In the vibrational photon echo approach, controlled optical pulse delays connect an initial frequency label with a later response. The 2D IR center-line slope helps summarize changing spectral geometry. Interpreting that observable as molecular memory requires the response assumptions and axis convention; a slope fitted to arbitrary data is not automatically $C(T)/C(0)$.

That distinction connects apparatus with theory. Pulse timing, the infrared spectroscopy response and a model of environmental fluctuations all enter the inverse problem. Measurements can constrain how a molecular environment changes while leaving multiple microscopic explanations possible. Michael D. Fayer's Stanford research description explicitly places theoretical analysis alongside the experiments.^[ref]

Missing-section ambiguity in volume reconstruction is nonuniqueness caused by unobserved material between acquired sections. Even perfect image registration cannot determine every structure inside that material. It is an inverse problem: different three-dimensional specimens can produce exactly the same retained image stack. The serial-section acquisition sequence makes the missing observation explicit: cut, collect and image sections, then align and stack the surviving observations.^[ref]

For an original exact counterexample, consider a binary column of three adjacent voxels, ordered along the cutting direction. A value 1 denotes a marked material and 0 denotes its absence. Two possible columns areIf the middle section is missing, both acquired records are $(1,1)$. Yet $A$ is one connected column and $B$ consists of two separated parts. The acquisition map $P$ that retains only the first and last entries satisfiesThis is an ambiguity in the observations, not a failure of a particular reconstruction program. Linear interpolation inserts the middle value 1 and selects $A$, but that choice uses a smoothness assumption rather than an additional measurement.

With expected section spacing $40\mathrm{nm}$, losing the section centered at $z=40\mathrm{nm}$ leaves acquired centers at $0$ and $80\mathrm{nm}$. Keep their physical positions: silently relabeling them as adjacent $40\mathrm{nm}$ observations also distorts depth. Additional physical observations or justified constraints may reduce the possible reconstructions. A visually smooth surface alone cannot distinguish the two columns above.

Recovering molecular geometry from rotational constants is an inverse problem. In the synthetic assignment, $B=5\times 10^9\mathrm{Hz}$ gives $I=h/(8{\pi }^{2}B)\simeq 1.6784\times 10^{-45}\mathrm{kg}{\mathrm{m}}^2$. For a diatomic molecule with known masses, rotational constant explains how to infer an effective bond length. For a longer linear molecule, one $I$ cannot determine every bond length: several different arrangements can have the same moment of inertia.

For a fixed three-atom linear model with positions $0,a,a+b$ and positive masses $m_1,m_2,m_3$, subtracting the centre of mass givesThis is one equation for two unknown bond lengths. Even fixing the order of the chemical atoms leaves a family of solutions. Additional isotope substitutions change known masses and supply independent equations, provided the effective molecular geometry can be treated as shared. For nonlinear molecules, three principal moments of inertia give frequency-valued constants $A=h/(8{\pi }^2{I}_a)$, $B=h/(8{\pi }^2{I}_b)$, and $C=h/(8{\pi }^2{I}_c)$; a single evenly spaced comb is generally insufficient.

Vibrational averaging, unresolved hyperfine structure, internal motion, and centrifugal distortion limit structural inference. A ground-vibrational fitted rotational constant is not automatically the equilibrium constant. Ordinary field-free rotational spectroscopy also gives identical spectra for opposite enantiomers under the usual parity-symmetric molecular Hamiltonian; assigning composition and conformation does not by itself determine handedness. Molecular identification from a rich fingerprint can therefore succeed even when a unique full molecular geometry cannot be reconstructed.

Non-identifiability of shape from height and mass means that different external geometries can produce the same measured height and mass. The observation map from a body shape to those two numbers is many-to-one. A formula using only those inputs cannot distinguish shapes that the observations themselves do not distinguish. This is an identifiability issue in an inverse problem, independently of measurement precision.

An original family of homogeneous cuboids supplies an explicit proof. Fix height $H$, mass $M$ and density $\rho$, and define horizontal area $A=M/(\rho H)$. LetEvery real $t$ gives $ab=A$, the same height and the same mass. Its external surface area, however, isThe square section $t=0$ minimizes this family's area. At $t=\log 2$, the sides have ratio four and the lateral surface area is 25% larger. Height and mass remain exactly unchanged. This is a geometric construction, not a proposed anatomical model of a person.

With stipulated $H=1.8$ m, $M=80$ kg and $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$, those two total areas are approximately $1.607$ and $1.986{\mathrm{m}}^2$. The Mosteller formula assigns both the same predicted body surface area, $2.00{\mathrm{m}}^2$, because its documented inputs are only height and mass.^[ref] No reduction of numerical rounding can make one identical prediction equal both different geometric areas.

Additional measurements and anatomical identifiability can change the problem. In this particular cuboid family, measuring the horizontal perimeter $P=2(a+b)$ identifies the area through $S=2A+HP$. That conclusion relies on the assumed geometry; a perimeter does not generally determine the full external surface of an arbitrary human body.

Validation must therefore ask whether the observations and model contain the shape information needed for the intended quantity. Repeated precise height and mass measurements constrain the inputs. Independent geometric measurements or a validated restriction on allowable shapes constrain the missing anatomical information.

A spatial charge density describing the electric charge inside an atomic nucleus. In a static description centered on the nucleus, normalize it bywhere $Z$ is the atomic number and $e>0$ the elementary charge. A nuclear charge distribution is an average over a nuclear quantum state, rather than a photograph of immobile protons. Even when the average has spherical symmetry, individual nucleons need not sit on spherical shells. The proton charge radius and the internal charge density of the neutron also contribute: a distribution of pointlike proton centers is an approximation to a nuclear charge distribution, not its definition.

The experimental attraction is that an electron probes electric charge through the electromagnetic interaction, without itself feeling the nuclear force. Varying its momentum transfer changes the spatial scale tested. Long wavelengths mainly constrain a nuclear charge radius; shorter wavelengths can distinguish a diffuse edge from a sharper one through diffraction. Recovering the entire nuclear charge distribution from a finite range of noisy scattering data remains an inverse problem.

The optical contribution of the nail bed is the portion of measured light associated with tissue beneath the nail plate, after accounting for illumination, the covering plate and collection geometry. The vascular tissue beneath the plate contributes to its ordinary pink appearance, while the compact plate itself consists of dead keratinocytes.^[ref] A colour measurement combines optical transport and anatomy rather than directly reporting a single blood quantity.

An original simplified reflectance model iswhere $P$ represents light returned by the plate and interfaces, $T$ is single-pass plate transmittance and $R$ is the fraction returned from the underlying region. The square represents outward and return passage through the plate. This deliberately omits detailed multiple scattering; it is an explanatory model, not a clinical calibration.

Even with known $P=0$, $y=0.32$ can arise from $(T,R)=(0.8,0.5)$ or $(\sqrt{0.4},0.8)$. One measurement cannot distinguish those pairs. For small changes with $P=0$,Thus a 1% relative change in plate transmittance contributes approximately a 2% relative signal change without any change in underlying reflectance. This is a concrete inverse problem.

A primary experiment separated one optical component physically: thirteen autopsy toenails were measured at fifteen wavelengths between 300 and 600 nm using a monochromator, integrating sphere and photodiode, with a correction for nail fluorescence.^[ref] Measuring an isolated plate constrains its transmission; it does not measure living nail-bed blood flow.

Plate-bed optical identifiability requires additional constraints: independent plate measurements, a specified spectral model or multiple collection geometries. More wavelengths help only if the added data constrain the unknowns; allowing an arbitrary new $T$ and $R$ at every wavelength preserves the ambiguity.

Orbit determination infers a body's position and velocity from observations and an equation of motion. It is an inverse problem; a residual of an inverse problem compares an observed readout with its modeled prediction.

Hadamard well-posedness requires existence of a solution for every admissible datum, uniqueness, and continuous dependence of the solution on the datum in the specified normed spaces. An inverse problem is ill-posed if any of these requirements fails. The choice of spaces and norms is part of this assertion: continuity in one norm does not establish continuity in a stronger norm.

For a bounded linear operator $A$, existence requires $y\in \mathrm{ran}A$, and uniqueness requires $\ker A=\{0\}$. Even when both hold for exact data, inversion can be unstable. For example, if $\Vert x_N{\Vert }_X=1$ but $\Vert A{x}_N{\Vert }_Y\to 0$, data $A{x}_N$ approach zero while their solutions do not approach zero. Thus $A^{-1}$ on its range is not continuous. Perturbed data can also leave the range entirely, so an exact solution need not exist after measurement error. These are distinct mechanisms of ill-posedness; the inverse scattering problem below exhibits unstable continuation even on data for which an outgoing continuation exists.

A forward problem evaluates a known model at a state: given $u$ and $K$, compute $f=Ku$. The associated inverse problem takes the data $f$ and the model $K$ and seeks a state $u$ with $Ku=f$.

In the sense of Hadamard well-posedness, a problem is well posed if a solution exists for every admissible datum, is unique, and depends continuously on the datum in the specified input and output topologies. An inverse problem is ill-posed if at least one of these requirements fails. Thus ill-posedness means failure of existence, uniqueness or stability. The topologies and the admissible data set matter: restricting the data to the range restores existence, but need not restore continuous dependence.

Use the Hilbert space interpretation of a linear inverse problem, with a bounded linear operator $K:\mathcal{U}\to \mathcal{V}$. A least-squares solution minimizes $\Vert Ku-f{\Vert }_{\mathcal{V}}$ over $u\in \mathcal{U}$. A minimum-norm least-squares solution additionally has the smallest $\mathcal{U}$-norm among all least-squares solutions. When the equation is consistent this is the smallest-norm exact solution.

The least-squares solutions form a nonempty closed affine subspace precisely whenThe closest point theorem in a Hilbert space then gives a unique nearest point to zero, hence a unique minimum-norm least-squares solution. Equivalently, it is the unique least-squares solution in $\mathcal{N}(K{)}^{\perp }$. If $\mathcal{R}(K)$ is closed, this exists for every datum; a nonclosed operator range can leave some data without any least-squares solution.

For example, on the l2 sequence space, take $(Ku{)}_n=u_n/n$ and $f_n=1/n$. The equation would require $u_n=1$ for every $n$, which is not square summable. However, setting the first $N$ coordinates of $u$ to one and the rest to zero givesThus the least-squares infimum is zero but is not attained, and no minimum-norm least-squares solution exists for these data.

A problem is well posed in the sense of Hadamard when a solution exists for every admissible datum, is unique, and depends continuously on the datum. It is ill posed if any one of these three properties fails.

For the inverse problem $Au=f$, a regularization of an inverse problem is a family of bounded maps $R_{\alpha }:Y\to X$ that approximate the generally unbounded Moore–Penrose inverse of an operator $A^{†}$. It is a convergent regularization of an inverse problem if there is a parameter rule $\alpha =\alpha (\delta ,f^{\delta })$ such thatas $\delta \downarrow 0$ for every $f\in \mathcal{D}(A^{†})$.

A detector too slow to resolve a femtosecond field can measure its integrated nonlinear response instead. Split the pulse, delay one replica, overlap both in a thin phase-matched crystal, and select the sum-frequency signal. For negligible depletion and sufficiently broad phase matching, the cross-generated field is proportional to $A(t)A(t-\tau )$; integrated signal gives the second-harmonic intensity autocorrelationA noncollinear geometry can reject each replica's individual second harmonic. The optical response supplies the fast gate while the detector integrates slowly. Crystal thickness, group-velocity mismatch, spectral acceptance, beam overlap and delay calibration limit resolution.

For a Gaussian pulse specified by its intensity profile, the autocorrelation's full width at half maximum is $\sqrt{2}$ times the intensity full width at half maximum. For a sech-squared intensity profile the ratio is about $1.543$. These are model-dependent conversion factors. An intensity autocorrelation is symmetric in delay and loses phase information; different pulses can share it. An interferometric autocorrelation retains rapid carrier fringes but still does not generally reconstruct an arbitrary complex field uniquely.

frequency-resolved optical gating also resolves the gated signal's spectrum, giving a two-dimensional delay-frequency record. For second-harmonic frequency-resolved optical gating,where $\omega$ is measured relative to the signal's carrier. An iterative retrieval seeks an envelope whose predicted record matches the measurement. Alternatively, polarization-gating frequency-resolved optical gating uses the optical Kerr effect to gate the probe by a delayed replica's intensity, with signal envelope proportional to $A(t)|A(t-\tau ){|}^2$. The original FROG work introduced nonlinear optical gating plus spectral resolution to characterize arbitrary femtosecond pulses.^[ref]

Absolute wave phase and absolute time origin remain undetermined by these intensity records. The SHG record also cannot distinguish $A(t)$ from $A^{\ast }(-t)$, so its time-direction ambiguity remains. Retrieval agreement, an independently measured spectrum, adequate signal dynamic range and gate-response calibration are practical consistency checks. A pulse train that changes between shots can yield an averaged record that no single pulse faithfully represents. Thus nonlinear gating solves the detector-speed problem but does not remove inverse problem ambiguities or experimental limitations.

What is measured. A monoenergetic electron beam is directed at a thin target. A movable magnetic spectrometer selects the scattering angle and outgoing momentum, separating elastic scattering from nuclear excitation and other backgrounds. Incident beam current and target thickness supply the normalization of the differential scattering cross-section. The Hofstadter nuclear electron scattering experiment is a classic implementation of this electromagnetic probe.

For a heavy, approximately spin-zero atomic nucleus, neglecting recoil and treating the electron in the Born approximation, the point-charge result is multiplied by the nuclear charge form factor:Here $q$ is the magnitude of the spatial momentum transfer in natural units. The cosine factor includes the electron spin. Dividing the elastic data by the point-charge factor isolates $|F(q){|}^2$ within these approximations. Recoil, magnetic terms for targets with spin, radiative corrections and particularly Coulomb distortion in electron scattering must be included when significant.

From scattering to size and shape. For a spherical nuclear charge distribution,Thus the nuclear charge form factor is a normalized Fourier transform. Its low-$q$ expansion,measures the nuclear charge radius $r_{\mathrm{ch}}=\sqrt{⟨r^2{⟩}_{\mathrm{ch}}}$. Larger momentum transfers resolve the density profile through diffraction maxima and minima. For a uniform sphere of sharp radius $R$, the uniform-sphere nuclear charge form factor is $3[\sin (qR)-qR\cos (qR)]/(qR{)}^3$; its first zero is at $qR\simeq 4.493$, while $r_{\mathrm{ch}}=\sqrt{3/5}R$.

A practical extraction fits a family of densities to all measured cross sections, rather than treating a few measured values of $|F{|}^2$ as a complete density reconstruction. Finite angular coverage, resolution and the missing sign information make this an inverse problem.

Observed pattern. Apart from the lightest nuclei, the central density is approximately constant; the size grows approximately as $R\simeq 1.2A^{1/3}\mathrm{fm}$ and there is a diffuse surface. The coefficient describes an equivalent sharp nuclear radius, rather than the root-mean-square radius itself. Heavy nuclear charge distributions are often described by the two-parameter Fermi nuclear charge distribution, $\rho (r)={\rho }_0/[1+\exp ((r-c)/a)]$, with half-density radius $c$ scaling roughly as $A^{1/3}$ and diffuseness $a$ around $0.5\mathrm{fm}$. The surface thickness is then about $4a\ln 3\simeq 2.2\mathrm{fm}$. Shell effects, deformation and the lightest systems give departures from this smooth pattern. These qualitative features are established by the electron scattering density measurements.^[ref]

Complementary probes. Transition energies in a muonic atom are especially sensitive to nuclear size because the heavy negative muon orbits close to the charge. Electronic isotope shifts, after separating nuclear-mass and finite-size contributions, measure changes in $⟨r^2⟩$. Both require atomic-structure calculations. Electromagnetic excitation and transition strengths constrain deformation and charge multipoles. Hadronic scattering can also constrain nuclear dimensions, but carries uncertainties from the strong interaction and probes matter rather than purely electric charge.

Hearing requires estimating sound-producing events from noisy mixtures at two ears. It involves mechanical detection, neural coding, source segregation and recognition. A frequency analysis alone cannot identify who or what produced the sound.

The physical input is an acoustic pressure fluctuation. At one ear a simplified room model iswhere source signals $s_j$ are filtered by propagation and reflection paths $h_j$, and $n$ is noise. The convolution model assumes approximately linear, time-invariant propagation over the analysis interval; a moving source or listener can violate that simplification. This is a many-to-one inverse problem: different source configurations can produce similar ear signals. Distance, reverberation and head orientation alter a source's waveform. Pressure level is logarithmic: sound pressure level is $L_p=20{\log }_{10}(p_{\mathrm{rms}}/p_0)$, with reference $p_0=20\mu \mathrm{Pa}$ in air. A 100 dB difference corresponds to a pressure ratio $10^{100/20}=10^5$. The system must remain sensitive to weak signals while avoiding saturation over a large operating range.

The middle ear helps transfer airborne vibrations into cochlear fluids despite an acoustic impedance mismatch. Within the cochlea, the graded mechanics of the basilar membrane produce frequency-dependent traveling-wave responses: high frequencies peak nearer the base, low frequencies farther toward the apex. This tonotopy provides a distributed spectral representation. Deflection of sensory hair cells drives mechanotransduction; inner hair cells provide most afferent signaling, while outer hair cells contribute to mechanical amplification and tuning. Detection therefore depends on an active biological filter, not merely a passive bank of resonators.

A revealing experiment by Mario A. Ruggero and colleagues measured chinchilla basilar-membrane vibrations using the Mössbauer technique or Doppler laser velocimetry before and after a reversible intervention that reduced cochlear transduction currents. Responses near the characteristic frequency were reduced particularly strongly for weak sounds, with smaller changes at high levels and away from that frequency.^[ref] This supports feedback-dependent cochlear gain. The intervention also affects other cochlear processes, so it does not by itself identify every cellular contribution or directly measure human loudness. The important methodological contrast is controlled sound input plus measured mechanical output and a reversible physiological perturbation.

The resulting neural code must preserve several attributes simultaneously. Place of excitation and neural phase locking can provide complementary frequency and temporal information; firing rates, adaptation and population responses encode level and envelopes. Individual neurons have limited dynamic ranges, and temporal precision has limits. Pitch is not simply the location of the strongest spectral component, and timbre depends on spectral shape and temporal development. Even physically, a signal containing harmonics at 200 and 300 Hz repeats every 10 ms although it contains no 100 Hz component: both components complete integer cycles in that interval. Recovering source periodicity is therefore different from selecting a component in a frequency spectrum. It does not follow that every such signal has one obligatory perceptual pitch.

Localization adds ambiguity. Ear separation yields interaural time differences; head shadow yields interaural level differences; outer-ear filtering supplies direction-dependent spectral cues. In an ideal far-field point-receiver model, $\Delta t\simeq d\sin \theta /c$, obtained by dividing the path-length difference by sound speed. Real head diffraction modifies this approximation. One time difference can be compatible with several spatial directions, and phase is ambiguous when several cycles fit the delay. Head-related transfer functions, head movements and environmental context help resolve these ambiguities. Reverberation also requires distinguishing useful early information from later reflections.

Competing sources create both energetic masking, in which the target is poorly represented against overlapping acoustic energy, and informational masking, involving source confusion or selection despite available acoustic cues. Common onsets, harmonic relations, coherent modulation, timbre and spatial cues help group components into auditory objects. The grouping must remain stable when streams overlap, but also update when a source changes. These are the problems of auditory scene analysis.

E. Colin Cherry's classic experiments used simultaneous spoken messages and separate messages delivered to the two ears. Asking listeners to repeat one message while ignoring the other made source selection experimentally measurable.^[ref] Ear and voice differences can aid selection; the result does not show that unattended information is never processed, nor that real rooms provide the clean separation of headphone channels.

Finally, meaningful recognition requires tolerance of speaker, rate and context variation. Richard M. Warren removed a speech segment and substituted a cough-like sound. Listeners could hear the utterance as continuous and had difficulty locating the replacement, whereas a silent interruption was more readily detected.^[ref] Phonemic restoration shows that recognition uses constraints beyond the momentary waveform; it does not mean missing acoustic information has been physically recovered. Successful auditory perception combines peripheral sensitivity with grouping and context-sensitive inference, while retaining the possibility of masking, ambiguity and systematic error.

Seismic surface-wave dispersion. Rayleigh waves combine vertical and longitudinal motion in the vertical plane of propagation; Love waves have horizontal transverse motion and can be guided by a low-shear-speed surface layer. A homogeneous elastic half-space supports a Rayleigh wave with a frequency-independent speed in the ideal elastic limit, but depth-dependent Earth structure makes the surface-wave eigenfunctions and their propagation speeds depend on frequency. This structural dispersion can occur even when each constituent material is itself nondispersive.

Distinguish phase velocity $c=\omega /k$ from group velocity $U=d\omega /dk$, which controls a narrow-band envelope in the appropriate weakly attenuating limit. Longer periods generally have deeper sensitivity, while shorter periods emphasize shallower structure. In a profile whose shear speed increases with depth this often makes long-period waves faster, but low-velocity zones, anisotropy and higher modes can change that trend. Frequency dependence need not be monotonic in every Earth model.

Measure phase delays or wave-packet arrival times across known source–receiver paths for many periods, accounting for source phase, cycle ambiguity, multipathing and mode identification. Compute predicted dispersion from trial elastic profiles and compare with observations. Linearization expresses the sensitivity schematically aswhere $\alpha ,\beta$ are compressional and shear speeds. A regularized inverse problem then estimates depth structure; many crossing paths additionally provide surface-wave tomography. Joint Love wave and Rayleigh wave measurements help distinguish shear-wave polarizations and constrain radial seismic anisotropy. Finite bandwidth, depth averaging and correlated sensitivities mean that one dispersion curve does not uniquely determine a detailed Earth profile.

Seismic ray theory, the eikonal equation and Fermat's principle. In a smoothly varying isotropic medium, insert a high-frequency ansatz $u=A(\mathbf{x})e^{i\omega [T(\mathbf{x})-t]}$ into the scalar wave equation. The leading powers of frequency giveThe travel-time surfaces are wavefronts and rays are their normal trajectories, with slowness vector $\mathbf{p}=\nabla T$. The next-order equation transports amplitude, including geometrical spreading; attenuation is an additional material effect.

Fermat's principle says that the ray makes $\int ds/c$ stationary under path variations with fixed endpoints. Stationarity, rather than always the global minimum, permits multiple arrivals. Its Euler–Lagrange equation, with arclength $s$ and unit tangent $\hat{\mathbf{t}}$, isFor horizontal layering, the horizontal slowness $p_h=\sin i/c$ is conserved, giving Snell's law at interfaces and smooth refraction between them. Here $i$ is measured from vertical. A turning point satisfies $p_{h}c(z)=1$; increasing speed with depth bends a downgoing ray back upwards. Travel-time tomography uses these ray integrals to infer spatial slowness variations.

The approximation requires wavelength to be small compared with structural variation scales. Caustics, turning-point amplitudes, diffraction, scattering and finite-frequency sensitivity require wave-based corrections. Thus seismic ray theory gives an economical geometric description and useful travel times, while surface-wave dispersion retains essential interference and boundary-mode physics.

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.