Local groundwater decline exacerbates response of dryland riparian woodlands to climatic drought.

Dryland riparian woodlands are considered to be locally buffered from droughts by shallow and stable groundwater levels. However, climate change is causing more frequent and severe drought events, accompanied by warmer temperatures, collectively threatening the persistence of these groundwater dependent ecosystems through a combination of increasing evaporative demand and decreasing groundwater supply. We conducted a dendro-isotopic analysis of radial growth and seasonal (semi-annual) carbon isotope discrimination (Δ13 C) to investigate the response of riparian cottonwood stands to the unprecedented California-wide drought from 2012 to 2019, along the largest remaining free-flowing river in Southern California. Our goals were to identify principal drivers and indicators of drought stress for dryland riparian woodlands, determine their thresholds of tolerance to hydroclimatic stressors, and ultimately assess their vulnerability to climate change. Riparian trees were highly responsive to drought conditions along the river, exhibiting suppressed growth and strong stomatal closure (inferred from reduced Δ13 C) during peak drought years. However, patterns of radial growth and Δ13 C were quite variable among sites that differed in climatic conditions and rate of groundwater decline. We show that the rate of groundwater decline, as opposed to climate factors, was the primary driver of site differences in drought stress, and trees showed greater sensitivity to temperature at sites subjected to faster groundwater decline. Across sites, higher correlation between radial growth and Δ13 C for individual trees, and higher inter-correlation of Δ13 C among trees were indicative of greater drought stress. Trees showed a threshold of tolerance to groundwater decline at 0.5 m year-1 beyond which drought stress became increasingly evident and severe. For sites that exceeded this threshold, peak physiological stress occurred when total groundwater recession exceeded ~3 m. These findings indicate that drought-induced groundwater decline associated with more extreme droughts is a primary threat to dryland riparian woodlands and increases their susceptibility to projected warmer temperatures.

Hourly potential evapotranspiration at 0.1° resolution for the global land surface from 1981-present.

Challenges exist for assessing the impacts of climate and climate change on the hydrological cycle on local and regional scales, and in turn on water resources, food, energy, and natural hazards. Potential evapotranspiration (PET) represents atmospheric demand for water, which is required at high spatial and temporal resolutions to compute actual evapotranspiration and thus close the water balance near the land surface for many such applications, but there are currently no available high-resolution datasets of PET. Here we develop an hourly PET dataset (hPET) for the global land surface at 0.1° spatial resolution, based on output from the recently developed ERA5-Land reanalysis dataset, over the period 1981 to present. We show how hPET compares to other available global PET datasets, over common spatiotemporal resolutions and time frames, with respect to spatial patterns of climatology and seasonal variations for selected humid and arid locations across the globe. We provide the data for users to employ for multiple applications to explore diurnal and seasonal variations in evaporative demand for water.

Groundwater dependence of riparian woodlands and the disrupting effect of anthropogenically altered streamflow.

Riparian ecosystems fundamentally depend on groundwater, especially in dryland regions, yet their water requirements and sources are rarely considered in water resource management decisions. Until recently, technological limitations and data gaps have hindered assessment of groundwater influences on riparian ecosystem health at the spatial and temporal scales relevant to policy and management. Here, we analyze Sentinel-2-derived normalized difference vegetation index (NDVI; n = 5,335,472 observations), field-based groundwater elevation (n = 32,051 observations), and streamflow alteration data for riparian woodland communities (n = 22,153 polygons) over a 5-y period (2015 to 2020) across California. We find that riparian woodlands exhibit a stress response to deeper groundwater, as evidenced by concurrent declines in greenness represented by NDVI. Furthermore, we find greater seasonal coupling of canopy greenness to groundwater for vegetation along streams with natural flow regimes in comparison with anthropogenically altered streams, particularly in the most water-limited regions. These patterns suggest that many riparian woodlands in California are subsidized by water management practices. Riparian woodland communities rely on naturally variable groundwater and streamflow components to sustain key ecological processes, such as recruitment and succession. Altered flow regimes, which stabilize streamflow throughout the year and artificially enhance water supplies to riparian vegetation in the dry season, disrupt the seasonal cycles of abiotic drivers to which these Mediterranean forests are adapted. Consequently, our analysis suggests that many riparian ecosystems have become reliant on anthropogenically altered flow regimes, making them more vulnerable and less resilient to rapid hydrologic change, potentially leading to future riparian forest loss across increasingly stressed dryland regions.

