The use of carbonate 'clumped isotope' thermometry as a geochemical technique to determine temperature of formation of a carbonate mineral is predicated on the assumption that the mineral has attained an internal thermodynamic equilibrium. If true, then the clumped isotope signature is dependent solely upon the temperature of formation of the mineral without the need to know the isotopic or elemental composition of coeval fluids. However, anomalous signatures can arise under disequilibrium conditions that can make the estimation of temperatures uncertain by several degrees Celsius. Here we use ab initio calculations to examine the potential disequilibrium mineral signatures that may arise from the incorporation of dissolved inorganic carbon (DIC) species (predominantly aqueous carbonate and bicarbonate ions) into growing crystals without full equilibration with the crystal lattice.We explore theoretically the nature of clumping in the individual DIC species and the composite DIC pool under varying pH, salinity, temperature, and isotopic composition, and speculate about their effects upon the resultant disequilibrium clumping of the precipitates. We also calculate equilibrium clumped signatures for the carbonate minerals calcite, aragonite, and witherite. Our models indicate that each DIC species has a distinct equilibrium clumped isotope signature such that, (witherite). We define the calcite clumped crossover pH as the pH at which the composite D 47 (DIC pool) = D 47 (equilibrium calcite). If disequilibrium D 47 (calcite) is misinterpreted as equilibrium D 47 (calcite), it is possible to overestimate or underestimate the growth temperature by small but significant amounts. Increases in salinity lower the clumped crossover pH and may cause larger effects. Extreme effects of pH, salinity, and temperature, such as between cold freshwater lakes at high latitudes to hot hypersaline environments, are predicted to have sizeable effects on the clumped isotope composition of aqueous DIC pools.In order to determine the most reliable and efficient modeling methods to represent aqueous dissolved inorganic carbon (DIC) species and carbonate minerals, we performed convergence and sensitivity testing on several different levels of theory. We used 4 different techniques for modeling the hydration of DIC: gas phase, implicit solvation (PCM and SMD), explicit solvation (ion with 3 water molecules) and supermolecular clusters (ion plus 21 to 32 water molecules with geometries generated by molecular dynamics). For each solvation technique, we performed sensitivity testing by combining different levels of theory (up to 8 ab initio/hybrid methods, each with up to 5 different sizes of basis sets) to understand the limits of each technique. We looked at the degree of convergence with the most complex (and accurate) models in order to select the most reliable and efficient modeling methods. The B3LYP method combined with the 6-311++G(2d,2p) basis set with supermolecular clusters worked well.

Clumped-isotope" thermometry is an emerging tool to probe the temperature history of surface and subsurface environments based on measurements of the proportion of 13 C and 18 O isotopes bound to each other within carbonate minerals in 13 C 18 O 16 O 2 2À groups (heavy isotope "clumps"). Although most clumped isotope geothermometry implicitly presumes carbonate crystals have attained lattice equilibrium (i.e., thermodynamic equilibrium for a mineral, which is independent of solution chemistry), several factors other than temperature, including dissolved inorganic carbon (DIC) speciation may influence mineral isotopic signatures. Therefore we used a combination of approaches to understand the potential influence of different variables on the clumped isotope (and oxygen isotope) composition of minerals.We conducted witherite precipitation experiments at a single temperature and at varied pH to empirically determine 13 C-18 O bond ordering (D 47 ) and d 18 O of CO 3 2À and HCO 3 À molecules at a 25°C equilibrium. Ab initio cluster models based on density functional theory were used to predict equilibrium 13 C- 18 O bond abundances and d 18 O of different DIC species and minerals as a function of temperature. Experiments and theory indicate D 47 and d 18 O compositions of CO 3 2À and HCO 3 À ions are significantly different from each other. Experiments constrain the D 47 -d 18 O slope for a pH effect (0.011 ± 0.001; 12 P pH P 7). Rapidly-growing temperate corals exhibit disequilibrium mineral isotopic signatures with a D 47 -d 18 O slope of 0.011 ± 0.003, consistent with a pH effect. Our theoretical calculations for carbonate minerals indicate equilibrium lattice calcite values for D 47 and d 18 O are intermediate between HCO 3 À and CO 3 2À. We analyzed synthetic calcites grown at temperatures ranging from 0.5 to 50°C with and without the enzyme carbonic anhydrase present. This enzyme catalyzes oxygen isotopic exchange between DIC species and is present in many natural systems. The two types of experiments yielded statistically indistinguishable results, and these measurements yield a calibration that overlaps with our theoretical predictions for calcite at equilibrium. The slow-growing Devils Hole calcite exhibits D 47 and d 18 O values consistent with lattice equilibrium.Factors influencing DIC speciation (pH, salinity) and the timescale for DIC equilibration, as well as reactions at the mineral-solution interface, have the potential to influence clumped-isotope signatures and the d 18 O of carbonate minerals. In fast-growing carbonate minerals, solution chemistry may be an important factor, particularly over extremes of pH and salinity. If a crystal grows too rapidly to reach an internal equilibrium (i.e., achieve the value for the temperature-dependent mineral lattice equilibrium), it may record the clumped-isotope signature of a DIC species (e.g., the temperature-dependent equilibrium of HCO 3 À ) or a mixture of DIC species, and hence record a disequilibrium mineral composition. For extremely slow-growing cry...

No abstract

No abstract

No abstract