Continental-scale isotope hydrology

Providing sustainable sources of fresh water for a growing population of 7 billion people is one of the grand challenges of the 21st century. This dissertation outlines several applications of isotope hydrology to address four previously unknown questions involving surface- and ground-water resources at regional- to continental-spatial scales over contemporary- to millennial-temporal scales. The four chapters in this dissertation investigate (1) the rate of plant transpiration, (2) the seasonality of groundwater recharge, (3) the climate of the last ice age, and (4) the chemistry of Ugandan waters. (1) Chapter one presents a new global compilation of lake water isotopic data, river isotopic data, stand-level transpiration rates, and water use efficiency measurements, and analyzes the newly synthesized data to show that plant transpiration is the largest water flux from Earth\u27s continents, exceeding both physical evaporation and continental runoff. (2) Chapter two presents a new global synthesis of rain, snow and groundwater isotopic compositions, and analyzes the paired precipitation-groundwater dataset to show that the percentage of precipitation that recharges aquifers is at a maximum during the winter (extra-tropics) and wet (tropics) seasons. (3) Chapter three presents a new global compilation of groundwater radiocarbon, tritium, and stable O and H isotopic data, and maps the isotopic shift of meteoric waters since the last ice age. The analysis shows that the majority (~90%) of precipitation during the last ice age had lower 18O/16O and 2H/1H ratios than the modern day, except in some exclusively coastal locations. We also show that current isotope-enabled general circulation models capture some, but not all, spatial variability in ice-age-to-late-Holocene 18O/16O and 2H/1H shifts, providing a new calibration tool that can be used to improve our understanding of glacial climate dynamics. (4) Chapter four presents isotopic and chemical analyses of Ugandan lake, river, rain, and ground water collected during a field expedition led in July of 2013. Analysis of this new dataset reveals new estimates of lake water balances across Uganda

Intensive rainfall recharges tropical groundwaters

Dependence upon groundwater to meet rising agricultural and domestic water needs is expected to increase substantially across the tropics where, by 2050, over half of the world's population is projected to live. Rare, long-term groundwater-level records in the tropics indicate that groundwater recharge occurs disproportionately from heavy rainfalls exceeding a threshold. The ubiquity of this bias in tropical groundwater recharge to intensive precipitation is, however, unknown. By relating available long-term records of stable-isotope ratios of O and H in tropical precipitation (15 sites) to those of local groundwater, we reveal that groundwater recharge in the tropics is near-uniformly (14/15 sites) biased to intensive monthly rainfall, commonly exceeding the ~70th intensity decile. Our results suggest that the intensification of precipitation brought about by global warming favours groundwater replenishment in the tropics. Nevertheless, the processes that transmit intensive rainfall to groundwater systems and enhance the resilience of tropical groundwater storage in a warming world, remain unclear

Stable isotope mass balance of the North American Laurentian Great Lakes

This thesis describes a method for calculating lake evaporation as a proportion of water inputs (E/I) for large surface water bodies, using stable isotope ratios of oxygen (18O/16O) and hydrogen (2H/1H) in water. Evaporation as a proportion of inflow (E/I) is calculated for each Laurentian Great Lake using a new dataset of 516 analyses of δ18O and δ2H in waters sampled from 75 offshore stations during spring and summer of 2007. This work builds on previous approaches by accounting for lake effects on the overlying atmosphere and assuming conservation of both mass and isotopes (18O and 2H) to better constrain evaporation outputs. Results show that E/I ratios are greatest for headwater Lakes Superior and Michigan and lowest for Lakes Erie and Ontario, controlled largely by the magnitude of hydrologic inputs from upstream chain lakes. For Lake Superior, stable isotopes incorporate evaporation over the past century, providing long-term insights to the lake’s hydrology that may be compared to potential changes under a future – expectedly warmer – climate. Uncertainties in isotopically derived E/I are comparable to conventional energy and mass balance uncertainties. Isotope-derived E/I values are lower than conventional energy and mass balance estimates for Lakes Superior and Michigan. The difference between conventional and isotope estimates may be explained by moisture recycling effects. The isotope-based estimates include only evaporated moisture that is also advected from the lake surface, thereby discounting moisture that evaporates and subsequently reprecipitates on the lake surface downwind as recycled precipitation. This shows an advantage of applying an isotope approach in conjunction with conventional evaporation estimates to quantify both moisture recycling and net losses by evaporation. Depth profiles of 18O/16O and 2H/1H in the Great Lakes show a lack of isotopic stratification in summer months despite an established thermocline. These results are indicative of very low over-lake evaporation during warm summer months, with the bulk of evaporation occurring during the fall and winter. This seasonality in evaporation losses is supported by energy balance studies. For Lakes Michigan and Huron, the isotope mass balance approach provides a new perspective into water exchange and evaporation from these lakes. This isotope investigation shows that Lake Michigan and Lake Huron waters are distinct, despite sharing a common lake level. This finding advocates for the separate consideration of Lake Michigan and Lake Huron in future hydrologic studies

