Ensemble of CMIP6 derived reference and potential evapotranspiration with radiative and advective components.

Assessing changes in future aridity requires an understanding of variations in the atmospheric demand for water. Such assessments are often driven by estimations of potential evapotranspiration (ETP) and/or reference evapotranspiration (ET0), yet no comprehensive and validated estimate of these climate metrics exists to date from the Coupled Model Intercomparison Project 6 (CMIP6). Here we describe the development and validation of a published dataset of global monthly estimates of the Penman-Monteith derived ET0, its advective and radiation components, Priestley-Taylor derived ETP, and vapor pressure deficit from 16 CMIP6 projections and four emissions scenarios. Historical validation of the ensemble of CMIP6 evaporative demand shows general agreement with observationally-derived baselines of ET0 and ETP from the Climate Research Unit (CRU) and ERA5-Land reanalysis products, with GCM biases driven primarily by regional differences in modeled humidity and advective contributions to ET0. Overall, evaporative demand is projected to increase across all emissions scenarios, with the largest increases over polar regions, and with a larger contribution from advection particularly for regions with higher baseline ET0.

Satellites reveal widespread decline in global lake water storage.

Climate change and human activities increasingly threaten lakes that store 87% of Earth's liquid surface fresh water. Yet, recent trends and drivers of lake volume change remain largely unknown globally. Here, we analyze the 1972 largest global lakes using three decades of satellite observations, climate data, and hydrologic models, finding statistically significant storage declines for 53% of these water bodies over the period 1992-2020. The net volume loss in natural lakes is largely attributable to climate warming, increasing evaporative demand, and human water consumption, whereas sedimentation dominates storage losses in reservoirs. We estimate that roughly one-quarter of the world's population resides in a basin of a drying lake, underscoring the necessity of incorporating climate change and sedimentation impacts into sustainable water resources management.

Bark beetle impacts on forest evapotranspiration and its partitioning.

Insect outbreaks affect forest structure and function and represent a major category of forest disturbance globally. However, the resulting impacts on evapotranspiration (ET), and especially hydrological partitioning between the abiotic (evaporation) and biotic (transpiration) components of total ET, are not well constrained. As a result, we combined remote sensing, eddy covariance, and hydrological modeling approaches to determine the effects of bark beetle outbreak on ET and its partitioning at multiple scales throughout the Southern Rocky Mountain Ecoregion (SRME), USA. At the eddy covariance measurement scale, 85 % of the forest was affected by beetles, and water year ET as a fraction of precipitation (P) decreased by 30 % relative to a control site, with 31 % greater reductions in growing season transpiration relative to total ET. At the ecoregion scale, satellite remote sensing masked to areas of >80 % tree mortality showed corresponding ET/P reductions of 9-15 % that occurred 6-8 years post-disturbance, and indicated that the majority of the total reduction occurred during the growing season; the Variable Infiltration Capacity hydrological model showed an associated 9-18 % increase in the ecoregion runoff ratio. Long-term (16-18 year) ET and vegetation mortality datasets extend the length of previously published analyses and allowed for clear characterization of the forest recovery period. During that time, transpiration recovery outpaced total ET recovery, which was lagged in part due to persistently reduced winter sublimation, and there was associated evidence of increasing late summer vegetation moisture stress. Overall, comparison of three independent methods and two partitioning approaches demonstrated a net negative impact of bark beetles on ET, and a relatively greater negative impact on transpiration, following bark beetle outbreak in the SRME.

Recent decreases in snow water storage in western North America.

Mountain snowpacks act as natural water towers, storing winter precipitation until summer months when downstream water demand is greatest. We introduce a Snow Storage Index (SSI), representing the temporal phase difference between daily precipitation and surface water inputs-sum of rainfall and snowmelt into terrestrial systems-weighted by relative magnitudes. Different from snow water equivalent or snow fraction, the SSI represents the degree to which the snowpack delays the timing and magnitude of surface water inputs relative to precipitation, a fundamental component of how snow water storage influences the hydrologic cycle. In western North America, annual SSI has decreased (p < 0.05) from 1950-2013 in over 25% of mountainous areas, as a result of substantially earlier snowmelt and rainfall in spring months, with additional declines in winter precipitation. The SSI and associated trends offer a new perspective on hydrologic sensitivity to climate change which have broad implications for water resources and ecosystems.

