Observed Changes in Climate and Streamflow in the Upper Rio Grande Basin

Observed streamflow and climate data are used to test the hypothesis that climate change is already affecting Rio Grande streamflow volume derived from snowmelt runoff in ways consistent with model‐based projections of 21st‐Century streamflow. Annual and monthly changes in streamflow volume and surface climate variables on the Upper Rio Grande, near its headwaters in southern Colorado, are assessed for water years 1958–2015. Results indicate winter and spring season temperatures in the basin have increased significantly, April 1 snow water equivalent (SWE) has decreased by approximately 25%, and streamflow has declined slightly in the April–July snowmelt runoff season. Small increases in precipitation have reduced the impact of declining snowpack on trends in streamflow. Changes in the snowpack–runoff relationship are noticeable in hydrographs of mean monthly streamflow, but are most apparent in the changing ratios of precipitation (rain + snow, and SWE) to streamflow and in the declining fraction of runoff attributable to snowpack or winter precipitation. The observed changes provide observational confirmation for model projections of decreasing runoff attributable to snowpack, and demonstrate the decreasing utility of snowpack for predicting subsequent streamflow on a seasonal basis in the Upper Rio Grande Basin.     

Effect of Snow Cover Conditions on the Hydrologic Regime: Case Study in a Pluvial‐Nival Watershed,Japan1

Abstract: Hydrologic monitoring in a small forested and mountainous headwater basin in Niigata Prefecture has been undertaken since 2000. An important characteristic of the basin is that the hydrologic regime contains pluvial elements year‐round, including rain‐on‐snow, in addition to spring snowmelt. We evaluated the effect of different snow cover conditions on the hydrologic regime by analyzing observed data in conjunction with model simulations of the snowpack. A degree‐day snow model is presented and applied to the study basin to enable estimation of the basin average snow water equivalent using air temperature at three representative elevations. Analysis of hydrological time series data and master recession curves showed that flow during the snowmelt season was generated by a combination of ground water flow having a recession constant of 0.018/day and diurnal melt water flow having a recession constant of 0.015/hour. Daily flows during the winter/snowmelt season showed greater persistence than daily flows during the warm season. The seasonal water balance indicated that the ratio of runoff to precipitation during the cold season (December to May) was about 90% every year. Seasonal snowpack plays an important role in defining the hydrologic regime, with winter precipitation and snowmelt runoff contributing about 65% of the annual runoff. The timing of the snowmelt season, indicated by the date of occurrence of the first significant snowmelt event, was correlated with the occurrence of low flow events. Model simulations showed that basin average snow water equivalent reached a peak around mid‐February to mid‐March, although further validation of the model is required at high elevation sites.     

SNOWPACK CHANGES RESULTING FROM TIMBER HARVEST INTERCEPTION,REDISTRIBUTION, AND EVAPORATION1

ABSTRACT: Three processes were examined as causing snowpack changes in forest clearings. Two of the three contribute to increases and one counteracts by reducing snowpack. The two that increase snowpack are redistribution and decreased loss to interception. Snow evaporation from a clearing counteracts snowpack increases. Research has indicated that as vegetation density increases, so too does the loss to interception. As snow in the canopy reaches the limit that the canopy can hold (the threshold amount) evaporation increases. Aerodynamics of the forest canopy were studied as well. As timber is cut, wind patterns are disturbed, creating disruptions in the wind velocity gradient depositing snow in openings. This redistribution leads to an increased snow water equivalent and augments runoff. Snow evaporation was shown to increase proportionally with opening size. Evaporation offsets the water yield gains derived from forest cut. It was found that this offset is inclusive to the measurements of water yield changes in experimental forests. An optimal size of harvest block may be five tree heights in width as suggested by numerous studies.     

Application of Satellite Microwave Images in Estimating Snow Water Equivalent1

