Paul L. Vaughan, M.Sc.
Without a good handle on its simple geometry, a seemingly complex time series can appear as a changeling yielding to the pressures of mysterious statistical manipulation.
For example, a fundamentally important seminal observation reported by Le Mouël, Blanter, Shnirman, & Courtillot (2010) revealed the quasistationary 11 year solar cycle in the rate of change of length of day (LOD’), but newcomers taking a preliminary look at daily resolution LOD’ are more likely to fixate on the 18.6 year lunisolar envelope.
Multiscale variance summaries highlight obvious envelopes:
Zooming in, a semi-annual envelope is also evident:
(WIDE GRAPH ABOVE –Click to view elongate graph^1 & then click again to magnify.)
(WIDE GRAPH ABOVE –Click to view elongate graph^2 & then click again to magnify.)
A parsimonious weekly-to-monthly timescale model of daily LOD’, explaining ~93% of the variance (r = 0.965), can be constructed using the following information (with model terms in bold italics):
| Year | Period (days) | Half-Period (days) | Defined by… |
| Tropical | 365.24219 | 182.621095 | equinoxes |
| Lunar Month | Period (days) | Half-Period (days) | Defined by… |
| Tropical | 27.321582 | 13.660791 | equator/equinoxes |
| Nodal or Draconic | 27.212221 | 13.6061105 | ecliptic |
| Anomalistic | 27.55455 | 13.777275 | apogee/perigee |
| Synodic | 29.530589 | 14.7652945 | new/full moon |
(27.321582)*(27.212221) / (27.321582 – 27.212221)
= 6798.410105 days = 18.61343046 years
(6798.410105)*(13.6061105) / (6798.410105 – 13.6061105)
= 13.63339592 days
(27.55455)*(13.660791) / (27.55455 + 13.660791)
= 9.132933018 days
Noteworthy envelopes apparent in the variance structure of LOD’ relate to:
1) lunar nodal cycle (LNC) = 18.6 years
2) lunar apse cycle (LAC) = 8.85 years
3) terrestrial year (1 year)
4) harmonics (e.g. 0.5 years & 4.42 years)
| Beat Period | (years) | Tropical | Nodal | Anomalistic | Synodic |
| 27.321582 | 27.212221 | 27.55455 | 29.530589 | ||
| Tropical | 27.321582 | – | 18.6134 | 8.8475 | 1.0000 |
| Nodal | 27.212221 | 18.6134 | – | 5.9970 | 0.9490 |
| Anomalistic | 27.55455 | 8.8475 | 5.9970 | – | 1.1274 |
| Synodic | 29.530589 | 1.0000 | 0.9490 | 1.1274 | – |
| Beat Period | (years) | Tropical/2 | Nodal/2 | Anomalistic/2 | Synodic/2 |
| 13.660791 | 13.6061105 | 13.777275 | 14.7652945 | ||
| Tropical/2 | 13.660791 | – | 9.3067 | 4.4238 | 0.5000 |
| Nodal/2 | 13.6061105 | 9.3067 | – | 2.9985 | 0.4745 |
| Anomalistic/2 | 13.777275 | 4.4238 | 2.9985 | – | 0.5637 |
| Synodic/2 | 14.7652945 | 0.5000 | 0.4745 | 0.5637 | – |
Beat Period = (A*B) / ( |A-B| )
| | indicates absolute value
The model:
| Relative | Cumulative | ||||
| Term | Period (days) | Amplitude | r^2 | r | Contribution |
| 1 | 13.660791 | 1 | 0.713 | 0.844 | | polarity | |
| 2 | 13.63339592 | 0.41 | 0.824 | 0.908 | LNC |
| 3 | 9.132950896 | 0.30 | 0.881 | 0.939 | LAC alternation |
| 4 | 27.55455 | 0.26 | 0.926 | 0.962 | LAC alternation |
| 5 | 14.7652945 | 0.08 | 0.931 | 0.965 | semi-annual |
(WIDE GRAPH ABOVE – Click to view elongate graph^3 & then click again to magnify.)
eLOD’ = estimated LOD’
The above tables & figures, while certainly nothing new to science, have been summarized here for the benefit of those striving to efficiently develop the foundations necessary to appreciate and build upon the recent seminal work of Le Mouël, Blanter, Shnirman, & Courtillot (2010). From their conclusions:
“The solid Earth behaves as a natural spatial integrator and time filter, which makes it possible to study the evolution of the amplitude of the semi-annual variation in zonal winds over a fifty-year time span. We evidence strong modulation of the amplitude of this lod spectral line by the Schwabe cycle (Figure 1a). This shows that the Sun can (directly or undirectly) influence tropospheric zonal mean-winds over decadal to multi-decadal time scales. Zonal mean-winds constitute an important element of global atmospheric circulation. If the solar cycle can influence zonal mean-winds, then it may affect other features of global climate as well […]”
[Typos: 1) “evidence” should read “observe”. 2) “undirectly” should read “indirectly”.]
Caution
Exclusive &/or excessive focus on the first moment (the mean) should not be at the expense of attention to higher moments (such as the variance), as the following graph should emphasize:
SOI = Southern Oscillation Index (an index of El Nino / La Nina)
[ ] indicates boxcar averaging [applied here to highlight interannual variability]
When studying the preceding graph, it is important to understand that the blue line is the normalized interannual average of the black line. (Take a minute to think about this carefully.)
To reinforce this point, here is another graph of the normalized mean at the semi-annual to annual timescale:
The occurrence of such patterns in the mean despite the maintenance of stationary variance limits suggests a need to carefully consider which equators (geographic, celestial, magnetic, meteorological, etc.) are relevant to the phenomena under study. (See for example Leroux (1993).)
Multimoment multiscale spatiotemporal integration reveals nonrandom harmonic pattern-summary discontinuities, exposing the comedy tragically advocated by deceitful &/or naive theoreticians who are in part constrained by a dominant culture that clings seemingly religiously to maladaptive traditions such as unjustifiable assumptions of randomness, independence, uniformity, linearity, etc. that are routinely misapplied (for example to conveniently render abstract conceptions mathematically tractable).
Bear in mind that for some phenomena, such as ice-jacking freeze/thaw cycles, the properties of the variance play a critically fundamental role in dynamics.
Conclusion
With awareness of key wavelengths and a solid conceptual understanding of the effect of integration across harmonics, we arrive at something truly simple: Earth, Sun, Moon.
Both of the ~11 year waves summarize the semi-annual wave, which summarizes biweekly & monthly LOD’ variations bounded by lunisolar limits.
While the magenta wave is isolated via complex wavelet methods, the sky-blue wave is accessible to any member of the general public with an understanding of this article, 5 minutes to spare, & a spreadsheet.
Acknowledgement
Tim Channon generously shared LOD’ models developed using his synthesizer software. Access to Tim’s models facilitated expeditious cross-checking of lunisolar theory, mainstream literature, & data.
Suggestion
I encourage responsible readers to download & archive daily LOD data. Scientifically-engaged citizens can keep a vigilant watch on potentially-arising future data vandalism.
Data
LOD
International Earth Rotation Service (IERS)
http://www.iers.org/IERS/EN/DataProducts/EarthOrientationData/eop.html
Related Reading
Li, G.-O.; & Zong, H.-F. (2007). 27.3-day and 13.6-day atmospheric tide. Science in China Series D – Earth Sciences 50(9), 1380-1395.
http://www.scichina.com:8080/sciDe/fileup/PDF/07yd1380.pdf
Sidorenkov, N.S. (2007). Long-term changes in the variance of the earth orientation parameters and of the excitation functions.
http://syrte.obspm.fr/journees2005/s3_07_Sidorenkov.pdf
Sidorenkov, N.S. (2005). Physics of the Earth’s rotation instabilities. Astronomical and Astrophysical Transactions 24(5), 425-439.
http://images.astronet.ru/pubd/2008/09/28/0001230882/425-439.pdf
Gross, R.S. (2007). Earth rotation variations – long period. In: Herring, T.A. (ed.), Treatise on Geophysics vol. 11 (Physical Geodesy), Elsevier, Amsterdam, in press, 2007.
http://geodesy.eng.ohio-state.edu/course/refpapers/Gross_Geodesy_LpER07.pdf
http://geodesy.geology.ohio-state.edu/course/refpapers/Gross_Geodesy_LpER07.pdf
Schwing, F.B.; Jiang, J.; & Mendelssohn, R. (2003). Coherency of multi-scale abrupt changes between the NAO, NPI, and PDO. Geophysical Research Letters 30(7), 1406. doi:10.1029/2002GL016535.
Maraun, D.; & Kurths, J. (2005). Epochs of phase coherence between El Nino-Southern Oscillation and Indian monsoon. Geophysical Research Letters 32, L15709. doi10.1029-2005GL023225.
http://www.cru.uea.ac.uk/~douglas/papers/maraun05a.pdf
Leroux, M. (1993). The Mobile Polar High: a new concept explaining present mechanisms of meridional air-mass and energy exchanges and global propagation of palaeoclimatic changes. Global and Planetary Change 7, 69-93.
http://ddata.over-blog.com/xxxyyy/2/32/25/79/Leroux-Global-and-Planetary-Change-1993.pdf
Trenberth, K.E.; Stepaniak, D.P.; & Smith, L. (2005). Interannual variability of patterns of atmospheric mass distribution. Journal of Climate 18, 2812-2825.
http://www.cgd.ucar.edu/cas/Trenberth/trenberth.papers/massEteleconnJC.pdf
Abarca del Rio, R.; Gambis, D.; & Salstein, D.A. (2000). Interannual signals in length of day and atmospheric angular momentum. Annals Geophysicae 18, 347-364.
http://hal-insu.archives-ouvertes.fr/docs/00/32/91/24/PDF/angeo-18-347-2000.pdf
Abarca del Rio, R.; Gambis, D.; Salstein, D.; Nelson, P.; & Dai, A. (2003). Solar activity and earth rotation variability. Journal of Geodynamics 36, 423-443.
http://www.cgd.ucar.edu/cas/adai/papers/Abarca_delRio_etal_JGeodyn03.pdf
Le Mouël, J.-L.; Blanter, E.; Shnirman, M.; & Courtillot, V. (2010). Solar forcing of the semi-annual variation of length-of-day. Geophysical Research Letters 37, L15307. doi:10.1029/2010GL043185.
Vaughan, P.L. (2010). Semi-annual solar-terrestrial power.
Technical Aside
For those interested in exploring LOD’ variance patterns that are not necessarily evident at first glance, another noteworthy envelope is the following:
(13.777275)*(13.63339592) / (13.777275 – 13.63339592)
= 1305.478517 days = 3.574281812 years
This polar-equatorial eclipse cycle is evident in the sequence of diagrams here:
http://eclipse.gsfc.nasa.gov/5MCLE/5MCLE-Figs-10.pdf (1733-2151)
From:
Espenak, F.; & Meeus, J. (2009). Five millennium canon of solar eclipses: -1999 to +3000 (2000 BCE to 3000 CE). NASA Technical Publication TP-2009-214172.
http://eclipse.gsfc.nasa.gov/SEpubs/5MCLE.html
h/t to WUWT commenter “lgl” for initially drawing attention to this pattern some time ago.
Earlier & Future Articles
I wrote the following articles before (a) acquiring access to Le Mouël, Blanter, Shnirman, & Courtillot (2010), (b) coming across Leroux (1993), and (c) re-reading Sidorenkov (2005) with consequently improved awareness:
1) http://wattsupwiththat.com/2010/08/18/solar-terrestrial-coincidence/
2) http://wattsupwiththat.com/2010/09/04/the-north-pacific-solar-cycle-change/
3) http://wattsupwiththat.com/2010/09/11/solar-cycle-length-its-rate-of-change-the-northern-hemisphere/
Related articles could have been written on All India Rainfall Index & other variables, but the audiences’ handle on the solar, lunisolar, & spatiotemporal nature of interannual variations was revealed to be inadequate in comments here:
4) http://wattsupwiththat.com/2010/10/11/atlantic-hurricanes-the-sun/
[Some audience members may benefit from careful consideration of issues raised by Tomas Milanovic at Dr. Judith Curry’s blog Climate Etc.]
Le Mouël, Blanter, Shnirman, & Courtillot’s (2010) game changing observation rendered earlier results much less mysterious:
For capable individuals striving to render these & related findings disgestible by a mainstream audience, I strongly recommend:
A) gleaning the primary point made by Schwing, Jiang, & Mendelssohn (2003) about the effect of windowing parameters on apparent phase, which can be reversed by spatial patterns, not just temporal evolution.
B) heeding the advice of Maraun & Kurths (2005) about “periods of coupling which are invisible to linear methods.”
Future posts in this series (if it continues) may draw attention to:
a) nonrandom relations between interannual terrestrial oscillations and interannual [not to be confused with decadal] rates of change of solar variables.
b) the guaranteed potential for naive investigators to be irrecoverably derailed by Simpson’s Paradox due to stubborn &/or blind adherence to seriously misguided conventional mainstream statistical inference paradigms & malpractices that rigidly & dogmatically insist on falsely assuming independence when none exists.
c) the [counterintuitive &/or paradoxical for some] influence of grain & extent – & aggregation criteria more generally – on summaries of spatiotemporal pattern.
“Grain” & “extent“?…
Grain is another term for spatiotemporal resolution. Important: Extent is a term which concisely encompasses the properties of spatiotemporal summary windows. The vast majority of mainstream researchers are either absolutely ignorant or insufficiently cognizant of the effect of extent on integrals across spatiotemporal harmonics (including the nonstationary variety). The consequences are serious: blindness and rejection of valid findings on nonsensical grounds.
Best Regards to All.