Development and Validation of an In Situ Groundwater Abstraction Sensor Network, Hydrologic Statistical Model, and Blockchain Trading Platform: A Demonstration in Solano County, California.

Megadrought in the western United States is jeopardizing water security. Groundwater regulations, such as California's Sustainable Groundwater Management Act (SGMA), aim to preserve groundwater resources in overdrafted basins. Water agencies must establish sufficient monitoring systems to measure local groundwater abstraction and devise plans to moderate groundwater use. However, few technologies are available to monitor and regulate groundwater abstraction spatially and temporally. In this study, we deployed satellite-connected electrical current sensors on 11 agricultural groundwater pumps in Solano County, California over 2019-2022. A high correlation (R 2 = 0.706) was found between the in situ sensors and in-line flow meters. We then combine in situ sensor data with a land surface model to develop a multiple linear regression model of groundwater abstraction and groundwater level. Using a 10-fold cross-validation, it is found that our predictive groundwater abstraction model has approximately a 3.5% bias and a mean absolute error of 1.21 acre-feet, while our predictive groundwater level model has approximately 4.2% bias and about 5.9 acre-feet mean absolute error. Finally, we integrated these data with a blockchain-based groundwater credit trading platform to demonstrate how such a tool could be used for SGMA compliance.

A review of simulation experiment techniques used to analyze wildfire effects on water quality and supply.

This review addresses the critical knowledge gap of techniques simulating combustion and heating characteristics present in natural wildfires and their use in assessing postfire impacts on water quality and quantity. Our assessment includes both laboratory and plot-scale techniques with burn and rainfall simulation components. Studies included focus on advancing understanding of changes in chemical and physical properties of soil, as well as subsequent runoff changes. Advantages of simulation experiments include: overcoming logistical challenges of collecting in situ wildfire data, reducing the high spatial variability observed in natural settings (i.e., the heterogeneity of burn intensity and the underlying vegetation and soil properties), and controlling the magnitude of key drivers of wildfire impacts. In sum, simulation experiments allow for more direct attribution of water quality and quantity responses to specific drivers than experiments conducted in situ. Drawbacks of simulation techniques include the limitation of observing only local-scale processes, the potential misrepresentation of natural settings (i.e., lack of spatial variability in vegetation, soil structure, burn intensity, etc.), uncertainty introduced through experimental error, and subsequent challenges in upscaling results to larger scales relevant for water management. This review focuses primarily on simulation techniques, with the goal of providing a foundation of knowledge for the design of future simulation experiments.

Growing impact of wildfire on western US water supply.

Streamflow often increases after fire, but the persistence of this effect and its importance to present and future regional water resources are unclear. This paper addresses these knowledge gaps for the western United States (WUS), where annual forest fire area increased by more than 1,100% during 1984 to 2020. Among 72 forested basins across the WUS that burned between 1984 and 2019, the multibasin mean streamflow was significantly elevated by 0.19 SDs (P < 0.01) for an average of 6 water years postfire, compared to the range of results expected from climate alone. Significance is assessed by comparing prefire and postfire streamflow responses to climate and also to streamflow among 107 control basins that experienced little to no wildfire during the study period. The streamflow response scales with fire extent: among the 29 basins where >20% of forest area burned in a year, streamflow over the first 6 water years postfire increased by a multibasin average of 0.38 SDs, or 30%. Postfire streamflow increases were significant in all four seasons. Historical fire-climate relationships combined with climate model projections suggest that 2021 to 2050 will see repeated years when climate is more fire-conducive than in 2020, the year currently holding the modern record for WUS forest area burned. These findings center on relatively small, minimally managed basins, but our results suggest that burned areas will grow enough over the next 3 decades to enhance streamflow at regional scales. Wildfire is an emerging driver of runoff change that will increasingly alter climate impacts on water supplies and runoff-related risks.

