ENSO Dominates NODC Ocean Heat Content Data

ENSO Dominates NODC Ocean Heat Content Data

Guest post by Bob Tisdale, BTW here is the current SST map. – Anthony

clickable global map of SST anomalies
Click for a larger image

The Royal Netherlands Meteorological Institute (KNMI) recently added the National Oceanographic Data Center (NODC) Ocean Heat Content (OHC) dataset to their Climate Explorer website, allowing users to download data based on user-defined global coordinates.

http://climexp.knmi.nl/selectfield_obs.cgi?someone@somewhere

This OHC dataset was presented in the Levitus et al (2009) paper “Global ocean heat content(1955-2008) in light of recent instrumentation problems” [GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L07608, doi:10.1029/2008GL037155, 2009]

ftp://ftp.nodc.noaa.gov/pub/data.nodc/woa/PUBLICATIONS/grlheat08.pdf

There are differences in the presentation of the data. The NODC illustrates their OHC data in 10^22 Joules, but KNMI presents the data on an area-averaged basis, in units of Gigajoules (10^9 Joules) per square meter. The data is the same; the units in which the data is presented are different. Also, the NODC provides the data on a quarterly basis; that is, the data is grouped in three-month averages. KNMI presents the NODC OHC data on a monthly basis by listing the quarterly data for each of the three months. This is why the OHC data appears to be squared off in the graphs of monthly raw data. This can be seen in Figure 1.

http://i32.tinypic.com/29de5ow.png

Figure 1

Figure 1 is a comparison graph of the Global OHC anomaly data (NODC), scaled NINO3.4 SST Anomalies (HADISST), and scaled Sato Index (GISS) data. This is the same format used in the graphs of the subsets illustrated in this post. The NINO3.4 SST anomalies are used to illustrate the timing of the El Nino-Southern Oscillation (ENSO) events. The Sato Index of Mean Optical Thickness at 500nm are provided to illustrate the timing of explosive volcanic eruptions. I’ve also smoothed the data for each OHC anomaly subset with a 13-month running-average filter, Figure 2. As you will see later, some of the subsets are noisy in their raw form.

http://i32.tinypic.com/sno57l.png

Figure 2

In the following, I’ve provided links to the graphs of the raw data, for those who are interested in seeing it in that form, but I have only posted the graphs of the data smoothed with a 13-month running-average filter. It’s much easier to see the step changes when the data is in that form.

TROPICS

The Tropical Pacific OHC anomaly data is illustrated in Figure 3. A number of things to note: The tropical Pacific OHC anomalies fall during El Nino events, but then recharge during the La Nina. For the most part, when the El Nino events occur at the same time as volcanic eruptions, the recharge does not return the OHC anomalies to the value they were at before the El Nino, but if the El Nino occurs without the influence of a volcanic eruption, the La Nina recharges the Tropical Pacific OHC anomalies to the pre-El Nino level. And it does it quickly. Note also how the 1972/73 El Nino event causes an upward step in the OHC anomalies of the Tropical Pacific. The OHC anomalies then decrease gradually, being influenced by the eruptions of El Chichon in 1982 and Mount Pinatubo in 1991, until they rise suddenly in 1995. In an earlier post, I illustrated how a shift in Tropical Pacific Total Cloud Amount may have caused the 1995 rise in Tropical Pacific OHC, providing fuel for the 1997/98 El Nino. Refer to my post Did A Decrease In Total Cloud Amount Fuel The 1997/98 El Nino?

http://i25.tinypic.com/wrz71x.png

Figure 3

http://i31.tinypic.com/2s96hd1.png

Figure 3 Raw

However, the Tropical Indian Ocean OHC anomaly data reveals a sudden decline in 1995. Did a shift of warm water from the Tropical Indian Ocean to the Tropical Pacific also fuel the 1997/98 El Nino? I’ll investigate this in a future post. Note how the Tropical Indian Ocean OHC anomalies correlate with NINO3.4 SST anomalies over a large portion of the term of the data, but after 1995, the amplitude of the variations changes drastically.

http://i25.tinypic.com/atkaa8.png

Figure 4

http://i30.tinypic.com/xfnk14.png

Figure 4 Raw