Abstract: Flood forecast and water resource management requires reliable estimates of snow pack properties [snow depth and snow water equivalent (SWE)]. This study focuses on application of satellite microwave images to estimate the spatial distribution of snow depth and SWE over the Great Lakes area. To estimate SWE, we have proposed the algorithm which uses microwave brightness temperatures (Tb) measured by the Special Sensor Microwave Imager (SSM/I) radiometer along with information on the Normalized Difference Vegetation Index (NDVI). The algorithm was developed and tested over 19 test sites characterized by different seasonal average snow depth and land cover type. Three spectral signatures derived from SSM/I data, namely T19V‐T37V (GTV), T19H‐T37H (GTH), and T22V‐T85V (SSI), were examined for correlation with the snow depth and SWE. To avoid melting snow conditions, we have used observations taken only during the period from December 1‐February 28. It was found that GTH, and GTV exhibit similar correlation with the snow depth/SWE and are most should be used over deep snowpack. In the same time, SSI is more sensitive to snow depth variations over a shallow snow pack. To account for the effect of dense forests on the scattering signal of snow we established the slope of the regression line between GTV and the snow depth as a function of NDVI. The accuracy of the new technique was evaluated through its comparison with ground‐based measurements and with results of SWE analysis prepared by the National Operational Hydrological Remote Sensing Center (NOHRSC) of the National Weather Service. The proposed algorithm was found to be superior to previously developed global microwave SWE retrieval techniques.     

ESTIMATING REGIONAL SNOW WATER EQUIVALENT WITH A SIMPLE SIMULATION MODEL1

A simple simulation model designed to monitor snow-packs of the central Sierra Nevada is described. The model estimates average snow water equivalent for rectangular subregions in the area. Static subregion characteristics, daily precipitation and mean and minimum air temperatures measured at three index stations are the only needed input values. A water balance technique simulates daily snowpack changes in each subregion. Reasonable basinwide water equivalent values are produced. The procedure should be useful for estimating snow water distribution in large mountainous watersheds.     

GENERAL-CIRCULATION-MODEL SIMULATIONS OF FUTURE SNOWPACK IN THE WESTERN UNITED STATES1

ABSTRACT: April 1 snowpack accumulations measured at 311 snow courses in the western United States (U.S.) are grouped using a correlation-based cluster analysis. A conceptual snow accumulation and melt model and monthly temperature and precipitation for each cluster are used to estimate cluster-average April 1 snowpack. The conceptual snow model is subsequently used to estimate future snowpack by using changes in monthly temperature and precipitation simulated by the Canadian Centre for Climate Modeling and Analysis (CCC) and the Hadley Centre for Climate Prediction and Research (HADLEY) general circulation models (GCMs). Results for the CCC model indicate that although winter precipitation is estimated to increase in the future, increases in temperatures will result in large decreases in April 1 snowpack for the entire western U.S. Results for the HADLEY model also indicate large decreases in April 1 snowpack for most of the western US, but the decreases are not as severe as those estimated using the CCC simulations. Although snowpack conditions are estimated to decrease for most areas of the western US, both GCMs estimate a general increase in winter precipitation toward the latter half of the next century. Thus, water quantity may be increased in the western US; however, the timing of runoff will be altered because precipitation will more frequently occur as rain rather than as snow.     

Hydrologic Calibration and Validation of SWAT in a Snow‐Dominated Rocky Mountain Watershed,Montana, U.S.A.1

Abstract: The Soil and Water Assessment Tool (SWAT) has been applied successfully in temperate environments but little is known about its performance in the snow‐dominated, forested, mountainous watersheds that provide much of the water supply in western North America. To address this knowledge gap, we configured SWAT to simulate the streamflow of Tenderfoot Creek (TCSWAT). Located in central Montana, TCSWAT represents a high‐elevation watershed with ～85% coniferous forest cover where more than 70% of the annual precipitation falls as snow, and runoff comes primarily from spring snowmelt. Model calibration using four years of measured daily streamflow, temperature, and precipitation data resulted in a relative error (RE) of 2% for annual water yield estimates, and mean paired deviations (Dv) of 36 and 31% and Nash‐Sutcliffe (NS) efficiencies of 0.90 and 0.86 for monthly and daily streamflow, respectively. Model validation was conducted using an additional four years of data and the performance was similar to the calibration period, with RE of 4% for annual water yields, Dv of 43% and 32%, and NS efficiencies of 0.90 and 0.76 for monthly and daily streamflow, respectively. An objective, regression‐based model invalidation procedure also indicated that the model was validated for the overall simulation period. Seasonally, SWAT performed well during the spring and early summer snowmelt runoff period, but was a poor predictor of late summer and winter base flow. The calibrated model was most sensitive to snowmelt parameters, followed in decreasing order of influence by the surface runoff lag, ground water, soil, and SCS Curve Number parameter sets. Model sensitivity to the surface runoff lag parameter reflected the influence of frozen soils on runoff processes. Results indicated that SWAT can provide reasonable predictions of annual, monthly, and daily streamflow from forested montane watersheds, but further model refinements could improve representation of snowmelt runoff processes and performance during the base flow period in this environment.     