A multi-sensor evaluation of precipitation uncertainty for landslide-triggering storm events.

Extreme precipitation can have profound consequences for communities, resulting in natural hazards such as rainfall-triggered landslides that cause casualties and extensive property damage. A key challenge to understanding and predicting rainfall-triggered landslides comes from observational uncertainties in the depth and intensity of precipitation preceding the event. Practitioners and researchers must select from a wide range of precipitation products, often with little guidance. Here we evaluate the degree of precipitation uncertainty across multiple precipitation products for a large set of landslide-triggering storm events and investigate the impact of these uncertainties on predicted landslide probability using published intensity-duration thresholds. The average intensity, peak intensity, duration, and NOAA-Atlas return periods are compared ahead of 177 reported landslides across the continental United States and Canada. Precipitation data are taken from four products that cover disparate measurement methods: near real-time and post-processed satellite (IMERG), radar (MRMS), and gauge-based (NLDAS-2). Landslide-triggering precipitation was found to vary widely across precipitation products with the depth of individual storm events diverging by as much as 296 mm with an average range of 51 mm. Peak intensity measurements, which are typically influential in triggering landslides, were also highly variable with an average range of 7.8 mm/h and as much as 57 mm/h. The two products more reliant upon ground-based observations (MRMS and NLDAS-2) performed better at identifying landslides according to published intensity-duration storm thresholds, but all products exhibited hit ratios of greater than 0.56. A greater proportion of landslides were predicted when including only manually verified landslide locations. We recommend practitioners consider low-latency products like MRMS for investigating landslides, given their near-real time data availability and good performance in detecting landslides. Practitioners would be well-served considering more than one product as a way to confirm intense storm signals and minimize the influence of noise and false alarms.

A continental-scale soil evaporation dataset derived from Soil Moisture Active Passive satellite drying rates.

This manuscript describes an observationally-based dataset of soil evaporation for the conterminous U.S. (CONUS), gridded to a 9 km resolution for the time-period of April 2015-March 2019. This product is termed E-SMAP (Evaporation-Soil Moisture Active Passive) in which soil evaporation is estimated from the surface layer, defined by the SMAP sensing depth of 50 mm, between SMAP overpass intervals that are screened on the basis of precipitation and SMAP quality control flags. Soil evaporation is estimated using a water balance of the surface soil that we show is largely dominated by SMAP-observed soil drying. E-SMAP soil evaporation is on average 0.72 mm day-1, which falls within the range of soil evaporation estimates (0.17-0.89 mm day-1) derived from operational land surface models and an alternative remote sensing product. E-SMAP is independent from existing soil evaporation estimates and therefore has the potential to improve understanding of evapotranspiration partitioning and model development.

PEMIP: Post-fire erosion model inter-comparison project.

Land managers often need to predict watershed-scale erosion rates after disturbance or other land cover changes. This study compared commonly used hillslope erosion models to simulate post-fire sediment yields (SY) at both hillslope and watershed scales within the High Park Fire, Colorado, U.S.A. At hillslope scale, simulated SY from four models- RUSLE, AGWA/KINEROS2, WEPP, and a site-specific regression model-were compared to observed SY at 29 hillslopes. At the watershed scale, RUSLE, AGWA/KINEROS2, and WEPP were applied to simulate spatial patterns of SY for two 14-16 km2 watersheds using different scales (0.5-25 ha) of hillslope discretization. Simulated spatial patterns were compared between models and to densities of channel heads across the watersheds. Three additional erosion algorithms were implemented within a land surface model to evaluate effects of parameter uncertainty. At the hillslope scale, SY was only significantly correlated to observed SY for the empirical model, but at the watershed scale, sediment loads were significantly correlated to observed channel head densities for all models. Watershed sediment load increased with the size of the hillslope sub-units due to the nonlinear effects of hillslope length on simulated erosion. SY's were closest in magnitude to expected watershed-scale SY when models were divided into the smallest hillslopes. These findings demonstrate that current erosion models are fairly consistent at identifying areas with low and high erosion potential, but the wide range of predicted SY and poor fit to observed SY highlight the need for better field observations and model calibration to obtain more accurate simulations.

