VARIABILITY OF MELTWATER AND SOLUTE FLUXES FROM HOMOGENEOUS MELTING SNOW AT THE LABORATORY SCALE

Four experiments were performed to examine the relationship between the meltwater flow field and ion release from melting snow. A 0.4 m³ volume of snow was placed in a Plexiglass box and melted from above using a heating plate. The meltwater and solute fluxes issuing from the bottom of the snow were monitored. In experiments with NaCl tracer added to the snow, the solute concentrations were generally lower in the flow fingers than in the background wetting front. Dye tracer experiments revealed contemporaneous areas of concentrated dye and dilute meltwater in flow fingers. This suggests that the meltwater in flow fingers is diluted by low concentration water from the top of the snowpack. Flow fingers contribute more meltwater flux primarily because the flow is maintained for a longer period of time than in the non-finger areas; however, the relative contribution of flow fingers to solute flux was apparently not as great as that of the background wetting front because of dilution of solute in the flow finger areas.     

Charred forests accelerate snow albedo decay: parameterizing the post‐fire radiative forcing on snow for three years following fire

As large, high‐severity forest fires increase and snowpacks become more vulnerable to climate change across the western USA, it is important to understand post‐fire disturbance impacts on snow hydrology. Here, we examine, quantify, parameterize, model, and assess the post‐fire radiative forcing effects on snow to improve hydrologic modelling of snow‐dominated watersheds having experienced severe forest fires. Following a 2011 high‐severity forest fire in the Oregon Cascades, we measured snow albedo, monitored snow, and micrometeorological conditions, sampled snow surface debris, and modelled snowpack energy and mass balance in adjacent burned forest (BF) and unburned forest sites. For three winters following the fire, charred debris in the BF reduced snow albedo, accelerated snow albedo decay, and increased snowmelt rates thereby advancing the date of snow disappearance compared with the unburned forest. We demonstrate a new parameterization of post‐fire snow albedo as a function of days‐since‐snowfall and net snowpack energy balance using an empirically based exponential decay function. Incorporating our new post‐fire snow albedo decay parameterization in a spatially distributed energy and mass balance snow model, we show significantly improved predictions of snow cover duration and spatial variability of snow water equivalent across the BF, particularly during the late snowmelt period. Field measurements, snow model results, and remote sensing data demonstrate that charred forests increase the radiative forcing to snow and advance the timing of snow disappearance for several years following fire. Copyright © 2016 John Wiley & Sons, Ltd.     

Modification and testing of a one‐dimensional energy and mass balance model for supraglacial snowpacks

A one‐dimensional energy and mass balance snow model (SNTHERM) has been modified for use with supraglacial snowpacks and applied to a point on Haut Glacier d'Arolla, Switzerland. It has been adapted to incorporate the underlying glacier ice and a site‐specific, empirically derived albedo routine. Model performance was tested against continuous measurements of snow depth and meltwater outflow from the base of the snowpack, and intermittent measurements of surface albedo and snowpack density profiles collected during the 1993 and 2000 melt seasons. Snow and ice ablation was simulated accurately. The timing of the daily pattern of meltwater outflow was well reproduced, although magnitudes were generally underestimated, possibly indicating preferential flow into the snowpack lysimeter. The model was used to assess the quantity of meltwater stored temporally within the unsaturated snowpack and meltwater percolation rates, which were found to be in agreement with dye tracer experiments undertaken on this glacier. As with other energy balance studies on alpine valley glaciers, the energy available for melt was dominated by net radiation (64%), with a sizable contribution from sensible heat flux (36%) and with a negligible latent heat flux overall, although there was more complex temporal variation on diurnal timescales. A basic sensitivity analysis indicated that melt rates were most sensitive to radiation, air temperature and snowpack density, indicating the need to accurately extrapolate/interpolate these variables when developing a spatially distributed framework for this model. Copyright © 2008 John Wiley & Sons, Ltd.     

Soil C and N response to changes in winter precipitation in a subalpine forest ecosystem,NW Italy