Multimoment multiscale spatiotemporal integration reveals nonrandom harmonic pattern-summary discontinuities
that’s uhm… EXACTLY what I was thinking as I was reading along.
Thank you for this, Paul.
Have you had a chance to review “Solar Minima, Earth’s rotation and Little Ice Ages in the past and in the future The North Atlantic–European case”?
I think you will find it quite interesting if you have not.
http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6VF0-4Y7P4NS-2&_user=10&_coverDate=07%2F31%2F2010&_rdoc=1&_fmt=high&_orig=gateway&_origin=gateway&_sort=d&_docanchor=&view=c&_searchStrId=1713037503&_rerunOrigin=google&_acct=C000050221&_version=1&_urlVersion=0&_userid=10&md5=591a5ac3b099dd4c0ac2ff0f9ec45a36&searchtype=a
Email me at sharkhearted@gmail.com and I will send you a copy.
Chris
Norfolk, VA, USA
Had me sufficiently confused/absorbed to keep reading till the end.
Wow, I read the whole thing. The graphics look good, but what is with the techno speak, I need a “Solar, Terrestrial, & Lunisolar Components of Rate of Change of Length of Day” for DUMMIES post.
I’d love to know how this LOD change on the surface interacts with the mantle and core. That surface change can’t possibly be affecting the iron core, I would think, so the mantle has to absorb the energy of the friction. I understand the core normally has a different period of rotation than does the surface so that has to create quite a bit of churn in the mantle. Illustrations of the core tend to show it smooth and round as a pool ball but I suspect it is actually quite irregular. I suspect too that regional viscosity/temperature variations of the mantle perturb the rotation of both the core and the surface. And I wonder too what happens at the poles where the mantle is less impacted by the relative surface speeds of the core and the surface.
Being old and curious about nature is frustrating given how little time is left to discover, but it’s still better than just being old.
Very interesting work and glad to see Leroux’s work in good place!
I have previously criticized the Le Moeul et al. paper and found reasons to reject it. First, instead of studying a proxy for the zonal winds, one should investigate the zonal wind directly. Second, the data is heavily smoothed [4-year sliding window which reduces the number of degrees of freedom enormously and makes the data points very dependent on each other]. Third, the cosmic ray data compared with are not correct. They seem to have been manipulated to improve the fit, as previously pointed out. Fourth, most of the time series used in sun-climate research have no spatial dimension, e.g. TSI, the sunspot number, or the galactic ray flux. Fifth, the mumbo-jumbo level is just too high. Sixth, the venom directed at mainstream scientists is inappropriate.
Fourth, most of the time series used in sun-climate research have no spatial dimension, e.g. TSI, the sunspot number, or the galactic ray flux. Fifth, the mumbo-jumbo level is just too high. Sixth, the venom directed at mainstream scientists is inappropriate.
These points are of course for the present posting.
So, if your spatialtemporal window on the data isn’t wide enough (or you purposely narrow it’s focus) then the results will be grainy and ignorant of the cyclic and/or harmonic nature of the data.
Manns Nature Trick purposely tacked on a narrow window of different data to the tree ring data. Had he left the recent tree ring data attached and run the full length of instrument data, it would not have been misleading. The two do NOT correlate so well, and that much would be instantly obvious. Had he averaged the overlap, he would have committed Simpsons Paradox.
I know that the length of day varies with wind and changing ocean currents and it is getting longer due to the moon, I also know that the variation in length of day is very small(fractions of a second).Do these very small changes alter the temperature of the earth?Some claim that the Sun alters wind patterns such as the jet stream if it does then it is obvious that it will have some correlation with the length of day.How do you fit solar cycles so precisely to the length of day ?(the lod measurements must have error bars).
Tim Channon, if you’re around & you have time, I think some people might be interested in seeing a pdf of an expanded model (one that doesn’t exclude any of the higher amplitude terms). I imagine people would also be interested in hearing the story of the synthesizer, its development, why it was developed, etc., if you are inclined to share.
Best Regards,
Paul.
I have a few questions for anyone who is willing to answer (to see if people are absorbing the message):
1) How would you crudely go about isolating the semi-annual + annual wave using nothing but simple (boxcar) averaging?
rbateman says, “So, if your spatialtemporal window on the data isn’t wide enough (or you purposely narrow it’s focus) then the results will be grainy and ignorant of the cyclic and/or harmonic nature of the data.”
Using the “wrong” extent doesn’t eliminate the cycles, but it can introduce phase reversals (& change the amplitude).
Leif Svalgaard, wrote [of Le Mouël, Blanter, Shnirman, & Courtillot (2010)], “Second, the data is heavily smoothed [4-year sliding window which reduces the number of degrees of freedom enormously and makes the data points very dependent on each other].”
As I’ve explained in a recent WUWT thread, I don’t think they chose their extent with proper cognizance of the properties of the time series. They did, however, get “close enough” to optimal to harness the pattern.
I’ve varied the wavelet extent parameter and found that there’s a fairly wide range around the optimum that will capture the pattern. Base on intuition, this is in absolutely no way surprising.
dp “Being old and curious about nature is frustrating given how little time is left to discover, but it’s still better than just being old.”
Love it.
Look on the bright side. One has more time and there is less distraction. You do really know what you want to do and you can get into it at 4.30 in the morning if you wish and snooze in the afternoon to recover. I like being old. And I like the feeling of focus and patience that comes after that afternoon nap.
It is not surprising that the SOI varies with some aspect of length of day. SOI relates directly to the differential pressure driving the major wind systems which impact the speed of rotation of the Earth and are the essence of the climate system as it evolves day by day basis.
But what does the following mean?
‘When studying the preceding graph, it is important to understand that the blue line is the normalized interannual average of the black line’.
Is this simply de-seasonalised data obtained by averaging over 12 months , or anomalies with respect to the monthly mean or what of?
The real interest is in what lies behind the SOI. I don’t think this analysis is going to help with that question. That’s the big question in my mind.
Thanks Paul,
Well and concisely put together, trying to make progress into the mechanics behind the drivers of climate, would be harder with out a good understanding of these basic Solar, Terrestrial, & Lunisolar harmonic interactions, and their further application to the longer periods of climate periodicities of oscillations.
You do the world a service in helping others to gain a good understanding of how the rest of the interplanetary harmonic interactions, that drive the decadal and longer period cycles, and to extend the thinking to the similarly inter-modulated Inner planet, solar, and moon harmonic period.
As opposed to the longer period cycles of synod conjunction of the outer planets, which jointly swings the sun and the inner planetary system, around the SSB which beats the two long period set together, the inner planets at a rate of once every 17.95 Years and the outer planets periodicity of ~179.5 years.
Leif Svalgaard wrote [of Le Mouël, Blanter, Shnirman, & Courtillot (2010)],“First, instead of studying a proxy for the zonal winds, one should investigate the zonal wind directly.”
…which is exactly why I extended their analysis soon after first reading it:
Vaughan, P.L. (2010). Semi-annual solar-terrestrial power.
http://wattsupwiththat.com/2010/12/23/confirmation-of-solar-forcing-of-the-semi-annual-variation-of-length-of-day/
Leif Svalgaard wrote [of Le Mouël, Blanter, Shnirman, & Courtillot (2010)], “Third, the cosmic ray data compared with are not correct. They seem to have been manipulated to improve the fit, as previously pointed out.”
NO manipulations have been performed on the cosmic ray data presented here:
Vaughan, P.L. (2010). Semi-annual solar-terrestrial power.
http://wattsupwiththat.com/2010/12/23/confirmation-of-solar-forcing-of-the-semi-annual-variation-of-length-of-day/
Direct links to the CR graphs:
1) http://wattsupwiththat.files.wordpress.com/2010/12/vaughn_lod_fig1a.png
2) http://wattsupwiththat.files.wordpress.com/2010/12/vaughn_lod_fig1b.png
Leif Svalgaard wrote, “Fourth, most of the time series used in sun-climate research have no spatial dimension, e.g. TSI, the sunspot number, or the galactic ray flux.”
EARTH is not uniform. See for example Trenberth, Stepaniak, & Smith (2005).
Paul, I am fascinated by your work but I can never understand it well enough to reach the point where I can begin to cross-check it, except to note that it looks as if there really are significant harmonic resonance patterns. If you want to avoid “corrupt science” taking it over, it would help if you made your work accessible to ordinary intelligent generalists (like Willis E). I think if you did this, your dark desire to rant against scientists would actually mostly evaporate.
Focus on straightforward communication – please!
The music of the spheres ?… including the harmonics? Did those old guys in their stone-age observatories come to the same conclusions?
BTW … is this what Leif means by ‘venom’?
“Multimoment multiscale spatiotemporal integration reveals nonrandom harmonic pattern-summary discontinuities, exposing the comedy tragically advocated by deceitful &/or naive theoreticians who are in part constrained by a dominant culture that clings seemingly religiously to maladaptive traditions such as unjustifiable assumptions of randomness, independence, uniformity, linearity, etc. that are routinely misapplied (for example to conveniently render abstract conceptions mathematically tractable).”
Richard Holle says:
April 10, 2011 at 11:29 pm
the inner planets at a rate of once every 17.95 Years and the outer planets periodicity of ~179.5 years.
Not technically correct, Jose when looking at the data only went back several hundred years. The 4 outer planets never really come back to the same position, but they roughly come back to a similar position on average about every 172 years, or the rough synodic period of Uranus & Neptune.
Geoff Sharp says:
April 11, 2011 at 1:38 am
~
Question for Geoff, bit off the topic here or not..
The Earth has a roid, asteroid that is orbiting in its zone.
They are saying that this roid has a “horseshoe” type orbit. Now how can that be, Geoff? Check out the orbit of this thing Geoff.
Astronomers Find Newly Discovered Asteroid is Earth’s Companion
http://www.arm.ac.uk/press/2011/aac_horseshoe_orbit.html
I know you like it when I show you things like this. . or not. .
One more thing if you know Geoff. It is particularly interesting to me, the way that the planetary magnetic dipoles are laid out from the sun to interstellar space. Do you know how closely we monitor the magnetic dipoles of the other planets? Or have any information about this? Saturn appears to be a far enough out of solar influence to vary more in its magnetic dipole locations. A cutoff point..
Now related to the topic or not, was wondering how GCR or ACR could stall vertical currents in Earths upper atmosphere. This would lower the earth electric potential thereby affecting…….
Carla says:
April 11, 2011 at 5:37 am
Carla, the queen of left field 🙂
Asteroids are numerous but have no real effect on planetary systems in this current age other than posing potential threats to life. The 4 outer planets only, have the clout to make solar system changes on big scales in my way of thinking, but leaving the door open re the inner planets on solar cycle length.
On the magnetic front, I see no evidence of any planetary magnetic feedback. The Tsunami like solar wind obliterates all.
What I’ve said:
“Regrettably, some mainstream climate science leaders fell victim to Simpson’s Paradox decades ago; the point is not to issue blame, but rather to suggest that we help pick up the pieces now to enable more efficient data exploration moving forward.” http://wattsupwiththat.com/2010/12/23/confirmation-of-solar-forcing-of-the-semi-annual-variation-of-length-of-day/
Leif Svalgaard wrote, “Sixth, the venom directed at mainstream scientists is inappropriate.”
What is inappropriate is your motion to reject a landmark conceptual finding that opens key doors for generations of researchers. Stalling the authors for a few hours with your editing demands would in no meaningful way affect their core empirical finding.
I would be very curious to see the authors respond to your criticisms here.