# Rapid microbial methanogenesis during CO2 storage in hydrocarbon reservoirs

R. L. Tyne, P. H. Barry, M. Lawson, D. J. Byrne, O. Warr, H. Xie, D. J.

Hillegonds, M. Formolo, Z. M. Summers, B. Skinner, et al.

► To cite this version:

R. L. Tyne, P. H. Barry, M. Lawson, D. J. Byrne, O. Warr, et al. Rapid microbial methanogenesis during CO2 storage in hydrocarbon reservoirs. Nature, 2021, 600, pp.670-674. 10.1038/s41586-021-04153-3. insu-03619882

HAL Id: insu-03619882

https://insu.hal.science/insu-03619882

Submitted on 25 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Rapid microbial methanogenesis during CO₂ storage in hydrocarbon reservoirs

https://doi.org/10.1038/s41586-021-04153-3

Received: 6 March 2021

Accepted: 14 October 2021

Published online: 22 December 2021

Open access

<figure><img src="image_2.png" /></figure>

Check for updates

R. L. Tyne<fnref n="1" />, P. H. Barry<fnref n="1" />,², M. Lawson<fnref n="3" />,⁹, D. J. Byrne⁴, O. Warr⁵, H. Xie⁶, D. J. Hillegonds<fnref n="1" />, M. Formolo⁷, Z. M. Summers⁸, B. Skinner⁷, J. M. Eiler⁶ & C. J. Ballentine<fnref n="1" />

Carbon capture and storage (CCS) is a key technology to mitigate the environmental impact of carbon dioxide (CO₂) emissions. An understanding of the potential trapping and storage mechanisms is required to provide confidence in safe and secure CO₂ geological sequestration<fnref n="1" />,². Depleted hydrocarbon reservoirs have substantial CO₂ storage potential<fnref n="1" />,<fnref n="3" />, and numerous hydrocarbon reservoirs have undergone CO₂ injection as a means of enhanced oil recovery (CO₂-EOR), providing an opportunity to evaluate the (bio)geochemical behaviour of injected carbon. Here we present noble gas, stable isotope, clumped isotope and gene-sequencing analyses from a CO₂-EOR project in the Olla Field (Louisiana, USA). We show that microbial methanogenesis converted as much as 13–19% of the injected CO₂ to methane (CH₄) and up to an additional 74% of CO₂ was dissolved in the groundwater. We calculate an in situ microbial methanogenesis rate from within a natural system of 73–109 millimoles of CH₄ per cubic metre (standard temperature and pressure) per year for the Olla Field. Similar geochemical trends in both injected and natural CO₂ fields suggest that microbial methanogenesis may be an important subsurface sink of CO₂ globally. For CO₂ sequestration sites within the environmental window for microbial methanogenesis, conversion to CH₄ should be considered in site selection.

A possible method for reducing current greenhouse gas emission rates is carbon capture and storage (CCS). Long-term carbon dioxide (CO₂) trapping mechanisms include structural or stratigraphic trapping in stable geological configurations, dissolution into pore fluids (solubility trapping)⁴, carbonate mineralization (precipitation and mineral trapping)⁵ or adsorption (for example, onto coal)⁶. Typically, during CO₂ enhanced oil recovery (CO₂-EOR), a proportion of the injected CO₂ remains within the reservoir post-injection<fnref n="3" />, providing an analogue for investigating and quantifying processes within CCS sites, over decadal timescales.

Here we investigate the behaviour of CO₂ within the Olla Oil Field, Louisiana, USA (Fig. 1), which was CO₂-EOR flooded in the 1980s, and compare its (bio)geochemical composition with the adjacent Nebo-Hemphill Oil Field, which was never subjected to CO₂-EOR (Methods). Approximately 9 × 10⁷ m<fnref n="3" /> (standard temperature and pressure (STP)) of injected CO₂ has been retained in the Olla reservoirs post-injection⁷. The injected CO₂ was sourced from the Black Lake Oil Field, which is located adjacent to the Sabine Island Complex, a basement high and likely conduit for mantle-derived fluids⁸. Although the Olla Oil Field has been geochemically characterized by several previous studies⁹⁻¹¹ (Methods), this study integrates noble gas and clumped isotope data with stable isotope and microbiological data to investigate the fate of the injected CO₂. Previous studies have suggested that microbial hydrocarbon degradation and methanogenesis may occur

within both the Olla and Nebo-Hemphill oil fields⁹⁻¹¹. Microbial methanogenesis takes acetate, methylated compounds, or hydrogen and CO₂ (hydrogenotrophic methanogenesis, equation (1)) as its starting point; the latter mechanism is most relevant to this study. Hydrogenotrophic methanogenesis results in ¹³C enrichment of residual CO₂ and the formation of methane (CH₄) initially ¹³C depleted compared with the source CO₂, which increases to the value of the consumed CO₂.

$$CO_2 + 4H_2 \rightarrow CH_4 + 2H_2O \quad (1)$$

Given their proximity and comparable geological histories (including hydrocarbon production from similar reservoirs; Methods), we assume that the Nebo-Hemphill and Olla oil fields initially had comparable geochemical compositions, before CO₂-EOR at Olla. For example, the Olla and Nebo-Hemphill oil fields have broadly overlapping pH, temperatures, salinities and ion chemistry (Extended Data Table 1). Despite CO₂ injection ceasing in 1986, CO₂ concentrations within the Olla Oil Field (18 ± 12 mol% (STP) across all reservoirs) are greater than breakthrough concentrations limits (10%)⁷ and significantly greater than those measured at Nebo-Hemphill (1.02 ± 0.69 mol% (STP)). We observe a higher δ¹³C<sub>VPDB</sub> of CO₂ at Olla (13.5 ± 3.4‰) compared with both Nebo-Hemphill (4.8 ± 4.9‰) and the injected CO₂ (0.85 ± 0.92‰)¹⁰ (Extended Data Table 2), where δ¹³C = [(¹³C/¹²C)<sub>sample</sub>/(¹³C/¹²C)<sub>standard</sub>] - 1 and VPDB is Vienna PeeDee Belemnite. Injection alone cannot account

<footnote n="1">¹Department of Earth Sciences, University of Oxford, Oxford, UK. ²Department of Marine Chemistry and Geochemistry, Woods Hole Oceanographic Institution, Woods Hole, MA, USA.</footnote>

<footnote n="3">³ExxonMobil Upstream Business Development, Spring, TX, USA. ⁴CRPG-CNRS, Université de Lorraine, Nancy, France. ⁵Department of Earth Sciences, University of Toronto, Toronto, Ontario, Canada. ⁶Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, USA. ⁷ExxonMobil Upstream Integrated Solutions, Spring, TX, USA. ⁸ExxonMobil Research and Engineering Co., Virginia, NJ, USA. ⁹Present address: Aker BP, Stavanger, Norway. ¹⁰e-mail: rebecca.tyne@earth.ox.ac.uk; pbarry@whoi.edu; miclawson@gmail.com; chris ballentine@earth.ox.ac.uk</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** | Map of study area showing locations of the Olla and Nebo-Hemphill oil fields as well as the Black Lake Oil Field, from which the injected CO₂ was sourced. Only the Olla Oil Field contains injected CO₂. The inset on the right</figcaption></figure>

shows an expanded view with individual sample locations (nearby urban areas are denoted by stars) and a stratigraphic column showing the relevant lithologic units⁹.

for the higher δ¹³C of CO₂ values at Olla; therefore, the integration of complementary datasets is required to evaluate additional in situ processes.

Owing to their inert nature, noble gases provide a powerful tool for tracing and quantifying physiochemical processes associated with CO₂ injection⁴,⁸,¹²⁻¹⁹. The average air-corrected helium isotope (³He/⁴He) ratio relative to the atmospheric ratio (R/R<sub>A</sub>, where R = ³He/⁴He<sub>sample</sub> and R<sub>A</sub> = ³He/⁴He<sub>air</sub>) was determined for the Olla and Nebo-Hemphill oil fields to be 1.76 ± 0.31 R<sub>A</sub> and 0.46 ± 0.12 R<sub>A</sub>, respectively (Extended Data Fig. 1, Extended Data Table 2). The elevated He isotope values in the Olla Oil

Field require an enhanced mantle-derived noble gases contribution, which is supported by additional noble gas data (Extended Data Fig. 1, Methods), associated with injected CO₂, which were sourced from the Black Lake Oil Field.

In addition to characterizing fluid origin, a combined noble gas and CO₂ isotope approach provides insight into processes associated with CO₂ trapping (for example, ref. ⁴) or methanogenesis following injection. ³He is inert and insoluble, with no significant sources within the crust. Thus, variations in CO₂/³He post-injection are directly attributable to the addition or removal of CO₂ within the system⁴,²⁰. Although

<figure><img src="image_4.png" /><figcaption>**Fig. 2** | The δ¹³C of CO₂ in the Olla (CO₂ injected field) samples. a, The δ¹³C of CO₂ as a function of the CO₂/³He ratio. The dashed lines show endmember methanogenesis and dissolution (pH 7) fractionation trajectories. The tick marks represent the total amount of CO₂ trapping within the system, relative to sample O5. The shaded region represents trapping by the combination of both microbial methanogenesis and dissolution. The upper and lower methanogenesis:dissolution ratios (M:D) are 0.33 and 0.19, respectively, showing that dissolution accounts for approximately three times more CO₂</figcaption></figure>

removal (M:D = 0.26) than microbial methanogenesis. The lines labelled 'consumption' show the portion of original injected CO₂ that has been removed by net microbial methanogenesis. b, The δ¹³C of CO₂ as a function of the δ¹³C of CH₄. The shaded region represents the CO₂ isotopic composition of the injectate from the Black Lake Oil Field into the Olla system. The Olla data are consistent with thermal re-equilibration with both reservoir temperatures (solid lines) and microbial methanogenesis (dotted lines). The 1σ level of uncertainty is within the symbol size.
two samples at Olla (samples O5 and O6) have CO₂/³He within the man-
tle range (2 ± 1 × 10⁹), consistent with the injection of mantle-derived
fluids, the majority of the ratios are lower and indicate CO₂ removal
(Fig. 2a, Extended Data Table 2). The difference between the highest
CO₂/³He (O5, considered the most pristine sample) and the remaining
samples provides a conservative estimate of post-injection CO₂ trap-
ping/consumption of between 39% and 89%.

A decrease in the CO₂/³He alongside increasing ⁴He and ²⁰Ne con-
centrations is observed within the Olla Oil Field (Extended Data Fig. 2,
Methods). In the subsurface, ⁴He is produced via radiogenic decay of
uranium and thorium and accumulates in formation water¹², which also
contains atmosphere-derived ²⁰Ne. An anticorrelation between CO₂/³He
and ⁴He and ²⁰Ne is observed in natural CCS analogues globally⁴, and
indicates that the extent of water contact controls the magnitude of
CO₂ trapping via dissolution or precipitation of CO₂ (refs. ⁴,¹⁵). Within
the observed pH range (7.0 ± 0.5) at Olla, dissolution will dominate over
precipitation as a CO₂ trapping mechanism⁴,¹⁰. However, as described
above, the δ¹³C of CO₂ in the Olla gases is more enriched than can be
explained by these processes alone (Fig. 2a).

Independent of the source of the injected CO₂, such positive δ¹³C values suggest significant modification from a starting composition that is unlikely to be heavier than 2‰ (refs. ²¹,²²). Notably, sorption to coal, hydrocarbon biodegradation to CO₂, and subsequent consumption and mixing processes could result in an increase in the δ¹³C of CO₂ (Methods). However, even from the most pristine sample (O5), it is unlikely that any of these processes could generate a fluid with the elevated concentrations of the isotopically enriched CO₂ observed at Olla.

Further constraints on the origins and processes affecting these
gases are provided by the molecular-average and clumped isotope
compositions of the associated CH₄. We observe a relatively high δ¹³C
of CO₂ (13.6 ± 3.2‰) and CH₄ (−47.7 ± 1.7‰) at Olla compared with
those observed at Nebo-Hemphill (3.5 ± 3.6‰ and −56.4 ± 3.6‰ for
CO₂ and CH₄, respectively) (Fig. 2b, Extended Data Fig. 3a, Extended
Data Table 2). The higher δ¹³C of CH₄ at Olla could be interpreted
as evidence that the Olla Oil Field contains a higher proportion of
thermogenic gas, or thermogenic gas of higher thermal maturity,
whereas the Nebo-Hemphill Oil Field contains more microbial gas or
thermogenic gas of lower thermal maturity. Using our multi-isotope
approach, we present a model that is consistent with the observed
higher δ¹³C values of CH₄ and an increase in the δ¹³C of residual CO₂.
In this model, Olla is assumed to have had a pre-CO₂-EOR composition
resembling Nebo-Hemphill, but following CO₂-EOR, microbial activity
converted significant amounts of injected CO₂ to CH₄ (Extended
Data Fig. 4). We arrive at this conclusion on the basis of the following
observations.

First, we report a difference between the δ¹³C of CH₄ and the δ¹³C of CO₂ of 53.7–64.4‰ at the Olla Oil Field and of 50.8–68.1‰ at the Nebo-Hemphill Oil Field (Fig. 2b). The observed carbon isotope fractionations between co-existing CO₂ and CH₄ are consistent with thermodynamic isotopic equilibrium at temperatures between 44.6 °C and 76.3 °C at Olla and between 29.8 °C and 85.3 °C at Nebo-Hemphill²³, which overlap the present-day reservoir temperatures within these fields (29.2–50.7 °C and 29.7–57.1 °C, respectively), which suggests that the systems are approaching isotopic equilibrium under current reservoir conditions (Methods). The isotopic approach to equilibrium under reservoir conditions appears to be a result of microbial cycling of carbon (that is, methanogenesis and anaerobic oxidation of methane (AOM); evidence for AOM is apparent in the clumped isotopologues, see below). Equilibrium with reservoir conditions between the δ¹³C of CO₂ and the δ¹³C of CH₄ has previously been observed under similar geological conditions²⁴.

