Widespread six degrees Celsius cooling on land during the Last Glacial Maximum.

The magnitude of global cooling during the Last Glacial Maximum (LGM, the coldest multimillennial interval of the last glacial period) is an important constraint for evaluating estimates of Earth's climate sensitivity1,2. Reliable LGM temperatures come from high-latitude ice cores3,4, but substantial disagreement exists between proxy records in the low latitudes1,5-8, where quantitative low-elevation records on land are scarce. Filling this data gap, noble gases in ancient groundwater record past land surface temperatures through a direct physical relationship that is rooted in their temperature-dependent solubility in water9,10. Dissolved noble gases are suitable tracers of LGM temperature because of their complete insensitivity to biological and chemical processes and the ubiquity of LGM-aged groundwater around the globe11,12. However, although several individual noble gas studies have found substantial tropical LGM cooling13-16, they have used different methodologies and provide limited spatial coverage. Here we use noble gases in groundwater to show that the low-altitude, low-to-mid-latitude land surface (45 degrees south to 35 degrees north) cooled by 5.8 ± 0.6 degrees Celsius (mean ± 95% confidence interval) during the LGM. Our analysis includes four decades of groundwater noble gas data from six continents, along with new records from the tropics, all of which were interpreted using the same physical framework. Our land-based result broadly supports a recent reconstruction based on marine proxy data assimilation1 that suggested greater climate sensitivity than previous estimates5-7.

Dissolved gases in the deep North Atlantic track ocean ventilation processes.

Gas exchange between the atmosphere and ocean interior profoundly impacts global climate and biogeochemistry. However, our understanding of the relevant physical processes remains limited by a scarcity of direct observations. Dissolved noble gases in the deep ocean are powerful tracers of physical air-sea interaction due to their chemical and biological inertness, yet their isotope ratios have remained underexplored. Here, we present high-precision noble gas isotope and elemental ratios from the deep North Atlantic (~32°N, 64°W) to evaluate gas exchange parameterizations using an ocean circulation model. The unprecedented precision of these data reveal deep-ocean undersaturation of heavy noble gases and isotopes resulting from cooling-driven air-to-sea gas transport associated with deep convection in the northern high latitudes. Our data also imply an underappreciated and large role for bubble-mediated gas exchange in the global air-sea transfer of sparingly soluble gases, including O2, N2, and SF6. Using noble gases to validate the physical representation of air-sea gas exchange in a model also provides a unique opportunity to distinguish physical from biogeochemical signals. As a case study, we compare dissolved N2/Ar measurements in the deep North Atlantic to physics-only model predictions, revealing excess N2 from benthic denitrification in older deep waters (below 2.9 km). These data indicate that the rate of fixed N removal in the deep Northeastern Atlantic is at least three times higher than the global deep-ocean mean, suggesting tight coupling with organic carbon export and raising potential future implications for the marine N cycle.

High 3He/4He in central Panama reveals a distal connection to the Galápagos plume.

It is well established that mantle plumes are the main conduits for upwelling geochemically enriched material from Earth's deep interior. The fashion and extent to which lateral flow processes at shallow depths may disperse enriched mantle material far (>1,000 km) from vertical plume conduits, however, remain poorly constrained. Here, we report He and C isotope data from 65 hydrothermal fluids from the southern Central America Margin (CAM) which reveal strikingly high 3He/4He (up to 8.9RA) in low-temperature (≤50 °C) geothermal springs of central Panama that are not associated with active volcanism. Following radiogenic correction, these data imply a mantle source 3He/4He >10.3RA (and potentially up to 26RA, similar to Galápagos hotspot lavas) markedly greater than the upper mantle range (8 ± 1RA). Lava geochemistry (Pb isotopes, Nb/U, and Ce/Pb) and geophysical constraints show that high 3He/4He values in central Panama are likely derived from the infiltration of a Galápagos plume-like mantle through a slab window that opened ∼8 Mya. Two potential transport mechanisms can explain the connection between the Galápagos plume and the slab window: 1) sublithospheric transport of Galápagos plume material channeled by lithosphere thinning along the Panama Fracture Zone or 2) active upwelling of Galápagos plume material blown by a "mantle wind" toward the CAM. We present a model of global mantle flow that supports the second mechanism, whereby most of the eastward transport of Galápagos plume material occurs in the shallow asthenosphere. These findings underscore the potential for lateral mantle flow to transport mantle geochemical heterogeneities thousands of kilometers away from plume conduits.

