Why Joe Bastardi sees red: A look at Sea Ice and GISTEMP and starting choices

AcuuWeather meteorologist Joe Bastardi has a question about two datasets and asks: If it is darn warm, how come there is so much sea ice?

click image to watch the video

Bastardi asks a simple question: how can we have above normal temperatures in the Arctic and the Antarctic and continue to add to the global sea ice trend? After all we’ve been told by media stories that both the Arctic and the Antarctic continue to melt at a frenetic pace. But it looks like this year we’ll see another Arctic recovery as we’ve seen in 2008 and 2009.

Bastardi also wonders about something we routinely ask about here at WUWT: data adjustments. GISS seems to be stuck with Arctic positive anomaly, yet the sea ice isn’t cooperating. Of course just having a positive temperature anomaly doesn’t guarantee melt, but members of the public who are less discerning, who look at red hot color presentations like GISS puts out, usually can’t tell the difference.

For reference here are the images Joe uses in his presentation. I’m going to help out a bit too with some simple comparisons.

First The GISS Dec-Feb 2010 Global Surface Anomaly as Joe presents it in his video:

click to enlarge

Source: http://data.giss.nasa.gov/cgi-bin/gistemp/do_nmap.py?year_last=2010&month_last=2&sat=4&sst=0&type=anoms&mean_gen=1203&year1=2010&year2=2010&base1=1951&base2=1980&radius=1200&pol=reg

Note that in the warmest places in the Arctic according to GISS, there are few if any land thermometers:

click to enlarge

Above: map of GHCN2 land stations (thanks to commenter Carrick at Lucia’s)

Note the cross section of the GISS data, most of the warmth is at the Arctic where there are no thermometers. The Antarctic comes in a close second, though it has a few thermometers at bases on the perimeter of the continent plus a couple at and near the center. Note the flat plateaus are each pole.

The effects of interpolation become clearer when you do a 250 km map instead of 1200 km:

click to enlarge

Source: http://data.giss.nasa.gov/cgi-bin/gistemp/do_nmap.py?year_last=2010&month_last=2&sat=4&sst=0&type=anoms&mean_gen=1203&year1=2010&year2=2010&base1=1951&base2=1980&radius=250&pol=reg

All of the sudden, the hot Arctic disappears. It disappears because there are no thermometers there as demonstrated by the cross section image which stops at about 80N.

Interestingly, the global surface anomaly also drops, from 0.80°C at 1200km of interpolation to 0.77°C with an interpolation of 250km.

One of the things that I and many other people criticize GISS for is the use of the 1951-1980 base period which they adopted as their “standard” base period. That period encompasses a lot of cool years, so anomalies plotted against that base period will tend to look warmer.

This famous GISS graph of surface temperatures from weather stations, shown worldwide in media outlets, is based on the 1951-1980 period:

Uncertainty bars (95% confidence limits) are shown for both the annual and five-year means, account only for incomplete spatial sampling of data.”]GISS doesn’t provide a utility to replot the graph above with a different base period on their webpage here http://data.giss.nasa.gov/gistemp/graphs/ but I can demonstrate what would happen to the GISS global maps using a different base period by using their plot selector here http://data.giss.nasa.gov/gistemp/maps/

Watch what happens when we use the same base period as the UAH satellite data, which is 1979-2009. The 1200km interpolated global temperature anomaly for Dec-Jan-Feb 2010 drops more than half to 0.38°C from 0.80°C. That number is not so alarming now is it? As for the graphic,  the flaming red is still there in the same places because the anomaly map colors always stay the same, no matter what the absolute temperature scale is. In the first map with the 1951-1980 base period, the max positive anomaly was 6.4°C for 1200km and  8.8°C  for 250km, while in the one below with the 1979-2009 base period the max positive anomaly of 7.1C  If colors were assigned to absolute temperatures, this map would look cooler than it’s counterpart with the 1951-1980 base period.

click to enlarge

Source: http://data.giss.nasa.gov/cgi-bin/gistemp/do_nmap.py?year_last=2010&month_last=2&sat=4&sst=0&type=anoms&mean_gen=1203&year1=2010&year2=2010&base1=1979&base2=2009&radius=1200&pol=reg

And here’s the 250km presentation, note that the global surface temp drops to 0.34°C

click to enlarge

Source: http://data.giss.nasa.gov/cgi-bin/gistemp/do_nmap.py?year_last=2010&month_last=2&sat=4&sst=0&type=anoms&mean_gen=1203&year1=2010&year2=2010&base1=1979&base2=2009&radius=250&pol=reg