Second, we observe a range in the two measurable clumped isotopo-
logues of CH₄, Δ¹³CH₃D and Δ¹²CH₃D₂, at Olla of 3.45–5.62‰ and 9.13–
12.4‰ (Fig. 3, Extended Data Table 2), respectively, which correspond
to apparent temperature ranges of 29⁻¹⁴°C to 128⁺¹⁷°C and 87⁻¹³°C to

<figure><img src="image_5.png" /></figure>

**Fig. 3 | Δ¹²CH₃D₂ versus Δ¹³CH₃D of the measured Olla and Nebo-Hemphill samples.** The clumped isotopologue space illustrates whether measured CH₄ is at internal thermodynamic equilibrium (black line) or not. The thermodynamic equilibrium curve is calculated following ref. ⁴¹. The shaded cross represents thermal equilibration to current reservoir temperatures in each isotopologue. The arrows represent the theoretical trends for methanogenesis (dark grey) and AOM (light grey) ²⁹. The 1σ level of uncertainty is shown on the measured samples.

132<sup>+30</sup><sub>-23</sub> °C, respectively, and are approaching equilibrium with current reservoir conditions. Clumped isotope compositions provide an independent constraint on the origin of CH₄ in petroleum systems<sup>25–30</sup>. Thermogenic CH₄ appears to dominantly form under internal isotopic equilibrium<sup>26,30</sup> whereas microbial CH₄ is highly variable in its clumped isotope compositions, with both non-equilibrium<sup>28,29</sup> and apparent equilibrium signatures<sup>24,26,28,29,31–34</sup>. For Olla, we expect thermogenic-CH₄-generation temperatures to have exceeded 163 ± 18 °C, from independent maturity constraints provided by biomarkers<sup>35</sup>. Two samples from Olla (O4 and O6) appear to be in internal isotopic equilibrium (about 125 °C and about 95 °C) but have lower apparent temperatures than expected for pure thermogenic CH₄. These temperatures are consistent with a component of microbial CH₄ formed at or close to equilibrium with current reservoir conditions. The remaining samples have Δ¹³CH₃D apparent temperatures of 28.7<sup>+13.6</sup><sub>-12.3</sub> °C to 74.0<sup>+16.8</sup><sub>-15.0</sub> °C, within error of the present-day reservoir temperature for these fluids (Methods), consistent with in situ microbial methanogenesis and AOM. Furthermore, these samples exhibit deficits in Δ¹²CH₂D₂, consistent with fluids that are dominated by microbial CH₄<sup>28,29,34</sup> in a system approaching equilibrium via active methanogenesis and AOM. In contrast, the Δ¹³CH₃D-based temperatures of 77.9<sup>+13.5</sup><sub>-12.3</sub> °C to 166<sup>+34</sup><sub>-29</sub> °C at Nebo-Hemphill suggest that these fluids are dominated by thermogenic CH₄, although, similarly to Olla, the Δ¹²CH₂D₂ clumped isotopologues also suggest a microbial CH₄ contribution.

Third, we identify evidence for the biodegradation of hydrocarbons in the molecular geochemistry of oils produced from both the Olla and the Nebo-Hemphill oil fields (Methods). In the case of Olla, we also report an elevated δ¹³C of propane (Extended Data Fig. 3b). The presence of oils co-existing with propane in these reservoirs preclude the possibility of high-temperature cracking as the cause of the positive δ¹³C values and suggests that the propane has been subjected to biodegradation³⁶, consistent with low-temperature microbial activity within the field.
Fourth, small subunit (SSU) ribosomal RNA gene sequencing of microbial communities identifies the presence of hydrogenotrophic methanogens, methanol and methylotrophic methanogens, and anaerobic methanotrophs (ANME) assigned sequences (Extended Data Fig. 5). These were sampled from the formation water in two wells within each field alongside previous microbiology analysis of these fields³⁷. Co-occurrence of both methanogenic and methanotrophic Archaea shows the potential for conversion of CO₂ to CH₄ and vice versa, consistent with isotopic signatures at intramolecular and intermolecular equilibrium. However, equilibration is probably rate limited by AOM, as the abundance of methanogens was 100 times higher than ANME, and AOM probably proceeds at a lower rate per cell.

Combined independent geochemical and microbiological data support our assertion that oil field microbes convert injected CO₂ to CH₄ at Olla. This study identifies the microbial processes associated with the injection of CO₂ into the deep subsurface operating on decadal timescales and the scale of the processing.

To quantify the impact of these microbial processes on the Olla Oil Field, we have constructed a model that considers both dissolution and microbial methanogenesis (Fig. 2a, Methods). From our most pristine sample (OS, highest CO₂/³He), we find that to match the isotopic composition of CO₂ with its corresponding CO₂/³He, dissolution must account for 3.1 ± 1.1 times more CO₂ removal than hydrogenotrophic methanogenesis. From this model, net microbial methanogenesis is estimated to have consumed a minimum of between 13% and 19% of the post-injection CO₂, depending on sample location, whereas dissolution is responsible for removal of as much as 74% more CO₂. This estimate is independently confirmed when CO₂ concentrations are compared with the δ¹³C of CO₂ values (Extended Data Fig. 6). On the basis of the elevated δ¹³C of CO₂, we recognize that our most pristine sample has probably also undergone a significant degree of modification, consistent with our estimate being conservative, and that CO₂ concentrations post-injection were likely to be higher and variable across the field.

Our combined noble gas and stable isotope approach allows estimates of in situ net microbial methanogenesis rates from within a natural system to be calculated. By extrapolating our results over the 29 years between the cessation of injection (1986) and sampling (2015), and assuming 13–19% microbial consumption of CO₂ since injection (from the hybrid methanogenesis–dissolution model, Fig. 2) with the remaining injected CO₂ volume, we calculate that a minimum of 1.15 × 10⁷–1.72 × 10⁷ m³ (STP) of microbial CH₄ has been produced at a minimum rate of 73–109 mmol CH₄ m⁻³ (STP) yr⁻¹. Previous estimates for CO₂ reduction following methanogenic oil degradation by hydrogenotrophic methanogens in lab microcosm incubations at similar temperatures are significantly lower (about 0.01–0.15 mmol CH₄ m⁻³ (STP) yr⁻¹)³⁸ by comparison.

This identification and quantification of CO₂-EOR-related microbial methanogenesis at Olla has significant implications, and warrants reconsideration of reservoir processes at similar sites worldwide. The fate of injected CO₂ in other oil fields (for example, Belridge Diatomite Formation, Lost Hills Oil Field, California, USA¹⁷,³⁹) has never been systematically explored using this integrated approach. Similarly to Olla, the Lost Hills Field has an elevated δ¹³C of dissolved inorganic carbon (DIC) (20.16–23.61‰)⁴⁰ compared with the injected δ¹³C of CO₂ (–30.1‰)³⁹. The fractionation between δ¹³C of DIC and CH₄ suggests that the Lost Hills system is approaching equilibrium with current reservoir conditions and the inverse correlation between the δ¹³C of the post-injection residual CO₂ and concentration is consistent with microbial methanogenesis. Furthermore, data from Nebo-Hemphill, as well as naturally CO₂-rich hydrocarbon systems in the Pannonian Basin, Hungary, are also consistent with significant microbial methanogenesis based on the CO₂/³He and the δ¹³C of CO₂ (Extended Data Fig. 7, Methods). These examples suggest that hydrogenotrophic methanogenesis probably occurs in both natural and injection-perturbed systems when reservoirs have favourable physiochemical conditions.

We conclude that although methanogenesis is not the dominant CO₂ sink, it can represent a substantial process within both natural and perturbed CO₂ fields and may be significantly more prevalent than previously considered. Even with the most conservative estimates of CO₂ conversion by methanogenesis in the Olla Oil Field, we find that, so far, as much as 13% to 19% of the emplaced CO₂ in sampled sections of the field has been consumed by methanogens. This is less than the amount of CO₂ that has been dissolved into water, but similarly occurs at significant rates on engineering timescales. CH₄ is less soluble and more mobile than CO₂ and therefore there is an enhanced risk of gas loss associated with CH₄ production due to microbial methanogenesis. Depleted hydrocarbon reservoirs have the second-largest CO₂ storage potential¹,³ after deep saline aquifers. If physicochemical and environmental conditions within these reservoirs (for example, low temperature, low salinity and high substrate availability) are conducive to microbial methanogenesis, we suggest that these microbial processes should be considered as criteria for future CCS site selection to ensure low-risk long-term storage. In addition, integrated studies like this may prove essential to effectively monitor CO₂ storage and the biogeochemical processes that result as a consequence of it.

## Online content

Any methods, additional references, Nature Research reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41586-021-04153-3.

1. Leung, D. Y. C., Caramanna, G. & Maroto-Valer, M. M. An overview of current status of carbon dioxide capture and storage technologies. *Renew. Sustain. Energy Rev.* **39**, 426–443 (2014).

2. Bui, M. et al. Carbon capture and storage (CCS): the way forward. *Energy Environ. Sci.* **11**, 1062–1176 (2018).

3. Stewart, R. J. & Haszeldine, R. S. *Carbon Accounting for Carbon Dioxide Enhanced Oil Recovery* (Scottish Carbon Capture & Storage, 2014); https://www.sccs.org.uk/images/expertise/misc/SCCS-CO2-EOR-JIP-Carbon-Balance.pdf

4. Gilfillan, S. M. V. et al. Solubility trapping in formation water as dominant CO₂ sink in natural gas fields. *Nature* **458**, 614–618 (2009).

5. Matter, J. M. et al. Rapid carbon mineralization for permanent disposal of anthropogenic carbon dioxide emissions. *Science* **352**, 1312–1314 (2016)

6. Busch, A., Gensterblum, Y., Krooss, B. M. & Littke, R. Methane and carbon dioxide adsorption–diffusion experiments on coal: upscaling and modeling. *Int. J. Coal Geol.* **60**, 151–168 (2004).

7. Moffat, C. I. Post-economic Evaluation 2800 foot Sand Pilot Project Olla Field, LaSalle Parish Louisiana (1990).

8. Byrne, D. J., Barry, P. H., Lawson, M. & Ballentine, C. J. The use of noble gas isotopes to constrain subsurface fluid flow and hydrocarbon migration in the East Texas Basin. *Geochim. Cosmochim. Acta* **268**, 186–208 (2020).

9. McIntosh, J. C., Warwick, P. D., Martini, A. M. & Osborn, S. G. Coupled hydrology and biogeochemistry of Paleocene–Eocene coal beds, northern Gulf of Mexico. *Geol. Soc. Am. Bull.* **122**, 1248–1264 (2010).

10. Shelton, J. L., McIntosh, J. C., Warwick, P. D. & Yi, A. L. Z. Fate of injected CO₂ in the Wilcox Group, Louisiana, Gulf Coast Basin: chemical and isotopic tracers of microbial–brine-rock–CO₂ interactions. *Appl. Geochem.* **51**, 155–169 (2014).

11. Shelton, J. L., McIntosh, J. C., Warwick, P. D. & McCray, J. E. Impact of formation water geochemistry and crude oil biodegradation on microbial methanogenesis. *Org. Geochem.* **98**, 105–117 (2016).

12. Ballentine, C. J. & Burnard, P. G. Production, release and transport of noble gases in the continental crust. *Rev. Mineral. Geochem.* **47**, 481–538 (2002).

13. LaForce, T., Ennis-King, J., Boreham, C. & Paterson, L. Residual CO₂ saturation estimate using noble gas tracers in a single-well field test: the CO₂CRC Otway project. *Int. J. Greenhouse Gas Control* **26**, 9–21 (2014).

14. Györe, D., Stuart, F. M., Gilfillan, S. M. V. & Waldron, S. Tracing injected CO₂ in the Cranfield enhanced oil recovery field (MS, USA) using He, Ne and Ar isotopes. *Int. J. Greenhouse Gas Control* **42**, 554–561 (2015).

15. Sathaye, K. J., Larson, T. E., Hesse, M. A. Noble gas fractionation during subsurface gas migration. *Earth Planet Sci. Lett.* **450**, 1–9 (2016).

16. Shelton, J. L. et al. Determining CO₂ storage potential during miscible CO₂ enhanced oil recovery: noble gas and stable isotope tracers. *Int. J. Greenhouse Gas Control* **51**, 239–253 (2016).

17. Barry, P. H. et al. Tracing enhanced oil recovery signatures in casing gases from the Lost Hills Oil Field using noble gases. *Earth Planet Sci. Lett.* **496**, 57–67 (2018).

18. Barry, P. H., Lawson, M., Meurer, W. P., Cheng, A. & Ballentine, C. J. Noble gases in deepwater oils of the U.S. Gulf of Mexico. *Geochem. Geophys. Geosyst.* **19**, 4218–4235 (2018).
19. Gilfillan, S. M. V. et al. The noble gas geochemistry of natural CO₂ gas reservoirs from the Colorado Plateau and Rocky Mountain provinces, USA. *Geochim. Cosmochim. Acta* **72**, 1174–1198 (2008).

20. Sherwood Lollar, B., Ballentine, C. J. & Onions, R. K. The fate of mantle-derived carbon in a continental sedimentary basin: integration of C-He relationships and stable isotope signatures. *Geochim. Cosmochim. Acta* **61**, 2295–2307 (1997).

