Eutrophication and Deoxygenation Drive High Methane Emissions from a Brackish Coastal System.

Coastal environments are a major source of marine methane in the atmosphere. Eutrophication and deoxygenation have the potential to amplify the coastal methane emissions. Here, we investigate methane dynamics in the eutrophic Stockholm Archipelago. We cover a range of sites with contrasting water column redox conditions and rates of organic matter degradation, with the latter reflected by the depth of the sulfate-methane transition zone (SMTZ) in the sediment. We find the highest benthic release of methane (2.2-8.6 mmol m-2 d-1) at sites where the SMTZ is located close to the sediment-water interface (2-10 cm). A large proportion of methane is removed in the water column via aerobic or anaerobic microbial pathways. At many locations, water column methane is highly depleted in 13C, pointing toward substantial bubble dissolution. Calculated and measured rates of methane release to the atmosphere range from 0.03 to 0.4 mmol m-2 d-1 and from 0.1 to 1.7 mmol m-2 d-1, respectively, with the highest fluxes at locations with a shallow SMTZ and anoxic and sulfidic bottom waters. Taken together, our results show that sites suffering most from both eutrophication and deoxygenation are hotspots of coastal marine methane emissions.

Anaerobic methanotrophy is stimulated by graphene oxide in a brackish urban canal sediment.

Anthropogenic activities are influencing aquatic environments through increased chemical pollution and thus are greatly affecting the biogeochemical cycling of elements. This has increased greenhouse gas emissions, particularly methane, from lakes, wetlands, and canals. Most of the methane produced in anoxic sediments is converted into carbon dioxide by methanotrophs before it reaches the atmosphere. Anaerobic oxidation of methane requires an electron acceptor such as sulphate, nitrate, or metal oxides. Here, we explore the anaerobic methanotrophy in sediments of three urban canals in Amsterdam, covering a gradient from freshwater to brackish conditions. Biogeochemical analysis showed the presence of a shallow sulphate-methane transition zone in sediments of the most brackish canal, suggesting that sulphate could be a relevant electron acceptor for anaerobic methanotrophy in this setting. However, sediment incubations amended with sulphate or iron oxides (ferrihydrite) did not lead to detectable rates of methanotrophy. Despite the presence of known nitrate-dependent anaerobic methanotrophs (Methanoperedenaceae), no nitrate-driven methanotrophy was observed in any of the investigated sediments either. Interestingly, graphene oxide stimulated anaerobic methanotrophy in incubations of brackish canal sediment, possibly catalysed by anaerobic methanotrophs of the ANME-2a/b clade. We propose that natural organic matter serving as electron acceptor drives anaerobic methanotrophy in brackish sediments.

Photocatalytic chlorine atom production on mineral dust-sea spray aerosols over the North Atlantic.

Active chlorine in the atmosphere is poorly constrained and so is its role in the oxidation of the potent greenhouse gas methane, causing uncertainty in global methane budgets. We propose a photocatalytic mechanism for chlorine atom production that occurs when Sahara dust mixes with sea spray aerosol. The mechanism is validated by implementation in a global atmospheric model and thereby explaining the episodic, seasonal, and location-dependent 13C depletion in CO in air samples from Barbados [J.E. Mak, G. Kra, T. Sandomenico, P. Bergamaschi, J. Geophys. Res. Atmos. 108 (2003)], which remained unexplained for decades. The production of Cl can also explain the anomaly in the CO:ethane ratio found at Cape Verde [K. A. Read et al., J. Geophys. Res. Atmos. 114 (2009)], in addition to explaining the observation of elevated HOCl [M. J. Lawler et al., Atmos. Chem. Phys. 11, 7617-7628 (2011)]. Our model finds that 3.8 Tg(Cl) y-1 is produced over the North Atlantic, making it the dominant source of chlorine in the region; globally, chlorine production increases by 41%. The shift in the methane sink budget due to the increased role of Cl means that isotope-constrained top-down models fail to allocate 12 Tg y-1 (2% of total methane emissions) to 13C-depleted biological sources such as agriculture and wetlands. Since 2014, an increase in North African dust emissions has increased the 13C isotope of atmospheric CH4, thereby partially masking a much greater decline in this isotope, which has implications for the interpretation of the drivers behind the recent increase of methane in the atmosphere.

