Deriving snow‐cover depletion curves for different spatial scales from remote sensing and snow telemetry data

During the melting of a snowpack, snow water equivalent (SWE) can be correlated to snow‐covered area (SCA) once snow‐free areas appear, which is when SCA begins to decrease below 100%. This amount of SWE is called the threshold SWE. Daily SWE data from snow telemetry stations were related to SCA derived from moderate‐resolution imaging spectroradiometer images to produce snow‐cover depletion curves. The snow depletion curves were created for an 80 000 km² domain across southern Wyoming and northern Colorado encompassing 54 snow telemetry stations. Eight yearly snow depletion curves were compared, and it is shown that the slope of each is a function of the amount of snow received. Snow‐cover depletion curves were also derived for all the individual stations, for which the threshold SWE could be estimated from peak SWE and the topography around each station. A station's peak SWE was much more important than the main topographic variables that included location, elevation, slope, and modelled clear sky solar radiation. The threshold SWE mostly illustrated inter‐annual consistency. Copyright © 2015 John Wiley & Sons, Ltd.     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

Spatio‐temporal variability of Great Lakes basin snow cover ablation events

In the Great Lakes basin of North America, annual run‐off is dominated by snowmelt. This snowmelt‐induced run‐off plays an important role within the hydrologic cycle of the basin, influencing soil moisture availability and driving the seasonal cycle of spring and summer lake levels. Despite this, relatively little is understood about the patterns and trends of snow ablation event frequency and magnitude within the Great Lakes basin. This study uses a gridded dataset of Canadian and United States surface snow depth observations to develop a regional climatology of snow ablation events from 1960 to 2009. An ablation event is defined as an interdiurnal snow depth decrease within an individual grid cell. A clear seasonal cycle in ablation event frequency exists within the basin and peak ablation event probability is latitudinally dependent. Most of the basin experiences peak ablation frequency in March, while the northern and southern regions of the basin experience respective peaks in April and February. An investigation into the interannual frequency of ablation events reveals ablation events significantly decrease within the northeastern and northwestern Lake Superior drainage basins and significantly increase within the eastern Lake Huron and Georgian Bay drainage basins. In the eastern Lake Huron and Georgian Bay drainage basins, larger ablation events are occurring more frequently, and a larger impact to the hydrology can be expected. Trends in ablation events are attributed primarily to changes in snowfall and snow depth across the region.     

Factors controlling seasonal groundwater and solute flux from snow‐dominated basins

Critical zone influences on hydrologic partitioning, subsurface flow paths and reactions along these flow paths dictate the timing and magnitude of groundwater and solute flux to streams. To isolate first‐order controls on seasonal streamflow generation within highly heterogeneous, snow‐dominated basins of the Colorado River, we employ a multivariate statistical approach of end‐member mixing analysis using a suite of daily chemical and isotopic observations. Mixing models are developed across 11 nested basins (0.4 to 85 km²) spanning a gradient of climatological, physical, and geological characteristics. Hydrograph separation using rain, snow, and groundwater as end‐members indicates that seasonal contributions of groundwater to streams is significant. Mean annual groundwater flux ranges from 12% to 33% whereas maximum groundwater contributions of 17% to 50% occur during baseflow. The direct relationship between snow water equivalent and groundwater flux to streams is scale dependent with a trend toward self‐similarity when basins exceed 5.5 km². We find groundwater recharge increases in basins of high relief and within the upper subalpine where maximum snow accumulation is coincident with reduced conifer cover and lower canopy densities. The mixing model developed for the furthest downstream site did not transfer to upstream basins. The resulting error in predicted stream concentrations points toward weathering reactions as a function of source rock and seasonal shifts in flow path. Additionally, the potential for microbial sulfate reduction in floodplain sediments along a low‐gradient, meandering portion of the river is sufficient to modify hillslope contributions and alter mixing ratios in the analysis. Soil flushing in response to snowmelt is not included as an end‐member but is identified as an important mechanism for release of solutes from these mountainous watersheds. End‐member mixing analysis used in combination with high‐frequency observations reveals important aspects of catchment hydrodynamics across scale.     