In Figure 5, I’ve combined the OHC anomaly data for the Tropical Indian and Pacific Oceans. The OHC anomaly data for this subset follows the base of the NINO3.4 SST anomalies remarkably well. The OHC anomalies of the Tropical Indian and Pacific Oceans follow the rise in NINO3.4 SST anomalies after the 1972/73 and 1997/98 El Nino events. In other words, like the Tropical Pacific, there also appears to be a 25-year decay after the upward step from the 1972/73 El Nino (also influenced by the 1982 and 1991 volcanic eruptions), until the 1997/98 El Nino causes another upward step.

http://i26.tinypic.com/2j60dfp.png

Figure 5

http://i28.tinypic.com/2a3y2a.png

Figure 5 Raw

The step changes in the Tropical Atlantic OHC anomalies are obvious. The first occurred three years after the peak of 1972/73 El Nino, as the NINO3.4 SST anomalies rose from the secondary minimum of the two-year La Nina event. The same thing occurred with the next significant El Nino that was strong enough to generate a La Nina that lasted through two ENSO seasons, and that was the 1997/98 El Nino. Note also how the OHC anomalies of the Tropical Atlantic have been dropping quickly since 2005. Click on the link to the raw data (Figure 6 Raw) to see just how precipitous that drop has been in recent years.

http://i28.tinypic.com/1jnp87.png

Figure 6

http://i28.tinypic.com/2a3y2a.jpg

Figure 6 Raw

MID-TO-HIGH LATITUDES

The North Pacific OHC anomalies are like no other OHC subset. In 1967, there was a sudden drop in the North Pacific OHC anomalies. Twenty plus years later North Pacific OHC anomalies rebounded. I’ll have to investigate this dataset further in a later post, to try to isolate where the majority of that variability takes place.

http://i28.tinypic.com/f56pfm.png

Figure 7

http://i29.tinypic.com/rwp8ut.png

Figure 7 Raw

As illustrated in Figure 8, the South Pacific OHC anomalies show a sharp upward step change following the 1997/98 El Nino. Between 1971 and 1996, the OHC anomalies oscillate at or near 0 GJ/sq meter. The cause of the small rise between the 1960s and 1970 is elusive, but it’s not a significant rise compared to the upward step after the 1997/98 El Nino.

http://i26.tinypic.com/xuhkn.png

Figure 8

http://i27.tinypic.com/25s5ta1.png

Figure 8 Raw

The South Indian Ocean OHC anomaly data, Figure 9, shows a decrease from 1955 until the late 1960s. Then the 1968/69/70 El Nino caused a minor rise in OHC anomalies. This was followed by a major upward step from the 1972/73 El Nino. OHC anomalies in the South Indian Ocean remained relatively flat until the eruption of Mount Pinatubo, when the OHC anomalies dipped. The upward step change after the 1997/98 El Nino is hard to miss. The decay until 2006 almost returned the South Indian Ocean OHC anomalies to the pre-1997/98 values, but the El Nino of 2006/07 bumped it back up again.

http://i31.tinypic.com/34jamtj.png

Figure 9

http://i31.tinypic.com/2dqvpfl.png

Figure 9 Raw

The North Atlantic OHC anomaly data, Figure 10, with its gradual climb, is clearly dominated by the Atlantic Multidecadal Oscillation. The impacts of ENSO events are visible, however. In a future post, I may detrend the North Atlantic OHC anomaly data to emphasize the ENSO impacts on this dataset.

http://i25.tinypic.com/2s17wpt.png

Figure 10

http://i32.tinypic.com/swa4xf.png

Figure 10 Raw

There is a clear step change in the South Atlantic OHC anomaly data, Figure 11, following the 1972/73 El Nino. In this case, however, the response appears to be lagged an extra couple of years. The response is so long, it appears to result from the lesser El Nino of 1976/77. The South Atlantic OHC anomalies remain relatively flat until they appear to respond to the 1997/98 El Nino with an upward step that starts again many years after the peak of the El Nino. Why so long?

http://i30.tinypic.com/2n9xsv6.png

Figure 11

http://i25.tinypic.com/2qdcx7l.png

Figure 11 Raw

Could the variations in the South Atlantic OHC anomalies simply be lagged responses to the Tropical Atlantic OHC anomalies, with surface and subsurface currents transporting the waters from the tropics to the mid-to-high latitudes of the South Atlantic? Refer to Figure 12.