Exploring groundwater and soil water storage changes across the CONUS at 12.5 km resolution by a Bayesian integration of GRACE data into W3RA.

Climate variability and change along with anthropogenic water use have affected the (re)distribution of water storage and fluxes across the Contiguous United States (CONUS). Available hydrological models, however, do not represent recent changes in the water cycle. Therefore, in this study, a novel Bayesian Markov Chain Monte Carlo-based Data Assimilation (MCMC-DA) approach is formulated to integrate Terrestrial Water Storage changes (TWSC) from the Gravity Recovery and Climate Experiment (GRACE) satellite mission into the W3RA water balance model. The benefit of this integration is its dynamic solution that uses GRACE TWSC to update W3RA's individual water storage estimates while rigorously accounting for uncertainties. It also down-scales GRACE data and provides groundwater and soil water storage changes at ~12.5 km resolution across the CONUS covering 2003-2017. Independent validations are performed against in-situ groundwater data (from USGS) and Climate Change Initiative (CCI) soil moisture products from the European Space Agency (ESA). Our results indicate that MCMC-DA introduces trends, which exist in GRACE TWSC, mostly to the groundwater storage and to a lesser extent to the soil water storage. Higher similarity is found between groundwater estimation of MCMC-DA and those of USGS in the southeastern CONUS. We also show a stronger linear trend in MCMC-DA soil water storage across the CONUS, compared to W3RA (changing from ±0.5 mm/yr to ±2 mm/yr), which is closer to independent estimates from the ESA CCI. MCMC-DA also improves the estimation of soil water storage in regions with high forest intensity, where ESA CCI and hydrological models have difficulties in capturing the soil-vegetation-atmosphere continuum. The representation of El Niño Southern Oscillation (ENSO)-related variability in groundwater and soil water storage are found to be considerably improved after integrating GRACE TWSC with W3RA. This new hybrid approach shows promise for understanding the links between climate and the water balance over broad regions.

Aridity is expressed in river topography globally.

It has long been suggested that climate shapes land surface topography through interactions between rainfall, runoff and erosion in drainage basins1-4. The longitudinal profile of a river (elevation versus distance downstream) is a key morphological attribute that reflects the history of drainage basin evolution, so its form should be diagnostic of the regional expression of climate and its interaction with the land surface5-9. However, both detecting climatic signatures in longitudinal profiles and deciphering the climatic mechanisms of their development have been challenging, owing to the lack of relevant global data and to the variable effects of tectonics, lithology, land surface properties and human activities10,11. Here we present a global dataset of 333,502 river longitudinal profiles, and use it to explore differences in overall profile shape (concavity) across climate zones. We show that river profiles are systematically straighter with increasing aridity. Through simple numerical modelling, we demonstrate that these global patterns in longitudinal profile shape can be explained by hydrological controls that reflect rainfall-runoff regimes in different climate zones. The most important of these is the downstream rate of change in streamflow, independent of the area of the drainage basin. Our results illustrate that river topography expresses a signature of aridity, suggesting that climate is a first-order control on the evolution of the drainage basin.

Isotopic evidence for mercury photoreduction and retention on particles in surface waters of Central California, USA.