Management Implications of Snowpack Sensitivity to Temperature and Atmospheric Moisture Changes in Yosemite National Park,CA

In order to investigate snowpack sensitivity to temperature increases and end‐member atmospheric moisture conditions, we applied a well‐constrained energy‐ and mass‐balance snow model across the full elevation range of seasonal snowpack using forcing data from recent wet and dry years. Humidity scenarios examined were constant relative humidity (high) and constant vapor pressure between storms (low). With minimum calibration, model results captured the observed magnitude and timing of snowmelt. April 1 snow water equivalent (SWE) losses of 38%, 73%, and 90% with temperature increases of 2, 4, and 6°C in a dry year centered on areas of greatest SWE accumulation. Each 2°C increment of warming also resulted in seasonal snowline moving upslope by 300 m. The zone of maximum melt was compressed upward 100–500 m with 6°C warming, with the range reflecting differences in basin hypsometry. Melt contribution by elevations below 2,000 m disappeared with 4°C warming. The constant‐relative‐humidity scenario resulted in 0–100 mm less snowpack in late spring vs. the constant‐vapor‐pressure scenario in a wet year, a difference driven by increased thermal radiation (+1.2 W/m²) and turbulent energy fluxes (+1.2 W/m²) to the snowpack for the constant‐relative‐humidity case. Loss of snowpack storage and potential increases in forest evapotranspiration due to warming will result in a substantial shift in forest water balance and present major challenges to land management in this mountainous region.     

OPERATIONAL APPLICATIONS OF SATELLITE SNOW COVER OBSERVATIONS1

ABSTRACT: Several federal and state water resources agencies and NASA have recently completed an Applications Systems Verification and Transfer (ASVT) project on the operational applications of satellite snow cover observations. When satellite snow cover data were tested in both empirical seasonal runoff estimation and short term modeling approaches, a definite potential for reducing forecast error was evident. Three years of testing in California resulted in reduction of seasonal stream flow forecast error was evident. Three years of testing in California resulted in reduction of seasonal stream flow forecast error from 15 percent to 10 percent on three study basins; and modeling studies on the Boise River basin in Idaho indicated that satellite snow cover could be used to reduce short term forecast error by up to 9.6 percent (5 day forecast). Potential benefits from improved satellite snow cover based predictions across the 11 western states total 10 million dollars for hydropower and 28 million dollars for irrigation annually. The truly operational application of the new technology in the West, however, will only be possible when the turnaround time for all data is reduced to 72 hours, and the water management agencies can be assured of a continuing supply of operational snow cover data from space.     

Snowmelt Runoff and Water Yield Along Elevation and Temperature Gradients in California’s Southern Sierra Nevada1

Hunsaker, Carolyn T., Thomas W. Whitaker, and Roger C. Bales, 2012. Snowmelt Runoff and Water Yield Along Elevation and Temperature Gradients in California’s Southern Sierra Nevada. Journal of the American Water Resources Association (JAWRA) 48(4): 667‐678. DOI: 10.1111/j.1752‐1688.2012.00641.x Abstract: Differences in hydrologic response across the rain‐snow transition in the southern Sierra Nevada were studied in eight headwater catchments – the Kings River Experimental Watersheds – using continuous precipitation, snowpack, and streamflow measurements. The annual runoff ratio (discharge divided by precipitation) increased about 0.1 per 300 m of mean catchment elevation over the range 1,800‐2,400 m. Higher‐elevation catchments have lower vegetation density, shallow soils with rapid permeability, and a shorter growing season when compared with those at lower elevations. Average annual temperatures ranged from 6.8°C at 2,400 m to 8.6 at 1,950 m elevation, with annual precipitation being 75‐95% snow at the highest elevations vs. 20‐50% at the lowest. Peak discharge lagged peak snow accumulation on the order of 60 days at the higher elevations and 20 to 30 days at the lower elevations. Snowmelt dominated the daily streamflow cycle over a period of about 30 days in higher elevation catchments, followed by a 15‐day transition to evapotranspiration dominating the daily streamflow cycle. Discharge from lower elevation catchments was rainfall dominated in spring, with the transition to evapotranspiration dominance being less distinct. Climate warming that results in a longer growing season and a shift from snow to rain would result in earlier runoff and a lower runoff ratio.     

