Abstract
Keywords
Introduction
Due to their annual resolution and absolute dating, tree-ring chronologies play an essential role in archaeological, climatic, and environmental research (Cook and Kairiukstis, 2013; Fritts, 1976). Further to dating relict wood remains (Bannister and Robinson, 1975; Tegel et al., 2022), dendro data are extensively used for reconstructing temperature and hydroclimate variability over centuries to millennia (e.g. Esper et al., 2016; Ljungqvist et al., 2020). Tree-ring records may also contain precise signatures of post-volcanic summer cooling (e.g. D’Arrigo et al., 2001; Esper et al., 2013), contribute to the interpretation of historical events (e.g. Büntgen et al., 2016; Oppenheimer et al., 2018), and refine the international radiocarbon (14C) calibration curves (IntCal; Reimer et al., 2020; Reinig et al., 2020). The quantity and quality of relict wood from different sources, however, decrease back in time, and limited tree-ring evidence is available for the early- and mid-Holocene (Bebchuk et al., 2024; Büntgen and Esper, 2025).
The Fenland region in eastern England comprises a vast, yet rapidly disappearing archive of subfossil wood (Figure 1; Bebchuk and Büntgen, 2025). Since the 1630s, this peat-rich flat lowland spreading at sea-level from Cambridge in the south to Lincoln in the north and from Peterborough in the west to the Wash Bay in the east, has been extensively drained for agricultural use (Darby, 1940). As peat thickness was reducing due to intensive farming, wind erosion and decomposition of organic materials, countless oak, yew, pine, willow, elm, birch, and ash trees came to the surface and have been excavated and subsequently burnt by generations of farmers (Astbury, 1958; Darby, 1983; Ennion, 1951; Godwin, 1978; Pryme, 1701; Pryor, 2019; Skertchly, 1877). Although the antiquity of trees was long recognized, their age remained unknown. English antiquary Abraham de la Pryme was the first to systematically describe these ‘subterraneous trees’ in 1701, arguing in his letter to the

