Introduction
The ability to rapidly measure the 14C content of sub-milligram quantities of biogenic carbonate without conversion to graphite has risen in popularity over the past decade (e.g., Lougheed et al. Reference Lougheed, Snowball, Moros, Kabel, Muscheler, Virtasalo and Wacker2012, Reference Lougheed, Metcalfe, Ninnemann and Wacker2018; Bush et al. Reference Bush, Santos, Xu, Southon, Thiagarajan, Hines and Adkins2013; Longworth et al. Reference Longworth, Robinson, Roberts, Beaupre, Burke and Jenkins2013; Wacker et al. Reference Wacker, Fahrni, Hajdas, Molnar, Synal, Szidat and Zhang2013a, Reference Wacker, Lippold, Molnár and Schulz2013b; Dominguez et al. Reference Dominguez, Kosnik, Allen, Hua, Jacob, Kaufman and Whitacre2016; Kosnik et al. Reference Kosnik, Hua, Kaufman, Kowalewski and Whitacre2017; Ritter et al. Reference Ritter, Erthal, Kosnik, Ciombra and Kaufman2017, Reference Ritter, Erthal, Kosnik, Kowalewski, Coimbra, Caron and Kaufman2023; Gottschalk et al. Reference Gottschalk, Szidat, Michel, Mazaud, Salazar, Battaglia, Lippold and Jaccard2018; Kowalewski et al. Reference Kowalewski, Casebolt, Hua, Whitacre, Kaufman and Kosnik2018; Tuna et al. Reference Tuna, Fagault, Bonvalot, Capano and Bard2018; Fagault et al. Reference Fagault, Tuna, Rostek and Bard2019; Lindauer et al. Reference Lindauer, Friedrich, van Gyseghem, Schöne and Hinderer2019; New et al. Reference New, Yanes, Cameron, Miller, Teixeira and Kaufman2019; Parker et al. Reference Parker, Yanes, Hernández, Hernández Marreno, Paris and Surge2019; Albano et al. Reference Albano, Hua, Kaufman, Tomašových, Zuschin and Agadi2020, Reference Albano, Sabbatini, Lattanzio, Päßler, Steger, Hua, Kaufman, Szidat, Zuschin and Negri2023a, Reference Albano, Hua, Kaufman and Zuschin2023b; Missiaen et al. Reference Missiaen, Wacker, Lougheed, Skinner, Hajdas, Nouet, Pichat and Waelbroeck2020; Dolman et al. Reference Dolman, Groeneveld, Mollenhauer, Ho and Laepple2021; Mollenhauer et al. Reference Mollenhauer, Grotheer, Gentz, Bonk and Hefter2021; Nawrot et al. Reference Nawrot, Berensmeier, Gallmetzer, Haselmair, Tomašových and Zuschin2022; Sanchez et al. Reference Sanchez, Yanes, Linstädter and Hutterer2022; Steger et al. Reference Steger, Bosnjak, Belmaker, Galil, Zuschin and Albano2022). Dating small carbonate samples has become increasingly popular in time averaging or population structure studies (e.g., Kowalewski et al. Reference Kowalewski, Casebolt, Hua, Whitacre, Kaufman and Kosnik2018; Nawrot et al. Reference Nawrot, Berensmeier, Gallmetzer, Haselmair, Tomašových and Zuschin2022; Ritter et al. Reference Ritter, Erthal, Kosnik, Kowalewski, Coimbra, Caron and Kaufman2023) and in the field of marine geochronology where a few dozen or even single large foraminifera can be dated (e.g., Lougheed et al. Reference Lougheed, Metcalfe, Ninnemann and Wacker2018; Mollenhauer et al. Reference Mollenhauer, Grotheer, Gentz, Bonk and Hefter2021). Two rapid 14C methods are currently employed: direct carbonate and gas ion source (GIS). The direct carbonate method produces carbon ions (C-) by directly sputtering powdered carbonate and metal binder (typically iron or niobium) mixtures with cesium atoms (Bush et al. Reference Bush, Santos, Xu, Southon, Thiagarajan, Hines and Adkins2013; Longworth et al. Reference Longworth, Robinson, Roberts, Beaupre, Burke and Jenkins2013; Hua et al. Reference Hua, Levchenko and Kosnik2019). Alternatively, C- ions are produced by the GIS method through the acid dissolution of carbonate and then sputtering the released CO2 gas with cesium atoms in the presence of titanium (Middleton Reference Middleton1984; Bronk and Hedges 1987). The GIS method has become increasingly popular with the development of the MIni CArbon DAting System (MICADAS) (Synal et al. Reference Synal, Stocker and Suter2007) coupled with a carbonate handling system (CHS) inlet to the GIS (Ruff et al. Reference Ruff, Wacker, Gäggeler, Suter, Synal and Szidat2007, Reference Ruff, Fahrni, Gäggeler, Hajdas, Suter, Synal, Szidat and Wacker2010; Wacker et al. Reference Wacker, Fahrni, Hajdas, Molnar, Synal, Szidat and Zhang2013a; Lindauer et al. Reference Lindauer, Friedrich, van Gyseghem, Schöne and Hinderer2019).
The Arizona Climate and Ecosystems (ACE) Isotope Laboratory at Northern Arizona University (NAU) brought a MICADAS system online in June 2021 (Ebert et al. Reference Ebert, Schuur, Kaufman, Brown, Propster, Kelley, Bright, Carbone, McKay and Koch2022). The ACE Lab includes elemental analyzers, a GIS, a CHS, and both automated graphitization equipment (AGE) and a manual cryogenic purification line.
This report compares the percent Modern Carbon (pMC) values and analytical precision from 48 samples of biogenic carbonate each analyzed as graphite, GIS, and as direct carbonate. The primary goal of our study is to compare the accuracy and reproducibility of carbonate pMC values from NAU’s MICADAS and GIS system against standard graphite and direct carbonate analyses.
Materials and Methods
This study features predominantly mollusk shells and fewer echinoid tests. Several of the standard graphite and direct carbonate analyses were included in an earlier study published prior to the arrival of the MICADAS at NAU (Bright et al. Reference Bright, Ebert, Kosnik, Southon, Whitacre, Albano, Flores, Frazer, Hua, Kowalewski, Martinelli, Oakley, Parker, Retelle, Ritter, Rivadeneira, Scarponi, Yanes, Zuschin and Kaufman2021). The compilation featured here comprises pMC values ranging from approximately 99 to 0.4 (radiocarbon ages approximately 0.1 to 45.1 ka BP), based on prior AMS analysis of samples archived at NAU (Bright et al. Reference Bright, Ebert, Kosnik, Southon, Whitacre, Albano, Flores, Frazer, Hua, Kowalewski, Martinelli, Oakley, Parker, Retelle, Ritter, Rivadeneira, Scarponi, Yanes, Zuschin and Kaufman2021). To fill in gaps in this dataset, new graphite and direct carbonate powders were processed and analyzed at NAU or at the W.M. Keck Carbon Cycle Accelerator Mass Spectrometer Facility at the University of California – Irvine (UCI), respectively. To obtain equivalent subsamples for the three analyses, each mollusk shell and echinoid test was sampled parallel to growth bands and/or sub-adjacent to each other. This approach minimizes the risk of sampling different aged shell material. For example, some mollusks are slow growing and can live for several hundred years (Moss et al. Reference Moss, Ivany, Judd, Cummings, Bearden, Kim, Artruc and Driscol2016).
All samples were sonicated briefly to remove loose debris before being exposed to ACS grade 2N hydrochloric acid which removed 30% of their mass. The leaching process was considered complete when effervescence stopped. All samples were rinsed three times with deionized reverse osmosis water (16.7 Mohm*cm) before being dried at 50°C in a convection oven.
Carbonate samples (0.5–1.5 mg; 60–180 μg C) destined for GIS 14C analysis were placed in baked (3 hr at 500° C) glass reaction vials sealed with screw-top caps containing a rubber septum (Exetainer #VW101). The vials were flushed with N2 gas and sealed until analyzed. The CHS automatically flushes each vial with helium using a double-walled needle, after which 85% phosphoric acid is injected to dissolve the shell material. The evolved CO2 is carried by helium through a water trap and then to the zeolite trap of the GIS. The CO2 is released from the zeolite at 450°C and transferred to a syringe where it is mixed with helium to a final CO2 concentration of 5% by volume, after which it is injected into the ion source. The zeolite trap is optimized for 80–100 μg C (Ruff et al. Reference Ruff, Fahrni, Gäggeler, Hajdas, Suter, Synal, Szidat and Wacker2010), which is equivalent to about 0.6–0.8 mg of calcium carbonate. Excess carbon is not analyzed. Detailed description of the MICADAS GIS and CHS analytical method is provided in Synal et al. (Reference Synal, Stocker and Suter2007) and Ruff et al. (Reference Ruff, Wacker, Gäggeler, Suter, Synal and Szidat2007, Reference Ruff, Fahrni, Gäggeler, Hajdas, Suter, Synal, Szidat and Wacker2010).
Ten new samples (0.3–0.5 mg; 36–60 μg C) for direct carbonate AMS analyses at UCI were manually ground to a fine powder with an agate mortar and pestle. The powdered shell was mixed with 6 to 7 mg of niobium powder (Alfa Aesar Puratronic, –325 mesh, 99.99%) in baked (3 hr at 500°C) Kimble borosilicate glass culture tubes (6 mm OD × 50 mm), flushed with nitrogen gas, and capped with Supelco plastic column caps (1/4” OD) until the carbonate-niobium mixture was manually pressed into pre-drilled (4.1 mm depth) aluminum targets before being sent to UCI.
Thirteen new samples for graphite analysis (8-10 mg; 960–1200 μg C) at NAU were placed in baked (3 hr at 500°C) acid-washed glass vials sealed with rubber septa from BD Vacutainer® plastic collection tubes (No. 366704). Ambient atmosphere was removed via vacuum before a small-bore needle added approximately 0.8 mL of ACS grade 85% phosphoric acid to each vial. Vials were warmed in a heating block set to 70°C until the shell material entirely dissolved. The evolved gas was removed via vacuum on a manual cryogenic purification line. Water vapor was removed by passing the gas through a mixture of ethanol and liquid nitrogen at approximately –80°C. Carbon dioxide was condensed to a solid using a liquid nitrogen bath and the remaining gases were drawn off. The purified CO2 was converted to graphite by reaction with iron powder (Alfa Aesar, -325 mesh, reduced 98%) in a hydrogen reducing environment at 550°C for 3 hr (Vogel et al. Reference Vogel, Southon, Nelson and Brown1984). The graphite-iron mixture was packed into targets using an automated press before being loaded onto the MICADAS.
Radiocarbon (14C) concentrations are reported as pMC following the conventions of Stuvier and Polach (Reference Stuvier and Polach1977). Sample preparation backgrounds have been subtracted based on measurements of 14C-free calcite processed in the same fashion as the unknowns. All 14C determinations have been corrected for isotopic fractionation according to the conventions of Stuvier and Polach (Reference Stuvier and Polach1977) with δ13C values being measured on the AMS. These δ13C values are superior when used to correct for fractionation but can differ from the actual value of the original material and are thus not reported. Machine performance was monitored by repeat analysis of the IAEA standards C1 (14C-dead marble; consensus pMC = 0.00 ± 0.02; Rozanski et al. Reference Rozanski, Stichler, Gofiantini, Scott, Buekens, Kromer and van der Plicht1992) and C2 (travertine; consensus pMC = 41.14 ± 0.03; Rozanski et al. Reference Rozanski, Stichler, Gofiantini, Scott, Buekens, Kromer and van der Plicht1992), the CAHI coral standard provided by UCI (pMC = 94.44 ± 0.19; Mollenhauer et al. Reference Mollenhauer, Grotheer, Gentz, Bonk and Hefter2021), and a Pliocene 14C-dead Tridacna shell provided by the Florida Museum of Natural History (specimen UF 143174) which accounts for possible matrix effects that the marble C1 blank might not. All GIS 14C determinations of unknowns have been calibrated by assigning a blank pMC value of 0 and forcing the IAEA C2 travertine standard and the CAHI coral standard runs in the same batch through their consensus values (e.g., Gottschalk et al. Reference Gottschalk, Szidat, Michel, Mazaud, Salazar, Battaglia, Lippold and Jaccard2018). The blank and standard compilations were analyzed using boxplot and whisker diagrams to identify outliers, which were removed, without iteration, from the mean and standard deviation calculations. The number of outliers is included in Table 1. Data from the MICADAS was processed using BATS software (version 4.0; Wacker et al. Reference Wacker, Bonani, Friedrich, Hajdas, Kromer, Nemec, Ruff, Suter, Synal and Vockenhuber2010). When reported, radiocarbon ages are uncalibrated years BP, without consideration of the marine reservoir effect, with BP meaning conventional radiocarbon years before AD 1950 (Stuvier and Polach Reference Stuvier and Polach1977).
N.A. – not applicable.
* Coral standard provided by UC-Irvine. Graphite reference value as cited in Mollenhauer et al. (Reference Mollenhauer, Grotheer, Gentz, Bonk and Hefter2021) from 294 analyses performed at the Keck Carbon Cycle Radiocarbon Laboratory at UC-Irvine. GIS “reference” value is based on 47 analyses reported in Table 3 of Mollenhauer et al. (Reference Mollenhauer, Grotheer, Gentz, Bonk and Hefter2021).
** Calibrated values.
The relationship between GIS and graphite or direct carbonate and graphite pMC values was evaluated using a reduced major axis regression (RMA) analysis, which minimizes the residual variation across both the X- and Y-axes (Quinn and Keough Reference Quinn and Keough2002; Smith Reference Smith2009). An RMA regression avoids assumptions about the dependent and independent variables between GIS and graphite or direct carbonate and graphite pMC values (Smith Reference Smith2009). The PAST 4.13 statistical program (Hammer et al. Reference Hammer, Harper and Ryan2001) was used for the RMA with 95% bootstrapped confidence intervals [n = 1999 replications]. The Akaike Information Criterion (AIC) (Akaike Reference Akaike1973; Cavanaugh and Neath Reference Cavanaugh and Neath2019) using the R language version 4.3.1 (R Core Team 2023) and R package “rcompanion” version 2.4.30 (Mangiafico Reference Mangiafico2023) was used to compare the performance of the GIS versus graphite and direct carbonate versus graphite models. Scores were calculated by AIC = 2K – 2ln(L), where K is the number of model parameters and ln(L) is the log-likelihood of the model. The lowest AIC score is considered the better fit.
Results and Discussion
ACE Laboratory Blank (IAEA C1 and Pliocene Tridacna Shell) and Holocene Standard (IAEA C2 and CAHI) Performance
The quality of our carbonate 14C determinations is monitored by repeat analysis of certified and internal reference blanks and standards (Table 1). The pMC values overlap published values at one standard deviation (Rozanski et al. Reference Rozanski, Stichler, Gofiantini, Scott, Buekens, Kromer and van der Plicht1992; Mollenhauer et al. Reference Mollenhauer, Grotheer, Gentz, Bonk and Hefter2021). At the ACE Lab, the C1 procedural blank yields graphite and GIS ages of 50.6 ± 1.9 and 40.2 ± 3.9 ka BP, respectively. The Tridacna procedural blank yields graphite and GIS ages of 49.2 ± 2.4 and 42.2 ± 4.3 ka BP, respectively.
GIS versus Graphite and Direct Carbonate versus Graphite FMC Determinations
There is considerable time savings when using the GIS method. Processing samples for the direct carbonate method requires manual powdering of the samples, followed by manual subsampling, weighing, and manipulation of that powder plus a metal binder into small borosilicate tubes. The powder and metal mixture is then carefully poured into manually drilled aluminum cathodes before being manually pressed into targets. In contrast, the GIS preparation method is more streamlined and only requires weighing the cleaned carbonate samples into glass reaction vials. Processing a shell for GIS analysis takes roughly half as much time as processing a shell for direct carbonate analysis. For example, an undergraduate worker typically produces a batch of 51 unknowns for direct carbonate analysis in approximately 12 hours. That same worker can produce a batch of 48 unknowns for GIS analysis in approximately five to six hours. Converting carbonate to graphite is more time consuming. It typically takes approximately three hours to process eight unknowns to where the purified CO2 is converted to graphite.
Forty-one blank-corrected pMC values range from 98.68 to 20.51 (Supple. Info.), or spanning 14C ages of approximately 0.1 to 12.6 ka BP. One mid-Pleistocene Rangia shell that yields measurable carbon was analyzed in triplicate using each of the three AMS methods. The triplicates were averaged into a single value for each method to provide a lower limit to the compilation at a pMC value of approximately 0.36 (Supple. Info.), or approximately 45.1 ka BP. Six additional 14C-dead mollusk and echinoid samples were analyzed using each AMS method to determine and compare the limits of detection (Supple. Info.).
The RMA regression of the graphite versus GIS pMC values yields a slope with a 95% confidence interval that includes 1.000 (Figure 1; Table 2), which outperforms the RMA regression of the graphite versus direct carbonate values (Figure 1; Table 2). The RMA regression of graphite versus GIS and graphite versus direct carbonate yield slopes that are statistically indistinguishable from each other (p-value same slope = 0.003). Analyzing the graphite versus GIS and graphite versus direct carbonate data as linear regressions using AIC, where graphite pMC values are the dependent variable, yields a lower AIC score for graphite versus GIS (102.5) than for graphite versus direct carbonate (147.2). The lower AIC score indicates that graphite values are better explained by a model in which GIS values are the predictor than a model in which direct carbonate values are the predictor (Akaike Reference Akaike1973; Cavanaugh and Neath Reference Cavanaugh and Neath2019).
* This study.
** Uncalibrated. Not normalized using a blank pMC value of 0.00 and IAEA C2 and UCI CAHI consensus values.
‡ Missiaen et al. (Reference Missiaen, Wacker, Lougheed, Skinner, Hajdas, Nouet, Pichat and Waelbroeck2020). Radiocarbon ages converted to pMC using the equation pMC = (e14C yr/−8033) × 100.
§ Wacker et al. (Reference Wacker, Lippold, Molnár and Schulz2013b).
An alternative method of data reduction that does not calibrate the unknown pMC values by forcing the IAEA C2 and CAHI standards in each batch through their consensus values does not appreciably change the outcome of the RMA comparison (Table 2). We conclude that GIS outperforms direct carbonate 14C determinations, acknowledging the relatively small number of paired analyses (n = 42) used in this comparison.
A larger study of 150 graphite versus direct carbonate pairs of mollusk, echinoid, and brachiopod carbonate yielded results that are more comparable to the graphite versus GIS comparison presented here (Table 2) (Bright et al. Reference Bright, Ebert, Kosnik, Southon, Whitacre, Albano, Flores, Frazer, Hua, Kowalewski, Martinelli, Oakley, Parker, Retelle, Ritter, Rivadeneira, Scarponi, Yanes, Zuschin and Kaufman2021). The RMA regression slope of graphite versus direct carbonate data in Bright et al. (Reference Bright, Ebert, Kosnik, Southon, Whitacre, Albano, Flores, Frazer, Hua, Kowalewski, Martinelli, Oakley, Parker, Retelle, Ritter, Rivadeneira, Scarponi, Yanes, Zuschin and Kaufman2021) and the graphite versus GIS regression slope in this study are not statistically distinguishable (p-value = 0.167). Additional comparisons using foraminifera have shown that graphite and GIS analyses also yield comparable results (Table 2) (Wacker et al. Reference Wacker, Lippold, Molnár and Schulz2013b; Gottschalk et al. Reference Gottschalk, Szidat, Michel, Mazaud, Salazar, Battaglia, Lippold and Jaccard2018; Missiaen et al. Reference Missiaen, Wacker, Lougheed, Skinner, Hajdas, Nouet, Pichat and Waelbroeck2020). Our mollusk and echinoid graphite versus GIS comparison yields slightly better results than the foraminifera comparisons, which we suspect reflects true age heterogeneity within the collections of analyzed foraminifera tests, which can be prone to mixing and reworking in some environments (Fagault et al. Reference Fagault, Tuna, Rostek and Bard2019; Dolman et al. Reference Dolman, Groeneveld, Mollenhauer, Ho and Laepple2021). Mollusk shells and echinoid tests have the advantage of being large enough to routinely date single specimens (e.g., Scarponi et al. Reference Scarponi, Kaufman, Amorosi and Kowalewski2013; Harnik et al. Reference Harnik, Torstenson and Williams2017; Kowalewski et al. Reference Kowalewski, Casebolt, Hua, Whitacre, Kaufman and Kosnik2018; Nawrot et al. Reference Nawrot, Berensmeier, Gallmetzer, Haselmair, Tomašových and Zuschin2022; Ritter et al. Reference Ritter, Erthal, Kosnik, Kowalewski, Coimbra, Caron and Kaufman2023) or individual growth lines (e.g., Lindauer et al. Reference Lindauer, Friedrich, van Gyseghem, Schöne and Hinderer2019; Towers Reference Towers2022).
The accuracy of the GIS 14C determinations outperforms the direct carbonate method when compared to graphite. Sixty percent (25/42) of the GIS determinations are within ±0.5 pMC of their graphite counterparts, compared to 26% (11/42) of the direct carbonate determinations (Figure 1). Furthermore, 86% (36/42) of the GIS determinations are within ±1.0 pMC of their graphite counterparts, compared to 67% (28/42) of the direct carbonate determinations (Figure 1). A larger study found that 39% (59/150) and 77% (116/150) of direct carbonate determinations were within ±0.5 pMC and ±1.0 pMC, respectively, of their graphite counterparts (Bright et al. Reference Bright, Ebert, Kosnik, Southon, Whitacre, Albano, Flores, Frazer, Hua, Kowalewski, Martinelli, Oakley, Parker, Retelle, Ritter, Rivadeneira, Scarponi, Yanes, Zuschin and Kaufman2021), indicating that GIS accuracy outperforms the direct carbonate method regardless of the number of samples analyzed. We suspect that part of the difference in performance is due to the smaller carbon content of the direct carbonate samples (36–60 μg C) compared to the GIS samples featured in this study (60–100 μg C), which may render the direct carbonate 14C determinations more vulnerable to any contamination introduced during the powdering process or with the addition of the metal binder. In addition, the direct carbonate 14C determinations were conducted over a roughly 5 yr span whereas most (38/42) of the GIS 14C determinations were run in a single day (Supple. Info.), thus, direct carbonate 14C measurements might have captured additional long-term variability in AMS performance. However, the average graphite and GIS pMC values over roughly two years of analysis for the IAEA C2 standard and the CAHI coral standard differ by 0.55 and 0.14 pMC, respectively (Table 1), suggesting the typically small differences between graphite and GIS pMC values measured in this study are robust and are not necessarily a function of the short analytical timeframe for the GIS samples.
Graphite 14C determinations have comparatively small, typically decadal-scale precision due in part to the long sputter times (approximately 75 minutes) and large number of carbon counts per analysis, which can be on the order of one million counts. In contrast, the GIS method analyzes a sample in roughly 15 minutes and generates a few thousand to a few tens of thousands of carbon counts, thus the lower precision is based on counting statistics. In our compilation, the GIS ages have precision on the order of ±60 to ±80 14C yrs back to about 6.5 ka BP. At around 9 ka BP, the precision decreases to about ±110 14C yrs, and decreases further to about ±130 14C yrs at 12.5 ka BP (Figure 2). The precision of the GIS is on par with that of the direct carbonate method and averages four to five times worse than graphite precision (Figure 2). Although the precision for the GIS and direct carbonate method both increase with age (Figure 2), the variation in the GIS precision (stdev = 17 14C years) is similar to that of graphite (stdev = 13 14C years), and both are better than for direct carbonate errors (stdev = 30 14C years). This is likely because the GIS and graphite methods spall C- ions from comparatively pure CO2 gas and graphite sources, respectively, which promotes more homogenous ion formation conditions (e.g., Middleton Reference Middleton1984; Bronk and Hedges 1987), whereas the direct carbonate method spalls C- ions from a mixture of powdered calcium carbonate and metal binder.
Six Plio-Pleistocene carbonate samples, determined to be 14C-dead by graphite analysis, yielded finite pMC values and 29.4 to 36.7 ka BP ages when analyzed using the direct carbonate method but yielded non-finite pMC values when analyzed by GIS (Supple. Info). Bush et al. (Reference Bush, Santos, Xu, Southon, Thiagarajan, Hines and Adkins2013) also reported a direct carbonate result from an old coral sample that was approximately 6 14C kyr younger than the paired graphite result. They attributed the difference to a variety of issues, including extended storage (five years) of their carbonate materials between the graphite and direct carbonate analyses, modern contamination of the carbonate powder during processing, and heterogeneity in the carbonate sample. They concluded that lower beam currents and lower count rates, coupled with modern contamination of the carbonate powders during processing led to less favorable results on samples with ages over 30 ka BP. In our study, there is a six-year difference between the direct carbonate analyses of the six 14C-dead samples and their graphite counterparts, whereas the GIS counterparts were processed five months after the graphite analyses. Thus, the more similar GIS and graphite pMC results on the 14C-dead shells could also be a function of reduced sampling time between the two analyses, as proposed by Bush et al. (Reference Bush, Santos, Xu, Southon, Thiagarajan, Hines and Adkins2013). The IAEA C1 blank yielded a pMC value of 0.67 ± 0.33 (n = 39) when run by GIS at ACE, whereas 188 direct carbonate analyses of IAEA C1 processed at ACE and analyzed at UCI between 2017 and 2021 yielded a pMC value of 1.67 ± 0.69. We suspect that the poorer performance of the direct carbonate analyses results from carbon contamination during powdering and the addition of a metal binder that is known to contain some carbon (Bush et al. Reference Bush, Santos, Xu, Southon, Thiagarajan, Hines and Adkins2013; Hua et al. Reference Hua, Levchenko and Kosnik2019). Similarly, 142 direct carbonate analyses of the IAEA C2 standard processed at ACE and analyzed at UCI between 2017 and 2021 yielded a pMC value of 40.81 ± 0.63, whereas 89 GIS analyses of C2 analyzed at ACE over the past two years yielded a pMC value of 41.16 ± 0.64 (Table 2). The two mean values are statistically distinguishable at a 95% confidence interval (p-value = 6 x 10−5) and the GIS average falls closer to the consensus value of 41.14 ± 0.03 (Rozanski et al. Reference Rozanski, Stichler, Gofiantini, Scott, Buekens, Kromer and van der Plicht1992).
Thus, collectively, we conclude that the GIS method is superior to the direct carbonate method, especially when it comes to older samples (e.g., > 20 ka BP).
Conclusions
This study compared 42 samples of biogenic carbonate (mollusk shells and echinoid tests) that were dated using a MICADAS plus gas interface system (GIS) at Northern Arizona University’s Arizona Climate and Ecosystems Lab, by direct carbonate methods, and by standard graphite.
-
Preparing carbonate samples for GIS analysis is roughly 50% less time-consuming than required for the direct carbonate method.
-
The GIS method yields pMC values that are virtually indistinguishable from graphite (slope = 1.011 ± 0.006, 95% CI [0.997–1.023]; R2 = 0.999) and outperforms the direct carbonate method (slope = 0.978 ± 0.009, 95% CI [0.959 – 0.999]; R2 = 0.997).
-
Sixty percent (25/42) and 86% (36/42) of the GIS determinations are within ±0.5 and ±1.0 pMC of their graphite counterparts, respectively, compared to 26% (11/42) and 67% (28/42) of the direct carbonate determinations, respectively.
-
Errors on the GIS determinations are on par with those generated by the direct carbonate method and are roughly four to five times worse than errors derived from standard graphite determinations.
-
Six Plio-Pleistocene shells, 14C-dead as determined by graphite analysis, yielded finite direct carbonate ages but yielded expected non-finite GIS ages.
-
Our evaluation of graphite, GIS, and direct carbonate 14C determinations reveals that GIS outperforms direct carbonate, especially for older samples (e.g., > 20 ka BP) and at smaller sample sizes.
-
MICADAS + GIS 14C determinations of small biogenic carbonate samples are an effective alternative to graphite 14C determinations, especially when the precision of graphite is not required. This rapid and reliable dating method is useful for a broad range of applications, including the dating of minute specimens (e.g., small mollusks, foraminifera tests, ostracode valves, fish otoliths, etc.), determining degrees of time averaging across environments and taxa, and inferring temporal changes in population dynamics from post-mortem age distributions.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/RDC.2024.45
Acknowledgments
We appreciate the efforts of Aibhlin Ryan, Fallon Born, and Michelle Olichney as undergraduate laboratory technicians at Northern Arizona University. We thank John Southon for the direct carbonate analyses, Dan Muhs for providing several mollusk shells for the comparative analysis, and Sheila Griffen and Ellen Druffel for sharing the CAHI standard. DSK and EAGS acknowledge funding from the National Science Foundation for the GIS-CHS system (EAR 1855381 and 1919506). Additional support from the National Science Foundation was provided to PGH (EAR CAREER 2041667), JWH (EAR CAREER 1650745), MK (EAR 2127623), and MR (OPP-1744433). PGH acknowledges funding from the Gulf Research Program of the National Academies of Science, Engineering, and Medicine (GRP 2000008410). CF thanks FONDECYT 11200953 and Millennium Nucleus UPWELL NCN19 153. The authors declare no competing interests with this manuscript.