Global environmental implications of atmospheric methane removal through chlorine-mediated chemistry-climate interactions.

Atmospheric methane is both a potent greenhouse gas and photochemically active, with approximately equal anthropogenic and natural sources. The addition of chlorine to the atmosphere has been proposed to mitigate global warming through methane reduction by increasing its chemical loss. However, the potential environmental impacts of such climate mitigation remain unexplored. Here, sensitivity studies are conducted to evaluate the possible effects of increasing reactive chlorine emissions on the methane budget, atmospheric composition and radiative forcing. Because of non-linear chemistry, in order to achieve a reduction in methane burden (instead of an increase), the chlorine atom burden needs to be a minimum of three times the estimated present-day burden. If the methane removal target is set to 20%, 45%, or 70% less global methane by 2050 compared to the levels in the Representative Concentration Pathway 8.5 scenario (RCP8.5), our modeling results suggest that additional chlorine fluxes of 630, 1250, and 1880 Tg Cl/year, respectively, are needed. The results show that increasing chlorine emissions also induces significant changes in other important climate forcers. Remarkably, the tropospheric ozone decrease is large enough that the magnitude of radiative forcing decrease is similar to that of methane. Adding 630, 1250, and 1880 Tg Cl/year to the RCP8.5 scenario, chosen to have the most consistent current-day trends of methane, will decrease the surface temperature by 0.2, 0.4, and 0.6 °C by 2050, respectively. The quantity and method in which the chlorine is added, its interactions with climate pathways, and the potential environmental impacts on air quality and ocean acidity, must be carefully considered before any action is taken.

Exploring the potential of Δ17O in CO2 for determining mesophyll conductance.

Mesophyll conductance to CO2 from the intercellular air space to the CO2-H2O exchange site has been estimated using δ18O measurements (gm18). However, the gm18 estimates are affected by the uncertainties in the δ18O of leaf water where the CO2-H2O exchange takes place and the degree of equilibration between CO2 and H2O. We show that measurements of Δ17O (i.e.Δ17O=δ17O-0.528×δ18O) can provide independent constraints on gm (gmΔ17) and that these gm estimates are less affected by fractionation processes during gas exchange. The gm calculations are applied to combined measurements of δ18O and Δ17O, and gas exchange in two C3 species, sunflower (Helianthus annuus L. cv. 'sunny') and ivy (Hedera hibernica L.), and the C4 species maize (Zea mays). The gm18 and gmΔ17 estimates agree within the combined errors (P-value, 0.876). Both approaches are associated with large errors when the isotopic composition in the intercellular air space becomes close to the CO2-H2O exchange site. Although variations in Δ17O are low, it can be measured with much higher precision compared with δ18O. Measuring gmΔ17 has a few advantages compared with gm18: (i) it is less sensitive to uncertainty in the isotopic composition of leaf water at the isotope exchange site and (ii) the relative change in the gm due to an assumed error in the equilibration fraction θeq is lower for gmΔ17 compared with gm18. Thus, using Δ17O can complement and improve the gm estimates in settings where the δ18O of leaf water varies strongly, affecting the δ18O (CO2) difference between the intercellular air space and the CO2-H2O exchange site.

Advancing understanding of land-atmosphere interactions by breaking discipline and scale barriers.

