Guest Post by Willis Eschenbach
When I was a kid I had the great fortune to be taught to splice rope by a grandson of Richard Henry Dana. He taught me how to do a long splice and a short splice … but I never learned how to do either the long or short data splice. For that, I had to wait for climate science.
Anthony recently highlighted a paper about reconstructed Greenland temperatures called “Past 4000 Years of Greenland Temperatures“. I don’t really have a dog in the fight, but I was greatly amused by their description of what they call “in situ observations”.
Here’s the top panel of their Figure 1.
Figure 1. ORIGINAL CAPTION FROM THE PAPER (top) Reconstructed Greenland snow surface temperatures for the past 4000 years and air temperature over the past 170 years (1840–2010) from three records. The thick blue line and blue band represents the reconstructed Greenland temperature and 1s error, respectively (this study). The reconstruction was made by two different methods before and after 1950. The “gas method” is as described in section 2, and the “forward model” is described by Kobashi et al. [2010]. Thick and thin black lines are the inversion‐adjusted reconstructed Summit annual air temperatures and 10‐year moving average temperatures, respectively [Box et al., 2009]. Thin and thick red lines are the inversion adjusted annual and 10‐year moving average AWS temperature records, respectively [Stearns and Weidner, 1991; Shuman et al., 2001; Steffen and Box, 2001; Vaarby‐Laursen, 2010]
I thought “hmmm, data from the summit of the Greenland ice cap, didn’t know about that” … so I went to find out what their data consists of, the famous “thin red line” in the Figure 1 above. Here’s what they say about the “in situ observations”:
5.1. Present Greenland Temperature
To place the Greenland temperature proxy reconstruction into a historical context, we incorporate two additional Summit temperature records. One record is obtained from a compilation of Summit Automatic Weather Station ∼2 m surface air temperature (SAT) observations (hereafter AWS or in‐situ record) that spans 23 years (1987–2010). The AWS were situated within 20 km of the GISP 2 coring site and within 25 m elevation of the ice sheet topographical summit (Figure 1, top, red line).
The series begins in May 1987 with Automatic Weather Station data after Stearns and Weidner [1991].
Shuman et al. [2001] merge this record with data from the Greenland Climate Network (GC‐Net) AWS data [Steffen and Box, 2001] to produce the first 12 years of this compilation.
Gaps before June 1996 are in‐filled using daily passive microwave emission brightness temperatures.
GC‐Net data then comprise the period spanning June 1996 to December 2003 with gaps in‐filled by Danish Meteorological Institute (DMI) Summit AWS data [Vaarby‐Laursen, 2010].
The DMI data exclusively form this data series from January 2004 through December 2010.
Man, I thought, that is a curious provenance for the summit temperature data. First, 12 years of observations, with gaps in-filled using satellite microwave observations. Then from 1999 to 2003 the summit temperature was estimated using surrounding stations. Gaps in the surrounding station data are infilled from the DMI Summit station data. Then from 2004 on, we have DMI summit station data only.
So I looked a bit deeper. It gets better. Remember the Shuman et al. observations that form the “first 12 years of this compilation”? It turns out that it was not one station, but three stations … and there was no overlap between the stations to compare temperatures. The Shuman paper is here. Their Figure 2 shows those first 12 years of the record, from three AWS stations: CATHY, GISP2, and SUMMIT.
Figure 2. Automated Weather Station records near and at the Greenland Summit.
Note the lack of overlap and the number of gaps, covering days, weeks, or months. Here’s Shuman’s description of what they did:
In order to complete a temperature record from the Greenland Summit (May 1987 to October 1999), it will be necessary to adjust the AWS Cathy temperature record to account for the difference in its location as well as to complete the AWS GISP2 and AWS Summit records across multiple periods of missing data .
The methods detailed in Shuman et al. (1996) or in Shuman et al. (1995), which rely on appropriately located and contemporaneous SSM/I brightness temperature data (Table 1), will be used to achieve a complete and consistent temperature record.
Inconveniently, SSM/I data are not available from 4 May 1987 to 10 July 1987 and are missing from 3 December 1987 to 13 January 1988 (Table 3). Smaller gaps in the SSM/I or AWS record of less than 5 days will be dealt with by interpolation.
So for the 23 year record we have a long data splice as follows:
2 years of AWS CATHY data,
7 years of AWS GSIP2 data,
3 years of AWS SUMMIT data,
5 years of estimated data, and
6 years of AWS SUMMIT data.
All of this is “infilled” from a couple of sources.
Now, remember that this data is what they are using to calibrate their algorithm that converts ∂O18 data into temperature data …
The authors of the study cite Ellen Vaarby-Larson 2010, which is here. It shows the Summit temperature for the period 1998- Feb 2010:
Figure 3. Air temperature from the 04416 Summit AWS station. From Vaarby-Larson 2010.
Note that even in the modern period there are gaps of up to about a year in the record. (Also of note is that there are only about six three-hour periods in the record where it is above freezing, and that the summer/winter spread is about 60°C [108°F]. Yikes!) Here’s an example of why the record contains so many gaps and spaces:
Figure 4. Condition of the Summit station 04416 during the 2007 annual visit. From Vaarby-Larson 2010
Let me conclude by looking at some other problems endemic with AWS records. I cannot improve on the words of Ms. Vaarby-Larson, who said (op. cit.):
Recommendations
Generally, great care should be taken when using the observations from station 04416 Summit, since the observations are influenced by:
• Extreme climatic exposure of the station measurement equipment (the extreme cold causing e.g. low availability of wind observations during winter)
• Non static barometer elevation above sea level, due to the Greenland Ice sheet flow-patterns
• Non static height of measurements above ground, due to burial of the station by snow falling
None the less, the observations of temperature, humidity and wind (see Figure 10 – Figure 16 ) show no great, obvious shift in level or variance, at the individual station relocations in 2005, 2006, 2007 and 2009.
Any existing bias due to change in measuring height above ground, might be difficult to identify since:
• 04416 Summit only issues synop every three hours.
• There’s no direct measurement of the ongoing change in actual measuring height above ground, in-between station relocations.
• Change in weather conditions might produce the same signal in the observations
• Natural variability shows e.g. great temperature variance during winter compared to summer, please confer with the early attempt at investigation of temperature variance in Figure 17.
Hmmmm …
In a recent post called “A Modest Proposal—Forget About Tomorrow“, I discussed a paper that showed why even perfect “fitted” or “calibrated” models may not have predictive (or reconstructive) capabilities. In the current case, they are reconstructing the Greenland temperature based on a collage of actual observations, nearby observations, estimates, and satellite microwave brightness.
Let me be clear that there is nothing inherently wrong with putting together a pastiche composed of a variety of estimates and local observations and satellite data. There may even be something to learn from such a collation of disparate elements.
My point is that it doesn’t engender confidence when said pastiche is used to calibrate an algorithm designed to transform an ice core ∂O18 record into a temperature reconstruction. Even with the best of data, that’s a tough sell.
My conclusion? Only have one.
The confidence intervals in the original paper on historical Greenland temperatures are way too narrow.
w.
PS—Ya gotta love the Google Earth image of the Greenlad Summit temperature measuring site (Station ID# 04416), found here:
When I first looked that that I thought “They left out the satellite part of the image” … then I realized that they hadn’t left anything out. Anthony, I do believe we’ve actually finally found a truly rural temperature station …
Why does the shape of the red line in figure 1 bear absolutely no resemblance to Figures 2 and 3?
Just asking… I really would like to understand.
Looks like a white out, I can’t see the station 😉
When I initially commented that the spliced in situ data at the end might be dodgy, I had no idea just how dodgy it would turn out to be. I wonder if any of this was ever caught and discussed in peer review?
James Reid says:
November 11, 2011 at 10:40 pm
Not sure, but I think they’ve used the thick red line from Figure 1 (10 year averages) rather than the annual (or monthly) data.
w.
One of my teachers once said that it is better to have an approximate answer to a good question than an exact answer to a poor question. He so saying this did so before climate “science” began fabricating and torturing data. If one doesn’t splice a rope correctly the consequences can be deadly. In climate “science” there are no consequences.
Thanks Willis.
Lemme see Willis, we’re using a coupla stations that have been relocated, to estimate the temperature at a station that has been wonky for at least some little time, from which we estimate a terminal velocity from temperature acknowledged to be even wonkier from time to time due to known spurious responses due to whatever might be discerned from “Natural variability shows e.g. great temperature variance during winter compared to summer”.
Got it..
Seems OK to me, but I did bump my head last night……..
re Google Map
“… it’s full of polar bears!”
Apologies to Clark & Kubrick (and Dave Bowman)
This is the satellite record for Central Greenland:
Monthly anomalies
http://climexp.knmi.nl/data/itlt_315-325E_70-75N_na.png
Annual anomalies
http://climexp.knmi.nl/data/itlt_315-325E_70-75N_n_mean1a.png
Pics will be active for few days.
So Figure 1 is what peer-reviewed “climate science” looks like?
Say no more.
“I don’t really have a dog in the fight”
Then I won’t really have to compare you to Michael Vick.
What a relief. Michael Vick is the lowest form of life on this planet.
summer/winter spring is about 60°C [108°F]. Yikes!)
Doesn’t that melt the snow ??
I can understand a typo of the missing minus ( this is apparently stand practice in temperature reporting) but when you did the fahrenheit conversion you should have clicked. 😉
Something is rotten in the state of Denmark and I suspect it is the Greenland Temperature reconstruction.
Take a look at the long sweeping curves in their Figure 2 Borehole temperature reconstruction down to 2,000 metres…
Now reconcile that to their jagged surface temperature graphics….
These sweeping curves look like an artefact of depth… not a temperature reconstruction.
Take a look at the temperature range in their Figure 2 Borehole temperature reconstruction down to 1,200 meters…
It ranges between -30.9 C and -31.9 C… one whole degree variability in 1,200 metres of ice…
I never realised weather and climate was sooooooooooooooo stable.
Try to find the temperature range in their Figure 2 Borehole temperature reconstruction below 1,950 meters…
It suddenly goes out of range on the graph… what is happening here…
Are they hiding an exponential artefact increase that kicks in around 1,600 meters?
The three hourly temperature records shows an annual temperature range of over 60 degrees C…
Its a very extreme environment…
Its an environment where the concept of an Average Temperature is pretty meaningless…
And I think that is the bottom line with these Greenland reconstructions: MEANINGLESS.
It looks exactly the same on Google Earth … There are some nice pics there of the Summit Station. If I knew how, I’d paste one here for you to look at.
Juraj V. says:
November 11, 2011 at 11:47 pm
Juraj, per the READ ME file from the UAH MSU data (emphasis mine):
Note that the Greenland Summit is at 3,200 metres …
Just sayin’ …
w.
malagaview says:
November 12, 2011 at 12:41 am
Actually, the solid body of the ice acts as an integrator, so the smoothness of the borehole curve is expected.
Whether it contains valuable information on surface temperature is another question. It depends on whether the thermal properties of the actual ice match up with the thermal properties of the theoretical ice. I discuss some of the issues with boreholes here.
w.
Willis:
It would be very interesting to read your take on the Greenland dating techniques.
Dating layers in Firn is one thing… at least there are layers that can be argued about…
But once you get down into ice it is another thing altogether… no layers…
This is where the computer models kick in…
Modelling “gaps” in the record
Modelling “periods of melting”
Modelling “rates of accumulation”
Modelling “flow of the ice cap”
Modelling the “annual melt rate” at the bottom of the ice cap
Modelling an ever thinning “notional annual layer” of ice in the flowing ice cap
Basically the models are circular logic… confirmation bias… call it what you like…
The input parameters to the models are a “temperature reconstruction”!
I suspect it is another case of models all the way down…
Especially as Greenland is called Greenland because it was once GREEN.
Thanks for the feedback….
Whether it contains valuable information on surface temperature is another question.
That’s the question I am trying to answer….
That’s the question I would ask before I wasted time and money on Greenland cores….
But then this is climate science 🙂
Willis
On the subject of splicing, you may be interested to look in detail at the “dance of the thermometers” in the Central England Temperature record that is presented as the oldest temperature record in the World but is no such thing, even if it does contain the oldest measurements. It is of interest because it is managed by the Hadley (Met Office) partner in the HadCRUT series and whereas the other partner Jones claims an insignificant UHI effect the CET series since 1959 has had adjustments for UHI from 0.1 to 0.3 oC depending on time of year.
From 1772 until 1852 it was a composite stitched together from seven London sites and then until 1877 from Oxford, from when until 1930 it was made up a an equally weighted composite from three sites in a triangle from the NW to SW and SE that have different climates. In 1931 another switch takes place with Cambridge tossed out and replaced by Rothamsted – one of the oldest agricultural research stations in the World with a continuous high quality temperature record going back 120 years. This mix lasted until 1958 (when the atmospheric carbon dioxide took off) when all but Rothamsted were tossed and three more introduced and then finally (so far) in 2007 one of these replaced. Details of the compilations and corrections are here:http://www.metoffice.gov.uk/hadobs/hadcet/Parker_etalIJOC1992_dailyCET.pdf
It should also be noted that the CET graphic used by the UK Government (http://ukclimateprojections.defra.gov.uk/content/view/751/500/) shows continued warming to 2010 whereas the original complied by Hadley (Download HadCET here / http://www.metoffice.gov.uk/hadobs/hadcet/ ) shows cooling substantial cooling after 2007/8!
I spent my youth near two of the current sites and can attest that they have been subject to a major increasing UHI effect – but so too have the other sites that were rural prewar and have continually urbanized. Take Rothamsted that had a population of 11000 in 1930 that grew to 15,000 in 1950 and about 27,000 now but where as far as I can tell there has been no correction for UHI. What is somewhat disturbing is that this raw data from Rothamsted (http://www.worldclimatereport.com/archive/previous_issues/vol4/v4n20/cutting.htm ) is converted by the statisticians at that site into this “hockey stick” with the “blade” from the time it joined the CET series:( http://www.rothamsted.bbsrc.ac.uk/aen/ecn/AirTemp.htm ) when “by eye” the series has three similar warming trends. A comparison of the four stations making up the recent CET record shows that they had a similar trends until around 1992 when they diverged with Rothamsted showing a major warming trend whilst Cambridge, the station it replaced went in the opposite direction and that Malvern that shows a cooling trend was replaced in 2007: http://www.anenglishmanscastle.com/archives/004438.html
Since the CET record is said to reflect the NAO, and as such can be taken as a proxy for the global record, the analysis of these series deserves more attention, especially the discrepancy between the partners in HadCRUT over corrections for UHI
malagaview says:
November 12, 2011 at 1:19 am
Dating is tricky. However, a variety of different methods are used and compared for agreement. There’s a paper here (PDF) whose abstract says (emphasis mine):
Since they give error measurements, that allows you to see how the dating error might affect any particular results that interest you. In any case, download the cited report, it will answer each and every one of your questions, and you can make up your own mind.
w.
Actually, the solid body of the ice acts as an integrator, so the smoothness of the borehole curve is expected.
The smooth curves indicate that ice acts as an integrator…
The smooth curves could also indicate that ice acts as an homogeniser…
In which case the smooth curves could just be an artefact of depth and pressure…
So I am wondering where are the studies that demonstrate that the argon and nitrogen ratios remain stable in “ice bubbles” when subjected to increasing pressure over,say, 4,000 years?
When you standardize a measurement the standard must be correct and accurate otherwise the end result will be wrong. To standardize an algorithm against a data set that is questionable, and some data from below surface level due to the surprise snow, is laughable.
Thanks for the link….
It’s the case of the exponentially disappearing ice…
As you get to the bottom of the ice core 98% of the ice is gone…
So I guess the bottom of the ice core is just “dust” and “bubbles”
Providing snow has been accumulating continuously for 161,313 years..
And no “ice” has melted for 161,313 years..
And the “bubbles” haven’t risen in 161,311 years…
And the “dust” hasn’t settled in 161,311…
And the ice has only “flowed” laterally for 161,311 years…
And I am happy that up to 98% of my sample has disappeared…
Then everything is looking good.
Malagaview
One of the ways they date the Greenland ice cores is by counting the individual annual layers. Kind of like counting tree rings. This has worked reasonably well down to ages of 50,000 years (or so) . Generally speaking layer counting works fairly well. The dated variations in ice and air bubble chemistry tie in quite nicely with other well dated climate sensitive proxies from ocean floor sediments, lake sediments and speleothems. Skepticism has its limits. Even a depths of 2 to 3 km and corresponding ages of up to 110,000 years the Greenland ice cores are each accurate to about plus or minus 2000 years, which might not be perfect, but its not bad either. The layer counted ages are archived and can be downloaded for free.
Why go to all this expense and set up only one thermometer? This is the difference between classroom scientists and engineers. There should be a dozen thermometers in an array that are correlated. You could even see when one of them went down or got buried. Jeesh, we have wasted 12 precious years and have to start over again. It takes 60 years before you might have a reliable record to correlate with O18 data. Also, the AWSs should be on platforms that could be raised to stay above the snow level – adjustments made for the change in elevation (hmm engineering again)
They are splicing snow surface temperatures to air temperatures. Anyone who has dabbled in the art and science of waxing cross country skis will tell you they can be surprisingly disconnected.
In tree rings, each line represents annual growth.
In ice snow, each line represents a snow event. How do we translate multiple lines into a single annual event? Is it a simple model calculation? ie total line count divided by the average typical snow events, in a Greenland year.
I feel, a little stupid, asking such a dumb question? I assume someone has radio dated organic matter in the ice layers, to confirm event layer conversion to annuals. Is there an issue here? GK