Replication #4: Tree Ring Site Selection Criteria

Mann et al. [2000] listed 5 quality control criteria, which were said to have been applied in MBH98 for the inclusion of tree ring sites. These criteria were again referred to in the MBH Corrigendum in July 2004. In total, the information in the SI to MBH98 showed the use of 394 sites, either as individual proxies or as a component of a principal component network. I have been able to identify original data at WDCP for 331 out of 333 non-Vaganov sites used in MBH98. 196 out of these 331 series failed one or more of the stated quality control tests. One Jacoby series failed the mean correlation test so spectacularly (mean correlation of only 0.045 between segments and archived chronology) that we contacted one of the originating authors (D’Arrigo) who requested WDCP to withdraw the data.

Summary: We examined 331 series used directly either as individual proxies or as a component of a principal component network. 196 of these series failed one or more of the stated quality control tests, most frequently the mean correlation, but occasionally each of the other tests (except for the end date of the series). Results of these calculations are shown in Appendices.

In our review of MBH98 sites, we noticed a material discrepancy between the series listed in the original SI and the series actually used. This is summarized in Table 1 below. In MM03, we had pointed out the discrepancy for the “miscellaneous” 5 series; after we obtained access to a previously unavailable directory at Mann’s FTP site in November 2003, we determined that 30 other series listed in the original SI were not used in calculations at Mann’s FTP site. We notified Nature of this and, in July 2004, these discrepancies, together with some, but not all of the problems affecting MBH98, were acknowledged in the Corrigendum, which provided a listing of the 35 series.

 The Corrigendum purported to provide an explanation for the discrepancy as follows:

These series, all of which come from the International Tree Ring Data Bank (ITRDB), met all the tests used for screening of the ITRDB data used in Mann et al. [1998] (see Mann et al. (2000)), except one—namely, that in 1997, either it could not be ascertained by the authors how these series had been standardized by the original contributors, or it was known that the series had been aggressively standardized, removing multidecadal to century-scale fluctuations.

The tests listed in Mann et al. (2000), and presumably the ones referred to in the Corrigendum, were as follows:

  • Reliable information on the methods used to remove biological trend was available;
  • The median length of the individual segments used to build the chronology was greater than 150 years;
  • The mean correlation of these individual segments with the site chronology was greater than 0.5;
  • The first year of the chronology was before AD 1626, and it contained at least 8 segments by 1680;
  • The last year was after 1970, and there were still 8 segments after 1960.

In this note, I report on verification attempts on the actual application of these criteria; I will examine the supposed explanation for the discrepancy in a forthcoming note.

Mann et al. use tree ring chronologies as part of 6 principal component networks, as individual proxies and indirectly in temperature and precipitation reconstructions from tree rings. Here I only consider quality control in the tree ring sites used in principal component networks and directly as proxies and, in another note, I will look at quality control for most of the tree ring sites used indirectly through precipitation and temperature reconstructions.

The original SI provided locations or ITRDB codes for 333 tree ring sites and reported the use of 61 Vaganov sites without providing locations. As shown in Table 2 below, out of the 333 sites with locations or codes, I was able to locate 331 series at WDCP/ITRDB, with one NOAMER series (ar045) and one Stahle.SWM series not being archived at ITRDB. I have not tried to determine the reasons for these discrepancies. Mann’s FTP site showed that 22 sites were used in the Stahle/SWM network rather than the 20 listed in the original SI. The Corrigendum attributed the 2 additional sites to Stahle, pers. comm., but unfortunately did not provide a location or code. Since the first 120-125 years of each series is identical to other series in the 20 identified sites, it appears that there has been a duplication or splicing. This remains unexplained. This affects retention policy in the 15th century. MBH98 stated that PC networks required a minimum of 7 series for calculation, except for the Stahle/SWM series,which, for some reason, only required 6 series. The unexplained duplications affect the Stahle/SWM network in the 15th century as, without the duplication, there are only 4 series. For final calculations, this probably does not “matter”, since nothing much “matters” in MBH98 except bristlecone pines, but it would be nice to reconcile the duplications.