Vegetation and atmosphere processes are coupled through a myriad of interactions linking plant transpiration, carbon dioxide assimilation, turbulent transport of moisture, heat and atmospheric constituents, aerosol formation, moist convection, and precipitation. Advances in our understanding are hampered by discipline barriers and challenges in understanding the role of small spatiotemporal scales. In this perspective, we propose to study the atmosphere-ecosystem interaction as a continuum by integrating leaf to regional scales (multiscale) and integrating biochemical and physical processes (multiprocesses). The challenges ahead are (1) How do clouds and canopies affect the transferring and in-canopy penetration of radiation, thereby impacting photosynthesis and biogenic chemical transformations? (2) How is the radiative energy spatially distributed and converted into turbulent fluxes of heat, moisture, carbon, and reactive compounds? (3) How do local (leaf-canopy-clouds, 1 m to kilometers) biochemical and physical processes interact with regional meteorology and atmospheric composition (kilometers to 100 km)? (4) How can we integrate the feedbacks between cloud radiative effects and plant physiology to reduce uncertainties in our climate projections driven by regional warming and enhanced carbon dioxide levels? Our methodology integrates fine-scale explicit simulations with new observational techniques to determine the role of unresolved small-scale spatiotemporal processes in weather and climate models.

Plastic photodegradation under simulated marine conditions.

Ocean plastic pollution is a problem of increasing magnitude; yet, the amount of plastic at the sea surface is much lower than expected. Solar ultraviolet (UV) radiation can induce photodegradation, but its importance in determining the longevity of floating plastic remains unconstrained. Here, we measured photodegradation rates of different plastic types slightly larger than microplastics (virgin polymers and floating plastic debris) under simulated marine conditions. UV irradiation caused all plastic types to leach dissolved organic carbon, and to a lesser degree carbon dioxide, carbon monoxide, methane, and other hydrocarbon gases. The release of photodegradation products translates to degradation rates of 1.7-2.3 % yr-1 of the tested plastic particles normalized to conditions as found in the subtropical surface ocean. Modelling the accumulation of floating plastic debris, our results show that solar UV radiation could already have degraded 7 to 22 % of all floating plastic that has ever been released to the sea.

Temperature dependence of isotopic fractionation in the CO2 -O2 isotope exchange reaction.

RATIONALE: Oxygen isotope exchange between O2 and CO2 in the presence of heated platinum (Pt) is an established technique for determining the δ17 O value of CO2 . However, there is not yet a consensus on the associated fractionation factors at the steady state. METHODS: We determined experimentally the steady-state α17 and α18 fractionation factors for Pt-catalyzed CO2 -O2 oxygen isotope exchange at temperatures ranging from 500 to 1200°C. For comparison, the theoretical α18 equilibrium exchange values reported by Richet et al. (1997) have been updated using the direct sum method for CO2 and the corresponding α17 values were determined. Finally, we examined whether the steady-state fractionation factors depend on the isotopic composition of the reactants, by using CO2 and O2 differing in δ18 O value from -66 ‰ to +4 ‰. RESULTS: The experimentally determined steady-state fractionation factors α17 and α18 are lower than those obtained from the updated theoretical calculations (of CO2 -O2 isotope exchange under equilibrium conditions) by 0.0024 ± 0.0001 and 0.0048 ± 0.0002, respectively. The offset is not due to scale incompatibilities between isotope measurements of O2 and CO2 nor to the neglect of non-Born-Oppenheimer effects in the calculations. There is a crossover temperature at which enrichment in the minor isotopes switches from CO2 to O2 . The direct sum evaluation yields a θ value of ~0.54, i.e. higher than the canonical range maximum for a mass-dependent fractionation process. CONCLUSIONS: Updated theoretical values of α18 for equilibrium isotope exchange are lower than those derived from previous work by Richet et al. (1997). The direct sum evaluation for CO2 yields θ values higher than the canonical range maximum for mass-dependent fractionation processes. This demonstrates the need to include anharmonic effects in the calculation and definition of mass-dependent fractionation processes for poly-atomic molecules. The discrepancy between the theory and the experimental α17 and α18 values may be due to thermal diffusion associated with the temperature gradient in the reactor.

Nanoplastics measurements in Northern and Southern polar ice.