Large contribution from anthropogenic warming to an emerging North American megadrought.

Severe and persistent 21st-century drought in southwestern North America (SWNA) motivates comparisons to medieval megadroughts and questions about the role of anthropogenic climate change. We use hydrological modeling and new 1200-year tree-ring reconstructions of summer soil moisture to demonstrate that the 2000-2018 SWNA drought was the second driest 19-year period since 800 CE, exceeded only by a late-1500s megadrought. The megadrought-like trajectory of 2000-2018 soil moisture was driven by natural variability superimposed on drying due to anthropogenic warming. Anthropogenic trends in temperature, relative humidity, and precipitation estimated from 31 climate models account for 47% (model interquartiles of 35 to 105%) of the 2000-2018 drought severity, pushing an otherwise moderate drought onto a trajectory comparable to the worst SWNA megadroughts since 800 CE.

Spatial variation of the rain-snow temperature threshold across the Northern Hemisphere.

Despite the importance of precipitation phase to global hydroclimate simulations, many land surface models use spatially uniform air temperature thresholds to partition rain and snow. Here we show, through the analysis of a 29-year observational dataset (n = 17.8 million), that the air temperature at which rain and snow fall in equal frequency varies significantly across the Northern Hemisphere, averaging 1.0 °C and ranging from -0.4 to 2.4 °C for 95% of the stations. Continental climates generally exhibit the warmest rain-snow thresholds and maritime the coolest. Simulations show precipitation phase methods incorporating humidity perform better than air temperature-only methods, particularly at relative humidity values below saturation and air temperatures between 0.6 and 3.4 °C. We also present the first continuous Northern Hemisphere map of rain-snow thresholds, underlining the spatial variability of precipitation phase partitioning. These results suggest precipitation phase could be better predicted using humidity and air temperature in large-scale land surface model runs.

Development of a gridded meteorological dataset over Java island, Indonesia 1985-2014.

We describe a gridded daily meteorology dataset consisting of precipitation, minimum and maximum temperature over Java Island, Indonesia at 0.125°×0.125° (~14 km) resolution spanning 30 years from 1985-2014. Importantly, this data set represents a marked improvement from existing gridded data sets over Java with higher spatial resolution, derived exclusively from ground-based observations unlike existing satellite or reanalysis-based products. Gap-infilling and gridding were performed via the Inverse Distance Weighting (IDW) interpolation method (radius, r, of 25 km and power of influence, α, of 3 as optimal parameters) restricted to only those stations including at least 3,650 days (~10 years) of valid data. We employed MSWEP and CHIRPS rainfall products in the cross-validation. It shows that the gridded rainfall presented here produces the most reasonable performance. Visual inspection reveals an increasing performance of gridded precipitation from grid, watershed to island scale. The data set, stored in a network common data form (NetCDF), is intended to support watershed-scale and island-scale studies of short-term and long-term climate, hydrology and ecology.

A spatially comprehensive, hydrometeorological data set for Mexico, the U.S., and Southern Canada 1950-2013.

A data set of observed daily precipitation, maximum and minimum temperature, gridded to a 1/16° (~6 km) resolution, is described that spans the entire country of Mexico, the conterminous U.S. (CONUS), and regions of Canada south of 53° N for the period 1950-2013. The dataset improves previous products in spatial extent, orographic precipitation adjustment over Mexico and parts of Canada, and reduction of transboundary discontinuities. The impacts of adjusting gridded precipitation for orographic effects are quantified by scaling precipitation to an elevation-aware 1981-2010 precipitation climatology in Mexico and Canada. Differences are evaluated in terms of total precipitation as well as by hydrologic quantities simulated with a land surface model. Overall, orographic correction impacts total precipitation by up to 50% in mountainous regions outside CONUS. Hydrologic fluxes show sensitivities of similar magnitude, with discharge more sensitive than evapotranspiration and soil moisture. Because of the consistent gridding methodology, the current product reduces transboundary discontinuities as compared with a commonly used reanalysis product, making it suitable for estimating large-scale hydrometeorologic phenomena.