A Statistical Graphical Model of the California Reservoir System

The recent California drought has highlighted the potential vulnerability of the state's water management infrastructure to multiyear dry intervals. Due to the high complexity of the network, dynamic storage changes in California reservoirs on a state-wide scale have previously been difficult to model using either traditional statistical or physical approaches. Indeed, although there is a significant line of research on exploring models for single (or a small number of) reservoirs, these approaches are not amenable to a system-wide modeling of the California reservoir network due to the spatial and hydrological heterogeneities of the system. In this work, we develop a state-wide statistical graphical model to characterize the dependencies among a collection of 55 major California reservoirs across the state; this model is defined with respect to a graph in which the nodes index reservoirs and the edges specify the relationships or dependencies between reservoirs. We obtain and validate this model in a data-driven manner based on reservoir volumes over the period 2003–2016. A key feature of our framework is a quantification of the effects of external phenomena that influence the entire reservoir network. We further characterize the degree to which physical factors (e.g., state-wide Palmer Drought Severity Index (PDSI), average temperature, snow pack) and economic factors (e.g., consumer price index, number of agricultural workers) explain these external influences. As a consequence of this analysis, we obtain a system-wide health diagnosis of the reservoir network as a function of PDSI

The Sensitivity of US Wildfire Occurrence to Pre-Season Soil Moisture Conditions Across Ecosystems

It is generally accepted that year-to-year variability in moisture conditions and drought are linked with increased wildfire occurrence. However, quantifying the sensitivity of wildfire to surface moisture state at seasonal lead-times has been challenging due to the absence of a long soil moisture record with the appropriate coverage and spatial resolution for continental-scale analysis. Here we apply model simulations of surface soil moisture that numerically assimilate observations from NASA's Gravity Recovery and Climate Experiment (GRACE) mission with the US Forest Service"TM"s historical Fire-Occurrence Database over the contiguous United States. We quantify the relationships between pre-fire-season soil moisture and subsequent-year wildfire occurrence by land-cover type and produce annual probable wildfire occurrence and burned area maps at 0.25-degree resolution. Cross-validated results generally indicate a higher occurrence of smaller fires when months preceding fire season are wet, while larger fires are more frequent when soils are dry. This result is consistent with the concept of increased fuel accumulation under wet conditions in the pre-season. These results demonstrate the fundamental strength of the relationship between soil moisture and fire activity at long lead-times and are indicative of that relationship's utility for the future development of national-scale predictive capability

A Statistical Graphical Model of the California Reservoir System

The recent California drought has highlighted the potential vulnerability of the state's water management infrastructure to multiyear dry intervals. Due to the high complexity of the network, dynamic storage changes in California reservoirs on a state-wide scale have previously been difficult to model using either traditional statistical or physical approaches. Indeed, although there is a significant line of research on exploring models for single (or a small number of) reservoirs, these approaches are not amenable to a system-wide modeling of the California reservoir network due to the spatial and hydrological heterogeneities of the system. In this work, we develop a state-wide statistical graphical model to characterize the dependencies among a collection of 55 major California reservoirs across the state; this model is defined with respect to a graph in which the nodes index reservoirs and the edges specify the relationships or dependencies between reservoirs. We obtain and validate this model in a data-driven manner based on reservoir volumes over the period 2003–2016. A key feature of our framework is a quantification of the effects of external phenomena that influence the entire reservoir network. We further characterize the degree to which physical factors (e.g., state-wide Palmer Drought Severity Index (PDSI), average temperature, snow pack) and economic factors (e.g., consumer price index, number of agricultural workers) explain these external influences. As a consequence of this analysis, we obtain a system-wide health diagnosis of the reservoir network as a function of PDSI

Separating Decadal Global Water Cycle Variability from Sea Level Rise

Under a warming climate, amplification of the water cycle and changes in precipitation patterns over land are expected to occur, subsequently impacting the terrestrial water balance. On global scales, such changes in terrestrial water storage (TWS) will be reflected in the water contained in the ocean and can manifest as global sea level variations. Naturally occurring climate-driven TWS variability can temporarily obscure the long-term trend in sea level rise, in addition to modulating the impacts of sea level rise through natural periodic undulation in regional and global sea level. The internal variability of the global water cycle, therefore, confounds both the detection and attribution of sea level rise. Here, we use a suite of observations to quantify and map the contribution of TWS variability to sea level variability on decadal timescales. In particular, we find that decadal sea level variability centered in the Pacific Ocean is closely tied to low frequency variability of TWS in key areas across the globe. The unambiguous identification and clean separation of this component of variability is the missing step in uncovering the anthropogenic trend in sea level and understanding the potential for low-frequency modulation of future TWS impacts including flooding and drought

Modeling groundwater levels in California's Central Valley by hierarchical Gaussian process and neural network regression

Modeling groundwater levels continuously across California's Central Valley (CV) hydrological system is challenging due to low-quality well data which is sparsely and noisily sampled across time and space. A novel machine learning method is proposed for modeling groundwater levels by learning from a 3D lithological texture model of the CV aquifer. The proposed formulation performs multivariate regression by combining Gaussian processes (GP) and deep neural networks (DNN). Proposed hierarchical modeling approach constitutes training the DNN to learn a lithologically informed latent space where non-parametric regression with GP is performed. The methodology is applied for modeling groundwater levels across the CV during 2015 - 2020. We demonstrate the efficacy of GP-DNN regression for modeling non-stationary features in the well data with fast and reliable uncertainty quantification. Our results indicate that the 2017 and 2019 wet years in California were largely ineffective in replenishing the groundwater loss caused during previous drought years.Comment: Submitted to Water Resources Researc