It has been established that various anthropogenic contaminants have already reached all the world's pristine locations, including the polar regions. While some of those contaminants, such as lead and soot, are decreasing in the environment, thanks to international regulations, other novel contaminants emerge. Plastic pollution has been shown as a durable novel pollutant, and, since recently, smaller and smaller plastics particles have been identified in various environments (air, water and soil). Considerable research already exists measuring the plastics in the 5 mm to micrometre size range (microplastics). However, far less is known about the plastics debris that fragmented to the sub-micrometre size (nanoplastics). As these small particles are light, it is expected that they have already reached the most remote places on Earth, e.g. transported across the globe by air movement. In this work, we used a novel method based on Thermal Desorption - Proton Transfer Reaction - Mass Spectrometry (TD-PTR-MS) to detect and measure nanoplastics of different types in the water sampled from a Greenland firn core (T2015-A5) and a sea ice core from Antarctica. We identify polyethylene (PE), polypropylene (PP), polyethylene terephthalate (PET), polystyrene (PS), polyvinyl chloride (PVC), and Tire wear nanoparticles in the 14 m deep Greenland firn core and PE, PP and PET in sea ice from Antarctica. Nanoplastics mass concentrations were on average 13.2 ng/mL for Greenland firn samples and 52.3 ng/mL for Antarctic sea ice. We further discuss the possible sources of nanoplastics that we found at these remote locations, which likely involve complex processes of plastic circulation (emission from both land and sea surface, atmospheric and marine circulation).

Observational Evidence of Large Contribution from Primary Sources for Carbon Monoxide in the South Asian Outflow.

South Asian air is among the most polluted in the world, causing premature death of millions and asserting a strong perturbation of the regional climate. A central component is carbon monoxide (CO), which is a key modulator of the oxidizing capacity of the atmosphere and a potent indirect greenhouse gas. While CO concentrations are declining elsewhere, South Asia exhibits an increasing trend for unresolved reasons. In this paper, we use dual-isotope (δ13C and δ18O) fingerprinting of CO intercepted in the South Asian outflow to constrain the relative contributions from primary and secondary CO sources. Results show that combustion-derived primary sources dominate the wintertime continental CO fingerprint (fprimary ∼ 79 ± 4%), significantly higher than the global estimate (fprimary ∼ 55 ± 5%). Satellite-based inventory estimates match isotope-constrained fprimary-CO, suggesting observational convergence in source characterization and a prospect for model-observation reconciliation. This "ground-truthing" emphasizes the pressing need to mitigate incomplete combustion activities for climate/air quality benefits in South Asia.

Nanoplastics transport to the remote, high-altitude Alps.

Plastic materials are increasingly produced worldwide with a total estimated production of >8300 million tonnes to date, of which 60% was discarded. In the environment, plastics fragment into smaller particles, e.g. microplastics (size < 5 mm), and further weathering leads to the formation of functionally different contaminants - nanoplastics (size <1 µm). Nanoplastics are believed to have entirely different physical (e.g. transport), chemical (e.g. functional groups at the surface) and biological (passing the cell membrane, toxicity) properties compared to the micro- and macroplastics, yet, their measurement in the environmental samples is seldom available. Here, we present measurements of nanoplastics mass concentration and calculated the deposition at the pristine high-altitude Alpine Sonnblick observatory (3106 MASL), during the 1.5 month campaigh in late winter 2017. The average nanoplastics concentration was 46.5 ng/mL of melted surface snow. The main polymer types of nanoplastics observed for this site were polypropylene (PP) and polyethylene terephthalate (PET). We measured significantly higher concentrations in the dry sampling periods for PET (p < 0.002) but not for PP, which indicates that dry deposition may be the preferential pathway for PET leading to a gradual accumulation on the snow surfaces during dry periods. Air transport modelling indicates regional and long-range transport of nanoplastics, originating preferentially from European urban areas. The mean deposition rate was 42 (+32/-25) kg km-2 year-1. Thus more than 2 × 1011 nanoplastics particles are deposited per square meter of surface snow each week of the observed period, even at this remote location, which raises significant toxicological concerns.

Leaf scale quantification of the effect of photosynthetic gas exchange on Δ47 of CO2.