Cache Creek (Coast Range, California) and the Yuba River (Sierra Nevada Foothills, California) are two river systems affected by extensive mercury (Hg) contamination due to legacy sources of Hg related to mining. Stable Hg isotope techniques have proven useful for elucidating the complex cycling of Hg within aquatic ecosystems, and we applied these techniques to improve understanding of Hg and methylmercury (MeHg) transformations in these watersheds. Total mercury (THg) concentrations and Hg stable isotope ratios were measured in filtered surface waters and suspended particulate matter collected from 14 sites within the Cache Cr. and Yuba R. watersheds. Filtered surface waters from both watersheds exhibited values of ∆199Hg (0.37‰ to 0.71‰), consistently elevated above those observed in sediments (∆199Hg average = 0.07‰). Associated suspended particulates from these surface water samples displayed a much greater range of values for ∆199Hg (-0.61‰ to 0.70‰), although suspended particulates from the Yuba R. exhibited mostly negative ∆199Hg values (-0.61‰ to 0.10‰). The relationship between ∆199Hg and ∆201Hg in the filtered surface waters and associated suspended particulates was calculated using a bivariate York regression, yielding a slope of 1.57 ±â€¯0.49 (±2SE) for the Yuba R. and 1.40 ±â€¯0.27 (±2SE) for Cache Cr., both within error of the previously reported experimentally-derived slopes for MeHg- and inorganic Hg(II)-photoreduction. This provides isotopic evidence that Hg photoreduction is occurring within these surface waters to a significant degree, and suspended particulate phases are retaining the reduced product of Hg photoreduction, particularly within the Yuba R. The isotopic compositions of filtered surface waters are consistent with the isotopic signatures recorded in biota at low trophic positions within these watersheds, suggesting that the reservoir of Hg incorporated within the biota of these systems is similar to the filter-passing Hg fraction in surface waters.

Methylmercury degradation and exposure pathways in streams and wetlands impacted by historical mining.

Monomethyl mercury (MMHg) and total mercury (THg) concentrations and Hg stable isotope ratios (δ(202)Hg and Δ(199)Hg) were measured in sediment and aquatic organisms from Cache Creek (California Coast Range) and Yolo Bypass (Sacramento Valley). Cache Creek sediment had a large range in THg (87 to 3870ng/g) and δ(202)Hg (-1.69 to -0.20‰) reflecting the heterogeneity of Hg mining sources in sediment. The δ(202)Hg of Yolo Bypass wetland sediment suggests a mixture of high and low THg sediment sources. Relationships between %MMHg (the percent ratio of MMHg to THg) and Hg isotope values (δ(202)Hg and Δ(199)Hg) in fish and macroinvertebrates were used to identify and estimate the isotopic composition of MMHg. Deviation from linear relationships was found between %MMHg and Hg isotope values, which is indicative of the bioaccumulation of isotopically distinct pools of MMHg. The isotopic composition of pre-photodegraded MMHg (i.e., subtracting fractionation from photochemical reactions) was estimated and contrasting relationships were observed between the estimated δ(202)Hg of pre-photodegraded MMHg and sediment IHg. Cache Creek had mass dependent fractionation (MDF; δ(202)Hg) of at least -0.4‰ whereas Yolo Bypass had MDF of +0.2 to +0.5‰. This result supports the hypothesis that Hg isotope fractionation between IHg and MMHg observed in rivers (-MDF) is unique compared to +MDF observed in non-flowing water environments such as wetlands, lakes, and the coastal ocean.

Hydrologic indicators of hot spots and hot moments of mercury methylation potential along river corridors.

