Search

The Cedars is an area in Northern California with a chain of highly alkaline springs resulting from CO2-charged meteorological water interacting with a peridotite body. Serpentinization resulting from this interaction at depth leads to the sequestration of various carbonate minerals into veins accompanied by a release of Ca2+ and OH– enriched water to the surface, creating an environment which promotes rapid precipitation of CaCO3 at surface springs. This environment enables us to apply the recently developed Δ47-Δ48 dual clumped isotope analysis to probe kinetic isotope effects (KIEs) and timescales of CO2 transformation in a region with the potential for geological CO2 sequestration. We analyzed CaCO3 recovered from various localities and identified significant kinetic fractionations associated with CO2 absorption in a majority of samples, characterized by enrichment in Δ47 values and depletion in Δ48 values relative to equilibrium. Surface floes exhibited the largest KIEs (ΔΔ47: 0.163‰, ΔΔ48: −0.761‰). Surface floe samples begin to precipitate out of solution within the first hour of CO2 absorption, and the dissolved inorganic carbon (DIC) pool requires a residence time of >100 h to achieve isotopic equilibria. The Δ48/Δ47 slope of samples from the Cedars (−3.223 ± 0.519) is within the range of published theoretical values designed to constrain CO2 hydrolysis-related kinetic fractionation (−1.724 to −8.330). The Δ47/δ18O slope (−0.009 ± 0.001) and Δ47/δ13C slope (−0.009 ± 0.001) are roughly consistent with literature values reported from a peridotite in Oman of −0.006 ± 0.002 and −0.005 ± 0.002, respectively. The consistency of slopes in the multi-isotope space suggests the Δ47-Δ48 dual carbonate clumped isotope framework can be applied to study CO2-absorption processes in applied systems, including sites of interest for geological sequestration.

The authors regret that in the original article, the kinetic fractionation factors (KFFs) related to the CO2 hydroxylation reaction (CO2(aq) + OH– → HCO3–) were calculated using atmospheric δ18OCO2 values of 0.0 ± 0.5 ‰ (VPDB) while the correct values are 10.25 ± 0.5 ‰ (VPDB). We erroneously employed the atmospheric δ18OCO2 values from the NOAA CO2 global network (Trolier et al., 1996) by not converting atmospheric δ18OCO2 from the VPDB-CO2 scale to the VPDB scale through the following expression (e.g. Srivastava and Verkouteren, 2018): [Formula presented] This correction affects the reported oxygen isotope KFFs for CO2 hydroxylation to form HCO3– prior to CaCO3 precipitation. We have corrected Table 7 and Fig. 13 accordingly (see below). The composition of atmospheric CO2 was also incorrectly plotted on the VPDB scale in Figs. 1, 8 and 15 of our original manuscript based on previous works (Fig. 8 in Clark and Fontes, 1990; Fig. 2 in Clark et al., 1992). Importantly, the modelling outputs and main conclusions of the original article are not affected by these corrections. Nevertheless, the updated KFFs resolve an apparent discrepancy between the laboratory and field studies in our original article. Moreover, the oxygen isotope KFF for CO2 hydroxylation now displays a positive temperature dependence over the 4–28 °C range, which was not statistically significant in the original article (Fig. 13c, r2 = 0.86, n = 10, p-value