Spatiotemporal Variability of Snow Depletion Curves Derived from SNODAS for the Conterminous United States, 2004‐2013

Assessment of water resources at a national scale is critical for understanding their vulnerability to future change in policy and climate. Representation of the spatiotemporal variability in snowmelt processes in continental‐scale hydrologic models is critical for assessment of water resource response to continued climate change. Continental‐extent hydrologic models such as the U.S. Geological Survey National Hydrologic Model (NHM) represent snowmelt processes through the application of snow depletion curves (SDCs). SDCs relate normalized snow water equivalent (SWE) to normalized snow covered area (SCA) over a snowmelt season for a given modeling unit. SDCs were derived using output from the operational Snow Data Assimilation System (SNODAS) snow model as daily 1‐km gridded SWE over the conterminous United States. Daily SNODAS output were aggregated to a predefined watershed‐scale geospatial fabric and used to also calculate SCA from October 1, 2004 to September 30, 2013. The spatiotemporal variability in SNODAS output at the watershed scale was evaluated through the spatial distribution of the median and standard deviation for the time period. Representative SDCs for each watershed‐scale modeling unit over the conterminous United States (n = 54,104) were selected using a consistent methodology and used to create categories of snowmelt based on SDC shape. The relation of SDC categories to the topographic and climatic variables allow for national‐scale categorization of snowmelt processes.     

Trends in Raindrop Kinetic Energy with Modeled Climate Warming in the Lake Tahoe Basin

Two means by which climate change may increase surface soil erosion in mountainous terrain are: (1) increasing the frequency of extreme rainfall events and (2) decreasing the duration of snow cover on bare soil. We used output from four general circulation models (GCMs) and two greenhouse gas trajectories to produce a suite of hydrologic variables at a daily time‐step for historic and projected 21st Century conditions. We statistically disaggregated the daily rainfall to hourly, using hourly rainfall from a network of nine weather stations in the Tahoe Basin, and filtered out rain falling on a snowpack. We applied published equations to convert hourly intensity to raindrop kinetic energy (KE) for each day and grid cell in the Basin, averaged across grid cells, and created time series of total annual and maximum annual hourly kinetic energy (TKE and MKE) on snow‐free ground. Using the Generalized Extreme Value distribution, we calculated the significance of long‐term trends in KE on snow‐free ground, and estimated energy levels for return periods of 2, 20, and 100 years. We then detrended the snowpack data and compared the resulting trends in KE with the trends resulting from changes in both rainfall energy and snowpack under two GCMs. Principal findings include (1) upward trends in MKE, (2) stronger upward trends in TKE; and (3) an effect of increasing rainfall intensities on KE in some cases, and a strong effect of reduced snowpack in all cases examined.     

SNOW: NATURE'S RESERVOIR1

ABSTRACT: Snow, one of Nature's greatest reservoirs, supplies most of the usable water in the Western United States. Reliable predictions of the quantity and timing of the release of this water are used in making management decisions involving irrigation, stock water and municipal water supplies, hydro-power generation, recreation, navigation, and pollution control Practically oriented research is vital for the proper development and management of this resource. In southwestern Idaho, the Northwest Watershed Research Center, ARS, USDA, is conducting intensive investigations for assessing snow Volumes, snow water content, and snow-melt over a watershed. Application of these research findings will result in better development and management of the water stored as snow in Nature's reservoir.     