Among the potential effects of climate change on subalpine forest ecosystems during the winter season, the shift in snowline towards higher altitudes and the increase in frequency of rain events on the snowpack are of particular interest. Here, we present the results of a 2‐year field experiment conducted in a forest stand (Larix decidua) in NW Italy at 2020 m a.s.l. From 2009 to 2011, we monitored soil physical characteristics (temperature and moisture), and soil and soil solution chemistry, in particular carbon (C) and nitrogen (N) forms and their change in time, as affected by simulated late snowpack accumulation and rain on snow events. Late snowpack accumulation determined a stronger effect on soil thermal and moisture regimes than rain on snow events. Also soil chemistry was significantly affected by late snowfall simulation. Although microbial biomass C and N were not reduced by soil freezing, soil contents of the more labile dissolved organic carbon and inorganic N increased when the soil was affected by mild/hard freezing. Variations in the soil solution were shifted with respect to those observed in soil, with an increase in N‐NO₃^? concentrations occurring during spring and summer. This study highlights the potential N loss in subalpine soils under changing environmental conditions driven by a changing climate. Copyright © 2013 John Wiley & Sons, Ltd.     

Estimating sublimation of intercepted and sub‐canopy snow using eddy covariance systems

Direct measurements of winter water loss due to sublimation were made in a sub‐alpine forest in the Rocky Mountains of Colorado. Above‐and below‐canopy eddy covariance systems indicated substantial losses of winter‐season snow accumulation in the form of snowpack (0·41 mm d^?1) and intercepted snow (0·71 mm d^?1) sublimation. The partitioning between these over and under story components of water loss was highly dependent on atmospheric conditions and near‐surface conditions at and below the snow/atmosphere interface. High above‐canopy sensible heat fluxes lead to strong temperature gradients between vegetation and the snow‐surface, driving substantial specific humidity gradients at the snow surface and high sublimation rates. Intercepted snowfall resulted in rapid response of above‐canopy latent heat fluxes, high within‐canopy sublimation rates (maximum = 3·7 mm d^?1), and diminished sub‐canopy snowpack sublimation. These results indicate that sublimation losses from the sub‐canopy snowpack are strongly dependent on the partitioning of sensible and latent heat fluxes in the canopy. This compels comprehensive studies of snow sublimation in forested regions that integrate sub‐canopy and over‐story processes. Copyright © 2007 John Wiley & Sons, Ltd.     

Incorporating topographic variability into a simple regional snowmelt model

General circulation models (GCMs), or stand‐alone models that are forced by the output from GCMs, are increasingly being used to simulate the interactions between snow cover, snowmelt, climate and water resources. The variation in snowpack extent, and hence albedo, through time in a cell is likely to be substantial, especially in mid‐latitude mountainous regions. As a consequence, the energy budget simulation by a GCM relies on a realistic representation of snowpack extent. Similarly, from a water resource perspective, the spatial extent of the pack is key in predicting meltwater discharges into rivers. In this paper a simple computationally efficient regional snow model has been developed, which is based on a degree‐day approach and simulates the fraction of the model domain covered by snow, the spatially averaged melt rate and the mean snowpack depth. Computational efficiency is achieved through a novel spatial averaging procedure, which relies on the assumptions that precipitation and temperature scale linearly with elevation and that the distribution of elevations in the domain can be modelled by a continuous function. The resulting spatially averaged model is compared with both observations of the duration of snow cover throughout Austria and with results from a distributed model based on the same underlying assumptions but applied at a fine spatial resolution. The new spatially averaged model successfully simulated the seasonal snow duration observations and reproduced the daily dynamics of snow cover extent, the spatially averaged melt rate and mean pack depth simulated by the distributed model. It, therefore, offers a computationally efficient and easily applied alternative to the current crop of regional snow models. Copyright © 2004 John Wiley & Sons, Ltd.     

Long‐term investigations of the snow cover in a subalpine semi‐forested catchment

To improve spring runoff forecasts from subalpine catchments, detailed spatial simulations of the snow cover in this landscape is obligatory. For more than 30 years, the Swiss Federal Research Institute WSL has been conducting extensive snow cover observations in the subalpine watershed Alptal (central Switzerland). This paper summarizes the conclusions from past snow studies in the Alptal valley and presents an analysis of 14 snow courses located at different exposures and altitudes, partly in open areas and partly in forest. The long‐term performance of a physically based numerical snow–vegetation–atmosphere model (COUP) was tested with these snow‐course measurements. One single parameter set with meteorological input variables corrected to the prevailing local conditions resulted in a convincing snow water equivalent (SWE) simulation at most sites and for various winters with a wide range of snow conditions. The snow interception approach used in this study was able to explain the forest effect on the SWE as observed on paired snow courses. Finally, we demonstrated for a meadow and a forest site that a successful simulation of the snowpack yields appropriate melt rates. Copyright © 2006 John Wiley & Sons, Ltd.     