Evaluating the functionality and streamflow impacts of explicitly modelling forest–snow interactions and canopy gaps in a distributed hydrologic model

Many plot‐scale studies have shown that snow‐cover dynamics in forest gaps are distinctly different from those in open and continuously forested areas, and forest gaps have the potential to alter the magnitude and timing of snowmelt. However, the watershed‐level impacts of canopy gap treatment on streamflows are largely unknown. Here, we present the first research that explicitly assesses the impact of canopy gaps on seasonal streamflows and particularly late‐season low flows at the watershed scale. To explicitly model forest–snow interactions in canopy gaps, we made major enhancements to a widely used distributed hydrologic model, distributed hydrology soil vegetation model, with a canopy gap component that represents physical processes of snowpack evolution in the forest gap separately from the surrounding forest on the subgrid scale (within a grid typically 10–150 m). The model predicted snow water equivalent using the enhanced distributed hydrology soil vegetation model showed good agreement (R² > 0.9) with subhourly snow water equivalent measurements collected from open, forested, and canopy gap sites in Idaho, USA. Compared with the original model that does not account for interactions between gaps and surrounding forest, the enhanced model predicted notably later melt in small‐ to medium‐size canopy gaps (the ratio of gap radius (r) to canopy height (h) ≤ 1.2), and snow melt rates exhibited great sensitivity to changing gap size in medium‐size gaps (0.5 ≤ r/h ≤ 1.2). We demonstrated the watershed‐scale implications of canopy gaps on streamflow in the snow‐dominated Chiwawa watershed, WA, USA. With 24% of the watershed drainage area (about 446 km²) converted to gaps of 60 m diameter, the mean annual 7‐day low flow was increased by 19.4% (i.e., 0.37 m³/s), and the mean monthly 7‐day low flows were increased by 13.5% (i.e., 0.26 m³/s) to 40% (i.e., 1.76 m³/s) from late summer through fall. Lastly, in practical implementation of canopy gaps with the same total gap areas, a greater number of distributed small gaps can have greater potential for longer snow retention than a smaller number of large gaps.     

Understanding hydrological processes in glacierized catchments: Evidence and implications of highly variable isotopic and electrical conductivity data

The spatial and temporal characterization of geochemical tracers over Alpine glacierized catchments is particularly difficult, but fundamental to quantify groundwater, glacier melt, and rain water contribution to stream runoff. In this study, we analysed the spatial and temporal variability of δ²H and electrical conductivity (EC) in various water sources during three ablation seasons in an 8.4‐km² glacierized catchment in the Italian Alps, in relation to snow cover and hydro‐meteorological conditions. Variations in the daily streamflow range due to melt‐induced runoff events were controlled by maximum daily air temperature and snow covered area in the catchment. Maximum daily streamflow decreased with increasing snow cover, and a threshold relation was found between maximum daily temperature and daily streamflow range. During melt‐induced runoff events, stream water EC decreased due to the contribution of glacier melt water to stream runoff. In this catchment, EC could be used to distinguish the contribution of subglacial flow (identified as an end member, enriched in EC) from glacier melt water to stream runoff, whereas spring water in the study area could not be considered as an end member. The isotopic composition of snow, glacier ice, and melt water was not significantly correlated with the sampling point elevation, and the spatial variability was more likely affected by postdepositional processes. The high spatial and temporal variability in the tracer signature of the end members (subglacial flow, rain water, glacier melt water, and residual winter snow), together with small daily variability in stream water δ²H dynamics, are problematic for the quantification of the contribution of the identified end members to stream runoff, and call for further research, possibly integrated with other natural or artificial tracers.     

Exploration of remotely sensed forest structure and ultrasonic range sensor metrics to improve empirical snow models

