The temporal and spatial continuity of spatially distributed estimates of snow-covered area (SCA) are limited by the availability of cloud-free satellite imagery; this also affects spatial estimates of snow water equivalent (SWE), as SCA can be used to define the extent of snow telemetry (SNOTEL) point SWE interpolation. In order to extend the continuity of these estimates in time and space to areas beneath the cloud cover, gridded temperature data were used to define the spatial domain of SWE interpolation in the Salt-Verde watershed of Arizona. Gridded positive accumulated degree-days (ADD) and binary SCA (derived from the Advanced Very High Resolution Radiometer (AVHRR)) were used to define a threshold ADD to define the area of snow cover. The optimized threshold ADD increased during snow accumulation periods, reaching a peak at maximum snow extent. The threshold then decreased dramatically during the first time period after peak snow extent owing to the low amount of energy required to melt the thin snow cover at lower elevations. The area having snow cover at this later time was then used to define the area for which SWE interpolation was done. The area simulated to have snow was compared with observed SCA from AVHRR to assess the simulated snow map accuracy. During periods without precipitation, the average commission and omission errors of the optimal technique were 7% and 11% respectively, with a map accuracy of 82%. Average map accuracy decreased to 75% during storm periods, with commission and omission errors equal to 11% and 12% respectively. The analysis shows that temperature data can be used to help estimate the snow extent beneath clouds and therefore improve the spatial and temporal continuity of SCA and SWE products. © 2004 John Wiley and Sons, Ltd.

Regression tree models have been shown to provide the most accurate estimates of distributed snow water equivalent (SWE) when intensive field observations are available. This work presents a comparison of regression tree models using different source digital elevation models (DEMs) and different combinations of independent variables. Different residual interpolation techniques are also compared. The analysis was performed in the 19.1 km2 Tokopah Basin, located in the southern Sierra Nevada of California. Snow depth, the dependent variable of the statistical models, was derived from three snow surveys (April, May and June 1997), with an average of 328 depth measurements per survey. Estimates of distributed SWE were derived from the product of the snow depth surfaces, the average snow density (54 measurements on average) and the fractional snow covered area (obtained from the Landsat Thematic Mapper and the Airborne Visible/Infrared Imaging Spectrometer). Independent variables derived from the standard US Geological Survey DEM yielded the lowest overall model deviance and lowest error in snow depth prediction. Simulations using the Shuttle Radar Topography Mission DEM and the National Elevation Dataset DEM were improved when northness was substituted for solar radiation in five of six cases. Co-kriging with maximum upwind slope and elevation proved to be the best method for distributing residuals for April and June, respectively. Inverse distance weighting was the best residual distribution method for May. Copyright © 2004 John Wiley & Sons, Ltd.

Abstract. We present results from snow-on and snow-off airborne-scanning LiDAR measurements over a 53 km2 area in the southern Sierra Nevada. We found that snow depth as a function of elevation increased approximately 15 cm per 100 m, until reaching an elevation of 3300 m, where depth sharply decreased at a rate of 48 cm per 100 m. Departures from the 15 cm per 100 m trend, based on 1 m elevation-band means of regression residuals, showed slightly less steep increases below 2050 m; steeper increases between 2050 and 3300 m; and less steep increases above 3300 m. Although the study area is partly forested, only measurements in open areas were used. Below approximately 2050 m elevation, ablation and rainfall are the primary causes of departure from the orographic trend. From 2050 to 3300 m, greater snow depths than predicted were found on the steeper terrain of the northwest and the less steep northeast-facing slopes, suggesting that ablation, aspect, slope and wind redistribution all play a role in local snow-depth variability. At elevations above 3300 m, orographic processes mask the effect of wind deposition when averaging over large areas. Also, terrain in this basin becomes less steep above 3300 m. This suggests a reduction in precipitation from upslope lifting and/or the exhaustion of precipitable water from ascending air masses. Our results suggest a cumulative precipitation lapse rate for the 2100–3300 m range of about 6 cm per 100 m elevation for the accumulation period of 3 December 2009 to 23 March 2010. This is a higher gradient than the widely used PRISM (Parameter-elevation Relationships on Independent Slopes Model) precipitation products, but similar to that from reconstruction of snowmelt amounts from satellite snow-cover data. Our findings provide a unique characterization of the consistent, steep average increase in precipitation with elevation in snow-dominated terrain, using high-resolution, highly accurate data and highlighs the importance of solar radiation, wind redistribution and mid-winter melt with regard to snow distribution.