(spoiler: If you have read Chapter 13 of Mathematical Geoenergy¹, you know where this is going)

The oceanic pole tide has historically been treated as a quasi-static equilibrium adjustment of the oceans to the Earth’s 433-day (~14 month) Chandler wobble. However, observed sea-level records present severe spatial discrepancies — most notably the extreme high-latitude amplification in semi-enclosed basins such as the Gulf of Bothnia in the northern Baltic Sea, where sea-level excursions exceed theoretical predictions by nearly an order of magnitude. This post reviews the historical discovery and mathematical derivation of the classical pole tide, examines why its classical tesseral potential ($\sin 2\theta$) fails to account for high-latitude maxima, and presents a physical and mathematical framework wherein the 14-month Chandler signal emerges as a stroboscopic, modular alias of the zonal ($m=0$) lunar draconic fortnightly tide gated by an annual impulse train. Not to say that the Chandler wobble and pole tide are not related; in fact, they are but the link is more to a common mechanism.

Discovery and Theoretical Origins

The discovery of the pole tide followed a classic sequence of astronomical deduction, geodynamic formulation, and subsequent tidal-gauge isolation.

In 1765, Leonhard Euler proved that any rotating, biaxial, rigid spheroid whose rotation axis departs slightly from its figure axis undergoes a free Eulerian nutation. For Earth’s flattening ratio:

$T_{\text{Euler}} = \frac{A}{C – A} \approx 305 \text{ days}$

where A and C are the equatorial and polar principal moments of inertia, respectively.

In 1891, the American astronomer Seth Carlo Chandler analyzed decades of precise stellar latitude observations, demonstrating that the true wobble had an amplitude of $\approx 0.7”$ and a period of approximately 433 days (~14 months), completely at odds with the rigid-body calculation.

In 1892, Simon Newcomb reconciled the 128-day gap by explaining that the Earth is not rigid. Instead, it yields elastically in its mantle and fluidly in its global oceans. As the instantaneous rotational axis wanders relative to the crust, the centrifugal potential shifts:

$\mathbf{\Omega}(t) = \Omega \begin{pmatrix} m_1(t) \\ m_2(t) \\ 1 + m_3(t) \end{pmatrix}$

The time-dependent centrifugal potential perturbation $\Phi_p$ acting on the fluid envelope is:

$\Phi_p = -\frac{1}{2} \Omega^2 r^2 \sin(2\theta) \left[ m_1(t)\cos\lambda + m_2(t)\sin\lambda \right]$

where:

• $\theta$ is geographic latitude,
• $\lambda$ is longitude,
• $m_1, m_2$ are the angular components of polar motion.

Newcomb termed this forced 14-month sea-surface deformation the pole tide.

Observational confirmation was not immediate as the equilibrium amplitude of this tide in the open ocean is exceedingly small:

$\zeta_{\text{eq}} \approx \frac{(1 + k_2 – h_2)\Phi_p}{g} \approx 5 \text{ to } 8 \text{ mm}$

where $h_2 \approx 0.61$ and $k_2 \approx 0.30$ are the Love numbers accounting for solid-Earth elastic yield and self-gravitation.

Despite the low signal-to-noise ratio, Christie (1895) in the United States and Bakhuyzen (1912) in the Netherlands successfully isolated a 14-month, 1-to-2 cm sea-level modulation from multi-decadal coastal tide records, establishing the pole tide as an empirical reality².

The wobble explanation of why the tide should vanish at the poles

The geometric signature of the classical pole tide is governed by the second-degree, first-order tesseral spherical harmonic:

$P_2^1(\sin\theta) = 3\sin\theta\cos\theta = \frac{3}{2}\sin(2\theta)$

The term $\sin(2\theta) = 2\sin\theta\cos\theta$ represents the product of two orthogonal projections, the z tilt projection and the $r_\perp$ lever arm:

              Z (Spin Axis)
              ▲
              │  Pole (θ = 90°): Lever arm $r_{\perp}$ = 0
              │  ──> No lateral centrifugal shift
              │
              │       • P(r, θ)
              │      /│
              │     / │
              │    /  │ z = r sin(θ)
              │   /   │ (Axis tilt projection)
              │  /    │
              │ /     │
              └───────┴──────► X, Y (Equator, θ = 0°)
              $r_{\perp}$ = r cos(θ)
              (Equator: z = 0 ──> No first-order tilt)

