Much of tidal analysis has been performed by Fourier analysis, whereby one can straightforwardly deduce the frequency components arising from the various lunar and solar orbital factors. In a perfectly linear world with only two ideal sinusoidal cycles, we would see the Fourier amplitude spectra of Figure 1.
But in a non-linear world such as ENSO where the tidal forces interact with the seasonal cycle via modulated feedbacks the picture is quite different. What happens is that the cycles interact and get folded multiple times until what originally were three cycles (yearly plus draconic and anomalistic lunar periods) end up appearing as Figure 2.
Probably over 90% of the scientists and engineers with any signal processing expertise will start off with the ENSO signal and do a Fourier analysis as a first step. Then they will try to make sense of the spectra and try to identify cyclic components from the peak locations. Eventually they will scratch their heads in not being able to make any headway due to the richness in the spectra, and possibly attribute it to a chaotic time series — or perhaps describe it as a red noise subject to random wind burst stimuli. Unfortunately, that’s a misguided or at least shortsighted view of how to do a comprehensive time-series analysis. It’s well known that non-linear equations such as the Mathieu equation cause multiple splitting of the fundamental spectral peaks resulting in spectra with rich profiles and complex phase space trajectories. In particular see this recent Schroedinger equation analysis of harmonic generation .
Fortunately, I already knew of these issues from the solid-state world when I started working on the ENSO topic a few years ago — yet it still didn’t make it any easier to get to my current understanding. As one can see from Figure 1, it’s almost impossible to intuit that something as basic as three interacting sinusoidal factors can lead to such a rich spectra, but that’s exactly what happens when solving the delayed/Mathieu differential equation. It just took this long to capture the dynamics, simply because there aren’t a lot of crutches like a Fourier analysis to help out along the way. Everyone has the same problem, but it’s partly luck and partly persistence of trying to find the right combination that creates the best fit.
 C. Zagoya, M. Bonner, H. Chomet, E. Slade, and C. F. de Morisson Faria, “Different time scales in plasmonically enhanced high-order-harmonic generation,” Physical Review A, vol. 93, no. 5, p. 053419, 2016.