Current methods to estimate snow accumulation and ablation at the plot and watershed levels can be improved as new technologies offer alternative approaches to more accurately monitor snow dynamics and their drivers. Here we conduct a meta‐analysis of snow and vegetation data collected in British Columbia to explore the relationships between a wide range of forest structure variables – obtained from Light Detection and Ranging (LiDAR), hemispherical photography (HP) and Landsat Thematic Mapper – and several indicators of snow accumulation and ablation estimated from manual snow surveys and ultrasonic range sensors. By merging and standardizing all the ground plot information available in the study area, we demonstrate how LiDAR‐derived forest cover above 0.5 m was the variable explaining the highest percentage of absolute peak snow water equivalent (SWE) (33%), while HP‐derived leaf area index and gap fraction (45° angle of view) were the best potential predictors of snow ablation rate (explaining 57% of variance). This study reveals how continuous SWE data from ultrasonic sensors are fundamental to obtain statistically significant relationships between snow indicators and structural metrics by increasing mean r² by 20% when compared to manual surveys. The relationships between vegetation and spectral indices from Landsat and snow indicators, not explored before, were almost as high as those shown by LiDAR or HP and thus point towards a new line of research with important practical implications. While the use of different data sources from two snow seasons prevented us from developing models with predictive capacity, a large sample size helped to identify outliers that weakened the relationships and suggest improvements for future research. A concise overview of the limitations of this and previous studies is provided along with propositions to consistently improve experimental designs to take advantage of remote sensing technologies, and better represent spatial and temporal variations of snow. Copyright © 2013 John Wiley & Sons, Ltd.     

Development and application of a spatially‐distributed Arctic hydrological and thermal process model (ARHYTHM)

A process‐based, spatially distributed hydrological model was developed to quantitatively simulate the energy and mass transfer processes and their interactions within arctic regions (arctic hydrological and thermal model, ARHYTHM). The model first determines the flow direction in each element, the channel drainage network and the drainage area based upon the digital elevation data. Then it simulates various physical processes: including snow ablation, subsurface flow, overland flow and channel flow routing, soil thawing and evapotranspiration. The kinematic wave method is used for conducting overland flow and channel flow routing. The subsurface flow is simulated using the Darcian approach. The energy balance scheme was the primary approach used in energy‐related process simulations (snowmelt and evapotranspiration), although there are options to model snowmelt by the degree‐day method and evapotranspiration by the Priestley–Taylor equation. This hydrological model simulates the dynamic interactions of each of these processes and can predict spatially distributed snowmelt, soil moisture and evapotranspiration over a watershed at each time step as well as discharge in any specified channel(s). The model was applied to Imnavait watershed (about 2·2 km²) and the Upper Kuparuk River basin (about 146 km²) in northern Alaska. Simulated results of spatially distributed soil moisture content, discharge at gauging stations, snowpack ablations curves and other results yield reasonable agreement, both spatially and temporally, with available data sets such as SAR imagery‐generated soil moisture data and field measurements of snowpack ablation, and discharge data at selected points. The initial timing of simulated discharge does not compare well with the measured data during snowmelt periods mainly because the effect of snow damming on runoff was not considered in the model. Results from the application of this model demonstrate that spatially distributed models have the potential for improving our understanding of hydrology for certain settings. Finally, a critical component that led to the performance of this modelling is the coupling of the mass and energy processes. Copyright © 2000 John Wiley & Sons, Ltd.     

Automatic mapping of snow cover depletion curves using optical remote sensing data under conditions of frequent cloud cover and temporary snow