Climate variability,snow, and physiographic controls on storm hydrographs in small forested basins,western Cascades,Oregon

Large floods are often attributed to the melting of snow during a rain event. This study tested how climate variability, snowpack presence, and basin physiography were related to storm hydrograph shape in three small (<1 km²) basins with old‐growth forest in western Oregon. Relationships between hydrograph characteristics and precipitation were tested for approximately 800 storms over a nearly 30‐year period. Analyses controlled for (1) snowpack presence/absence, (2) antecedent soil moisture, and (3) hillslope length and gradient. For small storms (<150 mm precipitation), controlling for precipitation, the presence of a snowpack on near‐saturated soil increased the threshold of precipitation before hydrograph rise, extended the start lag, centroid lag, and duration of storm hydrographs, and increased the peak discharge. The presence of a snowpack on near‐saturated soil sped up and steepened storm hydrographs in a basin with short steep slopes, but delayed storm hydrographs in basins with longer or more gentle slopes. Hydrographs of the largest events, which were extreme regional rain and rain‐on‐snow floods, were not sensitive to landform characteristics or snowpack presence/absence. Although the presence of a snowpack did not increase peak discharge in small, forested basins during large storms, it had contrasting effects on storm timing in small basins, potentially synchronizing small basin contributions to the larger basin hydrograph during large rain‐on‐snow events. By altering the relative timing of hydrographs, snowpack melting could produce extreme floods from precipitation events whose size is not extreme. Further work is needed to examine effects of canopy openings, snowpack, and climate warming on extreme rain‐on‐snow floods at the large basin scale. Copyright © 2008 John Wiley & Sons, Ltd.     

Tracer release in melting snow: Diurnal and seasonal patterns

Time sequences of tracer release from an alpine snowpack were investigated at Mammoth Mountain, California in 1989. Lysimeter discharge and conductivity were recorded at 30 minute intervals. Three separate applications of chemical tracers were added to the snow surface to provide an ionic signal with known origins in the snowpack. Grab samples of meltwater and snow from snow pits were analysed for chemical composition. There were three distinct discharge periods, each characterized by diurnal fluctuations in discharge and conductivity. An inverse relation between discharge and conductivity was interpreted as the combination of a concentrated signal from regions in the pack less subject to leaching and a relatively dilute signal from near the snow surface where the snow was actively melting Conductivity peaks were highest and diurnal changes greatest immediately following periods of freezing. Grab samples showed little correlation with either 30 minute or daily average conductivity. Relative concentrations of individual ions in meltwater were similar between samples. Non-systematic grab sampling of snowpack meltwater is shown to be potentially misleading because of multiple ionic pulses over the ablation season and strong diurnal fluctuations in chemical concentrations. Continuous measurements of discharge conductivity are a good indicator of diurnal and seasonal changes in the rate of ion release from the snowpack, and should be used to guide sampling. Composite, or time-integrated samples rather than grab samples may be required to estimate daily and weekly rates of ion release in melting snow.     

Seasonal snowpack dynamics and runoff in a cool temperate forest: lysimeter experiment in Niigata,Japan

Seasonal snowpack dynamics are described through field measurements under contrasting canopy conditions for a mountainous catchment in the Japan Sea region. Microclimatic data, snow accumulation, albedo and lysimeter runoff are given through the complete winter season 2002–03 in (1) a mature cedar stand, (2) a larch stand, and (3) a regenerating cedar stand or opening. The accumulation and melt of seasonal snowpack strongly influences streamflow runoff during December to May, including winter baseflow, mid‐winter melt, rain on snow, and diurnal peaks driven by radiation melt in spring. Lysimeter runoff at all sites is characterized by constant ground melt of 0·8–1·0 mm day⁻¹. Rapid response to mid‐winter melt or rainfall shows that the snowpack remains in a ripe or near‐ripe condition throughout the snow‐cover season. Hourly and daily lysimeter discharge was greatest during rain on snow (e.g. 7 mm h⁻¹ and 53 mm day⁻¹ on 17 December) with the majority of runoff due to rainfall passing through the snowpack as opposed to snowmelt. For both rain‐on‐snow and radiation melt events lysimeter discharge was generally greatest at the open site, although there were exceptions such as during interception melt events. During radiation melt instantaneous discharge was up to 4·0 times greater in the opening compared with the mature cedar, and 48 h discharge was up to 2·5 times greater. Perhaps characteristic of maritime climates, forest interception melt is shown to be important in addition to sublimation in reducing snow accumulation beneath dense canopies. While sublimation represents a loss from the catchment water balance, interception melt percolates through the snowpack and contributes to soil moisture during the winter season. Strong differences in microclimate and snowpack albedo persisted between cedar, larch and open sites, and it is suggested further work is needed to account for this in hydrological simulation models. Copyright © 2005 John Wiley & Sons, Ltd.     