So it is clear, with the GISS anomaly presentation, you can look at it many different ways, and get many different answers. Who decides then which maps and graphs with what base periods and interpolations get sent out in press releases? Jim? Gavin?, Reto? Consensus over coffee at Monks?

The answer as to what base period GISS chooses in temperature anomaly maps to present to the public is easily answered by looking at their main page here: http://data.giss.nasa.gov/gistemp/

Here’s a thumbnail of the page, and the full size version of the second graph from the top, note the caption on the top of the graph:

Clearly, they prefer the base period of 1951-1980 as the default base period for the public presentation [as well as 1200 km smoothing]  and by choosing that, the GISS results look a lot more alarming than they might be if a different base period was used, such as the 1979-2009 period used by UAH and RSS.

Anomalies can show anything you want based of choosing the base period. For example, a simple thought experiment. I could choose a base period from 11,000 years ago, during the last ice age, and plot maps and graphs of today’s temperatures against that base period. Would we see red? You betcha.

Here’s a graph that shows reconstructed northern hemisphere temps at the end of the last ice age 11k years ago, they were about 4.5°C cooler than today. Granted it’s not a global temp, but close enough for government work.

So if I used a 30 year slice of temperature 11,000 years before the present as a baseline period, our GISTEMP map would look something like this:

Obviously the map above is not an accurate representation, just a visual guesstimate. The more excitable who frequent here will likely cry foul at my abuse of the image. But it does illustrate how choices of colors and baseline periods can have a distinct effect on the final visual. Using a cold baseline period in the past (in this case 4.5°C globally cooler than the present) makes the present look broiling hot.

Anomalies are all about the starting choices made by people. Nature doesn’t give a hoot about anomalies. Generally, people don’t either. Imagine if your local TV weather forecaster gave tomorrow’s forecast in anomalies rather than absolute temperatures. He might say something like:

It’s going to be a hot one folks! Tomorrow we’ll have a high temperature that is 0.8C warmer than the 1951-1980 historical baseline for this city. Dress accordingly.

Useful isn’t it? Even more useful if he’s the weatherman in Svalbaard and people anticipating a heat wave go out in shorts and tank tops in mid February.

While anomalies are fine for illustrating many things, including temperature, bear in mind it’s all about the starting conditions chosen by the individuals doing the analysis. It’s all about choosing a baseline “normal”, which is subjective.

So when Joe Bastardi looks at the GISS map of the world, sees red, and wonders why we have a growing ice presence, the answer is in the choice of baseline and the choice of colors used to calculate and represent the anomaly.

0 0 votes
Article Rating

Discover more from Watts Up With That?

Subscribe to get the latest posts sent to your email.

185 Comments
Inline Feedbacks
View all comments
RR Kampen
March 25, 2010 1:59 am

Re: Richard M (07:18:19) :
Logic?????? As you can see from your chart the amount of 1st ice increased significantly in 2008. This became a large increase in 2nd year ice in 2009. Are you really claiming that we won’t see a huge increase in “older” ice this summer? If you are then I can easily ignore all of your future posts. Only someone completely devoid of logic would make such a statement.
You may be forgetting about the dynamic nature of the pack. So a continuing decrease in multiyear ice is quite possible even if the extent doesn’t decrease.
http://nsidc.org/sotc/images/npseaice_ssm_200902.png
You are also making a too wild induction, in thinking that one of my posts would render anything I ever say in future worth ignoring 🙂

March 25, 2010 2:34 am

>>Steve Dallas
>>Robinson which almost halves the area distortion at the poles.
>> http://en.wikipedia.org/wiki/Robinson_projection
Yes, now this would be a better projection to use. Much more true to reality.
.

PeterB in Indainapolis
March 25, 2010 6:48 am

In reply to Humanity Rules:
Why, precisely, do “trends” matter?
There were several 1000-2000 year long periods during the last 11,000 years in which it was warmer than it currently is now. This warming trend did not continue unabated and fry the globe, instead, it was followed by significant periods of cooling.
The current “trend” is no different. There are no “trends” in climate, only CYCLES.

Ryan
March 25, 2010 8:57 am

“multiple satellites measuring various climate data.”
Satellites sensors are only accurate to 0.5Celsius, and monitor average temp over a 45km footprint.

March 25, 2010 9:44 am