Snow cover depletion curves are required for several water management applications of snow hydrology and are often difficult to obtain automatically using optical remote sensing data owing to both frequent cloud cover and temporary snow cover. This study develops a methodology to produce accurate snow cover depletion curves automatically using high temporal resolution optical remote sensing data (e.g. Terra Moderate Resolution Imaging Spectroradiometer (MODIS), Aqua MODIS or National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR)) by snow cover change trajectory analysis. The method consists of four major steps. The first is to reclassify both cloud‐obscured land and snow into more distinct subclasses and to determine their snow cover status (seasonal snow cover or not) based on the snow cover change trajectories over the whole snowmelt season. The second step is to derive rules based on the analysis of snow cover change trajectories. These rules are subsequently used to determine for a given date, the snow cover status of a pixel based on snow cover maps from the beginning of the snowmelt season to that given date. The third step is to apply a decision‐tree‐like processing flow based on these rules to determine the snow cover status of a pixel for a given date and to create daily seasonal snow cover maps. The final step is to produce snow cover depletion curves using these maps. A case study using this method based on Terra MODIS snow cover map products (MOD10A1) was conducted in the lower and middle reaches of the Kaidu River Watershed (19 000 km²) in the Chinese Tien Shan, Xinjiang Uygur Autonomous Region, China. High resolution remote sensing data (charge coupled device (CCD) camera data with 19·5 m resolution of the China and Brazil Environmental and Resources Satellite (CBERS) data (19·5 m resolution), and the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data with 15 m resolution of the Terra) were used to validate the results. The study shows that the seasonal snow cover classification was consistent with that determined using a high spatial resolution dataset, with an accuracy of 87–91%. The snow cover depletion curves clearly reflected the impact of the variation of temperature and the appearance of temporary snow cover on seasonal snow cover. The findings from this case study suggest that the approach is successful in generating accurate snow cover depletion curves automatically under conditions of frequent cloud cover and temporary snow cover using high temporal resolution optical remote sensing data. Copyright © 2007 John Wiley & Sons, Ltd.     

Snowpack variability across various spatio‐temporal resolutions

High‐resolution snow depth (SD) maps (1 × 1 m) obtained from terrestrial laser scanner measurements in a small catchment (0.55 km²) in the Pyrenees were used to assess small‐scale variability of the snowpack at the catchment and sub‐grid scales. The coefficients of variation are compared for various plot resolutions (5 × 5, 25 × 25, 49 × 49, and 99 × 99 m) and eight different days in two snow seasons (2011–2012 and 2012–2013). We also studied the relation between snow variability at the small scale and SD, topographic variables, small‐scale variability in topographic variables. The results showed that there was marked variability in SD, and it increased with increasing scales. Days of seasonal maximum snow accumulation showed the least small‐scale variability, but this increased sharply with the onset of melting. The coefficient of variation (CV) in snowpack depth showed statistically significant consistency amongst the various spatial resolutions studied, although it declined progressively with increasing difference between the grid sizes being compared. SD best explained the spatial distribution of sub‐grid variability. Topographic variables including slope, wind sheltering, sub‐grid variability in elevation, and potential incoming solar radiation were also significantly correlated with the CV of the snowpack, with the greatest correlation occurring at the 99 × 99 m resolution. At this resolution, stepwise multiple regression models explained more than 70% of the variance, whereas at the 25 × 25 m resolution they explained slightly more than 50%. The results highlight the importance of considering small‐scale variability of the SD for comprehensively representing the distribution of snowpack from available punctual information, and the potential for using SD and other predictors to design optimized surveys for acquiring distributed SD data. Copyright © 2014 John Wiley & Sons, Ltd.     

Exploitation of Landsat data for snow zonation mapping in the Hindukush,Karakoram and Himalaya (HKH) region of Pakistan

Abstract

In the Hindukush, Karakoram and Himalaya (HKH) region of Pakistan, many glaciological variables are still not known due to the remoteness and harsh weather conditions of the area. A remote sensing technique is therefore applied to map the snow zonation in the HKH region. Landsat 7 ETM+ data for the year 2003 are used in this study. Image classification and image processing techniques are applied to map, for the first time, the major snow zones in the HKH region. Six classes are identified: the results show that the area covered by the highest-altitude snow (Snow I), lower-altitude snow (Snow II), bare ice, debris-covered ice, wet snow and shadow is 21 529.42, 22 472.58, 8696.41, 8038.75, 12 159.37 and 7322.30 km², respectively. The study also indicates that the equilibrium line altitude (ELA) lies between 5000 and 5500 m above sea level, with an accumulation area ratio (AAR) of 0.60.

