INTRODUCTION
The catastrophic Minoan eruption of Santorini (Thera) in the second millennium BCE provides a crucial chronological anchor for Bronze Age prehistory. The precise date of the eruption has been debated for decades (for a recent overview see Antiquity 2014). Radiocarbon (14C) dates play a major role in this discussion. These need to be calibrated to derive historical dates. Obviously, historical inference largely depends on the exact shape of the calibration curve. Recently, a new single year calibration record became available, its content indicating that the curve needs to be revised (Pearson et al. Reference Pearson, Brewer, Brown, Heaton, Hodgins, Jull, Lange and Salzer2018). This spawned major (re)dating efforts of dendrochronologically dated wood dating to the time of the Santorini eruption. Here, we present a short report on 18 such new dates from the Netherlands.
METHODS AND RESULTS
In 1996, four trees were found during infrastructural work by Holland Railconsult (now: Movares) in Bodegraven, the Netherlands. The approximate coordinates of the location are latitude 52.0823259, and longitude 4.7460844. The tree species is Quercus robur/petraea. The trees were analyzed and two of them were dated at the Netherlands Centre for Dendrochronology/RING Foundation in Amersfoort. One of the dated trees (coded BOF00041) has a growth pattern covering the time of the Santorini eruption. Hence, this wood has been selected presently for dating by 14C, for cross checking the calibration curve.
Dendrochronology
The tree sample BOF00041 contains a total of 186 rings. It does not contain sapwood, implying that at least a sapwood zone covering 26 ± 8 yr is lacking on the outside of the wood. The tree-rings were dated to 1737–1552 BCE (dataset developed by Hanraets and Jansma 1997, see Jansma et al. Reference Jansma, van Lanen, Brewer and Kramer2012). The reference chronologies used are from North Germany (Leuschner and Delorme Reference Leuschner and Delorme1988; Leuschner, unpublished data) and the Netherlands (‘NLPre_ZH; Jansma Reference Jansma1995). In general, narrow rings in the bog oaks from the Netherlands correspond very well with these from North Germany. Examples have been published in Leuschner et al. (Reference Leuschner, Sass-Klaassen, Jansma, Baillie and Spurk2003). The dendrochronological parameters are shown in Table 1.
The parameter %PV is percentage of parallel variaton, also known as “Gleichlaufigkeit”. It is the percentage of rings showing simultaneous increases and decreases of tree-ring width relative to the reference chronology (Jansma Reference Jansma1995). The parameter t is the result of a Student’s t-test for the Pearson’s cross-correlation coefficient between the tree-ring pattern and that of the reference chronology. The parameter P is the probability that the %PV value is coincidental, expressed as a fraction of 1.
Radiocarbon
The two series of tree rings were dated by AMS in Groningen. For 14C dating we selected 9 decadal- and 9 single-year wood samples, as indicated in Table 2. All the samples were pretreated to α-cellulose using the method of Groningen (Dee et al. Reference Dee, Palstra, Aerts-Bijma, Bleeker, de Bruijn, Ghebru, Jansen, Kuitems, Paul and Richie2020). In brief, the tree-rings (early- and late-wood not separated) were cut into much smaller fragments with a scalpel. Aliquots of ~50 mg (or all available) of the fragments were weighed into 12-mL test tubes. Then, the samples were exposed to a strong acid (HCl, 1.5 M, 80ºC, 20 min), strong base (NaOH, 17.5% w/vol, 60 min, RT), acid (HCl, 1.5 M, 80ºC, 20 min), and finally to strong oxidation (NaClO2, 1.5% w/vol in HCl, 0.06 M, 16 hr, 80ºC), with rinses using deionized, ultrapure water after each chemical step.
The alpha cellulose was combusted to CO2 by an elemental analyzer, connected to a stable isotope mass spectrometer (EA/IRMS, Elementar Vario Isotope Cube™/Isoprime 100™). The latter provides the stable isotope ratio δ13C (in ‰, relative to the VPDB standard; Mook Reference Mook2006). Part of the CO2 is transferred into graphite, by a reaction with H2 gas at a temperature of about 600ºC, using Fe powder as catalyst (Aerts-Bijma et al. Reference Aerts-Bijma, van der Plicht and Meijer2001). The graphite was pressed into target holders for the ion source of the AMS. The AMS is a MICADAS-17 (IonPlus®) (Mini Carbon Dating System; Synal et al. Reference Synal, Stocker and Suter2007) manufactured by IonPlus, installed in 2017. The present Groningen laboratory code is GrM.
The results are shown in Table 2. Listed are the GrM-number, the dendrochronological date in BCE, the ring width, the radiocarbon date in BP, its uncertainty (1 σ), and the δ13C value (in ‰) as determined by the IRMS.
The uncertainty of the latter is 0.15‰ (1 σ). These δ13C values are all consistent with those of the AMS.
The radiocarbon dates are reported by convention (van der Plicht and Hogg Reference van der Plicht and Hogg2006), using the oxalic acid reference, the conventional half-life and isotopic fractionation correction using the δ13C measured by the AMS. The uncertainties of the 14C dates are based on counting statistics and includes an estimate for internal laboratory error.
DISCUSSION AND CONCLUSION
The 14C dates (Table 2) of the selected samples from bog oak BOF00041 are shown in Figure 1a. Also shown is the relevant part of the calibration curve IntCal13 (Reimer et al. Reference Reimer, Bard, Bayliss, Beck, Blackwell, Bronk Ramsey, Buck, Cheng, Edwards and Friedrich2013). Also, the individual data used for construction of the curve are shown. These are high-precision radiometric dates from the laboratories Heidelberg (Hd), Belfast (UB) and Seattle (QL). The error bars plotted are all 1 σ. For the IntCal curve, the 1-σ envelope is shown in Figure 1a.
As visible in Figure 1a, the GrM-dates are fairly consistent with the IntCal13 curve. Indeed, most of our new data overlap within 1 σ with IntCal13 but a few data points are slightly increased compared to IntCal13. Therefore, our data confirm that the calibration curve needs to be raised for the time range around 1600 BCE, but not as dramatically as suggested by Pearson et al. (Reference Pearson, Brewer, Brown, Heaton, Hodgins, Jull, Lange and Salzer2018), see Figure 1b. Our data will be implemented in the new calibration curve IntCal20 (Reimer et al. Reference Reimer, Austin, Bard, Bayliss, Blackwell, Bronk Ramsey, Butzin, Edwards, Friedrich and Grootes2020 in this issue), together with denser (long single-year series) datasets which became available as well (this issue). Based on all available new data, the implications of calibration for the Santorini eruption will be discussed in more detail elsewhere (van der Plicht et al. Reference van der Plicht, Bronk Ramsey, Heaton, Scott and Talamo2020 in this issue).
ACKNOWLEDGMENTS
This work was supported by an European Research Council grant (ECHOES, nr. 714679).