Richard M (07:18:19) :
RR Kampen (01:29:02) :
“Maybe the Arctic sea ice ‘recovered’ to some extent because it is winter.
Also there was this peculiar temperature anomaly distribution associated with a record low NAO index.
Likely all Arctic sea ice older than two or three years will be gone by September 2010. Because that, in reality, is the trend.
http://nsidc.org/images/arcticseaicenews/20091005_Figure5_thumb.png
Logic?????? As you can see from your chart the amount of 1st ice increased significantly in 2008. This became a large increase in 2nd year ice in 2009. Are you really claiming that we won’t see a huge increase in “older” ice this summer? If you are then I can easily ignore all of your future posts. Only someone completely devoid of logic would make such a statement.

Actually only someone who thinks the Arctic ice is stationary and only melts in situ would make the statement you made.
If you look at the 2009 image in the figure linked above you’ll see that the 2nd year ice and some MY ice is situated north of the Fram Strait.
A MODIS image from yesterday shows that the ice in that region is substantially fragmented and being pushed out of the Fram to points south where it will melt later this year. By continuity that ice is being replaced by new ice upstream. Some really beautiful cloud streets in that image by the way.
http://rapidfire.sci.gsfc.nasa.gov/realtime/single.php?2010083/crefl1_143.A2010083134500-2010083135000.500m.jpg
Here’s the latest 6-day drift for the Arctic which shows the regions of MY and 2nd year ice being pushed into a strong transpolar drift and out of the Fram:
http://i302.photobucket.com/albums/nn107/Sprintstar400/Drift20100317-20100323.jpg
MODIS also shows the North water polynya is open with no ice bridge as in 2007 so that will also be a loss of MY ice.
http://rapidfire.sci.gsfc.nasa.gov/realtime/single.php?2010083/crefl1_143.A2010083202000-2010083202500.500m.jpg
This normally doesn’t happen until May/June.
It looks like 2010 is shaping up for a high loss of MY and 2nd year ice.

March 25, 2010 10:16 am

Ralph (02:34:27) :
>>Steve Dallas
>>Robinson which almost halves the area distortion at the poles.
>> http://en.wikipedia.org/wiki/Robinson_projection
Yes, now this would be a better projection to use. Much more true to reality.

OK here it is courtesy of GISS: http://www.giss.nasa.gov/tools/gprojector/
http://i302.photobucket.com/albums/nn107/Sprintstar400/Robinsongiss.jpg

Anu
March 25, 2010 11:50 am

@Willis Eschenbach (01:08:58) :
I don’t think CRU uses global average temperature anomalies for empty grids, I’d have to see a citation for that.
http://www.cru.uea.ac.uk/cru/data/temperature/
I’m pretty sure I found it off this page, if you’re interested. I don’t have my notes with me at this computer…
Nope. They use the Reynolds and Smith data, which doesn’t include the Arctic Ocean.
A 1993 paper ? Perhaps in the late 80’s GISS didn’t use satellite data in the Arctic, but a more recent citation from late 2007:
http://data.giss.nasa.gov/gistemp/references.html
Smith, T.M., R.W. Reynolds, Thomas C. Peterson, and Jay Lawrimore, 2008: Improvements to NOAA’s Historical Merged Land-Ocean Surface Temperature Analysis (1880-2006). Journal of Climate,21, 2283-2296.
http://www.ncdc.noaa.gov/oa/climate/research/sst/papers/SEA.temps08.pdf
shows that they use the Pathfinder satellite data from the NOAA, and the paper above shows they explicitly deal with sea ice affects on SST in the Arctic.
Pathfinder info can be found here:
http://www.nodc.noaa.gov/sog/pathfinder4km/userguide.html
Measuring the correlation of two datasets always includes time.
Not true. You could look at the correlation of height and weight among all males in a certain region on a single day, for instance. If you’re saying all measurements occur at an actual time, this is trivially true.
Again, I don’t think anyone uses “average global warming” to infill the Arctic.
They use only the data available and the rest is missing.
Yes, the hadcrut3gl.txt file shows the percentage of the planet that has been measured, e.g. 83% in January 2009.
http://www.cru.uea.ac.uk/cru/data/temperature/hadcrut3gl.txt
Calculating the temperature anomaly for the entire planet from 83% surface data, is equivalent to filling in the “missing” grids with planetary averages.
If you find differently, please mention it here.
So only GISS does this nonsense.
Yes, only GISS tried to guesstimate the missing data in the era before satellite coverage, because their climate models made predictions for every grid on the planet and they wanted to compare the two. Hadley could just report what was measured, and leave large gaps in their data.
Well, we don’t know, do we, that’s my point … but the issue is not the guess for the old, sparse years. It is the assumption that you can estimate trends in a radius of 1200 km.
Fine, you’d rather have no estimated data in the gaps. Your choice.
Since (as you point out below) the satellites only go to 82.5N, I’m not sure what you mean by this.
My point is that once they get surface temperature data for that final 7.5 deg Arctic circle, you can argue that “previous data” was not full coverage, so perhaps the warming in the Arctic was less than they estimated for their 82.5 deg coverage. Or, maybe it was even warmer.
As I say, measurements improve over time.
True … but what does that have to do with extrapolating trends out 1200 km from the nearest station?
It means as more money is spent on better measurement systems, the old techniques of guesstimating are no longer needed. Actual measurements trump proxies, optimal interpolation techniques, weighted averages, whatever.
Then, people will say “but you only have x years of this new, great data. Not enough to prove a longterm trend”.
You know the drill.