Influence of melt‐freeze‐cycles on the radionuclide transport in homogeneous laboratory snowpack

Radionuclides released to the environment and deposited with or onto snow can be stored over long time periods if ambient temperature stays low, particularly in glaciated areas or high alpine sites. The radionuclides will be accumulated in the snowpack during the winter unless meltwater runoff at the snow base occurs. They will be released to surface waters within short time during snowmelt in spring. In two experiments under controlled melting conditions of snow in the laboratory, radionuclide migration and runoff during melt‐freeze‐cycles were examined. The distribution of Cs‐134 and Sr‐85 tracers in homogeneous snow columns and their fractionation and potential preferential elution in the first meltwater portions were determined. Transport was associated with the percolation of meltwater at ambient temperatures above 0 °C after the snowpack became ripe. Mean migration velocities in the pack were examined for both nuclides to about 0.5 cm hr^?1 after one diurnal melt‐freeze‐cycle at ambient temperatures of ?2 to 4 °C. Meltwater fluxes were calculated with a median of 1.68 cm hr^?1. Highly contaminated portions of meltwater with concentration factors between 5 and 10 against initial bulk concentrations in the snowpack were released as ionic pulse with the first meltwater. Neither for caesium nor strontium preferential elution was observed. After recurrent simulated day‐night‐cycles (?2 to 4 °C), 80% of both radionuclides was released with the first 20% of snowmelt within 4 days. 50% of Cs‐134 and Sr‐85 were already set free after 24 hr. Snowmelt contained highest specific activities when the melt rate was lowest during the freeze‐cycles due to concentration processes in remaining liquids, enhanced by the melt‐freeze‐cycling. This implies for natural snowpack after significant radionuclide releases, that long‐time accumulation of radionuclides in the snow during frost periods, followed by an onset of steady meltwater runoff at low melt rates, will cause the most pronounced removal of the contaminants from the snow cover. This scenario represents the worst case of impact on water quality and radiation exposure in aquatic environments.     

A modelling framework to simulate field‐scale nitrate response and transport during snowmelt: The WINTRA model

Modelling nutrient transport during snowmelt in cold regions remains a major scientific challenge. A key limitation of existing nutrient models for application in cold regions is the inadequate representation of snowmelt, including hydrological and biogeochemical processes. This brief period can account for more than 80% of the total annual surface runoff in the Canadian Prairies and Northern Canada and processes such as atmospheric deposition, overwinter redistribution of snow, ion exclusion from snow crystals, frozen soils, and snow‐covered area depletion during melt influence the distribution and release of snow and soil nutrients, thus affecting the timing and magnitude of snowmelt runoff nutrient concentrations. Research in cold regions suggests that nitrate (NO₃) runoff at the field‐scale can be divided into 5 phases during snowmelt. In the first phase, water and ions originating from ion‐rich snow layers travel and diffuse through the snowpack. This process causes ion concentrations in runoff to gradually increase. The second phase occurs when this snow ion meltwater front has reached the bottom of the snowpack and forms runoff to the edge‐of‐the‐field. During the third and fourth phases, the main source of NO₃ transitions from the snowpack to the soil. Finally, the fifth and last phase occurs when the snow has completely melted, and the thawing soil becomes the main source of NO₃ to the stream. In this research, a process‐based model was developed to simulate hourly export based on this 5‐phase approach. Results from an application in the Red River Basin of southern Manitoba, Canada, shows that the model can adequately capture the dynamics and rapid changes of NO₃ concentrations during this period at relevant temporal resolutions. This is a significant achievement to advance the current nutrient modelling paradigm in cold climates, which is generally limited to satisfactory results at monthly or annual resolutions. The approach can inform catchment‐scale nutrient models to improve simulation of this critical snowmelt period.     

Snowfall interception in a deciduous Nothofagus forest and implications for spatial snowpack distribution

