skip to main content

Title: Enhancing Soil and Water Assessment Tool Snow Prediction Reliability with Remote-Sensing-Based Snow Water Equivalent Reconstruction Product for Upland Watersheds in a Multi-Objective Calibration Process
Precipitation occurs in two basic forms defined as liquid state and solid state. Different from rain-fed watershed, modeling snow processes is of vital importance in snow-dominated watersheds. The seasonal snowpack is a natural water reservoir, which stores snow water in winter and releases it in spring and summer. The warmer climate in recent decades has led to earlier snowmelt, a decline in snowpack, and change in the seasonality of river flows. The Soil and Water Assessment Tool (SWAT) could be applied in the snow-influenced watershed because of its ability to simultaneously predict the streamflow generated from rainfall and from the melting of snow. The choice of parameters, reference data, and calibration strategy could significantly affect the SWAT model calibration outcome and further affect the prediction accuracy. In this study, SWAT models are implemented in four upland watersheds in the Tulare Lake Basin (TLB) located across the Southern Sierra Nevada Mountains. Three calibration scenarios considering different calibration parameters and reference datasets are applied to investigate the impact of the Parallel Energy Balance Model (ParBal) snow reconstruction data and snow parameters on the streamflow and snow water-equivalent (SWE) prediction accuracy. In addition, the watershed parameters and lapse rate parameters-led equifinality is also more » evaluated. The results indicate that calibration of the SWAT model with respect to both streamflow and SWE reference data could improve the model SWE prediction reliability in general. Comparatively, the streamflow predictions are not significantly affected by differently lumped calibration schemes. The default snow parameter values capture the extreme high flows better than the other two calibration scenarios, whereas there is no remarkable difference among the three calibration schemes for capturing the extreme low flows. The watershed and lapse rate parameters-induced equifinality affects the flow prediction more (Nash-Sutcliffe Efficiency (NSE) varies between 0.2–0.3) than the SWE prediction (NSE varies less than 0.1). This study points out the remote-sensing-based SWE reconstruction product as a promising alternative choice for model calibration in ungauged snow-influenced watersheds. The streamflow-reconstructed SWE bi-objective calibrated model could improve the prediction reliability of surface water supply change for the downstream agricultural region under the changing climate. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Snowpack provides the majority of predictive information for water supply forecasts (WSFs) in snow-dominated basins across the western United States. Drought conditions typically accompany decreased snowpack and lowered runoff efficiency, negatively impacting WSFs. Here, we investigate the relationship between snow water equivalent (SWE) and April–July streamflow volume (AMJJ-V) during drought in small headwater catchments, using observations from 31 USGS streamflow gauges and 54 SNOTEL stations. A linear regression approach is used to evaluate forecast skill under different historical climatologies used for model fitting, as well as with different forecast dates. Experiments are constructed in which extreme hydrological drought years are withheld from model training, that is, years with AMJJ-V below the 15th percentile. Subsets of the remaining years are used for model fitting to understand how the climatology of different training subsets impacts forecasts of extreme drought years. We generally report overprediction in drought years. However, training the forecast model on drier years, that is, below-median years (P15,P57.5], minimizes residuals by an average of 10% in drought year forecasts, relative to a baseline case, with the highest median skill obtained in mid- to late April for colder regions. We report similar findings using a modified National Resources Conservation Servicemore »(NRCS) procedure in nine large Upper Colorado River basin (UCRB) basins, highlighting the importance of the snowpack–streamflow relationship in streamflow predictability. We propose an “adaptive sampling” approach of dynamically selecting training years based on antecedent SWE conditions, showing error reductions of up to 20% in historical drought years relative to the period of record. These alternate training protocols provide opportunities for addressing the challenges of future drought risk to water supply planning.

    Significance Statement

    Seasonal water supply forecasts based on the relationship between peak snowpack and water supply exhibit unique errors in drought years due to low snow and streamflow variability, presenting a major challenge for water supply prediction. Here, we assess the reliability of snow-based streamflow predictability in drought years using a fixed forecast date or fixed model training period. We critically evaluate different training protocols that evaluate predictive performance and identify sources of error during historical drought years. We also propose and test an “adaptive sampling” application that dynamically selects training years based on antecedent SWE conditions providing to overcome persistent errors and provide new insights and strategies for snow-guided forecasts.

    « less
  2. Abstract

    Snow dominated mountainous karst watersheds are the primary source of water supply in many areas in the western U.S. and worldwide. These watersheds are typically characterized by complex terrain, spatiotemporally varying snow accumulation and melt processes, and duality of flow and storage dynamics because of the juxtaposition of matrix (micropores and small fissures) and karst conduits. As a result, predicting streamflow from meteorological inputs has been challenging due to the inability of physically based or conceptual hydrologic models to represent these unique characteristics. We present a hybrid modeling approach that integrates a physically based, spatially distributed, snow model with a deep learning karst model. More specifically, the high‐resolution snow model captures spatiotemporal variability in snowmelt, and the deep learning model simulates the corresponding response of streamflow as influenced by complex surface and subsurface properties. The deep learning model is based on the Convolutional Long Short‐Term Memory (ConvLSTM) architecture capable of handling spatiotemporal recharge patterns and watershed storage dynamics. The hybrid modeling approach is tested on a watershed in northern Utah with seasonal snow cover and variably karstified carbonate bedrock. The hybrid models were able to simulate streamflow at the watershed outlet with high accuracy. The spatial and temporal rechargemore »and discharge patterns learned by the ConvLSTM model were then examined and compared with known hydrogeologic information. Results suggest that ConvLSTM simulates streamflow with higher accuracy than reference models for the study area and provides insight into spatially influenced hydrologic responses that are unavailable within lumped modeling approaches.

    « less
  3. Abstract

    Canopy‐snow unloading is the complex physical process of snow unloading from the canopy through meltwater drip, sublimation to the atmosphere, or solid snow unloading to the snowpack below. This process is difficult to parameterize due to limited observations. Time‐lapse photographs of snow in the canopy were characterized by citizen scientists to create a data set of snow interception observations at multiple locations across the western United States. This novel interception data set was used to evaluate three snow unloading parameterizations in the Structure for Unifying Multiple Modeling Alternatives (SUMMA) modular hydrologic modeling framework. SUMMA was modified to include a third snow unloading parameterization, termed Wind‐Temperature (Roesch et al., 2001,, which includes wind‐dependent and temperature‐dependent unloading functions. It was compared to a meltwater drip unloading parameterization, termed Melt (Andreadis et al., 2009,, and a time‐dependent unloading parameterization, termed Exponential‐Decay (Hedstrom & Pomeroy, 1998,<1611::AID-HYP684>3.0.CO;2-4). Wind‐Temperature performed well without calibration across sites, specifically in cold climates, where wind dominates unloading and rime accretion is low. At rime prone sites, Wind‐Temperature should be calibrated to account for longer interception events with less sensitivity to wind, otherwise Melt can be used without calibration. The absence of model physics in Exponential‐Decay requires local calibrationmore »that can only be transferred to sites with similar unloading patterns. The choice of unloading parameterization can result in 20% difference in SWE on the ground below the canopy and 10% difference in estimated average winter canopy albedo. These novel observations shed light on processes that are often overlooked in hydrology.

    « less
  4. Abstract

    The selection of calibration and validation time periods in hydrologic modelling is often done arbitrarily. Nonstationarity can lead to an optimal parameter set for one period which may not accurately simulate another. However, there is still much to be learned about the responses of hydrologic models to nonstationary conditions. We investigated how the selection of calibration and validation periods can influence water balance simulations. We calibrated Soil and Water Assessment Tool hydrologic models with observed streamflow for three United States watersheds (St. Joseph River of Indiana/Michigan, Escambia River of Florida/Alabama, and Cottonwood Creek of California), using time period splits for calibration/validation. We found that the choice of calibration period (with different patterns of observed streamflow, precipitation, and air temperature) influenced the parameter sets, leading to dissimilar simulations of water balance components. In the Cottonwood Creek watershed, simulations of 50‐year mean January streamflow varied by 32%, because of lower winter precipitation and air temperature in earlier calibration periods on calibrated parameters, which impaired the ability for models calibrated to earlier periods to simulate later periods. Peaks of actual evapotranspiration for this watershed also shifted from April to May due to different parameter values depending on the calibration period's winter airmore »temperatures. In the St. Joseph and Escambia River watersheds, adjustments of the runoff curve number parameter could vary by 10.7% and 20.8%, respectively, while 50‐year mean monthly surface runoff simulations could vary by 23%–37% and 169%–209%, depending on the observed streamflow and precipitation of the chosen calibration period. It is imperative that calibration and validation time periods are chosen selectively instead of arbitrarily, for instance using change point detection methods, and that the calibration periods are appropriate for the goals of the study, considering possible broad effects of nonstationary time series on water balance simulations. It is also crucial that the hydrologic modelling community improves existing calibration and validation practices to better include nonstationary processes.

    « less
  5. Abstract

    Summer streamflow predictions are critical for managing water resources; however, warming‐induced shifts from snow to rain regimes impact low‐flow predictive models. Additionally, reductions in snowpack drive earlier peak flows and lower summer flows across the western United States increasing reliance on groundwater for maintaining summer streamflow. However, it remains poorly understood how groundwater contributions vary interannually. We quantify recession limb groundwater (RLGW), defined as the proportional groundwater contribution to the stream during the period between peak stream flow and low flow, to predict summer low flows across three diverse western US watersheds. We ask (a) how do snow and rain dynamics influence interannual variations of RLGW contributions and summer low flows?; (b) which watershed attributes impact the effectiveness of RLGW as a predictor of summer low flows? Linear models reveal that RLGW is a strong predictor of low flows across all sites and drastically improves low‐flow prediction compared to snow metrics at a rain‐dominated site. Results suggest that strength of RLGW control on summer low flows may be mediated by subsurface storage. Subsurface storage can be divided into dynamic (i.e., variability saturated) and deep (i.e., permanently saturated) components, and we hypothesize that interannual variability in dynamic storage contribution tomore »streamflow drives RLGW variability. In systems with a higher proportion of dynamic storage, RLGW is a better predictor of summer low flow because the stream is more responsive to dynamic storage contributions compared to deep‐storage‐dominated systems. Overall, including RLGW improved low‐flow prediction across diverse watersheds.

    « less