With this analysis, I wanted to demonstrate the underlying order of the most concise SOI Model. This model characterizes the salient fit parameters:

- Two slightly offset forcing sinusoids which match the average QBO forcing cycle
- A forcing sinusoid that maps to the frequency of the Chandler wobble beat
- A Mathieu modulation perturbing the 2nd-order DiffEq with a periodicity of about 8 years

This set of four parameters was used to model both modern day records corresponding to the atmospheric pressure data describing the Southern Oscillation Index, as well as to proxy records of historical coral data. The parameters seem to match closely over widely separated time intervals (see Figure 5 in the latter link).

**Figure 1** is the modern-day SOI record, suitably filtered to show the multi-year excursions.

It is amazing that this erratic a waveform can be modeled by a limited set of parameters that actually make some physical sense, but that is nature for you and the idea behind "sloppy modeling" -- models that use just a few parameters to accurately describe a behavior. The simplified model strongly suggests that there is a hidden harmony acting to drive ENSO.

To get to the results, **Figure 2** is the Mathematica model of the second-order Mathieu model describing hydrodynamic sloshing of the equatorial Pacific waters:

So even though the modeled waveform is erratic, the forcing factors are anything but ! (see **Figure 3** below). The Chandler wobble waveform (cw) appears weak but it is closer to the characteristic frequency for sloshing of ~4.3 years so has a greater resonant impact than the ~2.3 year average QBO period.

This fit was accomplished primarily manually via adjusting the sinusoidal coefficients to achieve the highest correlation coefficient.

As substantiation, the Eureqa machine learning tool (differential solution guide here) finds the same general set of parameters for a higher complexity solution as shown in **Figure 4**.

**Figure 5** is one selected Eureqa fit:

**Figure 6** is the QBO factor that Eureqa discovers:

The less complex Eureqa solutions feature the solutions shown in **Figure 3**.

What is interesting is that the periods for QBO variability, when aliased, fold into the synodic lunar cycle in which the actual time between lunations may range from about 29.18 to about 29.93 days. In** Figure 3**, one QBO frequency corresponds to 29.3 days and one to about 29.4 days. The average synodic cycle is 29.53 days, and 13 sideral months alias the same as 12 synodic months. See this post on possible impacts that the lunar force may have on QBO.

It may be true as well that there is some essential harmony behind the Chandler wobble, QBO, and lunar/tidal effects which are manifestations of the same geophysical process relating to periodic changes in angular momentum.

## References

-- see index to other SOI posts

I'm having a hard time replicating the filter you use on the NCAR SOI data. I've got the data and I've tried both running means and LOESS smooth, which can both give "pretty close" but not exact replication to what you've posted.

I'm sure you've mentioned your filter in some previous post somewhere, but if you can point me to where it would be helpful.

Keith,

Recently I have been combining the Darwin and Tahiti data directly to get SOI.

In Mathematica, I will do this

which is a 7 month filter. It ends up looking very much like the NINO3.4 curve.

All that is needed is a filter to remove the fast month-to-month changes, which is what the 7 month filter accomplishes.

In this comprehensive and practical guide, the secrets of the ancient science of Tantra become available to a contemporary audience for the first time. Confined to small, hidden mystery schools for centuries, and often misunderstood and misinterprete