Put and . Acting on the Fourier mode, . Continuity and horizontal momentum give and . The buoyancy equation gives . Substitution into vertical momentum yields
The two signs describe the two frequency branches; when the intrinsic frequency is required positive, . For a stationary disturbance , so . With this requires and selects the minus sign, . The stationary dispersion condition is
Differentiate the unconstrained dispersion before imposing stationarity. The stationary mountain-wave group velocity is obtained from the ground-frame dispersion; before imposing stationarity it is
On the stationary locus this becomes
Differentiating after substituting as though all neighbouring modes were stationary would incorrectly erase the frequency gradients.
For upward energy propagation from the lower boundary, choose , since and . It requires real nonzero , hence
The downstream ray inclination above horizontal is
It is fixed for a given monochromatic forcing, with the two signs of giving the upward/downward rays.
Frame qualification for the requested phase angle: the normal intrinsic phase velocity is . Because , the boxed stationary group velocity is exactly : its angle to intrinsic phase propagation is , and it is perpendicular to the stationary phase lines. By contrast the intrinsic group velocity satisfies , so its angle to intrinsic phase propagation is . In the ground frame the phase velocity of a strictly stationary pattern is zero and has no nonzero propagation direction. For example , gives and . These distinctions prevent applying the usual zero-flow perpendicularity rule to two velocities expressed in different frames.
A mountain width is not an exactly specified sinusoidal wavelength. Taking its characteristic horizontal inverse length as with gives . Since , the stationary vertical wavenumber is
This is close to the long-horizontal-scale/hydrostatic estimate . If instead the quoted width is identified with a sinusoidal wavelength , then and . Thus several kilometres, roughly – under these two common width conventions, is the justified estimate; the terrain profile would fix its actual forcing spectrum.
For the inverse-width choice, both and are negative for upward energy radiation. The stationary phase lines tilt upstream with increasing height, while energy propagates upward and slightly downstream. Numerically , with ray inclination about . Intrinsic phase propagation points upstream and downwards; the ground pattern itself remains fixed. The original sketch below uses this explicit convention.
A ground observer may see stationary lenticular or wave clouds where ascent cools sufficiently moist air, with continuing air flow through a cloud whose position remains almost fixed. Descending regions can clear or warm. An aircraft crossing the pattern encounters alternating updraughts and downdraughts; gliders can exploit the lift. Strong amplitudes, shear or rotors may cause turbulence, but a linear gravity wave can also be smooth, and clear-air waves need not produce any visible cloud. Mountain height, humidity and nonlinear stability are not supplied, so neither a vertical-wind amplitude nor a definite atmospheric precipitation forecast follows.
Figure 1.
Stationary mountain-wave phase lines, intrinsic phase direction, energy rays and an upper-layer propagation cutoff
. Original schematic and dispersion calculation. The ridge height is illustrative; the phase sketch uses , while the upper-layer panel shows the cutoff for the reduced Brunt–Väisälä frequency. No cloud or turbulence amplitude is inferred. Details.
For constant , substitution of the vertical cosine and horizontal plane wave gives
The top rigid-lid condition is satisfied by . To select a discrete vertical mode one must also impose the corresponding no-flow, adiabatic condition at the flat bottom, . It requires . The mode is barotropic; the gravest nonzero baroclinic mode has . Hence
The lid condition alone is insufficient: also has zero derivative at the lid but violates the bottom condition. Constant depth motivates the additional physical boundary, which must be used to obtain the printed denominator.
Let and set . The zonal phase velocity and group velocity are
The phase curve is even in , with limiting value at zero and approach to from below at large . The group curve is also even: it starts at the same long-wave limit, equals at , reaches at , and approaches from above at large . The normalized curves appear in the bottom panels of the preceding figure. For , the phase speed means the limit ; there is no zonal crest propagation to define directly at exactly zero wavenumber.
Take a Northern-Hemisphere cyclonic vortex, with positive relative vorticity. For a fixed-depth incompressible layer, conservation of absolute vorticity gives
The anticlockwise flow on the east side moves parcels northward, increasing their planetary vorticity and creating negative relative vorticity anomalies. On the west side it moves parcels southward and creates positive relative vorticity anomalies. Initially,
so the positive anomaly grows on the western flank while the eastern flank loses positive vorticity. The cyclonic pattern consequently shifts westward. The induced anomalies themselves drive meridional flow and further relative vorticity anomalies, transmitting a Rossby wave pattern through the surrounding fluid. The restoring process is this exchange between planetary vorticity and relative vorticity, rather than a parcel executing a permanently westward translation.
The sketch shows the initial tendency, with the generated signs separated from the initial positive core. An isolated disturbance contains many wavevectors and disperses; westward phase velocity does not imply that all its wave energy travels westward. A Southern-Hemisphere cyclonic vortex reverses the signs of the initial and generated relative vorticities, while the same beta plane mechanism remains.
For and fixed , completing the square in the dispersion relation gives
These Rossby-wave isofrequency circles lie in the half-plane and touch the origin, where the dispersion relation is undefined. Their radii decrease as frequency increases. The zero-frequency line is .
The group velocity is normal to the constant-frequency curves, toward increasing frequency:
At , the group velocity is eastward; close to it is westward. For the positive-frequency circles its arrows point inward toward the circle centre. Its meridional component is southward in and northward in .
Figure 1.
Positive-frequency Rossby-wave isofrequency circles and group velocity directions in dimensionless wavenumber coordinates
. Details.
For , the zonal phase velocity, meaning propagation at fixed , is
It is westward everywhere it is defined. The zonal group velocity is eastward for , westward for , and zero on . There is no region of eastward zonal phase velocity. At , is undefined, though its continuous limiting expression is westward and . Using instead the normal-to-wavecrest phase velocity vector gives ; its zonal component is also westward for , but is not the same scalar as .
For a nonzero zonal wavenumber , stationarity requires
A real propagating solution therefore requires a westerly . For there are two meridional wavenumbers:
The source's assertion of two stationary solutions needs these qualifications: at they coincide, and at there is no real meridionally propagating solution. For example gives . Modes with are a separate zonally invariant degeneracy and do not have a defined zonal phase velocity.
Adding to the zonal group velocity and imposing stationarity gives
Thus a localized midlatitude disturbance radiates stationary Rossby wave activity downstream, along one poleward and one equatorward branch. For in the Northern Hemisphere, gives the northeastward branch and the southeastward branch. With strictly constant , these wave rays are straight. Changing latitude cannot bend them within that constant-coefficient model.
To discuss a planet, extend it locally to slowly varying and possibly . The conserved zonal wavenumber in a locally zonally uniform medium then has
For constant and northward propagation in the Northern Hemisphere, decreasing reduces until : the wave ray becomes zonal at a turning latitude and can refract equatorward. The equatorward branch initially encounters larger allowed and becomes more meridional. This local midlatitude description eventually needs replacement near the equator or wherever the background varies rapidly.
Orographic forcing by a mountain range or persistent spatially localized heating can supply such a disturbance. Its steady forcing selects the ground-stationary part of the response. The corresponding local total wavelength is
For fixed it increases toward the pole as decreases. Locally the zonal integer scale has upper bound , so fewer planetary-scale zonal waves fit around higher-latitude circles. This is not an unconditional trend in actual weather: changing , its curvature, stratification and the forcing spectrum can change the dominant Rossby wave wavelengths.
Put and . With and complex dependence , the northward momentum equation and continuity equation give
Thus . Write the signed phase velocity as , so . The eastward momentum equation is the transverse geostrophic balance . Hence
The decaying sign is selected by which side of contains the ocean. For an eastern boundary the ocean occupies , so decay offshore requires . Consequently the Kelvin wave moves northwards in the northern hemisphere and southwards in the southern hemisphere, with the coast respectively on its right and left. The positive offshore decay length is
The sign of reverses the propagation direction and the sign of relative to , but not the speed or offshore decay length. The real physical fields are the real parts of the complex solution. At the equator itself this local f-plane coastal solution loses its finite trapping scale.
Coastal Kelvin-wave radiation of long Rossby waves requires a common phase convention to match the two waves. The supplied Rossby wave formula conventionally multiplies , while the coastal expression uses . Since physical fields are real, matching may use the complex-conjugate representative; in either convention the boundary forcing fixes the absolute frequency and alongshore wavenumber. Thus and .
For , the shallow-water Rossby-wave dispersion relation reduces to
Choosing recovers the printed magnitude relation. Both the zonal phase velocity and zonal group velocity are approximately , so the long Rossby waves carry the coastal signal westwards into the basin.
The coastal condition alone does not guarantee a long zonal wavelength. Since , matching requires and hence the stronger consistency condition
For example, at a coastal wavenumber with gives , which is not asymptotically small even though the meridional wavelength is long. Thus the source's proposed excitation is possible at sufficiently low frequency; it is not an automatic consequence of its meridional long-wave assumption.
At N or S the beta plane parameter is
Taking a representative Pacific basin width m gives
A width of – m gives roughly – years. This is propagation near latitude, not an equatorial calculation. Slow Rossby wave adjustment supplies a basin-scale delay relevant to El Niño–Southern Oscillation: changing thermocline depth can affect later sea surface temperature and coupled wind forcing. Reflection and coupling may help produce interannual memory, but this travel-time estimate alone does not establish a climate oscillation or its period.
For group velocity from a period-dependent phase velocity, start with rather than equating group velocity to phase velocity. Differentiating,
For the stated inverse-frequency term, , so
The inequality assumes positive frequency and physical positive phase velocity, with . The symbol here has units of acceleration and is unrelated to the plate thickness in Question 1.
The phase velocity in Figure 1 is measured along the Earth's surface. Spherical geometry makes the local horizontal phase speed at radius equal to . In the ray approximation, squared radial shear wavenumber is
The spherical surface-wave evanescence threshold immediately below the discontinuity is therefore
Equality is the grazing limit. Using directly would omit the requested spherical correction. Immediately above the discontinuity the corresponding boundary is . Hence the interval specifically describes propagation just above the jump and evanescence just below it. Slower modes can already be evanescent above .
Intersecting the plotted branches with gives the following approximate upper periods. Below each upper period, that branch is evanescent immediately beneath the discontinuity:
For these overtones the displayed ranges thus start near and end at their respective crossings. Any still higher curve already above the threshold at the lower plotted period has no evanescent interval within the displayed range. Readings of the closely spaced higher curves are only accurate to about a second, and the plot does not establish their continuation to unshown periods. Nor does the local inequality prove that a mode remains evanescent at every greater depth in an unspecified Earth model.
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 from group velocity , 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 as
where 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 into the scalar wave equation. The leading powers of frequency give
The travel-time surfaces are wavefronts and rays are their normal trajectories, with slowness vector . The next-order equation transports amplitude, including geometrical spreading; attenuation is an additional material effect.
Fermat's principle says that the ray makes stationary under path variations with fixed endpoints. Stationarity, rather than always the global minimum, permits multiple arrivals. Its Euler–Lagrange equation, with arclength and unit tangent , is
For horizontal layering, the horizontal slowness is conserved, giving Snell's law at interfaces and smooth refraction between them. Here is measured from vertical. A turning point satisfies ; 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.
Figure 1.
Original illustrative Love-wave dispersion and depth sensitivity, alongside a turning ray
. Details.
For , phase velocity is the speed of a fixed phase, . The group velocity transports a narrow-band wave envelope. With ,
As with fixed, . Therefore
Both velocities tend to the shallow-water speed at leading order.
Set . The real part and imaginary part of the complex permittivity are
For positive frequency the loss part is nonnegative. The real part starts at , rises below resonance, falls steeply through the resonance region, and approaches one from below at high frequency. The imaginary part is a positive absorption peak near , tending to zero on both sides.
For a nonmagnetic medium write its complex refractive index as , with and optical extinction coefficient . In a transparent region with and ,
The real part controls phase velocity and optical dispersion; the loss part controls wave absorption and the cycle-averaged dissipated power density . Weak damping gives a narrow line but does not make wave absorption weak at its centre. The transparent-medium approximation applies away from the line.
Figure 1.
Single-resonance, multiple-resonance and free-electron optical responses
. Original dimensionless model plots. The first row shows one bound-electron resonance, the second shows three, and the third shows the Drude permittivity and its normal-incidence reflectance. Finite damping is retained; these are illustrative calculations, not material measurements. Details.
Assume a cold, unmagnetized, weakly collisional plasma, a small fractional bandwidth and . The dispersion relation is , with . The group velocity and propagation time are
The lower frequencies arrive later. Taking the magnitude of the frequency derivative over gives
For a nonuniform line of sight, replace by the dispersion measure . The light travels at the group velocity, not at the phase velocity.
Landau damping is the collisionless decay of a plasma wave through resonant energy exchange with particles whose velocity is close to the wave's phase velocity . Particles slightly slower than the wave tend to gain energy from it; particles slightly faster can give energy to it. For the decreasing resonant part of a Maxwell-Boltzmann distribution, the net transfer is from the wave to the electrons, so the wave amplitude decays.
In the linear Vlasov equation description this is collisionless phase mixing, not a requirement for frequent particle collisions. Fine-scale structure develops in the velocity dependence of the perturbation even while the macroscopic electric field decreases. The resonant contribution to the growth rate has the sign of at a positive wave phase velocity: a negative slope produces damping, whereas a sufficiently positive slope can drive an instability. A population comparison must be local to the resonance; merely counting all particles faster and slower than the wave is not the criterion.
The phase velocity is the speed of a constant-phase surface, . The group velocity is the envelope speed of a narrow packet, , requiring a local dispersion relation rather than only a single monochromatic wave. In an isotropic transparent medium with ,
For a vacuum wave, over the packet bandwidth, so both speeds are and . No superluminal signalling is implied by a phase speed in a dispersive material.
Strictly, the value at one frequency alone does not prove the product identity. For example gives but . The vacuum interpretation, or an equivalent nondispersive assumption, is necessary.
Write for the nonzero forcing frequency, to distinguish it from the vertical velocity . Put and . Then and , so it meets the rigid-wall conditions. Substituting into the scalar wave equation gives
or
The radiation sign for an inertial wave must exclude any disturbance arriving from . This means an radiation condition for propagating waves, or bounded decay for an evanescent wave.
For , put . The solutions are oscillatory. The inertial-wave dispersion relation is ; its horizontal group velocity is . Positive energy transport into the channel therefore requires . The solution is
For a positive frequency, the phase travels towards negative while the group velocity travels towards positive . Choosing the outgoing sign from the phase velocity alone would give the wrong result.
For , set . Excluding the exponentially growing solution gives
This disturbance decays away from the wavemaker and does not propagate horizontally. At the cutoff the ODE degenerates to ; zero frequency is not a nonzero forced harmonic mode with this vertical shape.
There is a second printed inconsistency: the stated velocity , if taken literally with no dependence and no transverse velocity, cannot support a nontrivial rotating inertial wave. Its transverse momentum equation would give , then continuity and the rigid walls would force . The scalar solutions above are the intended inertial-wave solutions when two-dimensional means dependence on while retaining the transverse component . For a mode, explicitly supplies the generally nonzero transverse velocity. Thus the scalar equation can be solved as requested, but the full vector field needs this correction to the printed two-component description.
For this density the Euler-Lagrange field equation with second derivatives gives
Insert a plane wave . The second and fourth spatial derivatives give and , respectively, so
For propagating modes, and the phase velocity has magnitude , which depends on the wave number. Thus the waves are dispersive. For , the squared frequency is negative: the corresponding solutions grow or decay exponentially rather than propagate with a real frequency. The uniform displacement is a zero-frequency shift mode.
Choose the rightgoing signed dispersion relation . The horizontal phase velocity and group velocity are
The leftgoing branch has their negatives. For a convention using positive frequency at both signs of , multiply each displayed speed by . For every nonzero , on the rightgoing branch; the horizontal group-speed bound for a confined internal wave is , approached as .
Pavel Cherenkov Created 2026-10-24 Updated 2026-11-05
Pavel Cherenkov investigated the visible radiation produced by fast charged particles in transparent materials. The emission's directional and polarization properties distinguished it from ordinary luminescence. Its interpretation as Cherenkov radiation relates the optical signal to a particle moving faster than the local phase velocity of light. His authored lecture connects the discovery with practical detectors, in which the presence and geometry of a light cone provide information about particle speed.[ref]