Revisiting the contribution of transpiration to global terrestrial evapotranspiration

Even though knowing the contributions of transpiration (T), soil and open water evaporation (E), and interception (I) to terrestrial evapotranspiration (ET=T+E+I) is crucial for understanding the hydrological cycle and its connection to ecological processes, the fraction of T is unattainable by traditional measurement techniques over large scales. Previously reported global mean T/(E+T+I) from multiple independent sources, including satellite-based estimations, reanalysis, land surface models, and isotopic measurements, varies substantially from 24% to 90%. Here we develop a new ET partitioning algorithm, which combines global evapotranspiration estimates and relationships between leaf area index (LAI) and T/(E+T) for different vegetation types, to upscale a wide range of published site-scale measurements. We show that transpiration accounts for about 57.2% (with standard deviation6.8%) of global terrestrial ET. Our approach bridges the scale gap between site measurements and global model simulations,and can be simply implemented into current global climate models to improve biological CO2 flux simulations

Predicting Spatial Patterns in Precipitation Isotope (δ2H and δ18O) Seasonality Using Sinusoidal Isoscapes

Understanding how precipitation isotopes vary spatially and temporally is important for tracer applications. We tested how well month‐to‐month variations in precipitation δ18O and δ2H were captured by sinusoidal cycles, and how well spatial variations in these seasonal cycles could be predicted, across Switzerland. Sine functions representing seasonal cycles in precipitation isotopes explained between 47% and 94% of the variance in monthly δ18O and δ2H values at each monitoring site. A significant sinusoidal cycle was also observed in line‐conditioned excess. We interpolated the amplitudes, phases, and offsets of these sine functions across the landscape, using multiple linear regression models based on site characteristics. These interpolated maps, here referred to as a sinusoidal isoscape, reproduced monthly observations with prediction errors that were smaller than or similar to those of other isoscapes. Sinusoidal isoscapes are likely broadly useful because they concisely describe seasonal isotopic behavior and can be estimated efficiently from sparse or irregular data. Plain Language Summary Naturally occurring isotopic variations in precipitation are used to trace water movement through landscapes and ecosystems. However, direct measurements are often unavailable, so many isotope‐based approaches to studying terrestrial processes require predicted isotopic inputs. We found that the isotopic composition of precipitation follows a predictable seasonal pattern. We developed a new approach for mapping precipitation isotope seasonality that will be useful in a wide range of fields

Global Sinusoidal Seasonality in Precipitation Isotopes