The clumped isotope composition (Δ47, the anomaly of the mass 47 isotopologue relative to the abundance expected from a random isotope distribution) of CO2 has been suggested as an additional tracer for gross CO2 fluxes. However, the effect of photosynthetic gas exchange on Δ47 has not been directly determined and two indirect/conceptual studies reported contradicting results. In this study, we quantify the effect of photosynthetic gas exchange on Δ47 of CO2 using leaf cuvette experiments with one C4 and two C3 plants. The experimental results are supported by calculations with a leaf cuvette model. Our results demonstrate the important roles of the Δ47 value of CO2 entering the leaf, kinetic fractionation as CO2 diffuses into, and out of the leaf and CO2-H2O isotope exchange with leaf water. We experimentally confirm the previously suggested dependence of Δ47 of CO2 in the air surrounding a leaf on the stomatal conductance and back-diffusion flux. Gas exchange can enrich or deplete the Δ47 of CO2 depending on the Δ47 of CO2 entering the leaf and the fraction of CO2 exchanged with leaf water and diffused back to the atmosphere, but under typical ambient conditions, it will lead to a decrease in Δ47.

A GC-IRMS method for measuring sulfur isotope ratios of carbonyl sulfide from small air samples.

A new system was developed for measuring sulfur isotopes δ 33S and δ 34S from atmospheric carbonyl sulfide (COS) on small air samples of several liters, using pre-concentration and gas chromatography - isotope ratio mass spectrometry (GC-IRMS). Measurements of COS isotopes provide a tool for quantifying the COS budget, which will help towards better understanding climate feedback mechanisms. For a 4 liter sample at ambient COS mixing ratio, ~500 parts per trillion (ppt), we obtain a reproducibility error of 2.1 ‰ for δ 33S and 0.4 ‰ for δ 34S. After applying corrections, the uncertainty for an individual ambient air sample measurement is 2.5 ‰ for δ 33S and 0.9 ‰ for δ 34S. The ability to measure small samples allows application to a global-scale sampling program with limited logistical effort. To illustrate the application of this newly developed system, we present a timeseries of ambient air measurements, during the fall and winter of 2020 and 2021 in Utrecht, the Netherlands. The observed background values were δ 33S = 1.0 ± 3.4 ‰ and δ 34S = 15.5 ± 0.8 ‰ (VCDT). The maximum observed COS mixing ratios was only 620 ppt. This, in combination with the relatively high δ 34S suggests that the Netherlands receives little COS-containing anthropogenic emissions. We observed a change in COS mixing ratio and δ 34S with different air mass origin, as modelled with HYSPLIT backward trajectory analyses. An increase of 40 ppt in mean COS mixing ratio was observed between fall and winter, which is consistent with the expected seasonal cycle in the Netherlands. Additionally, we present the results of samples from a highway tunnel to characterize vehicle COS emissions and isotopic composition. The vehicle emissions were small, with COS/CO 2 being 0.4 ppt/ppm; the isotopic signatures are depleted relatively to background atmospheric COS.

Publisher Correction: East Siberian Arctic inland waters emit mostly contemporary carbon.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

East Siberian Arctic inland waters emit mostly contemporary carbon.

Inland waters (rivers, lakes and ponds) are important conduits for the emission of terrestrial carbon in Arctic permafrost landscapes. These emissions are driven by turnover of contemporary terrestrial carbon and additional pre-aged (Holocene and late-Pleistocene) carbon released from thawing permafrost soils, but the magnitude of these source contributions to total inland water carbon fluxes remains unknown. Here we present unique simultaneous radiocarbon age measurements of inland water CO2, CH4 and dissolved and particulate organic carbon in northeast Siberia during summer. We show that >80% of total inland water carbon was contemporary in age, but pre-aged carbon contributed >50% at sites strongly affected by permafrost thaw. CO2 and CH4 were younger than dissolved and particulate organic carbon, suggesting emissions were primarily fuelled by contemporary carbon decomposition. Our findings reveal that inland water carbon emissions from permafrost landscapes may be more sensitive to changes in contemporary carbon turnover than the release of pre-aged carbon from thawing permafrost.