http://i28.tinypic.com/2uffyfr.png

Figure 12

ARCTIC AND SOUTHERN OCEANS

I’ve provided the Arctic and Southern Ocean OHC anomaly data in Figures 13 and 14, without commentary, for those who are interested in seeing what those curves look like.

http://i31.tinypic.com/23u23cz.png

Figure 13

http://i28.tinypic.com/wa0tu0.png

Figure 13 Raw

############

http://i28.tinypic.com/53ve2w.png

Figure 14

http://i28.tinypic.com/2niwilg.png

Figure 14 Raw

CLOSING

It is clear that significant El Nino events can and do cause upward step changes in Ocean Heat Content. This indicates that ENSO events do more than simply release heat from the tropical Pacific into the atmosphere. Apparently, El Nino events also cause changes in atmospheric circulation in ways that impact Ocean Heat Content. If and when GCMs are able to recreate the variations in atmospheric circulation that cause these changes in Ocean Heat Content, GCMs may have value in predicting future climate variability. At present, they do not.

SOURCES

The NINO3.4 SST anomaly data is based on HADISST data available through the KNMI Climate Explorer:

http://climexp.knmi.nl/selectfield_obs.cgi?someone@somewhere

Sato Index data is available through GISS:

http://data.giss.nasa.gov/modelforce/strataer/

Specifically:

http://data.giss.nasa.gov/modelforce/strataer/tau_line.txt

Get notified when a new post is published.
Subscribe today!
0 0 votes
Article Rating
72 Comments
Inline Feedbacks
View all comments
September 6, 2009 2:03 pm

hotrod: You wrote, “It seems to me like that sudden cooling in Tropical Indian Ocean heat content, with a double dip first at about 1997, and than the bigger dip in 1999-200, could be best explained by a sudden upwelling of deep cold water in the Indian ocean…”
I’ll look into those dips in the Indian Ocean in a follow up post. I believe we’ll find that those dips are strongest in the eastern tropical Indian Ocean. The second dip in 1997 should be part of the 1997/98 El Nino. That El Nino was so strong it drew warm water from the eastern Indian Ocean. (Remember, the Pacific Warm Pool is also known as the Indo-Pacific Warm Pool because it stretches into the Indian Ocean at times.) Here’s a video if SSH anomalies. Note how the eastern tropical Indian Ocean anomalies drop during the 1997/98 El Nino, about one minute into the video.

And I believe the rise in the western Indian Ocean during the 1997/98 El Nino is a response to the cooling in the east. It would have increased the temperature difference, which in turn would have kicked up the trade winds in the Indian Ocean, which would have pushed more waters to the west, draining more warm water from the east, and adding more positive feedback to the process.
But if the dips in 1995 and 1997 in OHC are located in the western or central tropics of the Indian Ocean, leaving the eastern Indian Ocean untouched, then we’ll have to look for other answers.
We shall see.

September 6, 2009 2:16 pm

lucklucky: You wrote, “How can anyone make a post in disparate geographic, time, instrument quality data AND still call it Science?”
You must be referring to the wiggle matching in my post and not the Levitus et al, HADISST, and Sato Index data. Apparently you only looked at the graphs and didn’t bother to read the post, because you would have noted that I discussed the reasons for including the NINO3.4 SST anomaly and Stratospheric Mean Optical Thickness data. Let me repeat it for you, “The NINO3.4 SST anomalies are used to illustrate the TIMING of the El Nino-Southern Oscillation (ENSO) events. The Sato Index of Mean Optical Thickness at 500nm are provided to illustrate the TIMING of explosive volcanic eruptions.” [Caps on the word TIMING added.] They’re simply time markers.
Have a nice day.

September 6, 2009 2:34 pm

TonyB: You wrote, “I am trying to find out the temperatures of arctic waters (which will vary considerably according to specific location) both close to the surface and under the ice. I want to compare it with historic records. I really need them to be location specfic rather than covering a huge area. Any ideas where to look?”
The HADISST data is presented in 1 deg geographic grids through the KNMI Climate Explorer, but that’s only surface. I have not found (also haven’t looked for) subsurface ocean temperatures at given depths presented on a webpage anywhere, but the NODC had to compile it for the OHC data.
Sorry that I couldn’t be more help.