21. Clark, I. D. & Fritz, P. *Environmental Isotopes in Hydrogeology* (CRC Press/Lewis Publishers, 1997).

22. Wycherley, H., Fleet, A. & Shaw, H. Some observations on the origins of large volumes of carbon dioxide accumulations in sedimentary basins. *Mar. Pet. Geol.* **16**, 489–494 (1999).

23. Bottinga, Y. Calculated fractionation factors for carbon and hydrogen isotope exchange in the system calcite–carbon dioxide–graphite–methane–hydrogen–water vapor. *Geochim. Cosmochim. Acta* **33**, 49–64 (1969).

24. Thiagarajan, N. et al. Isotopic evidence for quasi-equilibrium chemistry in thermally mature natural gases. *Proc. Natl Acad. Sci. USA* **117**, 3989–3995 (2020).

25. Stolper, D. A. et al. Formation temperatures of thermogenic and biogenic methane. *Science* **344**, 1500–1503 (2014).

26. Douglas, P. M. J. et al. Methane clumped isotopes: progress and potential for a new isotopic tracer. *Org. Geochem.* **113**, 262–282 (2017).

27. Stolper, D. A. et al. The utility of methane clumped isotopes to constrain the origins of methane in natural gas accumulations. *Geol. Soc. Lond. Spec. Publ.* **468**, 23–52 (2018).

28. Giunta, T. et al. Methane sources and sinks in continental sedimentary systems: new insights from paired clumped isotopologues $\Delta^{13}CH_3D$ and $\Delta^{13}CH_3D_2$. *Geochim. Cosmochim. Acta* **245**, 327–351 (2019).

29. Young, E. D. in *Deep Carbon* (eds Orcutt, B. N. et al.) 388–414 (Cambridge Univ. Press, 2019); https://doi.org/10.1017/9781108677950.013

30. Stolper, D. A. et al. Distinguishing and understanding thermogenic and biogenic sources of methane using multiply substituted isotopologues. *Geochim. Cosmochim. Acta* **161**, 219–247 (2015).

31. Wang, D. T. et al. Nonequilibrium clumped isotope signals in microbial methane. *Science* **348**, 428–431 (2015).

32. Gruen, D. S. et al. Experimental investigation on the controls of clumped isotopologue and hydrogen isotope ratios in microbial methane. *Geochim. Cosmochim. Acta* **237**, 339–356 (2018).

33. Cao, X., Bao, H. & Peng, Y. A kinetic model for isotopologue signatures of methane generated by biotic and abiotic CO₂ methanation. *Geochim. Cosmochim. Acta* **249**, 59–75 (2019).

34. Ono, S. et al. Clumped isotopolgue fractionation by microbial cultures performing the anaerobic oxidation of methane. *Geochim. Cosmochim. Acta* **293**, 70–85 (2021).

35. Sweeney, J. J. & Burnham, A. K. Evaluation of a simple model of vitrinite reflectance based on chemical kinetics. *AAPG Bull.* **74**, 1559–1570 (1990).

36. Head, I. M., Jones, D. M. & Larter, S. R. Biological activity in the deep subsurface and the origin of heavy oil. *Nature* **426**, 344–352 (2003).

37. Shelton, J. L. et al. Microbial community composition of a hydrocarbon reservoir 40 years after a CO₂ enhanced oil recovery flood. *FEMS Microbiol. Ecol.* **94**, fiy153 (2018).

38. Gray, N. D. et al. Biogenic methane production in formation waters from a large gas field in the North Sea. *Extremophiles* **13**, 511–519 (2009).

39. Perri, P. R. *Advanced Reservoir Characterization in the Antelope Shale to Establish the Viability of CO2 Enhanced Oil Recovery in California's Monterey Siliceous Shales* Annual Report (Chevron USA Production Company, 2003).

40. Gannon, R. S. et al. *Produced Water Chemistry Data for the Lost Hills and Fruitvale Study Areas, Southern San Joaquin Valley, California, 2016* (US Geological Survey, 2020); https://doi.org/10.5066/F7X929H9

41. Young, E. D. et al. The relative abundances of resolved $\Delta^{13}CH_3D$ and $\Delta^{13}CH_3D$ and mechanisms controlling isotopic bond ordering in abiotic and biotic methane gases. *Geochim. Cosmochim. Acta* **203**, 235–264 (2017).

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021
# Methods

## Modelling

Here we provide extra details of the δ¹³C-CO₂ fractionation modelling of methanogenesis and dissolution. In our models (Fig. 2a, Extended Data Fig. 6), we use the sample with the highest measured CO₂/³He (most pristine sample, O5) as a reference point to calculate the change in CO₂/³He and δ¹³C-CO₂ with CO₂ processing (that is, methanogenesis and dissolution). We assume there is no³He loss from the gas phase by dissolution or methanogenesis and therefore changes in the fraction of the CO₂ cause changes in the CO₂/³He.

Dissolution is modelled as open-system Rayleigh fractionation⁴.

$$ \delta^{13}C_{CO2} = \delta^{13}C_{CO2i} \times f^{(\alpha-1)}, \quad (2) $$

where f is the fraction of the original CO₂ remaining in the reservoir, δ¹³C<sub>CO2</sub> and δ¹³C<sub>CO2i</sub> are the carbon isotopic composition of the CO₂ for post-dissolution and in the initial CO₂ respectively, and α is the fractionation factor. The fractionation factor for dissolution (α<sub>d</sub> = 1.0038 ± 0.0012) was calculated based on the average measured pH¹⁰. Precipitation curves (Fig. 2a) can be modelled in the same way using fractionation factors (α<sub>p</sub> = 1.0086 ± 0.0012) calculated for the current reservoir temperature²¹.

From consideration of the carbon isotopes, which appear to be approaching equilibrium under current reservoir conditions, net methanogenesis has been modelled as a closed-system equilibrium process (equation (2), Fig. 2a). The fractionation resulting from this cycling, α<sub>m</sub>, is 0.9363 ± 0.003, which is consistent with thermal equilibration at current reservoir temperatures²³.

$$ \delta^{13}C_{CO2} = \delta^{13}C_{CO2i} \times f^{(\alpha-1)} - (\alpha_m - 1)(1-f) \times 1,000 \quad (3) $$

Within the Olla Oil Field, it is likely that there is a combination of both microbial cycling and dissolution. As a result, the net effect of microbial methanogenesis and dissolution can be described using equation (3).

$$ \delta^{13}C_{CO2} = \delta^{13}C_{CO2i} \times f^{(\alpha-1)} - (\alpha_m - 1)(1-f) \times 1,000 \quad (4) $$

where f is the fraction of CO₂ remaining after trapping/consumption and α<sub>A</sub> is the apparent fractionation factor for the system. α<sub>A</sub> depends on the degree of dissolution compared with methanogenesis and is described as⁴²:

$$ \alpha_A = F_d \alpha_d + (1 - F_d) \alpha_m \quad (5) $$

where F<sub>d</sub> is the fraction of CO₂ of the total CO₂ trapped/consumed due to dissolution.

## Sample collection and analysis

Samples for bulk gas composition, stable isotope, noble gas and clumped isotope analysis were collected directly from the wellhead using standard techniques following methods from refs.⁸,¹⁷,¹⁸,²⁴. Noble gas isotope determination was conducted in the Noble Laboratory at the University of Oxford, where the analysis of hydrocarbon gases is well established⁸,¹⁷,¹⁸,⁴³⁻⁴⁵. For bulk composition and C and H stable isotope analysis, gas cylinders were shipped to Isotech in Champaign, Illinois, USA. The standard procedures for these are described in ref.⁴⁶. CH₄ clumped isotopologues were determined at the California Institute of Technology following the procedures described in ref.⁴⁶.

Microbiological analysis was conducted at an internal ExxonMobil Research and Engineering Company Facility. The produced water samples were collected and filtered immediately from wells O1, O4, NHH3 and NHH6. Water was passed through 0.22-µm SVGPL10 RC Sterivex GP filters (EMD Milipore) until filters clogged in the following amounts: O1, 1,350 ml; O4, 1,350 ml; NHH3, 1,750 ml; NHH4, 300 ml. Filters were stored