The biogeochemical cycling of metals and other contaminants in river-floodplain corridors is controlled by microbial activity responding to dynamic redox conditions. Riverine flooding thus has the potential to affect speciation of redox-sensitive metals such as mercury (Hg). Therefore, inundation history over a period of decades potentially holds information on past production of bioavailable Hg. We investigate this within a Northern California river system with a legacy of landscape-scale 19th century hydraulic gold mining. We combine hydraulic modeling, Hg measurements in sediment and biota, and first-order calculations of mercury transformation to assess the potential role of river floodplains in producing monomethylmercury (MMHg), a neurotoxin which accumulates in local and migratory food webs. We identify frequently inundated floodplain areas, as well as floodplain areas inundated for long periods. We quantify the probability of MMHg production potential (MPP) associated with hydrology in each sector of the river system as a function of the spatial patterns of overbank inundation and drainage, which affect long-term redox history of contaminated sediments. Our findings identify river floodplains as periodic, temporary, yet potentially important, loci of biogeochemical transformation in which contaminants may undergo change during limited periods of the hydrologic record. We suggest that inundation is an important driver of MPP in river corridors and that the entire flow history must be analyzed retrospectively in terms of inundation magnitude and frequency in order to accurately assess biogeochemical risks, rather than merely highlighting the largest floods or low-flow periods. MMHg bioaccumulation within the aquatic food web in this system may pose a major risk to humans and waterfowl that eat migratory salmonids, which are being encouraged to come up these rivers to spawn. There is a long-term pattern of MPP under the current flow regime that is likely to be accentuated by increasingly common large floods with extended duration.

Isotopic Composition of Inorganic Mercury and Methylmercury Downstream of a Historical Gold Mining Region.

We measured total mercury (THg) and monomethyl mercury (MMHg) concentrations and mercury (Hg) isotopic compositions in sediment and aquatic organisms from the Yuba River (California, USA) to identify Hg sources and biogeochemical transformations downstream of a historical gold mining region. Sediment THg concentrations and δ(202)Hg decreased from the upper Yuba Fan to the lower Yuba Fan and the Feather River. These results are consistent with the release of Hg during gold mining followed by downstream mixing and dilution. The Hg isotopic composition of Yuba Fan sediment (δ(202)Hg = -0.38 ± 0.17‰ and Δ(199)Hg = 0.04 ± 0.03‰; mean ± 1 SD, n = 7) provides a fingerprint of inorganic Hg (IHg) that could be methylated locally or after transport downstream. The isotopic composition of MMHg in the Yuba River food web was estimated using biota with a range of %MMHg (the percent of THg present as MMHg) and compared to IHg in sediment, algae, and the food web. The estimated δ(202)Hg of MMHg prior to photodegradation (-1.29 to -1.07‰) was lower than that of IHg and we suggest this is due to mass-dependent fractionation (MDF) of up to -0.9‰ between IHg and MMHg. This result is in contrast to net positive MDF (+0.4 to +0.8‰) previously observed in lakes, estuaries, coastal oceans, and forests. We hypothesize that this unique relationship could be due to differences in the extent or pathway of biotic MMHg degradation in stream environments.

Floodplain ecohydrology: Climatic, anthropogenic, and local physical controls on partitioning of water sources to riparian trees.

Seasonal and annual partitioning of water within river floodplains has important implications for ecohydrologic links between the water cycle and tree growth. Climatic and hydrologic shifts alter water distribution between floodplain storage reservoirs (e.g., vadose, phreatic), affecting water availability to tree roots. Water partitioning is also dependent on the physical conditions that control tree rooting depth (e.g., gravel layers that impede root growth), the sources of contributing water, the rate of water drainage, and water residence times within particular storage reservoirs. We employ instrumental climate records alongside oxygen isotopes within tree rings and regional source waters, as well as topographic data and soil depth measurements, to infer the water sources used over several decades by two co-occurring tree species within a riparian floodplain along the Rhône River in France. We find that water partitioning to riparian trees is influenced by annual (wet versus dry years) and seasonal (spring snowmelt versus spring rainfall) fluctuations in climate. This influence depends strongly on local (tree level) conditions including floodplain surface elevation and subsurface gravel layer elevation. The latter represents the upper limit of the phreatic zone and therefore controls access to shallow groundwater. The difference between them, the thickness of the vadose zone, controls total soil moisture retention capacity. These factors thus modulate the climatic influence on tree ring isotopes. Additionally, we identified growth signatures and tree ring isotope changes associated with recent restoration of minimum streamflows in the Rhône, which made new phreatic water sources available to some trees in otherwise dry years. KEY POINTS: Water shifts due to climatic fluctuations between floodplain storage reservoirsAnthropogenic changes to hydrology directly impact water available to treesEcohydrologic approaches to integration of hydrology afford new possibilities.

