Carbonate clumped-isotope geothermometry is a tool used to reconstruct formation or (re)equilibration temperatures of carbonate bearing minerals, including carbonate groups substituted into apatite. It is based on the preference for isotopologues with multiple heavy isotopes (for example,13C16O218O2- groups) to be more abundant at equilibrium than would be expected if all isotopes were randomly distributed amongst all carbonate groups. Because this preference is only a function of temperature, excesses of multiply substituted species can be used to calculate formation temperatures without knowledge of the isotopic composition of water from which the mineral precipitated or other phases with which it may have equilibrated. However, the measured temperature could be modified after mineral growth if exchange of isotopes amongst carbonate groups within the mineral has occurred through internal isotope-exchange reactions. Because these exchange reactions occur through thermally activated processes, their rates depend on temperature and increase at higher temperatures. Thus internal isotope-exchange reactions could lead to effective reequilibration at high temperatures, overprinting the original temperatures recorded during mineral growth. We measured clumped-isotope temperatures in carbonate bearing minerals (including apatites) from several carbonatites to constrain the kinetics of these internal isotope-exchange reactions. We observe two key trends for clumped-isotope temperatures in carbonatites: (i) clumped-isotope temperatures of apatites and carbonate-bearing minerals decrease with increasing intrusion depth and (ii) apatites record lower clumped-isotope temperatures than carbonate minerals from the same intrusion. We additionally conducted heating experiments at different temperatures to derive the temperature dependence for the rate constants that describe the alteration of clumped-isotope temperatures with time in calcites and apatites. We find that calcites exhibit complex kinetics as has been seen in previous studies. To quantify these results, we constructed a model that incorporates both diffusion of isotopes through the crystal lattice and isotope-exchange reactions between adjacent carbonate groups. We tested this model through comparison to previous measurements of optical calcites and brachiopods and to samples with known cooling histories and find that the model is able to reasonably capture kinetic data from previous experiments and the observed clumped-isotope temperatures of calcites assuming geologically reasonable cooling rates. A similar model for apatite over-predicts the observed clumped-isotope temperatures found in natural samples; we hypothesize this discrepancy is the result of annealing of radiation damage in our experiments, which lowers the diffusivity and rate of isotope exchange of carbonate groups compared to damaged natural samples. Finally, we constructed models to explore how heating can alter recorded clumped-isotope temperatures. Our model predicts that samples change in clumped-isotope temperatures in two stages. The first stage changes the recorded clumped isotope temperatures by <1 °C if held at 75 °C for 100 million years and by up to ∼40 °C if held at 120 °C, but the clumped-isotope temperatures does not reach ambient values through this low-temperature mechanism. A second, slower change becomes effective at temperatures above 150 °C and can take the measured clumped-isotope temperature up to the true ambient temperature. This result implies that old (hundreds of million years) samples that have only experienced mild (<100-125 °C) thermal histories could exhibit small but measurable (order 10 °C) changes in their clumped-isotope temperatures. We compared this heating model to clumped-isotope measurements from paleosol samples from the Siwalik Basin in Nepal, which were buried up to 5 km and then rapidly exhumed to the surface; these samples often do not give reasonable surface temperatures. The modeled temperatures agree with measured temperatures of these samples, suggesting that partial re-equilibration during shallow crustal burial is responsible for their elevated clumpedisotope temperatures.

The history of atmospheric O₂ partial pressures (Po₂) is inextricably linked to the coevolution of life and Earth's biogeochemical cycles. Reconstructions of past Po₂ rely on models and proxies but often markedly disagree. We present a record of Po₂ reconstructed using O₂/N₂ ratios from ancient air trapped in ice. This record indicates that Po₂ declined by 7 per mil (0.7%) over the past 800,000 years, requiring that O₂ sinks were ~2% larger than sources. This decline is consistent with changes in burial and weathering fluxes of organic carbon and pyrite driven by either Neogene cooling or increasing Pleistocene erosion rates. The 800,000-year record of steady average carbon dioxide partial pressures (Pco₂) but declining Po₂ provides distinctive evidence that a silicate weathering feedback stabilizes Pco₂ on million-year time scales.

Exhumation and accompanying retrograde metamorphism alter the compositions and textures of metamorphic rocks through deformation, mineral–mineral reactions, water–rock reactions, and diffusion-controlled intra- and inter-mineral atomic mobility. Here, we demonstrate that these processes are recorded in the clumped- and single-isotope (δ13C and δ18O) compositions of marbles, which can be used to constrain retrograde metamorphic histories. We collected 27 calcite and dolomite marbles along a transect from the rim to the center of the metamorphic core-complex of Naxos (Greece), and analyzed their carbonate single- and clumped-isotope compositions. The majority of Δ47 values of whole-rock samples are consistent with exhumation- controlled cooling of the metamorphic complex. However, the data also reveal that water–rock interaction, deformation driven recrystallization and thermal shock associated with hydrothermal alteration may considerably impact the overall distribution of Δ47 values. We analyzed specific carbonate fabrics influenced by deformation and fluid–rock reaction to study how these processes register in the carbonate clumped-isotope system. Δ47 values of domains drilled from a calcite marble show a bimodal distribution. Low Δ47 values correspond to an apparent temperature of 260 °C and are common in static fabrics; high Δ47 values correspond to an apparent temperature of 200 °C and are common in dynamically recrystallized fabrics. We suggest that the low Δ47 values reflect diffusion-controlled isotopic reordering during cooling, whereas high Δ47 values reflect isotopic reordering driven by dynamic recrystallization. We further studied the mechanism by which dynamic recrystallization may alter Δ47 values by controlled heating experiments. Results show no significant difference between laboratory reactions rates in the static and dynamic fabrics, consistent with a mineral-extrinsic mechanism, in which slip along crystal planes was associated with atomic-scale isotopic reordering in the calcite lattice. An intrinsic mechanism (enhanced isotopic reordering rate in deformed minerals) is contraindicated by these experiments. We suggest that Δ47 values of dynamically recrystallized fabrics that form below the diffusion-controlled blocking-temperature for calcite constrain the temperature of deformation. We find that Δ47-based temperatures of static fabrics from Naxos marbles are ∼60–80 °C higher than commonly observed in slowly cooled metamorphic rocks, and would suggest cooling rates of ∼105°CMyr−1. A similar thermal history is inferred for dolomite marbles from the core vicinity, which preserve apparent temperatures up to 200 °C higher than a typical blocking temperature (∼300 °C). This finding could be explained by a hydrothermal event driving a brief thermal pulse and locally resetting Δ47 values. Rapid cooling of the core-complex region is consistent with a compilation of published cooling ages and a new apatite U–Th/He age, associating the thermal event with the emplacement of a granodiorite pluton at ∼12 Ma.