The radial distance to the rotational axis is $\cos\theta \propto r_\perp$. Centrifugal acceleration $F_c = \Omega^2 r_\perp$ scales with distance from the axis. Thus, at the geographic poles ($\theta = 90^\circ$), $r_\perp = 0$. Here, the spin axis merely tilts in place beneath the observer; with no perpendicular lever arm, the centrifugal perturbation is zero.

The coordinate projection along the axis is $\sin\theta \propto z$, in other words. the displacement of the tilted axis relative to the mean axis scales with vertical distance from the center of mass. In the equatorial plane ($\theta = 0^\circ$), z = 0, so that a pure tilt about the origin produces zero first-order variation in radial distance.

Because it requires both a non-zero height (z) and a non-zero lever arm ($r_\perp$), the perturbation achieves its maximum strictly at mid-latitudes:

$\frac{d}{d\theta}\sin(2\theta) = 2\cos(2\theta) = 0 \implies \theta = 45^\circ$

The classical theory dictates that the pole tide must decay toward zero at high latitudes.

The Bothnian anomaly defies the mechanism geometry

Tide-gauge networks throughout Northern Europe reveal an empirical distribution that directly inverts the classical prediction:

Empirically, instead of decaying as $\sin(2\theta)$ toward the Arctic, the amplitude expands monotonically up the Gulf of Bothnia, reaching its global coastal maximum in the sub-polar cul-de-sac of the Bothnian Bay in the northern Baltic.

Hydrodynamic and Meteorological Hypotheses

Physical oceanography has attempted to explain this discrepancy through local basin dynamics:

• Wind Setup Funneling: Southerly wind stress associated with the North Atlantic Oscillation (NAO) pushes water through the Baltic. As water enters the Bothnian Bay across the shallow Quarken sill ($h < 20 \text{ m}$), the one-dimensional momentum balance forces a steep surface slope: $\frac{\partial \eta}{\partial x} \approx \frac{\tau_w}{\rho g h}$

• Basin Seiche Antinodes: The Baltic acts as a long-wave resonator with an amphidromic nodal region near the Åland Islands and an antinode at the northern boundary.

While these local hydrodynamics account for local amplification factors of shallow water, they rely on an external weather forcing that happens to possess quasi-periodic power near 14 months. They do not explain why the fundamental forcing frequency itself matches the Chandler wobble across global decadal spectra.

Re-evaluating Causality: The Lunar Draconic Fortnightly Hypothesis

The severe shortcomings of the centrifugal $P_2^1$ potential at high latitudes invite a fundamental physical re-examination: What if the Chandler period itself is not an internal elastic resonance excited by random atmospheric noise, but rather an aliased response to an astronomical forcing that is naturally maximized at the poles?

First consider that direct gravitational tidal potentials do not share the $P_2^1$ mid-latitude cutoff if they are longitudinally symmetric (wavenumber m = 0), i.e. they are invariant across longitudes. This is because long-period tides generated by the Moon’s declinational excursions (the draconic or nodical tides) average over the diurnal cycle, eliminating longitudinal dependence. The resulting potential is governed by the Legendre polynomial of degree 2, order 0:

$P_2^0(\sin\theta) = \frac{1}{2}\left(3\sin^2\theta – 1\right)$

Evaluating $P_2^0(\sin\theta)$ across latitudes:

• Equator ($\theta = 0^\circ$): $P_2^0(0) = -0.5$
• Nodal zero-crossing $(\theta \approx 35.26^\circ): 3\sin^2\theta – 1 = 0$
• North Pole $(\theta = 90^\circ): P_2^0(1) = \mathbf{+1.0}$

A forcing with wavenumber m = 0 doubles in magnitude at the pole relative to the equator and increases monotonically across the sub-polar and polar seas.

The stroboscopic aliasing mechanism is the spoiler in explaining everything. Consider first that the Moon’s nodal crossing cycle—the draconic month—has a period of:

$T_d \approx 27.2122 \text{ days}$

Because tidal deformation depends on the square of the declination ($\sin^2\delta$) or is symmetric across the equatorial plane, the primary declinational forcing recurs every half-cycle, known as the draconic fortnightly cycle:

$T_{df} = \frac{T_d}{2} \approx 13.6061 \text{ days} \quad \implies \quad f_{df} = \frac{1}{T_{df}} \approx 0.0735 \text{ cycles/day (cpd)}$

