1. Introduction
The tropics are the most biologically diverse and important areas in the world (Ter Steege et al. Reference Ter Steege, Pitman, Killeen, Laurance, Peres, Guevara, Salomão, Castilho, Amaral, de Almeida Matos and de Souza Coelho2015), yet radiocarbon (14C) measurements on tropical tree rings are relatively sparse. Recent advances in tropical dendrochronology have shown that several tropical trees can produce recognizable annual bands due to cambial dormancy, driven by annual seasonality in environmental conditions (Brienen et al. Reference Brienen, Schöngart, Zuidema, Goldstein and Santiago2016). As wood cells develop and fix 14CO2 from the atmosphere during the growing season, tree rings can reflect the atmospheric 14C signal at the time the tissue was formed. Provided that annually resolved tree rings have been found and were validated by measuring 14C of selected calendar years after 1950 AD (which provides a high time resolution reference of ≤ 1 year—Fontana et al. Reference Fontana, López, Santos, Villalba, Hornink, Assis-Pereira, Roig and Tomazello-Filho2024; Hadad et al. Reference Hadad, Santos, Junent and Grainger2015; Haines et al. Reference Haines, Olley, English and Hua2018; Pearson et al. Reference Pearson, Hua, Allen and Bowman2011; Santos et al. Reference Santos, Granato-Souza, Barbosa, Oelkers and Andreu-Hayles2020, Reference Santos, Rodriguez, Barreto, Assis-Pereira, Barbosa, Roig and Tomazello-Filho2021), they can be further used to address observation-based atmospheric 14C records. This is relevant, due to the shortage of robust atmospheric 14C records across the tropics (Santos et al. Reference Santos, Granato-Souza, Barbosa, Oelkers and Andreu-Hayles2020).
Global atmospheric 14C concentrations post-AD 1950 are typically divided into five geographical zones (Hua et al. Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and Palmer2022), e.g., three across the Northern Hemisphere (NH1, NH2, and NH3), and two over the Southern Hemisphere (SH1-2 and SH3). As part of a major effort to produce and integrate high quality annual tree-ring records of past atmospheric 14C into future age calibration compilations (Hua et al. Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and Palmer2022), researchers have produced novel 14C records using dated tree rings from tropical and subtropical South America (Ancapichún et al. Reference Ancapichún, De Pol-Holz, Christie, Santos, Collado-Fabbri, Garreaud, Lambert, Orfanoz-Cheuquelaf, Rojas, Southon and Turnbull2021; Santos et al. Reference Santos, Linares, Lisi and Tomazello Filho2015, Reference Santos, Albuquerque, Barros, Ancapichún, Oelkers, Andreu-Hayles and De Faria2022, Reference Santos, Granato-Souza, Ancapichún, Oelkers, Haines, De Pol-Holz, Andreu-Hayles, Hua and Barbosa2024). Still, the methodology is intrinsically challenging, as it demands costly and laborious work. For example, 14C measurements using annual tree rings involve a sequence of i) sampling material from radial directions removed from cross-sections or increment cores in the field, ii) dating the trees and marking the selected tree rings for later wood subsampling, iii) careful isolation of tree rings with a scalpel, iv) chemical procedures to isolate inert carbon fractions (such as holocellulose and alpha-cellulose extracts), v) sample processing to CO2 and/or filamentous graphite, and finally vi) analysis by specific analytical instruments (such as, accelerator mass spectrometry—AMS).
Each of those steps can introduce uncertainties to the final results. Santos et al. (Reference Santos, Linares, Lisi and Tomazello Filho2015) and Turnbull et al. (Reference Turnbull, Mikaloff Fletcher, Brailsford, Moss, Norris and Steinkamp2017) reported differences in 14C data replication and attributed them to the difficulty of separating tree rings at their boundaries. Svarva et al. (Reference Svarva, Grootes, Seiler, Stene, Thun, Værnes and Nadeau2019) went a step further, adding that inhomogeneities in the tree-ring cellulose extracts can also lead to differences in 14C data replication.
In principle, 14C data variance can be extremely high on single tree-ring analysis when wood material has been mistakenly added, or accidentally lost during the cutting up of wood sections for further 14C analysis. Specifically, tree-ring subsampling from laths can experience both exclusion of portions of the desired growing season, or unintentional inclusion of more than one season. In the laboratory, caution is also required during single tree-ring reduction to particles (a maximum % yield recovery is desirable), and/or cellulose partitioning for the CO2 combustion step before the 14C-AMS measurements take place (use of homogenized fibers is advised). Radiocarbon variance among tree-ring replicated results tend to be higher during the increase of atmospheric 14C levels after nuclear bomb tests in the early 1950s (Enting Reference Enting1982), and its subsequent decline due to the 1963 moratorium on above ground testing (Levin and Hesshaimer, Reference Levin and Hesshaimer2000). During those periods, the relative level of 14C absorbed by wood material changes greatly from the onset to the end of wood cell formation. Thus, one should expect differences in 14C signatures when the wood material analyzed is imprecisely subsampled relative to the distinct lines (ring boundaries or earlywood/latewood boundaries) that demarcate a single tree ring (Turnbull et al. Reference Turnbull, Mikaloff Fletcher, Brailsford, Moss, Norris and Steinkamp2017). The effect of accidentally losing wood material during the early stages of sample handling can mimic the 14C results of selective earlywood or latewood measurements (e.g., Haines et al. Reference Haines, Hiscock, Palmer, Turney, Thomas, Cadd, Vohra and Marjo2023; Kudsk et al. Reference Kudsk, Olsen, Nielsen, Fogtmann-Schulz, Knudsen and Karoff2018).
As instrumentation has advanced, especially AMS systems used in 14C analysis (Bronk Ramsey Reference Bronk Ramsey2023), analytical uncertainties have decreased to well under 0.3% (Sookdeo et al. Reference Sookdeo, Kromer, Büntgen, Friedrich, Friedrich, Helle, Pauly, Nievergelt, Reinig, Treydte and Synal2020, for example). Further evaluations of sources of errors in tree-ring 14C data are still required. Notably, it is difficult to separate cause and effect when so many steps are involved in tree-ring sample processing for isotopic analysis, as aforementioned. Overall, differences in chemical extractions and spectrometers may not be the only factors causing data variability. Several researchers have compared 14C data from distinct chemical methods, e.g., holocellulose versus alpha-cellulose (Baker et al. Reference Baker, Santos, Gloor and Brienen2017; Capano et al. Reference Capano, Miramont, Guibal, Kromer, Tuna, Fagault and Bard2018; Cercatillo et al. Reference Cercatillo, Friedrich, Kromer, Paleček and Talamo2021; Michczyńska et al. Reference Michczyńska, Krąpiec, Michczyński, Pawlyta, Goslar, Nawrocka, Piotrowska, Szychowska-Krąpiec, Waliszewska and Zborowska2018; Němec et al. Reference Němec, Wacker, Hajdas and Gäggeler2010; Santos et al. Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023; Southon and Magana Reference Southon and Magana2010), independent chemical treatments of alpha-cellulose (Gaudinski et al. Reference Gaudinski, Dawson, Quideau, Schuur, Roden, Trumbore, Sandquist, Oh and Wasylishen2005; Santos et al. Reference Santos, Granato-Souza, Ancapichún, Oelkers, Haines, De Pol-Holz, Andreu-Hayles, Hua and Barbosa2024), extractions using different apparatus or devices (Fogtmann-Schulz et al. Reference Fogtmann-Schulz, Kudsk, Adolphi, Karoff, Knudsen, Loader, Muscheler, Trant, Østbø and Olsen2021; Santos et al. Reference Santos, Granato-Souza, Barbosa, Oelkers and Andreu-Hayles2020; Southon and Magana Reference Southon and Magana2010), as well as different types of mass spectrometers of laboratory analysis (Baker et al. Reference Baker, Santos, Gloor and Brienen2017; Kudsk et al. Reference Kudsk, Olsen, Hodgins, Molnár, Lange, Nordby, Jull, Varga, Karoff and Knudsen2021; Santos et al. Reference Santos, Albuquerque, Barros, Ancapichún, Oelkers, Andreu-Hayles and De Faria2022, Reference Santos, Granato-Souza, Ancapichún, Oelkers, Haines, De Pol-Holz, Andreu-Hayles, Hua and Barbosa2024; Sookdeo et al. Reference Sookdeo, Kromer, Büntgen, Friedrich, Friedrich, Helle, Pauly, Nievergelt, Reinig, Treydte and Synal2020). Regardless of differences in approach, 14C results have been somewhat similar within uncertainties (± 2σ). For example, holocellulose versus alpha-cellulose extractions of Holocene samples among distinct laboratories yielded pooled standard deviations of 0.25% or lower (n = 5 pairs/sets) (Santos et al. Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023); alpha-cellulose extractions from independent chemical procedures and apparatus yielded standard deviations of 0.3% (n = 27 pairs/sets; Santos et al. Reference Santos, Albuquerque, Barros, Ancapichún, Oelkers, Andreu-Hayles and De Faria2022); and different 14C-AMS facilities’ alpha-cellulose procedures and 14C analysis yielded standard deviations of 0.17% (n = 10 pairs/sets; Santos et al. Reference Santos, Granato-Souza, Ancapichún, Oelkers, Haines, De Pol-Holz, Andreu-Hayles, Hua and Barbosa2024). Here, pooled standard deviation values were calculated using the formulation shown in Cohen (Reference Cohen1988), as a way to assess data quality.
In order to separate spatial and temporal heterogeneity in the tree-ring samples, from chemical and/or analytical uncertainties, Santos et al. (Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023) focused on evaluating the very early stages of wood sample handling as well as cellulose extract homogeneity. Depending on its growth triggers and sites (Baker et al. Reference Baker, Santos, Gloor and Brienen2017), different tree species may experience different growing season lengths. Therefore, Santos et al. (Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023) emphasize the inclusion of earlywood and latewood from each full tree-ring sample analyzed. This can be done through the microscope inspection of tree-ring samples before samples are reduced to small particles (e.g., wood chips or slivers) for chemical processing. Variance-reduction strategies, as reported in Santos et al. (Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023), have then been used in the production of two post-AD 1950 atmospheric 14C records across the Amazon when using single absolutely dated trees (e.g., tree-ring 14C records produced from Central and Eastern Amazon Basin) (Santos et al. Reference Santos, Albuquerque, Barros, Ancapichún, Oelkers, Andreu-Hayles and De Faria2022 and Reference Santos, Granato-Souza, Ancapichún, Oelkers, Haines, De Pol-Holz, Andreu-Hayles, Hua and Barbosa2024, respectively). Nonetheless, just one dendrochronologically-dated tree was used per site studied.
The goal of this work is threefold: i) to reevaluate the aforementioned variance-reduction strategies in multiple tree species; ii) test data consistency among several trees on a given site (same species or different species); and iii) evaluate data consistency among trees based on their distances apart on a given site. For all of those, we present 14C tree-ring data of two post-AD 1950 Pantropical sites and their coexisting trees (same or distinct species). Those sites and samples pertain to a project to complete post-1950 14C records with the goal of mitigating important deficiencies of the current 14C atmospheric global map over the Tropical Low-Pressure Belt (TLPB; conventionally termed the Inter-Tropical Convergence Zone). Those records are timely (Reimer et al. Reference Reimer, Austin, Bard, Bayliss, Blackwell, Ramsey, Butzin, Cheng, Edwards, Friedrich and Grootes2020), given the dearth of datasets that can elucidate air transportation and climatic variability affecting atmospheric 14C concentrations across this region. Thus here, we took advantage of the 14C pairs/sets produced from the tree-ring calendar years associated with the post-1950 slopes, i.e., periods of major atmospheric 14C signature changes where small offsets can be highly detectable, to identify methods to improve precision and accuracy. These points are paramount for the assessment of the likely impact of our soon-to-be published Pantropical 14C records. As such, climatic patterns and geographic provenance of air-masses to sites are not pursued in depth in the present article, as they are not in the scope of this work. As already mentioned, to evaluate the reproducibility of the 14C data, we made use of the pair/sets of the overlapping calendar years, i.e., replicated data from each atmospheric 14C record at a given location. For accuracy data discussions, we use the current atmospheric zonal 14CO2 curves (such as SH1-2 and SH3, as defined in Hua et al. Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and Palmer2022), when using multiple tropical trees at single locations.
2. Materials and methods
2.1 Sites and samples
The 14C data discussed in this work are from samples collected at two localities in South America. Vegetation growth in South America is strongly regulated by the South American Monsoon System (SAMS), occurring between September and May (Novello et al. Reference Novello, Cruz, Vuille, Stríkis, Edwards, Cheng, Emerick, De Paula, Li, Barreto and Karmann2017). Ring growth, however, follows more specific site climatic conditions (e.g., the Köppen climate classification of each specific region – Figure S1). Thus, our comparisons explored i) tropical dry and moist forest conditions (Figure 1); ii) more than one tree and/or species per site, where trees were collected from dead (cross-section) and/or live (increment-core) individuals; and iii) source areas, as trees were separated by a range between 0.2 and 22.5 km (Figure 2).
Map of South America, detailing the two sample sites in Brazil (A), i.e., Juvenília, Minas Gerais (blue star) and Jamari National Forest, Rondônia (red star). Aerial view of the Juvenília site, and tree locations for the Cedrela fissilis samples, MTV001D (cross-section) and MTV114B (increment-core) (B). Aerial view of the Jamari National Forest site (C) tree positions for the Cedrela fissilis JAC09C and JAC21A, and Peltogyne paniculata JAR26C. Equator (black line) and the Dec/Jan/Feb-TLPB delimiting mean position (white dashed-line; adapted from Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and PalmerHua et al (2022)) is shown here for reference. Further details on wood samples characteristics are displayed in Figure 2.
Types of wood samples used in this study (A) from Juvenília, Minas Gerais, and (B) from Jamari National Forest, Rondônia dendrochronology field collections. Site coordinates, tree species, field deployable sampling method and sample identification are shown under each image. Examples of basic tree-ring sampling methods are shown in Figure S2, in supplementary material.
Trees from the Rio São Francisco Basin, specifically those from Juvenília, Minas Gerais, Brazil (blue star in Figure 1A) experience a long dry season followed by a wet season of heavy rainfall. Since Cedrela species’ tree growth and reproduction are normally synchronous with the onset of the rains, we used a climate classification diagram to determine the period of tree growth across this region (November to March; Figure S1).
Trees at the Juvenília site consisted of dendrochronologically correlated Cedrela fissilis trees (Pereira et al. Reference Pereira, Barbosa, Torbenson, Stahle, Granato-Souza, Santos and Barbosa2018). The series MTV001D and MTV114B (Figure 2), originating from MTV001 (1882 to 2015) and MTV114 (1994 to 2015) trees, respectively, were selected for 14C dating. In the first tree, ring widths were measured on four radii (cross section), while in the second, ring widths were measured on two radii (cores). The intercorrelation between the two series was r=0.692 (Table S1). Even so, MTV001 shows some evidence of wedging rings (Figure 2), protocols were followed to select the sample materials for optimal 14C results (i.e., laths that show rings with the earlywood and latewood boundaries clearly identified). The rationale of using both series, MTV001D and MTV114B, during the production of a complete post-1950 tree-ring 14C record was associated with the compressed section near the bark of the MTV001 tree (Figure 2). While the compressed section did not pose a problem during dendrochronological dating, it did not allow for precise sampling of wood for 14C analysis, as the rings in this section are far too narrow. Thus, we selected MTV114B series to complete the calendar years of the 14C record (not shown in this work). We then arranged to overlap the two wood sample sequences for 14C dating at the post-bomb slope section. This procedure assured us that high-quality data could be produced via several 14C replicates. The 15 overlapping calendar years selected for 14C analysis included duplicates of the single tree rings sampled from the cross-section MTV001D (the calendar years of 1964, 1966, and 1969), the increment core MTV114B (the calendar years of 1982, 1988, and 1989), and several others from both samples (nine calendar years between 1971 and 1980).
For the Jamari National Forest dendrochronology collections in the Southwestern Brazilian Amazon (red star in Figure 1A), three trees were selected for 14C analysis (Figure 2B). Trees at this locality experience a hot, humid, and seasonally wet type of climate (precipitation occurs between October and May; Figure S1). Trees species of the dendrochronologically dated collections of Peltogyne paniculata and Cedrela fissilis were selected for our study (Santos et al. Reference Santos, Rodriguez, Barreto, Assis-Pereira, Barbosa, Roig and Tomazello-Filho2021 and Ortega-Rodriguez et al. Reference Ortega-Rodriguez, Sánchez-Salguero, Hevia, Granato-Souza, Cintra, Hornink, Andreu-Hayles, Assis-Pereira, Roig and Tomazello-Filho2023, respectively). Here, we took advantage of those two annually-resolved tree species within a mutual territory to test for equal F14C means (Figure 1C), assuming that if they are subjected to the same 14CO2 air-masses they should have the same F14C value within error. In this case, we assume that P. paniculata tree species responses to the rainy season would be similar to those of Cedrelas at this location, even though their phenology has not been extensively examined yet.
From this site, 36 overlapping calendar years between 1955 and 1990 were tested. In all cases, at least two different trees were measured in replication, i.e., pairs/sets using same tree species (the C. fissilis cross-sections JAC09C and JAC21A, apart by at least 5 km), or different tree species (i.e., the C. fissilis cross-sections JAC09C and/or JAC21A, with the addition of a P. paniculata cross-section JAR26C located 0.2km from JAC09C, and ca. 5km from JAC21A).
The chronologies of both sites (BRA002 and BRA015) and the information of the tree-ring parameters of each sample can be accessed through the International Tree-Ring Data Bank (ITRDB, https://www.ncei.noaa.gov/products/paleoclimatology/tree-ring).
2.2 Sample handling for 14C analysis
For satisfactory results, especially regarding tree 14C reconstructions, dendrochronologically dated material sampled as single tree rings should contain the entire growing season. Thus, dendrochronology teams were advised to sample tree rings by cutting them along their delimited boundaries with great care without including neighboring rings. In Figure S3 (supplementary material) we illustrate suboptimal versus optimal tree-ring sampling for atmospheric 14C reconstructions. Researchers can determine ring boundaries due to anatomical features of growth ring formation in tropical hardwoods, which can include density variation of fiber cells, marginal parenchyma bands, alternating parenchyma and fiber bands, vessel cluster frequencies and vessel diameter sizes (Worbes Reference Worbes1989, Quintilhan et al. Reference Quintilhan, Santini, Rodriguez, Guillemot, Cesilio, Chambi-Legoas, Nouvellon and Tomazello-Filho2021, Ortega-Rodriguez et al. Reference Ortega-Rodriguez, Hevia, Sánchez-Salguero, Santini, de Carvalho, Roig and Tomazello-Filho2022). Two recommendations were further advised: i) to facilitate the isolation of individual samples from cross-sections, we emphasized first cutting a thin radial cross-section from it (see examples in Santos et al. Reference Santos, Rodriguez, Barreto, Assis-Pereira, Barbosa, Roig and Tomazello-Filho2021, Reference Santos, Granato-Souza, Ancapichún, Oelkers, Haines, De Pol-Holz, Andreu-Hayles, Hua and Barbosa2024; and Figure S3C); and ii) avoiding glue that is typically used to mount increment-cores for dendrochronology dating.
Once in the lab, before chemical extractions for 14C dating analysis took place, individual tree-ring samples were meticulously inspected under a microscope to ensure that they were not overcut or undercut, or if they needed to be re-shaped (Figure 3). Earlywood and latewood fractions of tree rings can have different 14C values (Haines et al. Reference Haines, Hiscock, Palmer, Turney, Thomas, Cadd, Vohra and Marjo2023; Kudsk et al. Reference Kudsk, Olsen, Nielsen, Fogtmann-Schulz, Knudsen and Karoff2018). Thus earlywood/latewood fractions should also be carefully evaluated for their ability to represent a full growth-season. Enlarged images of each tree species’ wood slabs or strips assisted in assessing the sampled tree-ring shapes. Undercut tree rings that potentially include a different year’s carbon, or those of irregular earlywood to latewood shapes, may be adjusted by using a transversal cutting technique. Overcut tree rings that showed missing portions of the earlywood or latewood fractions were replaced by samples that could represent the entire annual growing season.
Diagrammatic example of a misaligned tree ring received at the 14C-AMS laboratory (top right). Uneven tree rings with earlywood (EW) light-colored band somewhat shorter, than the latewood (LW). 14C-AMS laboratory and steps to trim the sample so the only potentially different dimension being sampled is the earlywood/latewood width. Proper transverse cutting will equalize the height and tangential dimensions of the ring block so only the ring widths would have different dimensions.
Here, all tree-ring samples were adequately evaluated for their quality regarding their earlywood and latewood fractions under the microscope. Later, they were weighed out for their total mass, then reduced to chips. Santos et al. (Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023) favor wood chip reduction by cutting radial slivers of wood from the tree-ring chips above glazed ceramic dishes of different sizes, as it allows for easy wood fraction recovery without contamination (refer to video in the supplementary material of Santos et al. Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023). Once the wood was reduced to slivers, the samples were weighed out again. This procedure aimed for high recovery yields. On average, we obtained a wood recovery yield of 99.6% (n = 30 independent wood samples) for the Juvenília site, and 99.9% (n = 89 independent wood samples) for Jamari National Forest.
Chemical extractions to alpha-cellulose were performed on an aliquot of tree-ring slivers weighing 10–30 mg, following the protocol of Santos et al. (Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023), which was adapted from a previous work of Southon and Magana (Reference Southon and Magana2010). The former chemical protocol produces pure alpha-cellulose extracts and well-mixed cellulose fibers. Chemical treatment can be summarized in the following steps. 1 N HCl and NaOH at 70°C cycling until the supernatant is clear. A bleaching step with 1 N HCl-NaClO2 at 70°C was employed for about 6 hr. The alpha-cellulose residual fraction was reached after 17.5% NaOH treatment for 1–2 hr at room temperature. The alpha-cellulose fraction was then treated with 1 N HCl at 70°C to remove atmospheric CO2 adsorbed during treatments and brought to a neutral pH. Cellulose fibers were homogenized in an ultra-sonic bath and dried for later processing.
For the purpose of 14C-AMS measurements, alpha-cellulose extracts were loaded into pre-baked quartz tubes with CuO, sealed after evacuation, and burned at 900°C to produce CO2. Later, CO2 samples were released in a vacuum line, cryogenically cleaned and transferred to pre-baked vessels loaded with Zn/TiH2 and Fe for the reduction of CO2 to filamentous graphite (Santos and Xu et al. Reference Santos and Xu2017). Graphite-Fe mixtures were pressed into target holders and underwent measurement by AMS, using an in-house modified NEC 0.5MV 1.5SDH-1 spectrometer equipped with online δ13C-AMS capabilities (Beverly et al. Reference Beverly, Beaumont, Tauz, Ormsby, von Reden, Santos and Southon2010). This instrument typically allows for precision better than 0.3%, based on measurements of 6-7 targets of Oxalic acid I (OX-I or NIST HOxI SRM 4990B), as the primary standard. Combustible reference materials, such as OX-II (oxalic acid II) and ANU (sucrose), were used to confirm instrumentation precision levels. For quality control of alpha-cellulose chemical processing and for background, the reference materials of FIRI-J (barley mash) and AVR-07-PAL-37 (wood blank) were used, respectively. Radiocarbon measurements were performed on graphite targets of 0.4 to 0.9mgC. Those are considered regular-sized targets on our in-house modified NEC 0.5MV 1.5SDH-1 spectrometer, as no significant changes in beam current and/or uncertainties (precision) were observed. The data are reported as a fraction of modern carbon (F14C), following the conventions shown in Stuiver and Polach (Reference Stuiver and Polach1977) and Reimer et al. (Reference Reimer, Brown and Reimer2004) when dealing with post-1950 AD samples.
3. Results and discussion
3.1 Assessment of 14C results
To assess potential sources of error (such as, 14C sample processing to graphite, measurement instrument bias errors, etc.) we measured combustible samples of OX-I, OX-II or ANU along with cellulose fibers. For alpha-cellulose analysis, we processed and measured duplicates of reference materials, such as FIRI-J barley. In all instances, precision was on the order of 0.25% or lower. Precision fitness from tropical cellulose fibers was evaluated by the weighted mean F14C and propagated errors between 14C results of pairs and/or sets of tree rings that shared the same calendar years of a given site. To assess alpha-cellulose homogeneity and tree-related variables influencing 14C values, we compared 14C results using several pair/set combinations. The legends in Figures 4 and 5 contain information about those combinations.
Data accuracy was appraised through the difference between the weighted mean F14C detected from tree rings for each site and the datasets associated with the Dec/Jan/Feb-TLPB-influenced regions, i.e., of SH1-2 or SH3 (Hua et al. Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and Palmer2022). Since atmospheric 14C compilations still lack data in SH lower and middle latitudes, we posit that small 14C differences are expected. For example, current SH1-2 and/or SH3 curves do not fully reflect latitudinal/longitudinal air mass circulation, terrestrial carbon isoflux effects, and/or local contributions of small nuclear tests after the 1963 Limited Test Ban Treaty around the world. In addition, Hua et al. (Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and Palmer2022) assume that vegetation over SH grows during Dec/Jan/Feb, hence precise fine-tuning cannot be expected (see discussions in sections 3.1.1 and 3.1.2).
To visually compare data, we made use of a 2-dimensional scatter plot where dendrochronological dates (or calendar years AD) were directly compared to 14C calibrated ages by CALIBomb software (http://calib.org/CALIBomb/). When corresponding dates from the two datasets were equal to each other, they fell exactly on the identity line 1:1 scatter plot. In order to convert the post-1950 AD F14C values obtained to 14C calibrated ages, we first calculated the weighted mean of pairs/sets and their respective propagated and standard error. We used the greatest of two given errors of the weighted mean F14C values as inputs in CALIBomb to produce ≈95% (2-sigma errors; P< 0.05). As expected, when converting 14C ages to 14C calibrated ages after 1950 AD, each weighted mean F14C and error resulted in two possible 14C calibrated age ranges. Thus, we selected the calibrated age ranges that were in closer proximity to the sequence of calendar ages expected for the dendrochronological dates.
To adjust the dendrochronological dates to the time when trees photosynthesized atmospheric 14CO2 and form tree-ring cellulose, we made use of the length of the growing season at each site, which is associated with the rainy season (Figure S1). Even though annual variations in vegetation growth may occur, small differences from the growing-season length should not be sufficient to affect our conclusions. A decimal date was then centered in the middle of the growing season, while the length of the growth season was considered as the error-bar. Thus, for the Juvenília site all calendar dates were adjusted to 0.002 ± 0.166 (November to March; 4 months in total; Figure S1A), while for Jamari we used 0.042 ± 0.291 (October to May; 7 months in total; Figure S1B).
In the supplementary material, we also provide the F14C associated value of each calendar year’s individual tree rings alongside the current atmospheric post-1950 AD 14C calibration curves for each site (see Figures S4 and S5), as well as the 14C data per dated tree (see Tables S2 and S3).
3.1.1 Juvenília samples 14C data precision/accuracy performance
In Figure 4 we show the 14C data precision/accuracy performance associated with the trees of the Juvenília site. A total of 15 replicate 14C analyses were produced based on 30 measurements from this site.
In Figure 4A, we compared the data reproducibility based on propagated errors calculated for each weighted mean F14C result of the tree rings pairs. For the 14C results/pairs that belong solely to the cross-section MTV001D and to the increment core MTV114B (cellulose homogeneity), and to those using both individual trees > 20km apart (Figure 1), we used the different symbols following the scheme shown in the figure 4A legend. Although our calculations are comparing just two data points per calendar year, the context and initial differences on how they were sampled, handled, and produced are enough to justify them. Here, 93% of replicated data (as cellulose fibers) yielded precision that is < 0.3% (n = 14 pairs/sets), while 67% are notably below 0.2% (n = 10 pairs/sets). Individual 14C results are shown in Table S2 in the supplementary material. A standard pooled deviation, using all measured pairs, yielded 0.23% (n = 15 pairs).
Radiocarbon data accuracy was evaluated by using the description in section 3.1, and it is displayed in Figure 4B. Due to the location of the Juvenília site, outside of the TLPB boundary on South America (Figure 1A), we chose the SH1-2 curve (Hua et al. Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and Palmer2022) to estimate the 14C calibrated ages based on the F14C values derived from the tree rings. Note that the fragmented nature of the 14C datasets in the SH1-2 led this compilation to be dominated by 14C signatures from sites at higher southern latitudes in Australia and New Zealand. Nonetheless, besides the calendar year of 1988, our data in figure 4B show an excellent agreement with 14C calibrated ages based on the SH1-2 curve (Hua et al. Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and Palmer2022) for all other Juvenília assigned calendar years (as homogenized alpha-cellulose fibers). The extremely high 14C data consistency detected here on post-1950 C. fissilis tree-ring samples shows that field wood sample collection techniques (cross-sections versus increment core) as well as tree spatial distribution are not relevant for data heterogeneity overall (Figure 4B, dashed line = 1:1 match, solid line = best-fit linear regression, R2 =0.998).
Radiocarbon data performance associated with trees of the Juvenília site. (A) Data reproducibility among pairs for different tree samples for calendar years ranging between 1964 and 1989. The legend/diagram describes the combinations tested. Solid arrows represent pairs within a tree and dashed arrows between trees. (B) Dendrochronological assigned calendar years (x-axis) were adjusted to 0.002 (middle of the growing season, as described in text) versus 14C-calibrated calendar years (y-axis) of averaged pairs/sets. The error for dendrochronological dates has been estimated based on the length of the growing season (i.e., ± 0.166; 4 months in total) and is hidden by symbols. The 14C-calibrated calendar year errors were calculated based on age ranges derived by CALIBomb. The black-dashed line shows the 1:1 ratio that is expected if all 14C-calibrated calendar years agree exactly with dendrochronological dates. Detailed results are provided in supplementary materials (Figure S4 and Table S2).
3.1.2 Jamari National Forest samples 14C data precision/accuracy performance
In Figure 5, we show the 14C data precision/accuracy performance associated with the trees of the Jamari National Forest site. A total of 36 replications were performed based on 89 individual 14C measurements.
Data reproducibility in Figure 5A represents the error associated with the weighted mean F14C result of the tree-rings pairs/sets corresponding to each calendar year between 1955 and 1990. Reproducibility tests were more complex at this site, i.e., we explored possible differences between same individual trees/same species, pair/sets involving all tree species, and pairs/sets using two distinct tree species (legend in Figure 5A). The consistency of the 14C data is rather similar among pairs/sets regardless of distances between trees (Figure 1), differences in tree species and/or pairs/sets combinations (Figure 2). Overall, most of the tree-ring calendar-year sets analyzed yielded a precision below 0.3% (83%), with the majority below 0.2% (64%). Reproducibility that fell between 0.3% and 0.8% (17%) are from pair/sets that involved P. paniculata and one or both of the C. fissilis samples (Figure 5A).
Despite the noticeable differences in wood anatomy, density, and chemistry profiles between P. paniculata and C. fissilis, both species can produce climatically sensitive growth rings, well-defined by parenchyma bands (Ortega-Rodriguez et al. Reference Ortega-Rodriguez, Hevia, Sánchez-Salguero, Santini, de Carvalho, Roig and Tomazello-Filho2022). Thus, possible errors due to the delimitation of the ring boundaries can be completely ruled out. As recovery yield during laboratory wood chipping was close to 100%, problems regarding this step can also be ruled out. Third, we assumed that growth of P. paniculata and C. fissilis tree species in the Jamari National Forest occurred during similar time frames, as mentioned earlier. Still, maximum variance was observed just for pairs/sets related to the calendar years of 1962, 1963, and 1965, but not 1964 (Table S3). Throughout all those years, atmospheric 14C changes during a single year were enormous, but that should include the 1964 calendar year as well. Thus, we suspect that small differences in tree life-cycles are not a major factor either. P. paniculata, also known as “purpleheart,” has a higher wood density than C. fissilis, i.e., 0.80 versus 0.55 g/cm3 (Zanne et al. Reference Zanne, Lopez-Gonzalez, Coomes, Ilic, Jansen, Lewis, Miller, Swenson, Wiemann and Chave2009). In addition, the wood of P. paniculata has a straight or wavy grain (sometimes interlocked) and is loaded with resins (Galletly Reference Galletly1889). Those characteristics make wood working with P. paniculata more challenging, including sawing, planning, turning, drilling and shaping (ITTO, 2024; Wood Database). Thus, we suspect this was the main reason for the small variabilities (0.3% < x < 0.8%) observed for the 14C pairs/sets of 1962, 1963 and 1965, but not to 1964. As mentioned before, isolation of whole single tree rings from denser and resinous woods by hand requires more skill, time and effort. Still, the 14C results of multiple tree species, when combined as propagated error, did not exceed 0.8%. A pooled standard deviation yielded 0.35% by computing all measured pairs/sets (n = 36 pairs/sets). Given the difficulty of working with denser and/or resinous woods, our results are extremely encouraging. As annually-resolved rings in quantities for dendrochronology dating are not always available at key sites of interest in the tropics, overlapping tree species can be a viable option for 14C reconstructions, when master chronologies are not feasible. However, sufficient 14C data replication will be obviously required to assure data quality, as we demonstrated here.
For this site, we also examined 14C data accuracy by using the estimated 14C calibrated ages averaging the F14C values derived from the tree-rings pairs/sets and CALIBomb. Since Jamari National Forest seems to be inside the TLPB boundary in South America (Figure 1A), we chose SH3 curve to calibrate our data for this location (Hua et al. Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and Palmer2022). We observed that all Jamari National Forest assigned calendar years are in good (if not excellent) agreement with expected 14C calibrated ages, regardless of the fact that 86% of these data make use of mixed tree species (Figure 5B, dashed line = 1:1 match, solid line = best-fit linear regression, R2 =0.998). Here, slight differences in 14C were mostly concentrated around the early 1960s and the 1970s. The first can be attributed to asymmetry among site longitudes and terrestrial carbon fluxes, as evidenced by Santos et al. (Reference Santos, Granato-Souza, Ancapichún, Oelkers, Haines, De Pol-Holz, Andreu-Hayles, Hua and Barbosa2024), when comparing 14C records of the Eastern Amazon Basin, Brazil (1°S, 53°W; 1940–2016) and Muna Island (5°S, 122°E; 1951-1979; the backbone of the SH3 record—Hua et al. Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and Palmer2022). Regarding the 1970s, after the 1963 Limited Test Ban Treaty, we attribute differences in 14C between our site and SH3 to China’s Lop Nor nuclear tests (42.4°N, 88.3°E; 1964–1980). These nuclear tests have been detected in the Muna Island and the Huangzhong district, China (36.3°N, 101.7°E; Xiong et al Reference Xiong, Zhou, Hou, Cheng, Du, Zhao, Wu, Hou, Lu and Fu2021) 14C records, but not at the Eastern Amazon Basin 14C record (Santos et al. Reference Santos, Granato-Souza, Ancapichún, Oelkers, Haines, De Pol-Holz, Andreu-Hayles, Hua and Barbosa2024). As the Jamari site is located in the Southwestern Brazilian Amazonia (Figure 1A), we suspect that this site would not capture Lop Nor nuclear tests as well (Figure S5). China’s Lop Nor nuclear tests were about 2-5 Mt TNT in magnitude (Xiong et al Reference Xiong, Zhou, Hou, Cheng, Du, Zhao, Wu, Hou, Lu and Fu2021), and therefore would have a small influence in the upper atmosphere. Nonetheless, high consistency 14C results obtained here on post-1950 tree-ring samples from distinct tree species (P. paniculata and C. fissilis) show that tree species diversity is not the primary factor influencing heterogeneity of results.
Radiocarbon data performance associated with the trees of the Jamari site. (A) Data reproducibility among sets of distinct tree samples of same or different tree species for calendar years between 1955 and 1990. The inset legend/diagram indicates the combinations tested. Open-square symbols compare reproducibility of two C. fissilis trees, open diamond symbols compare one C. fissilis and one P. paniculata tree, and solid diamond symbols compare all three trees. (B) Dendrochronologically determined calendar years (x-axis) were adjusted to 0.042 (middle of the growing season, as described in text) versus 14C calibrated calendar years (y-axis) of averaged pairs/sets. The error for dendrochronological dates has been estimated based on the length of the growing season (i.e., ± 0.291; 7 months in total), and are within the bounds of the symbol area. The 14C calibrated calendar years errors were calculated based on age ranges derived by CALIBomb. The black-dashed line shows the 1:1 ratio that is expected if all 14C calibrated calendar years agree exactly with dendrochronological dates. Detailed results are provided in supplementary materials (Figure S5 and Table S3).
3.2 Relevance of findings
Physical and chemical tree-ring preparation are essential prerequisites for precise and accurate results during 14C analysis. However, the many involved steps during sample handling and processing can add unwanted errors. Santos et al. (Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023) laid out steps on how to prevent and reduce errors when analyzing single trees for atmospheric 14C reconstructions and applied those techniques to produce tropical tree-ring 14C records using parenchyma-rich Hymenolobium petraeum and Cedrela odorata tree species from the Amazon Basin (Santos et al. Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023 and Reference Santos, Granato-Souza, Ancapichún, Oelkers, Haines, De Pol-Holz, Andreu-Hayles, Hua and Barbosa2024, respectively). However, questions remained regarding other tropical tree species and the possibility of combining several trees or tree species from a given region to derive fine-tuned atmospheric 14C information.
Here, we tested a number of Cedrela fissilis samples (i.e., multiple trees) from two sites in South America, from dry and wet forests. Trees from each dendrochronological collection were separated by 5 or 22.5 km and were from cross-sections or increment cores. Our 14C results also represent samples from species with distinctly different wood density, i.e., C. fissilis (0.55 g/cm3) and P. paniculata (0.80 g/cm3) (Zanne et al. Reference Zanne, Lopez-Gonzalez, Coomes, Ilic, Jansen, Lewis, Miller, Swenson, Wiemann and Chave2009). In sum, 15 pairs of 14C results from the Juvenília site and another 36 from the Jamari site, which include multiple tree species combinations are shown. In all instances, 14C data appear precise and accurate within error, and therefore external factors (such as microhabitat-related conditions affecting trees) do not appear to be relevant for bomb-induced 14C variability in the tropics.
Precision frequency distributions are shown in Figure 6 for the pair/sets presented in this study and those of duplicated tree 14C data using Araucaria aungustifolia (n = 9 pairs; Santos et al. Reference Santos, Linares, Lisi and Tomazello Filho2015), Pinus radiata and Agathis australis (n = 14 pairs; Turnbull et al. Reference Turnbull, Mikaloff Fletcher, Brailsford, Moss, Norris and Steinkamp2017), and Polylepis tarapacana (n = 21 pairs; Ancapichún et al. Reference Ancapichún, De Pol-Holz, Christie, Santos, Collado-Fabbri, Garreaud, Lambert, Orfanoz-Cheuquelaf, Rojas, Southon and Turnbull2021) that have been included in the global post-1950 map with replicated data (Hua et al. Reference Hua, Turnbull, Santos, Rakowski, Ancapichún, De Pol-Holz, Hammer, Lehman, Levin, Miller and Palmer2022). Except for the later work (i.e., Ancapichún et al. Reference Ancapichún, De Pol-Holz, Christie, Santos, Collado-Fabbri, Garreaud, Lambert, Orfanoz-Cheuquelaf, Rojas, Southon and Turnbull2021, which used more than one tree of the same species in their replication analysis), the former two studies reported 14C results of pairs of the same trees. Note that we chose not to add the Campbell Island tree 14C record of Turney et al. (Reference Turney, Palmer, Maslin, Hogg, Fogwill, Southon, Fenwick, Helle, Wilmshurst, McGlone, Bronk Ramsey, Thomas, Lipson, Beaven, Jones, Andrews and Hua2018) in this evaluation. Due to low tree availability across this region, this record was constructed with multiple distinct single trees that bear large calendar date errors. Thus, some inherent variability in the Campbell Island tree 14C record is expected.
Poor performance in tree 14C data is sometimes attributed to chemical procedures adopted by laboratories, earlywood nonstructural carbon transfer or other confounding factors that are difficult to demonstrate. Here, we evaluate tree 14C data fitness using a wider range of factors, including dry and wet forests, individual tree distances, and distinct tree species together with wood’s variable properties (density). Regardless of the combinations assessed (Figures 4 and 5), our data show a large proportion of 14C results with minimum variance (close to 0.3% or lower). This outcome is similar or even better than that of other studies (Figure 6). We can conclude that careful wood and cellulose processing (prior 14C measurements) played important roles in precision of data.
Frequency distribution versus 14C precision of pair/sets from several sites and studies. Number of replicated pairs/sets were variable, i.e., 14 pairs from Bearing Head (Turnbull et al. Reference Turnbull, Mikaloff Fletcher, Brailsford, Moss, Norris and Steinkamp2017), 9 pairs from Camanducaia (Santos et al. Reference Santos, Linares, Lisi and Tomazello Filho2015), 21 pairs from Altiplano (Ancapichún et al. Reference Ancapichún, De Pol-Holz, Christie, Santos, Collado-Fabbri, Garreaud, Lambert, Orfanoz-Cheuquelaf, Rojas, Southon and Turnbull2021), and 51 pairs/sets from this study (15 pairs from Juvenília and 36 from Jamari sites). Precision bins are shown as percentages (%).
The findings reported here are significant, as by improving 14C data consistency among several trees/sites, researchers can fill several spatial and temporal gaps and/or improve atmospheric 14C data in a more comprehensive manner. Moreover, even though the Amazon Basin hosts over 16,000 tree species (Ter Steege et al. Reference Ter Steege, Pitman, Killeen, Laurance, Peres, Guevara, Salomão, Castilho, Amaral, de Almeida Matos and de Souza Coelho2015), not all those tree species show annual growth rings. Plus, those that show annual growth may not occur within the same micro-environment or may not be long-living species (Linares et al. Reference Linares, Santos, Brandes, Barros, Lisi, Balieiro and de Faria2017). Thus, combining tree species for a given site to derive atmospheric 14C information may be required. Here, we demonstrated that this is feasible when we overlapped C. fissilis and P. Paniculata tree species. In addition, mapping present-day 14CO2 distribution in trees in places where air-sampling stations are mostly scarce (such as the tropics) is relevant to better evaluating future changes, e.g., how atmospheric circulation may change as the climate warms (TLPB expansion or displacement), and how trees will cope under extreme climate conditions (physiological responses). While our efforts to improve 14C data quality focus on tropical trees in general during the post-1950 AD period, we believe that the steps adopted here can also be used on extratropical tree species, other 14C age-ranges, and applications.
Achieving robust 14C data consistency can be then attributed to the coordinated procedures of careful handling of whole tree rings during physical individual sampling from wood laths and laboratory steps, a straightforward chemical extraction procedure (Santos et al. Reference Santos, Komatsu, Renteria, Brandes, Leong, Collado-Fabbri and De Pol-Holz2023, Griffin et al. Reference Griffin, Santos, Nguyen, Rodriguez, Pereira, Jaén-Barrios, Assis-Pereira, de Oliveira Barreto, Brandes, Barbosa and Groenendijk2024), as well as homogenization of cellulose once chemical extractions reach completion. The latter is useful when cellulose extracts per individual tree ring are significantly large, as just a fraction of it will be used for further 14C sample processing.
4. Conclusions and further recommendations
A fundamental limitation to any scientific question is normally associated with the overall data uncertainty. In tree-ring 14C analysis, which involves considerable sample handling and processing, data uncertainty can sometimes become overwhelming. Hence, important questions can no longer be answered. While the overall data uncertainty can never be reduced to zero due to analytical errors inherent in 14C-AMS measurements, the implementation of careful methodologies while handling tree-ring samples can help reduce sample processing errors.
In this work, we reported post-1950 AD 14C replicated results from tropical tree rings sampled from dendrochronological collections containing Peltogyne paniculata and multiple Cedrela fissilis tree species. The tree rings tested showed different anatomical and chemical characteristics and were collected from two separate sites on South America under different precipitation regimes. Moreover, replicated data involved trees spaced apart by significant distances and/or combined different tree species. Regardless of the criteria differences among distinctive types of trees/sites tested, variance and accuracy among 14C results were fully reproducible. Thus, the most important criteria for 14C data consistency were credited to careful tree-ring sampling and handling at both laboratories (tree ring and 14C), followed by thorough homogenization of cellulose fibers before 14C sample processing and analysis took place. By following those strategies, whole rings representing the intake 14C for a full growing season were properly detected and reproduced within uncertainties.
For most cases, the sum of all variances in the final tree-ring 14C measurements were 0.3% or lower (88% of the data) regardless of the distances between trees, especially when the same tree species with similar density levels were used. The overall high-performance for our post-1950 AD results, which involved a comprehensive number of tree-ring 14C measurements (n = 119 individual measurements in total), allowed us to obtain a more accurate picture of the 14C concentration of those tree rings. Our results suggest that several trees from the tropics with continuous annual growth patterns that can be discerned by dendrochronology methods, can also be incorporated into fine-tuned atmospheric 14C reconstructions, when variance-reduction strategies (as those reported here) are applied. Further recommendations for future 14C datasets would include sufficient data replication per 14C records produced, for a better understanding of data and the factors affecting them.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/RDC.2024.105
Acknowledgments
G.M.S. thanks the support of the United States National Science Foundation (Grant# AGS-1903690). N.O.B., D.R.O.R, and G.A. thanks for the partial support provided by FAPESP (Grants # 2023/14668-5; 2023/07753-6; 2020/01378-0; 2019/26350-4; 2018/22914-8; 2017/50085-3; 2009/53951-7). A.C.B. was supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq: PQ 313129/2022-3) and Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG: APQ-01544-22). We thank the anonymous reviewers for their careful reading of our manuscript and their many insightful comments and suggestions, as well as the journal editor Dr. Jull. A special thanks to Dr. Leavitt for meaningful suggestions to our work that have boosted clarity.
Open access