on ice immediately and stored at -80 °C until DNA extraction. DNA was extracted from the Sterivex filters using the DNeasy PowerWater Sterivex Kit (Qiagen). The SSU rRNA gene was amplified with V4V5 specific primers (https://vamps2.mbl.edu/) targeting Bacteria. The SSU rRNA gene was amplified with specific primers targeting Bacteria: forward primers (967F) CTAACCGANGAACCTYACC, CNACGCGAAG AACCTTANC, CAACGCGMARAACCTTACC and ATACGCGARGAACCTT ACC; reverse primer (1064R) CGACRCCATGCANCACCT. A second set of amplifications was performed using specific primers targeting Archaea: (Arch2A519F) CAGCMGCCGCGGTAA and (Arch1017R) GGCC ATGCACCWCTCTC. Illumina MiSeq 2x300bp paired-end sequencing was performed by Mr. DNA. Microbiome bioinformatics were performed with QIIME 2.2017.4⁴⁷. Raw sequence data were demultiplexed and quality filtered using the q2-demux plugin followed by denoising with DADA2⁴⁸. Amplicon sequence variants were aligned with mafft and used to construct a phylogeny with fasttree²⁴⁹,⁵⁰. Taxonomy was assigned to amplicon sequence variants using the q2-feature-classifier⁵¹ classify-sklearn Naïve Bayes taxonomy classifier against the Greengenes 13_8 99% operational taxonomic unit (OTU) reference sequences⁵².

## Geological and production history

The Olla Oil Field is located in La Salle Parish, Louisiana, in the northern Gulf of Mexico and was CO₂ injected between 1983 and 1986 for enhanced oil recovery (CO₂-EOR). The Nebo-Hemphill Oil Field lies 20 km to the southeast of the Olla Oil Field and provides a control field for this study, having had no CO₂ injected (Fig. 1). Both fields are producing from the Palaeocene-Eocene Wilcox Group, which in northern Louisiana is composed of 600-1,500 m of shallow marine clastic sediments to fluvial deltaic deposits and up to 6-m-thick coal beds⁵³,⁵⁴. The group is bound either side by the Palaeocene Midway Group and the Eocene Clairborne Group confining units⁵⁴. Typically, the Wilcox Group is subdivided into the lower (sandstone with coal beds), middle (sandstone, lignite and clay) and upper (deltaic deposits) intervals, with the 'Big Shale' informal unit separating the middle and upper groups⁵³,⁵⁵. Late Cretaceous intrusive igneous bodies in northern Louisiana, west-central Mississippi and southeastern Arkansas are thought to have affected the regional heat flow, which may have caused an increase in the thermal maturity of the Wilcox Group during the early Palaeocene⁵⁶,⁵⁷.

Gas/oil ratios are variable within the field, and do not seem to be associated with proximity to the known gas cap and higher gas/oil ratios may instead reflect gas coming out of solution as bottom hole pressure decreases.

CO₂ is present in low concentrations (1.02 ± 0.69 mol%) within the Nebo-Hemphill Oil Field, which suggests that there is some natural background CO₂ present in both fields. Between March 1983 and April 1986, the 2,800-ft sandstone (within the middle Wilcox Group) in the Olla Oil Field underwent a pilot project for CO₂ flooding. The CO₂ injected into the field was immiscible and injected in the gas phase (that is, not supercritical). It is assumed that a negligible amount of CH₄ was injected alongside the CO₂. Previous studies on the Olla Oil Field have assumed that all the injected CO₂ was retained within the 2,800-ft sandstone⁹,¹⁰. However, the initial CO₂ injection report⁷ states that "mole percentages in excess of 10% more than likely indicate production of injected CO₂". Only 1 sample outside of the 2,800-ft sandstone has CO₂ concentrations less than 10% (that is, all other CO₂ concentrations are greater than 10%). Notably, concentrations greater than 10% in other producing zones have previously been reported¹⁰. These signatures are consistent with injected CO₂ migration across the different reservoirs and remaining within the formations upon sampling. If such high concentrations of CO₂ were present in this field before injection, CO₂ would not have been extracted and transported from the Black Lake Oil Field and would have instead have been extracted from the producing fluids of CO₂ rich fluids at Olla. In addition, we see a clear mantle-derived noble gas signature (see below),
which is most evident in the He and neon (Ne) isotopes (Extended Data Fig. 1), in the Olla Oil Field, which shows a clear incursion of injection CO₂ into the proximal lithologies.

A total of 2.2 × 10⁸ m³ of CO₂ was injected into the Olla Oil Field and was sourced from the nearby Black Lake Oil Field (produced from the Lower Cretaceous carbonates of the Petit Formation) in Natchitoches Parish, Louisiana (Fig. 1). Approximately one-third of the injected CO₂ remained within the reservoir after the cessation of the EOR project⁷. It is important to note that given the age of this CO₂-EOR project, we were unable to obtain constraints such as the ratio of injected CO₂ to reservoir pore volume that are sometimes cited to quantify the scale of the EOR project and the distribution it achieves within the subsurface. Further, we anticipate that the injected CO₂ is not uniformly distributed within the pore volume of the Olla Oil Field due to different phases of EOR and gas cycling.

In 1984, water flooding also started within the field. In contrast, Nebo-Hemphill has never undergone any type of EOR.

## Molecular geochemistry

Whole-oil gas chromatography (WOGC) was performed to provide a hydrocarbon fingerprint of each oil, and open-column liquid chromatography was utilized to separate and quantify the saturate, aromatic, resins and asphaltene fractions of each oil. In general, we see similar fluid compositions between the two fields. The WOGC shows evidence of marine-sourced oils with a smooth distribution of *n*-alkanes, pristine/phytane ratios of approximately 2 and low wax content. The notable differences are probably due to varying levels of biodegradation illustrated by the loss of *n*-alkanes in some samples and American Petroleum Institute (API) gravity ranging from 21 to 42. The similarities in the geochemical signatures of these oils provide a basis to use the Nebo-Hemphill Oil Field as a control site for any biogeochemical differences that may arise as a result of the injection of CO₂ at Olla.

## Source facies

Gas chromatography mass spectroscopy (GC/MS) was utilized to analyse the saturate and aromatic fractions to provide the biomarker distributions of each oil. The interpretation of a dominantly marine source is based on the following evidence: (1) a smooth distribution of tricyclic terpanes with C₂₃ being the maximum along with extended hopanes (C₃₁–C₃₅) decreasing in order is a classic marine source signature; (2) the high presence of C₃₀ hopane compared with C₂₉ hopane is indicative of a marine clastic source; (3) a higher abundance of C₂₈ bisnorhopane compared with oleanane suggest a higher marine input compared with the above source⁵⁸; (4) a greater abundance of C₂₆ tricyclic terpane as compared to C₂₅ tricyclic terpane provides evidence of some terrigenous input. Consistent with previous interpretations, we consider this source to be a part of the Wilcox Group oils, more specifically Wilcox Group I as originally described by refs.⁵⁹,⁶⁰.

## Thermal maturity

We adopted standard saturate and aromatic biomarker ratios to develop apparent thermal maturity estimates for each of the oil samples analysed. We discuss these in the context of apparent thermal maturities because we recognize these estimates are based on the bulk weighted average composition of the fluid. The interpretation of the GC/MS analyses of the 12 oils yielded a thermal maturity range of R₀ equivalent of 0.9–1.25%, which is classified as main-to-late-stage oil. The interpretation of this thermal maturity range is based on the following: (1) the presence of higher concentrations of high-molecular-weight compared with low-molecular-weight triaromatic steriods, described as % ΔTAS (20–57%); (2) the isomerization of C₂₉ ααα-20S relative to C₂₉ ααα-20R steranes (0.45–0.48)⁶¹; (3) the isomerization of 1,2,5-trimethylnaphthalene to 1,3,7-trimethylnaphthalene represented as ratio TMNr (0.48–0.67) as defined by ref.⁶²; (4) the ratio of C₂₃ tricyclic terpane to C₃₀ hopane, which is calculated as T23/(T23 + H30) (0.10–0.37), where the increased presence of C₂₃ tricyclic represents higher maturity.

## Alteration

The use of WOGC and GC/MS analyses help provide a picture of the alteration history of these oils. The shallow Nebo-Hemphill NHH3 fluid is the most biodegraded end member of the samples tested, based on the apparent lack of *n*-alkanes, an API gravity of 21 and the presence of demethylated hopanes. This sample would be classified as severely biodegraded based on the classification of ref. ⁶³. The remaining Nebo-Hemphill samples have APIs of 32–41 and are classified as unbiodegraded to very slight biodegradation⁶³. The Olla Oil Field fluids show evidence of reduced light *n*-alkanes, display similar *n*-alkane-to-isoprenoid ratios and have APIs of 33 to 37, suggesting similar levels of biodegradation between the fluids. Therefore, the Olla fluids would be classified as having very slight biodegradation⁶³.

## Summary of previous findings

The formation water geochemistry, molecular oil composition and produced gas composition (including δ¹³C) in the Olla and Nebo-Hemphill oil fields have been previously characterized⁹⁻¹¹. Microbial methanogenesis was first hypothesised for the Olla Oil Field based on its enriched δ¹³C of CO₂ (ref. ⁹). However, owing to the differences in the isotopic composition between the 2,800-ft sandstone and the other producing formations with the injected material, they concluded that the injected CO₂ was not present in the producing zones¹⁰ and consequently any microbial methanogenesis was independent of injection and is probably associated with microbial oxidation of hydrocarbons. The enhanced methanogenesis within the Olla Oil Field compared with other fields produced from the Wilcox Group was thought to be due to the ideal geochemical conditions (that is, suitable range of pH, temperature and salinity, and lack of alternative electron acceptors)¹⁰. They conclude that the injected CO₂ was most likely dissolved and trapped in a gas phase¹⁰. Further investigation of microbial methanogenesis found that there was little correlation between crude oil biodegradation and methanogenesis; however, the extent of oil biodegradation was correlated with the temperature, salinity and depth of samples¹¹. In addition, there has also been an investigation of the microbial community composition within the Olla reservoir³⁷, which found methanogens present within the system, and no significant difference in microbial communities between the 2,800-ft sandstone and other producing zones³⁷. With the addition of the inert noble gas isotopes and clumped isotopes, we are able to re-investigate the fate of CO₂ within the Olla reservoirs.

## Formation water chemistry in the Olla Oil Field

The pH and alkalinity of the formation waters within the Olla Oil Field are between 6.5 and 7.5 and between 16.3 meq kg⁻¹ and 57.7 meq kg⁻¹, respectively¹⁰. The most dominant anion was Cl⁻ (1,122–1,621 mmol l⁻¹) and the most dominant cation measured was Na⁺ (1,108–1,630 mmol l⁻¹) (Extended Data Table 1), which are typical of basinal brines⁶⁴ and in agreement with previous measurements¹⁰. Other cations present include: Mg²⁺ (119–329 ppm), Ca²⁺ (80–500 ppm), K⁺ (223–524 ppm), Ba²⁺ (95–165 ppm), Si⁴⁺ (10.1–12.5 ppm) and Sr²⁺ (77–215 ppm). Measured iron concentrations are less than 0.1 ppm. With the exception of sample O1 where SO₄²⁻ concentrations are 62 ppm, SO₄²⁻ concentrations are less than 2 ppm. Both the iron and SO₄²⁻ concentrations are in agreement with those previously reported¹⁰.

## Identification of injected material

CO₂ concentrations in the injected Olla Oil Field (18 ± 12 mol%) are significantly greater than the Nebo-Hemphill Oil Field (1.02 ± 0.69 mol%). CH₄ concentrations at Olla are 77 ± 13 mol% (STP) compared with 92.6 ± 8.2 mol% (STP) at Nebo-Hemphill. However, the δ¹³C of CO₂ of the injected material (0.85 ± 0.91‰) is isotopically distinct from the δ¹³C of CO₂ now within Olla (13.6 ± 3.2‰) (Fig. 2b). The injected CO₂ was sourced from the Black Lake Oil Field, which is located close to the
Sabine Island Complex. The Sabine Island Complex forms part of the Sabine Uplift basement high and has been suggested to be a conduit for mantle-derived fluids⁸. However, the measured δ¹³C of CO₂ values at Black Lake (0.85 ± 0.91‰)¹⁰ are elevated relative to mantle CO₂ (−2‰ to −5%; refs. ²¹,²²), probably reflecting some degree of biological modification when the reservoir was shallower and at temperatures favourable to microbial activity. An alternative interpretation for the elevated δ¹³C of CO₂ at Black Lake is that it is a mixture of mantle-derived and crustal-derived CO₂; however, this would increase the CO₂/³He above that expected in the mantle, which is not observed in the Olla samples.

As the stable isotopes of carbon are frequently modified by chemical or biological processes, the inert nature of noble gases means that they are ideal tracers for injected CO₂ in such scenarios. As injectate samples were not collected during CO₂-EOR, we do not have access to the noble gas compositions of the material injected into the Olla Oil Field. Nevertheless, it is reasonable that the noble gases record pre-injection conditions at the nearby Nebo-Hemphill Oil Field, and that these are subsequently modified at the Olla Oil Field through mixing of the injected material during the CO₂ injection.

The ⁴He/²⁰Ne ratio is between 3,210 and 20,900 in both fields, which is significantly greater than the atmospheric value of 0.32 and thus there is no significant atmospheric contribution to these samples. The air-corrected He isotopic ratio (³He/⁴He), reported relative to the atmospheric ratio (R<sub>C</sub>/R<sub>A</sub> where air is 1R<sub>A</sub>), in the Olla Oil Field is 1.84 ± 0.10R<sub>A</sub> and 0.46 ± 0.08R<sub>A</sub> at Nebo-Hemphill (Extended Data Table 2). These are both higher than a typical hydrocarbon system where the majority of He is derived from crustal radiogenic production (³He/⁴He = 0.02R<sub>A</sub>; ref. ¹²), and higher than other fields within the Gulf of Mexico region⁸,¹⁸,⁴³. The elevated ³He/⁴He values are the result of a resolvable mantle contribution. Using a simple two endmember mixing model, with subcontinental lithospheric mantle (³He/⁴He = 6.1; ref. ⁶⁵) and crustal endmembers, we resolve an average 7.8% mantle contribution of mantle He in Nebo-Hemphill Oil Field. It is likely that the mantle fluid contribution is coming from the Little Creek collapse structure, which is situated between the two fields⁶⁶. At Olla, the average mantle contribution to ³He is 30%, which suggests that there is an additional source of mantle fluids to this field, which is reasonably explained by the injectate sourced from the Black Lake Oil Field. Ne isotope ratios are elevated at Olla (²⁰Ne/²²Ne 10.3 ± 0.2 and ²¹Ne/²²Ne 0.0399 ± 0.002) compared with Nebo-Hemphill (9.85 ± 0.09 and 0.0335 ± 0.001, respectively) and air (9.8 and 0.029) (Extended Data Fig. 1, Extended Data Table 3)⁶⁷,⁶⁸,⁶⁹, again suggesting that the injected material contains mantle-derived fluids.

The ⁴⁰Ar/³⁶Ar values at Nebo-Hemphill are significantly greater than the atmospheric value (298.56; ref. ⁷⁰), with measured ratios of 370 ± 3 to 604 ± 5. The Olla Oil Field has a higher ⁴⁰Ar/³⁶Ar (1,476 ± 13 to 2,766 ± 24) (Extended Data Table 3). Elevated ⁴⁰Ar/³⁶Ar ratios can be a result of both an excess of radiogenic ⁴⁰Ar and mantle-derived fluids¹²,⁷¹. The difference in the ⁴⁰Ar/³⁶Ar ratios could be due to variable amounts of radiogenic production or fundamental differences in the extent of contact with formation waters between the two fields. However, as all the wells sampled are within the same geological formation, the two fields are in relatively close proximity and they have indistinguishable ⁴He concentrations, it seems unlikely that this would be the cause. Thus, the difference is probably a result of the injected CO₂, which probably had mantle-derived ⁴⁰Ar/³⁶Ar excesses, concordant with the observed mantle-derived He and Ne isotopes. Finally, whereas xenon (Xe) isotope ratios at Nebo-Hemphill are indistinguishable from air, the Olla Oil Field showed excesses in some Xe isotope ratios relative to air. For example, the ¹²⁹Xe/¹³⁰Xe ratio is 6.59 ± 0.01 compared with an atmospheric value of 6.49 (Extended Data Table 3); this excess in ¹²⁹Xe relative to ¹³⁰Xe is consistent with a contribution from mantle fluids with an end member ¹²⁹Xe/¹³⁰Xe value of 8.2 (ref. ⁷¹). The Olla Oil Field therefore shows a clear mantle component that is not present within the Nebo-Hemphill control field.

The detection of injected fluids within the Olla Oil Field using the conservative noble gas tracers differs from the findings of the previous studies that used the δ¹³C of CO₂ to conclude that injected CO₂ is no longer or never was present in the majority of the sampled wells (for example, ref. ¹⁰). The subsurface behaviour of the injected CO₂ within the system can therefore be re-examined using a combination of noble gases and stable isotopes.

## Dissolution within the Olla Oil Field

CO₂ is highly soluble in water under reservoir conditions. The most likely trapping mechanisms for free-phase subsurface CO₂ are dissolution or precipitation as carbonate minerals mediated by dissolution⁴,²¹. A decrease in the CO₂/³He correlates with increasing ⁴He and ²⁰Ne concentrations within the Olla Oil Field (Extended Data Fig. 2). ²⁰Ne (and often ⁴He)⁴ are derived from the formation water, suggesting that CO₂ trapping in the Olla Oil Field is also proportional to the amount of formation water that has been contacted and degassed, consistent with dissolution within the field⁴. Dissolution (solubility trapping) has previously been found to be a major sink for CO₂ (ref. ⁴) and imparts a slight pH-dependent fractionation towards a light isotopic composition²¹. However, this fractionation effect is marginal compared with the amount of CO₂ that can be trapped and the fractionation that can be imparted by methanogenesis-AOM cycling (Fig. 2a).

Given the presence of discontinuous coal beds within the Wilcox Group⁵⁷,⁷², it is possible that some of the injected CO₂ may have been sorbed onto these coal beds as CO₂ has a much greater affinity to coal than CH₄ (ref. ⁹). The resulting isotopic fractionation imparted by this process remains poorly constrained, with some studies suggesting a greater than 5‰ shift⁷³,⁷⁴ whereas others suggest negligible effects⁷⁵⁻⁷⁷ on the δ¹³C of CO₂. Coal sorption would cause fractionation in the opposition direction of the CH₄ isotopic signature. Given the heavier δ¹³C of CO₂ and δ¹³C of CH₄ isotopic signatures in the injected Olla Oil Field, it is unlikely that CO₂ sorption to coal is a major CO₂ sink in the field, in agreement with previous findings¹⁰.

## Subsurface behaviour of carbon at Olla

The Olla samples are enriched in both the δ¹³C of CO₂ and the δ¹³C of CH₄ compared with Nebo-Hemphill and the injected CO₂ values. The offset between the δ¹³C of CH₄ and the δ¹³C of CO₂ in the Olla Oil Field are close to those expected under thermally driven isotope exchange equilibrium for the observed temperatures (Fig. 2b). Carbon isotopes of CO₂ and CH₄ approaching equilibrium in hydrocarbon systems with current reservoir conditions have previously been observed²⁴. However, the thermal equilibrium of CO₂ and CH₄ is kinetically limited at the current temperatures (29.2–50.7 °C) in the Olla reservoirs and would not occur on the timescale of injection, without a process to enhance the reaction rate⁷⁸. A combination of both microbial methanogenesis and AOM active within the field would allow for the system to be approaching equilibrium under the current thermal conditions within the reservoir. Further observations detailed below provide evidence for this conclusion and allow for the net microbial methanogenesis to be quantified.

The clumped isotopologues of CH₄ (Δ¹³CH₃D and Δ¹²CH₂D₂) can provide independent constraints on its origin within petroleum systems²⁵⁻²⁹. Thermogenic CH₄ is expected to have a composition consistent with the source-rock maturation temperature, which is typically at or above 157 °C (ref. ²⁵). At Olla, we expect CH₄ generation to have occurred at 163 ± 18 °C (based on biomarker maturity estimates following ref. ³⁵). Microbial CH₄ is variable in its clumped isotopic compositions with both non-equilibrium²⁸,²⁹ and apparent equilibrium signatures reported in samples obtained from the deep subsurface²⁶,²⁹. Microbial methanogenesis forms CH₄ characterized by low Δ¹³CH₃D and extremely low Δ¹²CH₂D₂ (refs. ²⁸,²⁹), whereas AOM tends to drive the isotopologue composition of CH₄ towards equilibrium, with equilibrium expected to be reached at a faster rate in Δ¹³CH₃D compared
with in $\Delta^{12}\text{CH}_2\text{D}_2$ (refs. 29,30,34,79-81). In the Olla Oil Field, we see $\Delta^{13}\text{CH}_3\text{D}$ of 3.45-5.62‰ (apparent temperature of $29^{+14}_{-12}^\circ\text{C}$ to $128^{+17}_{-16}^\circ\text{C}$) and $\Delta^{12}\text{CH}_2\text{D}_2$ of 9.13-12.4‰ (apparent temperature of $87^{+13}_{-11}^{+13}_{-12}^{+30}_{-23}^\circ\text{C}$) (Fig. 3, Extended Data Table 2); these temperatures preclude a thermogenic origin for the CH$_4$ and are instead consistent with a microbial contribution. Samples O4 and O6 are in internal clumped isotopologue equilibrium at temperatures lower than expected for the thermogenic CH$_4$ and appear to be approaching equilibrium with current reservoir conditions. The remaining samples from Olla have already reached apparent equilibrium with the current reservoir conditions in $\Delta^{13}\text{CH}_3\text{D}$ but have yet to reach apparent equilibrium in $\Delta^{12}\text{CH}_2\text{D}_2$ probably due to kinetic effects. The clumped isotopologues at Olla are approaching equilibrium at reservoir temperatures, which would be accounted for by active AOM within the Olla Oil Field. This supports the conclusion that there is active microbial methanogenesis and AOM driving the carbon isotopes to approach equilibrium under current thermal conditions in the Olla reservoirs.

There is microbial degradation of the higher alkane gases (Extended Data Fig. 3) within the Olla Oil Field. Propane is preferentially consumed during microbial degradation and thus in biodegraded gases, propane will see the most pronounced $\delta^{13}\text{C}$ enrichment$^{36}$. This $\delta^{13}\text{C}$ enrichment is clear within the Olla gases (Extended Data Table 3) and supports active microbial hydrocarbon degradation within the field.

The presence of methanogenic activity and the idea of a methanogenic 'hotspot' in the Olla reservoir has been previously suggested$^{9,10}$. The environmental conditions for which methanogenesis can occur are relatively well defined. Within nutrient-poor environments such as oil reservoirs or coal beds, methanogens can survive at temperatures up to 80 °C, with the largest growth rates and thus greatest amount of methanogenic activity between 40 °C and 50 °C (refs. 36,82-84). Methanogenic activity is greatest at near neutral pH and is inhibited at pH <4 and >9 (ref. 84). Formation water salinity can also inhibit methanogenesis above $\text{Cl}^-$ concentrations of 1,500 mM, with greatest growth between 200 nM and 600 mM (refs. 82,83,85). The presence of ferric iron and sulfate at concentrations greater than 1 mM can also limit methanogenesis (as a result of the dominance of other bacteria)$^{86-89}$. Notably, the Olla Oil Field exhibits ideal environmental conditions for methanogenesis, in terms of pH (6.5-7.5)$^{10}$, temperature (29.2-50.7 °C), salinity (1,094-1,622 mM), sulfate (<0.6 mM $\text{I}^{-1}$)$^{10}$ and iron concentrations (≤0.1 mM $\text{I}^{-1}$) (Methods).

In addition, geochemical and isotopic parameters associated with microbial methanogenesis in organic-rich settings are also well established$^{10,11,83,88,90}$. Traditionally, carbon isotopes of CH$_4$ have been used to determine the presence of biogenic CH$_4$ (typically less than -60‰)$^{91-93}$ (Extended Data Fig. 3a). However, methanogenesis can result in heavy biogenic CH$_4$ as the residual CO$_2$ pool gets increasingly elevated, consistent with observations at Olla (Fig. 2b). In addition, where methanogenesis has occurred, formation water alkalinity is greater than 10 meq kg$^{-1}$, the acetate concentrations are lower than 1 mM, the Ca/Mg ratios are less than 1.5 and the $\delta^{13}\text{C}$-DIC is greater than +20‰ (refs. 87,88,90). Within the Olla Oil Field, the alkalinity ranges between 21.4 meq kg$^{-1}$ and 57.7 meq kg$^{-1}$, the average Ca/Mg ratio is 0.2 and the $\delta^{13}\text{C}$-DIC is 14.2-28.7‰, again indicating that methanogenesis has occurred within the field.

SSU rRNA gene-sequencing analysis of the microbial communities within the Olla Oil Field (samples O1 and O4) indicate that within the archaeal sequences, hydrogenotrophic methanogens, methanol and methylotrophic methanogens, and ANME assigned sequences are present (Extended Data Fig. 5). This microbial array would facilitate both microbial methanogenesis and AOM within the field.

During the cycling of carbon between methanogenesis and AOM, we would not expect a change in the isotopic composition of hydrogen. It is possible that as a result of microbial activity or the reaction kinetics that extracellular hydrogen is added to the system. This hydrogen could be sourced from the water and free hydrogen from the degradation of the higher hydrocarbon gases or from the oil. As a result, we are unable

to determine whether the hydrogen isotopes within the system are also approaching equilibrium under current reservoir temperatures.

It is possible that the CO$_2$ and CH$_4$ isotopes in the Olla system are reflecting purely the kinetic process of methanogenesis. The difference in the $\delta^{13}\text{C}$ of CO$_2$ and the $\delta^{13}\text{C}$ of CH$_4$ at Olla (53.7-64.4‰) is consistent with that expected from microbial methanogenesis (about 58‰; ref. 94). In such a scenario, all the CH$_4$ at Olla would be microbial in origin, thus requiring any thermogenic CH$_4$ component originally in the system to be overwhelmed. The high concentration and mixed thermogenic-microbial composition of CH$_4$ within the Nebo-Hemphill Oil Field (Fig. 3, Methods) is inconsistent with such a scenario. In addition, the clumped isotopologues at Olla show a clear need for AOM within the field and 16S rRNA gene sequencing shows active methanotrophs, which is not compatible with a scenario of purely microbial methanogenesis dominating the system. Therefore, future studies investigating the processes affecting CO$_2$ in potential CCS storage targets will benefit from a complete geological, geochemical and microbiological characterization of the field study site.

## Subsurface behaviour of carbon in the Nebo-Hemphill Oil Field

The carbon isotopic composition of CO$_2$ and CH$_4$ at Nebo-Hemphill is depleted compared with the Olla Oil Field (Fig. 2b). Despite this, the difference in the $\delta^{13}\text{C}$ of CO$_2$ and the $\delta^{13}\text{C}$ of CH$_4$ at Nebo-Hemphill overlaps with that of Olla (50.8-68.1‰ and 53.7-64.4‰, respectively). The degree of isotopic fractionation at Nebo-Hemphill is consistent with thermal equilibrium at 29.8-85.3 °C, in agreement with present-day reservoir temperatures (29.7-57.1 °C). In addition, the $\delta^{13}\text{C}$ of CO$_2$ at Nebo-Hemphill is enriched compared with that typically expected for crustal or mantle-sourced CO$_2$ (less than 2‰)$^{21,22}$, consistent with a small amount of microbial methanogenesis within the field. We propose that there is also microbial methanogenesis (in agreement with previous studies$^{10}$) and AOM within the Nebo-Hemphill Oil Field; however, unlike at Olla, this has not been enhanced by any injection to the field.

Microbiological analysis of NHH3 and NHH6 show that hydrogenotrophic methanogens, methanol and methylotrophic methanogens, and ANME assigned sequences are present within Nebo-Hemphill, in agreement with both active methanogenesis and AOM within the field (Extended Data Fig. 5). Clumped isotopologues of CH$_4$ within the Nebo-Hemphill Oil Field give apparent temperatures for $\Delta^{13}\text{CH}_3\text{D}$ of $78^{+14}_{-12}^{+166}_{-29}^\circ\text{C}$ and apparent temperatures for $\Delta^{12}\text{CH}_2\text{D}_2$ that are greater than expected for thermogenic CH$_4$ in the field (163 ± 18 °C, Fig. 3). They suggest that the CH$_4$ is likely to be dominantly thermogenic with some microbial contributions. Samples NHH1, NHH2 and NHH8 have apparent temperatures below that expected for thermogenic CH$_4$ and appear to be approaching low-temperature equilibrium in $\Delta^{13}\text{CH}_3\text{D}$, consistent with active methanogenesis and AOM cycling within the field.

## Modelling of the Pannonian Basin

The unequivocal identification and quantification of CO$_2$-related microbial methanogenesis at Olla has globally significant implications, and warrants the reconsideration of reservoir processes at other sites worldwide. Re-examination of previously published datasets allows the elucidation of the significance of methanogenesis in different reservoir environments, especially in systems that have or had naturally high CO$_2$ concentrations, such as the Pannonian Basin. A significant proportion of the CO$_2$ found in natural gas fields within the Pannonian Basin is of mantle origin$^{20}$. The initial, $\text{CO}_2/\text{He}$ in the Szegholm North, Szegholm South, Ebes and Hajduszoboszlo gas fields was estimated to be $7.9 \pm 5.4 \times 10^9$, consistent with the European subcontinental lithospheric mantle range ($0.6-40 \times 10^9$) and the mid-ocean ridge basalt mantle range ($2 \pm 1 \times 10^9$)$^{95-97}$ (Extended Data Fig. 7a). Since emplacement, significant CO$_2$ trapping in many fields was identified from lower observed $\text{CO}_2/\text{He}$ ratios (for example, the Szegholm North and South, Ebes and Hajduszoboszlo gas fields) and speculation that some of the
CO₂ may have been converted into CH₄ (ref. ²⁰), although no viable mechanism was identified.

Using an average geothermal gradient of 45 °C km⁻¹ the temperature within these fields is predicted to be between 40 °C and 105 °C (with temperatures being up to 20 °C cooler at the time of CO₂ emplacement)⁹⁸, and thus within the range of microbially driven methanogenesis. In addition, the presence of thermogenic hydrocarbons would also ensure H₂ availability, similar to the Olla and Nebo-Hemphill oil fields. From our work in Olla, we surmise that the Pannonian Basin may be analogous to Olla and that hydrogenotrophic microbial methanogenesis active within the field may provide the missing mechanism for CO₂ removal, which was previously unidentified, and point to this mechanism as an important subsurface process operating on systems naturally enriched in CO₂ as well as anthropogenically CO₂ injected oil and gas fields.

In the absence of supporting clumped isotopologue data, if we assume that the system is methanogenesis driven we can nevertheless construct similar methanogenesis–dissolution models using published data from the Szegholm South Gas Field (where multiple CO₂ isotope data are available). By comparing the observed CO₂/³He ratios against the CO₂ concentration, a clear correlation emerges (Extended Data Fig. 7a). Attributing this to the effect of dissolution (a pre-requisite for microbial methanogenesis), it is possible to extrapolate back to an initial CO₂/³He end member value of 13.3 × 10⁹, consistent with previous estimates.

Our model assumes mantle-derived CO₂ in the field with a typical isotopic composition of δ¹³C of CO₂ = -5‰ (Extended Data Fig. 7b). Taking this as a starting point, 16.4 ± 0.8 times more dissolution than methanogenesis is required in the Szegholm South Gas Field to produce the measured δ¹³C of CO₂. As such, proportionally less methanogenesis, relative to the Olla Oil Field, would be expected here as the gas fields are not being constantly perturbed via water injection and therefore nutrient availability is probably much more limited. The presence of biogenic CH₄ (through CO₂ conversion) can also be clearly identified from the measured C₁/Cₙ in the Szegholm South Gas Field (0.946 ± 0.004), which is greater than what would be expected in a purely thermogenic system (0.909). Using the excess in the C₁/Cₙ ratios, we determine that 41 ± 5% of the CH₄ currently present in the Szegholm South Gas Field part of the basin is biogenic. Similar C₁/Cₙ excesses are seen in the other gas fields in the Pannonian Basin⁹⁹. In addition, other fields such as Szegholm North²⁰,⁹⁹, Ebes²⁰,⁹⁹, Hajduszoboszło²⁰,⁹⁹, Kismarja²⁰,⁹⁹, Répcelak¹⁰⁰ and Mihàly¹⁰⁰ (Pannonian Basin), Subei¹⁰¹ (China) and JM Brown Bassett⁹⁵ (Permian Basin) also show geochemical signatures consistent with the occurrence of methanogenesis. Thus, we conclude that although methanogenesis is not the dominant CO₂ loss mechanism, it may play a substantial role that has not previously been considered in many natural systems.

## Data availability

The geochemical data that support the findings of this study are available in the NERC EDS National Geoscience Data Centre at https://doi.org/10.5285/a4070f5d-2064-4caf-a82c-79a786d6af9e. The microbial SRA and biosample data can be found at https://www.ncbi.nlm.nih.gov/bioproject/PRJNA744568. Source data are provided with this paper.

42. Van Breukelen, B. M. Extending the Rayleigh equation to allow competing isotope fractionating pathways to improve quantification of biodegradation. *Environ. Sci. Technol.* **41**, 4004–4010 (2007).

43. Byrne, D. J., Barry, P. H., Lawson, M. & Ballentine, C. J. Determining gas expulsion vs retention during hydrocarbon generation in the Eagle Ford Shale using noble gases. *Geochim. Cosmochim. Acta* **241**, 240–254 (2018).

44. Barry, P. H. et al. Noble gases solubility models of hydrocarbon charge mechanism in the Sleipner Vest gas field. *Geochim. Cosmochim. Acta* **194**, 291–309 (2016).

45. Tyne, R. L. et al. A novel method for the extraction, purification, and characterization of noble gases in produced fluids. *Geochem. Geophys. Geosyst.* **20**, 5588–5597 (2019).

46. Thiagarajan, N. et al. Identifying thermogenic and microbial methane in deep water Gulf of Mexico reservoirs. *Geochim. Cosmochim. Acta* **275**, 188–208 (2020).

47. Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIME 2. Nat. Biotechnol. **37**, 852–857 (2019).

48. Callahan, B. J. et al. DADA2: high-resolution sample inference from Illumina amplicon data. *Nat. Methods* **13**, 581–583 (2016).

49. Katoh, K., Misawa, K., Kuma, K. & Miyata, T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. *Nucleic Acids Res.* **30**, 3059–3066 (2002).

50. Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2—approximately maximum-likelihood trees for large alignments. *PLoS ONE* **5**, e9490 (2010).

51. Bokulich, N. A. et al. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2's q2-feature-classifier plugin. *Microbiome* **6**, 90 (2018).

52. McDonald, D. et al. An improved Greenoxes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. *ISME J.* **6**, 610–618 (2012).

53. Tye, R. S., Moslow, T. F., Kimbrell, W. C. & Wheeler, C. W. Lithostratigraphy and production characteristics of the Wilcox Group (Paleocene–Eocene) in central Louisiana. *AAPG Bull.* **75**, 1675–1713 (1991).

54. Johnston, J. E. III, Heinrich, P. V., Lovelace, J. K., McCulloh, R. P. & Zimmerman, R. K. *Stratigraphic Charts of Louisiana* (2000).

55. Echols, D. J. & Malkin, D. S. Wilcox (Eocene) stratigraphy, a key to production. *AAPG Bull.* **32**, 11–33 (1948).

56. Heydari, E., Byerly, G. R. & Henry, D. J. Contact metamorphism and over maturation of organic matter associated with an igneous intrusion in the Smackover Formation, Northeastern Louisiana. *Gulf Coast Assoc. Geol. Soc. Trans.* **47**, 201–213 (1997).

57. Warwick, P. D., Breland, F. C. & Hackley, P. C. Biogenic origin of coalbed gas in the northern Gulf of Mexico Coastal Plain, U.S.A. *Int. J. Coal Geology* **76**, 119–137 (2008).

58. Wenger, L. M., Sassen, R. & Shurmacher, D. S. in *Gulf Coast Oils and Gases: Their Characteristics, Origin, Distribution, and Exploration and Production Significance* (eds Schumacher, D. & Perkins, B. F.) Ch. 6 (GCSSEPM, 1990).

59. Sassen, R. Lower Tertiary and Upper Cretaceous source rocks in Louisiana and Mississippi: implications to Gulf of Mexico crude oil. *AAPG Bull.* **74**, 857–878 (1990).

60. Sassen, R., Tye, R. S., Chinn, E. W. & Lemoine, R. C. Origin of crude oil in the Wilcox trend of Louisiana and Mississippi: evidence of long-range migration. *Gulf Coast Assoc. Geol. Soc. Trans.* **38**, 27–34 (1988).

61. Seifert, W. K. & Moldowan, J. M. in *Methods in Geochemistry and Geophysics* Vol. 24 (ed. Johns, R. B.) 261–290 (Elsevier, 1986).

62. van Aarssen, B. G. K., Bastow, T. P., Alexander, R. & Kagi, R. I. Distributions of methylated naphthalenes in crude oils: indicators of maturity, biodegradation and mixing. *Org. Geochem.* **30**, 1213–1227 (1999).

63. Wenger, L. M., Davis, C. L. & Isaksen, G. H. Multiple controls on petroleum biodegradation and impact on oil quality. *SPE Reserv. Eval. Eng.* **5**, 375–383 (2002).

64. Kharaka, Y. K. & Hanor, J. S. Deep fluids in the continents: I. Sedimentary basins. *Treatise Geochem.* **5**, 499–540 (2003).

65. Day, J. M. D. et al. The helium flux from the continents and ubiquity of low-³He/³He recycled crust and lithosphere. *Geochim. Cosmochim. Acta* **153**, 116–133 (2015).

66. McCulloh, R. P. Little Creek structure, northwestern La Salle Parish, Louisiana. In *Geological Society of America South Central Section 50th Annual Meeting Paper no. 15-3* (Geological Society of America, 2016); https://doi.org/10.1130/abs/2016SC-271780

67. Ballentine, C. J., Marty, B., Lollar, B. S. & Cassidy, M. Neon isotopes constrain convection and volatile origin in the Earth's mantle. *Nature* **433**, 33–38 (2005).

68. Kennedy, B. M., Hiyagon, H. & Reynolds, J. H. Crustal neon: a striking uniformity. *Earth Planet. Sci. Lett.* **98**, 277–286 (1990).

69. Clarke, W. B., Jenkins, W. J. & Top, Z. Determination of tritium by mass spectrometric measurement of ³He. *Int. J. Appl. Radiat. Isot.* **27**, 515–522 (1976).

70. Lee, J.-Y. et al. A redeterminaton of the isotopic abundances of atmospheric Ar. *Geochim. Cosmochim. Acta* **70**, 4507–4512 (2006).

71. Graham, D. W. Noble gas isotope geochemistry of mid-ocean ridge and ocean island basalt: characterization of mantle source reservoirs. *Rev. Mineral. Geochem.* **47**, 247–317 (2002).

72. Warwick, P. D., Breland, F. C. J., Ratchford, M. E., Ingram, S. L. S. & Pashin, J. C. Coal gas resource potential of Cretaceous and Paleogene coals of the eastern Gulf of Mexico Coastal Plain. In *Southeast Section Meeting* (Geological Society of America, 2004).

73. Prinzhofer, A. A. & Huc, A. Y. Genetic and post-genetic molecular and isotopic fractionations in natural gases. *Chem. Geol.* **126**, 281–290 (1995).

74. Prinzhofer, A. & Pernaton, É. Isotopically light methane in natural gas: bacterial imprint or diffusive fractionation? *Chem. Geol.* **142**, 193–200 (1997).

75. Fuex, A. N. Experimental evidence against an appreciable isotopic fractionation of methane during migration. *Phys. Chem. Earth* **12**, 725–732 (1980).

76. Zhang, T., Wang, X. & Chen, J. Chemical composition of gases as a geochemical tracer of natural gas migration. *Acta Sedimentol. Sinica* **17**, 627–632 (1999).

77. Xia, X. & Tang, Y. Isotope fractionation of methane during natural gas flow with coupled diffusion and adsorption/desorption. *Geochim. Cosmochim. Acta* **77**, 489–503 (2012).

78. Giggenbach, W. F. Carbon-13 exchange between CO₂ and CH₄ under geothermal conditions. *Geochim. Cosmochim. Acta* **46**, 159–165 (1982).

79. 50. Ash, J. L. et al. Exchange catalysis during anaerobic methanotrophy revealed by Δ¹⁵CH₃D and Δ¹⁵CH₃D in methane. *Geochem. Persp. Lett.* **10**, 26–30 (2019).

80. Holler, T. et al. Carbon and sulfur back flux during anaerobic microbial oxidation of methane and coupled sulfate reduction. *Proc. Natl Acad. Sci. USA* **108**, E1484–E1490 (2011).

81. Yoshinaga, M. Y. et al. Carbon isotope equilibration during sulphate-limited anaerobic oxidation of methane. *Nat. Geosci* **7**, 190–194 (2014).

82. Doerfert, S. N., Reichlen, M., Iyer, P., Wang, M. & Ferry, J. G. *Methanolobus zinderi* sp. nov., a methylotrophic methanogen isolated from a deep subsurface coal seam. *Int. J. Syst. Evol. Microbiol.* **59**, 1064–1069 (2009).

83. Larter, S. R. et al. Biodegradation, gas destruction and methane generation in deep subsurface petroleum reservoirs: an overview. *Geol. Soc. Lond. Pet. Geol. Conf. Ser.* **6**, 633–639 (2005).
84. Zinder, S. H. in *Methanogenesis: Ecology, Physiology, Biochemistry & Genetics* (ed Ferry, J. G.) 128–206 (Springer, 1993); https://doi.org/10.1007/978-1-4615-2391-8_4

85. Oren, A. Thermodynamic limits to microbial life at high salt concentrations. *Environ. Microbiol.* **13**, 1908–1923 (2011).

86. Capone, D. G. & Kiene, R. P. Comparison of microbial dynamics in marine and freshwater sediments: contrasts in anaerobic carbon catabolism. *Linnol. Oceanogr.* **33**, 725–749 (1988).

87. Lovley, D. R. & Klug, M. J. Sulfate reducers can outcompete methanogens at freshwater sulfate concentrations. *Appl. Environ. Microbiol.* **45**, 187–192 (1983).

88. Rice, D. D. & Claypool, G. E. Generation, accumulation, and resource potential of biogenic gas. *AANG Bull.* **65**, 5–25 (1981).

89. van Bodegom, P. M., Scholten, J. C. M. & Stams, A. J. M. Direct inhibition of methanogenesis by ferric iron. *FEMS Microbiol. Ecol.* **49**, 261–268 (2004).

90. McIntosh, J. C. & Martini, A. M. Hydrogeochemical indicators for microbial methane in fractured organic-rich shales: case studies of the Antrim, New Albany, and Ohio shales. **13**, 162–174 (2008).

91. Etiope, G. & Lollar, B. S. Abiotic methane on Earth. *Rev. Geophys.* **51**, 276–299 (2013).

92. Chung, H. M., Gormly, J. R. & Squires, R. M. Origin of gaseous hydrocarbons in subsurface environments: theoretical considerations of carbon isotope distribution. *Chem. Geol.* **71**, 97–104 (1988).

93. Whiticar, M. J., Faber, E. & Schoell, M. Biogenic methane formation in marine and freshwater environments: CO₂ reduction vs. acetate fermentation—isotope evidence. *Geochim. Cosmochim. Acta* **50**, 693–709 (1986).

94. Conrad, R. Quantification of methanogenic pathways using stable carbon isotopic signatures: a review and a proposal. *Org. Geochem.* **36**, 739–752 (2005).

95. Ballentine, C. J., Schoell, M., Coleman, D. & Cain, B. A. 300-Myr-old magmatic CO₂ in natural gas reservoirs of the west Texas Permian Basin. *Nature* **409**, 327–331 (2001).

96. Holland, G. & Gilfillan, S. in *The Noble Gases as Geochemical Tracers* (ed. Burnard, P.) 177–223 (Springer, 2013); https://doi.org/10.1007/978-3-642-28836-4_8

97. Kämpf, H., Geissler, W. H. & Bäuer, K. in *Mantle Plumes: A Multidisciplinary Approach* (eds Ritter J.R.R. & Christensen U.R.) 127–158 (Springer, 2007).

98. Royden, L., Horváth, F., Nagymarosy, A. & Stegena, L. Evolution of the Pannonian Basin system: 2. Subsidence and thermal history. *Tectonics* **2**, 91–137 (1983).

99. Sherwood Lollar, B., O’Nions, R. K. & Ballentine, C. J. Helium and neon isotope systematics in carbon dioxide-rich and hydrocarbon-rich gas reservoirs. *Geochim. Cosmochim. Acta* **58**, 5279–5290 (1994).

100. Palcsu, L. et al. In-reservoir mixing of mantle-derived CO₂ and metasedimentary CH₄–N₂ fluids—noble gas and stable isotope study of two multistacked fields (Pannonian Basin System, W-Hungary). *Mar. Pet. Geol.* **54**, 216–227 (2014).

101. Sheng, X., Nakai, S., Wakita, H., Yongchang, X. & Wang, X. Carbon isotopes of hydrocarbons and carbon dioxide in natural gases in China. *J. Asian Earth Sci.* **15**, 89–101 (1997).

**Acknowledgements** R.L.T. was supported by a Natural Environment Research Council studentship (grant reference NE/LO02612/1). C.J.B. and P.H.B. acknowledge A. Regberg and B. Meurer for their support of the project and help with sample collection. C.J.B. was part supported by an Earth4D CIFAR fellowship. P.H.B. was supported by NSF awards 1923915 and 2015789. O.W. was supported by Natural Sciences and Engineering Research Council of Canada Discovery and Accelerator grants awarded to the Sherwood Lollar research group and acknowledges B. Sherwood Lollar’s support for the project. Z.M.S. acknowledges J. Biddle and G. Christman for their help in generating the microbial data.

**Author contributions** The project was conceived by R.L.T., P.H.B., M.L. and C.J.B. R.L.T. performed the noble gas isotopic analysis, managed the project and prepared the first draft of the manuscript. P.H.B. collected the samples. Z.M.S., H.X. and B.S. conducted the microbial, clumped isotope and molecular geochemical analysis, respectively. Clumped isotope modelling was developed by H.X. and J.M.E. Integrated noble gas and stable isotope modelling was developed by R.L.T., P.H.B., M.L., C.J.B., O.W. and D.J.B. All authors contributed to the final manuscript.

**Competing interests** M.F. and B.S. are employed by/an employee of ExxonMobil Upstream Integrated Solutions Company. Z.M.S. is employed by/an employee of ExxonMobil Research and Engineering Company. When the manuscript was first submitted, M.L. was employed by/ an employee of ExxonMobil Upstream Integrated Solutions Company. M.L. is now employed by/an employee of Aker BP. The views expressed are those of the author(s) and not necessarily those of ExxonMobil Upstream Integrated Solutions Company, ExxonMobil Research and Engineering Company, or Aker BP.

**Additional information**

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41586-021-04153-3.

**Correspondence and requests for materials** should be addressed to R. L. Tyne, P. H. Barry, M. Lawson or C. J. Ballentine.

**Peer review information** *Nature* thanks Werner Aeschbach, Jenna Shelton and Marc Strous for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permissions information** is available at http://www.nature.com/reprints.
<figure><img src="image_7.png" /></figure>

**Extended Data Fig. 1** | Relationship between noble gas isotopic ratios in the Olla (blue circles) and Nebo-Hemphill (orange triangles) oil fields. Both fields show a resolvable mantle component that is significantly greater in the Olla Oil Field. One standard deviation error bars are within symbol size. a, He isotopic ratios versus ⁴He/²⁰Ne ratio. Air-corrected He isotopes are reported relative to air ($R_A = 1.38 \times 10^{-6}$ (ref. ⁶⁹)). All samples are in excess of typical radiogenic production (0.02$R_A$), requiring a mantle contribution. High ⁴He/²⁰Ne

indicates negligible air contamination in the samples. b, Three neon isotope plot. The crustal production line⁶⁸ reflects the observed trend in typical crustal fluids. The mantle-air line represents mixing between air and a MORB-like endmember⁶⁷. Data are consistent with mixing between air and a mantle-rich component. In both He and Ne isotopes, the mantle signal is stronger in the Olla Oil Field.
<figure><img src="image_8.png" /></figure>

**Extended Data Fig. 2** | CO₂/³He variation vs. a, b, ⁴He (a) and ²⁰Ne (b)

abundances in the Olla (blue circles) and Nebo-Hemphill (orange triangles) oil fields. 1σ level of uncertainty is shown on the measured samples. A negative correlation between ⁴He, which accumulates in formation waters as a function

<figure><img src="image_9.png" /></figure>

of time and ²⁰Ne (which is groundwater derived) and CO₂/³He, has been interpreted in other CO₂-rich natural gas fields⁴ to support the importance of the formation water in controlling CO₂ removal.
<figure><img src="image_10.png" /><figcaption>**Extended Data Fig. 3** | Hydrocarbon isotopic composition in the Olla (blue circles) and Nebo-Hemphill (orange triangles) oil fields.</figcaption></figure>

a, Comparison of δ¹³C-CH₄ and δD-CH₄, with typical CH₄ provenances after Whiticar⁹³. One standard deviation error bars are within symbol size. In isolation, the Olla samples (blue circles, this study; unfilled circles, from ref. ¹⁰) are consistent with a mainly thermogenic origin, whereas Nebo-Hemphill (orange triangles, this study; unfilled triangles, from ref. ¹⁰) appears to be a

microbial-thermogenic mix. Clumped isotopologues nevertheless provide clear evidence that the Olla CH₄ is microbial in origin (Fig. 3). b, Plot of δ¹³C of Cₙ versus 1/Cₙ variation in natural gases⁹². One standard deviation error bars are within symbol size. Both the Olla (blue circles) and Nebo-Hemphill (orange triangles) fields show an enrichment in their isotopic signature especially in propane; however, this is most pronounced within the Olla Oil Field.
<figure><img src="image_11.png" /></figure>

**Extended Data Fig. 4** | Schematic of the processes occurring within the Olla Oil Field and resulting changes in the isotopic composition within the field. The first panel shows CO₂ injection into the 2,800-ft sandstone (sst) and subsequent lateral and vertical migration (yellow arrows) to all the producing formations within the field. The middle panel demonstrate the different processes (microbial methangonesis, AOM and dissolution), which are

occurring in the field and their effects on the isotopic composition.

Yellow circles represent CO₂ dissolution and green circles are the net CH₄ production from microbial activity. Hydrogen for microbial methanogenesis is sourced from the hydrocarbons and water. The final panel shows the current state of the reservoir.
<figure><img src="image_12.png" /><figcaption>**Extended Data Fig. 5** | Microbial community analysis of most abundant taxa from the SSU rRNA gene-sequence region. Bar charts represent the most resolved taxa from amplicon sequence variants generated in QIIME</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>2⁴⁹,⁵⁰,⁵¹ for archaeal primers (a) and bacterial primers (b). Hydrogenotrophic methanogenesis-capable taxa are represented with bold red font.</figcaption></figure>
<figure><img src="image_14.png" /></figure>

**Extended Data Fig. 6** | Carbon isotopic composition of the Olla samples as a function of CO₂ concentration for data from this study (dark blue circles) and those collected without ³He analysis¹⁰ (unfilled circles). One standard deviation error bars are within symbol size. Dashed lines show endmember methanogenesis and dissolution fractionation (at pH 7) trajectories. Tick marks represent the total amount of CO₂ trapping within the system, relative to sample O5. The shaded region represents trapping by the combination of both microbial methanogenesis and dissolution. The same combination of

dissolution and microbial methanogenesis has been modelled (grey region, cf. Fig. 2). Lines labelled ‘consumption’ at 5%, 13% and 19% show the proportion of original CO₂ that has been removed by net microbial methanogenesis. One sample from ref. ¹⁰ has a CO₂ concentration greater than our ‘most pristine’ CO₂/³He sample (O5). This illustrates that by using sample O5 as our least altered composition we calculate a conservative (minimum) amount of CO₂ consumption, and that initial CO₂ concentrations were probably greater.
<figure><img src="image_15.png" /></figure>

**Extended Data Fig. 7** | Carbon composition of the Pannonian Basin samples. a, Szegholm South<sup>20,99</sup> (red circles), Szegholm North<sup>20,99</sup> (purple triangles), Ebes<sup>20,99</sup> (blue diamonds) and Hajduszoboszlo<sup>20,99</sup> (grey squares) gas field samples as a function of their CO<sub>2</sub> concentrations and CO<sub>2</sub>/<sup>3</sup>He. Dashed lines represent the projection of CO<sub>2</sub> trapping in the fields. b, The Szegholm South field had multiple CO<sub>2</sub> isotope data and this field has been used to investigate the processes controlling this pattern using the carbon isotopic composition as a function of CO<sub>2</sub>/<sup>3</sup>He. Dashed lines show endmember

<figure><img src="image_16.png" /></figure>

methanogenesis and dissolution fractionation (at pH 7) trajectories. Tick marks represent the total amount of CO₂ trapping within the system. The shaded region represents trapping by the combination of both microbial methanogenesis and dissolution. The upper and lower methanogenesis: dissolution ratios (M:D) are 0.050 and 0.060, respectively, showing dissolution accounts for approximately 16 times (M:D = 0.055) more CO₂ removal than microbial methanogenesis. Lines labelled ‘consumption’ show the proportion of original CO₂ that has been removed by net microbial methanogenesis.
<table><caption>Extended Data Table 1 | Geochemistry and isotopic composition of the formation waters</caption><thead><tr><th>Sample</th><th>Latitude, Longitude</th><th>Producing zone</th><th>Depth m</th><th>Cl mmol/l</th><th>Na mmol/l</th><th>Mg ppm</th><th>Ca ppm</th><th>K ppm</th><th>Ba ppm</th><th>Si ppm</th><th>Sr ppm</th><th>Fe ppm</th><th>SO<sub>4</sub><sup>2-</sup> ppm</th><th>Bicarbonate ppm</th></tr></thead><tbody><tr><td>O1</td><td>31.83128, -92.2095</td><td>L Cent</td><td>677</td><td>1435</td><td>1409</td><td>222</td><td>181</td><td>223</td><td>26.7</td><td>11.3</td><td>83</td><td>N.D.</td><td>62</td><td>3655</td></tr><tr><td>O2</td><td>31.83015, -92.2093</td><td>L Cent</td><td>656</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>O3</td><td>31.82453, -92.2137</td><td>2800 ft</td><td>835</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>O4</td><td>31.80899, -92.2297</td><td>U Cent</td><td>663</td><td>1122</td><td>1108</td><td>329</td><td>80.2</td><td>228</td><td>95</td><td>10.1</td><td>77</td><td>N.D.</td><td>&lt;1</td><td>3912</td></tr><tr><td>O5</td><td>31.80119, -92.2224</td><td>2800 ft</td><td>842</td><td>1357</td><td>1340</td><td>305</td><td>500</td><td>420</td><td>165</td><td>10.6</td><td>215</td><td>N.D.</td><td>&lt;1</td><td>3164</td></tr><tr><td>O6</td><td>31.76926, -92.2133</td><td>H2-SD</td><td>1051</td><td>1621</td><td>1630</td><td>119</td><td>343</td><td>524</td><td>158</td><td>12.5</td><td>207</td><td>N.D.</td><td>2</td><td>3486</td></tr><tr><td>NHH1</td><td>31.60037, -92.1668</td><td>Wilcox B2</td><td>1211</td><td>1538</td><td>1502</td><td>581</td><td>761</td><td>212</td><td>198</td><td>12.0</td><td>265</td><td>N.D.</td><td>&lt;1</td><td>812</td></tr><tr><td>NHH2</td><td>31.59793, -92.1468</td><td>Wilcox A2</td><td>1027</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NHH3</td><td>31.58830, -92.1752</td><td>T Wilcox</td><td>753</td><td>1497</td><td>1452</td><td>495</td><td>924</td><td>157</td><td>188</td><td>7.2</td><td>267</td><td>2.1</td><td>&lt;1</td><td>108</td></tr><tr><td>NHH4</td><td>31.58895, -92.1775</td><td>Wilcox C3</td><td>1099</td><td>1579</td><td>1507</td><td>548</td><td>862</td><td>188</td><td>210</td><td>10.4</td><td>266</td><td>2.3</td><td>&lt;1</td><td>308</td></tr><tr><td>NHH5</td><td>31.58907, -92.1666</td><td>Wilcox F3a &amp; F3b</td><td>1244</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NHH6</td><td>31.58238, -92.1487</td><td>Wilcox F2</td><td>1202</td><td>1633</td><td>1562</td><td>474</td><td>1122</td><td>195</td><td>216</td><td>10.4</td><td>270</td><td>7.5</td><td>&lt;1</td><td>38</td></tr><tr><td>NHH7</td><td>31.57380, -92.1512</td><td>Wilcox H2</td><td>1273</td><td>1786</td><td>1662</td><td>481</td><td>1618</td><td>242</td><td>203</td><td>10.4</td><td>324</td><td>20.2</td><td>&lt;1</td><td>162</td></tr><tr><td>NHH8</td><td>31.56586, -92.1557</td><td>Wilcox A1</td><td>1023</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
Extended Data Table 2 | Sample location, reservoir temperature, major gas species, carbon isotopes and clumped isotope data

<table><thead><tr><th rowspan="2">Sample</th><th>Reservoir temp*</th><th>CO<sub>2</sub><sup>†</sup></th><th>CH<sub>4</sub><sup>†</sup></th><th>δ<sup>13</sup>C-<br>CO<sub>2</sub></th><th>δ<sup>13</sup>C-<br>CH<sub>4</sub></th><th>δ<sup>13</sup>C-CO<sub>2</sub><br>-CH<sub>4</sub></th><th>CO<sub>2</sub> <sup>2</sup>He</th><th><sup>3</sup>He/<sup>4</sup>He</th><th><sup>4</sup>He</th><th>Δ<sup>13</sup>CH<sub>3</sub>D</th><th>Apparent<br>Δ<sup>13</sup>CH<sub>3</sub>D<br>Temperature</th><th>Δ<sup>12</sup>CH<sub>2</sub>D<sub>2</sub></th><th>Apparent<br>Δ<sup>13</sup>CH<sub>2</sub>D<sub>2</sub><br>Temperature</th></tr><tr><th>°C</th><th>mol<br>%(STP)</th><th>mol<br>%(STP)</th><th>%<sub>VDPB</sub></th><th>%<sub>VDPB</sub></th><th>%<sub>VDPB</sub></th><th>x10<sup>2</sup></th><th></th><th></th><th>x10<sup>-4</sup><br>cm<sup>3</sup>/cm<sup>3</sup>(STP)</th><th>°C</th><th></th><th>°C</th></tr></thead><tbody><tr><td>O1</td><td>39.1</td><td>11</td><td>86</td><td>16.4</td><td>-47.6</td><td>64</td><td>2.77 ± 0.06</td><td>1.78 ± 0.04</td><td>23.7 ± 0.2</td><td>4.89 ± 0.36</td><td>56<sup>+17</sup><sub>-19</sub></td><td>12.41 ± 1.29</td><td>88<sup>+15</sup><sub>-14</sub></td></tr><tr><td>O2</td><td>39.7</td><td>6.7</td><td>89</td><td>13.7</td><td>-49.3</td><td>63</td><td>1.92 ± 0.04</td><td>1.79 ± 0.04</td><td>20.7 ± 0.2</td><td>4.47 ± 0.34<sup>‡</sup></td><td>74<sup>+16</sup><sub>-14</sub></td><td></td><td></td></tr><tr><td>O3</td><td>44.4</td><td>9.4</td><td>88</td><td>15.1</td><td>-49.3</td><td>64.4</td><td>2.47 ± 0.06</td><td>1.93 ± 0.04</td><td>20.9 ± 0.2</td><td>5.62 ± 0.31<sup>‡</sup></td><td>30<sup>+10</sup><sub>-10</sub></td><td>9.13 ± 1.60</td><td>132<sup>+30</sup><sub>-33</sub></td></tr><tr><td>O4</td><td>39.4</td><td>18</td><td>73</td><td>15.8</td><td>-48.2</td><td>64</td><td>6.60 ± 0.15</td><td>1.72 ± 0.04</td><td>16.0 ± 0.2</td><td>3.93 ± 0.52</td><td>100<sup>+29</sup><sub>-28</sub></td><td>12.19 ± 1.39</td><td>90<sup>+17</sup><sub>-15</sub></td></tr><tr><td>O5</td><td>44.5</td><td>38</td><td>57</td><td>7.05</td><td>-46.7</td><td>53.75</td><td>17.3 ± 0.4</td><td>1.98 ± 0.04</td><td>11.8 ± 0.1</td><td>5.66 ± 0.39</td><td>29<sup>+14</sup><sub>-12</sub></td><td>12.46 ± 1.07</td><td>87<sup>+13</sup><sub>-11</sub></td></tr><tr><td>O6</td><td>50.7</td><td>27</td><td>67</td><td>13.3</td><td>-44.8</td><td>58.1</td><td>10.6 ± 0.2</td><td>1.81 ± 0.04</td><td>15.1 ± 0.2</td><td>3.45 ± 0.26<sup>‡</sup></td><td>128<sup>+17</sup><sub>-16</sub></td><td>9.82 ± 1.53</td><td>121<sup>+25</sup><sub>-21</sub></td></tr><tr><td>NHH1</td><td>51</td><td>1.9</td><td>97</td><td>13.9</td><td>-46</td><td>59.9</td><td></td><td></td><td></td><td>3.77 ± 0.37</td><td>109<sup>+24</sup><sub>-21</sub></td><td>4.27 ± 1.29</td><td>252<sup>+64</sup><sub>-44</sub></td></tr><tr><td>NHH2</td><td>50</td><td>0.28</td><td>96</td><td>0.33</td><td>-49.8</td><td>49.5</td><td>1.44 ± 0.06</td><td>0.38 ± 0.04</td><td>5.41 ± 0.04</td><td>4.38 ± 0.28<sup>‡</sup></td><td>78<sup>+13</sup><sub>-12</sub></td><td>3.85 ± 1.81</td><td>270<sup>+118</sup><sub>-65</sub></td></tr><tr><td>NHH3</td><td>41.9</td><td>0.77</td><td>98</td><td>7.1</td><td>-61</td><td>68.1</td><td>1.79 ± 0.05</td><td>0.39 ± 0.01</td><td>11.8 ± 0.1</td><td>3.17 ± 0.35</td><td>146<sup>+22</sup><sub>-20</sub></td><td>5.03 ± 1.24</td><td>224<sup>+48</sup><sub>-36</sub></td></tr><tr><td>NHH4</td><td>52.1</td><td>1.7</td><td>94</td><td>9.4</td><td>-54.5</td><td>63.9</td><td>9.97 ± 0.12</td><td>0.52 ± 0.01</td><td>2.92 ± 0.03</td><td>3.16 ± 0.54</td><td>147<sup>+32</sup><sub>-39</sub></td><td>1.56 ± 1.67</td><td>443<sup>+118</sup><sub>-142</sub></td></tr><tr><td>NHH5</td><td>56</td><td>1.8</td><td>76</td><td>3.08</td><td>-57</td><td>60.08</td><td>12.2 ± 0.3</td><td>0.58 ± 0.01</td><td>2.74 ± 0.03</td><td>3.28 ± 0.40</td><td>139<sup>+34</sup><sub>-29</sub></td><td>1.80 ± 1.13</td><td>413<sup>+222</sup><sub>-24</sub></td></tr><tr><td>NHH6</td><td>55</td><td>0.99</td><td>84</td><td>2.52</td><td>-57.5</td><td>60.02</td><td>2.57 ± 0.6</td><td>0.46 ± 0.01</td><td>8.93 ± 0.07</td><td>2.90 ± 0.42</td><td>166<sup>+34</sup><sub>-29</sub></td><td>0.01 ± 1.63</td><td>12700<sup>+n/a</sup><sub>-12200</sub></td></tr><tr><td>NHH7</td><td>57.1</td><td>0.39</td><td>99</td><td>-0.93</td><td>-59.1</td><td>58.17</td><td>0.53 ± 0.02</td><td>0.50 ± 0.01</td><td>15.6 ± 0.2</td><td>3.33 ± 0.57</td><td>135<sup>+44</sup><sub>-35</sub></td><td>2.02 ± 1.38</td><td>390<sup>+257</sup><sub>-98</sub></td></tr><tr><td>NHH8</td><td>49.8</td><td>0.35</td><td>97</td><td>3.31</td><td>-56.2</td><td>59.51</td><td>0.84 ± 0.03</td><td>0.39 ± 0.01</td><td>9.64 ± 0.10</td><td>3.64 ± 0.35</td><td>116<sup>+25</sup><sub>-20</sub></td><td>5.14 ± 1.10</td><td>211<sup>+41</sup><sub>-31</sub></td></tr></tbody></table>

<footnote n="10">10 errors are given. Location is given as latitude and longitude in decimal degrees. The <sup>3</sup>He/<sup>4</sup>He ratio, R, is shown relative to the <sup>3</sup>He/<sup>4</sup>He ratio in air, R<sub>A</sub>. *Adiabatic corrected. ′Air corrected. ′Corrected Δ18 from prototype Ultra).</footnote>
<table><caption>Extended Data Table 3 | Noble gas concentrations and ratios for the sampled wells (with 1σ errors), as well as measured carbon isotopic compositions of alkanes</caption><thead><tr><th rowspan="2">Sample</th><th><sup>26</sup>Ne</th><th><sup>26</sup>Ar</th><th><sup>54</sup>Kr</th><th><sup>133</sup>Xe</th><th><sup>20</sup>Ne/<sup>22</sup>Ne</th><th><sup>21</sup>Ne/<sup>22</sup>Ne</th><th><sup>40</sup>Ar/<sup>26</sup>Ar</th><th><sup>85</sup>Kr/<sup>84</sup>Kr</th><th><sup>123</sup>Xe/<sup>130</sup>Xe</th><th>δ<sup>13</sup>C-C<sub>2</sub></th><th>δ<sup>13</sup>C-C<sub>3</sub></th><th>δ<sup>13</sup>C-nC<sub>4</sub></th></tr><tr><th>x10<sup>-3</sup> cm<sup>3</sup>/cm<sup>3</sup><br>(STP)</th><th>x10<sup>-3</sup> cm<sup>3</sup>/cm<sup>3</sup><br>(STP)</th><th>x10<sup>-3</sup> cm<sup>3</sup>/cm<sup>3</sup><br>(STP)</th><th>x10<sup>-10</sup> cm<sup>3</sup>/cm<sup>3</sup><br>(STP)</th><th></th><th></th><th></th><th></th><th></th><th>%<sub>VPDB</sub></th><th>%<sub>VPDB</sub></th><th>%<sub>VPDB</sub></th></tr></thead><tbody><tr><td>O1</td><td>25.8 ± 0.5</td><td>84.9 ± 0.8</td><td>3.21 ± 0.09</td><td>3.46 ± 0.04</td><td>10.25 ± 0.03</td><td>0.0384 ± 0.0004</td><td>1910 ± 17</td><td>0.302 ± 0.002</td><td>6.59 ± 0.04</td><td></td><td></td><td></td></tr><tr><td>O2</td><td>26.9 ± 0.5</td><td>80.8 ± 0.8</td><td>3.66 ± 0.10</td><td>3.78 ± 0.05</td><td>10.18 ± 0.03</td><td>0.0371 ± 0.0003</td><td>1797 ± 16</td><td>0.302 ± 0.002</td><td>6.61 ± 0.04</td><td></td><td></td><td></td></tr><tr><td>O3</td><td>20.6 ± 0.4</td><td>83.3 ± 0.8</td><td>1.81 ± 0.05</td><td>2.04 ± 0.03</td><td>10.42 ± 0.03</td><td>0.0405 ± 0.0004</td><td>2344 ± 20</td><td>0.302 ± 0.002</td><td>6.60 ± 0.04</td><td>-27.3</td><td>-17.74</td><td>-22.54</td></tr><tr><td>O4</td><td>10.5 ± 0.2</td><td>66.0 ± 0.6</td><td>3.20 ± 0.09</td><td>9.04 ± 0.10</td><td>10.19 ± 0.03</td><td>0.0385 ± 0.0003</td><td>1475 ± 13</td><td>0.302 ± 0.002</td><td>6.57 ± 0.04</td><td></td><td></td><td></td></tr><tr><td>O5</td><td>10.1 ± 0.2</td><td>45.6 ± 0.4</td><td>2.07 ± 0.06</td><td>2.77 ± 0.03</td><td>10.42 ± 0.03</td><td>0.0409 ± 0.0004</td><td>2484 ± 22</td><td>0.302 ± 0.002</td><td>6.61 ± 0.04</td><td></td><td></td><td></td></tr><tr><td>O6</td><td>8.83 ± 0.17</td><td>46.0 ± 0.4</td><td>2.17 ± 0.05</td><td>2.83 ± 0.03</td><td>10.64 ± 0.03</td><td>0.0459 ± 0.0004</td><td>2766 ± 24</td><td>0.302 ± 0.002</td><td>6.59 ± 0.04</td><td>-28.14</td><td>-20.22</td><td>-21.58</td></tr><tr><td>NHH1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NHH2</td><td>10.3 ± 0.2</td><td>44.5 ± 0.4</td><td>2.42 ± 0.06</td><td>2.12 ± 0.03</td><td>9.93 ± 0.03</td><td>0.0330 ± 0.0003</td><td>435.7 ± 3.8</td><td>0.303 ± 0.002</td><td>6.51 ± 0.04</td><td>-24.19</td><td>-25.87</td><td>-26.00</td></tr><tr><td>NHH3</td><td>36.7 ± 0.7</td><td>176 ± 2</td><td>5.95 ± 0.16</td><td>7.39 ± 0.09</td><td>9.83 ± 0.03</td><td>0.0310 ± 0.0003</td><td>369.7 ± 3.2</td><td>0.302 ± 0.002</td><td>6.52 ± 0.04</td><td></td><td></td><td></td></tr><tr><td>NHH4</td><td>3.90 ± 0.07</td><td>66.0 ± 0.7</td><td>4.09 ± 0.11</td><td>6.04 ± 0.07</td><td>9.76 ± 0.03</td><td>0.0342 ± 0.0003</td><td>398.6 ± 3.5</td><td>0.302 ± 0.002</td><td>6.50 ± 0.04</td><td>-35.37</td><td>-26.52</td><td>-27.58</td></tr><tr><td>NHH5</td><td>3.49 ± 0.07</td><td>47.2 ± 0.4</td><td>5.23 ± 0.13</td><td>10.53 ± 0.11</td><td>9.87 ± 0.03</td><td>0.0353 ± 0.0003</td><td>573.1 ± 5.0</td><td>0.304 ± 0.002</td><td>6.47 ± 0.04</td><td>-36.96</td><td>-28.08</td><td>-26.62</td></tr><tr><td>NHH6</td><td>14.9 ± 0.3</td><td>79.1 ± 0.7</td><td>8.14 ± 0.22</td><td>6.15 ± 0.08</td><td>9.84 ± 0.03</td><td>0.0335 ± 0.0003</td><td>494.6 ± 4.3</td><td>0.304 ± 0.002</td><td>6.50 ± 0.04</td><td></td><td></td><td></td></tr><tr><td>NHH7</td><td>21.6 ± 0.4</td><td>112 ± 1</td><td>4.70 ± 0.12</td><td>3.99 ± 0.05</td><td>10.02 ± 0.03</td><td>0.0333 ± 0.0003</td><td>603.8 ± 5.3</td><td>0.304 ± 0.002</td><td>6.47 ± 0.04</td><td></td><td></td><td></td></tr><tr><td>NHH8</td><td>7.08 ± 0.13</td><td>78.6 ± 0.8</td><td>2.32 ± 0.06</td><td>4.81 ± 0.06</td><td>9.77 ± 0.03</td><td>0.0308 ± 0.0003</td><td>379.4 ± 3.3</td><td>0.303 ± 0.002</td><td>6.51 ± 0.04</td><td></td><td></td><td></td></tr><tr><td>Air</td><td></td><td></td><td></td><td></td><td>9.81</td><td>0.0290</td><td>298.56</td><td>0.303</td><td>6.49</td><td></td><td></td><td></td></tr></tbody></table>