Samples characteristics. (a) Photos of sampling sites with subfossil yews (left) and oaks (middle and right). (b) Sampling locations of subfossil yew (red) and oak (green) wood. (c) Temporal distribution of 100 yew tree-ring width (TRW, red bars) series provisionally radiocarbon (14C) dated and sorted by their first ring. Small dots indicate sample distribution when sorted by end dates. The 12 samples used for stable oxygen isotope (δ18O) analysis are shown in blue, of which eight passed the ISODATE statistical thresholds (solid line) and four did not pass (dashed line). Black dots mark the samples used for 14C analysis. Data are in calibrated years Before Present (BP) at an uncertainty of 7 years at the 95.4% probability. The green bar shows the span of the reference oak δ18O chronology.
Calendar dates of wood samples are usually established via cross-dating – a fundamental technique in dendrochronology that associates ring width patterns with absolute calendar dates (Douglass, 1941; Wigley et al., 1987). While the method works reliably at climatically extreme sites, it may fail where tree growth is not sufficiently stressed by a single dominant climate factor. The sparse distribution of subfossil wood finds worldwide, along with the limited number of contributing tree species further constrains the potential to develop absolutely dated tree-ring records over several millennia and from different regions and species (Tegel et al., 2022). To varying extents, most tree-ring chronologies for the early- and mid-Holocene therefore rely on provisional radiocarbon dating, which has resulted in important, yet floating tree-ring chronologies being developed for England (Batchelor et al., 2020), Greece (Christopoulou et al., 2024), Albania and North Macedonia (Hafner et al., 2021), Turkey (Pearson et al., 2020), Egypt (Kuniholm et al., 2014), and China (Shi et al., 2025), for example.
Tree-ring stable oxygen isotopes (δ18O) have been recently employed as an alternative dating tool (e.g. McCarroll et al., 2019; Nakatsuka et al., 2020; Römer et al., 2023). The δ18O in tree-ring cellulose is captured in wood without the need for plants to be climatically stressed as it is primarily controlled by the isotopic composition of source water and evaporative enrichment during transpiration (McCarroll and Loader, 2004; Siegwolf et al., 2022). As a result, a common δ18O signal is usually recorded even under favourable growth conditions at non-extreme sites with mild climates (Hartl-Meier et al., 2015; Loader et al., 2008). Another advantage of δ18O over traditional TRW cross-dating is that a relatively small sample replication and short measurement series can be sufficient to develop a reliable δ18O chronology (Loader et al., 2019). These benefits have facilitated cross-dating of δ18O records for different species (Loader et al., 2021) and broad geographical regions (Haneca et al., 2025; Shi et al., 2025), with promising applications for archaeological dating (Nayling et al., 2024a, 2024b; Sano et al., 2022).
Here, we present a TRW chronology from a 100 subfossil yew (
Data and methods
Between 2020 and 2024, we collected cross-sectional discs from over 400 subfossil yew trees in the Fenland region of eastern England (Figure 1a and b) and measured TRW along 2−5 radii per sample (for details, see Bebchuk et al., 2024). Visual cross-dating of all TRW measurements and statistical verification were performed using the TSAP-Win and COFECHA software, respectively (Holmes, 1983; Rinn, 1996). To provisionally date the yew record, 11 blocks of 10−20 tree rings from 8 samples were selected for 14C analysis. The samples were processed via the base-acid-base-acid method followed by bleaching and graphitization, using a high precision MICADAS facility at ETH in Zurich, Switzerland (Figure 1c; Němec et al., 2010; Sookdeo et al., 2020; Wacker et al., 2010). All 14C dates were wiggle-matched (Bronk Ramsey et al., 2001) and statistically resolved against IntCal20 (Reimer et al., 2020) using the OxCal v.4.4 software (Bronk Ramsey, 1995, 2021). The resulting 456-year-long TRW chronology comprises 100 yew samples (Supplemental Table S1) and spans from 4634 to 4179 ±7 calibrated years BP at the 95.4% probability (Figure 1c). Since no yew TRW chronology extends continuously from the present back into the mid-Holocene, we compared our record against an absolutely dated oak TRW chronology from the same region (Figure 1a), which was developed at Queen’s University in Belfast (QUB) in the 1980s (Baillie and Brown, 1988). We also compared our floating yew TRW chronology against English pine and German oak TRW chronologies (Billamboz, 2002, 2003; Boswijk, 1998). However, no statistically significant match was found with any of the existing chronologies.
To overcome these dating constraints, we analysed δ18O values from the tree rings of 12 yew samples (Figure 1c). Annual rings, including both early- and latewood, were split under a stereo microscope, with alpha-cellulose extracted following the modified Jayme-Wise isolation method (Boettger et al., 2007). Approximately 1 mg of the material was pyrolysed to carbon monoxide (CO) at 1450°C using a varioPYROcube elemental analyser (Elementar Analysensysteme, Germany), and the stable oxygen isotopes in the CO gas were measured using an ISOPRIME100 continuous-flow Isotope Ratio Mass Spectrometer (IRMS; IsoPrime, Manchester, UK) at the Global Change Research Institute in Brno, Czech Republic. Prior to each set of analyses, the IRMS was tuned and tested for signal stability (standard deviation ⩽ 0.04‰ over 10 pulses of reference gas) and linearity (⩽0.03‰/nA) across the expected ion current range obtained from the measurements of the test samples. Precision was maintained with standard deviations ⩽0.10‰ based on 5 consecutive measurements of the same alpha-cellulose sample. Isotopic values were calibrated against certified reference materials from the International Atomic Energy Agency (IAEA, Vienna, Austria). The δ18O values were referenced to benzoic acids (IAEA-601 and IAEA-602), and are reported in permil (‰) relative to Vienna Standard Mean Ocean Water (VSMOW; Coplen, 1995).
The individual yew δ18O series were cross-dated using the novel ISODATE software package (Davies et al., 2025). After indexation with a 9-year rectangular filter by subtraction, the series were first compared against each other and then against a reference oak δ18O chronology that is currently under development at Swansea University, UK.
To date, the reference oak δ18O chronology is anchored using ring width dendrochronology and has a near continuous coverage between 2982 and 1954 BCE (datum considers the year zero; Büntgen and Oppenheimer, 2020). It consists of 15 oak samples, out of which 11 form part of the absolutely dated TRW chronology developed at QUB by Baillie and Brown (1988), and the other four samples were recently included to the chronology to increase its sample replication (Supplemental Figure S1). These samples were cross-dated with the QUB samples using ring-width or stable oxygen isotope dendrochronology (Supplemental Table S2). The chronology’s coverage is being constantly extended as further samples from the QUB English chronology, that had precise dates previously assigned to them, are processed and added to the developing record. All reported dates are currently dependent on the dendrochronological dating of the QUB English archive. Once complete, the third millennium BCE isotope reference chronology will be made publicly available. This will include all cross-matching statistics between the isotope series used in its construction. Additionally, the series will be made publicly accessible via ISODATE (Davies et al., 2025) for dating purposes. For the stable oxygen isotope analysis, the α-cellulose was extracted only from the latewood component of each oak sample, homogenized and freeze-dried (Loader et al., 1997; Wieloch et al., 2011). It was then weighed into silver capsules and pyrolized to carbon monoxide gas at 1400°C. The δ18O were measured on an isotope ratio mass spectrometer at Swansea University, UK. The measurements are reported in permil (‰) relative to the VSMOW standard (Coplen, 1995), with a typical analytical precision (σn−1) of ±0.3 ‰ (Loader et al., 2013).The 15 oak δ18O series were indexed using a 9-year rectangular filter, and once securely cross-dated, averaged to produce the dating reference record. Following the statistical framework proposed by Loader et al. (2019), the agreement between the δ18O series was assessed using Pearson’s correlation coefficients
Results
Independent cross-dating of the 12 yew δ18O series against each other positions 8 of them in agreement with the original TRW-based cross-dating results. The 8 series correlate with each other with Pearson’s

Cross-dating characteristics. (a) Eight individual yew δ18O series between which a statistically significant match is found, and that are used to develop a mean yew δ18O chronology, all detrended with a 9-year rectangular filter. Note that the gaps in data are due to a loss of samples at the lab. (b) Yew (blue) and reference oak (green) δ18O chronologies, both
Intriguingly, our study presents the first example of oxygen isotopic agreement between oak and yew wood. We further emphasize that tree-ring stable isotopes can help overcome challenges associated with traditional TRW cross-dating, especially if a species has an irregular growth pattern, or tree growth is influenced by multiple climatic drivers and species exhibit different climate sensitivities. This approach can significantly advance archaeological and palaeoclimatic research, which are strongly limited by the availability of relict material.
Our absolute dating supports the results of the initial radiocarbon analysis (Figure 3, Supplemental Table S4). The calendar dates newly assigned to the 11 samples that were used for 14C measurements are within the 99.7% probability range given by wiggle-matching (Bronk Ramsey et al., 2001). Although the 14C values are slightly higher than the IntCal20 curve, especially across the plateau between circa 2450 and 2350 BCE, this discrepancy likely arises from both, samples integrating over 10–15 years and the calibration curve itself lacking a high temporal precision. Our study emphasizes the importance of increasing the total number of absolutely dated and annually resolved Holocene tree-ring chronologies, which may further help constrain radiocarbon calibration curves (Pearson et al., 2022).

Radiocarbon and absolute dating. Raw results of radiocarbon dating are listed in the first column of the inset table (for full statistics, see Supplemental Table S4). The individual 14C dates are resolved against the most recent radiocarbon calibration curve for the Northern Hemisphere IntCal20 (in blue; Reimer et al., 2020) with their probability distributions shown in light grey and calibrated ages BCE given in the second column. To further reduce the temporal uncertainty, the dates are wiggle-matched using dendrochronological information (dark grey shadings and white dots for mean values, third column). Note that the 14C dates are integrated over 10–15 years, and this uncertainty is not taken into account in OxCal calibration (Bronk Ramsey, 2001, 2021). The calibrated dates are reported at the 95.4% probability. Absolute cross-dating of the yew δ18O chronology against the reference oak δ18O chronology assigns a calendar year BCE to each sample (red dots, last column). The difference between wiggle-matched and calendar dates (white and red dots, respectively) is 17 years.
The statistically significant inter-species isotopic agreement now places our yew TRW chronology precisely from 2668 to 2213 years BCE (Figure 4a). By establishing an absolute date, we secure a benchmark for further improvement of Fenland dendrochronology, add a new species into the limited list of high-resolution Holocene chronologies, and enable additional investigations of long tree-ring chronologies from Ireland (Baillie and Brown, 1988; Brown et al., 1986), Fennoscandia (Eronen et al., 2002; Grudd et al., 2002; Helama et al., 2008), and western Europe (Boettger et al., 2007; Eckstein et al., 2009; Leuschner et al., 2002).

Chronology characteristics and applications. (a) Absolutely dated yew TRW chronology in red and its sample replication in black. EPS is the expressed population signal and
Discussion
Spanning the transition from the Late Neolithic to the Early Bronze Age, our absolutely dated yew chronology offers new opportunities for refining archaeological interpretation (Figure 4b). Systematic archaeological surveys in eastern England have unveiled vast and well-preserved evidence of early human activity (Godwin, 1997; Hall and Coles, 1994). Numerous houses, causewayed enclosures, trackways, and burial round borrows have been found (Medlycott 2011). Moreover, the period is characterized by the spread of metal-work and a remarkably wide range of pottery styles across the Fenland attributed to the Peterborough, Grooved Ware, and Beaker cultures (Bullmore, 2022; Glazebrook, 1997; Yates and Bradley, 2010). Despite a seemingly uniform landscape, regional diversity in archaeological findings has always been recognized, and developing a more precise spatiotemporal framework has consistently been identified as a research priority (Brown and Glazebrook, 2000; Cooper, 2024; Glazebrook, 1997; Medlycott, 2011. The calendar-dated yew TRW record now offers an opportunity for precise archaeological attribution.
The mid-Holocene hydroclimate reconstruction recently developed from Fenland yew wood (Bebchuk et al., 2025) is now anchored to the period 2645−2228 BCE. Further to shedding new light onto the hydro-geomorphological and biogeographic consequences of sea-level changes in the North Sea (Figure 4b) and their possible impacts on human migration and societal transformation, the record secures a timescale for high-resolution palaeoclimatic research on the still debated climate anomaly around 4200 years ago, the so-called ‘4.2 ka event’ (Helama, 2024; Weiss et al., 1993). Associated with an extreme drought in central Asia, its effect and extent across the North Atlantic remain unclear (Bradley and Bakke, 2019; McKay et al., 2024). Several lines of evidence point to increased storminess (Goslin et al., 2019; Sorrel et al., 2012) and a sea-level rise in the North Sea (Bebchuk et al., 2024, 2025; Shennan et al., 2018; Waller, 1994), which could have shifted (or be triggered by the shift of) the Inter-Tropical Convergence Zone southwards and contributed to drying across parts of Asia (Yan and Liu, 2019) and the European Alps (Arosio et al., 2025). We expect our yew archive to serve as a reference in advancing the quest for the causes and mechanisms of the 4.2 ka event.
Furthermore, our new absolutely dated yew chronology may help establish calendar dates for Icelandic volcanic eruptions, such as those of Katla and Hekla 4 that are currently estimated at 2388 ± 8 and 2375 ± 8 BCE, respectively (Figure 4b; Davies et al., 2024). This far distant dating might be possible through wood anatomical investigations of ‘blue rings’ that are known to occur in coniferous species and to be indicative of abrupt summer cooling often caused by large, sulphur-rich volcanic eruptions (Büntgen et al., 2022; Piermattei et al., 2015). We also expect many absolutely dated and annually resolved radiocarbon dates from our still growing yew archive to contribute to the refinement of the next IntCal curve (Reimer et al., 2020), and thus help reconstruct changes in solar activity during the mid-Holocene (Figure 4b). In the light of these opportunities, we hope that our study will stimulate measurements of tree-ring stable isotopes in relict wood across different species and regions and showcase their potential for advancing archaeological and paleoenvironmental research worldwide.
Supplemental Material
sj-docx-1-hol-10.1177_09596836251407634 – Supplemental material for An absolutely dated mid-Holocene English yew chronology offers new opportunities for archaeological and palaeoenvironmental research
Supplemental material, sj-docx-1-hol-10.1177_09596836251407634 for An absolutely dated mid-Holocene English yew chronology offers new opportunities for archaeological and palaeoenvironmental research by Tatiana Bebchuk, Darren Davies, Neil J Loader, Otmar Urban, Tito Arosio, Lukas Wacker, Natálie Pernicová, Josef Čáslavský, Miroslav Trnka, Alexander V Kirdyanov, David Brown, Jan Esper, Clive Oppenheimer and Ulf Büntgen in The Holocene
Footnotes
Acknowledgements
The authors thank Inna Roshka for her work at the Czech isotope laboratory, Danny McCarroll and the QUERCUS project team for support, Adamson Low and Martin Bridge for sample provision and dating, and Iona Robinson Zeki and David Matzliach from the Cambridge Archaeological Unit for providing the photo of wood remains at the archaeological site.
Author contribution(s)
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: TB was funded by the Hill Foundation Cambridge Trust, DD and NJL are supported by UKRI EP/X025098/1, TA by SNSF P500PN_225686, and JE and UB received funding from the Czech Science Foundation (# 23-08049S; Hydro8) and the ERC Advanced Grant (# 882727; Monostar). UB received the ERC Synergy Grant (# 101118880; Synergy-Plague). Stable isotope analyses were performed using analytical infrastructure that is part of AnaEE Czech Republic and AnaEE-ERIC networks (
), supported by the Ministry of Education, Youth, and Sports of the Czech Republic through the CzeCOS programme, grant number LM2023048.
Data availability
Raw yew tree-ring width (TRW) measurements are available at the International Tree-Ring DataBase (ITRDB) at https://www.ncei.noaa.gov/access/paleo-search/study/41020, and raw yew tree-ring stable oxygen isotope (TRSOI) measurements are published at the NOAA for Palaeoclimatology data repository at https://www.ncei.noaa.gov/access/paleo-search/study/41021. Oak stable isotope measurements form part of an evolving 7000-year chronology for the UK, which is currently being developed at the Swansea University, UK. Once complete, the data will be publicly available. Radiocarbon dates are provided in the Supplementary material. The ISODATE software is freely available at
.
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