Native Nothofagus forests in the midlatitude region of the Andes Cordillera are notorious biodiversity hot spots, uniquely situated in the Southern Hemisphere such that they develop in snow‐dominated reaches of this mountain range. Spanning a smaller surface area than similar ecosystems, where forests and snow coexist in the Northern Hemisphere, the interaction between vegetation and snow processes in this ecotone has received lesser attention. We present the first systematic study of snow–vegetation interactions in the Nothofagus forests of the Southern Andes, focusing on how the interplay between interception and climate determines patterns of snow water equivalent (SWE) variability. The Valle Hermoso experimental catchment, located in the Nevados de Chillán vicinity, was fitted with eight snow depth sensors that provided continuous measurements at varying elevations, aspect, and forest cover. Also, manual measurements of snow properties were obtained during snow surveys conducted during end of winter and spring seasons for 3 years, between 2015 and 2017. Each year was characterized by distinct climatological conditions, with 2016 representing one of the driest winters on record in this region. Distance to canopy, leaf area index, and total gap area were measured at each observational site. A regression model was built on the basis of statistical analysis of local parameters to model snow interception in this kind of forest. We find that interception implied a 23.2% reduction in snow accumulation in forested sites compared with clearings. The interception in these deciduous trees represents, on average, 23.6% of total annual snowfall, reaching a maximum measured interception value of 13.8‐mm SWE for all snowfall events analysed in this research.     

Subgrid variability of snow water equivalent at operational snow stations in the western USA

The spatial distribution of snow water equivalent (SWE) is a key variable in many regional‐scale land surface models. Currently, the assimilation of point‐scale snow sensor data into these models is commonly performed without consideration of the spatial representativeness of the point data with respect to the model grid‐scale SWE. To improve the understanding of the relationship between point‐scale snow measurements and surrounding areas, we characterized the spatial distribution of snow depth and SWE within 1‐, 4‐ and 16‐km² grids surrounding 15 snow stations (snowpack telemetry and California snow sensors) in California, Colorado, Wyoming, Idaho and Oregon during the 2008 and 2009 snow seasons. More than 30 000 field observations of snowpack properties were used with binary regression tree models to relate SWE at the sensor site to the surrounding area SWE to evaluate the sensor representativeness of larger‐scale conditions. Unlike previous research, we did not find consistent high biases in snow sensor depth values as biases over all sites ranged from 74% overestimates to 77% underestimates. Of the 53 assessments, 27 surveys indicated snow station biases of less than 10% of the surrounding mean observed snow depth. Depth biases were largely dictated by the physiographic relationship between the snow sensor locations and the mean characteristics of the surrounding grid, in particular, elevation, solar radiation index and vegetation density. These scaling relationships may improve snow sensor data assimilation; an example application is illustrated for the National Operational Hydrologic Remote Sensing Center National Snow Analysis SWE product. The snow sensor bias information indicated that the assimilation of point data into the National Operational Hydrologic Remote Sensing Center model was often unnecessary and reduced model accuracy. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling snow water equivalent and spring runoff in a boreal watershed,James Bay,Canada

The hydrology of boreal regions is strongly influenced by seasonal snow accumulation and melt. In this study, we compare simulations of snow water equivalent (SWE) and streamflow by using the hydrological model HYDROTEL with two contrasting approaches for snow modelling: a mixed degree‐day/energy balance model (small number of inputs, but several calibration parameters needed) and the thermodynamic model CROCUS (large number of inputs, but no calibration parameter needed). The study site, in Northern Quebec, Canada was equipped with a ground‐based gamma ray sensor measuring the SWE continuously for 5 years in a small forest clearing. The first simulation of CROCUS showed a tendency to underestimate SWE, attributable to bias in the meteorological inputs. We found that it was appropriate to use a threshold of 2 °C to separate rain and snow. We also applied a correction to account for snowfall undercatch by the precipitation gauge. After these modifications to the input dataset, we noticed that CROCUS clearly overestimated the SWE, likely as a result of not including loss in SWE because of blowing snow sublimation and relocation. To correct this, we included into CROCUS a simple parameterisation effective after a certain wind speed threshold, after which the thermodynamic model performed much better than the traditional mixed degree‐day/energy balance model. HYDROTEL was then used to simulate streamflow with both snow models. With CROCUS, the main peak flow could be captured, but the second peak because of delayed snowmelt from forested areas could not be reproduced due to a lack of sub‐canopy radiation data to feed CROCUS. Despite the relative homogeneity of the boreal landscape, data inputs from each land cover type are needed to generate satisfying simulation of the spring runoff. Copyright © 2013 John Wiley & Sons, Ltd.     

Optimal design of snow stake networks to estimate snow depth in an alpine mountain range