Solubility Equilibrium Isotope Effects of Noble Gases in Water: Theory and Observations.

The abundance and isotopic composition of noble gases dissolved in water have many applications in the geosciences. In recent years, new analytical techniques have opened the door to the use of high-precision measurements of noble gas isotopes as tracers for groundwater hydrology, oceanography, mantle geochemistry, and paleoclimatology. These analytical advances have brought about new measurements of solubility equilibrium isotope effects (SEIEs) in water (i.e., the relative solubilities of noble gas isotopes) and their sensitivities to the temperature and salinity. Here, we carry out a suite of classical molecular dynamics (MD) simulations and employ the theoretical method of quantum correction to estimate SEIEs for comparison with experimental observations. We find that classical MD simulations can accurately predict SEIEs for the isotopes of Ar, Kr, and Xe to order 0.01‰, on the scale of analytical uncertainty. However, MD simulations consistently overpredict the SEIEs of Ne and He by up to 40% of observed values. We carry out sensitivity tests at different temperatures, salinities, and pressures and employ different sets of interatomic potential parameters and water models. For all noble gas isotopes, the TIP4P water model is found to reproduce observed SEIEs more accurately than the SPC/E and TIP4P/ice models. Classical MD simulations also accurately capture the sign and approximate magnitude of temperature and salinity sensitivities of SEIEs for heavy noble gases. We find that experimental and modeled SEIEs generally follow an inverse-square mass dependence, which implies that the mean-square force experienced by a noble gas atom within a solvation shell is similar for all noble gases. This inverse-square mass proportionality is nearly exact for Ar, Kr, and Xe isotopes, but He and Ne exhibit a slightly weaker mass dependence. We hypothesize that the apparent dichotomy between He-Ne and Ar-Kr-Xe SEIEs may result from atomic size differences, whereby the smaller noble gases are more likely to spontaneously fit within cavities of water without breaking water-water H-bonds, thereby experiencing softer collisions during translation within a solvation shell. We further speculate that the overprediction of simulated He and Ne SEIEs may result from the neglection of higher-order quantum corrections or the overly stiff representation of van der Waals repulsion by the widely used Lennard-Jones 6-12 potential model. We suggest that new measurements of SEIEs of heavy and light noble gases may represent a novel set of constraints with which to refine hydrophobic solvation theories and optimize the set of interatomic potential models used in MD simulations of water and noble gases.

Terrestrial amplification of past, present, and future climate change.

Terrestrial amplification (TA) of land warming relative to oceans is apparent in recent climatic observations. TA results from land-sea coupling of moisture and heat and is therefore important for predicting future warming and water availability. However, the theoretical basis for TA has never been tested outside the short instrumental period, and the spatial pattern and amplitude of TA remain uncertain. Here, we investigate TA during the Last Glacial Maximum (LGM; ~20 thousand years) in the low latitudes, where the theory is most applicable. We find remarkable consistency between paleotemperature proxies, theory, and climate model simulations of both LGM and future climates. Paleoclimate data thus provide crucial new support for TA, refining the range of future low-latitude, low-elevation TA to [Formula: see text] (95% confidence interval), i.e., land warming ~40% more than oceans. The observed data model theory agreement helps reconcile LGM marine and terrestrial paleotemperature proxies, with implications for equilibrium climate sensitivity.

Ultrahigh-precision noble gas isotope analyses reveal pervasive subsurface fractionation in hydrothermal systems.