Sources of methane to sedimentary environments are commonly identified and quantified using the stable isotopic compositions of methane. The methane "clumped-isotope geothermometer", based on the measurement of multiply substituted methane isotopologues (¹³CH3D and ¹²CH2D2), shows promise in adding new constraints to the sources and formational environments of both biogenic and thermogenic methane. However, questions remain about how this geothermometer behaves in systems with mixtures of biogenic and thermogenic gases and different biogenic environments. We have applied the methane clumped-isotope thermometer to a mixed biogenic-thermogenic system (Antrim Shale, USA) and to biogenic gas from gas seeps (Santa Barbara and Santa Monica Basin, USA), a pond on the Caltech campus, and methanogens grown in pure culture. We demonstrate that clumped-isotope based temperatures add new quantitative constraints to the relative amounts of biogenic vs. thermogenic gases in the Antrim Shale indicating a larger proportion (~50%) of thermogenic gas in the system than previously thought. Additionally, we find that the clumped-isotope temperature of biogenic methane appears related to the environmental settings in which the gas forms. In systems where methane generation rates appear to be slow (e.g., the Antrim Shale and gas seeps), microbial methane forms in or near both internal isotopic equilibrium and hydrogen-isotope equilibrium with environmental waters. In systems where methane forms rapidly, microbial methane is neither in internal isotopic equilibrium nor hydrogen-isotope equilibrium with environmental waters. A quantitative model of microbial methanogenesis that incorporates isotopes is proposed to explain these results.

The Athel silicilyte is an enigmatic, hundreds of meters thick, finely laminated quartz deposit, in which silica precipitated in deep water (>~100-200 m) at the Ediacaran-Cambrian boundary in the South Oman Salt Basin. In contrast, Meso-Neoproterozoic sinks for marine silica were dominantly restricted to peritidal settings. The silicilyte is known to contain sterane biomarkers for demosponges, which today are benthic, obligately aerobic organisms. However, the basin has previously been described as permanently sulfidic and time-equivalent shallow-water carbonate platform and evaporitic facies lack silica. The Athel silicilyte thus represents a unique and poorly understood depositional system with implications for late Ediacaran marine chemistry and paleoecology. To address these issues, we made petrographic observations, analyzed biomarkers in the solvent-extractable bitumen, and measured whole-rock iron speciation and oxygen and silicon isotopes. These data indicate that the silicilyte is a distinct rock type both in its sedimentology and geochemistry and in the original biology present as compared to other facies from the same time period in Oman. The depositional environment of the silicilyte, as compared to the bounding shales, appears to have been more reducing at depth in sediments and possibly bottom waters with a significantly different biological community contributing to the preserved biomarkers. We propose a conceptual model for this system in which deeper, nutrient-rich waters mixed with surface seawater via episodic mixing, which stimulated primary production. The silica nucleated on this organic matter and then sank to the seafloor, forming the silicilyte in a sediment-starved system. We propose that the silicilyte may represent a type of environment that existed elsewhere during the Neoproterozoic. These environments may have represented an important locus for silica removal from the oceans.

Methane is a potent greenhouse gas, and there are concerns that its natural emissions from the Arctic could act as a substantial positive feedback to anthropogenic global warming. Determining the sources of methane emissions and the biogeochemical processes controlling them is important for understanding present and future Arctic contributions to atmospheric methane budgets. Here we apply measurements of multiply-substituted isotopologues, or clumped isotopes, of methane as a new tool to identify the origins of ebullitive fluxes in Alaska, Sweden and the Arctic Ocean. When methane forms in isotopic equilibrium, clumped isotope measurements indicate the formation temperature. In some microbial methane, however, non-equilibrium isotope effects, probably related to the kinetics of methanogenesis, lead to low clumped isotope values. We identify four categories of emissions in the studied samples: thermogenic methane, deep subsurface or marine microbial methane formed in isotopic equilibrium, freshwater microbial methane with non-equilibrium clumped isotope values, and mixtures of deep and shallow methane (i.e., combinations of the first three end members). Mixing between deep and shallow methane sources produces a non-linear variation in clumped isotope values with mixing proportion that provides new constraints for the formation environment of the mixing end-members. Analyses of microbial methane emitted from lakes, as well as a methanol-consuming methanogen pure culture, support the hypothesis that non-equilibrium clumped isotope values are controlled, in part, by kinetic isotope effects induced during enzymatic reactions involved in methanogenesis. Our results indicate that these kinetic isotope effects vary widely in microbial methane produced in Arctic lake sediments, with non-equilibrium δ18 values spanning a range of more than 5‰.