September 6, 2009 2:40 pm

Paul Vaughan: You wrote, “Bob, if warm Indian water moved to the Pacific in a ‘bubble’ or burst in the 90s, I’ll be curious to see (in a future post) which path(s) you figure it might have taken. Politics & economics aside, there is something fundamental about nature to be learned here.”
As I just noted to hotrod…
The second dip in 1997 should be part of the 1997/98 El Nino. That El Nino was so strong it drew warm water from the eastern Indian Ocean. (Remember, the Pacific Warm Pool is also known as the Indo-Pacific Warm Pool because it stretches into the Indian Ocean at times.) Here’s a video if SSH anomalies. Note how the eastern tropical Indian Ocean anomalies drop during the 1997/98 El Nino, about one minute into the video.

And I believe the rise in the western Indian Ocean during the 1997/98 El Nino is a response to the cooling in the east. It would have increased the temperature difference (between east and west), which in turn would have kicked up the trade winds in the Indian Ocean, which would have pushed more waters to the west, draining more warm water from the east, and adding more positive feedback to the process.
What happened in 1995 is beyond me at the present time. I still haven’t narrowed down the location of the dip. I’ve assumed it was in the eastern Indian Ocean.

September 6, 2009 2:55 pm

matt v.: You wrote, “How difficult would be for you to run Southern Hemisphere Heat Content and HADCRUT3? Also where would I find Southern Pacific SST anomalies table data . I have found a close correlation between Southern Hemisphere SST and Hadcrut3gl”
All of those datasets you seek are on the KNMI Climate Explorer.
http://climexp.knmi.nl/selectfield_obs.cgi?someone@somewhere
It’s easy. Select a dataset. On the next page, enter the coordinates. Next page, scroll down to the third graph (anomalies) and click on “raw data”. Convert the data to a format useable by your spreadsheet.

September 6, 2009 3:15 pm

Lindsay H: You asked, “any theory as to why there is a coincidence with the 13 month smoothing , it seems to show a considerable lag of up to two years to recover from events…”
Sorry, Lindsay, I just found the apparent relationship between significant ENSO events and OHC. No theories yet. I agree that two (plus) years appears to be a long time before an apparent response from some datasets, like the Atlantic. Is it the magnitude of the La Nina that follows the significant El Nino that shifts cloud cover and causes what appears to be a delayed response to the El Nino? That is, are those shifts responses to the significant La Ninas, not the El Ninos themselves? Dunno.

Stephen Wilde
September 6, 2009 3:36 pm

“Roy Spencer (10:34:43)
“warming in the pipeline” refers to a global radiative imbalance which has not yet been relieved by warming because the heat capacity of the system is so large. It exists in the theory of AGW, but has not actually been measured because it is so small (Hansen says somewhat less than 1 W/m2). If the climate system is insensitive, though, then any radiative imbalances are fairly rapidly eliminated by temperature change.
Which brings up a point of confusion…a very sensitive climate system takes longer to equilibrate to a new temperature, but it’s time rate of change of temperature is actually larger than for an insensitive system. Let’s see if anyone here can figure out how that can be true.”
Reply:
If, say, there is a rapid and proportionate response to a temperature forcing then that would be an indication of a system that is sensitive to forcings.
If the system is sensitive in that way then it will delay a temperature rise by immediately and proportionately initiating a neutralisation process and in due course cancel that temperature rise so in that separate sense the system is insensitive to such forcings. It would take a very powerful forcing to make any significant temperature difference for any length of time.
That is exactly how a change in the speed of the hydrological cycle would work and it operates through changes in the latitudinal position of all the air circulation systems.
Sorry if I’m getting to be a bore.

kim
September 6, 2009 8:01 pm

Roy Spencer 10:34:43
That’s easy unless I’m wrong; the sensitive system cycles above and below the equilibration temperature, before settling, and the insensitive system goes straight to it
====================================

Paul Vaughan
September 6, 2009 11:02 pm

Re: Bob Tisdale (14:40:04)
Thanks for the comments Bob.
If you find an ACC (Antarctic Circumpolar Current) link, that will be very interesting.