Anu
March 25, 2010 2:46 pm

@Willis Eschenbach (01:18:11) :
And yet the comparison to a baseline value varies with different baselines:
0.281 million sq km below one baseline, such as in the graph you gave, might be 0.600 million sq km below a different baseline.
Wasn’t somebody just saying that different baselines make the temperature anomalies look more extreme ?
I’m glad you see the value of “up to date” research – I wouldn’t want to fixate on 1987 or 1993 papers, for example…
The chart I gave showed how the average is affected by including the recent, more extreme Arctic melting. This allows more perspective for the other NSIDC charts such as:
http://nsidc.org/data/seaice_index/images/daily_images/N_stddev_timeseries.png
which use the 1979-2000 averages.
But quite right, 30 years and more is the time period of interest in climate studies. Now that they have 2009 data, perhaps future plots will all be made with the 1979-2009 baseline.
Well, you know, that’s why I provide citations to my sources, so you can examine them and answer these kinds of questions rather than read my (perhaps mistaken) interpretation of them, and make your own judgements on them.
Perhaps you intended to give a link to the paper ” Polyakov et al.:”, and gave a link to a textfile of data instead:
http://www.frontier.iarc.uaf.edu/~igor/research/data/composite_sat_slp
Since Polyakov has published many papers, I don’t know where this particular set of temperatures from 1875 onwards comes from, but you seem to trust the data, so I thought you would mention more about where they were from. I am curious how these old measurements could be better than anything GISS has – perhaps Russia has some data they are not sharing.
I can’t say if the GISS record is “better” if I have no idea where your data came from.
Cheers.

Editor
March 26, 2010 10:12 pm

Anu (14:46:07)

@Willis Eschenbach (01:18:11) :
And yet the comparison to a baseline value varies with different baselines:
0.281 million sq km below one baseline, such as in the graph you gave, might be 0.600 million sq km below a different baseline.
Wasn’t somebody just saying that different baselines make the temperature anomalies look more extreme ?

If you are concerned with how it looks, fine, your point is well taken. I’m concerned with how it is, not how some particular baseline makes it look.

Well, you know, that’s why I provide citations to my sources, so you can examine them and answer these kinds of questions rather than read my (perhaps mistaken) interpretation of them, and make your own judgements on them.

Perhaps you intended to give a link to the paper ” Polyakov et al.:”, and gave a link to a textfile of data instead:
http://www.frontier.iarc.uaf.edu/~igor/research/data/composite_sat_slp
Since Polyakov has published many papers, I don’t know where this particular set of temperatures from 1875 onwards comes from, but you seem to trust the data, so I thought you would mention more about where they were from. I am curious how these old measurements could be better than anything GISS has – perhaps Russia has some data they are not sharing.
I can’t say if the GISS record is “better” if I have no idea where your data came from.

All you have to do is go up two steps in that link and you’ll find out … but you knew that already, I suspect. Or perhaps not, but that’s where the info is.

Editor
March 26, 2010 10:35 pm

Anu (11:50:50)

@Willis Eschenbach (01:08:58) :

I don’t think CRU uses global average temperature anomalies for empty grids, I’d have to see a citation for that.

http://www.cru.uea.ac.uk/cru/data/temperature/
I’m pretty sure I found it off this page, if you’re interested. I don’t have my notes with me at this computer…

I didn’t find it there.

Nope. They use the Reynolds and Smith data, which doesn’t include the Arctic Ocean.

