the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 17 May 2021
Evidence for old carbon contamination in 14C wiggle-match age series for the 946 CE eruption of Changbaishan volcano
Abstract. Volcanic eruptions that are not historically attested are commonly radiocarbon dated by "wiggle matching" sequential 14C measurements of the rings of trees killed by the eruption against an accepted calibration curve. It is generally assumed that carbon laid down in the wood is uncontaminated by 14C-free ("old") carbon, although evidence for contamination is well documented. Often, ill-fitting ring ages are excluded from analysis. The ‘Millennium Eruption' of Changbaishan volcano on the China-DPR Korea border offers a valuable case study in wiggle match dating, since several independent groups reported age estimates before the determination and acceptance of a precise eruption year of 946 CE. Some of the discrepancies and incompatibilities between published dates were attributed to old carbon effects. Here, we apply a new methodology to correct for contamination levels of up to 4.5% old carbon to eight wiggle match date series for the Millennium Eruption. Without discarding ring ages, we find agreement indices as high as, or higher than, those for the published dates, and five of the eight date series yielded high-agreement-index eruption dates closer to 946 CE than the published dates. None of the five yield a best result at zero contamination. Differences between the eruption dates reveal a weak association with the direction of the sampled tree from the caldera, but no relationship with distance. Our results suggest that old carbon contamination is possible over a wide area, potentially leading to over-estimation of eruption ages by years, decades or more, cautioning against over-reliance on wiggle-match ages that are not corroborated by other lines of evidence. Our revised protocol that accounts for contamination offers a way forward in the application of wiggle match dating of eruptions and provides a platform for understanding discrepancies that exist when comparing wiggle match series.
- Preprint
(2460 KB) - Metadata XML
-
Supplement
(2564 KB) - BibTeX
- EndNote
Status: closed
-
RC1: 'Comment on gchron-2021-13', Anonymous Referee #1, 01 Aug 2021
Summary: This is a generally well written manuscript which makes interesting use of the calendar dating of the Changbaishan Millennium Eruption to attempt to test assumptions behind previous wigglematching of radiocarbon samples for the same event in order to evaluate the potential impacts of “old” carbon contamination. It is an interesting use of these published data to address an important question - the potential impact of volcanic contamination on wigglematching - and also a good way to highlight the differences in approach of using radiocarbon in tree rings to date volcanoes. There is some nice detailed content on the general topic of volcanic old carbon contamination. For all these reasons I'd really like to see this further developed towards publication, but I do have a series of fairly significant concerns that mean that I just can't recommend it for publication at this stage.
These are my main concerns:
- The published data points are only described in terms of volcanic contamination and this is not sufficient because there are many other interlaboratory or sampling based factors (numbers of rings in each sample, length of time represented, overall age of trees as represented by tree-rings present, changes in the shape of callibration curve relative to these) which could similarly explain apparent outliers. This at the very least should be addressed / discussed and considered within the methodological / decision making processes involved. Without this the whole aim as stated in the title seems fundamentally flawed. How can this study provide "Evidence for old carbon contamination in 14C wiggle-match age series for the 946 CE eruption of Changbaishan volcano" if these other possible explanations are not discussed and ruled out in the manuscript.
- Figure 2 "Confidence intervals (95.4%) for wiggle match dates for the Changbaishan Millennium eruption, in relation to the
946 CE date for the eruption (vertical broken line)" is very nice and a logical step in the analysis, but my interpretation of this is that in fact the majority of the previous wigglematches look pretty good with the known calendar date. Yes there are some probability ranges which do have much older possibilities, and some slight off-sets towards older but on the whole this figure gives me confidence in the wigglematch method which can not hope to offer the sort of precision and accuracy that securing 14C on a cosmogenic anchor point can. Generally most studies got somewhere quite close to, or in overlap with, the actual date. The whole 95.4% range needs to be considered. The slightly older age extents for some of the samples could also be explained by interlab variations and differences in models or calibration curves. On page 15, line 12, the authors state: "For example, an earlier wiggle match date of 938 +8/-5 CE, had been favoured from external evidence from dendrochronology (937–938 CE) and varves (912–972 CE) (Nakamura et al., 2007a) until further wiggle match series were measured (Xu et al., 2013).." - Making the point I think that the earlier estimate was somehow wrong, but this fails to recognize that 938 +8/-5 CE is in fact spot on for the later confirmed calendar date of 946 CE.... For such very slight differences in carbon that appear to being evaluated here there are a multitude of possible explanations, I would have thought, based on what we know of volcanic contamination in other studies, the smoking gun for volcanic contamination would be older dates of c100 or so year or more, not more subtle variations, but here we seem to be focusing on possibilities that are very slight.(What particularly stands out in the figure to me are the younger outliers which are not explained by old carbon contamination. I like that these are addressed in the text in terms of uptake of younger carbon, but here again but a wide range of other possibilities (like the geological interpretation of the context for the sample) are not really explored (though it seems like the later spatial information might also be very relevant here)... But this is a side issue... something to think about if reframing and refocusing this text for resubmission)
- All other comments aside, I would not feel that this approach could be applied with confidence in other studies especially given that the authors note that it can only be applied where a consistent level of contamination can be assumed over the life of the tree. Volcanic eruptions and the potential for volcanic gas release and other factors which might impact storage of carbon in sequential tree-ring series overtime are highly variable and very site/eruption specific. Also factors such as growth season of the trees, the period of time covered by the tree-ring sequence involved and the fact that old carbon may be released in some years and not in others all combines to complicate the application of a blanket correction factor for a given study. This sort of blanket correction could equally apply to correct a lab-off set and this would be nothing new.
A few other more specific comments to help with revisions
P1: Lines 1-2: Suggest simplification from: "increasing accuracy as the technologies developed, more refined calibration curves promulgated (Hogg et al., 2020; Reimer et al., 2020), and sophistication of statistical analyses (Crema and Bevan, 2020) to: "increasing accuracy with the use of more refined calibration curves (Hogg et al., 2020; Reimer et al., 2020), and sophisticated statistical analyses (Crema and Bevan, 2020)"
P2: Lines 5-6: Suggest change from " as the most likely to give a calendar date " to " as the most likely to give a close to
calendar date " (calendar date implies actual calendar date which is only possible with a 14C excursion such as the 774/5 CE event, not through conventional radiocarbon wigglematching, so either make this more subtle change or expand a bit further on making the distinction between the two different ways of using radiocarbon to date samples from eruption contexts, this may be helpful to the reader at this stage as the study rather depends on people appreciating the distinction.)P4: Line 27: "cosmogenically-tuned date" - change to "cosmogenically secured" important for readers to understand that this 'radiocarbon date' based on the 774/5 CE event IS a calendar date not subject to variations like wigglematches always will be with new iterations of the curve. Also on line 27 you introduce 3 possible options for why previous radiocarbon dates might be older than the date for the eruption given by Oppenheimer / Hakozaki and then go into a nicely comprehensive list of possible volcanic causes... but you do not introduce the significant range of other possibilities to explain this - for example, were Sun et al's dates based on a sequence of older to younger samples from a single tree? If they were, how old was the tree in total and were the samples blocks of years? Half the dates being older in a consecutive sequence from a tree would be predicted, especially if that tree was longer lived. Across how many years is the 14C of a single sample averaged to get the values for the wigglematched blocks and how much does the calibration curve vary (or how much could it vary around the time of the event (i.e. curve plateaus or slopes) - these are all really important to address and provide a thorough framing for the 'contamination' discussion..
P6: Figure 2, line 3, 774 not 74
Citation: https://doi.org/10.5194/gchron-2021-13-RC1 -
AC1: 'Reply on RC1', Richard Holdaway, 29 Sep 2021
Reviewer 1
Summary: This is a generally well written manuscript which makes interesting use of the calendar dating of the Changbaishan Millennium Eruption to attempt to test assumptions behind previous wigglematching of radiocarbon samples for the same event in order to evaluate the potential impacts of “old” carbon contamination. It is an interesting use of these published data to address an important question - the potential impact of volcanic contamination on wigglematching - and also a good way to highlight the differences in approach of using radiocarbon in tree rings to date volcanoes. There is some nice detailed content on the general topic of volcanic old carbon contamination. For all these reasons I'd really like to see this further developed towards publication, but I do have a series of fairly significant concerns that mean that I just can't recommend it for publication at this stage.
***We appreciate the detailed, thoughtful and overall positive review of our MS, and are grateful for the time and effort involved. We respond to their concerns in detail below.
These are my main concerns:
- The published data points are only described in terms of volcanic contamination and this is not sufficient because there are many other interlaboratory or sampling based factors (numbers of rings in each sample, length of time represented, overall age of trees as represented by tree-rings present, changes in the shape of callibration curve relative to these) which could similarly explain apparent outliers. This at the very least should be addressed / discussed and considered within the methodological / decision making processes involved. Without this the whole aim as stated in the title seems fundamentally flawed. How can this study provide "Evidence for old carbon contamination in 14C wiggle-match age series for the 946 CE eruption of Changbaishan volcano" if these other possible explanations are not discussed and ruled out in the manuscript.
***We agree that the absence of discussion of other potential problems with the ages in wiggle match series is a weakness in our MS and will rectify that omission in a revised version. We note in passing that such concerns are rarely dealt with in conventional wiggle match analyses: we agree that they should be addressed in any such study. Our aim was not to explain “apparent outliers” based on the assumption that no contamination affects wiggle match series, but to demonstrate that such series may better fit other sections of a calibration curve if some level of contamination – whatever the source but especially magmatic/geologic carbon – was present. We note that the effects listed also affect conventional age series but are seldom discussed except when invoking a method for identifying and rejecting ages that do not fit the assumed section of the relevant curve. Standard analyses concentrate on forcing a fit (a statistical Procrustes process) to what may not be an appropriate section of the curve.
- Figure 2 "Confidence intervals (95.4%) for wiggle match dates for the Changbaishan Millennium eruption, in relation to the 946 CE date for the eruption (vertical broken line)" is very nice and a logical step in the analysis, but my interpretation of this is that in fact the majority of the previous wigglematches look pretty good with the known calendar date. Yes there are some probability ranges which do have much older possibilities, and some slight off-sets towards older but on the whole this figure gives me confidence in the wigglematch method which can not hope to offer the sort of precision and accuracy that securing 14C on a cosmogenic anchor point can. Generally most studies got somewhere quite close to, or in overlap with, the actual date. The whole 95.4% range needs to be considered.
***We agree with the reviewer’s point, but wish to emphasise that our Fig. 2 is a post hoc representation of the extremely unusual situation of multiple wiggle matches and an independently derived date for an eruption. The main aim of our paper is to show that there can be other, even better, fits of wiggle match age series to a calibration curve which need to be considered if the usual assumption of no contamination from any source is set aside. We agree that standard wiggle matches cannot hope to isolate an exact eruption date that other, independent, dating might provide. We contend that it is possible to get better agreement (within error bars) between wiggle matches series and a known date, when some level of contamination is allowed. In the (usual) absence of a known eruption date, our method of examining the effects of hypothesised non-zero contamination should allow better agreement between wiggle match ages. It provides a systematic methodology for exploration of possible contamination from any source and as a basis for explaining disagreement in ages..
The slightly older age extents for some of the samples could also be explained by interlab variations and differences in models or calibration curves. On page 15, line 12, the authors state: "For example, an earlier wiggle match date of 938 +8/-5 CE, had been favoured from external evidence from dendrochronology (937–938 CE) and varves (912–972 CE) (Nakamura et al., 2007a) until further wiggle match series were measured (Xu et al., 2013).." - Making the point I think that the earlier estimate was somehow wrong, but this fails to recognize that 938 +8/-5 CE is in fact spot on for the later confirmed calendar date of 946 CE....
***We agree that the Xu et al date is “spot on”, we will emphasize this in the new text.But we will also emphasize that this is the last in a series of wiggle match attempts, some of which, although accepted for publication, can now be seen to have been decidedly inaccurate.
For such very slight differences in carbon that appear to being evaluated here there are a multitude of possible explanations, I would have thought, based on what we know of volcanic contamination in other studies, the smoking gun for volcanic contamination would be older dates of c100 or so year or more, not more subtle variations, but here we seem to be focusing on possibilities that are very slight.
***We agree with the reviewer, larger eruptions are more likely to show large systematic offsets (Holdaway et al. 2018). However we are unaware of any such eruptions with an independently known date and multiple wiggle match trees by which to test the hypothesis. The discussion is only possible at Changbaishan because multiple wiggle match series for the eruption and the systematic methodology shows onlyslight effects which might be expected in an open forest on a windy mountainside. Our intention was to introduce a methodology grounded on a well-constrained example that can be used by the community to explore the possibility of alternative, perhaps better, fits to a calibration curve. We feel that it is premature to set an arbitrary “smoking gun” value such as 100 years (c. 1.2% old carbon contamination) when such contamination can conceivably affect fits at any level above 0%.
***Contamination from any source can conceivably vary from 0 through “very slight” (0.25 – 1%?) to considerable (>1%) which can occur at different times during the life of the tree leading up to the time of the eruption which is inferred to have killed the tree. Changing levels of contamination could be enough to distort a wiggle match series and create a false match and this is a possibility that must be faced. At present, the standard approach is to seek out and remove ages that compromise the fit at the assumed (no contamination) section of the contamination curve.
(What particularly stands out in the figure to me are the younger outliers which are not explained by old carbon contamination. I like that these are addressed in the text in terms of uptake of younger carbon, but here again but a wide range of other possibilities (like the geological interpretation of the context for the sample) are not really explored (though it seems like the later spatial information might also be very relevant here)... But this is a side issue... something to think about if reframing and refocusing this text for resubmission)
***We agree that the “young” eruption dates are problematic and will expand the discussion on that issue. We note that these are not “outliers” in the conventional sense: each whole age series provides a match with the then current calibration curve that was accepted before later wiggle matches were attempted. The “young dates” resulted in the former name “Millennial Eruption”.
- All other comments aside, I would not feel that this approach could be applied with confidence in other studies especially given that the authors note that it can only be applied where a consistent level of contamination can be assumed over the life of the tree.
***We do not suggest that the method can be applied only where a consistent level of contamination is assumed. We propose the method recognising that by applying a constant contamination is better than applying no contamination at all (a fact which is not recognised in conventional wiggle match dating) as an essential first step in assessing a wiggle match series to see if the assumption of no contamination, consistent or otherwise, is valid. Where there are better matches using the constant contamination assumption then surely the possibility of contamination should be considered. We say that, and will make it clearer in a revised MS that we simplify the analysis by assuming a constant rate. This is analogous to the present situation where no contamination is assumed and only a single fit sought. The number of wiggle match ages is important here, as series of only 4 or 5 ages can conceivably fit several parts of a curve. By introducing the method with this assumption we hope more studies follow that explore non constant contamination models and suspect these will improve dating accuracy further..
Volcanic eruptions and the potential for volcanic gas release and other factors which might impact storage of carbon in sequential tree-ring series overtime are highly variable and very site/eruption specific. Also factors such as growth season of the trees, the period of time covered by the tree-ring sequence involved and the fact that old carbon may be released in some years and not in others all combines to complicate the application of a blanket correction factor for a given study. This sort of blanket correction could equally apply to correct a lab-off set and this would be nothing new.
***We agree, but this is seldom if ever done at the moment, and suggest our methodology would be equally useful for identifying and adjusting for interlab off set.. The comment serves to emphasise what our MS proposes, that wiggle match dating of eruptions (or any other event where systematic contamination/offset might occur) is not simple and investigation of potential methods of refining the analyses are necessary, rather than a statistically Procrustean forced fit by ad hoc exclusion of ages that compromise a no-contamination fit. It is certainly an uncomfortable prospect but one which cannot be ignored.
A few other more specific comments to help with revisions
P1: Lines 1-2: Suggest simplification from: "increasing accuracy as the technologies developed, more refined calibration curves promulgated (Hogg et al., 2020; Reimer et al., 2020), and sophistication of statistical analyses (Crema and Bevan, 2020) to: "increasing accuracy with the use of more refined calibration curves (Hogg et al., 2020; Reimer et al., 2020), and sophisticated statistical analyses (Crema and Bevan, 2020)"
***We will consider the suggested simplification.
P2: Lines 5-6: Suggest change from " as the most likely to give a calendar date " to " as the most likely to give a close to calendar date " (calendar date implies actual calendar date which is only possible with a 14C excursion such as the 774/5 CE event, not through conventional radiocarbon wigglematching, so either make this more subtle change or expand a bit further on making the distinction between the two different ways of using radiocarbon to date samples from eruption contexts, this may be helpful to the reader at this stage as the study rather depends on people appreciating the distinction.)
***We take the reviewers’ point and will expand the text to make the distinction clearer.
P4: Line 27: "cosmogenically-tuned date" - change to "cosmogenically secured" important for readers to understand that this 'radiocarbon date' based on the 774/5 CE event IS a calendar date not subject to variations like wigglematches always will be with new iterations of the curve. Also on line 27 you introduce 3 possible options for why previous radiocarbon dates might be older than the date for the eruption given by Oppenheimer / Hakozaki and then go into a nicely comprehensive list of possible volcanic causes... but you do not introduce the significant range of other possibilities to explain this - for example, were Sun et al's dates based on a sequence of older to younger samples from a single tree? If they were, how old was the tree in total and were the samples blocks of years? Half the dates being older in a consecutive sequence from a tree would be predicted, especially if that tree was longer lived. Across how many years is the 14C of a single sample averaged to get the values for the wigglematched blocks and how much does the calibration curve vary (or how much could it vary around the time of the event (i.e. curve plateaus or slopes) - these are all really important to address and provide a thorough framing for the 'contamination' discussion..
***We will explore the implications as to other possible causes (a point raised also by Reviewer 2) and will alter the text accordingly.
P6: Figure 2, line 3, 774 not 74
***We will correct the typo.
Citation: https://doi.org/10.5194/gchron-2021-13-AC1
-
AC1: 'Reply on RC1', Richard Holdaway, 29 Sep 2021
-
RC2: 'Comment on gchron-2021-13', Anonymous Referee #2, 08 Sep 2021
This is a really interesting article that also raises an interesting point that there is at least the potential that trees being dated by Dendro/14C wiggle matching may be influenced by gological dead 14C. As such the paper should be published. It is an important but fundamentally overlooked point. I do, however, have a number of cocernes that need in my mind to be addressed.
My first issue is one that is also raised by the first refereee, in that there are many sources of contamination that can cuase problems with 14C ages, and my not be removed by standard techniques. There is clearly a problem that is unlikely to be related to volcanic CO2 in the ages that are too young, as opposed to being too old. This needs to be discussed in more detail. In fact I think the paper needs to more be an examination of the problems that can be encountered with 14C ages and placed in the context of wiggle match dating of eruptions. In this unusual case we know the age very specifically and that allows a discussion of these problems. Some of the peroblems with the dates relates to contamination by older 14C, which may be from volcanic degassing, but could be from a range of other issues. Some dates have been contaminated by young 14C and this is unlikely to be a volcanic origin. One way of looking for potential contamination sources is through delta 13C measurements. In cases where 14C and 13C of volcanic emissions have been considered on modern trees there is evidence of a shift in both isotopes. The fact that it is pretty difficult to find 13C measurements for most of the cases used in this study makes the point that there are a whole host of quality assurance issues in raduiocarbon dating of volcanic eruptions that need to be considered. The known age of the B-TM makes this a great paper around which to have that wider discussion. This is moreso becuase I do not see how, unless we know the exact age of the eruption, we can use this approach, to finding old carbon. Also even is we can find the percentage of old carbon in a sample, withouot additional measurements such as 13C we can't start to suggest the source of the old carbon with much reliability.
I also wasn't sure why the outlier function in OxCal wasn't used instead as this has largely superceedeed the use of AI's in OxCal outlier detection studies. The AI is simply the area of overlap between posterior and prior. I would like to see at the very least the analyses done again with appropriate outlier detection.
Citation: https://doi.org/10.5194/gchron-2021-13-RC2 -
AC2: 'Reply on RC2', Richard Holdaway, 29 Sep 2021
Reviewer 2
This is a really interesting article that also raises an interesting point that there is at least the potential that trees being dated by Dendro/14C wiggle matching may be influenced by gological dead 14C. As such the paper should be published. It is an important but fundamentally overlooked point. I do, however, have a number of cocernes that need in my mind to be addressed.
***We appreciate the reviewer’s view that our paper deals with “an important but fundamentally overlooked point” and are pleased with their recommendation. We respond to their concerns in detail below.
My first issue is one that is also raised by the first refereee, in that there are many sources of contamination that can cuase problems with 14C ages, and my not be removed by standard techniques. There is clearly a problem that is unlikely to be related to volcanic CO2 in the ages that are too young, as opposed to being too old. This needs to be discussed in more detail. In fact I think the paper needs to more be an examination of the problems that can be encountered with 14C ages and placed in the context of wiggle match dating of eruptions.
***As in our response to this point as raised by Reviewer 1, we will expand our discussion of other potential sources of contamination while confirming our position that, for wiggle match series measured on trees within the influence of a volcanic edifice, geologic/magmatic carbon is an unappreciated but potentially dominant source. And one which can provide enough contaminating carbon to actually change the position of best fit on a calibration curve.
In this unusual case we know the age very specifically and that allows a discussion of these problems. Some of the peroblems with the dates relates to contamination by older 14C, which may be from volcanic degassing, but could be from a range of other issues. Some dates have been contaminated by young 14C and this is unlikely to be a volcanic origin. One way of looking for potential contamination sources is through delta 13C measurements. In cases where 14C and 13C of volcanic emissions have been considered on modern trees there is evidence of a shift in both isotopes. The fact that it is pretty difficult to find 13C measurements for most of the cases used in this study makes the point that there are a whole host of quality assurance issues in raduiocarbon dating of volcanic eruptions that need to be considered.
***We agree that δ13C is a potentially useful marker of the source of any contaminating carbon. We have explored, and continue to explore, the possibility of δ13C evidence for the source(s) of contamination. As the reviewer notes, such data are rarely reported and we very much agree that there are issues of quality assurance in radiocarbon dating of eruptions.
The known age of the B-TM makes this a great paper around which to have that wider discussion. This is moreso becuase I do not see how, unless we know the exact age of the eruption, we can use this approach, to finding old carbon. Also even is we can find the percentage of old carbon in a sample, withouot additional measurements such as 13C we can't start to suggest the source of the old carbon with much reliability.
***We agree and certainly hope that our paper could be a focus for a wider, constructive, discussion on the issues which undoubtedly surround the wiggle match process. We propose our method as a way of identifying the potential presence of old carbon in a wiggle match tree, by exploring the possibility of the existence ofother, conceivably better, matches if zero contamination is not assumed.
I also wasn't sure why the outlier function in OxCal wasn't used instead as this has largely superceedeed the use of AI's in OxCal outlier detection studies. The AI is simply the area of overlap between posterior and prior. I would like to see at the very least the analyses done again with appropriate outlier detection.
***We investigated but did not pursue the use of the “outlier” function in OxCal. We think that this is basically another way of isolating ages that don’t fit with an assumed position of fit on the curve. A different way of applying the Procrustean approach to a wiggle match fit. As such, we do not believe that its application would be appropriate here. We used AIs because they are recognised numerical method of assessing overall fit of a wiggle match series to a calibration curve. They allow fits at different sections of the curve to be compared objectively.
Citation: https://doi.org/10.5194/gchron-2021-13-AC2
-
AC2: 'Reply on RC2', Richard Holdaway, 29 Sep 2021
Status: closed
-
RC1: 'Comment on gchron-2021-13', Anonymous Referee #1, 01 Aug 2021
Summary: This is a generally well written manuscript which makes interesting use of the calendar dating of the Changbaishan Millennium Eruption to attempt to test assumptions behind previous wigglematching of radiocarbon samples for the same event in order to evaluate the potential impacts of “old” carbon contamination. It is an interesting use of these published data to address an important question - the potential impact of volcanic contamination on wigglematching - and also a good way to highlight the differences in approach of using radiocarbon in tree rings to date volcanoes. There is some nice detailed content on the general topic of volcanic old carbon contamination. For all these reasons I'd really like to see this further developed towards publication, but I do have a series of fairly significant concerns that mean that I just can't recommend it for publication at this stage.
These are my main concerns:
- The published data points are only described in terms of volcanic contamination and this is not sufficient because there are many other interlaboratory or sampling based factors (numbers of rings in each sample, length of time represented, overall age of trees as represented by tree-rings present, changes in the shape of callibration curve relative to these) which could similarly explain apparent outliers. This at the very least should be addressed / discussed and considered within the methodological / decision making processes involved. Without this the whole aim as stated in the title seems fundamentally flawed. How can this study provide "Evidence for old carbon contamination in 14C wiggle-match age series for the 946 CE eruption of Changbaishan volcano" if these other possible explanations are not discussed and ruled out in the manuscript.
- Figure 2 "Confidence intervals (95.4%) for wiggle match dates for the Changbaishan Millennium eruption, in relation to the
946 CE date for the eruption (vertical broken line)" is very nice and a logical step in the analysis, but my interpretation of this is that in fact the majority of the previous wigglematches look pretty good with the known calendar date. Yes there are some probability ranges which do have much older possibilities, and some slight off-sets towards older but on the whole this figure gives me confidence in the wigglematch method which can not hope to offer the sort of precision and accuracy that securing 14C on a cosmogenic anchor point can. Generally most studies got somewhere quite close to, or in overlap with, the actual date. The whole 95.4% range needs to be considered. The slightly older age extents for some of the samples could also be explained by interlab variations and differences in models or calibration curves. On page 15, line 12, the authors state: "For example, an earlier wiggle match date of 938 +8/-5 CE, had been favoured from external evidence from dendrochronology (937–938 CE) and varves (912–972 CE) (Nakamura et al., 2007a) until further wiggle match series were measured (Xu et al., 2013).." - Making the point I think that the earlier estimate was somehow wrong, but this fails to recognize that 938 +8/-5 CE is in fact spot on for the later confirmed calendar date of 946 CE.... For such very slight differences in carbon that appear to being evaluated here there are a multitude of possible explanations, I would have thought, based on what we know of volcanic contamination in other studies, the smoking gun for volcanic contamination would be older dates of c100 or so year or more, not more subtle variations, but here we seem to be focusing on possibilities that are very slight.(What particularly stands out in the figure to me are the younger outliers which are not explained by old carbon contamination. I like that these are addressed in the text in terms of uptake of younger carbon, but here again but a wide range of other possibilities (like the geological interpretation of the context for the sample) are not really explored (though it seems like the later spatial information might also be very relevant here)... But this is a side issue... something to think about if reframing and refocusing this text for resubmission)
- All other comments aside, I would not feel that this approach could be applied with confidence in other studies especially given that the authors note that it can only be applied where a consistent level of contamination can be assumed over the life of the tree. Volcanic eruptions and the potential for volcanic gas release and other factors which might impact storage of carbon in sequential tree-ring series overtime are highly variable and very site/eruption specific. Also factors such as growth season of the trees, the period of time covered by the tree-ring sequence involved and the fact that old carbon may be released in some years and not in others all combines to complicate the application of a blanket correction factor for a given study. This sort of blanket correction could equally apply to correct a lab-off set and this would be nothing new.
A few other more specific comments to help with revisions
P1: Lines 1-2: Suggest simplification from: "increasing accuracy as the technologies developed, more refined calibration curves promulgated (Hogg et al., 2020; Reimer et al., 2020), and sophistication of statistical analyses (Crema and Bevan, 2020) to: "increasing accuracy with the use of more refined calibration curves (Hogg et al., 2020; Reimer et al., 2020), and sophisticated statistical analyses (Crema and Bevan, 2020)"
P2: Lines 5-6: Suggest change from " as the most likely to give a calendar date " to " as the most likely to give a close to
calendar date " (calendar date implies actual calendar date which is only possible with a 14C excursion such as the 774/5 CE event, not through conventional radiocarbon wigglematching, so either make this more subtle change or expand a bit further on making the distinction between the two different ways of using radiocarbon to date samples from eruption contexts, this may be helpful to the reader at this stage as the study rather depends on people appreciating the distinction.)P4: Line 27: "cosmogenically-tuned date" - change to "cosmogenically secured" important for readers to understand that this 'radiocarbon date' based on the 774/5 CE event IS a calendar date not subject to variations like wigglematches always will be with new iterations of the curve. Also on line 27 you introduce 3 possible options for why previous radiocarbon dates might be older than the date for the eruption given by Oppenheimer / Hakozaki and then go into a nicely comprehensive list of possible volcanic causes... but you do not introduce the significant range of other possibilities to explain this - for example, were Sun et al's dates based on a sequence of older to younger samples from a single tree? If they were, how old was the tree in total and were the samples blocks of years? Half the dates being older in a consecutive sequence from a tree would be predicted, especially if that tree was longer lived. Across how many years is the 14C of a single sample averaged to get the values for the wigglematched blocks and how much does the calibration curve vary (or how much could it vary around the time of the event (i.e. curve plateaus or slopes) - these are all really important to address and provide a thorough framing for the 'contamination' discussion..
P6: Figure 2, line 3, 774 not 74
Citation: https://doi.org/10.5194/gchron-2021-13-RC1 -
AC1: 'Reply on RC1', Richard Holdaway, 29 Sep 2021
Reviewer 1
Summary: This is a generally well written manuscript which makes interesting use of the calendar dating of the Changbaishan Millennium Eruption to attempt to test assumptions behind previous wigglematching of radiocarbon samples for the same event in order to evaluate the potential impacts of “old” carbon contamination. It is an interesting use of these published data to address an important question - the potential impact of volcanic contamination on wigglematching - and also a good way to highlight the differences in approach of using radiocarbon in tree rings to date volcanoes. There is some nice detailed content on the general topic of volcanic old carbon contamination. For all these reasons I'd really like to see this further developed towards publication, but I do have a series of fairly significant concerns that mean that I just can't recommend it for publication at this stage.
***We appreciate the detailed, thoughtful and overall positive review of our MS, and are grateful for the time and effort involved. We respond to their concerns in detail below.
These are my main concerns:
- The published data points are only described in terms of volcanic contamination and this is not sufficient because there are many other interlaboratory or sampling based factors (numbers of rings in each sample, length of time represented, overall age of trees as represented by tree-rings present, changes in the shape of callibration curve relative to these) which could similarly explain apparent outliers. This at the very least should be addressed / discussed and considered within the methodological / decision making processes involved. Without this the whole aim as stated in the title seems fundamentally flawed. How can this study provide "Evidence for old carbon contamination in 14C wiggle-match age series for the 946 CE eruption of Changbaishan volcano" if these other possible explanations are not discussed and ruled out in the manuscript.
***We agree that the absence of discussion of other potential problems with the ages in wiggle match series is a weakness in our MS and will rectify that omission in a revised version. We note in passing that such concerns are rarely dealt with in conventional wiggle match analyses: we agree that they should be addressed in any such study. Our aim was not to explain “apparent outliers” based on the assumption that no contamination affects wiggle match series, but to demonstrate that such series may better fit other sections of a calibration curve if some level of contamination – whatever the source but especially magmatic/geologic carbon – was present. We note that the effects listed also affect conventional age series but are seldom discussed except when invoking a method for identifying and rejecting ages that do not fit the assumed section of the relevant curve. Standard analyses concentrate on forcing a fit (a statistical Procrustes process) to what may not be an appropriate section of the curve.
- Figure 2 "Confidence intervals (95.4%) for wiggle match dates for the Changbaishan Millennium eruption, in relation to the 946 CE date for the eruption (vertical broken line)" is very nice and a logical step in the analysis, but my interpretation of this is that in fact the majority of the previous wigglematches look pretty good with the known calendar date. Yes there are some probability ranges which do have much older possibilities, and some slight off-sets towards older but on the whole this figure gives me confidence in the wigglematch method which can not hope to offer the sort of precision and accuracy that securing 14C on a cosmogenic anchor point can. Generally most studies got somewhere quite close to, or in overlap with, the actual date. The whole 95.4% range needs to be considered.
***We agree with the reviewer’s point, but wish to emphasise that our Fig. 2 is a post hoc representation of the extremely unusual situation of multiple wiggle matches and an independently derived date for an eruption. The main aim of our paper is to show that there can be other, even better, fits of wiggle match age series to a calibration curve which need to be considered if the usual assumption of no contamination from any source is set aside. We agree that standard wiggle matches cannot hope to isolate an exact eruption date that other, independent, dating might provide. We contend that it is possible to get better agreement (within error bars) between wiggle matches series and a known date, when some level of contamination is allowed. In the (usual) absence of a known eruption date, our method of examining the effects of hypothesised non-zero contamination should allow better agreement between wiggle match ages. It provides a systematic methodology for exploration of possible contamination from any source and as a basis for explaining disagreement in ages..
The slightly older age extents for some of the samples could also be explained by interlab variations and differences in models or calibration curves. On page 15, line 12, the authors state: "For example, an earlier wiggle match date of 938 +8/-5 CE, had been favoured from external evidence from dendrochronology (937–938 CE) and varves (912–972 CE) (Nakamura et al., 2007a) until further wiggle match series were measured (Xu et al., 2013).." - Making the point I think that the earlier estimate was somehow wrong, but this fails to recognize that 938 +8/-5 CE is in fact spot on for the later confirmed calendar date of 946 CE....
***We agree that the Xu et al date is “spot on”, we will emphasize this in the new text.But we will also emphasize that this is the last in a series of wiggle match attempts, some of which, although accepted for publication, can now be seen to have been decidedly inaccurate.
For such very slight differences in carbon that appear to being evaluated here there are a multitude of possible explanations, I would have thought, based on what we know of volcanic contamination in other studies, the smoking gun for volcanic contamination would be older dates of c100 or so year or more, not more subtle variations, but here we seem to be focusing on possibilities that are very slight.
***We agree with the reviewer, larger eruptions are more likely to show large systematic offsets (Holdaway et al. 2018). However we are unaware of any such eruptions with an independently known date and multiple wiggle match trees by which to test the hypothesis. The discussion is only possible at Changbaishan because multiple wiggle match series for the eruption and the systematic methodology shows onlyslight effects which might be expected in an open forest on a windy mountainside. Our intention was to introduce a methodology grounded on a well-constrained example that can be used by the community to explore the possibility of alternative, perhaps better, fits to a calibration curve. We feel that it is premature to set an arbitrary “smoking gun” value such as 100 years (c. 1.2% old carbon contamination) when such contamination can conceivably affect fits at any level above 0%.
***Contamination from any source can conceivably vary from 0 through “very slight” (0.25 – 1%?) to considerable (>1%) which can occur at different times during the life of the tree leading up to the time of the eruption which is inferred to have killed the tree. Changing levels of contamination could be enough to distort a wiggle match series and create a false match and this is a possibility that must be faced. At present, the standard approach is to seek out and remove ages that compromise the fit at the assumed (no contamination) section of the contamination curve.
(What particularly stands out in the figure to me are the younger outliers which are not explained by old carbon contamination. I like that these are addressed in the text in terms of uptake of younger carbon, but here again but a wide range of other possibilities (like the geological interpretation of the context for the sample) are not really explored (though it seems like the later spatial information might also be very relevant here)... But this is a side issue... something to think about if reframing and refocusing this text for resubmission)
***We agree that the “young” eruption dates are problematic and will expand the discussion on that issue. We note that these are not “outliers” in the conventional sense: each whole age series provides a match with the then current calibration curve that was accepted before later wiggle matches were attempted. The “young dates” resulted in the former name “Millennial Eruption”.
- All other comments aside, I would not feel that this approach could be applied with confidence in other studies especially given that the authors note that it can only be applied where a consistent level of contamination can be assumed over the life of the tree.
***We do not suggest that the method can be applied only where a consistent level of contamination is assumed. We propose the method recognising that by applying a constant contamination is better than applying no contamination at all (a fact which is not recognised in conventional wiggle match dating) as an essential first step in assessing a wiggle match series to see if the assumption of no contamination, consistent or otherwise, is valid. Where there are better matches using the constant contamination assumption then surely the possibility of contamination should be considered. We say that, and will make it clearer in a revised MS that we simplify the analysis by assuming a constant rate. This is analogous to the present situation where no contamination is assumed and only a single fit sought. The number of wiggle match ages is important here, as series of only 4 or 5 ages can conceivably fit several parts of a curve. By introducing the method with this assumption we hope more studies follow that explore non constant contamination models and suspect these will improve dating accuracy further..
Volcanic eruptions and the potential for volcanic gas release and other factors which might impact storage of carbon in sequential tree-ring series overtime are highly variable and very site/eruption specific. Also factors such as growth season of the trees, the period of time covered by the tree-ring sequence involved and the fact that old carbon may be released in some years and not in others all combines to complicate the application of a blanket correction factor for a given study. This sort of blanket correction could equally apply to correct a lab-off set and this would be nothing new.
***We agree, but this is seldom if ever done at the moment, and suggest our methodology would be equally useful for identifying and adjusting for interlab off set.. The comment serves to emphasise what our MS proposes, that wiggle match dating of eruptions (or any other event where systematic contamination/offset might occur) is not simple and investigation of potential methods of refining the analyses are necessary, rather than a statistically Procrustean forced fit by ad hoc exclusion of ages that compromise a no-contamination fit. It is certainly an uncomfortable prospect but one which cannot be ignored.
A few other more specific comments to help with revisions
P1: Lines 1-2: Suggest simplification from: "increasing accuracy as the technologies developed, more refined calibration curves promulgated (Hogg et al., 2020; Reimer et al., 2020), and sophistication of statistical analyses (Crema and Bevan, 2020) to: "increasing accuracy with the use of more refined calibration curves (Hogg et al., 2020; Reimer et al., 2020), and sophisticated statistical analyses (Crema and Bevan, 2020)"
***We will consider the suggested simplification.
P2: Lines 5-6: Suggest change from " as the most likely to give a calendar date " to " as the most likely to give a close to calendar date " (calendar date implies actual calendar date which is only possible with a 14C excursion such as the 774/5 CE event, not through conventional radiocarbon wigglematching, so either make this more subtle change or expand a bit further on making the distinction between the two different ways of using radiocarbon to date samples from eruption contexts, this may be helpful to the reader at this stage as the study rather depends on people appreciating the distinction.)
***We take the reviewers’ point and will expand the text to make the distinction clearer.
P4: Line 27: "cosmogenically-tuned date" - change to "cosmogenically secured" important for readers to understand that this 'radiocarbon date' based on the 774/5 CE event IS a calendar date not subject to variations like wigglematches always will be with new iterations of the curve. Also on line 27 you introduce 3 possible options for why previous radiocarbon dates might be older than the date for the eruption given by Oppenheimer / Hakozaki and then go into a nicely comprehensive list of possible volcanic causes... but you do not introduce the significant range of other possibilities to explain this - for example, were Sun et al's dates based on a sequence of older to younger samples from a single tree? If they were, how old was the tree in total and were the samples blocks of years? Half the dates being older in a consecutive sequence from a tree would be predicted, especially if that tree was longer lived. Across how many years is the 14C of a single sample averaged to get the values for the wigglematched blocks and how much does the calibration curve vary (or how much could it vary around the time of the event (i.e. curve plateaus or slopes) - these are all really important to address and provide a thorough framing for the 'contamination' discussion..
***We will explore the implications as to other possible causes (a point raised also by Reviewer 2) and will alter the text accordingly.
P6: Figure 2, line 3, 774 not 74
***We will correct the typo.
Citation: https://doi.org/10.5194/gchron-2021-13-AC1
-
AC1: 'Reply on RC1', Richard Holdaway, 29 Sep 2021
-
RC2: 'Comment on gchron-2021-13', Anonymous Referee #2, 08 Sep 2021
This is a really interesting article that also raises an interesting point that there is at least the potential that trees being dated by Dendro/14C wiggle matching may be influenced by gological dead 14C. As such the paper should be published. It is an important but fundamentally overlooked point. I do, however, have a number of cocernes that need in my mind to be addressed.
My first issue is one that is also raised by the first refereee, in that there are many sources of contamination that can cuase problems with 14C ages, and my not be removed by standard techniques. There is clearly a problem that is unlikely to be related to volcanic CO2 in the ages that are too young, as opposed to being too old. This needs to be discussed in more detail. In fact I think the paper needs to more be an examination of the problems that can be encountered with 14C ages and placed in the context of wiggle match dating of eruptions. In this unusual case we know the age very specifically and that allows a discussion of these problems. Some of the peroblems with the dates relates to contamination by older 14C, which may be from volcanic degassing, but could be from a range of other issues. Some dates have been contaminated by young 14C and this is unlikely to be a volcanic origin. One way of looking for potential contamination sources is through delta 13C measurements. In cases where 14C and 13C of volcanic emissions have been considered on modern trees there is evidence of a shift in both isotopes. The fact that it is pretty difficult to find 13C measurements for most of the cases used in this study makes the point that there are a whole host of quality assurance issues in raduiocarbon dating of volcanic eruptions that need to be considered. The known age of the B-TM makes this a great paper around which to have that wider discussion. This is moreso becuase I do not see how, unless we know the exact age of the eruption, we can use this approach, to finding old carbon. Also even is we can find the percentage of old carbon in a sample, withouot additional measurements such as 13C we can't start to suggest the source of the old carbon with much reliability.
I also wasn't sure why the outlier function in OxCal wasn't used instead as this has largely superceedeed the use of AI's in OxCal outlier detection studies. The AI is simply the area of overlap between posterior and prior. I would like to see at the very least the analyses done again with appropriate outlier detection.
Citation: https://doi.org/10.5194/gchron-2021-13-RC2 -
AC2: 'Reply on RC2', Richard Holdaway, 29 Sep 2021
Reviewer 2
This is a really interesting article that also raises an interesting point that there is at least the potential that trees being dated by Dendro/14C wiggle matching may be influenced by gological dead 14C. As such the paper should be published. It is an important but fundamentally overlooked point. I do, however, have a number of cocernes that need in my mind to be addressed.
***We appreciate the reviewer’s view that our paper deals with “an important but fundamentally overlooked point” and are pleased with their recommendation. We respond to their concerns in detail below.
My first issue is one that is also raised by the first refereee, in that there are many sources of contamination that can cuase problems with 14C ages, and my not be removed by standard techniques. There is clearly a problem that is unlikely to be related to volcanic CO2 in the ages that are too young, as opposed to being too old. This needs to be discussed in more detail. In fact I think the paper needs to more be an examination of the problems that can be encountered with 14C ages and placed in the context of wiggle match dating of eruptions.
***As in our response to this point as raised by Reviewer 1, we will expand our discussion of other potential sources of contamination while confirming our position that, for wiggle match series measured on trees within the influence of a volcanic edifice, geologic/magmatic carbon is an unappreciated but potentially dominant source. And one which can provide enough contaminating carbon to actually change the position of best fit on a calibration curve.
In this unusual case we know the age very specifically and that allows a discussion of these problems. Some of the peroblems with the dates relates to contamination by older 14C, which may be from volcanic degassing, but could be from a range of other issues. Some dates have been contaminated by young 14C and this is unlikely to be a volcanic origin. One way of looking for potential contamination sources is through delta 13C measurements. In cases where 14C and 13C of volcanic emissions have been considered on modern trees there is evidence of a shift in both isotopes. The fact that it is pretty difficult to find 13C measurements for most of the cases used in this study makes the point that there are a whole host of quality assurance issues in raduiocarbon dating of volcanic eruptions that need to be considered.
***We agree that δ13C is a potentially useful marker of the source of any contaminating carbon. We have explored, and continue to explore, the possibility of δ13C evidence for the source(s) of contamination. As the reviewer notes, such data are rarely reported and we very much agree that there are issues of quality assurance in radiocarbon dating of eruptions.
The known age of the B-TM makes this a great paper around which to have that wider discussion. This is moreso becuase I do not see how, unless we know the exact age of the eruption, we can use this approach, to finding old carbon. Also even is we can find the percentage of old carbon in a sample, withouot additional measurements such as 13C we can't start to suggest the source of the old carbon with much reliability.
***We agree and certainly hope that our paper could be a focus for a wider, constructive, discussion on the issues which undoubtedly surround the wiggle match process. We propose our method as a way of identifying the potential presence of old carbon in a wiggle match tree, by exploring the possibility of the existence ofother, conceivably better, matches if zero contamination is not assumed.
I also wasn't sure why the outlier function in OxCal wasn't used instead as this has largely superceedeed the use of AI's in OxCal outlier detection studies. The AI is simply the area of overlap between posterior and prior. I would like to see at the very least the analyses done again with appropriate outlier detection.
***We investigated but did not pursue the use of the “outlier” function in OxCal. We think that this is basically another way of isolating ages that don’t fit with an assumed position of fit on the curve. A different way of applying the Procrustean approach to a wiggle match fit. As such, we do not believe that its application would be appropriate here. We used AIs because they are recognised numerical method of assessing overall fit of a wiggle match series to a calibration curve. They allow fits at different sections of the curve to be compared objectively.
Citation: https://doi.org/10.5194/gchron-2021-13-AC2
-
AC2: 'Reply on RC2', Richard Holdaway, 29 Sep 2021
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 804 | 423 | 49 | 1,276 | 132 | 31 | 31 |
- HTML: 804
- PDF: 423
- XML: 49
- Total: 1,276
- Supplement: 132
- BibTeX: 31
- EndNote: 31
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