CLIMATE,SNOW COVER,GLACIERS, AND RUNOFF IN THE TIEN SHAN,CENTRAL ASIA1

ABSTRACT: The Pica Shan, a mountainous region located on the northern periphery of central Asia, has a wide range of climatic and hydrological conditions. On the basis of long term data from 348 meteorological and glaciological stations, the annual distribution of precipitation in different regions and elevational zones of the Tien Shan was calculated. Major climatic features are the entrance of moisture during spring-summer, small winter precipitation, decrease of precipitation towards the east and the center of the mountains or with distance up valleys, and increase of precipitation with altitude up to crest-lines of ranges. Annual total evaporation from snow can be 50–60 mm per year, reaching 30 percent of snow accumulation. Four main groups of rivers were identified: rivers with mainly snow nourishment, rivers with mainly glacial nourishment, rivers with mainly rain nourishment, and rivers with mainly ground water nourishment. Coefficient of runoff variation in Tien Shan's rivers is about 0.20, and coefficient of glacial runoff variation is about 0.15. Glacial runoff is 15–20 percent of the total volume of river runoff.     

EFFECTS OF CLIMATE CHANGE ON HYDROLOGY AND WATER RESOURCES IN THE COLUMBIA RIVER BASIN1

ABSTRACT: As part of the National Assessment of Climate Change, the implications of future climate predictions derived from four global climate models (GCMs) were used to evaluate possible future changes to Pacific Northwest climate, the surface water response of the Columbia River basin, and the ability of the Columbia River reservoir system to meet regional water resources objectives. Two representative GCM simulations from the Hadley Centre (HC) and Max Planck Institute (MPI) were selected from a group of GCM simulations made available via the National Assessment for climate change. From these simulations, quasi-stationary, decadal mean temperature and precipitation changes were used to perturb historical records of precipitation and temperature data to create inferred conditions for 2025, 2045, and 2095. These perturbed records, which represent future climate in the experiments, were used to drive a macro-scale hydrology model of the Columbia River at 1/8 degree resolution. The altered streamflows simulated for each scenario were, in turn, used to drive a reservoir model, from which the ability of the system to meet water resources objectives was determined relative to a simulated hydrologic base case (current climate). Although the two GCM simulations showed somewhat different seasonal patterns for temperature change, in general the simulations show reasonably consistent basin average increases in temperature of about 1.8–2.1°C for 2025, and about 2.3–2.9°C for 2045. The HC simulations predict an annual average temperature increase of about 4.5°C for 2095. Changes in basin averaged winter precipitation range from -1 percent to + 20 percent for the HC and MPI scenarios, and summer precipitation is also variously affected. These changes in climate result in significant increases in winter runoff volumes due to increased winter precipitation and warmer winter temperatures, with resulting reductions in snowpack. Average March 1 basin average snow water equivalents are 75 to 85 percent of the base case for 2025, and 55 to 65 percent of the base case by 2045. By 2045 the reduced snowpack and earlier snow melt, coupled with higher evapotranspiration in early summer, would lead to earlier spring peak flows and reduced runoff volumes from April-September ranging from about 75 percent to 90 percent of the base case. Annual runoff volumes range from 85 percent to 110 percent of the base case in the simulations for 2045. These changes in streamflow create increased competition for water during the spring, summer, and early fall between non-firm energy production, irrigation, instream flow, and recreation. Flood control effectiveness is moderately reduced for most of the scenarios examined, and desirable navigation conditions on the Snake are generally enhanced or unchanged. Current levels of winter-dominated firm energy production are only significantly impacted for the MPI 2045 simulations.     

iTree‐Hydro: Snow Hydrology Update For The Urban Forest Hydrology Model1

