Guest Post by Willis Eschenbach
In my last post, which was about the Mauna Loa Observatory (MLO) in Hawaii, Dr. Richard Keen and others noted that for a good comparison, there was a need to remove the variations due to El Nino. Dr. Keen said that he uses the Multivariate ENSO Index (MEI) for such removal.
And what is the MEI when it is at home? Here’s the description from NOAA:
El Niño/Southern Oscillation (ENSO) is the most important coupled ocean-atmosphere phenomenon to cause global climate variability on interannual time scales. Here we attempt to monitor ENSO by basing the Multivariate ENSO Index (MEI) on the six main observed variables over the tropical Pacific. These six variables are:
sea-level pressure (P),
zonal (U) and
meridional (V) components of the surface wind,
sea surface temperature (S),
surface air temperature (A),
and total cloudiness fraction of the sky (C).
Now me, I’m a bit wary of the MEI, because of the possibility of it sharing a variable with something that I’m investigating. For example, in the post I did on the Mauna Loa Observatory in Hawaii, cloudiness was a variable because I was looking at solar radiation. However, I used it because Dr. Keen used it, and because for the purposes of my post it turned out the considerations didn’t matter.
The effects of the El Nino don’t happen immediately, of course. In general, the effect of the El Nino on the global temperature lags the El Nino by a couple of months. You can determine the lag by using a “cross-correlation analysis”, which shows the correlation between the El Nino and the variable of interest over a wide range of lags.
So imagine my surprise when I did the cross-correlation between the MEI and the Mauna Loa Observatory temperature and got the following result:
Figure 1. Cross-correlation between Mauna Loa Observatory (MLO) temperatures and the Multivariate ENSO Index (MEI).
Zowie … I was expecting a two or three-month lag, but the peak correlation is not lagged a couple months, a couple of quarters, or even a full year. Instead, peak correlation is at no less than a fifteen-month lag.
Now the joy of science is in the surprises. When I get surprised, I don’t sleep right until I learn more about what it was that surprised me. I couldn’t figure out how it was that Hawaii got basically no correlation for six months, and then after that, the correlation kept increasing until it peaked at fifteen months.
So I made an investigation of the correlation of the MEI with the individual 1° latitude x 1° longitude gridcells of the planetary surface. As you might imagine, at a lag of one month you have the strongest correlation between the MEI and the tropical Pacific. Here’s that map:
Figure 2. Correlation, MEI and 1°x1° gridcells. The dark blue lines outline the areas where the correlation is less than minus 0.2. Red outlines the areas where the correlation is greater than plus 0.2.
You can see that the area of the central Equatorial Pacific has the highest correlation with the MEI. The light blue rectangle shows the NINO 3.4 area, which is used in the same way as the MEI is used, to diagnose the state of the El Nino. So the high correlation there makes sense.
Figure 2 also reveals why the correlation with Hawaii is so low at the one-month lag. It is because Hawaii (black dot above the left side of the light blue rectangle) is very near the edge between the red and the blue areas, where the correlation is small.
To investigate the longer lags, I decided to make a movie so I could understand the evolution of the El Nino variations as they spread out and affected other parts of the world. Here’s that movie. It shows the correlation of the MEI and the individual gridcells at periods from one month to 24 months and then back down again to one month.

Again, more surprises. The correlation dies away quickly in some areas, but in Hawaii it builds until about fifteen months, and then decreases after that.
.
How amazing is that? If you want to know what the temperatures at the Mauna Loa Observatory will be doing fifteen months from now, you can look at the MEI today.
To demonstrate this odd fact, here are the MLO temperatures compared to the Multivariate ENSO Index lagged by 15 months.
And that’s why I love climate science … because there is so much to learn about it.
[UPDATE] I got to thinking that I should do the same thing using the NINO34 Index … here is that result. As you can see, it is extremely similar to the MEI graph above.

