Archive for the ‘Paleoclimatology’ Category
The Chronicle of Higher Education has a very nice review of Bill Ruddiman’s hypothesis that land use change by Holocene human populations had a detectable influence on climate and atmospheric chemistry prior to the Industrial Revolution. What I particularly like about this article is that it appears to successfully collect in a readable format the diverse array of evidence, opinions, and personalities involved in proposing, testing, and modifying this rather interesting hypothesis. Kudos to Josh Fischman for an enjoyable article.
Via Jeff at his blog The Air Vent, in the midst of a post that includes a nontrivial quantity of ALL CAPS and multiple exclamation points (no <blink> tag?), the following statement (since that time stated in a more reasonable tone here):
For some reason EVERY RCS CORRECTION Briffa can conceive of refuses to turn upward to fit the ACTUAL data. This lack of flexibility in the RCS curves is what creates the HS.
I’m not sure why he doesn’t take my post on growth curves at Yamal seriously, but I’ve already explained why this is probably incorrect. The problem is an intermingling of the signal and the noise. As I showed previously, when I constructed an RCS curve using only the recent, living trees from Yamal, because these trees were growing roughly at the same time, and at that during a period of increasing temperatures, estimating a growth curve from them potentially intermingles the climate signal with the geometric growth curve.
One way to look at this is to consider only those trees that ceased growing prior to 1900. I’ve extracted these from the full Yamal set and estimated the regional curve standardization using a 67% spline. Here it is:
Let me be as clear as I can be, there is no sign that I can detect that it is old trees that increase their growth at Yamal (even if identified, this phenomenon would require some hypothesis as to the cause), At Yamal, a portion of the old trees are those that were growing together during a period of climate warming. If you examine the raw ring width, there are a few fossil series that have rapid increases toward the end. If Jeff’s hypothesis were correct, we’d expect these to be the oldest, right? In fact, the seven subfossil samples I identified as having rapidly increasing growth in their later years, six had a wide range of ages from 90 to 180 years (this comes with the caveat that we don’t know the exact pith age).
Furthermore, and as my earlier posts have shown, and as I’ll show again below, using a 67% spline (so that the later part of the regional growth curve bends slightly upward) you still get an increase in the chronology values during the 20th century.
Finally, what of the claim that combining the mean-detrended series demonstrates that the RCS method is invalid? One way we can test this is by first aligning all the ring widths by age (as was done here), observing where the curve of the juvenile grow trend flattens in the majority of trees (eyeballing it at around 175 to 200 years), throwing away the ring widths for the time when every tree was between 1 and 199 years of age (you get essentially the same result if you use 175 years), then realigning by time (year A.D.), and removing the mean. Lets look at the truncated growth curves in time and the mean of these series:
What immediately jumps out at you is that the mean for the most recent century and half is noticeably higher than the earlier part of the millennium. This doesn’t a priori indicate warmer temperature, but as I explained to Jeff here, once again there is the potential for this approach to remove climate signal in the process of detrending. For the most part, the mean following truncation is lower for the non-living trees at Yamal, but for the living 17 trees, the mean is similar or in some cases higher after truncation. Why? At least partially because after truncation the mean of the recent, living trees is of a large portion of the growth occurring during 20th century warming.
So, what do the chronologies look like if we use three different types of detrending? Let’s try constructing an RCS chronology using a 67% spline, negative exponential, and a generalized negative exponential curve:
The top shows the regional curve fit to all the raw Yamal data, the bottom panel zooms in on the part of the chronology since A.D. 1700 (heavy lines are 20 year low pass Butterworth filtered values). The regional curve fit that I think Jeff would endorse (67% spline) to fit his ‘U shape’ hypothesis actually results in slightly higher chronology values following the Little Ice Age, and slightly higher values at the end of the chronology in the mid 1990s.
My point is not to indicate a ‘correct’ method here — that goes well beyond ‘Blog Science’. Rather my point is this: detrending and standardization is one of the most challenging tasks in accurately estimating past low frequency climate variability from tree rings. Divergence is a serious challenge worthy of further study. What I don’t understand, however, is why methods or analyses (mean detrending, identifying growth curves from a small number of simultaneously growing trees) that can be shown to have the potential biases in specific instances that I’ve demonstrated here and in a previous post can still touted by those with a visceral dislike of paleoclimatologists as proof that another method is incorrect. If this debate was a collegial one that might be one thing, but there is nothing collegial — or scientific — about the language and tone from the other side. Too bad. To endorse once again Rob Wilson’s comment from Climate Audit:
“In fact, the fatal flaw in this blog and what keeps it from being a useful tool for the palaeoclimatic and other communities is its persistent and totally unnecessary negative tone and attitude, and the assumption that our intention is faulty and biased, which keeps real discourse from taking place.”
UPDATE: Minor grammar corrections
Keith Briffa and Tom Melvin have posted an interesting and thorough examination of the Yamal data here:
This now supersedes much if not all additional analysis I had considered for possible future posts.
Via Deep Climate, I found this post by Jeff Id at The Air Vent. Comments there and elsewhere lead me to believe there is some confusion about the related question of regional curve standardization and the reason for the importance of sample size in dendrochronology — dendroclimatology in particular — and while this post is only an indirect commentary on Jeff Id’s post, hopefully it will be more broadly useful or stimulate some interesting technical discussion. For more information on regional curve standardization, this book chapter [PDF] is currently your best bet.
Jeff Id fits two separate exponential growth curves to the most recent 12 trees in the Yamal chronology and to the full Yamal series, and notes that they are different. Let’s emulate this here. Let me note first of all that this is an emulation — the published Yamal series uses a time-varying spline fit that I haven’t integrated in my own code.
What I’ve done is align the full Yamal set (blue) and the most recent 12 ring width series (red) by age, assuming no pith offset (that is, assuming the innermost ring in the core or cross section was the innermost ring in the tree). The heavy lines are the mean regional curves. The black line is the Khadyta River mean regional curve.
There are few interesting features, some of them I believe are noted by Jeff in his post. The more recent Yamal trees had a somewhat lower growth rate when they were young than the average of the full set of living and subfossil trees; however, it is not at all out of the range of the full Yamal population. The regional curve for just these twelve is therefore lower than that of the regional curve for the full population. Another feature to note is that the red line (the regional curve for the 12 trees alone) rises at 100 years of age and again near 300 years of age, since these represent the ages of the most recent wider ring widths of several of the individual samples in this small set. Finally, note that the Khadyta River trees as a whole are relatively young, and their growth falls for the most part lower than the regional curves for the Yamal full and recent subsets.
UPDATE: I’m adding here the spline curve fits
So what is the consequence of performing a separate RCS on the recent Yamal series only vs. the full Yamal set?
The 12-series only chronology is somewhat noisier overall, since it also excludes 5 other tree ring series that come into the second half of the 20th century but not all the way into the 1990s. The influence of sample size on the chronology variance can also be seen in the 1600s, probably, when the year-to-year variability is reproduced but the small number of series (prior to 1660 or so, there are only 3 cores) influences the variance. The influence of the different regional curves — shown above — is more difficult to detect, since it is intermingled with the influence of the loss of the other 5 cores, but slightly higher levels in the 12-tree only chronology in the 1600s and parts of the 1700s might reflect it. The most notable difference is therefore perhaps that the recent-tree only chronology is slightly lower than the full Yamal chronology starting in the early 1980s.
So what is going on? In fact, you are witnessing the importance of overall sample size in the specific case of Regional Curve Standardization. It is important to understand the importance of sample replication for two different (but of course complementary) purposes in dendrochronology, specifically when applying Regional Curve Standardization:  Adequate sample replication overall so as to accurately estimate the ‘true’ regional growth curve, and  sample replication through time adequate to estimate the transient climate signal. Remember that the goal of regional curve standardization is to remove a common age-related growth trend while preserving low frequency climate variability — to have any hope of estimating this you need a large number of trees whose actual period of growth was well-distributed over time. The reason for this is that you need to avoid intermingling your climate signal of interest with your age-related growth trend. You can imagine an age-related growth trend estimated by trees of more or less the same age that grew more or less at the same time could intermingle the time-related environmental signals with the age-related geometric growth patterns. On the face of it then, Yamal is a good candidate for RCS since it has a large number of total trees whose actual time of growth is well-distributed over the length of the chronology. Isolating the 12 most recent trees, however, runs the risk of intermingling recent patterns of temperature variability with the trees’ common growth signal. The ‘regional curve’ from just these twelve trees is quite unlikely to be very representative of some significant fraction of the mean regional growth pattern associated with tree age.
The full Yamal regional growth curve is therefore likely to be a much better estimate of the ‘true’ regional growth curve common to trees from the region than a growth curve from a small number of trees growing over a period of anthropogenic climate change, because the climate signal of interest is a common feature of the growth of many of the trees. The lower chronology values in the recent-only chronology is red above is a consequence of at least part of the temperature signal being subtracted because it is intermingled when the regional curve is calculated over only a few trees growing, at the end of their life, in a warming world. Jeff’s post is a little hard to parse in places (for one thing, he keeps referring to ‘climatology’, but I think he mean ‘climatologists’ or the ‘climatology community’), but reading carefully, I think he might recognize this as a potential problem.
Now, the other important part of having multiple samples from the same site is maximize the signal to noise ratio (for our purposes, the signal is climate) at any given time. Dendrochronologists have ways of traditionally estimated whether their chronology is sufficiently well-replicated and contains a common signal, including the Expressed Population Signal or Subsample Signal Strength. Using 20 year windows with 10 years of overlap, the Expressed Population Signal for the 12-series only RCS chronology is consistently above the (arbitrary but historical) 0.85 level back to the 18th century (and, indeed, most of the way back to the earlier parts of the chronology before sample size is reduced to a few cores). For 10 year windows with a 5 year overlap (not something I would consider particularly stable, but it allows us to look at very small slice of time), the EPS exceeds 0.85 from at least 1990 back to beginning of the chronology with only two decades in the 19th and 18th centuries with low interseries correlation. Note that these windows are shorter than we normally use.
My take home message is this and it is intended to be general: it is important to understand the two complementary parts of the importance of samples size in developing RCS chronologies for climate reconstruction. Lots of ring width series are necessary to develop an accurate regional curve. The number of chronologies needed at any given point in time to capture the transient climate signal can be estimated using EPS. Strong average interseries correlation between cores can mean that even relative few trees collectively capture a significant portion — again, as estimated from established metrics — of the climate variance and allow for adequate signal to noise ratios in the mean chronology. Replication gives us increased confidence in the value of the mean chronology, but a strong common signal is an important part of the equation.
A fair number of people have been visiting this humble blog in the last few days, presumably not for my excellent tips on where to find the best huevos rancheros in the world. My previous posts on Yamal are long and perhaps technical in places, so I wanted to provide a quick summary.
in Yamal I, I emulated the modifications of the Yamal chronology by Steve McIntyre. As I noted in that post:
Adding the Khadyta River series reduces the the level of the chronology though the 1970s and 1980s and into the early 1990s, when those data end. But if one includes both data sets, the series terminates similarly to the original Yamal chronology, of course (because the last few years are only present in the modern trees from Yamal). These changes are potentially important, and the actual scientific questions are interesting
I also noted that the manner in which the series are low pass filtered can change the impression of the impact of these changes. Smoothing and filtering series — particularly dealing with end points — is a complicated issue, so care is always necessary when designing, deploying, and subsequently interpreting smoothed series.
In doing the emulation, I noted that there are differences in the archived Yamal, McIntyre emulation, and my own. These probably arise due to the different curve fitting techniques (time varying spline for Yamal, nonlinear least squares for McIntyre’s R code, and standard spline or Supersmoother fitting for my emulation). In terms of identifying relative differences due to adding the Khadyta River data, though, these don’t seem to matter too much. Apparently, McIntyre noted that his RCS routine created a greater rise than the archived Yamal chronology (again, this is probably differences in curve fitting, if I had to guess), but I can’t locate the statement now.
I also see, via Deep Climate, that McIntyre has posted graphs showing the chronology from the full data set (including the complete set of trees from Yamal and Khadyta River). We seem to differ on showing the period from 1991 to 1996, apparently — I think if you are going to make rather strong statements about using all the data, you should use all the data, but I don’t expect us to agree on this.
So, more interesting perhaps is my second post, Yamal II, where I looked at the raw data from Yamal and Khadyta River. Here, I concluded:
Khadyta River does display the divergence problem, in that it ceases to track temperatures as it did from 1883 to the 1960s.
Relatively little more can be said about the specific case of Yamal at this point — I’ll leave that to the scientists working in this part of the world — but even my quick review of these data here shows that including Khadyta River raw data in the Yamal chronology does not result in a more accurate nor precise understanding of past temperatures in the region. This isn’t to say that some time in the past that Yamal didn’t experience divergence (this after all is a large part of the concern about divergence), but we can clearly see that Khadyta River does exhibit modern divergence.
That’s where I stand now, and I’m not sure what further insights will come from solely the data at hand. However, reading blog posts about this topic, I’ve become aware that there persist misunderstandings in how field and lab dendrochronology is actually done. So, hopefully there will a post about this from me in the near future.