Quantifying seasonal variations in precipitation δ2H and δ18O is important for many stable isotope applications, including inferring plant water sources and streamflow ages. Our objective is to develop a data product that concisely quantifies the seasonality of stable isotope ratios in precipitation. We fit sine curves defined by amplitude, phase, and offset parameters to quantify annual precipitation isotope cycles at 653 meteorological stations on all seven continents. At most of these stations, including in tropical and subtropical regions, sine curves can represent the seasonal cycles in precipitation isotopes. Additionally, the amplitude, phase, and offset parameters of these sine curves correlate with site climatic and geographic characteristics. Multiple linear regression models based on these site characteristics capture most of the global variation in precipitation isotope amplitudes and offsets; while phase values were not well predicted by regression models globally, they were captured by zonal (0–30∘ and 30–90∘) regressions, which were then used to produce global maps. These global maps of sinusoidal seasonality in precipitation isotopes based on regression models were adjusted for the residual spatial variations that were not captured by the regression models. The resulting mean prediction errors were 0.49 ‰ for δ18O amplitude, 0.73 ‰ for δ18O offset (and 4.0 ‰ and 7.4 ‰ for δ2H amplitude and offset), 8 d for phase values at latitudes outside of 30∘, and 20 d for phase values at latitudes inside of 30∘. We make the gridded global maps of precipitation δ2H and δ18O seasonality publicly available. We also make tabulated site data and fitted sine curve parameters available to support the development of regionally calibrated models, which will often be more accurate than our global model for regionally specific studies

Late-Glacial to Late-holocene Shifts in Global Precipitation Delta(sup 18)O

Reconstructions of Quaternary climate are often based on the isotopic content of paleo-precipitation preserved in proxy records. While many paleo-precipitation isotope records are available, few studies have synthesized these dispersed records to explore spatial patterns of late-glacial precipitation delta(sup 18)O. Here we present a synthesis of 86 globally distributed groundwater (n 59), cave calcite (n 15) and ice core (n 12) isotope records spanning the late-glacial (defined as 50,000 to 20,000 years ago) to the late-Holocene (within the past 5000 years). We show that precipitation delta(sup 18)O changes from the late-glacial to the late-Holocene range from -7.1% (delta(sup 18)O(late-Holocene) > delta(sup 18)O(late-glacial) to +1.7% (delta(sup 18)O(late-glacial) > delta(sup 18)O(late-Holocene), with the majority (77) of records having lower late-glacial delta(sup 18)O than late-Holocene delta(sup 18)O values. High-magnitude, negative precipitation delta(sup 18)O shifts are common at high latitudes, high altitudes and continental interiors

Modeling the isotopic evolution of snowpack and snowmelt : Testing a spatially distributed parsimonious approach

This work was funded by the NERC/JPI SIWA project (NE/M019896/1) and the European Research Council ERC (project GA 335910 VeWa). The Krycklan part of this study was supported by grants from the Knut and Alice Wallenberg Foundation (Branch-points), Swedish Research Council (SITES), SKB and Kempe foundation. The data and model code is available upon request. Authors declare that they have no conflict of interest. We would like to thank the three anonymous reviewers for their constructive comments that improved the manuscript.Peer reviewedPublisher PD

Characterizing the diurnal patterns of errors in the prediction of evapotranspiration by several land-surface models : an NACP analysis

Land-surface models use different formulations of stomatal conductance and plant hydraulics, and it is unclear which type of model best matches the observed surface-atmosphere water flux. We use the North American Carbon Program data set of latent heat flux (LE) measurements from 25 sites and predictions from 9 models to evaluate models' ability to resolve subdaily dynamics of transpiration. Despite overall good forecast at the seasonal scale, the models have difficulty resolving the dynamics of intradaily hysteresis. The majority of models tend to underestimate LE in the prenoon hours and overestimate in the evening. We hypothesize that this is a result of unresolved afternoon stomatal closure due to hydrodynamic stresses. Although no model or stomata parameterization was consistently best or worst in terms of ability to predict LE, errors in model-simulated LE were consistently largest and most variable when soil moisture was moderate and vapor pressure deficit was moderate to limiting. Nearly all models demonstrate a tendency to underestimate the degree of maximum hysteresis which, across all sites studied, is most pronounced during moisture-limited conditions. These diurnal error patterns are consistent with models' diminished ability to accurately simulate the natural hysteresis of transpiration. We propose that the lack of representation of plant hydrodynamics is, in part, responsible for these error patterns

Illuminating hydrological processes at the soil-vegetation-atmosphere interface with water stable isotopes

Funded by DFG research project “From Catchments as Organised Systems to Models based on Functional Units” (FOR 1Peer reviewedPublisher PDFPublisher PD