skip to main content


Title: Isotopes in the Water Cycle: Regional- to Global-Scale Patterns and Applications
Stable isotope ratios of hydrogen and oxygen have been applied to water cycle research for over 60 years. Over the past two decades, however, new data, data compilations, and quantitative methods have supported the application of isotopic data to address large-scale water cycle problems. Recent results have demonstrated the impact of climate variation on atmospheric water cycling, provided constraints on continental- to global-scale land-atmosphere water vapor fluxes, revealed biases in the sources of runoff in hydrological models, and illustrated regional patterns of water use and management by people. In the past decade, global isotopic observations have spurred new debate over the role of soils in the water cycle, with potential to impact both ecological and hydrological theory. Many components of the water cycle remain underrepresented in isotopic databases. Increasing accessibility of analyses and improved platforms for data sharing will refine and grow the breadth of these contributions in the future. ▪ Isotope ratios in water integrate information on hydrological processes over scales from cities to the globe. ▪ Tracing water with isotopes helps reveal the processes that govern variability in the water cycle and may govern future global changes. ▪ Improvements in instrumentation, data sharing, and quantitative analysis have advanced isotopic water cycle science over the past 20 years.  more » « less
Award ID(s):
1759730 1802880 1565128 1241286
NSF-PAR ID:
10107109
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Annual Review of Earth and Planetary Sciences
Volume:
47
Issue:
1
ISSN:
0084-6597
Page Range / eLocation ID:
453 to 479
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract. Reconstructions of global hydroclimate during the Common Era (CE; the past ∼2000 years) are important for providing context for current and future global environmental change. Stable isotope ratios in water are quantitative indicators of hydroclimate on regional to global scales, and these signals are encoded in a wide range of natural geologic archives. Here we present the Iso2k database, a global compilation of previously published datasets from a variety of natural archives that record the stable oxygen (δ18O) or hydrogen (δ2H) isotopic compositions of environmental waters, which reflect hydroclimate changes over the CE. The Iso2k database contains 759 isotope records from the terrestrial and marine realms, including glacier and ground ice (210); speleothems (68); corals, sclerosponges, and mollusks (143); wood (81); lake sediments and other terrestrial sediments (e.g., loess) (158); and marine sediments (99). Individual datasets have temporal resolutions ranging from sub-annual to centennial and include chronological data where available. A fundamental feature of the database is its comprehensive metadata, which will assist both experts and nonexperts in the interpretation of each record and in data synthesis. Key metadata fields have standardized vocabularies to facilitate comparisons across diversearchives and with climate-model-simulated fields. This is the firstglobal-scale collection of water isotope proxy records from multiple typesof geological and biological archives. It is suitable for evaluatinghydroclimate processes through time and space using large-scale synthesis,model–data intercomparison and (paleo)data assimilation. The Iso2k databaseis available for download at https://doi.org/10.25921/57j8-vs18 (Konecky and McKay, 2020) and is also accessible via the NOAA/WDS Paleo Datalanding page: https://www.ncdc.noaa.gov/paleo/study/29593 (last access: 30 July 2020). 
    more » « less
  2. This dataset contains monthly average output files from the iCAM6 simulations used in the manuscript "Enhancing understanding of the hydrological cycle via pairing of process-oriented and isotope ratio tracers," in review at the Journal of Advances in Modeling Earth Systems. A file corresponding to each of the tagged and isotopic variables used in this manuscript is included. Files are at 0.9° latitude x 1.25° longitude, and are in NetCDF format. Data from two simulations are included: 1) a simulation where the atmospheric model was "nudged" to ERA5 wind and surface pressure fields, by adding an additional tendency (see section 3.1 of associated manuscript), and 2) a simulation where the atmospheric state was allowed to freely evolve, using only boundary conditions imposed at the surface and top of atmosphere. Specific information about each of the variables provided is located in the "usage notes" section below. Associated article abstract: The hydrologic cycle couples the Earth's energy and carbon budgets through evaporation, moisture transport, and precipitation. Despite a wealth of observations and models, fundamental limitations remain in our capacity to deduce even the most basic properties of the hydrological cycle, including the spatial pattern of the residence time (RT) of water in the atmosphere and the mean distance traveled from evaporation sources to precipitation sinks. Meanwhile, geochemical tracers such as stable water isotope ratios provide a tool to probe hydrological processes, yet their interpretation remains equivocal despite several decades of use. As a result, there is a need for new mechanistic tools that link variations in water isotope ratios to underlying hydrological processes. Here we present a new suite of “process-oriented tags,” which we use to explicitly trace hydrological processes within the isotopically enabled Community Atmosphere Model, version 6 (iCAM6). Using these tags, we test the hypotheses that precipitation isotope ratios respond to parcel rainout, variations in atmospheric RT, and preserve information regarding meteorological conditions during evaporation. We present results for a historical simulation from 1980 to 2004, forced with winds from the ERA5 reanalysis. We find strong evidence that precipitation isotope ratios record information about atmospheric rainout and meteorological conditions during evaporation, but little evidence that precipitation isotope ratios vary with water vapor RT. These new tracer methods will enable more robust linkages between observations of isotope ratios in the modern hydrologic cycle or proxies of past terrestrial environments and the environmental processes underlying these observations.   Details about the simulation setup can be found in section 3 of the associated open-source manuscript, "Enhancing understanding of the hydrological cycle via pairing of process‐oriented and isotope ratio tracers." In brief, we conducted two simulations of the atmosphere from 1980-2004 using the isotope-enabled version of the Community Atmosphere Model 6 (iCAM6) at 0.9x1.25° horizontal resolution, and with 30 vertical hybrid layers spanning from the surface to ~3 hPa. In the first simulation, wind and surface pressure fields were "nudged" toward the ERA5 reanalysis dataset by adding a nudging tendency, preventing the model from diverging from observed/reanalysis wind fields. In the second simulation, no additional nudging tendency was included, and the model was allowed to evolve 'freely' with only boundary conditions provided at the top (e.g., incoming solar radiation) and bottom (e.g., observed sea surface temperatures) of the model. In addition to the isotopic variables, our simulation included a suite of 'process-oriented tracers,' which we describe in section 2 of the manuscript. These variables are meant to track a property of water associated with evaporation, condensation, or atmospheric transport. Metadata are provided about each of the files below; moreover, since the attached files are NetCDF data - this information is also provided with the data files. NetCDF metadata can be accessed using standard tools (e.g., ncdump). Each file has 4 variables: the tagged quantity, and the associated coordinate variables (time, latitude, longitude). The latter three are identical across all files, only the tagged quantity changes. Twelve files are provided for the nudged simulation, and an additional three are provided for the free simulations: Nudged simulation files iCAM6_nudged_1980-2004_mon_RHevap: Mass-weighted mean evaporation source property: RH (%) with respect to surface temperature. iCAM6_nudged_1980-2004_mon_Tevap: Mass-weighted mean evaporation source property: surface temperature in Kelvin iCAM6_nudged_1980-2004_mon_Tcond: Mass-weighted mean condensation property: temperature (K) iCAM6_nudged_1980-2004_mon_columnQ: Total (vertically integrated) precipitable water (kg/m2).  Not a tagged quantity, but necessary to calculate depletion times in section 4.3 (e.g., Fig. 11 and 12). iCAM6_nudged_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_nudged_1980-2004_mon_d18Oevap_0: Mass-weighted mean evaporation source property - d18O of the evaporative flux (e.g., the 'initial' isotope ratio prior to condensation), (‰ VSMOW) iCAM6_nudged_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_nudged_1980-2004_mon_dexevap_0: Mass-weighted mean evaporation source property - deuterium excess of the evaporative flux iCAM6_nudged_1980-2004_mon_lnf: Integrated property - ln(f) calculated from the constant-fractionation d18O tracer (see section 3.2). iCAM6_nudged_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. iCAM6_nudged_1980-2004_mon_residencetime: Mean atmospheric water residence time (in days). iCAM6_nudged_1980-2004_mon_transportdistance: Mean atmospheric water transport distance (in km). Free simulation files iCAM6_free_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_free_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_free_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. 
    more » « less
  3. Abstract

    Because of the pervasive role of water in the Earth system, the relative abundances of stable isotopologues of water are valuable for understanding atmospheric, oceanic, and biospheric processes, and for interpreting paleoclimate proxy reconstructions. Isotopologues are transported by both large‐scale and turbulent flows, and the ratio of heavy to light isotopologues changes due to fractionation that can accompany condensation and evaporation processes. Correctly predicting the isotopic distributions requires resolving the relationships between large‐scale ocean and atmospheric circulation and smaller‐scale hydrological processes, which can be accomplished within a coupled climate modeling framework. Here we present the water isotope‐enabled version of the Community Earth System Model version 1 (iCESM1), which simulates global variations in water isotopic ratios in the atmosphere, land, ocean, and sea ice. In a transient Last Millennium simulation covering the 850–2005 period, iCESM1 correctly captures the late‐twentieth‐century structure of δ18O and δD over the global oceans, with more limited accuracy over land. The relationship between salinity and seawater δ18O is also well represented over the observational period, including interbasin variations. We illustrate the utility of coupled, isotope‐enabled simulations using both Last Millennium simulations and freshwater hosing experiments with iCESM1. Closing the isotopic mass balance between all components of the coupled model provides new confidence in the underlying depiction of the water cycle in CESM, while also highlighting areas where the underlying hydrologic balance can be improved. The iCESM1 is poised to be a vital community resource for ongoing model development with both modern and paleoclimate applications.

     
    more » « less
  4. This Chapter considers triple oxygen isotope variations and their 4 Gyr temporal evolution in bulk siliciclastic sedimentary rocks and in granites. The d18O and D'17O values provide new insights into weathering in the modern and ancient hydrosphere and coeval crustal petrogenesis. We make use of the known geological events and processes that affect the rock cycle: supercontinent assembly and breakup that influence continent-scale and global climate, the fraction of the exposed crust undergoing weathering, and isotopic values of precipitation. New data from a 5000 m Texas drillhole into the Oligocene Frio Formation demonstrate minimal isotopic shifts from mudrocks to shales during diagenesis, mostly related to expulsion of water from smectite-rich loosely cemented sediment and its conversion to illite-rich shale. Inversion of triple oxygen isotope fractionations return isotopic values and temperatures along the hole depth that are more consistent with weathering conditions in the Oligocene and modern North America (d18O = -7 to -15‰, and T of +15 to +45°C) rather than d18O from 8 to 10‰ diagenetic water in the drill hole at 175-195°C. More precise T and d18Owater are obtained where the chemical index of alteration (CIA) based detrital contribution is subtracted from these sediments. Triple oxygen isotopes from suspended sediments in major world rivers record conditions (T and d18Ow) of their watersheds, and not the composition of bedrock because weathering is water-dominated. In parallel, the Chapter presents new analyses of 100 granites, orthogneisses, migmatites, tonalite-trondhjemite-granodiorite (TTG), and large-volume ignimbrites from around the world that range in age from 4 Ga to modern. Most studied granites are orogenic and anatectic in origin and represent large volume remelting/assimilation of shales and other metasediments; the most crustal and high-d18O of these are thus reflect and record the average composition of evolving continental crust. Granites also develop a significant progressive increase in d18O values from 6-7‰ (4-2.5 Ga) to 10-13‰ (~1.8-1.2 Ga) after which d18O stays constant or even decreases. More importantly, we observe a moderate -0.03‰ step-wise decrease in D'17O between 2.1 and 2.5 Ga, which is about half of the step-wise decrease observed in shales over this time interval. We suggest that granites, as well as shales, record the significant advent and greater volumetric appearance of low-D'17O, high-d18O weathering products (shales) altered by meteoric waters upon rapid emergence of large land masses at ~2.4 Ga, although consider alternative interpretations. These weathering products were incorporated into abundant 2.0-1.8 Ga orogens around the world, where upon remelting, they passed their isotopic signature to the granites. We further observe the dichotomy of high-D'17O Archean shales, and unusually low-D'17O Archean granites. We attribute this to greater contribution from shallow crustal hydrothermal contribution to shales in greenstone belts, while granites in the earliest 3.0-4.0 Ga crust and TTGs require involvement of hydrothermal products with lower-D'17O signatures at moderately high-d18O, which we attribute to secondary silicification of their protoliths before partial melting. The Chapter further discusses evolution of the shale record through geologic history and discusses the step-wise change in d18O and D'17O values at Archean/Proterozoic transition. Denser coverage for shales in the past 1 billion years permits investigation of the rocks and their weathering in the last supercontinent cycle, with observed lighter d18O values, characteristic for the mid-Phanerozoic at the initiation of Gondwana breakup. The continuing increase in d18O values of the shales since 4 Ga is interpreted to reflect accumulation of weathering products via shale accretion to continents, as low-density and buoyant shales tend to not subduct back into the mantle. The rock cycle passes triple oxygen isotopic signatures from precipitation to sedimentary, metasedimentary, and finally to anatectic igneous rocks. Continental crust became progressively heavier in d18O, lighter in D'17O due to incremental accumulation of high-d18O sediments in accretionary wedges. Second-order trends in d18O and D'17O are due to supercontinent cycles and glacial episodes. 
    more » « less
  5. Abstract

    The hydrologic cycle is a fundamental component of the climate system with critical societal and ecological relevance. Yet gaps persist in our understanding of water fluxes and their response to increased greenhouse gas forcing. The stable isotope ratios of oxygen and hydrogen in water provide a unique opportunity to evaluate hydrological processes and investigate their role in the variability of the climate system and its sensitivity to change. Water isotopes also form the basis of many paleoclimate proxies in a variety of archives, including ice cores, lake and marine sediments, corals, and speleothems. These records hold most of the available information about past hydrologic variability prior to instrumental observations. Water isotopes thus provide a ‘common currency’ that links paleoclimate archives to modern observations, allowing us to evaluate hydrologic processes and their effects on climate variability on a wide range of time and length scales. Building on previous literature summarizing advancements in water isotopic measurements and modeling and describe water isotopic applications for understanding hydrological processes, this topical review reflects on new insights about climate variability from isotopic studies. We highlight new work and opportunities to enhance our understanding and predictive skill and offer a set of recommendations to advance observational and model-based tools for climate research. Finally, we highlight opportunities to better constrain climate sensitivity and identify anthropogenically-driven hydrologic changes within the inherently noisy background of natural climate variability.

     
    more » « less