This is a concise derivation for a model of QBO by reducing Laplace's tidal equations *along the equator*.

## 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

as

The terms are in longitude so that we can use a separation of variables approach and create a spatial standing wave for QBO, 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 .

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. Not including this approximation leaves us with an indeterminate set of equations.

Now consider that the QBO itself is precisely the term - the horizontal longitudinal velocity of the fluid, i.e. the wind in other words - which can be derived from the above by applying the solution to Laplace's third tidal equation in simplified form above:

So the *acceleration *of wind, not the velocity, is what obeys the Sturm-Liouville equation. The acceleration of wind is simply a *F=ma* response to a tractive gravitational forcing. A derivative *preserves the periods* of the Fourier components, but not the amplitudes, so what we see is a differently shaped envelope for QBO — i.e. one that is more spiky due to the time derivative applied.

## Part 2 : Deriving the lunar forcing periods

We apply a seasonal aliasing of the lunar gravitational pull to generate the latitudinal tractive forcing terms above. This turns into a set of harmonics which we derive next.

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 describe the tractive forcing.

## Application

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.

To create a multiple linear regression fit, we consider all the aliased Fourier factors that are known from the lunisolar tidal cycles. What the regression does is determine the weighting of the terms, across the set of lunar terms.

In the regression fit shown above, the training region is 1970 to 1982, and the rest is extrapolated. To get back the wind speed QBO, the curves are integrated:

Again the model fitting is only conducted on the interval from 1970 to 1982. Outside of that interval the model doesn't track every peak but enough of the fine detail is captured that it is quite obvious that the model has predictive power.

## Conclusions

This *acceleration* version of the QBO is spikier than the more smoothly oscillating conventional wind speed description of QBO. Yet this spikiness makes sense as it is physically representative of the seasonally modulated lunar forcing. A straightforward time integration of the acceleration model will recover the wind speed version of the QBO.

Besides the excursion-limiting behavior imposed by a modulation, 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.

By carefully applying these factors within the tidal formulation, we can generate a concise solution that has all the characteristics of the measured QBO. This needs to be adopted as the *de facto* QBO model which replaces the original formulation proposed by Lindzen in the late 1960's [1].

[1] R. S. Lindzen and J. R. Holton, “A theory of the quasi-biennial oscillation,” Journal of the Atmospheric Sciences, vol. 25, no. 6, pp. 1095–1107, 1968.