skip to main content


Title: Enhancing Understanding of the Hydrological Cycle via Pairing of Process‐Oriented and Isotope Ratio Tracers
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.

 
more » « less
Award ID(s):
1954660 1954663
NSF-PAR ID:
10447447
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
13
Issue:
10
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Abstract

    The isotopic composition of precipitation is used to trace water cycling and climate change, but interpretations of the environmental information recorded in central Andean precipitation isotope ratios are hindered by a lack of multi‐year records, poor spatial distribution of observations, and a predominant focus on Rayleigh distillation. To better understand isotopic variability in central Andean precipitation, we present a three‐year record of semimonthly δ18Opand δ2Hpvalues from 15 stations in southern Peru and triple oxygen isotope data, expressed as ∆′17Op, from 32 precipitation samples. Consistent with previous work, we find that elevation correlates negatively with δ18Opand that seasonal δ18Opvariations are related to upstream rainout and local convection. Spatial δ18Opvariations and atmospheric back trajectories show that both eastern‐ and western‐derived air masses bring precipitation to southern Peru. Seasonal d‐excesspcycles record moisture recycling and relative humidity at remote moisture sources, and both d‐excesspand ∆′17Opclearly differentiate evaporated and non‐evaporated samples. These results begin to establish the natural range of unevaporated ∆′17Opvalues in the central Andes and set the foundation for future paleoclimate and paleoaltimetry studies in the region. This study highlights the hydrologic understanding that comes from a combination of δ18Op, d‐excessp, and ∆′17Opdata and helps identify the evaporation, recycling, and rainout processes that drive water cycling in the central Andes.

     
    more » « less
  4. Atmospheric rivers (ARs), and frontal systems more broadly, tend to exhibit prominent “V” shapes in time series of stable isotopes in precipitation. Despite the magnitude and widespread nature of these “V” shapes, debate persists as to whether these shifts are driven by changes in the degree of rainout, which we determine using the Rayleigh distillation of stable isotopes, or by post-condensation processes such as below-cloud evaporation and equilibrium isotope exchange between hydrometeors and surrounding vapor. Here, we present paired precipitation and water vapor isotope time series records from the 5–7 March 2016, AR in Bodega Bay, CA. The stable isotope composition of surface vapor along with independent meteorological constraints such as temperature and relative humidity reveal that rainout and post-condensation processes dominate during different portions of the event. We find that Rayleigh distillation controls during peak AR conditions (with peak rainout of 55%) while post-condensation processes have their greatest effect during periods of decreased precipitation on the margins of the event. These results and analyses inform critical questions regarding the temporal evolution of AR events and the physical processes that control them at local scales. 
    more » « less
  5. 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