Laplace developed his namesake tidal equations to mathematically explain the behavior of tides by applying straightforward Newtonian physics. In their expanded form, known as the primitive equations, Laplace's starting formulation is used as the basis of almost all detailed climate models. Since that's what they are designed to do, this post provides the details for solving Laplace's tidal equations in the context of the El Nino Southern Oscillation (ENSO) of the equatorial Pacific ocean. The derivation and results shown below essentially describe the framework of my presentation at this month's AGU meeting: Biennial-Aligned Lunisolar-Forcing of ENSO: Implications for Simplified Climate Models
The concise derivation for a model of ENSO depends on reducing Laplace's tidal equations along the equator. I could not find anyone taking a similar approach anywhere in the literature, even though it appears to be routinely obvious: (1) solve Laplace's tidal equations in a simplified context, then (2) apply the known tidal forcing and observe if the result correlates or matches the ENSO time series. In fact, it does, as I have shown before (and for QBO as well); but this is the first time that I have worked out the details in full for ENSO. Below is a two part solution.
Part 1 : Deriving a closed-form solution
For a fluid sheet of average thickness , the vertical tidal elevation , as well as the horizontal velocity components and (in the latitude and longitude directions).
This is the set of Laplace's tidal equations (wikipedia). Along the equator, for at zero we can reduce this.
where is the angular frequency of the planet's rotation, is the planet's gravitational acceleration at the mean ocean surface, is the planetary radius, and is the external gravitational tidal-forcing potential.
The main candidates for removal due to the small-angle approximation along the equator are the second terms in the second and third equation. The plan is to then substitute the isolated and terms into the first equation, after taking another derivative of that equation with respect to .
Taking another derivative of the first equation:
Next, on the bracketed pair we invert the order of derivatives (and pull out the constant )
Notice now that the bracketed terms can be replaced by the 2nd and 3rd of Laplace's equations
The terms are in longitude so that we can use a separation of variables approach and create a spatial standing wave for ENSO (i.e. the sloshing model), where is a wavenumber.
That is essentially close to a variation of a Sturm-Liouville equation (see right). To solve this for a plausible set of boundary conditions, we make a connection between a change in latitudinal forcing with a temporal change:
After properly applying the chain rule, this reduces the equation to a function of and , along with a constant A. The A subsumes the wavenumber SW(s) portion so there will be multiple solutions for the various standing waves, which will be used in fitting the model to the data.
So if we fix to a periodic function with a long-term mean of zero
to describe the tractive latitudinal displacement terms near the equator, then the solution is the following:
where is an aggregate of the constants of the differential equation and represents the fixed phase offset necessary for aligning on a seasonal peak. This approximation of a tangential tractive force as a cyclic displacement for is a subtle yet very effective means of eliminating a big unknown in the dynamics ( this is essentially similar to a Berry phase applied as a cyclic adiabatic process, with the phase being the internal sin modulation ). Not including this approximation leaves us with an indeterminate set of equations.
Now consider that the ENSO itself is precisely the term - the horizontal longitudinal acceleration of the fluid, i.e. leading to the observed sloshing in the thermocline - which can be derived from the above by applying the solution to Laplace's third tidal equation in simplified form above:
There is also a cosh solution for when A is negative
Essentially this result is simply a F=ma response to a tractive gravitational forcing which obeys Laplace's tidal equations. The cos solutions show a negative feedback while the cosh solutions are the positive feedback variation.
Part 2 : Deriving the lunar forcing periods
One of the features of the ENSO time-series is a strong biennial component. To model this characteristic, we apply a seasonal aliasing of the lunar gravitational pull to generate the terms needed as a forcing stimulus. This turns into a set of harmonics which we can fit the data to.
The starting premise is that a known lunar tidal forcing signal is periodic
The seasonal signal is likely a strong periodic delta function, which peaks at a specific time of the year. This can be approximated as a Fourier series of period .
For now, the exact form of this doesn't matter, as what we are trying to show is how the aliasing comes about.
The forcing is then a combination of the lunar cycles amplified in some way by the strongly cyclically peaked seasonal signal .
Multiplying this out, and pulling the lunar factor into the sum
then with the trig identity
Expanding the lower frequency difference terms and ignoring the higher frequency additive terms
Thus we can now understand how the high frequency lunar tidal terms get reduced in frequency by multiples of , until it nears the period of the seasonal cycle. These are the Fourier harmonics of the aliased lunar cycles that comprise the forcing. The aliasing explains why we do not see the ENSO cycles at the monthly scale, but instead observed cycles at the multi-year scale.
In a more precise fashion, we can apply the known gravitational forcing from the lunar orbit and the interaction with the sun's yearly (lunar-perturbed) cycle. This works very effectively, and the closer one can get to the precise orbital path, the better resolved the fit.
An interesting mathematical consideration is how best to adopt the biennial modulation. Either of these mechanisms work equivalently to some degree: (1) alternate the sign of the biennial forcing impulse, (2) apply a 2-year Mathieu modulation to the DiffEq, or (3) apply a 1 year-delay as a delayed differential equation. These are equivalent in the sense that similar results can be obtained if only one of these are used or they are all used, but each scaled by 1/3. For the Laplace formulation, approach (1) provides the most convenient approach because it does not require a rework to the closed-form solution.
The Part 1 derivation provides the closed-form natural response and Part 2 provides the boundary condition forcing terms due to the lunar nodal cycle. As an example of a typical fit with this approach, we apply the gravitational forcing described here:
where a(t) and d(t) are the composed monthly+fortnightly anomalistic and draconic lunar cycles. This generates a rich set of harmonics that expand as a Fourier series used as input to the Laplace solution.
In the regression fit shown above, the training region is 1880 to 2000, and the rest is extrapolated.
An excursion-limiting behavior is imposed by a modulation, and this formulation can also show amplitude folding at the positive and negative extremes. In other words, if the amplitude is too large, the outer modulation starts to shrink the excursion, instead of just limiting it. When the amplitude starts folding and bifurcating in this fashion (as A grows larger), it falls into the behavioral category identified as wave breaking. As a likely application of this breaking, we can substitute a larger A factor in Laplace's solution to model the more rapid cycles observed in the ENSO SOI time-series; these would arise from a differing wavenumber (likely larger index) standing wave solution to the sloshing. The results of this interpretation were shown in the last post on a HiRes ENSO model . This enables a fit over a wide dynamic range of time-scales without adding many additional degrees of freedom, i.e. an altered A term and phase.
As with QBO, by carefully applying these factors within the tidal formulation, we can generate a concise solution that has all the characteristics of the measured ENSO. The main challenge in store for others desiring to duplicate this work is the necessary meticulous bookkeeping of the orbital cycles, keeping track of all the second-order perturbations to the lunar cycles. This can also be validated against sensitive LOD measurements of the earth's slight angular rotation variations.
This book  provides an algorithm for the Anomalistic Earth-Moon distance perturbation (shown below), where an accurate mapping to the valley excursions has the strongest gravitational impact and thus is most sensitive to a good fit.The ENSO model based on Laplace's tidal equation needs to be carefully considered as a physically-sound alternative to the chaotic network model published by Tsonis in a highly influential physics journal . As is the case of Lindzen's model of the QBO, a generation of climate science has taken Tsonis' lead down a rabbit hole of unnecessary complexity in understanding ENSO.
 Astronomical Algorithms, Jean Meeus, 2nd edition (December 1998), Willmann-Bell, ISBN: 0943396638
 Topology and Predictability of El Niño and La Niña Networks, Anastasios A. Tsonis and Kyle L. Swanson, Physical Review Letters 100, 228502 – Published 5 June 2008