At Mann’s FTP site, there are latitudes and longitudes and three-letter codes for the Vaganov sites. I spot-checked some Vaganov sites against ITRDB sites in similar locations and didn’t have much success in matching and accordingly excluded the Vaganov network from this quality control verification.

The quality control checks, said to have applied to the sites listed in the original SI, were checked for each of the 331 sites, which could be identified at WDCP/ITRDB. I found that 196 of the 331 sites failed one or more of the quality control criteria, said to have been applied, as summarized in Table 3 below. In some cases, there were differences between the WDCP Version and MBH version

Some Follow-Up on Quality Checks

One QC failure was so spectacular (a mean correlation between chronology and individual segments of only 0.045) that we followed up with the original author, Rosanne D’Arrigo, as follows:

The *.crn series at WDCP for the site cana158 has extremely and atypically low correlations to the individual trees used in the calculation. The *.crn series also extends to 1982, while the *rwl data goes only to 1978. By any chance, is there an error in the labelling in the series submitted to WDCP?

Also I note that several treeline *.crn series have been updated to 1990-92 (ak031, ak032 and cana177), but the *.rwl series have not been updated commensurately. Perhaps this is an oversight and, if so, perhaps you could remedy the matter.

This prompted the originating author to request the withdrawal of the chronology from WDCP.

Other QC Issues

During the course of this QC exercise, the following issues with certain WDCP/ITRDB datasets were noticed.

  •  46 tree series in 7 South American sites had missing periods in the individual tree series, with the largest gap being 89 years for tree PP1993 at site arge015. All of these series were collected several decades ago and WDCP had no information on the reasons for the gaps (Bruce Bauer, pers. comm.). MBH provided no information on how they handled these gaps in their own quality control. All 7 sites affected by gaps in individual tree series were used in MBH98.
  •  The Polar Urals site and Tornetrask site also have many gaps. In these cases, the ring widths were apparently not measured prior to densitometric calculations, and gaps were lost re-orienting the core for densitrometric measurements (P.D. Jones, pers. comm.) but the number of rings were counted.
  •  arge015 has a misprint in the rwl file in which tree PP1982 is dated one millennium too early. The above was corrected in these calculations

Appendices

Quality Control Testing QC.XLSX

  1. Stahle/OK
  2. Stahle/SWM
  3. NOAMER
  4. SOAMER
  5. AUSTRAL
  6. Jacoby
  7. Miscellaneous

 

Bring the Proxies Up to Date!!

I will make here a very simple suggestion: if IPCC or others want to use “multiproxy” reconstructions of world temperature for policy purposes, stop using data ending in 1980 and bring the proxies up-to-date. I would appreciate comments on this note as I think that I will pursue the matter with policymakers. Continue reading

Graybill and Idso [1993] Online

A lot of people who are interested in the issue of the effect of bristlecone pines on MBH98 do not have access to Graybill and Idso [1993], the underlying reference, which we referred to in both MM05 articles and which I mentioned in this post . Here’s a pdf version , courtesy of the U.S. Water Conservation Lab. I was really amazed when I connected the MBH98 PC1 to the Graybill sites and, when you read Graybill and Idso [1993], I think that you will be as well. As we’ve mentioned before, there is extraordinary irony that Sherwood Idso’s data should now be relied upon by Mann et al. as essential to their temperature reconstruction. Every single expedient presented by them, whether couched in Preisendorfer’s Rule N or other lofty talk about “significance”, has one single purpose: to get Idso’s bristlecone pine growth series into their calculation. pdf version

Replication #3: What if a step is not replicable?

This is a short discussion of the issues arising when a calculation step in a multi-step study cannot be replicated. Continue reading

Sci-AM: Confirming the importance of MBH 99

A new profile of Michael Mann appears in the March edition of Scientific American (on the grounds that you can never have too many profiles) Update: the author of the profile is David Appell. Isn’t it a small world? You’ll have to pay real money to see all of it, but here’s a rather revealing extract (my emphasis added):

