THE EARLY YEARS
Not long after James Arnold and Willard Libby published the “Curve of Knowns” as preliminary proof that radiocarbon (14C) dating worked at least for the past 5000 years (Arnold and Libby Reference Arnold and Libby1949), researchers began to evaluate the validity of two of the fundamental assumptions of the method. Hessel de Vries, who pioneered radiocarbon dating at the University of Groningen, was the first to convincingly demonstrate, using tree rings of known age, that the radiocarbon concentration in the atmosphere was not constant but varied by up to 1% over the last 400 years (de Vries Reference de Vries1958). This study was expanded on by an intercomparison exercise undertaken by the laboratories of Cambridge, Copenhagen, and Heidelberg on 25 sequoia tree rings sampled at 50-year intervals from AD 659 to AD 1859 (Willis et al. Reference Willis, Tauber and Munnich1960). They reported a variation of 1.5% and observed that there appeared to be a periodicity of 150–200 years superimposed on a longer period (Figure 1). That same year the University of Pennsylvania laboratory published results for a number of dendrochronologically dated tree rings (mostly Douglas fir) and some well-constrained archaeological samples relative to the average value of 39 oak tree rings ranging from AD 1747 to 1849 (Ralph and Stuckenrath Reference Ralph and Stuckenrath1960). They found three time periods (AD 1687, 1297 BC, and 1925 BC) with significant deviations from the oak average. Stuiver (Reference Stuiver1961) noted a correlation between solar activity and the 14C of known age tree rings which was strengthened by further measurements (Stuiver Reference Stuiver1965; Suess Reference Suess1965). Stuiver and Suess (Reference Stuiver and Suess1966) analyzed known age tree rings from European oak, North American Douglas fir and sequoia at La Jolla and Yale relative to the oxalic acid standard provided by the U.S. National Bureau of Standards (NBS). They published tables of radiocarbon corrections and what amounts to the first radiocarbon calibration curve extending from AD 1000 to 1800 (Figure 2).
Precise 14C measurements on tree rings from Europe, North and South America, South Africa, and India were made at the University of Groningen to investigate the possibility of geographic differences in atmospheric 14C (Lerman et al. Reference Lerman, Mook, Vogel and Olsson1970). They found that the Northern Hemisphere values were similar for all locations, which supported the radiocarbon dating assumption of a well-mixed atmosphere. However, there measurements showed that Southern Hemisphere samples from around AD 1835 were 4.5 ± 1‰ depleted compared to Northern Hemisphere samples. They also provided the first Southern Hemisphere calibration curve from measurements of Argentine trees from AD 1400–1950.
The construction of a bristlecone pine chronology at the University of Arizona Laboratory of Tree Ring Research (Ferguson Reference Ferguson1968) made it possible to extend radiocarbon calibration. Suess published calibration corrections from 4100 BC to 1500 BC based on bristlecone pine measured at La Jolla (Suess Reference Suess1967). The calibration was extended with further measurements of bristlecone pine from 5400 BC to the present (Suess Reference Suess1970) and was widely used by archaeologists (Ottaway and Ottaway Reference Ottaway and Ottaway1972). As noted in comments by Renfrew, calibration had a tremendous impact on the interpretation of the dates of prehistoric Europe (Bermingham and Renfrew Reference Bermingham and Renfrew1972).
Suess’s calibration curve was not overwhelming accepted as a number of researchers did not believe that the fine structure or “wiggles” were real (Burleigh et al. Reference Burleigh, Longworth and Wainwright1972; Switsur Reference Switsur1973). The fact that Suess said he had drawn the curves with “cosmic schwung” to demonstrate the “most probable character” of the curve did not instil confidence (Adams Reference Adams1973; Olsson Reference Olsson2009). Switsur (Reference Switsur1973) compared the calibration tables presented by the University of Arizona (Damon et al. Reference Damon, Ferguson, Long and Wallick1974) and the University of Pennsylvania (Ralph et al. Reference Ralph, Michael and Han1974) at the 8th International Conference on Radiocarbon Dating at Lower Hutt City, New Zealand, and found that the difference in corrected ages was fairly small. He therefore provided calibration corrections from the averages of the Arizona and Pennsylvania measurements. These corrections had what was felt to be an advantage over the Suess calibration corrections as they produced only one calibrated age for each radiocarbon measurement. Suess argued “nature is not necessarily as simple as we would wish and an attempt to draw a smooth curve amongst the measured points is probably the wrong approach” (Switsur Reference Switsur1973).
Renfrew and Clark (Reference Renfrew and Clark1974) reviewed all the curves that had been published to date, as well as a formula for calibration (Wendland and Donley Reference Wendland and Donley1971) and found these were all statistically unsatisfactory. Clark (Reference Clark1975) then examined all the data and excluded duplicates as well as some data that did not use the NBS oxalic acid standard. He found that the variability between replicates of the same tree rings in different laboratories was in excess of the reported errors which generally only included counting statistics. Clark used a spline function to construct a curve from the data with the standard error determined from the replicates and a smoothing parameter determined by a cross-validation technique. The resulting curve was smoother than that of Suess (Reference Suess1970) but not as smooth as Switsur’s (Reference Switsur1973). Indeed, Clark concluded that the “kinks and wobbles” in Suess’s calibration curve are “not justified statistically.” Suess responded that there was no statistical justification necessary for kinks, whereas Clark responded that his method had been tested and would detect “kinks” if they were there and not swamped by measurement errors (Suess and Clark Reference Suess and Clark1976). No international agreement on a calibration curve could be reached at the Ninth International Conference Los Angeles and La Jolla in 1976. It is no wonder that archaeologists began to doubt the validity of radiocarbon dating and calibration (Pilcher and Baillie Reference Pilcher and Baillie1978).
Meanwhile German and Irish oak dendrochronologies, though mostly still floating (i.e., not dendrochronologically linked to known calendar ages), were under development (Becker Reference Becker1979; Pilcher and Baillie Reference Pilcher and Baillie1978). Suess presented measurements of floating German oak and California bristlecone pine to show that the 14C “wiggles” in them were synchronous and of the same amplitude within uncertainty (Suess Reference Suess1979). High-precision gas proportional counters were developed in Groningen with an overall precision of 1.3‰ for modern carbon (Tans and Mook Reference Tans and Mook1978) and in Seattle with a precision of 1.5–2‰ for a Pacific coast tree-ring series dating back to 1200 AD (Stuiver et al. Reference Stuiver, Robinson and Yang1979). A 500-yr-long section of the Southern German Neolithic floating chronology oak was measured in Groningen (de Jong et al. Reference de Jong, Mook and Becker1979) which confirmed the existence of the Suess-type wiggles. Pearson (Reference Pearson1979) developed high-precision liquid scintillation counters and made a thorough study of the source of errors in the data. With careful controls he was able to produce 14C data on samples of Irish oak with an overall uncertainty of 2.5‰ (ca. 20 14C yr).
In 1979, a meeting was held at Tucson, Arizona, with the primary objective being to construct a radiocarbon calibration correction table using only 14C analyses funded by the U.S. National Science Foundation (Damon et al. Reference Damon, Lerman, Long, Bannister, Klein and Linick1980). A method for producing the radiocarbon calibrations was agreed based on multistage linear regression previously presented by Ralph and Klein (Reference Ralph and Klein1978). The workshop report also stated that “the most important accomplishment of the Workshop was the participants’ universal agreement to accept a single calibration curve based on the composite data set.” The participants also recommended establishing an International Calibration Committee to extend the workshop to include all available data worldwide. Statistical analysis results from the workshop were reported at the 10th International Radiocarbon Conference in Bern (Klein et al. Reference Klein, Lerman, Damon and Timothy1980). The calibration curve and tables constructed as agreed were published two years later (Klein et al. Reference Klein, Lerman, Damon and Ralph1982).
A number of new radiocarbon measurements were also presented at the 10th International Radiocarbon conference. Pearson (Reference Pearson1980) reported dates on 20-ring blocks of Irish and Scottish oak trees from AD 1510 to AD 1840 using liquid scintillation counters which were in agreement with those from Stuiver (Reference Stuiver1978) on North American Douglas fir measured using gas proportional counters. A visual match of an older floating Irish tree-ring sequence with that of the Suess bristlecone pine provided an approximate age range of 1810–3535 BC. The data exhibited a saw tooth pattern with a period of 110 to 180 years, similar to that seen in the German oak (de Jong et al. Reference de Jong, Mook and Becker1979) with which there was a short overlap. Pearson concluded, however, that the variability in the bristlecone data (Suess Reference Suess1979) was too high to identify wiggles. Bruns et al. (Reference Bruns, Munnich and Becker1980) measured 1–3 ring samples of Southern German Oak from AD 200-800 with gas counters at the University of Heidelberg. Using Fourier analysis, they determined a periodicity of 150 to 180 years similar to that seen by de Jong et al. (Reference de Jong, Mook and Becker1979) and de Jong and Mook (Reference de Jong and Mook1980) in tree rings from several millennia earlier. Likewise, after removing a sine curve with a period of 12,100 years and applying a smoothing spline, Fourier analysis showed a significant 200-yr period in the bristlecone pine data (Suess Reference Suess1980).
The theoretical framework for the Suess “wiggles” was set out by Stuiver and Quay (Reference Stuiver and Quay1980). They used a carbon reservoir model to calculate the 14C production rate from precise 14C measurements of Pacific coast tree rings. The calculated production rates were then compared to those derived from neutron flux measurements and geomagnetic indices and to historical records of sunspot numbers and auroral indices. This work showed that the natural atmospheric 14C changes were related to Earth’s geomagnetic field variation and changes in solar wind shielding properties. In later work, ocean circulation and climate changes were also shown to influence atmospheric 14C (Stuiver et al. Reference Stuiver, Braziunas, Becker and Kromer1991).
The Klein et al. (Reference Klein, Lerman, Damon and Ralph1982) calibration curve was based on dendrochronologically dated trees, primarily bristlecone pine and sequoia of 20 or fewer rings measured at the Universities of Groningen, California at La Jolla, Pennsylvania, and Yale. No other selection criteria were applied and no adjustment was made for systematic differences between laboratories (Klein et al. Reference Klein, Lerman, Damon and Timothy1980), however, the uncertainty from the interlaboratory standardization was included in the calibration uncertainty. The data were weighted by their measurement uncertainty and an additional term to account for increasing uncertainty with sample age and scaled logarithmically. A polynomial regression was done to derive the long period of the curve on the order of a few thousand years. Next a piecewise regression was done on the residuals around the polynomial in intervals of 500 years (“shingles”) overlapping by 250 years and Fourier analysis was done to determine the short period variability. A composite calibration function with uncertainty was determined from these. A section of the curve from AD 600–1950 is shown in Figure 3 with 90% confidence limits for radiocarbon errors of 20, 100, 200, and 300 14C yr.
Stuiver produced a high-precision decadal calibration curve for AD 1–1950 based on Douglas Fir from the Pacific Northwest and Sequoia from California (Stuiver Reference Stuiver1982). Measurement errors were ±12 14C yr, which were increased to ±18 14C yr with an error multiplier derived from comparison of a set of 30 pairs of contemporaneous samples from different trees. The difference between the Seattle data and that from Belfast (Pearson Reference Pearson1980) for this period was 4.4 14C years with mean age difference of 2 ± 3 14C yr between the complete data sets of 53 sample pairs (AD 955–1840). Offsets varied between 27 for German oak measured and La Jolla and Douglas fir in Seattle, between 37 and 55 14C yr for La Jolla bristlecone pine versus Seattle sequoia and Seattle Douglas fir, respectively. Preliminary measurements from Heidelberg on German oak were offset 58 14C yr from Seattle sequoia.
At the business meeting of the 11th International Radiocarbon Conference in Seattle, 1982, the International Calibration Committee set up by Damon et al. (Reference Damon, Lerman, Long, Bannister, Klein and Linick1980) “was discharged with thanks for their useful work and a smaller committee with W. G. Mook as chairman was re-constituted for the 1982–1985 interval” (Stuiver Reference Stuiver1983).
THE STUIVER YEARS
The 1986 Calibration Curves
The 12th International Radiocarbon Conference, Trondheim, 1985, saw a proliferation of calibration datasets and curves. The bi-decadal curves from 2500 BC to the present (Pearson and Stuiver Reference Pearson and Stuiver1986; Stuiver and Pearson Reference Stuiver and Pearson1986) were officially recommended at the business meeting (Mook Reference Mook1986), the first time the radiocarbon community had come to an agreement on radiocarbon calibration. These curves, which were based on measurements of Irish oak (Pilcher et al. Reference Pilcher, Baillie, Schmidt and Becker1984) and trees from the western U.S., were published in a special calibration issue that also included a decadal curve for the same interval based on the Seattle measurements only (Stuiver and Becker Reference Stuiver and Becker1986). Additional measurements on German oak, both dendrochronologically dated and floating, and bristlecone pine were published in the calibration issue (Stuiver Reference Stuiver1986 and references therein). Stuiver (Reference Stuiver1986) also suggested the possibility of extending the calibration curve with varved sediment data from Lake of the Clouds.
These decadal and bidecadal data, with the exception of the 14C data from varved sediment and the older floating German oak chronology, were combined with a simple weighted average of the data within each bidecade, rather than using more complicated statistical processes, to create the curve atm20 for calibration of terrestrial samples from 7210 cal BC–AD 1950 using the computer program CALIB (Stuiver and Reimer Reference Stuiver and Reimer1986). An ocean-atmosphere box model developed by Oeschger et al. (Reference Oeschger, Siegenthaler, Schotterer and Gugelmann1975) was used to construct a “global” marine radiocarbon calibration curve using the atmospheric curve as input age (Stuiver et al. Reference Stuiver, Pearson and Braziunas1986). The concept of ΔR, the difference between the “global” ocean radiocarbon and that of a specific region was introduced and a table of ΔR values given based on known age, pre-bomb samples for use with the marine curve.
A comparison of different calibration methods and programs, using the 1986 datasets, was undertaken by Aitchinson et al. (1989). The theoretical basis for combining archaeological information and radiocarbon dates through a Bayesian calibration model was developed by Buck et al. (Reference Buck, Kenworthy, Litton and Smith1991) and the first versions of the software packages OxCal (Bronk Ramsey Reference Bronk Ramsey1995) and BCal (Buck et al. Reference Buck, Christen and James1999) were made available in due course. Both BCal and OxCal utilize the ratified IntCal curves as do age-depth modeling programs such as BChron (Parnell et al. Reference Parnell, Haslett, Allen, Buck and Huntley2008) and BACON (Blaauw and Christen Reference Blaauw and Christen2011).
The 1993 Update
The calibration curves were corrected and extended and a probability method of calibration (Stuiver and Reimer Reference Stuiver and Reimer1989) was added to the CALIB program (Stuiver and Reimer Reference Stuiver and Reimer1993). A correction to earlier Seattle data had to be made due to a change in lab procedure that introduced radon into the samples which produced alpha particles during decay (Stuiver and Becker Reference Stuiver and Becker1993). The radon correction was based on the difference between the first and last day of counting for hundreds of samples. The Belfast data was also adjusted to account for “the variation of efficiency with time” in standard counts (Pearson and Qua Reference Pearson and Qua1993). Additional data of German oak measured by Heidelberg and floating German pine with a tentative match to the oak (Kromer and Becker Reference Kromer and Becker1993) were also included in the bidecadal curve. Adjustments for laboratory offsets were also made. For the first time, U-Th and 14C dated coral (Bard et al. Reference Bard, Arnold, Fairbanks and Hamelin1993) were included with a 400-yr reservoir offset to extend the curve with a smoothing spline to nearly 22,000 cal BP. A ratification vote for the bidecadal IntCal93 and the corresponding Marine93 curves was not sought and the curves were not generally accepted in Europe (F. G. McCormac, pers. comm.)
IntCal98 and Marine98
In 1996 a workshop was held at the Wissenschaftsforum, Heidelberg, Germany, to address aspects of high-precision radiocarbon calibration (Kromer et al. Reference Kromer, Ambers, Baillie, Damon, Hesshaimer, Hofmann, Jöris, Levin, Manning, Mccormac, van der Plicht, Spurk, Stuiver and Weninger1996). With regards to calibration data, two main points were made. An intercomparison of the Hohenheim and the Gottingen German oak chronology (Leuschner and Delorme Reference Leuschner and Delorme1988; Becker Reference Becker1993) revealed that 41 years were missing from the published Hohenheim chronology (Spurk et al. Reference Spurk, Friedrich, Hofmann, Remmele, Frenzel, Leuschner and Kromer1998). The previous link between the German oak and Preboreal Pine Chronology (Becker Reference Becker1993; Kromer and Becker Reference Kromer and Becker1993) was also now discarded so the pine chronology was left floating. Gerry McCormac presented evidence to support the original Irish oak 14C data (Pearson et al. Reference Pearson, Pilcher, Baillie, Corbett and Qua1986) as opposed to the revised dates (Pearson et al. Reference Pearson, Becker and Qua1993) used in the 1993 calibration curves. Further potential calibration datasets were presented at the 16th International Radiocarbon Conference in Groningen in 1997. Among these were U-Th dated coral datasets from Barbados, Mururoa and Tahiti (Bard et al. Reference Bard, Arnold, Hamelin, Tisnerat-Laborde and Cabioch1998) and from Vanuatu (Burr et al. Reference Burr, Beck, Taylor, Recy, Edwards, Cabioch, Correge, Donahue and O’Malley1998) and foraminifera from varved Cariaco Basin sediments (Hughen et al. Reference Hughen, Overpeck, Lehman, Kashgarian and Southon1998). The agreement of the high resolution Cariaco Basin data (Hughen et al. Reference Hughen, Overpeck, Lehman, Kashgarian and Southon1998) with the Barbados coral data (Bard et al. Reference Bard, Arnold, Hamelin, Tisnerat-Laborde and Cabioch1998) convinced Minze Stuiver to include varve data in the calibration curve for the first time. At the conference a few scientists expressed their discontent to me about having their data included in the calibration curve without direct involvement. At my suggestion they went to talk to Minze Stuiver about it and he agreed that they would become co-authors on the resulting publication. The international part of the curve name “IntCal” became a reality.
A few other corrections were made to previously published data. Comparison of tree-ring samples remeasured in the Seattle lab showed that the radon correction made in 1993 needed to be reduced by 50% (Stuiver et al. Reference Stuiver, Reimer and Braziunas1998b). A weak section was found on re-examination of the older part of the German oak chronology and additional trees were found to bridge this section resulting in a shift of 54 yr to the older trees (Spurk et al. Reference Spurk, Friedrich, Hofmann, Remmele, Frenzel, Leuschner and Kromer1998). In addition, a 14C match of the floating German pine chronology to the German oak provided a tentative link and extended the chronology to 11,857 cal BP with an error estimate for the link of 20 cal years (Kromer and Spurk Reference Kromer and Spurk1998). Other new tree-ring data included in the IntCal98 curve were German oak and pine, Sitka spruce and Douglas fir (Stuiver et al. Reference Stuiver, Reimer, Bard, Beck, Burr, Hughen, Kromer, McCormac, van der Plicht and Spurk1998a), Irish oak (McCormac et al. Reference McCormac, Hogg, Higham, Baillie, Palmer, Xiong, Pilcher, Brown and Hoper1998a, Reference McCormac, Hogg, Higham, Lynch-Stieglitz, Broecker, Baillie, Palmer, Xiong, Pilcher, Brown and Hoper1998b) and previously published data (Kromer et al. Reference Kromer, Rhein, Bruns, Schochfischer, Munnich, Stuiver and Becker1986; Pearson et al. Reference Pearson, Becker and Qua1993; Vogel and van der Plicht Reference Vogel and van der Plicht1993).
Previously, 14C ages on dendrodated wood covering only a few years (sub-decadal) were disregarded, but a different approach was used for IntCal98. The “decadal” values were obtained by averaging all full decadal and sub-decadal results. Samples measured on 20-yr blocks of wood were also included by assigning this age to two decades with the standard deviation multiplied by 1.4. Laboratory error multipliers (k values) and offsets were estimated for each tree-ring dataset based on comparison of results from samples of identical calendar age.
Marine reservoir offsets (R) were estimated for the coral data from Bard et al. (Reference Bard, Arnold, Hamelin, Tisnerat-Laborde and Cabioch1998), Burr et al. (Reference Burr, Beck, Taylor, Recy, Edwards, Cabioch, Correge, Donahue and O’Malley1998), and Edwards et al. (Reference Edwards, Beck, Burr, Donahue, Chappell, Bloom, Druffel and Taylor1993) by comparing the marine 14C ages to overlapping tree-ring data. The weighted mean R value was 414 ± 31 14C yr (n =12) for 10,000–8000 cal BP and 509 ± 25 14C yr (n = 21) for 12,000–10,000 cal BP. While the possibility of reservoir age changes over time was acknowledged it was decided to use a value of 500 14C yr for the Late Glacial tropical surface ocean. The coral data were found to have higher 14C variability than the tree rings. To compensate for this the standard deviation of coral 14C ages was based on a 2σ error in the coral 14C measurement and an error multiplier of 1.3 was applied. Foraminifera 14C ages from the Cariaco Basin decreased by 500 and 400 14C yr for, respectively, 12–10 and 10–8 ka cal BP were matched to the IntCal98 tree-ring data to update the floating varve chronology. The observed scatter in the 14C ages was similar to that of the tree rings with an error multiplier of k = 1.3. The IntCal98 curve from 12,000 to 24,000 was generated from the reservoir corrected marine data using a spline with minimum smoothing (Stuiver et al Reference Stuiver, Reimer, Bard, Beck, Burr, Hughen, Kromer, McCormac, van der Plicht and Spurk1998a).
Marine98 curve, which extended from 0 to 12,000 cal BP, was generated from the IntCal98 curve using an atmosphere-ocean box diffusion model (Stuiver et al. Reference Stuiver, Reimer and Braziunas1998b) with model parameters as discussed in Stuiver and Braziunas (Reference Stuiver and Braziunas1993).
THE INTCAL WORKING GROUP
2004 Calibration Curves
When Minze Stuiver retired, I continued the calibration effort with encouragement from Gerry McCormac. Together we received a Leverhulme Trust networking grant to form the IntCal Working Group (IWG) in 2002. The grant funded meetings in Belfast and Woods Hole and a research assistant. Previous and potential data providers, as well as an archaeologist and a statistician, were invited to join the group. At the first meeting in Belfast in 2002 the IWG set out criteria for inclusion of data into the planned IntCal04 calibration curve (Reimer et al. Reference Reimer, Hughen, Guilderson, McCormac, Baillie, Bard, Barratt, Beck, Buck, Damon, Friedrich, Kromer, Ramsey, Reimer, Remmele, Southon, Stuiver and van der Plicht2002). It was agreed that a new statistical method was needed, in part to get rid of what Ron Reimer christened the “pig-in-the-python” around 15,000 cal BP (Figure 4). At the second meeting in Woods Hole in 2003, a random walk model was chosen to account for the uncertainty in the calendar ages of samples as well as the uncertainty in their 14C ages to calculate the curve (Buck and Blackwell Reference Buck and Blackwell2004). Constant regional offsets were used to incorporate the marine data. The IntCal04 and Marine04 curves were extended to 26,000 cal BP (Hughen et al. Reference Hughen, Baillie, Bard, Beck, Bertrand, Blackwell, Buck, Burr, Cutler, Damon, Edwards, Fairbanks, Friedrich, Guilderson, Kromer, Mccormac, Manning, Ramsey, Reimer, Reimer, Remmele, Southon, Stuiver, Talamo, Taylor, van der Plicht and Weyhenmeyer2004; Reimer et al. Reference Reimer, Baillie, Bard, Bayliss, Beck, Bertrand, Blackwell, Buck, Burr, Cutler, Damon, Edwards, Fairbanks, Friedrich, Guilderson, Hogg, Hughen, Kromer, McCormac, Manning, Ramsey, Reimer, Remmele, Southon, Stuiver, Talamo, Taylor, van der Plicht and Weyhenmeyer2004) but no recommendation was given beyond that because of the disparities between datasets. However, a predictive curve called “NotCal04” (to signify that it was not for calibration purposes) was plotted based on a random effects model (van der Plicht et al. Reference van der Plicht, Beck, Bard, Baillie, Blackwell, Buck, Friedrich, Guilderson, Hughen, Kromer, Mccormac, Ramsey, Reimer, Reimer, Remmele, Richards, Southon, Stuiver and Weyhenmeyer2004). This paper also demonstrated the erroneous conclusions that could result from attempts to calibrate radiocarbon ages in this time frame with alternative datasets, such as a smooth curve proposed by van Andel (Reference van Andel1998) based on geomagnetic records, and caution was urged. Nonetheless the NotCal04 curve was used for calibration of Palaeolithic radiocarbon dates such as those from Chauvet cave (e.g., Mellars Reference Mellars2006). van Andel (Reference van Andel2005) argued that alternative datasets and calibration methods such as CalPal (Jöris and Weninger Reference Jöris and Weninger1998) should not be discouraged due to the need to calibrate radiocarbon ages beyond the 26,000 cal BP. He also disputed the importance of the IntCal consensus calibration to the community. Fortunately, this issue was resolved by further updates to the IntCal curves extending to close to the limits of radiocarbon dating and increasing community involvement in the process. In addition to IntCal04 and Marine04, a preliminary calibration curve for the Southern Hemisphere atmosphere, SHCal02 (McCormac et al. Reference McCormac, Reimer, Hogg, Higham, Baillie, Palmer and Stuiver2002) was extended to 11,000 cal BP using the tree-ring portion of the IntCal04 curve with a modeled offset between the Northern and Southern Hemisphere (McCormac et al. Reference McCormac, Hogg, Blackwell, Buck, Higham and Reimer2004).
2009 Updates
In 2009 the IntCal09 and Marine09 calibration curves (Reimer et al. Reference Reimer, Baillie, Bard, Bayliss, Beck, Blackwell, Ramsey, Buck, Burr, Edwards, Friedrich, Grootes, Guilderson, Hajdas, Heaton, Hogg, Hughen, Kaiser, Kromer, Mccormac, Manning, Reimer, Richards, Southon, Talamo, Turney, van der Plicht and Weyhenmeyer2009) were extended to 50,000 cal BP using new data from U-Th dated coral (Fairbanks et al. Reference Fairbanks, Mortlock, Chiu, Cao, Kaplan, Guilderson, Fairbanks, Bloom, Grootes and Nadeau2005) and from marine sediment with timescales transferred through climatic correlation to the U-Th dated Hulu Cave speleothem (Wang et al. Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001). These “tie-pointed” records included the non-varved portion of Cariaco Basin (Hughen et al. Reference Hughen, Southon, Lehman, Bertrand and Turnbull2006) and the Iberian margin (Bard et al. Reference Bard, Ménot-Combes and Rostek2004). The Marine09 curve was taken from the Marine04 curve for 0–12,500 cal BP whereas the older portion was derived from the IntCal09 curve with a constant 405 14C yr offset. A Bayesian evaluation of the interhemispheric offset was undertaken but SHCal04 was not updated at this time (Hogg et al. Reference Hogg, Ramsey, Turney and Palmer2009). From a statistical perspective, the random walk model for curve construction was considerably developed for the 2009 update. The IntCal09 extension was the first time a fully Bayesian approach, achieved via Markov Chain Monte Carlo (MCMC), had been used to estimate the calibration curve (Blackwell and Buck Reference Blackwell and Buck2008; Heaton et al. Reference Heaton, Blackwell and Buck2009).
2013 Calibration Curves
For the 2013 curves, three scientists were elected by the community to the IntCal Oversight Committee to provide independent expertise and assessment. New tree-ring measurements were included in the IntCal13 curve most importantly from the floating Late Glacial Pine chronology (Schaub et al. Reference Schaub, Buntgen, Kaiser, Kromer, Talamo, Andersen and Rasmussen2008) anchored by a 14C wiggle-match as described by Hua et al. (Reference Hua, Barbetti, Fink, Kaiser, Friedrich, Kromer, Levchenko, Zoppi, Smith and Bertuch2009) which extended to 14,200 cal BP. The density of data beyond the tree rings was greatly increased compared to IntCal09 due to the incorporation of 14C measurements from speleothems from the Bahamas (Beck et al. Reference Beck, Richards, Edwards, Silverman, Smart, Donahue, Hererra-Osterheld, Burr, Calsoyas, Jull and Biddulph2001; Hoffmann et al. Reference Hoffmann, Beck, Richards, Smart, Singarayer, Ketchmark and Hawkesworth2010) and Hulu Cave, China (Southon et al. Reference Southon, Noronha, Cheng, Edwards and Wang2012) and from terrestrial plant macrofossils from Lake Suigetsu, Japan (Kitagawa and van der Plicht Reference Kitagawa and van der Plicht1998; Bronk Ramsey et al. Reference Bronk Ramsey, Staff, Bryant, Brock, Kitagawa, van der Plicht, Schlolaut, Marshall, Brauer, Lamb, Payne, Tarasov, Haraguchi, Gotanda, Yonenobu, Yokoyama, Tada and Nakagawa2012). There were also new coral data from Tahiti (Durand et al. Reference Durand, Deschamps, Bard, Hamelin, Camoin, Thomas, Henderson, Yokoyama and Matsuzaki2013) and a “tie-pointed” record from the Pakistan margin (Bard et al. Reference Bard, Ménot, Rostek, Licari, Böning, Edwards, Cheng, Wang and Heaton2013; Heaton et al. Reference Heaton, Bard and Hughen2013). Marine reservoir ages were determined for each carbonate dataset from the offset with overlapping tree-ring data, where possible, and kept constant beyond the end of the tree rings. Comparison of Cariaco Basin and Barbados data to overlapping tree-ring data suggested an abrupt drop in reservoir ages during both the Younger Dryas and Heinrich Stadial, therefore those data were not used in IntCal13 and Marine13.
The random walk model for curve construction was developed further (Niu et al. Reference Niu, Heaton, Blackwell and Buck2013) to account for the dependencies in the calendar age estimates of the samples within both the Lake Suigetsu and the “tie-pointed” records; and the wiggle-matched tree-ring measurements. Discrepancies between the diverse datasets used to construct the curve between 13,900–50,000 cal BP were accommodated by allowing for a potential additive increase in the uncertainty of the data. This resulted in a smoother character for the IntCal13 mean curve in this older time period. The increase in the volume of data for 2013, together with the complex dependencies in the calendar age chronologies of the constituent records, meant that the random walk approach had reached its practical limit. It was decided that, for the future, a new method that enabled more flexibility and greater potential to investigate the effect of specific modeling and dataset choices was needed.
The SHCal13 curve was constructed with additional Holocene tree-ring data including a Tasmanian Huon pine floating chronology from 12,679 to 12,073 cal BP (Hua et al. Reference Hua, Barbetti, Fink, Kaiser, Friedrich, Kromer, Levchenko, Zoppi, Smith and Bertuch2009). For portions of the curve without Southern Hemisphere tree-ring measurements, the IntCal13 curve with an atmospheric radiocarbon offset between the Northern and Southern hemispheres (N-S offset) of 43 ± 23 yr was used, which extended SHCal13 to 50,000 cal BP (Hogg et al. Reference Hogg, Hua, Blackwell, Niu, Buck, Guilderson, Heaton, Palmer, Reimer, Reimer, Turney and Zimmerman2013).
2020 Calibration Curves
In 2013 IntCal focus groups were instigated to promote wider engagement with the IWG and provide additional expertise on various topics. As requested by the IWG, the statistical approach to curve construction was completely redesigned for the 2020 calibration curves (Heaton et al. Reference Heaton, Blaauw, Blackwell, Bronk Ramsey, Reimer and Scott2020a). The new methodology utilized Bayesian splines. MCMC was still used to provide the curves but the new Bayesian spline approach allowed much faster implementation and more reliable curve estimation than had been possible with the previous random walk.
The spline method of construction creates a calibration curve based upon a trade-off between the quality of the curve’s fit to the calibration samples, and its roughness/wiggliness. This penalty on roughness aims to prevent over-fitting and the introduction of potential spurious variability. Fit to the calibration data was assessed in F14C where measurement uncertainties are symmetric, while roughness was penalized in Δ14C space. The method was refined to include an additive error term to account for any potential additional variability in observed tree-ring 14C measurements and to explicitly recognize rapid atmospheric 14C change events such as the AD 774–775 event (Miyake et al. Reference Miyake, Nagaya, Masuda and Nakamura2012).
In the older portion of the curve, from 13,900–55,000 cal BP, much work was done by the group to resolve the discrepancies between the various speleothem, ocean, and lacustrine datasets (Butzin et al. Reference Butzin, Köhler and Lohmann2017; Heaton et al. Reference Heaton, Blaauw, Blackwell, Bronk Ramsey, Reimer and Scott2020a; Hughen and Heaton Reference Hughen and Heaton2020). To deal with any discrepancies that may remain, “heavy-tailed” errors were used to reduce the influence of possible outliers.
A total of 220 Holocene tree-ring datasets were considered for possible inclusion in the IntCal20 curve and screened for suitability (Reimer et al. Reference Reimer, Austin, Bard, Bayliss, Blackwell, Bronk Ramsey, Butzin, Cheng, Edwards, Friedrich, Grootes, Guilderson, Hajdas, Heaton, Hogg, Hughen, Kromer, Manning, Muscheler, Palmer, Pearson, van der Plicht, Reimer, Richards, Scott, Southon, Turney, Wacker, Adolphi, Büntgen, Capano, Fahrni, Fogtmann-Schulz, Friedrich, Köhler, Kudsk, Miyake, Olsen, Reinig, Sakamoto, Sookdeo and Talamo2020). Improvements in the linkages of the European Preboreal Pine and the Swiss chronologies were reinforced by the inclusion of single-year subfossil pine trees from the French Alps and allowed a fully atmospheric record to ca. 13,910 cal BP. The Hulu Cave speleothem dataset, which had been extended to 53,900 cal BP (Cheng et al. Reference Cheng, Edwards, Southon, Matsumoto, Feinberg, Sinha, Zhou, Li, Li and Xu2018), provided the backbone for the curve beyond the dendrochronologically dated and 14C wiggle-matched tree-ring measurements ca. 14,000 cal BP. Floating tree rings from Northern Italy and kauri from New Zealand provided higher resolution details to the curves further back in time. The Lake Suigetsu macrofossil data on an extended and improved varve timescale were also included (Bronk Ramsey et al. Reference Bronk Ramsey, Heaton, Schlolaut, Staff, Bryant, Brauer, Lamb, Marshall and Nakagawa2020). Because many of the U-Th coral datasets older than 25,000 cal BP exhibited rather high variability, these data were excluded from IntCal20. For marine data, rather than using a constant offset from the atmospheric curve as had been done in the past, the Hamburg Large Scale Geostrophic model was used to estimate the large variations in the marine reservoir offset over time for each region (Butzin et al. Reference Butzin, Heaton, Köhler and Lohmann2020) using the Hulu 14C record as input. For further details of corrections and adjustments to the datasets, see Reimer et al. (Reference Reimer, Austin, Bard, Bayliss, Blackwell, Bronk Ramsey, Butzin, Cheng, Edwards, Friedrich, Grootes, Guilderson, Hajdas, Heaton, Hogg, Hughen, Kromer, Manning, Muscheler, Palmer, Pearson, van der Plicht, Reimer, Richards, Scott, Southon, Turney, Wacker, Adolphi, Büntgen, Capano, Fahrni, Fogtmann-Schulz, Friedrich, Köhler, Kudsk, Miyake, Olsen, Reinig, Sakamoto, Sookdeo and Talamo2020) and references therein. The IntCal20 curve is more detailed than IntCal13, in part because of increased data coverage and in part due to the efforts made by the group to develop the methods that synthesized the diverse data. In particular, the work to resolve the differences between the various datasets from 13,900–55,000 cal BP has meant we are hopefully better able to discern common signal in this period. While there are some other differences between IntCal20 and IntCal13 in the Holocene, the major changes are from the Late Glacial through to the end of the curve at 55,000 cal BP. For example, radiocarbon ages calibrated with the IntCal20 curve between 34,000 and 42,000 cal BP calibrated ages may be up to 700 years older than with IntCal13 and between 42,000 to 50,000 cal BP these may be more than 1000 years younger. This has major consequences for scientific studies in this time range, for example examining the duration of the overlap between Neanderthals and anatomically modern humans in Europe (Bard et al. Reference Bard, Heaton, Talamo, Kromer, Reimer and Reimer2020).
The SHCal20 curve was reinforced with a number of new tree-ring datasets (Hogg et al. Reference Hogg, Heaton, Hua, Palmer, Turney, Southon, Bayliss, Blackwell, Boswijk, Bronk Ramsey, Pearson, Petchey, Reimer, Reimer and Wacker2020). For time periods where direct Southern Hemisphere measurements were not available, the IntCal20 curve was used with an N-S offset model (Heaton et al. Reference Heaton, Blaauw, Blackwell, Bronk Ramsey, Reimer and Scott2020a).
For the Marine20 surface ocean calibration curve the BICYCLE global carbon cycle model, which incorporates a number of components of the carbon cycle (Köhler et al. Reference Köhler, Muscheler and Fischer2006), was revised to allow atmospheric CO2 to be specified by an ice core-based reconstruction and atmospheric Δ14C by the IntCal20 curve (Heaton et al. Reference Heaton, Köhler, Butzin, Bard, Reimer, Austin, Bronk Ramsey, Grootes, Hughen, Kromer, Reimer, Adkins, Burke, Cook, Olsen and Skinner2020b). The resulting model produced a curve with variable marine reservoir offsets similar to the Hamburg Large Scale Geostrophic model. Marine20 is intended to represent the non-polar surface ocean as it has no provision for ice cover, however, Marine20 should be suitable for calibration of polar samples during the Holocene. It is worth noting that the previous model used for marine calibration curves (e.g., Marine86 through Marine13) also did not include ice cover. Because the parameters used in the model differ from those used in the previous marine curves, the average offset between the ocean and the atmosphere is larger. ΔR values used for regional corrections must be calculated with Marine20. Recalculated known age, pre-bomb ΔR values are given in the Marine Radiocarbon Reservoir Database (calib.org/marine) whereas ΔR values for paired marine and terrestrial samples can be calculated with the deltar software (Reimer and Reimer Reference Reimer and Reimer2017) available on the website.
OUTLOOK
The IntCal Working Group, which is now led by Christopher Bronk Ramsey, has already begun discussions on the next update. Kauri measurements through the Laschamp geomagnetic excursion ca. 41–42 ka cal BP have recently been published (Cooper et al. Reference Cooper, Turney, Palmer, Hogg, Mcglone, Wilmshurst, Lorrey, Heaton, Russell, Mccracken, Anet, Rozanov, Friedel, Suter, Peter, Muscheler, Adolphi, Dosseto, Faith, Fenwick, Fogwill, Hughen, Lipson, Liu, Nowaczyk, Rainsley, Bronk Ramsey, Sebastianelli, Souilmi, Stevenson, Thomas, Tobler and Zech2021) and there are new single year tree-ring data coming out soon. It is expected that the statistical methods will not have major modifications for the next update.
In the longer term, methods for calibrating sets of radiocarbon dates with a large subset of the possible realizations generated by the MCMC method are likely to improve “wiggle-matched” dates. As more tree-ring data from different regions with low measurement uncertainty become available it may be possible to provide region-specific terrestrial calibration curves. For the marine curve, the global carbon cycle model could be improved by incorporating time-varying changes to climatic parameters based on palaeoceanographic and palaeoclimatic records, and eventually, with sufficient data, we may be able provide location-specific curves, in particular for the polar regions (Heaton et al. Reference Heaton, Köhler, Butzin, Bard, Reimer, Austin, Bronk Ramsey, Grootes, Hughen, Kromer, Reimer, Adkins, Burke, Cook, Olsen and Skinner2020b).
ACKNOWLEDGMENTS
I would like to thank the Leverhulme Trust for a networking grant which enabled the initial IntCal Working Group meetings, the UK Natural Environment Research Council NE/E018807/1 for funding research leading to IntCal09 and IntCal13, and IGBP PAGES (Past Global Changes) for providing funding for two palaeoscientists on the IntCal Oversight Committee to participate in meetings leading to IntCal13. Pieter Grootes provided valuable suggestions on the manuscript. I am grateful to Tim Heaton for input on the statistical methods used in IntCal09 through IntCal20.