Citation Butt, M.J., 2013. Exploitation of Landsat data for snow zonation mapping in the Hindukush, Karakoram and Himalaya (HKH) region of Pakistan. Hydrological Sciences Journal, 58 (5), 1088–1096.     

Snow in Lebanon: a preliminary study of snow cover over Mount Lebanon and a simple snowmelt model / Etude préliminaire du couvert neigeux et modèle de fonte des neige pour le Mont Liban

Abstract

The physical properties of snow, including apparent density, snow cover distribution and snowmelt in the Nahr El Kelb basin (Mount Lebanon), were studied in order to design a simple empirical snowmelt model. In February 2001, snow covered an area of 1600 km² on Mount Lebanon, representing a water equivalent of 1.1 x 10⁹ m³. The snow surface area was calculated by combining TM5 images with a digital elevation model, and field observations made every three days, from 1400 to 2300 m altitude. The depletion of snow cover was measured from the end of December 2000 to the end of June 2001. The snowmelt was measured from surface depletion on a degree-day basis. A simple model relating the daily snowmelt to the product of wind speed and average positive daily air temperature, is presented and discussed. For Mount Lebanon, this model gave a better approximation of snowmelt than a simple degree-day model.     

Snow accumulation distribution inferred from time‐lapse photography and simple modelling

The spatial and temporal distribution of snow accumulation is complex and significantly influences the hydrological characteristics of mountain catchments. Many snow redistribution processes, such as avalanching, slushflow or wind drift, are controlled by topography, but their modelling remains challenging. In situ measurements of snow accumulation are laborious and generally have a coarse spatial or temporal resolution. In this respect, time‐lapse photography shows itself as a powerful tool for collecting information at relatively low cost and without the need for direct field access. In this paper, the snow accumulation distribution of an Alpine catchment is inferred by adjusting a simple snow accumulation model combined with a temperature index melt model to match the modelled melt‐out pattern evolution to the pattern monitored during an ablation season through terrestrial oblique photography. The comparison of the resulting end‐of‐winter snow water equivalent distribution with direct measurements shows that the achieved accuracy is comparable with that obtained with an inverse distance interpolation of the point measurements. On average over the ablation season, the observed melt‐out pattern can be reproduced correctly in 93% of the area visible from the fixed camera. The relations between inferred snow accumulation distribution and topographic variables indicate large scatter. However, a significant correlation with local slope is found and terrain curvature is detected as a factor limiting the maximal snow accumulation. Copyright © 2010 John Wiley & Sons, Ltd.     

Measuring snow ablation rates in alpine terrain with a mobile multioffset ground‐penetrating radar system

Ground‐penetrating radar (GPR) has become a promising technique in the field of snow hydrological research. It is commonly used to measure snow depth, density, and water equivalent over large distances or along gridded snow courses. Having built and tested a mobile lightweight set‐up, we demonstrate that GPR is capable of accurately measuring snow ablation rates in complex alpine terrain. Our set‐up was optimized for efficient measurements and consisted of a multioffset radar with four pairs of antennas mounted to a plastic sled, which was small enough to permit safe and convenient operations. Repeated measurements at intervals of 2 to 7 days were taken during the 2014/2015 winter season along 10 profiles of 50 to 200 m length within two valleys located in the eastern Swiss Alps. Resulting GPR‐based data of snow depth, density, and water equivalent, as well as their respective change over time, were in good agreement with concurrent manual measurements, in particular if accurate alignment between repeated overpasses could be achieved. Corresponding root‐mean‐square error (RMSE) values amounted to 4.2 cm for snow depth, 17 mm for snow water equivalent, and 22 kg/m³ for snow density, with similar RMSE values for corresponding differential data. With this performance, the presented radar set‐up has the potential to provide exciting new and extensive datasets to validate snowmelt models or to complement lidar‐based snow surveys.     

Role of temporal resolution of meteorological inputs for process‐based snow modelling