Yang, Yang, Theodore A. Endreny, and David J. Nowak, 2011. iTree‐Hydro: Snow Hydrology Update for the Urban Forest Hydrology Model. Journal of the American Water Resources Association (JAWRA) 47(6):1211–1218. DOI: 10.1111/j.1752‐1688.2011.00564.x Abstract: This article presents snow hydrology updates made to iTree‐Hydro, previously called the Urban Forest Effects—Hydrology model. iTree‐Hydro Version 1 was a warm climate model developed by the USDA Forest Service to provide a process‐based planning tool with robust water quantity and quality predictions given data limitations common to most urban areas. Cold climate hydrology routines presented in this update to iTree‐Hydro include: (1) snow interception to simulate the capture of snow by the vegetation canopy, (2) snow unloading to simulate the release of snow triggered by wind, (3) snowmelt to simulate the solid to liquid phase change using a heat budget, and (4) snow sublimation to simulate the solid to gas phase via evaporation. Cold climate hydrology routines were tested with research‐grade snow accumulation and weather data for the winter of 1996‐1997 at Umpqua National Forest, Oregon. The Nash‐Sutcliffe efficiency for open area snow accumulation was 0.77 and the Nash‐Sutcliffe efficiency for under canopy was 0.91. The USDA Forest Service offers iTree‐Hydro for urban forest hydrology simulation through http://www.iTreetools.org .     

AN ELEVATIONAL CONTROL OF PEAK SNOWPACK VARIABILITY1

ABSTRACT: The records of the seasonal peaks in snow accumulation on 24 snow courses show that the relative variability (measured by the coefficient of variation) is inversely related to elevation in the San Juan Mountains, southwestern Colorado. Analysis on an annual basis shows that this is due to a tendency for the peak snowpack at high elevations to be closer to the long-term mean while that at low levels is further from the normal. This is true in both above-normal and below-normal accumulation seasons. Extrapolation of 20 annual elevation-accumulation trends suggests that mean accumulation is reached at 4107 ± 186 m elevation. This is approximately the height of the topographic harrier of the San Juan Mountains suggesting that the pattern of variability is a partial function of an atmospheric-topographic interaction.     

ESTIMATED RUNOFF FROM MAN-MADE SNOW1

ABSTRACT: The consumptive loss from man-made snowmaking at six Colorado ski areas is calculated. The focus of the procedures in this investigation is on the consumptive loss that occurs to man-made snow particles during the period they reside on or in the snowpack until spring snowmelt (termed the watershed loss). Calculated watershed losses under a variety of precipitation and temperature conditions at six ski areas varied from 7 to 33 percent. These calculations were made using the calibrated Subalpine Water Balance Simulation Model (Leaf and Brink, 1973a, 1973b). The watershed loss of 7 to 33 percent indicates the range of likely watershed losses that can be expected at Colorado ski areas. A previous paper by the authors (Eisel et al., 1988) provided estimates of the mean consumptive loss during the snowmaking process (termed initial loss) for conditions existing at Colorado ski areas to be 6 percent of water applied. Therefore, based on the mean initial loss, the total consumptive loss from man-made snowmaking under conditions found at Colorado ski areas could be expected to range from 13 to 37 percent. These results demonstrate the range of total consumptive losses that could be expected in various years and for various watershed conditions. These total percentage losses cannot be extrapolated directly to other specific sites because the total consumptive loss is dependent on temperature during actual snowmaking, temperature and precipitation throughout the winter at the specific ski area, and watershed conditions at the ski area. Consumptive losses to man-made snow for a specific ski area should be estimated using the handbook procedures developed especially for this purpose (Colorado Ski Country USA, 1986b).     

Defining the Diurnal Pattern of Snowmelt Using a Beta Distribution Function

Snow is an important component of the hydrologic cycle for many regions worldwide. In addition to vital water resources, snowmelt can be important for forest ecosystem dynamics and flood risk. However, standard design events in the United States lack a design snowmelt event, including only precipitation events, though snowmelt has been shown to be larger than rainfall. In this article, we present a method using hourly snow water equivalent data to develop and test a function for representing the diurnal pattern of snowmelt. A two‐parameter beta distribution function is modified for the purposes of this study and found to fit the pattern of snowmelt well with a root mean squared error of 0.008. Soil moisture sensors were additionally utilized to assess the timing of the snowmelt water outflow from the base of the snowpack that supports the shape of the function, but suggests that the timing of losses recorded on snow pillows lag as much as 3 h. Further testing of the function showed the shape of the function to be accurate. The methods developed and tested in this paper can be applied for design purposes comparing snowmelt and rainfall events or to improve hydrological models investigating processes such as streamflow or groundwater recharge.