skip to main content


Title: Effects of spatial and temporal variability in surface water inputs on streamflow generation and cessation in the rain–snow transition zone
Abstract. Climate change affects precipitation phase, which can propagate into changes in streamflow timing and magnitude. This study examines how the spatial and temporal distribution of rainfall and snowmelt affects discharge in rain–snow transition zones. These zones experience large year-to-year variations in precipitation phase, cover a significant area of mountain catchments globally, and might extend to higher elevations under future climate change. We used observations from 11 weather stations and snow depths measured from one aerial lidar survey to force a spatially distributed snowpack model (iSnobal/Automated Water Supply Model) in a semiarid, 1.8 km2 headwater catchment. We focused on surface water input (SWI; the summation of rainfall and snowmelt on the soil) for 4 years with contrasting climatological conditions (wet, dry, rainy, and snowy) and compared simulated SWI to measured discharge. A strong spatial agreement between snow depth from the lidar survey and model (r2 = 0.88) was observed, with a median Nash–Sutcliffe efficiency (NSE) of 0.65 for simulated and measured snow depths at snow depth stations for all modeled years (0.75 for normalized snow depths). The spatial pattern of SWI was consistent between the 4 years, with north-facing slopes producing 1.09–1.25 times more SWI than south-facing slopes, and snowdrifts producing up to 6 times more SWI than the catchment average. Annual discharge in the catchment was not significantly correlated with the fraction of precipitation falling as snow; instead, it was correlated with the magnitude of precipitation and spring snow and rain. Stream cessation depended on total and spring precipitation, as well as on the melt-out date of the snowdrifts. These results highlight the importance of the heterogeneity of SWI at the rain–snow transition zone for streamflow generation and cessation, and emphasize the need for spatially distributed modeling or monitoring of both snowpack and rainfall dynamics.  more » « less
Award ID(s):
1653998
NSF-PAR ID:
10337000
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Hydrology and Earth System Sciences
Volume:
26
Issue:
10
ISSN:
1607-7938
Page Range / eLocation ID:
2779 to 2796
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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 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. 
    more » « less
  2. The goal of this project is to characterize and constrain the physical mechanisms that control snowmelt delivery to streams in headwater basins. This project leverages new observation and modeling techniques to quantify and simulate the snow distribution, water holding capacity, snowmelt production, and dynamic flowpaths. This is achieved through state-of-the-science observation techniques including ground penetrating radar (GPR), Terrestrial LiDAR Scanning (TLS), global positioning system (GPS) instrumentation, a network of sensor nodes continuously measuring soil moisture and snow depth, and a weir to monitor streamflow. Finally, hydrologic modeling will be conducted with the Structure for Unified Multiple Modeling Alternatives (SUMMA) model to assess the impact of modeling decisions and the ability to simulate snowmelt dynamics. The overarching research question of this project is: How do snowpack liquid water storage and through-snow hydrologic flowpaths affect hillslope-stream connectivity, and how do these processes evolve throughout the snowmelt season? This research question will be investigated in a snow-dominated headwater catchment. This work will observe and simulate the spatially and temporally variable snowmelt season to complete the following project objectives: O1) Map the dynamics of catchment snow water equivalent (SWE) using TLS surveys, GPR surveys, a network of sensor nodes, and manual observations. O2) Monitor the spatial and temporal progression of snowpack liquid water content and transport using combined TLS and GPR surveys, automated GPS signal attenuation, soil moisture sensors, and catchment streamflow response. O3) Evaluate the skill of hydrologic models to simulate the observed dynamics of the snowpack, soil, and streamflow response by systematically analyzing multiple model representations of hydrologic processes and scaling behavior. The work builds upon decades of local research in hydrology, biogeochemistry, and ecological processes. 
    more » « less
  3. The Yukon River basin encompasses over 832,000 km2 of boreal Arctic Alaska and northwest Canada, providing a major transportation corridor and multiple natural resources to regional communities. The river seasonal hydrology is defined by a long winter frozen season and a snowmelt-driven spring flood pulse. Capabilities for accurate monitoring and forecasting of the annual spring freshet and river ice breakup (RIB) in the Yukon and other northern rivers is limited, but critical for understanding hydrologic processes related to snow, and for assessing flood-related risks to regional communities. We developed a regional snow phenology record using satellite passive microwave remote sensing to elucidate interactions between the timing of upland snowmelt and the downstream spring flood pulse and RIB in the Yukon. The seasonal snow metrics included annual Main Melt Onset Date (MMOD), Snowoff (SO) and Snowmelt Duration (SMD) derived from multifrequency (18.7 and 36.5 GHz) daily brightness temperatures and a physically-based Gradient Ratio Polarization (GRP) retrieval algorithm. The resulting snow phenology record extends over a 29-year period (1988–2016) with 6.25 km grid resolution. The MMOD retrievals showed good agreement with similar snow metrics derived from in situ weather station measurements of snowpack water equivalence (r = 0.48, bias = −3.63 days) and surface air temperatures (r = 0.69, bias = 1 day). The MMOD and SO impact on the spring freshet was investigated by comparing areal quantiles of the remotely sensed snow metrics with measured streamflow quantiles over selected sub-basins. The SO 50% quantile showed the strongest (p < 0.1) correspondence with the measured spring flood pulse at Stevens Village (r = 0.71) and Pilot (r = 0.63) river gaging stations, representing two major Yukon sub-basins. MMOD quantiles indicating 20% and 50% of a catchment under active snowmelt corresponded favorably with downstream RIB (r = 0.61) from 19 river observation stations spanning a range of Yukon sub-basins; these results also revealed a 14–27 day lag between MMOD and subsequent RIB. Together, the satellite based MMOD and SO metrics show potential value for regional monitoring and forecasting of the spring flood pulse and RIB timing in the Yukon and other boreal Arctic basins. 
    more » « less
  4. Abstract

    Quantitative estimations of ecohydrological water partitioning into evaporation and transpiration remains mostly based on plot‐scale investigations that use well‐instrumented, small‐scale experimental catchments in temperate regions. Here, we attempted to upscale and adapt the conceptual tracer‐aided ecohydrology model STARRtropics to simulate water partitioning, tracer, and storage dynamics over daily time steps and a 1‐km grid larger‐scale (2565 km2) in a sparsely instrumented tropical catchment in Costa Rica. The model was driven by bias‐corrected regional climate model outputs and was simultaneously calibrated against daily discharge observations from 2 to 30 years at four discharge gauging stations and a 1‐year, monthly streamwater isotope record of 46 streams. The overall model performance for the best discharge simulations ranged in KGE values from 0.4 to 0.6 and correlation coefficients for streamflow isotopes from 0.3 to 0.45. More importantly, independent model‐derived transpiration estimates, point‐scale residence time estimates, and measured groundwater isotopes showed reasonable model performance and simulated spatial and temporal patterns pointing towards an overall model realism at the catchment scale over reduced performance in the headwaters. The simulated catchment system was dominated by low‐seasonality and high precipitation inputs and a marked topographical gradient. Climatic drivers overrode smaller, landcover‐dependent transpiration fluxes giving a seemingly homogeneous rainfall‐runoff dominance likely related to model input bias of rainfall isotopes, oversimplistic Potential Evapotranspiration (PET) estimates and averaged Leaf Area Index (LAI). Topographic influences resulted in more dynamic water and tracer fluxes in the headwaters that averaged further downstream at aggregated catchment scales. Modelled headwaters showed greater storage capacity by nearly an order of magnitude compared to the lowlands, which also favoured slightly longer residence times (>250 days) compared to superficially well‐connected groundwater contributing to shorter streamflow residence times (<150 days) in the lowlands. Our findings confirm that tracer‐aided ecohydrological modelling, even in the data‐scarce Tropics, can help gain a first, but crucial approximation of spatio‐temporal dynamics of how water is partitioned, stored and transported beyond the experimental catchment scale of only a few km2.

     
    more » « less
  5. Abstract

    Despite its location in the rain shadow of the southern Sierra Nevada, the Panamint Range hosts a complex mountain groundwater system supporting numerous springs which have cultural, historical, and ecological importance. The sources of recharge that support these quintessential desert springs remain poorly quantified since very little hydrogeological research has been completed in the Panamint Range. Here we address the following questions: (i) what is the primary source of recharge that supports springs in the Panamint Range (snowmelt or rainfall), (ii) where is the recharge occurring (mountain‐block, mountain‐front, or mountain‐system) and (iii) how much recharge occurs in the Panamint Range? We answer questions (i) and (ii) using stable isotopes measured in spring waters and precipitation, and question (iii) using a chloride mass‐balance approach which is compared to a derivation of the Maxey–Eakin equation. Our dataset of the stable isotopic composition (δ18O andδ2H) of precipitation is short (1.5 years), but analyses on spring water samples indicate that high‐elevation snowmelt is the dominant source of recharge for these springs, accounting for 57 (±9) to 79 (±12) percent of recharge. Recharge from rainfall is small but not insignificant. Mountain‐block recharge is the dominant recharge mechanism. However, two basin springs emerging along the western mountain‐front of the Panamint Range in Panamint Valley appear to be supported by mountain‐front and mountain‐system recharge, while Tule Spring (a basin spring emerging at the terminus of the bajada on the eastern side of the Panamint Range) appears to be supported by mountain‐front recharge. Calculated recharge rates range from 19 mm year−1(elevations < 1000 mrsl) to 388 mm year−1(elevations > 1000 mrsl). The average annual recharge is approximately 91 mm year−1(equivalent to 19.4 percent of total annual precipitation). We infer that the springs in the Panamint Range (and their associated ecosystems) are extremely vulnerable to changes in snow cover associated with climate change. They are heavily dependent on snowmelt recharge from a relatively thin annual snowpack. These findings have important implications for the vulnerability of desert springs worldwide.

     
    more » « less