Accurate snow accumulation and melt simulations are crucial for understanding and predicting hydrological dynamics in mountainous settings. As snow models require temporally varying meteorological inputs, time resolution of these inputs is likely to play an important role on the model accuracy. Because meteorological data at a fine temporal resolution (~1 hr) are generally not available in many snow‐dominated settings, it is important to evaluate the role of meteorological inputs temporal resolution on the performance of process‐based snow models. The objective of this work is to assess the loss in model accuracy with temporal resolution of meteorological inputs, for a range of climatic conditions and topographic elevations. To this end, a process‐based snow model was run using 1‐, 3‐, and 6‐hourly inputs for wet, average, and dry years over Boise River Basin (6,963 km²), which spans rain dominated (≤1,400 m), rain–snow transition (>1,400 and ≤1,900 m), snow dominated below tree line (>1,900 and ≤2,400 m), and above tree line (>2,400 m) elevations. The results show that sensitivity of the model accuracy to the inputs time step generally decreases with increasing elevation from rain dominated to snow dominated above tree line. Using longer than hourly inputs causes substantial underestimation of snow cover area (SCA) and snow water equivalent (SWE) in rain‐dominated and rain–snow transition elevations, due to the precipitation phase mischaracterization. In snow‐dominated elevations, the melt rate is underestimated due to errors in estimation of net snow cover energy input. In addition, the errors in SCA and SWE estimates generally decrease toward years with low snow mass, that is, dry years. The results indicate significant increases in errors in estimates of SCA and SWE as the temporal resolution of meteorological inputs becomes coarser than an hour. However, use of 3‐hourly inputs can provide accurate estimates at snow‐dominated elevations. The study underscores the need to record meteorological variables at an hourly time step for accurate process‐based snow modelling.     

Effects of land‐cover changes on the hydrological response of interior Columbia River basin forested catchments

The topographically explicit distributed hydrology–soil–vegetation model (DHSVM) is used to simulate hydrological effects of changes in land cover for four catchments, ranging from 27 to 1033 km², within the Columbia River basin. Surface fluxes (stream flow and evapotranspiration) and state variables (soil moisture and snow water equivalent) corresponding to historical (1900) and current (1990) vegetation are compared. In addition a sensitivity analysis, where the catchments are covered entirely by conifers at different maturity stages, was conducted. In general, lower leaf‐area index (LAI) resulted in higher snow water equivalent, more stream flow and less evapotranspiration. Comparisons with the macroscale variable infiltration capacity (VIC) model, which parameterizes, rather than explicitly represents, topographic effects, show that runoff predicted by DHSVM is more sensitive to land‐cover changes than is runoff predicted by VIC. This is explained by model differences in soil parameters and evapotranspiration calculations, and by the more explicit representation of saturation excess in DHSVM and its higher sensitivity to LAI changes in the calculation of evapotranspiration. Copyright © 2002 John Wiley & Sons, Ltd.     

Dynamics of snow ablation in a small Alpine catchment observed by repeated terrestrial laser scans

The spatio‐temporal distribution of snow in a catchment during ablation reflects changes in the total amount of snow water equivalent and is thus a key parameter for the estimation of melt water run‐off. This study explores possible rules behind the spatial variability of snow depth during the ablation season in a small Alpine catchment with complex topography. The snow depth observations are based on more than 160 000 terrestrial laser scanner data points with a spatial resolution of 1 m, which were obtained from 11 scanning campaigns of two consecutive ablation seasons. The analysis suggests that for estimating cumulative snow melt dynamics from the catchment investigated, assessing the initial snow distribution prior to the melt season is more important than addressing spatial differences in the melt behaviour. Snow volume and snow‐covered area could be predicted well using a conceptual melt model assuming spatially uniform melt rates. However, accurate results were only obtained if the model was initialized with a pre‐melt snow distribution that reflected measured mean and standard deviation. Using stratified melt rates on the other hand did not improve the model results. At least for sites with similar meteorological and topographical conditions, the model approach presented here comprises an efficient way to estimate snow depletion dynamics, especially if persistent snow accumulation pattern between years facilitate the characterization of the initial snow distribution prior to the melt. Copyright © 2011 John Wiley & Sons, Ltd.     