A 1993 paper ? Perhaps in the late 80’s GISS didn’t use satellite data in the Arctic, but a more recent citation from late 2007:
http://data.giss.nasa.gov/gistemp/references.html
Smith, T.M., R.W. Reynolds, Thomas C. Peterson, and Jay Lawrimore, 2008: Improvements to NOAA’s Historical Merged Land-Ocean Surface Temperature Analysis (1880-2006). Journal of Climate,21, 2283-2296.
http://www.ncdc.noaa.gov/oa/climate/research/sst/papers/SEA.temps08.pdf
shows that they use the Pathfinder satellite data from the NOAA, and the paper above shows they explicitly deal with sea ice affects on SST in the Arctic.
Pathfinder info can be found here:
http://www.nodc.noaa.gov/sog/pathfinder4km/userguide.html

Still nope. They are concerned with sea surface temperature. They use the satellites to figure out ice concentration, and adjust the temperature towards the freezing point as there is more ice. When the ice is above 90%, they use -1.8°C (freezing point of sea water) for the sea surface temperature.
But that tells us nothing about what we want to know, which is not the water temperature below the ice, but the air temperature above the ice.

Measuring the correlation of two datasets always includes time.

Not true. You could look at the correlation of height and weight among all males in a certain region on a single day, for instance. If you’re saying all measurements occur at an actual time, this is trivially true.

Oh, please. This is meaningless nitpicking. We are discussing temperature datasets which include time.

Again, I don’t think anyone uses “average global warming” to infill the Arctic.

They use only the data available and the rest is missing.
Yes, the hadcrut3gl.txt file shows the percentage of the planet that has been measured, e.g. 83% in January 2009.
http://www.cru.uea.ac.uk/cru/data/temperature/hadcrut3gl.txt
Calculating the temperature anomaly for the entire planet from 83% surface data, is equivalent to filling in the “missing” grids with planetary averages.
If you find differently, please mention it here.

Nope. If you have missing grids, you don’t have a planetary average, just an average for some parts of the planet. So you couldn’t be infilling with the planetary average. If I average the height of ten men, am I infilling the rest of the missing men on the planet with the “planetary average height”?

So only GISS does this nonsense.

Yes, only GISS tried to guesstimate the missing data in the era before satellite coverage, because their climate models made predictions for every grid on the planet and they wanted to compare the two. Hadley could just report what was measured, and leave large gaps in their data.

Nope. Look at the top of the page. See the gray areas? Large gaps in their data … GISS infills up to 1200 km out, and they also infill out to 250 km out. Either way leaves gaps.

Well, we don’t know, do we, that’s my point … but the issue is not the guess for the old, sparse years. It is the assumption that you can estimate trends in a radius of 1200 km.

Fine, you’d rather have no estimated data in the gaps. Your choice.

If GISS used a reasonable estimation method, I’d have much less of a problem with it. As it is, they assume that since temperatures are correlated out to 1200 km, they can use the trends out to 1200 km. Bad math, no cookies. And you still haven’t commented on that …

Since (as you point out below) the satellites only go to 82.5N, I’m not sure what you mean by this.

My point is that once they get surface temperature data for that final 7.5 deg Arctic circle, you can argue that “previous data” was not full coverage, so perhaps the warming in the Arctic was less than they estimated for their 82.5 deg coverage. Or, maybe it was even warmer.
As I say, measurements improve over time.

The issue is not better data, or improved measurements. It is that for some parts of the earth we have no data at all. I object to filling in those blank parts of the planet with red “here there be dragons”.

True … but what does that have to do with extrapolating trends out 1200 km from the nearest station?

It means as more money is spent on better measurement systems, the old techniques of guesstimating are no longer needed. Actual measurements trump proxies, optimal interpolation techniques, weighted averages, whatever.
Then, people will say “but you only have x years of this new, great data. Not enough to prove a longterm trend”.
You know the drill.

I do know the drill. For me, it’s to use the data we have, and not pretend that we know what’s happening were there is no data. This is as true in time as in space. If we only have x years of data, that’s all we have. If it’s not enough to show a longterm trend, well, that’s the way it is. Should we extrapolate that data backwards in time, just like GISS extrapolates data outwards in space? By your reckoning I guess we should, can’t just leave those early years blank any more than we can leave the Arctic blank, by your lights we should put in our best guess.
The field of climate science is polluted by the inability of many climate scientists to say three simple words … “We don’t know”. We don’t know what the temperatures are like where we don’t have thermometers. We don’t know what happened before the start of the records. We don’t know what’s happening in the blank gridcells. We don’t know a whole host of things about the climate … and yet everyone keeps saying “The science is settled”.
Riiiight.

1 6 7 8