Mantle-derived noble gases in volcanic gases are powerful tracers of terrestrial volatile evolution, as they contain mixtures of both primordial (from Earth's accretion) and secondary (e.g., radiogenic) isotope signals that characterize the composition of deep Earth. However, volcanic gases emitted through subaerial hydrothermal systems also contain contributions from shallow reservoirs (groundwater, crust, atmosphere). Deconvolving deep and shallow source signals is critical for robust interpretations of mantle-derived signals. Here, we use a novel dynamic mass spectrometry technique to measure argon, krypton, and xenon isotopes in volcanic gas with ultrahigh precision. Data from Iceland, Germany, United States (Yellowstone, Salton Sea), Costa Rica, and Chile show that subsurface isotope fractionation within hydrothermal systems is a globally pervasive and previously unrecognized process causing substantial nonradiogenic Ar-Kr-Xe isotope variations. Quantitatively accounting for this process is vital for accurately interpreting mantle-derived volatile (e.g., noble gas and nitrogen) signals, with profound implications for our understanding of terrestrial volatile evolution.

Groundwater residence time estimates obscured by anthropogenic carbonate.

Groundwater is an important source of drinking and irrigation water. Dating groundwater informs its vulnerability to contamination and aids in calibrating flow models. Here, we report measurements of multiple age tracers (14C, 3H, 39Ar, and 85Kr) and parameters relevant to dissolved inorganic carbon (DIC) from 17 wells in California's San Joaquin Valley (SJV), an agricultural region that is heavily reliant on groundwater. We find evidence for a major mid-20th century shift in groundwater DIC input from mostly closed- to mostly open-system carbonate dissolution, which we suggest is driven by input of anthropogenic carbonate soil amendments. Crucially, enhanced open-system dissolution, in which DIC equilibrates with soil CO2, fundamentally affects the initial 14C activity of recently recharged groundwater. Conventional 14C dating of deeper SJV groundwater, assuming an open system, substantially overestimates residence time and thereby underestimates susceptibility to modern contamination. Because carbonate soil amendments are ubiquitous, other groundwater-reliant agricultural regions may be similarly affected.

Deglacial water-table decline in Southern California recorded by noble gas isotopes.

Constraining the magnitude of past hydrological change may improve understanding and predictions of future shifts in water availability. Here we demonstrate that water-table depth, a sensitive indicator of hydroclimate, can be quantitatively reconstructed using Kr and Xe isotopes in groundwater. We present the first-ever measurements of these dissolved noble gas isotopes in groundwater at high precision (≤0.005‰ amu-1; 1σ), which reveal depth-proportional signals set by gravitational settling in soil air at the time of recharge. Analyses of California groundwater successfully reproduce modern groundwater levels and indicate a 17.9 ± 1.3 m (±1 SE) decline in water-table depth in Southern California during the last deglaciation. This hydroclimatic transition from the wetter glacial period to more arid Holocene accompanies a surface warming of 6.2 ± 0.6 °C (±1 SE). This new hydroclimate proxy builds upon an existing paleo-temperature application of noble gases and may identify regions prone to future hydrological change.

Helium, inorganic and organic carbon isotopes of fluids and gases across the Costa Rica convergent margin.

In 2017, fluid and gas samples were collected across the Costa Rican Arc. He and Ne isotopes, C isotopes as well as total organic and inorganic carbon concentrations were measured. The samples (n = 24) from 2017 are accompanied by (n = 17) samples collected in 2008, 2010 and 2012. He-isotopes ranged from arc-like (6.8 RA) to crustal (0.5 RA). Measured dissolved inorganic carbon (DIC) δ13CVPDB values varied from 3.55 to -21.57‰, with dissolved organic carbon (DOC) following the trends of DIC. Gas phase CO2 only occurs within ~20 km of the arc; δ13CVPDB values varied from -0.84 to -5.23‰. Onsite, pH, conductivity, temperature and dissolved oxygen (DO) were measured; pH ranged from 0.9-10.0, conductivity from 200-91,900 µS/cm, temperatures from 23-89 °C and DO from 2-84%. Data were used to develop a model which suggests that ~91 ± 4.0% of carbon released from the slab/mantle beneath the Costa Rican forearc is sequestered within the crust by calcite deposition with an additional 3.3 ± 1.3% incorporated into autotrophic biomass.