Estimating snow distribution over a large area and its application for water resources

The spatial and temporal distribution of the snow water equivalent (SWE), snow density and snow depth were estimated by a method combining remote sensing technology and degree‐day techniques over a study area of 370 000 km². The advantages of this simulation model are its simplicity and the availability of degree‐day parameters, which can be successively evaluated by referring to snow area maps created from satellite images. This simulation worked very well for estimating SWE and helped to separate the areas of thin snow cover from heavier snowfall. However, shallow snow in warm regions led to some misjudgments in the snow area maps because of the time lag between when the satellite image was acquired and the simulation itself. Vulnerable areas, where a large variation in the amount of snow affects people's life, could be identified from the differences between heavy and light snow years. This vulnerability stems from a predicted lack of irrigation water for rice production caused by future climate change. The model developed in this study has the potential to contribute to water management activities and decision‐making processes when considering necessary adaptations to future climate change. Copyright © 2007 John Wiley & Sons, Ltd.     

Improvement of distributed snowmelt energy balance modeling with MODIS‐based NDSI‐derived fractional snow‐covered area data

Describing the spatial variability of heterogeneous snowpacks at a watershed or mountain‐front scale is important for improvements in large‐scale snowmelt modelling. Snowmelt depletion curves, which relate fractional decreases in snow‐covered area (SCA) against normalized decreases in snow water equivalent (SWE), are a common approach to scale‐up snowmelt models. Unfortunately, the kinds of ground‐based observations that are used to develop depletion curves are expensive to gather and impractical for large areas. We describe an approach incorporating remotely sensed fractional SCA (FSCA) data with coinciding daily snowmelt SWE outputs during ablation to quantify the shape of a depletion curve. We joined melt estimates from the Utah Energy Balance Snow Accumulation and Melt Model (UEB) with FSCA data calculated from a normalized difference snow index snow algorithm using NASA's moderate resolution imaging spectroradiometer (MODIS) visible (0·545–0·565 µm) and shortwave infrared (1·628–1·652 µm) reflectance data. We tested the approach at three 500 m² study sites, one in central Idaho and the other two on the North Slope in the Alaskan arctic. The UEB‐MODIS‐derived depletion curves were evaluated against depletion curves derived from ground‐based snow surveys. Comparisons showed strong agreement between the independent estimates. Copyright © 2010 John Wiley & Sons, Ltd.     

A comparison of wetness indices for the prediction of observed connected saturated areas under contrasting conditions

For lack of other widely available spatial information, topography is often used to predict water fluxes and water quality in mesoscale watersheds. Such data have however proven to be misleading in many environments where large and flat valley bottoms and/or highly conducive soil covers determine water storage and water transport mechanisms. Also, the focus is generally on the prediction of saturation areas regardless of whether they are connected to the catchment hydrographic network or rather present in isolated topographic depressions. Here soil information was coupled with terrain data towards the targeted prediction of connected saturated areas. The focus was on the 30 km² Girnock catchment (Cairngorm Mountains, northeast Scotland) and its 3 km² sub‐catchment, Bruntland Burn in which seven field surveys were done to capture actual maps of connected saturated areas in both dry and humid conditions. The 1 km² resolution UK Hydrology of Soil Types (HOST) classification was used to extract relevant, spatially variable, soil parameters. Results show that connected saturated areas were fairly well predicted by wetness indices but only in wet conditions when they covered more than 30% of the whole catchment area. Geomorphic indices including information on terrain shape, steepness, aspect, soil texture and soil depth showed potential but generally performed poorly. Indices based on soil and topographic data did not have more predictive power than those based on topographic information only: this was attributed to the coarse resolution of the HOST classification. Nevertheless, analyses provided interesting insights into the scale‐dependent water storage and transport mechanisms in both study catchments. Copyright © 2013 John Wiley & Sons, Ltd.