Paul Vaughan
September 6, 2009 11:11 pm

Bob, I’ve just been thinking about methods of visualization that might help with this mystery that has arisen. It will be great if someone designs a tool that will allow a user-defined hovmollering reference-line (via mouse-trace, for example) rather than the usual fixed- lat or long. I’ll continue pondering other alternatives…

September 7, 2009 3:31 am

Paul Vaughan: I’ve been trying to create SSH anomaly Hovmollers like so…
http://www.aviso.oceanobs.com/en/data/data-access-services/live-access-server-las/lively-data/2003/dec-22-2003-longitude-time-diagram-2/index.html
…at the AVISO site, but I keep getting error messages. Here’s their Live Access Server link:
http://www.aviso.oceanobs.com/en/data/data-access-services/live-access-server-las/index.html
Maybe you’ll have better luck.
Regards

coaldust
September 7, 2009 7:29 am

Bob Tisdale (04:42:15)
I was refering to the first graphic, although I wrote graph titled “NOAA/NESDIS SST Anomaly”. There is an area of relatively cool water in that graphic at about 118 W 5 N. I’m wondering if it is common during a develpoing El Nino to have an area of relatively cool water such as that. Perhaps someone can enlighten me.

coaldust
September 7, 2009 7:30 am

The first sentence in my previous post should read “I was refering to the first graphic titled ‘NOAA/NESDIS SST Anomaly’, although I wrote graph”.

September 7, 2009 10:41 am

Great stuff – even if it is complex and raises more questions than it answers! At least these questions are being approached. I have tried, as a general ecologist, to get a handle on the relationship of ENSO, the PDO, NAO and Arctic cycles from the point of view of mass movements of warm water generated in the equatorial regions by what appear to be pulses in the short-wave flux to the ocean surface. If you look at Hadley and NOAA (Scripps Inst., papers by Lyman and others) where they subtract various years of upper ocean heat content to get a sense of change according to the surface grid, it is clear that accumulated heat is primarily stored at depth in northern waters – in the main gyres of the north Pacific and north Atlantic. If we can figure out how these heat stores are built up over time and then how they are depleted over time, we may find our way to understanding the long term natural cycles that are missing from AGW models.
Cloud cover correlations are not easy to find. ISCCP data contradicts ship-board data, Earthshine data gives another picture. The ISCCP low level cloud data does show a 4% decline globally from 1983-2000 and this is certainly enough to account for the ocean warming (and land warming is a consequence of heat transfer).
With regard to Levitus’s new analysis – I am not sure how reliable it is. Prior to this 2008 paper, his previous estimates (before 2005) were used by climatologists to confirm their models (with much trumpeting) – but three further analyses (Gouretski and Koltermann, Palmer and others at Hadley, and Domingues and team) all held that these previous estimates were out by a factor of two – which is rather a lot when considering the planetary budget – this gave rise to Pielke’s conclusion that there was no more ‘warming in the pipeline’. Obviously, Levitus has responded to the critiques – at least to Gouretski’s and he cuts by half their previous cut – it is still lower and still a problem for the modellers than used the early analysis as validation of their models. But we might expect Levitus not to agree with the other three, and must wait until they respond. As you note, Bob, it is a complex field and not one for the generalist to second-guess.
I am surprised, however, at the size of the 21st century rise to 2000-2003, though not by the plateau thereafter, or the fall in the northern Atlantic. As many comments point out, it is this upper ocean heat content that is the only real measure of global warming or cooling that we have – and despite the vast array of monitors, there is still great uncertainty.
The spatial distribution of the deeper warm water pools holds the key to future patterns – I would like to see how these respond to the shifting jetstream as they are both in its track. Lower cloud in northern latitudes will lead to loss of heat and in equatorial latitudes, heat gain. More cloud in the north will insulate the waters and delay release, whereas in the tropics it will lower the heat gain.
Right now, I think the massive shift in the jetstream of the past 3 summers has depleted the Atlantic warm water pool (if SST is anything to go by) – the heat has been dumped in the UK and western Europe as rain! Hadley’s SSTs are showing a growing cold patch in the central northern Atlantic – this is the NAO on the turn, and the reason some modellers are now predicting cooling for the next decade. Hadley has such a model but they seem rather shy about, only publicising what they call the more certain outcomes of their long-term model to 2050-2080 (i.e. warming).
Once again, great work – a tribute to this blog and to the quest for truth!