Now consider the coupling of this continuous zonal carrier to an annual impulse train or gating envelope (e.g., seasonal stratification, solar insolation, annual sea-ice formation, or annual air-mass exchange):

$T_a \approx 365.2422 \text{ days} \quad \implies \quad f_a = \frac{1}{T_a} \approx 0.00274 \text{ cpd}$

Modular arithmetic of the interaction generates the aliasing. The ratio of annual to draconic fortnightly cycles defines the phase progression:

$N = \frac{T_a}{T_{df}} = \frac{365.2422}{13.6061} \approx 26.844 \text{ cycles/year}$

When sampled or gated stroboscopically at the annual interval $T_a$, frequencies fold modulo the annual sampling frequency. The fractional residual $\Delta \phi$ is:

$\Delta \phi = N \pmod 1 = 26.844 – 26 = 0.844 \text{ cycles/year}$

Computing the beat period corresponding to this residual:

$T_{\text{beat}} = \frac{T_a}{\Delta \phi} = \frac{365.2422 \text{ days}}{0.84406} \approx \mathbf{432.72 \text{ days}}$

This result is mathematically identical to the observed Chandler wobble period ($\sim 433d$).

Comparison
The hypothesis that the Chandler period and its corresponding sea-level signal originate via modular aliasing of the lunar draconic tide resolves the high-latitude amplification paradox:

As with all these models, the best way to proceed is to compare to actual data. Using the approach as described in Chapter 12 Mathematical Energy on Wave Energy, we set up the convolution of the tidal sinusoids with the annual impulse comb to generate a forcing manifold. Then apply the LTE modulation to generate a type of winding response.

The following is for the Kemi MSL station in northern Finland at the upper reaches of the Bothnian Gulf

Oulu in Finland

Most of these stations seem to be in phase, so the set of Baltic time-series can be averaged and fit again.

Note that in the above I have been using a multiple-leg fitting approach, whereby I take a common forcing manifold for NAO and the MSL data. The windings are essentially 180 degrees apart which makes sense from the inverted barometric effect but it may be that the dipole measurement is reversed here. Whichever it may be, both the atmospheric NAO and the tidal stations do indeed have a common mechanism.

The amplitude of the windings act essentially as the scaling of the tidal values. Both the fundamental winding value and the harmonics show the same latitudinal dependence as cited elsewhere.

a

The Arctic Oscillation fit is similar to the NAO, with the same set of windings. Only the phase is reversed below, with the NAO having a low-pass filter applied, so that the faster windings come through on AO.

The logical next connection to test is an index at the other pole, so here is the southern annular mode (SAM)

The windings for NAO are quite apparent, below from the winding scalogram applied:

Conclusion

It may be all tidal forcing after all.

The designation of a pole tide demonstrates the challenge of explicating the observations from what is expected form spatial potential analysis. While the Chandler wobble itself undoubtedly tilts the centrifugal potential in a $P_2^1$ harmonic, the inability of that harmonic to explain high-latitude sea-level data points toward another mechanism. I suppose it could be amplification from wind funneling or seiche antinodes, but the common-mode mechanism of draconic lunar torque driving both CW and the MSL cycle is more parsimonious and possibly more plausible.

Consider that the modular aliasing of the $m=0$ lunar draconic fortnightly tide across the annual solar cycle provides both:

• An exact mathematical derivation of the 433-day period ($365.2422 / 0.84406 \approx 432.7 d$).
• A zonal $P_2^0$ potential geometry whose intrinsic maximum resides at the poles, resolving the Bothnian amplification paradox from first physical principles

Perhaps it could still be called a pole tide, but if the polar wobble doesn’t cause it the naming would now be misguided. Consider the perhaps absurd notion that the Chandler wobble is on the same scale as the LOD changes signifying variation of angular motion — this is orders of magnitude larger than the CW kick, yet there is no such thing as an “LOD tide” (except by me of course, see Chapter 12 in the reference below, where it explains ENSO and AMO) .

References

1. P. Pukite, D. Challou, and D.Coyne, Mathematical Geoenergy, (Wiley, 2018), https://agupubs.onlinelibrary.wiley.com/doi/book/10.1002/9781119434351 ↩︎
2. S. R. Dickman, The self-consistent dynamic pole tide in global oceans, Geophysical Journal International, Volume 81, Issue 1, April 1985, Pages 157–174, https://doi.org/10.1111/j.1365-246X.1985.tb01356.x https://watermark02.silverchair.com/81-1-157.pdf ↩︎