Monitoring and estimation of snow depth in alpine catchments is needed for a proper assessment of management alternatives for water supply in these water resources systems. The distribution of snowpack thickness is usually approached by using field data that come from snow samples collected at a given number of locations that constitute the monitoring network. Optimal design of this network is required to obtain the best possible estimates. Assuming that there is an existing monitoring network, its optimization may imply the selection of an optimal network as a subset of the existing one (if there are no funds to maintain them) or enlarging the existing network by one or more stations (optimal augmentation problem). We propose an optimization procedure that minimizes the total variance in the estimate of snowpack thickness. The novelty of this work is to treat, for the first time, the problem of snow observation network optimization for an entire mountain range rather than for small catchments as done in the previous studies. Taking into account the reduced data available, which is a common problem in many mountain ranges, the importance of a proper design of these observation networks is even larger. Snowpack thickness is estimated by combining regression models to approach the effect of the explanatory variables and kriging techniques to consider the influence of the stakes location. We solve the optimization problems under different hypotheses, studying the impacts of augmentation and reduction, both, one by one and in pairs. We also analyse the sensitivity of results to nonsnow measurements deduced from satellite information. Finally, we design a new optimal network by combining the reduction and augmentation methods. The methodology has been applied to the Sierra Nevada mountain range (southern Spain), where very limited resources are employed to monitor snowfall and where an optimal snow network design could prove critical. An optimal snow observation network is defined by relocating some observation points. It would reduce the estimation variance by around 600 cm² (15%).     

Intercomparison of measurements of bulk snow density and water equivalent of snow cover with snow core samplers: Instrumental bias and variability induced by observers

Manually collected snow data are often considered as ground truth for many applications such as climatological or hydrological studies. However, there are many sources of uncertainty that are not quantified in detail. For the determination of water equivalent of snow cover (SWE), different snow core samplers and scales are used, but they are all based on the same measurement principle. We conducted two field campaigns with 9 samplers commonly used in observational measurements and research in Europe and northern America to better quantify uncertainties when measuring depth, density and SWE with core samplers. During the first campaign, as a first approach to distinguish snow variability measured at the plot and at the point scale, repeated measurements were taken along two 20 m long snow pits. The results revealed a much higher variability of SWE at the plot scale (resulting from both natural variability and instrumental bias) compared to repeated measurements at the same spot (resulting mostly from error induced by observers or very small scale variability of snow depth). The exceptionally homogeneous snowpack found in the second campaign permitted to almost neglect the natural variability of the snowpack properties and focus on the separation between instrumental bias and error induced by observers. Reported uncertainties refer to a shallow, homogeneous tundra-taiga snowpack less than 1 m deep (loose, mostly recrystallised snow and no wind impact). Under such measurement conditions, the uncertainty in bulk snow density estimation is about 5% for an individual instrument and is close to 10% among different instruments. Results confirmed that instrumental bias exceeded both the natural variability and the error induced by observers, even in the case when observers were not familiar with a given snow core sampler.     

Evaluation of temporal consistency of snow depth drivers of a Rocky Mountain watershed in southern Alberta

Collecting spatially representative data over large areas is a challenge within snow monitoring frameworks. Identifying consistent trends in snow accumulation properties enables increased sampling efficiency by minimizing field collection time and/or remote sensing costs. Seasonal snowpack depth estimations during mid-winter and melt onset conditions were derived from airborne Lidar over the West Castle Watershed in the southern Canadian Rockies on three dates. Each dataset was divided into five sets of snow depth driver classes: elevation, aspect, topographic position index, canopy cover and slope. Datasets were quality controlled by eliminating snow depth values above the 99th percentile value, which had a negligible effect on average snow depths. Consistent trends were observed among driver classes with peak snow accumulation occurring within the treeline ecotone, north-facing aspects, open canopies, topographic depressions and areas with low slope angle. Although mid-winter class trends for each driver were similar and watershed-scale snow depth distributions were significantly correlated (0.76, p < .01), depth distributions within the same driver class of the three datasets were not correlated due to recent snowfall events, redistribution and settling processes. Trends in driver classes during late season melt onset were similar to mid-winter conditions but watershed scale distribution correlation results varied with seasonality (0.68 mid-winter 2014 and melt onset 2016; 0.65 mid-winter 2017 and melt onset 2016, p < .1). This is due to the differing stages of accumulation or ablation and the upward migration in the 0°C isotherm during spring, when snow depth can be declining in valley bottoms while still increasing at higher elevations. The observed consistency in depth driver controls can be used to guide future integrated snow monitoring frameworks.