To construct the hockey-stick plot, Mann, Raymond S. Bradley of the University of Massachusetts Amherst and Malcolm K. Hughes of the University of Arizona analyzed paleoclimatic data sets such as those from tree rings, ice cores and coral, joining historical data with thermometer readings from the recent past. In 1998 they obtained a "reconstruction" of Northern Hemisphere temperatures going back 600 years; by the next year they had extended their analysis to the past 1,000 years. In 2003 Mann and Philip D. Jones of the University of East Anglia in England used a different method to extend results back 2,000 years.

In each case, the outcome was clear: global mean temperature began to rise dramatically in the early 20th century. That rise coincided with the unprecedented release of carbon dioxide and other heat-trapping gases into the earth’s atmosphere, leading to the conclusion that industrial activity was boosting the world’s mean temperature. Other researchers subsequently confirmed the plot.

Continue reading

Replication #2: Selection of Gridcells

The Corrigendum SI stated that:

MBH98 made use of all nearly continuous monthly gridpoint surface temperature records (no single gap greater than 24 months, and no more than 10 years of total missing data.

I checked these criteria against the temperature dataset archived at the Corrigendum SI and found that 242 out of 1082 gridcells selected for the 1902-1993 period failed one or both criteria, while 229 gridcells in the 1854-1993 ( rather than the 219 shown in MBH98) met both criteria. Some curious differences between newer and older HadCRU versions are noticed in passing. More

Replication #1: MBH98 Temperature Dataset

I have a considerable inventory of material on replication issues pertaining to MBH98, which does not really fit into academic journal formats. I’ll probably do about 25 of these notes, which may interest a few people and will illustrate the obstacles to replicating MBH98 without a close examination of source code. I’ll start first with the temperature dataset used in MBH98, identifying the correct reference rather than the one provided in the Corrigendum SI. Continue reading

Categories

I’ve updated to the Categories on the right bar to distinguish comments on MBH98, Moberg [2005], News and Commentary, Post-1980 Proxies etc. , which cover the main current themes on this blog.

More Commentary #4: Wall Street Journal Editorial

Further to the lengthy article by Antionio Regalado on Feb. 14, 2005, the Wall Street Journal on Feb. 18, 2005 had an editorial referring to the present matter:

Mr. Mann’s chart was both a scientific and political sensation. It contradicted a body of scientific work suggesting a warm period early in the second millennium, followed by a "Little Ice Age" starting in the 14th century. It also provided some visually arresting scientific support for the contention that fossil-fuel emissions were the cause of higher temperatures. Little wonder, then, that Mr. Mann’s hockey stick appears five times in the Intergovernmental Panel on Climate Change’s landmark 2001 report on global warming, which paved the way to this week’s global ratification–sans the U.S., Australia and China–of the Kyoto Protocol….

… as the Journal’s Antonio Regalado reported Monday, Mr. Mann tried to shut down debate by refusing to disclose the mathematical algorithm by which he arrived at his conclusions. All the same, Mr. Mann was forced to publish a retraction of some of his initial data, and doubts about his statistical methods have since grown. Statistician Francis Zwiers of Environment Canada (a government agency) notes that Mr. Mann’s method "preferentially produces hockey sticks when there are none in the data." Other reputable scientists such as Berkeley’s Richard Muller and Hans von Storch of Germany’s GKSS Center essentially agree…

It says something that it took two non-climate scientists to bring Mr. Mann’s errors to light.

But the important point is this: The world is being lobbied to place a huge economic bet–as much as $150 billion a year–on the notion that man-made global warming is real. Businesses are gearing up, at considerable cost, to deal with a new regulatory environment; complex carbon-trading schemes are in the making. Shouldn’t everyone look very carefully, and honestly, at the science before we jump off this particular cliff?

Link

More Commentary #3: Melanie Phillips

Melanie Phillips says:

so many people have made huge reputations from all this rubbish. What an astonishing scientific scandal this is, and growing by the day.

Link