Snow water equivalent along elevation gradients in the Merced and Tuolumne River basins of the Sierra Nevada
|
|
- Martin McDowell
- 5 years ago
- Views:
Transcription
1 Manuscript for Water Resource Research Rice et al. May 2011 Snow water equivalent along elevation gradients in the Merced and Tuolumne River basins of the Sierra Nevada Robert Rice Sierra Nevada Research Institute University of California, Merced Roger C. Bales Sierra Nevada Research Institute University of California, Merced Thomas H. Painter Jet Propulsion Laboratory California Institute of Technology Jeff Dozier Bren School of Environmental Science & Management University of California, Santa Barbara Corresponding Author and Address Robert Rice Sierra Nevada Research Institute University of California 5200 N. Lake Road Merced, CA95343 Phone: 209/
2 Abstract We used daily remotely sensed fractional snow-covered area (SCA) at 500-m resolution to estimate snow water equivalent (SWE) across the Upper Merced and Tuolumne River basins of the Sierra Nevada of California for 2004 (dry and warm) and 2005 (wet and cool). From 1800 to 3900 m, each successively higher 300 m elevation band consistently melts out 2-3 weeks later than the one below it. We compared two methods of estimating SWE from SCA: (i) blending the fractional SCA with SWE interpolated from snow-pillow measurements; and (ii) retrospectively estimating cumulative snowmelt based on a degree-day calculation after the snow disappeared. The interpolation approach estimates a lower snowmelt volume above 3000 m and a higher snowmelt contribution at elevations between m. Snowmelt timing from the depletion approach matches observed streamflow timing much better than snowmelt estimated by the interpolation method. The snow-pillow sites used in the interpolation method do not cover the highest elevations and melted out several weeks before the basin itself was free of snow. Middle elevations ( m) contributed 40-60% of the annual snowmelt in both basins, the lower elevations ( m) 10-15%, and elevations above 3000 m the remaining 30-40%. The presence of snow in the highest elevations highlights their critical buffering effect in accumulating snow every year. Variability in lower-elevation snow illustrates its susceptibility to climate variability and change Index terms: snowmelt (0740), snow and ice (1863), estimation and forecasting (1816), remote sensing (1855)
3 Introduction In mid-latitude montane regions, snow is critical for seasonal storage of water, releasing winter precipitation in spring and summer to provide the soil moisture and streamflow needed to sustain ecosystems and human populations. In the semi-arid western United States, snowmelt runoff accounts for up to 80% of annual streamflow [Daly et al., 2001]. Historically, most forecasts of seasonal runoff across this region have been based on statistical relationships between historical runoff and snow water equivalent (SWE) measured at manual snow courses and automatically telemetered sites, but these relationships have been developed during a period when climate was changing so their reliability will degrade as climate continues to change. Moreover, these index sites (i.e. snow courses) often fail to provide spatially representative measures of SWE and do not capture the physiographic variability across a basin [Molotch and Bales, 2005; Dressler et al., 2006; Rice and Bales, 2010]. Snowfall dominates precipitation in the mountains of the western U.S.; and snow accumulation varies because of topography, vegetation cover, and larger-scale synoptic processes [Bales et al., 2006], resulting in snow measurements that exhibit considerable variability, even between sites that are close together [Carroll et al., 1999]. Therefore, accurate estimates of the actual volume of snow and its spatial distribution are needed to provide quantitative estimates for more robust water-supply estimates, flood forecasts, and resource-management decisions. Blending spatial SWE estimates developed from interpolation of ground-based sensors with satellite-derived estimates of snow-covered area (SCA) potentially provides a more representative estimate of SWE across a basin than interpolation alone, because the snow-cover information constrains the interpolation [Bales et al., 2008; Bavera and De Michele, 2009; Harshburger et al., 2010]. Both the SCA interpolation and the surface measurements can be 1
4 made available in near real time, offering a promising way to provide spatial snow data to drive hydrologic forecast models. In this approach, snowmelt is computed as the difference between successive interpolated spatial SWE estimates. However, these SWE interpolations have not been evaluated at the basin scale, owing to the lack of spatially distributed ground truth estimates. Fractional SCA from NASA s Moderate Resolution Imaging Spectroradiometer (MODIS), and in the future with the Visible Infrared Imaging Radiometer Suite (VIIRS), makes possible a real-time product for operational hydrology [Dozier and Frew, 2009; Painter et al., 2009]. Martinec and Rango [1981] developed an alternate approach for estimating the spatial distribution of SWE, by combining a retrospective time series of snowmelt estimates with a time series of spatial SCA, and back-calculating the amount of SWE that was present earlier in the melt season. Because snowfall sometimes occurs during the snowmelt season, the SWE can be estimated only back to the last significant snowfall. This depletion approach provides an independent spatial SWE and snowmelt time series, which, although available only after the snow is gone, can be compared to the interpolation approach. Subsequent investigations using energy-balance snowmelt models have validated the accuracy of the method [Cline et al., 1998; Liston, 1999], and Molotch et al.[2004] point out that the energy-balance estimates are sensitive to accurate measurements of snow albedo. Homan et al. [2011] recently showed that an energybalance snowmelt model, together with a simple estimate of fractional SCA from MODIS based on a normalized difference snow index, can reproduce single-pixel values of SCA and SWE observed in the field over a snowmelt season, supporting this approach for estimating basin-wide snowmelt depletion curves. In an earlier study, Shamir and Georgakakos [2006] developed 2
5 basin-wide depletion curves for the American River using a temperature-index model and showed that model results compare well with a binary SCA product from MODIS. The aims of the research reported in this paper are to: (i) determine the patterns of snowcovered area in two Sierra Nevada basins based on a fractional SCA product, (ii) evaluate the adequacy of blending ground-based measurements with the fractional SCA to produce estimates of spatially distributed SWE and snowmelt; and (iii) determine the relative importance of different elevation zones within a basin to seasonally integrated snowmelt generation. Methods and data We carried out the study for 2004 and 2005 in the upper Tuolumne and Merced River basins, located in the Sierra Nevada above the foothill reservoirs of Don Pedro (Tuolumne) and Lake McClure (Merced) (Figure 1). These two years represent one below- and one above-normal year. In 2004, SWE measured at snow courses in the Tuolumne (Merced) was 83% (84%) of the historical April 1 records, where as in 2005 the values were 163% of normal in both basins. The Tuolumne basin s area is 4182 km 2 with elevations of m, while the Merced s is 2844 km 2 with elevations of m. Both rivers drain into the San Joaquin River in the Central Valley. Further, both basins are largely free from current human influences like dams, diversions, and major land-use changes, and the basins steep slopes and shallow soils make the hydrology straightforward [Slack and Landwehr, 1992; Jeton et al., 1996]. We restricted this analysis to seasonally snow-covered areas, which start at the rain/snow transition of 1500 m, representing 58% (2420 km 2 ) and 50% (1403 km 2 ) of the Tuolumne and Merced basins above their foothill reservoirs. Much of the snowfall lies within Yosemite National Park. Using topographic data from the Shuttle Radar Topography Mission at 30-m spatial sampling [Farr et 3
6 al., 2007], the study area was partitioned into eight elevation bands of 300-m increments, beginning at 1500 m and extending to 4000 m (Table 1). Daily fractional SCA maps at 500-m resolution were generated from the MODIS surfacereflectance product MOD09. SCA is retrieved with a spectral mixing model, through which the snow, soil, and vegetation fractions in each grid cell on each day are estimated [Dozier and Painter, 2004; Painter et al., 2009]. The MODIS snow-cover products were produced from the daily Terra satellite with a morning overpass about 10:30 am local time. Because cloud cover and low viewing angle limit the availability of quality data in some scenes, an interpolation scheme was applied to develop a spatially and temporally continuous snow product [Dozier et al., 2008]. No adjustment of snow cover was made for the occlusion by tree canopy, so the SCA values represent projected snow cover, i.e. snow cover that the satellite s sensor can detect. The fractional SCA threshold was set at 0.15, i.e. values below 15% were set to zero to prevent identification of spurious snow. Spatially distributed snowmelt was estimated daily for 2004 and 2005 by both methods described in the Introduction: interpolation and snow-cover depletion. The interpolation involved point SWE measurements from snow-pillows, operated by the California Department of Water Resources and U.S. Natural Resource Conservation Service, followed by masking with the MODIS fractional SCA product [Fassnacht et al., 2003; Bales et al., 2008]. Masking involves multiplying the pixel-by-pixel MODIS SCA product by the interpolated SWE product. SWE for each 500-m grid cell was interpolated using a linear regression between elevation and SWE for all snow-pillows within a 200-km radius of that pixel, including those outside the basin. A residual was obtained at each grid cell where an observing snow pillow was located by removing the observed value from the analysis (i.e. jack-knifing) and subtracting the observed 4
7 SWE from the computed SWE. Elevation-dependent biases in the residuals were removed by regressing residuals to a fixed datum of 5000 m using a lapse rate. Once regressed to the common datum, the lapsed residuals were spatially distributed using inverse distance squared weighting. The gridded residual surface was then regressed back to the basin surface using the same lapse rate and subtracted from the hypsometrically derived SWE grid in order to derive the SWE surface, thus preserving the SWE observation at each station. Daly et al. [2001] used a similar approach but computed one hypsometric relationship for each sub-basin instead of using a moving search radius to compute the hypsometric relationship at each pixel. Daily snowmelt was calculated by the daily grid cell differences in SWE. The alternate depletion approach to estimating snowmelt combined remote sensing, groundbased SWE and temperatures, and a temperature-index model [Anderson, 1968; Rango and Martinec, 1995; Ohmura, 2001]: M = a (T a T b ) (1) M is daily snowmelt, a is a degree-day coefficient (mm deg -1 day -1 ), T a is average daily temperature, T b = 0 o C; when T a <T b, no melt occurs. In this method, the ground-based SWE data were used to estimate degree-day coefficients. The approach provides an index of the average energy flux, but does not explicitly consider the individual fluxes and controlling factors that influence snowmelt (e.g. solar radiation, albedo, topography, turbulent-energy exchanges). Since temperature and SWE vary from year to year, producing different interannual daily melt rates, a daily degree-day coefficient for each year was calculated for both the Tuolumne and Merced using each of the nine co-located snow-pillow and temperature measurements (Figure 1); data were from the California Data Exchange Center (CDEC, Daily average temperature and SWE 5
8 were calculated from hourly and daily CDEC data for each of the nine snow pillows. The daily degree-day coefficient was calculated as the ratio of the daily decrease in snowmelt and the daily degree day [Anderson, 1968]. Our method of parameterizing the daily degree-day coefficient is similar to such semidistributed runoff models as the Snowmelt Runoff Model [Kustas et al., 1994] and Snow-17 [Anderson, 1973]. The degree-day coefficient increased over the ablation period in both years (Figure 2). As the two basins are adjacent and have similar physiographic characteristics, the same set of daily values was used for both. However, since the snow covered part of the basin extends over more than 2000 m, daily coefficients were calculated for 4 of the 8 elevations bands between m, where station data were available (Figure 3). Stations were grouped by elevation band and a linear trend fit to the values. In addition, in order to reduce the effects of possible site-specific differences (i.e. local shading, vegetation) snow pillow-stations within ±80 m of an elevation band were included to increase the sample size and reduce the bias. Below 1800 m and above 3000 m, daily coefficients from the adjacent elevation band were applied; above 2700 m no station data are available in the Merced, so the Tuolumne data were used. Two snow pillows were not used because the data were incomplete in The temperature-index calculation used daily average temperature for each 300-m elevation band. Stations within an elevation band were used to estimate degree days for that band. A simple lapse rate based on a linear fit to the daily station data was used to estimate the temperature in elevation bands with no station data. The calculated fixed lapse rates across the full elevation range for both basins were -6.0 and -7.0 o C per km for 2004 and We recognize that a simple lapse rate may not fully describe the spatial structure of temperature across complex terrain; a more involved interpolation scheme can improve spatially distributed 6
9 models incorporating microclimates [Lundquist and Cayan, 2007; Lundquist et al., 2008]. However, for whole-basin calculations integrating the full range of physiographic variability, a simple fixed lapse rate should not bias results. Daily snowmelt, M, was calculated for each elevation band using equation (1) and the average fractional SCA for each band was applied to correct for snow-free area. That is, if an elevation band had snow, we assumed it contributed a melt volume equal to the snowmelt rate times SCA. Results The degree-day coefficient, a, increased over the ablation season, with values for 2004 (2005) ranging from 1.2 (3.4) mm day -1 deg -1 at the onset of snowmelt to 7.8 (6.1) mm day -1 deg -1 by the end of the ablation season, but showed a high degree of daily variability (Figure 2). The different patterns for the degree-day coefficients in 2004 versus 2005 reflect the differing timing of melt. Snowmelt at the snow-pillow sites started much earlier in 2004 (early March) and was largely complete in mid-may, a month earlier than in Following Anderson [1968], the degree-day coefficient increased up to the summer solstice (June 21), with values afterward decreasing. In 2004 the degree-day coefficients show a distinct increase in elevation along the 4 elevation bands between m, where the snow-pillow sites are located (Figure 3). In 2005, the degree-day coefficients differed across the 4 elevation bands, but the elevation dependence was not as strong as it was in 2004, with the m elevation band showing a lower coefficient than the m elevation band. The daily degree-day coefficients were used to back calculate snowmelt at each snow-pillow site to determine whether these values were representative of the snowmelt. At all snow-pillow sites the elevation-dependent degree-day coefficient, a, overestimated snowmelt for both 2004 and 2005 (Figure 4). The mean absolute 7
10 errors (MAE) were 130 and 210 mm in 2004 and 2005, respectively. Relative differences between measured and calculated values at each site averaged 15 and 18%. On April 1, 2004 and 2005, the snow courses, which provide manual monthly snow measurements and are near the snow pillows, measured more SWE as compared to the snow pillows in those respective years. The MAE between the snow pillows and snow courses throughout the Tuolumne and Merced River basins were 120 and 240 mm in 2004 and 2005 (18 and 25% average relative differences). Therefore part of the error in our snowmelt calculations might be attributed to uncertainty in the snow-pillow measurements. The fractional SCA just prior to snowmelt in each 300-m elevation band was similar for 2004 and 2005, ranging from near 1.0 at the highest elevations bands to under 0.4 in the lowest band (Figure 5). The lowest band showed large fluctuations associated with snowfall and melt through the winter and spring, while at higher elevations these variations are smaller, but still indicate periods of snowmelt and snowfall. Near depletion of snow-covered area in 2004, identified here as SCA<0.1, occurred by March 15 in the lowest band, and occurred an average of 2-3 weeks later in each successively higher 300-m band. In 2005 SCA depletion occurred on average about five weeks later than in 2004 at the same elevation. There was little difference between the two basins in the timing of snow depletion, except in 2005 at the two highest elevation bands where the Merced melted out 3-5 weeks later than the Tuolumne. Interpolation of SWE from snow pillows and unconstrained by SCA gives estimates across elevation bands that during the accumulation period are fairly close together, with values becoming different during snowmelt (Figure 6). Note that melt at the snow pillows started in early March in both years, with further March accumulation in The spread in SWE values during snowmelt reflects greater melt at the lower snow pillows. 8
11 Masking interpolated SWE with the fractional SCA results in differences in SWE with elevation that reflect the SCA with elevation (Figure 7). Masking accounts for the interpolated SWE for areas with no snow, which is especially important at the lower to mid elevations. Note that while decreases in SCA indicate that snowmelt does occur early in the year at lower elevations, this is not reflected in the interpolated SWE of Figure 6 because there are no snow pillows at the lower elevations. Integrating the decline from peak SWE over time for each of the elevation bands in Figure 7 provides an estimate of the cumulative snowmelt volume (Figure 8). By this method, in 2004 snowmelt began in January, with more significant contributions beginning in February, and snowmelt was nearly complete by mid-may. Using the depletion approach to back-calculate snowmelt from the SCA patterns in Figure 5 gives the second estimate of snowmelt by elevation band (Figure 8). In 2004, this method shows no snowmelt in either the Tuolumne or Merced until March. That is, even though SCA decreased in the lower elevations early in the year, daily average temperature was below zero and thus snowmelt was calculated to be zero. Basin-total January-September snowmelt volumes estimated by the interpolation method were greater than that estimated by depletion: 42-44% higher in the Merced, 30-34% higher in the Tuolumne. The depletion method accounts for snow at elevations above that of all pillows, so it continues to estimate snowmelt above 2700 m in July and estimates a greater contribution to basin-wide melt from these higher elevations than does interpolation, whereas interpolated estimates of snowmelt contributions by elevation were approximately proportional to basin areas. Compared to the interpolation method, Table 2 and Figure 9 show that the depletion method shifts snowmelt contributions to higher elevations, especially in the drier year
12 Using values from the depletion method, elevations m contributed about 20% (24%) of the March 2004 snowmelt for the Tuolumne (Merced), declining to 19% (15%) and near zero in April (Table 3). Mid elevations ( m) contributed about 54% (70%) of March snowmelt, increasing to 61% (75%) in April, declining to 54% (66%) in May and dropping to about 25% (31%) by June and nearly disappearing by July. Elevations above 3000 m contributed only 22% (10%) of the monthly total in March, increasing to 43% (32%) in May and 95% in July. Trends were similar in March 2005, however, by April 1 the pattern of relative snowmelt contributions was about 1 month later than in 2004, with elevations m providing about 32-68% of the April snowmelt. Elevations m contributed 29% to 74% of the April, May and June snowmelt, declining to about 10% in August. At elevations above 3000 m snowmelt volume increased as the lower-elevation snow melted; and in July and August these upper elevations accounted for most of the monthly snowmelt (Table 3). For comparison with streamflow at Pohono Bridge in Yosemite National Park in the Merced basin (Figure 10), the per unit area potential snowmelt values were adjusted for vapor losses (sublimation + evapotranspiration) of 0.54 m yr -1 [Dettinger et al., 2004; Christensen et al., 2008]. Vapor losses were adjusted based on the distinct seasonality, nearly zero in November, increasing in February, and peaking in June. Figure 11 shows the comparison of the cumulative March-August streamflow averaged by depth over the basin (0.37 m in 2004 and 0.85 m in 2005) with the interpolation and depletion methods. The interpolation method overestimated the snowmelt by 9.4% in 2004 and 11% in In the calculation, snowmelt was nearly complete by May 1, 2004 and June 1, 2005, almost 2 months before the stream gauge at Pohono Bridge returned to baseflow levels. The depletion method overestimated the total snowmelt by 4.3% in 10
13 and underestimated by 8.8% in 2005, with the end of snowmelt delayed until mid-june in 2004 and July in Compared to traditional operational models, which use of regression models to predict seasonal streamflow volumes, the depletion method provides reasonable estimates, given that it is not calibrated against streamflow. The forecasts from the California Department of Water Resources in the same years over-predicted the April-July streamflow by 8.9% in 2004 and under-predicted it by 2.2% in It should be noted that CDWR had forecasted the total streamflow volume on April 1 for the 4 months (April-July); achieving this level of accuracy in forecast mode is impressive and should be acknowledged. The forecasts from 1990-present are available each year from February through May at Discussion Some possible reasons for differences in snowmelt estimates by the two methods include uncertainty in the degree-day factors, lack of representative ground-based measurements on which to base estimates, and uncertainties in SCA values for forested terrain. The degree-day factor (Figure 3), a, generally increases over the melt period because of seasonal increases in solar radiation [Granger and Male, 1978]. Previous studies have shown that in most circumstances the degree-day factor ranges between 1.8 and 6.5 mm deg -1 day -1 [Kuusisto, 1980; Kustas et al., 1994; Seidel and Martinec, 2004]. In this study, however, the degree-day factor showed significant day-to-day variations, likely a result of daily variations in the radiation flux. In 2004 snowmelt started on or before March 1; unseasonably high degreeday coefficients were calculated, as the average daily temperatures were low compared to the daily snowmelt rates. Since solar radiation, albedo, and turbulent-energy fluxes dominate snowmelt, day-to-day variations in the energy balance resulting from topography and landcover 11
14 can make the degree-day approach an unreliable estimate of daily melt, and can contribute to daily variability in the coefficient. However, over time the approach can provide a good measure of the average energy flux and melt, smoothing the day-to-day variations [Anderson, 1968; Ohmura, 2001]. Where sufficient data are available, the depletion method should be based on an energy-balance approach [Cline et al., 1998; Molotch et al., 2004], which can account for early season melt above treeline and topographic effects over the full snow season. The SCA values show clear increases with elevation during winter and steady decreases at all elevations during snowmelt despite small day-to-day variations resulting from patchy cloud cover, satellite view angle, occasional instrument noise from the MODIS sensor, and occasional late-season snowfalls that melt quickly (Figure 5). Although data from the sparse network of snow pillows contain less noise, they also contain less information about basin-scale distribution patterns of snow. Snow pillows were placed in elevations that are representative of the waterproducing regions of a watershed, in order to provide indices of streamflow for statistical watersupply outlooks [Farnes, 1967], not to calculate the water volume in a basin s snowpack. Locations are on flat or nearly flat ground, and thus do not represent the range of physiographic conditions in the surrounding catchments. Because of the multi-variant relationships between the geophysical variables and snowpack properties, interpolation schemes generally explain about only 50% of the SWE variability [Erickson et al., 2005; Molotch and Bales, 2006]. The snow pillows fail to sample the higher elevations (Figure 1). The MODIS images showed that the upper elevations were not snow free until August 2004 and never fully depleted of snow in 2005 (Figure 5), whereas the snow pillows showed complete snow depletion by June 2004 and July 2005 (Figure 6). Although the interpolation method uses observed SCA in the snowmelt estimations, all the input snow-pillow values are zero so the snow-volume estimate is zero. 12
15 Therefore, a spatial statistical model that uses the snow-pillow data will miss late-season melt unless one guesses at the SWE in the elevations above the snow-free pillows. Streamflow measured at Pohono Bridge on the Merced River in Yosemite Valley showed a snowmelt response beginning in early March in both years (Figure 10) with a return to baseflow by early July and August for 2004 and 2005, respectively. This suggests that snowmelt in both years began around the beginning of March. Comparing the snowmelt for just the March to September period, the interpolation and depletion methods give basin totals that are much closer, within about 5-25% (Figure 11; Tuolumne data not shown). Although snowmelt continues through July and August at the higher elevations, as shown by the depletion approach, the volume contribution was small. This lateseason contribution is, however, important for streamflow, soil moisture and ecosystem services. Because the MODIS fractional snow-cover product estimates the projected snow area, which is the snow not hidden by the forest canopy, during winter the snow cover is underestimated in the middle and lower elevation bands, but the effects of the trees diminish as snow melts. Discrepancies between SCA and canopy openings and the ability of the satellite to determine the amount of snow underneath the canopy [Liu et al., 2004], especially in highly forested areas, is problematic at elevations below 2100 m. At higher elevations, steep slopes keep SCA below 1.0. These issues of forest canopy and steep slopes apply the same bias to both the interpolation and depletion methods because SCA is used in both. To reduce this bias, other studies have used a vegetation gap fraction correction to account for more SCA than the satellite observes [Dressler et al., 2006; Durand et al., 2008]. If the fractional snow-covered area is greater than the gap fraction, then the satellite is likely measuring snow in the canopy, whereas fractional SCA less than the gap fraction may indicate that the canopy has hidden the snow from the sensor 13
16 (Figure 12). It should be noted that there is good correlation between canopy openings and fractional SCA as the R 2 was 0.91 for both the Tuolumne and Merced in 2004 and During the accumulation seasons of 2004 and 2005, vegetation may have caused underestimations of SWE by 18-30% at elevations below 2100 m, but when considering the total snowmelt volume, the underestimate represents only 2-4% of the total snowmelt volume. While the majority of the total snowmelt (50-60%) was derived from the m elevation range, after June the main source of snowmelt was above 3000 m. In 2004, less than 10% of the total snowmelt was derived from elevations of m, while elevations above 3000 m contributed more than 33% of the basin-wide snowmelt. During 2005, the m elevations contributed 11-17% of the overall snowmelt volumes, with elevations above 3000 m contributing a smaller fraction of the total, about 36% in the Tuolumne and 25% in the Merced. However, the actual volumes of snowmelt from the higher elevations were basically the same in both years, for the Merced and Tuolumne combined, equivalent to an area-weighted value of about 2.0 m in 2004 versus 2.2 m in 2005 (calculated from Table 1). Much of the difference between the total snowpack volume in the drier (2004) versus wetter (2005) year was in the m elevation range. For the two basins combined in this elevation range, this translated to an area-weighted average of about 0.60 m in 2004 versus 1.73 m in Note that 2005 was 47 degree days warmer than 2004 in January-February, so the difference in snow accumulation likely reflects a difference in precipitation, not a difference in temperature. Finally, it is important to consider the elevation dependence of the various methods. The depletion method gives a relatively steep elevation increase for snowmelt, about 1.5 m per 1000 m of elevation in the elevation range (Figure 13). The interpolation method gives SWE increases with elevation of about 0.1 m per 1000 m elevation (Figure 6); when masked 14
17 with SCA, this difference becomes about 0.2 m per 1000 m of elevation (Figure 7). In contrast, the commonly used PRISM data show essentially no increase in elevation for October to March precipitation (Figure 13). PRISM precipitation values are higher than snowmelt estimated by depletion. We have not used the snowmelt estimates to back calculate precipitation, as this would require an estimation of rainfall, which is currently not independently measured in these basins. Conclusions The MODIS fractional SCA product at 500-m resolution provides a consistent estimate of basin-scale snow-distribution patterns over mountainous terrain. Noise introduced by off-nadir viewing angles and cloud cover is generally small relative to trends, allowing for a spatially and temporally continuous snow product. Analysis of snowcover by 300-m elevation bands, which averages pixels with a range of slope, aspect and thus energy balance and melt rate, shows decreasing snow cover during the snowmelt season in each elevation band. The 2-3 weeks delay in melt in each successively higher elevation band represents a 2 o C decline in average temperature (based on average ground-surface lapse rate of 6.5 o C per 1000 m). Comparing the interpolated time series with that from the depletion calculation highlights the weakness in the ground-based SWE measurement network, particularly above 3000 m. Despite their smaller area, the high elevations will become more important as seasonal storage decreases at lower elevations with warming of the climate. For estimating a snowmelt time series, it is better to use the snow-pillow data to estimate degree-day factors for melt, i.e. basin-wide daily melt rates, than as absolute measures of SWE in a basin. The resulting back-calculated snowmelt time series using the depletion method provides basin-scale estimates of snowmelt that are consistent with spatially distributed observations of snow depletion. In the two basins 15
18 studied, about 55% of the snowmelt resulted from the m elevation bands in both drier (2004) and wetter (2005) years. The approximate one-third of the snowmelt from elevations above 3000 m, whose magnitude is not accurately estimated from ground measurements, contributes significantly to late-season runoff, baseflow and basin recharge. The 10-20% of seasonal snowmelt from below 2100-m elevation is most susceptible to falling as rain rather than snow and melting earlier in a warmer climate. These lower elevations also contribute to the inconsistency in January-February melt between the two methods, i.e. too little melt in the depletion method and too much melt in the interpolation approach. Though MODIS is a research instrument, and past its design lifetime, it is well suited to provide operational snow products with daily temporal resolution and large spatial coverage, and sets the stage for future snow products from VIIRS and GOES-R Advanced Baseline Imager. It is crucial that accurate fractional SCA products are available from these instruments, and the use of the MODSCAG product demonstrates the ability to estimate the spatial and temporal variability of snowmelt at research and operational scales for hydrologic modeling. Acknowledgements Support for this research was provided by NASACooperative AgreementNNG04GC52A and NSF Grant EAR We also acknowledge the contributions of X. Meng and P. Slaughter. Comments from the four anonymous reviewers greatly improved this manuscript
19 References Anderson, E. A. (1968), Development and testing of snow pack energy balance equations, Water Resour. Res., 4, 19-37, doi: /WR004i001p Anderson, E. A. (1973), National Weather Service River Forecast System snow accumulation and ablation model, Technical Memorandum NWS HYDRO-17, 217 pp, NOAA. Bales, R. C., N. P. Molotch, T. H. Painter, M. D. Dettinger, R. Rice, and J. Dozier (2006), Mountain hydrology of the western United States, Water Resour. Res., 42, W08432, doi: /2005WR Bales, R. C., K. A. Dressler, B. Imam, S. R. Fassnacht, and D. J. Lampkin (2008), Fractional snow cover in the Colorado and Rio Grande basins, , Water Resour. Res., 44, W01425, doi: /2006WR Bavera, D., and C. De Michele (2009), Snow water equivalent estimation in the Mallero basin using snow gauge data and MODIS images and fieldwork validation, Hydrol. Proc., 23, , doi: /hyp Carroll, S. S., T. R. Carroll, and R. W. Poston (1999), Spatial modeling and prediction of snow water equivalent using ground-based, airborne, and satellite snow data, J. Geophys. Res., 104, , doi: /1999JD Christensen, L., C. L. Tague, and J. S. Baron (2008), Spatial patterns of simulated transpiration response to climate variability in a snow dominated mountain ecosystem, Hydrol. Proc., 22, , doi: /hyp Cline, D. W., R. C. Bales, and J. Dozier (1998), Estimating the spatial distribution of snow in mountain basins using remote sensing and energy balance modeling, Water Resour. Res., 34, , doi: /97WR
20 Daly, S. F., R. Davis, E. Ochs, and T. Pangburn (2001), An approach to spatially distributed snow modelling of the Sacramento and San Joaquin basins, California, Hydrol. Proc., 14, , doi: / ( )14:18<3257::AID-HYP199>3.0.CO;2-Z. Dettinger, M. D., D. R. Cayan, M. K. Meyer, and A. E. Jeton (2004), Simulated hydrologic responses to climate variations and change in the Merced, Carson, and American River Basins, Sierra Nevada, California, , Climatic Change, 62, Dozier, J., and T. H. Painter (2004), Multispectral and hyperspectral remote sensing of alpine snow properties, Ann. Rev. Earth Planet. Sci., 32, , doi: /annurev.earth Dozier, J., T. H. Painter, K. Rittger, and J. E. Frew (2008), Time-space continuity of daily maps of fractional snow cover and albedo from MODIS, Adv. Water Res., 31, , doi: /j.advwatres Dozier, J., and J. Frew (2009), Computational provenance in hydrologic science: A snow mapping example, Phil. Trans. Royal Soc. A, 367, , doi: /rsta Dressler, K. A., G. H. Leavesley, R. C. Bales, and S. R. Fassnacht (2006), Evaluation of gridded snow water equivalent and satellite snow cover products for mountain basins in a hydrologic model, Hydrol. Proc., 20, , doi: /hyp Durand, M., N. P. Molotch, and S. A. Margulis (2008), A Bayesian approach to snow water equivalent reconstruction, J. Geophys. Res., 113, D20117, doi: /2008JD Erickson, T. A., M. W. Williams, and A. Winstral (2005), Persistence of topographic controls on the spatial distribution of snow in rugged mountain terrain, Colorado, United States, Water Resour. Res., 41, W04014, doi: /2003WR
21 Farnes, E. F. (1967), Criteria for determining mountain snow-pillow sites, Proc. Western Snow Conf., 35, Farr, T. G., P. A. Rosen, E. Caro, R. Crippen, R. Duren, S. Hensley, M. Kobrick, M. Paller, E. Rodriguez, L. Roth, D. Seal, S. Shaffer, J. Shimada, J. Umland, M. Werner, M. Oskin, D. Burbank, and D. Alsdorf (2007), The Shuttle Radar Topography Mission, Rev. Geophys., 45, RG2004, doi: /2005RG Fassnacht, S. R., K. A. Dressler, and R. C. Bales (2003), Snow water equivalent interpolation for the Colorado River Basin from snow telemetry (SNOTEL) data, Water Resour. Res., 39, 1208, doi: /2002WR Granger, R. J., and D. H. Male (1978), Melting of a prairie snowpack, J. Appl. Meteorol., 17, , doi: / (1978)017<1833:MOAPS>2.0.CO;2. Harshburger, B. J., K. S. Humes, V. P. Walden, T. R. Blandford, B. C. Moore, and R. J. Dezzani (2010), Spatial interpolation of snow water equivalency using surface observations and remotely sensed images of snow-covered area, Hydrol. Proc., 24, , doi: /hyp Jeton, A. E., M. D. Dettinger, and J. L. Smith (1996), Potential effects of elimate change on streamflow, eastern and western slopes of the Sierra Nevada, California and Nevada, Water Resources Investigations Rep , U. S. Geological Survey. Kustas, W., A. Rango, and R. Uijlenhoet (1994), A simple energy budget algorithm for the snowmelt runoff model, Water Resour. Res., 30, , doi: /94WR Kuusisto, E. (1980), On the values and variability of degree-day melting factor in Finland, Nordic Hydrol., 11, , doi: /nh
22 Liston, G. E. (1999), Interrelationships among snow distribution, snowmelt, and snow cover depletion: Implications for atmospheric, hydrologic, and ecologic modeling, J. Appl. Meteorol., 44, doi: / (1999)038<1474:IASDSA>2.0.CO;2. Liu, J., R. A. Melloh, C. E. Woodcock, R. E. Davis, and E. S. Ochs (2004), The effect of viewing geometry and topography on viewable gap fractions through forest canopies, Hydrol. Proc., 18, , doi: /hyp Lundquist, J. D., and D. R. Cayan (2007), Surface temperature patterns in complex terrain: Daily variations and long-term change in the central Sierra Nevada, California, J. Geophys. Res., 112, D11124, doi: /2006JD Lundquist, J. D., N. Pepin, and C. Rochford (2008), Automated algorithm for mapping regions of cold-air pooling in complex terrain, J. Geophys. Res., 113, D22107, doi: /2008JD Martinec, J., and A. Rango (1981), Areal distribution of snow water equivalent evaluated by snow cover monitoring, Water Resour. Res., 17, , doi: /WR017i005p Molotch, N. P., T. H. Painter, R. C. Bales, and J. Dozier (2004), Incorporating remotely sensed snow albedo into spatially distributed snowmelt modeling, Geophys. Res. Lett., 31, L03501, doi: /2003GL Molotch, N. P., and R. C. Bales (2005), Scaling snow observations from the point to the grid element: implications for observation network design, Water Resour. Res., 41, W11421, doi: /2005WR Molotch, N. P., and R. C. Bales (2006), SNOTEL representativeness in the Rio Grande headwaters on the basis of physiographics and remotely sensed snow cover persistence, Hydrol. Proc., 20, , doi: /hyp
23 Ohmura, A. (2001), Physical basis for the temperature-based melt-index method, J. Appl. Meteorol., 40, , doi: / (2001)040<0753:PBFTTB>2.0.CO;2. Painter, T. H., K. Rittger, C. McKenzie, R. E. Davis, and J. Dozier (2009), Retrieval of subpixel snow-covered area, grain size, and albedo from MODIS, Remote Sens. Environ., 113, , doi: /j.rse Rango, A., and J. Martinec (1995), Revisiting the degree-day method for snowmelt computations, Water Resour. Bull., 31, , doi: /j tb03392.x. Rice, R., and R. C. Bales (2010), Embedded sensor network design for snowcover measurements around snow-pillow and snow-course sites in the Sierra Nevada of California, Water Resour. Res., 46, W03537, doi: /2008WR Seidel, K., and J. Martinec (2004), Remote Sensing in Snow Hydrology: Runoff Modeling, Effect of Climate Change, 172 pp., Springer-Verlag, Berlin. Shamir, E., and K. P. Georgakakos (2006), Distributed snow accumulation and ablation modeling in the American River basin, Adv. Water Res., 29, , doi: /j.advwatres Slack, J. R., and J. M. Landwehr (1992), Hydro-climatic data detwork (HCDN): A U.S. Geological Survey streamflow data set of the United States for the study of climate variations, , Open-File Rep , U. S. Geological Survey
24 488 Tables Table 1. Areas of Tuolumne and Merced River basins above 1500 m, and snowmelt estimates a Merced Tuolumne elevation band, m Avg snowmelt, m Avg snowmelt, m area, km 2 area, km 2 midpoint range a Snowmelt values based on end-of-season volumes from Figure 8 by the depletion method. 22
25 Table 2. Snowmelt contributions in elevation bands by interpolation and depletion methods Elevation band, m Catchment area, % Contributions by interpolation (depletion) method, % Tuolumne Merced Tuolumne Merced (8) 28 (17) 28 (7) 31 (11) (43) 55 (47) 62 (57) 61 (64) > (49) 17 (36) 10 (36) 8 (25) 23
26 492 Table 3. Monthly snowmelt contribution by elevation band, calculated by depletion method Elevation band, m Tuolumne Contribution for 2004 (2005), % March April May June July August 24 (54) 19 (68) 3 (16) 0 (3) 0 (0) 0 (0) 54 (42) 61 (29) 54 (60) 25 (59) 5 (25) 5 (25) > (4) 20 (3) 43 (24) 75 (38) 95 (75) 95 (75) Merced (27) 15 (32) 2 (11) 0 (2) 0 (0) 0 (0) 70 (70) 75 (65) 66 (74) 31 (69) 5 (29) 5 (29) > (3) 10 (3) 32 (15) 69 (29) 95 (71) 95 (71)
27 Figure captions Figure 1. Upper Tuolumne and Merced River basins with the MODIS fractional snow-covered area for April 1, 2005 at 500-m resolution. Elevations of the 9 snow sensors are also shown. Figure 2. Degree-day coefficients for sites in Figure 1. Figure 3. Elevation dependent degree-day coefficients across 4 of 8 elevation bands in 2004 and Values are shown for period of snowmelt at each elevation. Figure 4. Observed snowmelt as measured at the snow-pillow and the back-calculated snowmelt in 2004 and 2005 using the elevation-dependent degree-day coefficients, a. Figure 5. Daily snow-covered area in the eight 300-m elevation bands. Figure 6. Daily interpolated snow water equivalent in 8 elevation bands in the Tuolumne and Merced River basins for 2004 and Figure 7. Interpolated snow water equivalent masked with fractional snow-covered area in the eight elevation bands. Figure 8. Cumulative snowmelt for January-September in 2004 and 2005; solid lines represent the interpolation method and dotted lines the depletion method. Note that each curve is cumulative over both time and elevation. Figure 9. Distribution of catchment area and elevation contributions of snowmelt for Tuolumne and Merced basins for January- September, calculated by both the interpolation and depletion methods. Figure 10. Streamflow for Water Years 2004 and 2005 at Pohono Bridge on the Merced River in Yosemite National Park measured by the United States Geological Survey (site ). Figure 11. Cumulative streamflow depth at Pohono Bridge on the Merced River in Yosemite National Park and the cumulative per unit area potential snowmelt depth from the interpolation 25
28 and depletion methods minus sublimation and evapotranspiration (ET) for March September 2004 and Figure 12. Mean snow-covered area across the eight elevation bands during the accumulation period plotted with the mean canopy openings across for each elevation band, averaged across all pixels in that elevation band. The error bars show the standard deviation from the mean for both SCA and canopy openings. Vegetation data from J. Van Wagtendonk, USGS, based on aerial photography gridded at 10m. Figure 13. Comparison of the PRISM 4 km precipitation (Oct-Mar) across the 8 elevation bands in the Merced and Tuolumne River basins versus the estimated depth of SWE derived from the depletion method. 26
29 Figure 1. Upper Tuolumne and Merced River basins with the MODIS fractional snow-covered area for April 1, 2005 at 500-m resolution. Elevations of the 9 snow sensors are also shown. 27
30 531 Figure 2. Degree-day coefficients for sites in Figure 1 28
31 Figure 3. Elevation dependent degree-day coefficients across 4 of 8 elevation bands in 2004 and Values are shown for period of snowmelt at each elevation
32 Figure 4. Observed snowmelt as measured at the snow-pillow and the back-calculated snowmelt in 2004 and 2005 using the elevation dependent degree-day coefficients, a. 30
33 533 Figure 5. Daily snow-covered area in the eight 300-m elevation bands. 31
34 Figure 6. Daily interpolated SWE in 8 elevation bands in the Tuolumne and Merced River basins for 2004 and
35 Figure 7. Interpolated snow water equivalent masked with fractional snow-covered area in the eight elevation bands. 33
36 Figure 8. Cumulative snowmelt for January-September in 2004 and 2005; solid lines represent the interpolation method and dotted lines the depletion method. Note that each curve is cumulative over both time and elevation. 34
37 Figure 9. Distribution of catchment area and elevation contributions of snowmelt for Tuolumne and Merced basins for January- September, calculated by both the interpolation and depletion methods
38 Figure 10. Streamflow for Water Years 2004 and 2005 at Pohono Bridge on the Merced River in Yosemite National Park measured by the United States Geological Survey (site )
39 Figure 11. Cumulative streamflow depth at Pohono Bridge on the Merced River in Yosemite National Park and the cumulative per unit area potential snowmelt depth from the interpolation and depletion methods minus sublimation and evapotranspiration (ET) for March September 2004 and
40 Figure 12. Mean snow-covered area across the eight elevation bands during the accumulation period plotted with the mean canopy openings across for each elevation band, averaged across all pixels in that elevation band. The error bars show the standard deviation from the mean for both SCA and canopy openings. Vegetation data from J. Van Wagtendonk, USGS, based on aerial photography gridded at 10m
Real Time Snow Water Equivalent (SWE) Simulation March 29, 2015 Sierra Nevada Mountains, California
Real Time Snow Water Equivalent (SWE) Simulation March 29, 2015 Sierra Nevada Mountains, California Introduction We have developed a real-time SWE estimation scheme based on historical SWE reconstructions
More informationSnowcover along elevation gradients in the Upper Merced and Tuolumne River basin of the Sierra Nevada of California from MODIS and blended ground data
Snowcover along elevation gradients in the Upper Merced and Tuolumne River basin of the Sierra Nevada of California from MODIS and blended ground data Robert Rice a, Roger Bales a, Thomas H. Painter b,
More informationReal Time Snow Water Equivalent (SWE) Simulation February 5, 2014 Sierra Nevada Mountains, California
Real Time Snow Water Equivalent (SWE) Simulation February 5, 2014 Sierra Nevada Mountains, California Introduction We have developed a real-time SWE estimation scheme based on historical SWE reconstructions
More informationReal Time (RT) Snow Water Equivalent (SWE) Simulation May 11, 2014 Sierra Nevada Mountains, California
Real Time (RT) Snow Water Equivalent (SWE) Simulation May 11, 2014 Sierra Nevada Mountains, California Abstract On May 11 th, percent of average SWE values for this date have shown an increase (see note
More informationFigure 2. April 8, 2011 spatial SWE shown for the Weber, UT watershed.
Real-Time Snow Water Equivalent (SWE) Simulation April 8, 2011 Intermountain West Noah Molotch, Dominik Schneider and Leanne Lestak University of Colorado, Institute of Arctic and Alpine Research Noah.Molotch@colorado.edu,
More informationSurface & subsurface processes in mountain environments
Surface & subsurface processes in mountain environments evapotranspiration snowmelt precipitation infiltration Roger Bales Martha Conklin Robert Rice Fengjing Liu Peter Kirchner runoff sublimation ground
More informationQuenching the Valley s thirst: The connection between Sierra Nevada snowpack & regional water supply
Quenching the Valley s thirst: The connection between Sierra Nevada snowpack & regional water supply Roger Bales, UC Merced Snow conditions Snow & climate change Research directions Sierra Nevada snow
More informationMODIS True Color Terra image from from NASA WorldView (worldview.earthdata.nasa.gov)
Estimating the snowpack with remote sensing MODIS True Color Terra image from 2016-09-07 from NASA WorldView (worldview.earthdata.nasa.gov) Edward (Ned) Bair, Ph.D. Earth Research Institute, University
More informationSierra Nevada Hydroclimatology: An Experimental Prediction of Maximum Daily Snowmelt Discharge in 2005
Sierra Nevada Hydroclimatology: An Experimental Prediction of Maximum Daily Snowmelt Discharge in 2005 David Peterson, Jessica Lundquist, Iris Stewart, Noah Knowles, Madeline Solomon and Stephen Hager
More informationSNOW WATER EQUIVALENT ESTIMATES IN THE HINDU KUSH AND THE SIERRA NEVADA USING PASSIVE MICROWAVE AND RECONSTRUCTION ABSTRACT
SNOW WATER EQUIVALENT ESTIMATES IN THE HINDU KUSH AND THE SIERRA NEVADA USING PASSIVE MICROWAVE AND RECONSTRUCTION Edward H. Bair 1,2, Jeff Dozier 3, Karl Rittger 4, Carrie Voyuvich 1, and Robert E. Davis
More informationWater information system advances American River basin. Roger Bales, Martha Conklin, Steve Glaser, Bob Rice & collaborators UC: SNRI & CITRIS
Water information system advances American River basin Roger Bales, Martha Conklin, Steve Glaser, Bob Rice & collaborators UC: SNRI & CITRIS Opportunities Unprecedented level of information from low-cost
More informationABSTRACT INTRODUCTION
INTER- AND INTRA-ANNUAL VARIABILITY IN SNOW ALBEDO, GRAIN SIZE, AND SNOW COVERED AREA FROM THE AIRBORNE SNOW OBSERVATORY DURING LOW SNOW YEARS (2013-2015), TUOLUMNE RIVER BASIN, CA S. McKenzie Skiles 1,
More informationHyMet Company. Streamflow and Energy Generation Forecasting Model Columbia River Basin
HyMet Company Streamflow and Energy Generation Forecasting Model Columbia River Basin HyMet Inc. Courthouse Square 19001 Vashon Hwy SW Suite 201 Vashon Island, WA 98070 Phone: 206-463-1610 Columbia River
More informationESTIMATING SNOWMELT CONTRIBUTION FROM THE GANGOTRI GLACIER CATCHMENT INTO THE BHAGIRATHI RIVER, INDIA ABSTRACT INTRODUCTION
ESTIMATING SNOWMELT CONTRIBUTION FROM THE GANGOTRI GLACIER CATCHMENT INTO THE BHAGIRATHI RIVER, INDIA Rodney M. Chai 1, Leigh A. Stearns 2, C. J. van der Veen 1 ABSTRACT The Bhagirathi River emerges from
More informationPrediction of Snow Water Equivalent in the Snake River Basin
Hobbs et al. Seasonal Forecasting 1 Jon Hobbs Steve Guimond Nate Snook Meteorology 455 Seasonal Forecasting Prediction of Snow Water Equivalent in the Snake River Basin Abstract Mountainous regions of
More informationUse of Snow Data from Remote Sensing in Operational Streamflow Prediction
Use of Snow Data from Remote Sensing in Operational Streamflow Prediction Stacie Bender 1, Thomas H. Painter 2, Paul Miller 1, Michelle Stokes 1 1 NOAA/National Wear Service Colorado Basin River Forecast
More informationEvaluation of gridded snow water equivalent and satellite snow cover products for mountain basins in a hydrologic model
Evaluation of gridded snow water equivalent and satellite snow cover products for mountain basins in a hydrologic model K.A. Dressler 1, G. H. Leavesley 2, R.C. Bales 3, S.R. Fassnacht 4 1 Penn State Institutes
More informationPreliminary Runoff Outlook February 2018
Preliminary Runoff Outlook February 2018 Prepared by: Flow Forecasting & Operations Planning Water Security Agency General Overview The Water Security Agency (WSA) is preparing for 2018 spring runoff including
More informationStreamflow estimation from hydrologic model updates of remotely sensed snow information in snowmelt dominated basins
Streamflow estimation from hydrologic model updates of remotely sensed snow information in snowmelt dominated basins K.A. Dressler 1, G. H. Leavesley 2, R.C. Bales 3, S.R. Fassnacht 4 1 Department of Hydrology
More informationNIDIS Intermountain West Regional Drought Early Warning System February 7, 2017
NIDIS Drought and Water Assessment NIDIS Intermountain West Regional Drought Early Warning System February 7, 2017 Precipitation The images above use daily precipitation statistics from NWS COOP, CoCoRaHS,
More informationClimate Change and Water Supplies in the West. Michael Dettinger, USGS
Climate Change and Water Supplies in the West Michael Dettinger, USGS During the past 30 years--and in the near-term future--natural climate variability has been augmented by warming trends associated
More informationDirection and range of change expected in the future
Direction and range of Air Temperature Over the past 30 years, air Across the greater PNW and temperature has been Columbia Basin, an ensemble increasing an average of forecast from ten of the best 0.13
More informationCentral Asia Regional Flash Flood Guidance System 4-6 October Hydrologic Research Center A Nonprofit, Public-Benefit Corporation
http://www.hrcwater.org Central Asia Regional Flash Flood Guidance System 4-6 October 2016 Hydrologic Research Center A Nonprofit, Public-Benefit Corporation FFGS Snow Components Snow Accumulation and
More informationRELATIVE IMPORTANCE OF GLACIER CONTRIBUTIONS TO STREAMFLOW IN A CHANGING CLIMATE
Proceedings of the Second IASTED International Conference WATER RESOURCE MANAGEMENT August 20-22, 2007, Honolulu, Hawaii, USA ISGN Hardcopy: 978-0-88986-679-9 CD: 978-0-88-986-680-5 RELATIVE IMPORTANCE
More informationEmbedded sensor network design for snow cover measurements around snow pillow and snow course sites in the Sierra Nevada of California
Click Here for Full Article WATER RESOURCES RESEARCH, VOL. 46,, doi:10.1029/2008wr007318, 2010 Embedded sensor network design for snow cover measurements around snow pillow and snow course sites in the
More informationJanuary 2011 Calendar Year Runoff Forecast
January 2011 Calendar Year Runoff Forecast 2010 Runoff Year Calendar Year 2010 was the third highest year of runoff in the Missouri River Basin above Sioux City with 38.8 MAF, behind 1978 and 1997 which
More informationColorado s 2003 Moisture Outlook
Colorado s 2003 Moisture Outlook Nolan Doesken and Roger Pielke, Sr. Colorado Climate Center Prepared by Tara Green and Odie Bliss http://climate.atmos.colostate.edu How we got into this drought! Fort
More informationClimate Change in Colorado: Recent Trends, Future Projections and Impacts An Update to the Executive Summary of the 2014 Report
Climate Change in Colorado: Recent Trends, Future Projections and Impacts An Update to the Executive Summary of the 2014 Report Jeff Lukas, Western Water Assessment, University of Colorado Boulder - Lukas@colorado.edu
More informationA Review of the 2007 Water Year in Colorado
A Review of the 2007 Water Year in Colorado Nolan Doesken Colorado Climate Center, CSU Mike Gillespie Snow Survey Division, USDA, NRCS Presented at the 28 th Annual AGU Hydrology Days, March 26, 2008,
More informationScaling snow observations from the point to the grid element: Implications for observation network design
WATER RESOURCES RESEARCH, VOL. 41, W11421, doi:10.1029/2005wr004229, 2005 Scaling snow observations from the point to the grid element: Implications for observation network design Noah P. Molotch Cooperative
More informationUsing MODIS imagery to validate the spatial representation of snow cover extent obtained from SWAT in a data-scarce Chilean Andean watershed
Using MODIS imagery to validate the spatial representation of snow cover extent obtained from SWAT in a data-scarce Chilean Andean watershed Alejandra Stehr 1, Oscar Link 2, Mauricio Aguayo 1 1 Centro
More informationHydrologic Analysis and Process-Based Modeling for the Upper Cache La Poudre Basin
Hydrologic Analysis and Process-Based Modeling for the Upper Cache La Poudre Basin Stephanie K. Kampf Eric E. Richer October 2010 Completion Report No. 222 Acknowledgements The authors thank the Colorado
More informationLake Tahoe Watershed Model. Lessons Learned through the Model Development Process
Lake Tahoe Watershed Model Lessons Learned through the Model Development Process Presentation Outline Discussion of Project Objectives Model Configuration/Special Considerations Data and Research Integration
More informationUpper Missouri River Basin February 2018 Calendar Year Runoff Forecast February 6, 2018
Upper Missouri River Basin February 2018 Calendar Year Runoff Forecast February 6, 2018 Calendar Year Runoff Forecast Explanation and Purpose of Forecast U.S. Army Corps of Engineers, Northwestern Division
More informationStreamflow estimation from hydrologic model updates of remotely sensed snow information in snowmelt dominated basins
Streamflow estimation from hydrologic model updates of remotely sensed snow information in snowmelt dominated basins K.A. Dressler 1, G. H. Leavesley 2, R.C. Bales 3, S.R. Fassnacht 4 1 Penn State Institutes
More informationUpper Missouri River Basin December 2017 Calendar Year Runoff Forecast December 5, 2017
Upper Missouri River Basin December 2017 Calendar Year Runoff Forecast December 5, 2017 Calendar Year Runoff Forecast Explanation and Purpose of Forecast U.S. Army Corps of Engineers, Northwestern Division
More informationGreat Lakes Update. Volume 199: 2017 Annual Summary. Background
Great Lakes Update Volume 199: 2017 Annual Summary Background The U.S. Army Corps of Engineers (USACE) tracks and forecasts the water levels of each of the Great Lakes. This report is primarily focused
More informationDEVELOPMENT OF A LARGE-SCALE HYDROLOGIC PREDICTION SYSTEM
JP3.18 DEVELOPMENT OF A LARGE-SCALE HYDROLOGIC PREDICTION SYSTEM Ji Chen and John Roads University of California, San Diego, California ABSTRACT The Scripps ECPC (Experimental Climate Prediction Center)
More informationANALYZING THE IMPACT OF CLIMATE CHANGE ON MONTHLY RIVER FLOWS IN CALIFORNIA S SIERRA NEVADA AND SOUTHERN CASCADE MOUNTAIN RANGES. Gary J.
ANALYZING THE IMPACT OF CLIMATE CHANGE ON MONTHLY RIVER FLOWS IN CALIFORNIA S SIERRA NEVADA AND SOUTHERN CASCADE MOUNTAIN RANGES Gary J. Freeman 1 ABSTRACT The impact of climate change on monthly river
More informationINTERPOLATING SURFACE AIR TEMPERATURE FOR USE IN SEMI-DISTRIBUTED SNOWMELT RUNOFF MODELS
INTERPOLATING SURFACE AIR TEMPERATURE FOR USE IN SEMI-DISTRIBUTED SNOWMELT RUNOFF MODELS Best Poster Award, WSC 2005 - T. R. Blandford, B. J. Harshburger, K. S. Humes, B. C. Moore, V. P. Walden 1 ABSTRACT
More informationThe Documentation of Extreme Hydrometeorlogical Events: Two Case Studies in Utah, Water Year 2005
The Documentation of Extreme Hydrometeorlogical Events: Two Case Studies in Utah, Water Year 2005 Tim Bardsley1*, Mark Losleben2, Randy Julander1 1. USDA, NRCS, Snow Survey Program, Salt Lake City, Utah.
More informationOregon Water Conditions Report April 17, 2017
Oregon Water Conditions Report April 17, 2017 Mountain snowpack continues to maintain significant levels for mid-april. By late March, statewide snowpack had declined to 118 percent of normal after starting
More informationJanuary 25, Summary
January 25, 2013 Summary Precipitation since the December 17, 2012, Drought Update has been slightly below average in parts of central and northern Illinois and above average in southern Illinois. Soil
More informationNIDIS Intermountain West Drought Early Warning System February 6, 2018
NIDIS Intermountain West Drought Early Warning System February 6, 2018 Precipitation The images above use daily precipitation statistics from NWS COOP, CoCoRaHS, and CoAgMet stations. From top to bottom,
More informationNIDIS Intermountain West Drought Early Warning System April 18, 2017
1 of 11 4/18/2017 3:42 PM Precipitation NIDIS Intermountain West Drought Early Warning System April 18, 2017 The images above use daily precipitation statistics from NWS COOP, CoCoRaHS, and CoAgMet stations.
More informationNIDIS Intermountain West Drought Early Warning System February 12, 2019
NIDIS Intermountain West Drought Early Warning System February 12, 2019 Precipitation The images above use daily precipitation statistics from NWS COOP, CoCoRaHS, and CoAgMet stations. From top to bottom,
More informationCanadian Prairie Snow Cover Variability
Canadian Prairie Snow Cover Variability Chris Derksen, Ross Brown, Murray MacKay, Anne Walker Climate Research Division Environment Canada Ongoing Activities: Snow Cover Variability and Links to Atmospheric
More informationOregon Water Conditions Report May 1, 2017
Oregon Water Conditions Report May 1, 2017 Mountain snowpack in the higher elevations has continued to increase over the last two weeks. Statewide, most low and mid elevation snow has melted so the basin
More informationAn investigation of sampling efficiency using historical data. Patrick Didier Advisor: Justine Blanford
An investigation of sampling efficiency using historical data Patrick Didier Advisor: Justine Blanford Overview Introduction to Airborne Snow Survey program Goals of the Capstone Project Methods and Expected
More informationSouris River Basin Spring Runoff Outlook As of March 15, 2018
Souris River Basin Spring Runoff Outlook As of March 15, 2018 Prepared by: Flow Forecasting & Operations Planning Water Security Agency Basin Conditions Summer rainfall in 2017 in the Saskatchewan portion
More informationNIDIS Drought and Water Assessment
NIDIS Drought and Water Assessment PRECIPITATION The images above use daily precipitation statistics from NWS COOP, CoCoRaHS, and CoAgMet stations. From top to bottom, and left to right: most recent 7-days
More informationOperational Perspectives on Hydrologic Model Data Assimilation
Operational Perspectives on Hydrologic Model Data Assimilation Rob Hartman Hydrologist in Charge NOAA / National Weather Service California-Nevada River Forecast Center Sacramento, CA USA Outline Operational
More informationNovel Snotel Data Uses: Detecting Change in Snowpack Development Controls, and Remote Basin Snow Depth Modeling
Novel Snotel Data Uses: Detecting Change in Snowpack Development Controls, and Remote Basin Snow Depth Modeling OVERVIEW Mark Losleben and Tyler Erickson INSTAAR, University of Colorado Mountain Research
More informationMapping the extent of temperature-sensitive snowcover and the relative frequency of warm winters in the western US
Mapping the extent of temperature-sensitive snowcover and the relative frequency of warm winters in the western US Anne Nolin Department of Geosciences Oregon State University Acknowledgements Chris Daly,
More informationNIDIS Intermountain West Drought Early Warning System November 21, 2017
NIDIS Drought and Water Assessment NIDIS Intermountain West Drought Early Warning System November 21, 2017 Precipitation The images above use daily precipitation statistics from NWS COOP, CoCoRaHS, and
More informationTHE ROLE OF MICROSTRUCTURE IN FORWARD MODELING AND DATA ASSIMILATION SCHEMES: A CASE STUDY IN THE KERN RIVER, SIERRA NEVADA, USA
MICHAEL DURAND (DURAND.8@OSU.EDU), DONGYUE LI, STEVE MARGULIS Photo: Danielle Perrot THE ROLE OF MICROSTRUCTURE IN FORWARD MODELING AND DATA ASSIMILATION SCHEMES: A CASE STUDY IN THE KERN RIVER, SIERRA
More informationStorm and Runoff Calculation Standard Review Snowmelt and Climate Change
Storm and Runoff Calculation Standard Review Snowmelt and Climate Change Presented by Don Moss, M.Eng., P.Eng. and Jim Hartman, P.Eng. Greenland International Consulting Ltd. Map from Google Maps TOBM
More informationModeling of peak inflow dates for a snowmelt dominated basin Evan Heisman. CVEN 6833: Advanced Data Analysis Fall 2012 Prof. Balaji Rajagopalan
Modeling of peak inflow dates for a snowmelt dominated basin Evan Heisman CVEN 6833: Advanced Data Analysis Fall 2012 Prof. Balaji Rajagopalan The Dworshak reservoir, a project operated by the Army Corps
More informationPresentation Overview. Southwestern Climate: Past, present and future. Global Energy Balance. What is climate?
Southwestern Climate: Past, present and future Mike Crimmins Climate Science Extension Specialist Dept. of Soil, Water, & Env. Science & Arizona Cooperative Extension The University of Arizona Presentation
More informationUpper Missouri River Basin January 2018 Calendar Year Runoff Forecast January 3, 2018
Calendar Year Runoff Forecast Upper Missouri River Basin January 2018 Calendar Year Runoff Forecast January 3, 2018 Explanation and Purpose of Forecast U.S. Army Corps of Engineers, Northwestern Division
More information2003 Moisture Outlook
2003 Moisture Outlook Nolan Doesken and Roger Pielke, Sr. Colorado Climate Center Prepared by Tara Green and Odie Bliss http://climate.atmos.colostate.edu Through 1999 Through 1999 Fort Collins Total Water
More informationWater, geology, climate & soils: Sierra Nevada, Merced River & Central Valley connections
Water, geology, climate & soils: Sierra Nevada, Merced River & Central Valley connections Geography & geology Merced River basin Climate & hydrology Soils & soil formation Roger Bales UC Merced Water,
More informationNIDIS Intermountain West Drought Early Warning System October 17, 2017
NIDIS Drought and Water Assessment NIDIS Intermountain West Drought Early Warning System October 17, 2017 Precipitation The images above use daily precipitation statistics from NWS COOP, CoCoRaHS, and
More informationDrought and Future Water for Southern New Mexico
Drought and Future Water for Southern New Mexico David S. Gutzler University of New Mexico gutzler@unm.edu * The big drought of 2018 * Longer term challenges for water supply * Forecasting streamflow Elephant
More informationNIDIS Intermountain West Drought Early Warning System May 23, 2017
NIDIS Drought and Water Assessment NIDIS Intermountain West Drought Early Warning System May 23, 2017 Precipitation The images above use daily precipitation statistics from NWS COOP, CoCoRaHS, and CoAgMet
More informationMapping Temperature across Complex Terrain
Mapping Temperature across Complex Terrain Jessica Lundquist 1, Nick Pepin 2, Phil Mote 3 1 Assistant Professor, Civil and Environmental Engineering, University of Washington 2 Lecturer, Department of
More informationNIDIS Intermountain West Drought Early Warning System January 15, 2019
NIDIS Drought and Water Assessment NIDIS Intermountain West Drought Early Warning System January 15, 2019 Precipitation The images above use daily precipitation statistics from NWS COOP, CoCoRaHS, and
More informationFlood Risk Assessment
Flood Risk Assessment February 14, 2008 Larry Schick Army Corps of Engineers Seattle District Meteorologist General Assessment As promised, La Nina caused an active winter with above to much above normal
More informationThe Importance of Snowmelt Runoff Modeling for Sustainable Development and Disaster Prevention
The Importance of Snowmelt Runoff Modeling for Sustainable Development and Disaster Prevention Muzafar Malikov Space Research Centre Academy of Sciences Republic of Uzbekistan Water H 2 O Gas - Water Vapor
More informationNIDIS Intermountain West Drought Early Warning System December 11, 2018
NIDIS Drought and Water Assessment NIDIS Intermountain West Drought Early Warning System December 11, 2018 Precipitation The images above use daily precipitation statistics from NWS COOP, CoCoRaHS, and
More informationS. R. FASSNACHT. Watershed Science, Colorado State University, Fort Collins, Colorado Z.-L. YANG
298 J O U R N A L O F H Y D R O M E T E O R O L O G Y VOLUME 7 Effects of Averaging and Separating Soil Moisture and Temperature in the Presence of Snow Cover in a SVAT and Hydrological Model for a Southern
More informationCARFFG System Development and Theoretical Background
CARFFG Steering Committee Meeting 15 SEPTEMBER 2015 Astana, KAZAKHSTAN CARFFG System Development and Theoretical Background Theresa M. Modrick, PhD Hydrologic Research Center Key Technical Components -
More informationP1.34 MULTISEASONALVALIDATION OF GOES-BASED INSOLATION ESTIMATES. Jason A. Otkin*, Martha C. Anderson*, and John R. Mecikalski #
P1.34 MULTISEASONALVALIDATION OF GOES-BASED INSOLATION ESTIMATES Jason A. Otkin*, Martha C. Anderson*, and John R. Mecikalski # *Cooperative Institute for Meteorological Satellite Studies, University of
More informationModelling snow accumulation and snow melt in a continuous hydrological model for real-time flood forecasting
IOP Conference Series: Earth and Environmental Science Modelling snow accumulation and snow melt in a continuous hydrological model for real-time flood forecasting To cite this article: Ph Stanzel et al
More informationForecast Challenges for the Colorado Basin River Forecast Center
Forecast Challenges for the Colorado Basin River Forecast Center John Lhotak, Development and Operations Hydrologist Improving Sub-Seasonal to Seasonal (S2S) Workshop May 16th, 2018 1 River Forecast Centers
More informationPRECIPITATION. Last Week Precipitation:
12/3/2014 NIDIS Drought and Water Assessment PRECIPITATION The images above use daily precipitation statistics from NWS COOP, CoCoRaHS, and CoAgMet stations. From top to bottom, and left to right: most
More informationAppendix A Calibration Memos
Appendix A Calibration Memos University of Washington Department of Civil and Environmental Engineering TO: Joe Dvorak FROM: Margaret Hahn and Richard Palmer RE: DVSVM Calibration DATE: June 28, 21 This
More informationEVALUATION AND MONITORING OF SNOWCOVER WATER RESOURCES IN CARPATHIAN BASINS USING GEOGRAPHIC INFORMATION AND SATELLITE DATA
EVALUATION AND MONITORING OF SNOWCOVER WATER RESOURCES IN CARPATHIAN BASINS USING GEOGRAPHIC INFORMATION AND SATELLITE DATA Gheorghe Stancalie, Simona Catana, Anisoara Iordache National Institute of Meteorology
More informationAnalysis of meteorological measurements made over three rainy seasons in Sinazongwe District, Zambia.
Analysis of meteorological measurements made over three rainy seasons in Sinazongwe District, Zambia. 1 Hiromitsu Kanno, 2 Hiroyuki Shimono, 3 Takeshi Sakurai, and 4 Taro Yamauchi 1 National Agricultural
More informationSnowcover interaction with climate, topography & vegetation in mountain catchments
Snowcover interaction with climate, topography & vegetation in mountain catchments DANNY MARKS Northwest Watershed Research Center USDA-Agricultural Agricultural Research Service Boise, Idaho USA RCEW
More informationHua Zhang, Steven M. Gorelick, Nicolas Avisse, Amaury Tilmant, Deepthi Rajsekhar, and Jim Yoon
Remote Sens. 216,, 735; doi:1.339/rs9735 S1 of S5 Supplementary Materials: A New Temperature- Vegetation Triangle Algorithm with Variable Edges (VE) for Satellite-Based Actual Evapotranspiration Estimation
More informationThe National Operational Hydrologic Remote Sensing Center Operational Snow Analysis
The National Operational Hydrologic Remote Sensing Center Operational Snow Analysis World Meteorological Organization Global Cryosphere Watch Snow-Watch Workshop Session 3: Snow Analysis Products Andrew
More informationDaily, Monthly and Yearly Norm Temperature ( ): TnormD8110, TnormM8110 and TnormY8110
Federal Department of Home Affairs FDHA Federal Office of Meteorology and Climatology MeteoSwiss Documentation of MeteoSwiss Grid-Data Products Daily, Monthly and Yearly Norm Temperature (1981-2010): TnormD8110,
More informationEffects of forest cover and environmental variables on snow accumulation and melt
Effects of forest cover and environmental variables on snow accumulation and melt Mariana Dobre, William J. Elliot, Joan Q. Wu, Timothy E. Link, Ina S. Miller Abstract The goal of this study was to assess
More informationComparing different MODIS snow products with distributed simulation of the snowpack in the French Alps
Comparing different MODIS snow products with distributed simulation of the snowpack in the French Alps Luc Charrois 1, Marie Dumont 1,*, Pascal Sirguey 2, Samuel Morin 1, Matthieu Lafaysse 1 and Fatima
More informationClimate Impacts to Southwest Water Sector. Dr. Dave DuBois New Mexico State Climatologist
Climate Impacts to Southwest Water Sector Dr. Dave DuBois New Mexico State Climatologist Presented at EPA Extreme Events and Climate Adaptation Planning Workshop June 3, 2015 Water Year Precip. % of Average
More informationSubgrid variability of snow water equivalent at operational snow stations in the western USA
HYDROLOGICAL PROCESSES Hydrol. Process. (212) Published online in Wiley Online Library (wileyonlinelibrary.com) DOI: 1.12/hyp.9355 Subgrid variability of snow water equivalent at operational snow stations
More informationGreat Lakes Update. Volume 193: 2015 January through June Summary. Vol. 193 Great Lakes Update August 2015
Great Lakes Update Volume 193: 2015 January through June Summary The U.S. Army Corps of Engineers (USACE) monitors the water levels of each of the Great Lakes. This report provides a summary of the Great
More informationInvestigating Mechanisms of Cool Season Upper Colorado River Basin Precipitation
Investigating Mechanisms of Cool Season Upper Colorado River Basin Precipitation Benjamin Hatchett 1,2, Nina Oakley 1,2, John Abatzoglou 3, and Jon Rutz 4 1 Division of Atmospheric Sciences, Desert Research
More informationSupplementary Materials for
advances.sciencemag.org/cgi/content/full/3/12/e1701169/dc1 Supplementary Materials for Abrupt shift in the observed runoff from the southwestern Greenland ice sheet Andreas P. Ahlstrøm, Dorthe Petersen,
More informationStudying snow cover in European Russia with the use of remote sensing methods
40 Remote Sensing and GIS for Hydrology and Water Resources (IAHS Publ. 368, 2015) (Proceedings RSHS14 and ICGRHWE14, Guangzhou, China, August 2014). Studying snow cover in European Russia with the use
More informationSensor placement strategies for snow water equivalent (SWE) estimation in the American River basin
WATER RESOURCES RESEARCH, VOL. 49, 891 903, doi:10.1002/wrcr.20100, 2013 Sensor placement strategies for snow water equivalent (SWE) estimation in the American River basin Stephen C. Welch, 1 Branko Kerkez,
More informationEvaluation of a MODIS Triangle-based Algorithm for Improving ET Estimates in the Northern Sierra Nevada Mountain Range
Evaluation of a MODIS Triangle-based Algorithm for Improving ET Estimates in the Northern Sierra Nevada Mountain Range Kyle R. Knipper 1, Alicia M. Kinoshita 2, and Terri S. Hogue 1 January 5 th, 2015
More informationModelling runoff from large glacierized basins in the Karakoram Himalaya using remote sensing of the transient snowline
Remote Sensing and Hydrology 2000 (Proceedings of a symposium held at Santa Fe, New Mexico, USA, April 2000). IAHS Publ. no. 267, 2001. 99 Modelling runoff from large glacierized basins in the Karakoram
More informationPrecipitation. Standardized Precipitation Index. NIDIS Intermountain West Regional Drought Early Warning System January 3, 2017
1/3/2017 NIDIS Drought and Water Assessment NIDIS Intermountain West Regional Drought Early Warning System January 3, 2017 Precipitation The images above use daily precipitation statistics from NWS COOP,
More informationSnow Melt with the Land Climate Boundary Condition
Snow Melt with the Land Climate Boundary Condition GEO-SLOPE International Ltd. www.geo-slope.com 1200, 700-6th Ave SW, Calgary, AB, Canada T2P 0T8 Main: +1 403 269 2002 Fax: +1 888 463 2239 Introduction
More informationA BRIEF HISTORY OF EVALUATIONS PERFORMED ON THE OPERATIONAL KINGS RIVER WINTER OROGRAPHIC CLOUD SEEDING PROGRAM. Don A. Griffith and David P.
AUGUST 2014 GRIFFITH ET AL. 29 A BRIEF HISTORY OF EVALUATIONS PERFORMED ON THE OPERATIONAL KINGS RIVER WINTER OROGRAPHIC CLOUD SEEDING PROGRAM Don A. Griffith and David P. Yorty North American Weather
More informationSouris River Basin Spring Runoff Outlook As of March 1, 2019
Souris River Basin Spring Runoff Outlook As of March 1, 2019 Prepared by: Flow Forecasting & Operations Planning Water Security Agency Basin Conditions Summer rainfall in 2018 over the Saskatchewan portion
More informationRegional offline land surface simulations over eastern Canada using CLASS. Diana Verseghy Climate Research Division Environment Canada
Regional offline land surface simulations over eastern Canada using CLASS Diana Verseghy Climate Research Division Environment Canada The Canadian Land Surface Scheme (CLASS) Originally developed for the
More informationDrought in Southeast Colorado
Drought in Southeast Colorado Nolan Doesken and Roger Pielke, Sr. Colorado Climate Center Prepared by Tara Green and Odie Bliss http://climate.atmos.colostate.edu 1 Historical Perspective on Drought Tourism
More information