Micro- and Nanoplastics in Alpine Snow: A New Method for Chemical Identification and (Semi)Quantification in the Nanogram Range.

We present a new method for chemical characterization of micro- and nanoplastics based on thermal desorption-proton transfer reaction-mass spectrometry. The detection limit for polystyrene (PS) obtained is <1 ng of the compound present in a sample, which results in 100 times better sensitivity than those of previously reported by other methods. This allows us to use small volumes of samples (1 mL) and to carry out experiments without a preconcentration step. Unique features in the high-resolution mass spectrum of different plastic polymers make this approach suitable for fingerprinting, even when the samples contain mixtures of other organic compounds. Accordingly, we got a positive fingerprint of PS when just 10 ng of the polymer was present within the dissolved organic matter of snow. Multiple types of microplastics (polyethylene terephthalate (PET), polyvinyl chloride, and polypropylene carbonate), were identified in a snowpit from the Austrian Alps; however, only PET was detected in the nanometer range for both snowpit and surface snow samples. This is in accordance with other publications showing that the dominant form of airborne microplastics is PET fibers. The presence of nanoplastics in high-altitude snow indicates airborne transport of plastic pollution with environmental and health consequences yet to be understood.

Global 3-D Simulations of the Triple Oxygen Isotope Signature Δ17O in Atmospheric CO2.

The triple oxygen isotope signature Δ17O in atmospheric CO2, also known as its "17O excess," has been proposed as a tracer for gross primary production (the gross uptake of CO2 by vegetation through photosynthesis). We present the first global 3-D model simulations for Δ17O in atmospheric CO2 together with a detailed model description and sensitivity analyses. In our 3-D model framework we include the stratospheric source of Δ17O in CO2 and the surface sinks from vegetation, soils, ocean, biomass burning, and fossil fuel combustion. The effect of oxidation of atmospheric CO on Δ17O in CO2 is also included in our model. We estimate that the global mean Δ17O (defined as Δ 17 O = ln ( δ 17 O + 1 ) - λ RL · ln ( δ 18 O + 1 ) with λ RL = 0.5229) of CO2 in the lowest 500 m of the atmosphere is 39.6 per meg, which is ∼20 per meg lower than estimates from existing box models. We compare our model results with a measured stratospheric Δ17O in CO2 profile from Sodankylä (Finland), which shows good agreement. In addition, we compare our model results with tropospheric measurements of Δ17O in CO2 from Göttingen (Germany) and Taipei (Taiwan), which shows some agreement but we also find substantial discrepancies that are subsequently discussed. Finally, we show model results for Zotino (Russia), Mauna Loa (United States), Manaus (Brazil), and South Pole, which we propose as possible locations for future measurements of Δ17O in tropospheric CO2 that can help to further increase our understanding of the global budget of Δ17O in atmospheric CO2.

Determination of the triple oxygen and carbon isotopic composition of CO2 from atomic ion fragments formed in the ion source of the 253 Ultra high-resolution isotope ratio mass spectrometer.

RATIONALE: Determination of δ17 O values directly from CO2 with traditional gas source isotope ratio mass spectrometry is not possible due to isobaric interference of 13 C16 O16 O on 12 C17 O16 O. The methods developed so far use either chemical conversion or isotope equilibration to determine the δ17 O value of CO2 . In addition, δ13 C measurements require correction for the interference from 12 C17 O16 O on 13 C16 O16 O since it is not possible to resolve the two isotopologues. METHODS: We present a technique to determine the δ17 O, δ18 O and δ13 C values of CO2 from the fragment ions that are formed upon electron ionization in the ion source of the Thermo Scientific 253 Ultra high-resolution isotope ratio mass spectrometer (hereafter 253 Ultra). The new technique is compared with the CO2 -O2 exchange method and the 17 O-correction algorithm for δ17 O and δ13 C values, respectively. RESULTS: The scale contractions for δ13 C and δ18 O values are slightly larger for fragment ion measurements than for molecular ion measurements. The δ17 O and Δ17 O values of CO2 can be measured on the 17 O+ fragment with an internal error that is a factor 1-2 above the counting statistics limit. The ultimate precision depends on the signal intensity and on the total time that the 17 O+ beam is monitored; a precision of 14 ppm (parts per million) (standard error of the mean) was achieved in 20 hours at the University of Göttingen. The Δ17 O measurements with the O-fragment method agree with the CO2 -O2 exchange method over a range of Δ17 O values of -0.3 to +0.7‰. CONCLUSIONS: Isotope measurements on atom fragment ions of CO2 can be used as an alternative method to determine the carbon and oxygen isotopic composition of CO2 without chemical processing or corrections for mass interferences.