Here in the forest, after doing most of the work in this post, I was moving my computer files into new computer folders yesterday evening and I managed to destroy about half of them … and my last backup was three weeks ago.
So I was up until 2 AM saying bad words and reconstructing lost functions that went into the bit bucket. Grrrr … then I spent all day today beating my files back into submission so that I could recreate the work that I’d already done on this post. However, it allowed me to clean up some poorly written functions, and I suppose anything worth doing is worth doing twice.
It was also another demonstration of my rule of thumb gained from living about 20 years on tropical Pacific islands, which is:
The Universe doesn’t really give a damn what I think should happen next.
Factcheck: True … the good news is in the corollary to that rule of thumb, which is:
I do have a choice in the matter: I can dig it or whine about it.
Endeavoring to follow my own maxims re my monumental computer blunder, I remain,
Yr. Obdt. Svt.,
w.
I Know This Gets Old But: When you comment, please quote the exact words you are discussing. Misunderstandings abound on the intarwebs. Clarity about your subject can minimize them.



This is worth writing a paper about.
Thanks, Anthony. I was quite impressed with the accuracy of the 15-month prediction … perhaps someone with a PhD would like to assist me with writing it up. My problem is that my writing style is so far from the desired journal style that I need to get out an icepick and give myself a lobotomy to write that way.
If someone wants to assist me to do the writing and submission, I’m glad to give them first author status.
w.
This paper gives an ~8 month lag:
http://iopscience.iop.org/article/10.1088/1748-9326/aa9c5b#erlaa9c5bs4
Willis,
NEVER change your style! You are WUWT’s top communicator and are probably having some positive effect on others. Then your thinking “outside the box” is …. well …. outside the box and you can never see outside the box unless you look there. Keep up the good work!
Willis: your writing style makes it easy for the reader to understand why you started the project, see your data, see where it comes from, follow your arguments and understand how you get to your conclusions. This makes it better than 90 percent of the scientific literature out there. (at least). In climate science the number is more like 99 percent.
Another thing, reading a Willis post never makes me feel drowsy and/or sends my eyes drifting to clickbait around the edges of the screen. Can’t say that about a lot of scientific literature, even in my own field of economic geology, where a lot of what I read is familiar, and relevant to my work.
You may be right, you may be wrong, but by golly, there’s never any doubt about what you’re saying. Please don’t try and sound more “scientific”.
Willis,
Your ‘telling a story as it unfolded’ narrative style is a very effective communications technique! It ‘leads’ the reader down the investigation path, highlights why certain branches on the path were followed and others passed by, and develops a deeper understanding and appreciation of the results thereby in the target audience.
Please don’t change it….
J Mac
I don’t know; like every good writer, Willis is cognizant of his audience. Most rarely in human beings, he is also cognizant of his limitations, Speaking for myself, I could use a dose of that self-awareness from time to time. Kudos to Willis for reaching out for assistance.
Willis:
Extend all numbers to three decimal places.
Blame everything on global warming.
Claim your data are proof
of a coming climate catastrophe
Submit to any scientific journal for publication.
This is how “modern climate science” works.
My climate change blog:
http://www.elOnionBloggle.Blogspot.com
Read this tomorrow, I am obviously unable to understand it tonight.
I saw the head as MLO and ME! Expecting pictures of Willis on the mountain…
AWESOME! The new WordPress has an edit function for the first five minutes. Brilliant.
w.
Willis, here’s wishing you the best henceforth at the new, new wuwt! May this be the first of many enjoyable posts for you (especially so in the comment sections)…
Edit!! Yes!!!! Noticed that yesterday. Haven’t tried it, but the concept is wonderful.
Edit facility is only good for those who take the trouble to re-read their comment. My (despondent) bet is that many clangers will remain because those making them don’t have a clue about good grammar, manners, or sometimes, even facts.
Back on topic, another excellent dissection of the data, Willis, & I also love your writing style.
Willis, this is yet another fascinating piece of work, thank you. It seems to me that it is also troubling from the perspective of the old chestnut of trying to arrive at a meaningful global average surface temperature and/or the same expressed as an anomaly. There seem to be strong variations even within relatively small areas of the globe, whereas you might have expected the ENSO effects to be more consistent over larger areas, especially in the oceans. This, it seems to me, highlights the importance of comprehensiveness of data coverage and thus the degree to which we can rely on past figures – or indeed not!
Anyhow, a related question: whilst I appreciate that these maps show correlations with a calculated index, rather than temperatures, is it possible to infer what the difference between a dark red area and a light blue one represents in terms of temperature?
I would love to see an overlay of the shipping lanes on this map. It would give an indication how ‘accurate’ the historic sea surface temperatures might have been, or how sensitive they were to this phenomenon.
Is the lag simply due to the speed of the north pacific current? It moves at about 0.05 m/s and would take
about 15 months to move the warm water from the western pacific to Hawaii.
Willis, how might this feed into your mystery question of why MLO temps didn’t dip during el chichon & pinatubo?
No clue, fonz. It’s a separate question. I already adjusted for the effect in my post on el chichon & pinatubo.
w.
Is the lag not just due to the speed of the north pacific current? It would take about that long for
the warm water to move from the western pacific to Hawaii
Willis……thanks, as always, for posting your analyses here with such candour, humility and sense of wonderment. We are all enriched and better-informed thanks to your efforts and I, for one, am really intrigued by the stuff you have been researching in this series of recent posts. Absolutely fascinating stuff and true scientific inquiry. Brilliant! Thanks!
Great work Willis. Your mastery of R for pulling out this sort of graphical representation of the data is impressive.
From the animation it seems that Hawaii is in an island of positive correlation throughout the 24mo lag range. This does not seem consistent with MLO being negative up to six months. Is this just a different between air temp at MLO ( top of a mountain ) and the 1 degree cells predominantly SST (?) around Hawaii.
What was the data you were using? You do not seem to state anything more precise than ”
gridcells of the planetary surface”. Good practice would suggest you should at least state what dataset you are using here. An oversight you may like to correct in a note.
It’s also interesting that continental land areas seem to be weakly negatively correlated throughout.
Interesting as always. Thanks.
Greg June 1, 2018 at 11:57 pm
Thanks, Greg. The pictures are what make sense of the numbers.
I suspect that’s because there is slightly different data being used in Figure 1 and the movie. See below.
True, and good practice would suggest that I shouldn’t have erased a good chunk of my computer programs, so I’d have more time to write and spend less time programming.
The movie is made entirely by comparing gridcell data from the CERES satellite dataset with the MEI. Figure 1, on the other hand, shows the correlation between the MEI and the Mauna Loa Observatory dataset.
I think that the answer to that was provided by Percy Jackson in this thread, who said:
I believe Percy has hit on the explanation, that we’re seeing the movement of water masses in response to the El Nino.
My pleasure.
w
Willis,
It’s a very interesting correlation plot and movie. Is it temperature that you are correlating, or anomaly? If temperature, then seasonality is an issue. MEI doesn’t have much, but MLO does, and ENSO correlates with seasons.
The lagged correlations of MEI with itself are interesting too.
Interesting Nick, do you have an ACF of MEI to point to?
Greg,
No, I was just looking at the region in Willis’ movie.
OK that’s CCF with “surface” , not sure exactly what that is yet.
Apparently there is a 60 mo periodicity in MEI but I have not found a non paywalled copy of the paper claiming that.
I see the later graphs of MLO are anomaly, so I guess that is the answer. So no seasonality problem.
Nick Stokes June 2, 2018 at 12:07 am
Good question. I’ve removed the seasonality from both the MLO data and the CERES gridcell data before doing the analysis.
w.
If El Nuno events are supposed to have an effect on world temperatures then why are there areas where little or no correlation exists. Have I missed something basic here?
I have been trying to read as much as possible about climate and climate change so perhaps all of those facts and figures bouncing around inside my small mind have become so scrambled that all that I have managed to do is confuse the hell out of myself.
I know exactly how you feel Crowcane!
Seems like the more you find out, the less you know and there are far more questions than before you started your journey.
Also makes you even more cynical about politicians and bureaucrats as they cannot possibly understand the basics of climate science yet act and sound with confidence as if they really do.
PS El Nuno – sounds better!
You have a temperature, MLO. If you want to correct for El Nino, why not use another temperature, Nino 3.4 SST?
I wonder what is the correlation between Nino 3.4 SST and MEI.
Thanks, Bob. I used the Nino34 anomaly for the final plot. There’s good correlation between the NINO34 index and the MEI, hang on … correlation = 0.90.
w.
I missed that.
Given the lag I assume, as do several others, that the answer is ocean currents. As I googled the currents I was struck by the differences between different images. link link The second link indicates that we’re still learning new stuff about ocean currents. Willis, I wonder if your observation will promote a better understanding of the currents involved.
Hi Willis,
Slightly different but closely related:
I get only a 4-month lag of Global LT temperature after the Nino34 anomaly. Here is a previous post.
Regards, Allan
Post Script:
My plots used to display – at least for me – in the old server. Now they do not – only my Facebook reference appears. Can the rest of you at least access the plots via Facebook?
https://wattsupwiththat.com/2018/05/clouds-and-el-nino/#comment-2349718
The Nino34 Area Sea Surface Temperature (the blue line in the upper plot), adjusted by the Sato Global Mean Optical Depth Index (for major volcanoes – the yellow line), correlates quite well with the Global UAH LT temperature four months later (the red line).
https://www.facebook.com/photo.php?fbid=1527601687317388&set=a.1012901982120697.1073741826.100002027142240&type=3&theater
UAH LT temperature also correlates well with Equatorial Atmospheric Water Vapor of one month earlier (the GREEN line in the lower plot). This mechanism ties to Willis’ clouds – down to the hourly level – the rising Sun drives water vapor off the oceans into the atmosphere, clouds form, shade the ocean, and that is the natural regulator of temperature.
https://www.facebook.com/photo.php?fbid=1665255773551978&set=a.1012901982120697.1073741826.100002027142240&type=3&theater
This close correlation of Nino34 SST’s with global temperatures four months later either means that this small Nino34 area of the Pacific Ocean drives global temperature or other tropical oceanic areas have a similar temperature signature – that they correlate with Nino34 SST’s.
How does CO2 play into this equation? Atmospheric CO2 is increasing at a “base rate” of about 2ppm/year, probably due to fossil fuel combustion, deforestation, etc. However, the rate of change dCO2/dt changes ~contemporaneously with average global temperature, and its integral, the trend of atm. CO2 changes ~9 months after global temperature, so it is clear that global temperature drives atm. CO2 much more than atm. CO2 drives temperature, and so the sensitivity of temperature to atm. CO2 must be relatively small, and far too small for a real global warming crisis to exist.
It appears to be just about that simple! And the very-scary global warming climate models do not reflect this mechanism at all. I suggest we can do a much better job of modelling climate with the simple closed-form solutions presented in these plots, maybe with a bit more refinement – and they will show NO GLOBAL WARMING CRISIS.
More details here
https://wattsupwiththat.com/2018/05/is-climate-alarmist-consensus-about-to-shatter/#comment-2343879
and here
https://wattsupwiththat.com/2018/04/solar-activity-flatlines-weakest-solar-cycle-in-200-years/#comment-2341336
Are you sure the files you lost haven’t been inadvertently moved somewhere else? Windows Explorer (if that is what you were using) is notorious for doing this kind of thing. I avoid using it. There are better alternatives like Altap Salamander.
A search for a filename known to be missing might throw some light.
Ian Macdonald June 2, 2018 at 1:10 am
Thanks for the good thoughts, Ian, but sadly, what I did was to overwrite current versions with archived older versions …
Grrrr …
w.
UAH Global Temperature Update for May, 2018: +0.18 deg. C [for LT = Lower Troposphere]
http://www.drroyspencer.com/2018/06/uah-global-temperature-update-for-may-2018-0-18-deg-c/
Thank you again Roy Spencer and John Christy.
My note:
This down slightly from 0.21C in April 2018 but is still 0.18C higher than my prediction of ~4 months ago of 0.0C.
Isn’t it annoying when the planet just refuses to cooperate perfectly with one’s prediction? 🙂
My longer post made last night is still in the spam filter, so I am posting this one.
Hi Willis,
Slightly different but closely related:
I get only a FOUR-month lag of Global LT temperature after the Nino34 anomaly. Here is a previous post from May 2018 – my other similar posts date back to 2016.
Regards, Allan
https://wattsupwiththat.com/2018/05/clouds-and-el-nino/#comment-2349674
Really good Willis. Here is some stuff posted on wattsup over previous months/years that may help tie it all together.
Best, Allan
https://www.facebook.com/photo.php?fbid=1527601687317388&set=a.1012901982120697.1073741826.100002027142240&type=3&theater
https://www.facebook.com/photo.php?fbid=1665255773551978&set=a.1012901982120697.1073741826.100002027142240&type=3&theater
The correct mechanism is described as follows (approx.):
Equatorial Pacific Sea Surface Temperature up –> Equatorial Atmospheric Water Vapor up 3 months later –> Equatorial Temperature up -> Global Temperature up one month later -> Global Atmospheric dCO2/dt up (contemporaneous with Global Temperature) -> Atmospheric CO2 trends up 9 months later
What drives Equatorial Pacific Sea Surface Temperature? In sub-decadal timeframes, El Nino and La Nina (ENSO); longer term, probably [edit: some function of] the Integral of Solar Activity.
The base CO2 increase of ~2ppm/year could have many causes, including fossil fuel combustion, deforestation, etc, but it has a minor or insignificant impact on global temperatures.
Hi Allan, re slight decline on May temp – I’ve commented back when Tisdale was writing about ENSO during the unusual build-up of the last El Nino, that I felt Global SSTs were beginning to “depart” from a close relationship with ENSO. We had the stationary ‘hot blob’ off the West Coast of NA and other hot blobs in the temperate zones that appear to have modulated the relationship.
Today, with ENSO moving toward threshold El Nino we have large ‘cold blobs’ N and S of the equatorial zones. Instead of the usual strong upwelling of cold water in the eastern eq area, the equatorial band is being fed by cold water slanting equatorward from these large cold masses. I believe we will continue to cool even with ENSO venturing into warming. The warm water moving eastward at depth seems to be having cold water piled over it and the western warm pool has similarly been invaded by cold water. It is happening in the atlantic as well.
Mr Eschenbach, that is very impressive programming and the finding of the time lag and how it varies world wide demonstrates just how complex the Earth’s Climate System is.
It also demonstrates that the Oceans rule the Atmosphere.
Great post as always, Willis! I just love your way of thinking! Which triggers some questions:
I wonder how the movie looks from the opposite side of the world, centered at the prime meridian and going on even a little longer. In the North Atlantic the area of positive correlation seems to be still growing at 24 months lag. (It is non-existent at 12 months lag.) Which part of the world has the longest lag? Is there a an area with a 3 year lag?
What about the opposite direction? Which areas of the planet can indicate the development of the NINO34 area one, two or even three years in advance?
Is there an anti-NINO34 area on the planet? An area that when cooling has a lagged correlation to many other parts of the world?
Or more broadly: which areas, cooling or warming, have the largest lagged correlation with other parts of the planet?
If NINO34 has by far the largest impact, then it is even more likely that your governor theory is describing the most important mechanism for a stable climate on an ocean planet.
Thanks again for all the interesting posts and all the beautiful colors…
Re. The article. That’s neat. Predictions that actually work. IMHO, Climate research needs a whole lot more of that and a whole lot less Old Time Religion pretending to be “science”.
Re: Data loss. If it’s a computer, you are going to lose/overwrite files every now and then. And yes, it is super-annoying at best. If you are using Unix or a Mac, you might want to look into Rsync. It’s not easy to configure, but once set up, you can do backups every day if you want to without a lot of hassle. Personally, I do them every other day. A Windows equivalent seems to be Deltacopy — which I have no experience at all with.
Don, the Mac makes it simple because it has the Time Machine built into it. Unfortunately, I moved my office three weeks ago and neglected to plug my Time Machine backup disk into the system, so my a good chunk of my last three weeks work was lost.
w.
If there is a multi-month lag between El Nino releasing heat into the atmosphere, why is there no lag between ENSO and interannual variations in length of day?
http://lh5.googleusercontent.com/-Z7oqT0Wm3HE/VI11FCPzRSI/AAAAAAAAAj8/ZSBUiI0Iexs/s800/LOD-vs-ENSO_II.png
The effect of EL Nino on the inter-annual LOD signal is supposedly caused by conservation of angular momentum, when the atmosphere inflates due to warming. So why is there no lag in LOD if there is a delay in the release of heat from the ocean surface?
Why does LOD change immediately?
I believe the seawater expands contributing to delta LOD more quickly
While weather smoothly crosses the International Date Line, weather forecasts are dated so they cannot. A crust of ENSO forecast model is always building up behind the dateline and successive runs become more densely packed, piling up with increasing friction. A Global Calastrophic Event occurs when those with 12 month calendars are unable to reconcile with others’ 18 month calendars. Your 15 month lag is the gestation period when 12 month calendars breed with 18s, and it is no coincidence that when the remainder is paired with a human gestation period there are two full cycles. People are often conceived when there is nothing written on one’s calendar.
Are there any Hawaiian weather forecasters out there who were aware of this lagged correlation? If they weren’t aware of it then it seems they might owe Willis a drink.
“In general, the effect of the El Nino on the global temperature lags the El Nino by a couple of months.”
An AMO lag of ~8 months. E.g. peaking warm in August 1998, August 2010, and August-September 2016. That will effect CO2 uptake rates in the North Atlantic, and drive regional precipitation changes, effecting land CO2 uptake rates, like El Nino does to the Amazon region.
https://www.esrl.noaa.gov/psd/data/correlation/amon.us.data
I published an article “Forecasting the onset of southwest monsoon over Kerala”. In this I presented dates of onset collected by me from Indian Daily Weather, Reports, weekly weather reports and monthly weather reports published by Indian Meteorological Department. for 1921 to 1975 These were extended backwards upto 1921 to 1901 using the method developed by us for answering a query from the Indian Parliament on Delhi normal date of onset.
Also collected the Singapore 50 mb level zonal winds and compared with the dates of onset over Kerala for 1964 to 1975. When the winds were easterlies at 50 mb in the month of May, the onset will be late and when the winds are westerlies the onset will be early. The winds at 50 mb are very systematic: westerly regime, it is always 12 and 24 months but easterly regime is 12 and 6 months.
In fact I collected the data for Singapore, Bogota, Malakal and Lima for 100, 50 and 30 mb. Except 50 mb of Singapore, none of them showed good relation with the onset date. Few decades later some modelling groups from IITM,Pune and Ahemdabad tried to through dust on the work. They used estimated dates of onset and compared with other stations data [which I did not get any relation] and said my paper is not correct. When I sent my counter, the journal did not published my version as the editorial board is under their control. These are the people involved in IPCC.
Dr. S. Jeevananda Reddy
Very interesting. Like Anthony noted, it is important enough to write a paper about – not so much just MLO, but the global / regional lag correlations witht the ENSO.
Note that there is a Pacific Ocean circulation system and some of the ENSO circulation leaks north into the Kuroshio current which would then change MLO temperatures as the westward winds move energy into the atmosphere. It takes time for warm El Niño waters to move all the way to the North Pacific through the Kuroshio. The PDO index and the circulation it represents is an example of how these lags develop.
El Nino suppresses CO2 outgassing from the ocean, and is followed by ‘a larger (and lagged) response from the terrestrial component’.
(No mention of a lagged CO2 uptake reduction in the North Atlantic though)
https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/20180000620.pdf