Enduring legacy of a toxic fan via episodic redistribution of California gold mining debris.

The interrelationships between hydrologically driven evolution of legacy landscapes downstream of major mining districts and the contamination of lowland ecosystems are poorly understood over centennial time scales. Here, we demonstrate within piedmont valleys of California's Sierra Nevada, through new and historical data supported by modeling, that anthropogenic fans produced by 19th century gold mining comprise an episodically persistent source of sediment-adsorbed Hg to lowlands. Within the enormous, iconic Yuba Fan, we highlight (i) an apparent shift in the relative processes of fan evolution from gradual vertical channel entrenchment to punctuated lateral erosion of fan terraces, thus enabling entrainment of large volumes of Hg-laden sediment during individual floods, and (ii) systematic intrafan redistribution and downstream progradation of fan sediment into the Central Valley, triggered by terrace erosion during increasingly long, 10-y flood events. Each major flood apparently erodes stored sediment and delivers to sensitive lowlands the equivalent of ~10-30% of the entire postmining Sierran Hg mass so far conveyed to the San Francisco Bay-Delta (SFBD). This process of protracted but episodic erosion of legacy sediment and associated Hg is likely to persist for >10(4) y. It creates, within an immense swath of river corridor well upstream of the SFBD, new contaminated floodplain surfaces primed for Hg methylation and augments/replenishes potential Hg sources to the SFBD. Anticipation, prediction, and management of toxic sediment delivery, and corresponding risks to lowland ecology and human society globally, depend on the morphodynamic stage of anthropogenic fan evolution, synergistically coupled to changing frequency of and duration extreme floods.

Sediment-adsorbed total mercury flux through Yolo Bypass, the primary floodway and wetland in the Sacramento Valley, California.

The fate and transport of mercury are of critical concern in lowland floodplains and wetlands worldwide, especially those with a history of upstream mining that increases the mobility of both dissolved and sediment-bound Hg in watersheds. A mass budget of total mercury (THg) quantifies sources and storage for particular areas - knowledge that is required for understanding of management options in lowland floodplains. In order to assess contaminant risk in the largest flood-control bypass, prime wetland, and restoration target in the Sacramento River basin, we estimated empirical relationships between THg, suspended sediment concentration (SSC), and streamflow (Q) for each of the major inputs and outputs using data from various publicly available sources. These relationships were improved by incorporating statistical representations of the dynamics of seasonal and intra-flood exhaustion (hysteresis) of sediment and mercury. Using continuous records of Q to estimate SSC suspended sediment flux and SSC to estimate THg flux, we computed the net transfer of sediment-adsorbed mercury through the Yolo Bypass over a decade, 1993-2003. Flood control weirs spilling Sacramento River floodwaters into the bypass deliver ~75% of the water and ~50% of the river's suspended sediment load, while one Coast Range tributary of the bypass, Cache Creek, contributes twice the THg load of the mainstem Sacramento. Although estimated sediment flux entering Yolo Bypass is balanced by efflux to the Sacramento/San Francisco Bay-Delta, there is much evidence of deposition and remobilization of sediment in Yolo Bypass during flooding. These factors point to the importance of the bypass as sedimentary reservoir and as an evolving substrate for biogeochemical processing of heavy metals. The estimates of mercury flux suggest net deposition of ~500 kg in the 24,000 ha floodway over a decade, dominated by two large floods, representing a storage reservoir for this important contaminant.