Measurement of 18 O18 O and 17 O18 O in atmospheric O2 using the 253 Ultra mass spectrometer and applications to stratospheric and tropospheric air samples.

RATIONALE: The doubly substituted isotopologues (e.g., 18 O18 O, 17 O18 O) in atmospheric O2 are potential tracers for ozone photochemistry and atmospheric temperatures. Their low abundances and isobaric interference are the major analytical challenges. The 253 Ultra high-resolution stable isotope ratio mass spectrometer is suitable for resolving isobaric interferences. METHODS: O2 from air is purified using gas chromatography on a packed column filled with molecular sieve 5 Å and cooled to -78°C. The δ17 O, δ18 O, Δ17 O, Δ35 and Δ36 values are measured on the extracted O2 with the 253 Ultra at medium mass resolution (M/ΔM ~10000) using Faraday detectors for the singly substituted isotopologues and ion counters for the doubly substituted isotopologues. RESULTS: Interferences from isobars, mainly 35 Cl for 17 O18 O and H35 Cl and 36 Ar for 18 O18 O, are sufficiently resolved to enable high-precision determination of Δ35 and Δ36 . The Δ35 and Δ36 values of O2 after photochemical isotope equilibration at -63°C and heating to 850°C agree with the theoretical prediction. The stratospheric Δ35 and Δ36 values are close to isotopic equilibrium at the ambient temperatures. However, the values for tropospheric O2 differ from those expected at equilibrium. CONCLUSIONS: The 253 Ultra allows interference-free clumped isotope measurements of O2 at medium mass resolution. The Δ35 and Δ36 signatures in atmospheric O2 are mainly governed by O3 photochemistry, temperature and atmospheric transport. Tropospheric O2 is isotopically well mixed and retains a significant stratospheric signature.

H2 clumped isotope measurements at natural isotopic abundances.

RATIONALE: Molecular hydrogen (H2 ) is an important gas for atmospheric chemistry, and an indirect greenhouse gas due to its reaction with OH. The isotopic composition of H2 (δD) has been used to investigate its atmospheric budget; here we add a new observable, the clumped isotopic signature ΔDD, to the tools that can be used to study the global cycle of H2 . METHODS: A method for determining ΔDD in H2 was developed using the high-resolution MAT 253-Ultra isotope ratio mass spectrometer (Thermo Fisher). The HH, HD and DD abundances are quantified at medium resolution (M/ΔM ≈ 6000), which is sufficient for HD+ and DD+ to be distinguished from H3 + and H2 D+ , respectively. The method involves sequential measurement of isotopologues, and DD is measured using an ion counter. For verification, catalytic ΔDD equilibration experiments were performed at temperatures of up to 850°C. RESULTS: The typical precision obtained for ΔDD is 2-6‰, close to the theoretical counting statistics limit, and adequate for detecting the expected natural variations. Compatibility and medium-term reproducibility are consistent with the precision values. The method was validated using temperature equilibration experiments, which showed a dependence of ΔDD on temperature as expected form theoretical calculations. CONCLUSIONS: We have established a method for determining ΔDD in H2 at natural isotopic abundances, with a precision that is adequate for observing the expected variations in atmospheric and other natural H2 . This method opens the road to new research on the natural H2 cycle.