High-Tide Floods and Storm Surges During Atmospheric Rivers on the US West Coast

Amospheric rivers (ARs) effect inland hydrological impacts related to extreme precipitation. However, little is known about the possible coastal hazards associated with these storms. Here we elucidate high-tide floods (HTFs) and storm surges during ARs through a statistical analysis of data from the US West Coast during 1980-2016. HTFs and landfalling ARs co-occur more often than expected from random chance. Between 10%-63% of HTFs coincide with landfalling ARs, depending on location. However, only 2%-15% of ARs coincide with HTFs, suggesting that ARs typically must co-occur with anomalously high tides or mean sea levels to cause HTFs. Storm surges during ARs are interpretable in terms of local wind, pressure, and precipitation forcing. Meridional wind and barometric pressure are the primary drivers of the storm surge. This study highlights the relevance of ARs to coastal impacts, clarifies the drivers of storm surge during ARs, and identifies future research directions

Emerging Trends in Global Freshwater Availability

Freshwater availability is changing worldwide. Here we quantify 34 trends in terrestrial water storage (TWS) observed by the Gravity Recovery and Climate Experiment (GRACE) satellites during 2002-2016 and categorize their drivers as natural interannual variability, unsustainable groundwater consumption, or climate change. Several of these trends had been lacking thorough investigation and attribution, including massive changes in northwestern China and the Okavango delta. Others are consistent with climate model predictions. This observation-based assessment of how the world's water landscape is responding to human impacts and climate variations provides a blueprint for evaluating and predicting emerging threats to water and food security

Observation-Driven Estimation of the Spatial Variability of 20th Century Sea Level Rise

Over the past two decades, sea level measurements made by satellites have given clear indications of both global and regional sea level rise. Numerous studies have sought to leverage the modern satellite record and available historic sea level data provided by tide gauges to estimate past sea level rise, leading to several estimates for the 20th century trend in global mean sea level in the range between 1 and 2 mm/yr. On regional scales, few attempts have been made to estimate trends over the same time period. This is due largely to the inhomogeneity and quality of the tide gauge network through the 20th century, which render commonly used reconstruction techniques inadequate. Here, a new approach is adopted, integrating data from a select set of tide gauges with prior estimates of spatial structure based on historical sea level forcing information from the major contributing processes over the past century. The resulting map of 20th century regional sea level rise is optimized to agree with the tide gauge-measured trends, and provides an indication of the likely contributions of different sources to regional patterns. Of equal importance, this study demonstrates the sensitivities of this regional trend map to current knowledge and uncertainty of the contributing processes

GRACE storage-runoff hystereses reveal the dynamics of regional watersheds

We characterize how regional watersheds function as simple, dynamic systems through a series of hysteresis loops using measurements from NASA's Gravity Recovery and Climate Experiment (GRACE) satellites. These loops illustrate the temporal relationship between runoff and terrestrial water storage in three regional-scale watersheds (> 150 000 km²) of the Columbia River Basin, USA and Canada. The shape and size of the hysteresis loops are controlled by the climate, topography, and geology of the watershed. The direction of the hystereses for the GRACE signals moves in opposite directions from the isolated groundwater hystereses. The subsurface water (soil moisture and groundwater) hystereses more closely resemble the storage-runoff relationship of a soil matrix. While the physical processes underlying these hystereses are inherently complex, the vertical integration of terrestrial water in the GRACE signal encapsulates the processes that govern the non-linear function of regional-scale watersheds. We use this process-based understanding to test how GRACE data can be applied prognostically to predict seasonal runoff (mean Nash-Sutcliffe Efficiency of 0.91) and monthly runoff during the low flow/high demand month of August (mean Nash-Sutcliffe Efficiency of 0.77) in all three watersheds. The global nature of GRACE data allows this same methodology to be applied in other regional-scale studies, and could be particularly useful in regions with minimal data and in trans-boundary watersheds.For a previous discussion paper please see: http://hdl.handle.net/1957/57160. This is the publisher’s final pdf. The published article is copyrighted by the author(s) and published by Copernicus Publications on behalf of the European Geosciences Union. The published article can be found at: http://www.hydrology-and-earth-system-sciences.net/. The Supplement related to this article is available onlineat doi:10.5194/hess-19-3253-2015-supplement

GRACE storage-streamflow hystereses reveal the dynamics of regional watersheds

We characterize how regional watersheds function as simple, dynamic systems through a series of hysteresis loops. These loops illustrate the temporal relationship between runoff and terrestrial water storage using measurements from NASA's Gravity Recovery and Climate Experiment (GRACE) satellites in three regional-scale watersheds (>150 000 km² ) of the Columbia River Basin, USA and Canada. The direction of the hystereses for the GRACE signal move in opposite directions from the isolated groundwater hystereses, suggesting that regional scale watersheds require soil water storage to reach a certain threshold before groundwater recharge and peak runoff occur. While the physical processes underlying these hystereses are inherently complex, the vertical integration of terrestrial water in the GRACE signal encapsulates the processes that govern the non-linear function of regional-scale watersheds. We use this process-based understanding to test how GRACE data can be applied prognostically to predict seasonal runoff (mean R² of 0.91) and monthly runoff (mean R² of 0.77) in all three watersheds. The global nature of GRACE data allows this same methodology to be applied in other regional-scale studies, and could be particularly useful in regions with minimal data and in trans-boundary watersheds.This discussion paper has been under review for the journal Hydrology and Earth System Sciences (HESS). Please refer to the corresponding final paper in HESS. The published article is copyrighted by the author(s) and published by Copernicus Publications on behalf of the European Geosciences Union. The published article can be found at: http://hdl.handle.net/1957/5716