September 7, 2009 11:45 am

coaldust: Sorry. I forgot that Anthony put up the map that wasn’t part of my post. Here’s another map (Sept 6, 2009) that show the little cool spot just east of the NINO3.4 region.
http://weather.unisys.com/surface/sst_anom.html
That cool spot also shows up in the subsurface anomalies:
http://www.cpc.ncep.noaa.gov/products/analysis_monitoring/enso_update/wkxzteq.shtml
It’s definitely not a stong El Nino. There was a preliminary pool of warm surface water that appears to be breaking up. However, as shown in the animation above and in the following map of SSH, there’s still some more warm water in the equatorial Pacific that hasn’t made it’s way to the surface.
http://sealevel.jpl.nasa.gov/science/jason1-quick-look/2009/images/20090821G.jpg
The archive of those SSH maps is here:
http://sealevel.jpl.nasa.gov/science/jason1-quick-look/archive.html

Paul Vaughan
September 7, 2009 1:18 pm

Re: Bob Tisdale (03:31:33)
The software is very slow, the series only seem to go back to 2001, and frames are forced to be parallel & perpendicular to lines of longitude & latitude (i.e. it won’t even allow angled-straight-lines, let alone curved reference-lines & area-averaged hovmollers (as an integrated alternative to line-transects)), but I encounter few problems otherwise plotting from here:
http://las.aviso.oceanobs.com/las/servlets/dataset
Online software like this is a step in the right direction, but there are miles to go.
Data visualization needs to move up in the climate funding hierarchy. (The complex stuff we are going after isn’t going to show up in computer fantasies if no human has ever detected or imagined it.)

Paul Vaughan
September 7, 2009 2:23 pm


Indian-Pacific Hovmoller 2001-2009

Note the temporal variation in the pattern at ~75E. Sometimes what appears to be a wave from the South Pacific summer keeps propagating within the annual band, but in other instances (2 in particular) it appears to slow down (&/or split) thereby spilling into the next year. This is the sort of thing we need to be looking for to understand interannual variations (such as ENSO).
Clearly a longer temporal window (and more control over the geometry of the reference area/line) would be helpful.

Paul Vaughan
September 7, 2009 4:06 pm


Indian-Pacific Hovmoller 11S 2001-2009

Looking 5 degrees further south (at 11S), reveals 1, 2, & 4 year waves. It looks like there are also 8 year waves, but the record is short. Note particularly how the 1 year waves cross the 2 years waves south of the tip of India (~75-80E) in some years.
Thanks for pointing out another useful site Bob.

Here’s a look at 61S.

Note the structure just east of NZ and off the tip of S.America. Considered in conjunction with the Antarctic Circumpolar Current (ACC), the Antarctic Circumpolar Wave (ACW), the Antarctic Oscillation (AAO), and Bob’s findings about the Southern Ocean, this continues becoming more & more interesting.

kim
September 7, 2009 5:17 pm

Thanks, Bob, for this extremely helpful thread, and thanks, Peter, Paul, and all the other Harrys, for the tips about ocean heat content. You have it, as did Arnd Baerndt:
Climate is the continuation of the oceans by other means.
====================================

beng
September 8, 2009 8:38 am

I think Stephen Wilde & Kim above have the right answer to Roy Spencer’s question above. The sensitive system will oscillate around the mean (“ring”) for an extended time, while the insensitive will be “damped” and in a short time quickly return to the mean. Essentually the sensitive system “stores” energy (in the pipeline, which then takes a while to come out), while the insensitive does not store much if any.
We see this action in the global temp records when a global impulse (+ or -) occurs, like Pinatubo or the 1998 El Nino (the time period for recovery is less than a yr!) We’ve even seen it last month on a much shorter time scale (because the impulse was prb’ly not global & more localized) when temps spiked up and in less than a month lost most of the “spike”.
IOW, little to no heat energy in the pipeline.

September 8, 2009 4:47 pm

Anthony: As always, thanks.

October 23, 2009 5:53 am

Independent Forensic Investigation If Visa deems necessary, an independent forensic investigation must be conducted by a QIRA. ,