The cold rain‐on‐snow event of June 2013 in the Canadian Rockies — characteristics and diagnosis

The June 2013 flood in the Canadian Rockies featured rain‐on‐snow (ROS) runoff generation at alpine elevations that contributed to the high streamflows observed during the event. Such a mid‐summer ROS event has not been diagnosed in detail, and a diagnosis may help to understand future high discharge‐producing hydrometeorological events in mountainous cold regions. The alpine hydrology of the flood was simulated using a physically based model created with the modular cold regions hydrological modelling platform. The event was distinctive in that, although at first, relatively warm rain fell onto existing snowdrifts inducing ROS melt; the rainfall turned to snowfall as the air mass cooled and so increased snowcover and snowpacks in alpine regions, which then melted rapidly from ground heat fluxes in the latter part of the event. Melt rates of existing snowpacks were substantially lower during the ROS than during the relatively sunny periods preceding and following the event as a result of low wind speeds, cloud cover and cool temperatures. However, at the basin scale, melt volumes increased during the event as a result of increased snowcover from the fresh snowfall and consequent large ground heat contributions to melt energy, causing snowmelt to enhance rainfall–runoff by one fifth. Flow pathways also shifted during the event from relatively slow sub‐surface flow prior to the flood to an even contribution from sub‐surface and fast overland flow during and immediately after the event. This early summer, high precipitation ROS event was distinctive for the impact of decreased solar irradiance in suppressing melt rates, the contribution of ground heat flux to basin scale snowmelt after precipitation turned to snowfall, the transition from slow sub‐surface to fast overland flow runoff as the sub‐surface storage saturated and streamflow volumes that exceeded precipitation. These distinctions show that summer, mountain ROS events should be considered quite distinct from winter ROS and can be important contributors to catastrophic events. Copyright © 2016 John Wiley & Sons, Ltd.     

Spatio‐temporal controls of snowmelt and runoff generation during rain‐on‐snow events in a mid‐latitude mountain catchment

A network of 30 standalone snow monitoring stations was used to investigate the snow cover distribution, snowmelt dynamics, and runoff generation during two rain‐on‐snow (ROS) events in a 40 km² montane catchment in the Black Forest region of southwestern Germany. A multiple linear regression analysis using elevation, aspect, and land cover as predictors for the snow water equivalent (SWE) distribution within the catchment was applied on an hourly basis for two significant ROS flood events that occurred in December 2012. The available snowmelt water, liquid precipitation, as well as the total retention storage of the snow cover were considered in order to estimate the amount of water potentially available for the runoff generation. The study provides a spatially and temporally distributed picture of how the two observed ROS floods developed in the catchment. It became evident that the retention capacity of the snow cover is a crucial mechanism during ROS. It took several hours before water was released from the snowpack during the first ROS event, while retention storage was exceeded within 1 h from the start of the second event. Elevation was the most important terrain feature. South‐facing terrain contributed more water for runoff than north‐facing slopes, and only slightly more runoff was generated at open compared to forested areas. The results highlight the importance of snowmelt together with liquid precipitation for the generation of flood runoff during ROS and the large temporal and spatial variability of the relevant processes. Copyright © 2015 John Wiley & Sons, Ltd.     

Spatial and temporal diversity for debris‐flow meteorological control in subarctic oceanic periglacial environments in Iceland

Slopes in fjord environments of Iceland are prone to debris‐flow initiation, responding to a wide variety of meteorological triggering factors, such as rain on snow, rapid snowmelt, long‐lasting rainfall or intense rainfall. If all fjord regions have similar debris flows with regards to their magnitude, their meteorological control is diverse both in space and in time. Debris flows in Northwest Iceland are triggered mostly by rain‐on‐snow and long‐lasting rainfall, while snowmelt is more characteristic in North Iceland, and rainfall has a clear impact in East Iceland. Most debris‐flow events occur on a single slope, and only a few are recorded at the same time in different regions. Observations of the threshold values underline the diversity of debris‐flow initiation, occurring with huge amounts of sudden water supply as well as with very moderate ones. Copyright © 2007 John Wiley & Sons, Ltd.     

Rain on snow induced urban floods in cold maritime climate: Risk,indicators and trends

Urban floods pose a societal and economical risk. This study evaluated the risk and hydro-meteorological conditions that cause pluvial flooding in coastal cities in a cold climate. Twenty years of insurance claims data and up to 97 years of meteorological data were analysed for Reykjavík, Iceland (64.15°N; <100 m above sea level). One third of the city's wastewater collection system is combined, and pipe grades vary from 0.5% to 10%. Results highlight semi-intensive rain (<7 mm/h; ≤3 year return period) in conjunction with snow and frozen ground as the main cause for urban flood risk in a climate which undergoes frequent snow and frost cycles (avg. 13 and 19 per season, respectively). Floods in winter were more common, more severe and affected a greater number of neighbourhoods than during summer. High runoff volumes together with debris remobilized with high winds challenged the capacity of wastewater systems regardless of their age or type (combined vs. separate). The two key determinants for the number of insurance claims were antecedent frost depth and total precipitation volume per event. Two pluvial regimes were particularly problematic: long duration (13–25 h), late peaking rain on snow (RoS), where snowmelt enhanced the runoff intensity, elongated and connected independent rainfall into a singular, more voluminous (20–76 mm) event; shorter duration (7–9 h), more intensive precipitation that evolved from snow to rain. Closely timed RoS and cooling were believed to trigger frost formation. A positive trend was detected in the average seasonal snow depth and volume of rain and snowmelt during RoS events. More emphasis, therefore, needs to be placed on designing and operating urban drainage infrastructure with regard to RoS co-acting with frozen ground. Furthermore, more detailed, routine monitoring of snow and soil conditions is important to predict RoS flood events.     

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.     

Aspect control of water movement on hillslopes near the rain–snow transition of the Colorado Front Range

In the Colorado Front Range, forested catchments near the rain–snow transition are likely to experience changes in snowmelt delivery and subsurface water transport with climate warming and associated shifts in precipitation patterns. Snowpack dynamics are strongly affected by aspect: Lodgepole pine forested north‐facing slopes develop a seasonal snowpack, whereas Ponderosa pine‐dotted south‐facing slopes experience intermittent snow accumulation throughout winter and spring. We tested the degree to which these contrasting water input patterns cause different near‐surface hydrologic response on north‐facing and south‐facing hillslopes during the snowmelt period. During spring snowmelt, we applied lithium bromide (LiBr) tracer to instrumented plots along a north–south catchment transect. Bromide broke through immediately at 10‐ and 30‐cm depths on the north‐facing slope and was transported out of soil waters within 40 days. On the south‐facing slope, Br^? was transported to significant depths only during spring storms and remained above the detection limit throughout the study. Modelling of unsaturated zone hydrologic response using Hydrus‐1D corroborated these aspect‐driven differences in subsurface transport. Our multiple lines of evidence suggest that north‐facing slopes are dominated by connected flow through the soil matrix, whereas south‐facing slope soils experience brief periods of rapid vertical transport following snowmelt events and are drier overall than north‐facing slopes. These differences in hydrologic response were largely a function of energy‐driven differences in water supply, emphasizing the importance of aspect and climate forcing when considering contributions of water and solutes to streamflow in catchments near the snow line. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Testing a spatially distributed tracer‐aided runoff model in a snow‐influenced catchment: Effects of multicriteria calibration on streamwater ages

Integrating stable isotope tracers into rainfall‐runoff models allows investigation of water partitioning and direct estimation of travel times and water ages. Tracer data have valuable information content that can be used to constrain models and, in integration with hydrometric observations, test the conceptualization of catchment processes in model structure and parameterization. There is great potential in using tracer‐aided modelling in snow‐influenced catchments to improve understanding of these catchments' dynamics and sensitivity to environmental change. We used the spatially distributed tracer‐aided rainfall‐runoff (STARR) model to simulate the interactions between water storage, flux, and isotope dynamics in a snow‐influenced, long‐term monitored catchment in Ontario, Canada. Multiple realizations of the model were achieved using a combination of single and multiple objectives as calibration targets. Although good simulations of hydrometric targets such as discharge and snow water equivalent could be achieved by local calibration alone, adequate capture of the stream isotope dynamics was predicated on the inclusion of isotope data in the calibration. Parameter sensitivity was highest, and most local, for single calibration targets. With multiple calibration targets, key sensitive parameters were still identifiable in snow and runoff generation routines. Water ages derived from flux tracking subroutines in the model indicated a catchment where runoff is dominated by younger waters, particularly during spring snowmelt. However, resulting water ages were most sensitive to the partitioning of runoff sources from soil and groundwater sources, which was most realistically achieved when isotopes were included in the calibration. Given the paucity of studies where hydrological models explicitly incorporate tracers in snow‐influenced regions, this study using STARR is an important contribution to satisfactorily simulating snowpack dynamics and runoff generation processes, while simultaneously capturing stable isotope variability in snow‐influenced catchments.     

The June 2013 Alberta Catastrophic Flooding Event: Part 1—Climatological aspects and hydrometeorological features

In June 2013, excessive rainfall associated with an intense weather system triggered severe flooding in southern Alberta, which became the costliest natural disaster in Canadian history. This article provides an overview of the climatological aspects and large‐scale hydrometeorological features associated with the flooding event based upon information from a variety of sources, including satellite data, upper air soundings, surface observations and operational model analyses. The results show that multiple factors combined to create this unusually severe event. The event was characterized by a slow‐moving upper level low pressure system west of Alberta, blocked by an upper level ridge, while an associated well‐organized surface low pressure system kept southern Alberta, especially the eastern slopes of the Rocky Mountains, in continuous precipitation for up to two days. Results from air parcel trajectory analysis show that a significant amount of the moisture originated from the central Great Plains, transported into Alberta by a southeasterly low level jet. The event was first dominated by significant thunderstorm activity, and then evolved into continuous precipitation supported by the synoptic‐scale low pressure system. Both the thunderstorm activity and upslope winds associated with the low pressure system produced large rainfall amounts. A comparison with previous similar events occurring in the same region suggests that the synoptic‐scale features associated with the 2013 rainfall event were not particularly intense; however, its storm environment was the most convectively unstable. The system also exhibited a relatively high freezing level, which resulted in rain, rather than snow, mainly falling over the still snow‐covered mountainous areas. Melting associated with this rain‐on‐snow scenario likely contributed to downstream flooding. Furthermore, above‐normal snowfall in the preceding spring helped to maintain snow in the high‐elevation areas, which facilitated the rain‐on‐snow event. Copyright © 2016 John Wiley & Sons, Ltd.     

A comparison of snowmelt‐derived streamflow from temperature‐index and modified‐temperature‐index snow models

Reliable estimation of the volume and timing of snowmelt runoff is vital for water supply and flood forecasting in snow‐dominated regions. Snowmelt is often simulated using temperature‐index (TI) models due to their applicability in data‐sparse environments. Previous research has shown that a modified‐TI model, which uses a radiation‐derived proxy temperature instead of air temperature as its surrogate for available energy, can produce more accurate snow‐covered area (SCA) maps than a traditional TI model. However, it is unclear whether the improved SCA maps are associated with improved snow water equivalent (SWE) estimation across the watershed or improved snowmelt‐derived streamflow simulation. This paper evaluates whether a modified‐TI model produces better streamflow estimates than a TI model when they are used within a fully distributed hydrologic model. It further evaluates the performance of the two models when they are calibrated using either point SWE measurements or SCA maps. The Senator Beck Basin in Colorado is used as the study site because its surface is largely bedrock, which reduces the role of infiltration and emphasizes the role of the SWE pattern on streamflow generation. Streamflow is simulated using both models for 6 years. The modified‐TI model produces more accurate streamflow estimates (including flow volume and peak flow rate) than the TI model, likely because the modified‐TI model better reproduces the SWE pattern across the watershed. Both models also produce better performance when calibrated with SCA maps instead of point SWE data, likely because the SCA maps better constrain the space‐time pattern of SWE.     

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.     

Identifying run‐off contributions during melt‐induced run‐off events in a glacierized alpine catchment

We analysed contributions to run‐off using hourly stream water samples from seven individual melt‐induced run‐off events (plus one rainfall event) during 2011, 2012 and 2013 in two nested glacierized catchments in the Eastern Italian Alps. Electrical conductivity and stable isotopes of water were used for mixing analysis and two‐component and three‐component hydrograph separation. High‐elevation snowmelt, glacier melt and autumn groundwater were identified as major end‐members. Discharge and tracers in the stream followed the diurnal variations of air temperature but markedly reacted to rainfall inputs. Hysteresis patterns between discharge and electrical conductivity during the melt‐induced run‐off events revealed contrasting loop directions at the two monitored stream sections. Snowmelt contribution to run‐off was highest in June and July (up to 33%), whereas the maximum contribution of glacier melt was reached in August (up to 65%). The maximum groundwater and rainfall contributions were 62% and 11%, respectively. Run‐off events were generally characterized by decreasing snowmelt and increasing glacier melt fractions from the beginning to the end of the summer 2012, while run‐off events in 2013 showed less variable snowmelt and lower glacier melt contributions than in 2012. The results provided essential insights into the complex dynamics of melt‐induced run‐off events and may be of further use in the context of water resource management in alpine catchments. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaluating the SWAT's snow hydrology over a Northern Quebec watershed

The snow treatment becomes an important component of Soil and Water Assessment Tool (SWAT)’s hydrology when spring flows are dominated by snow melting. However, little is known about SWAT's snow hydrology performance because most studies using SWAT were conducted in rainfall‐driven catchments. To fill this gap, the present study aims to evaluate the ability of SWAT in simulating snow‐melting‐dominated streamflow in the Outardes Basin in Northern Quebec. SWAT performance in simulating snowmelt is evaluated against observed streamflow data and compared to simulations from the operationally used Streamflow Synthesis and Reservoir Regulation (SSARR) model over that catchment. The SWAT 5‐year calibration showed a satisfactory performance at the daily and seasonal time scales with low volume biases. The SWAT validation was conducted over two (17‐year and 15‐year) periods. Performances were similar to the calibration period in simulating the daily and seasonal streamflows again with low model biases. The spring‐snowmelt‐generated peak flow was accurately simulated by SWAT both in magnitude and timing. When SWAT's results are compared to SSARR, similar performances in simulating the daily discharges were observed. SSARR simulates more accurately streamflow generated at the snowmelt onset whereas SWAT better predicts streamflow in summer, fall and winter. SWAT provided reasonable streamflow simulations for our snow‐covered catchment, but refinement of the process‐driven baseflow during the snowmelt onset could improve spring performances. Therefore, SWAT becomes an attractive tool for evaluating water resources management in Nordic environments when a distributed model is preferred or when water quality information (e.g. temperature) is required. Copyright © 2013 John Wiley & Sons, Ltd.     

Changes in seasonality and timing of peak streamflow in snow and semi‐arid climates of the north‐central United States, 1910–2012

Changes in the seasonality and timing of annual peak streamflow in the north‐central USA are likely because of changes in precipitation and temperature regimes. A source of long‐term information about flood events across the study area is the U.S. Geological Survey peak streamflow database. However, one challenge of answering climate‐related questions with this dataset is that even in snowmelt‐dominated areas, it is a mixed population of snowmelt/spring rain generated peaks and summer/fall rain generated peaks. Therefore, a process was developed to divide the annual peaks into two populations, or seasons, snowmelt/spring, and summer/fall. The two series were then tested for the hypotheses that because of changes in precipitation regimes, the odds of summer/fall peaks have increased and, because of temperature changes, snowmelt/spring peaks happen earlier. Over climatologically and geographically similar regions in the north‐central USA, logistic regression was used to model the odds of getting a summer/fall peak. When controlling for antecedent wet and dry conditions and geographical differences, the odds of summer/fall peaks occurring have increased across the study area. With respect to timing within the seasons, trend analysis showed that in northern portions of the study region, snowmelt/spring peaks are occurring earlier. The timing of snowmelt/spring peaks in three regions in the northern part of the study area is earlier by 8.7– 14.3 days. These changes have implications for water interests, such as potential changes in lead‐time for flood forecasting or changes in the operation of flood‐control dams. Copyright © 2015 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.     

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.     

On the role of spatial resolution on snow estimates using a process‐based snow model across a range of climatology and elevation

Hydrological processes in mountainous settings depend on snow distribution, whose prediction accuracy is a function of model spatial scale. Although model accuracy is expected to improve with finer spatial resolution, an increase in resolution comes with modelling costs related to increased computational time and greater input data and parameter information. This computational and data collection expense is still a limiting factor for many large watersheds. Thus, this work's main objective is to question which physical processes lead to loss in model accuracy with regard to input spatial resolution under different climatic conditions and elevation ranges. To address this objective, a spatially distributed snow model, iSnobal, was run with inputs distributed at 50‐m—our benchmark for comparison—and 100‐m resolutions and with aggregated (averaged from the fine to the large resolution) inputs from the 50‐m model to 100‐, 250‐, 500‐, and 750‐m resolution for wet, average, and dry years over the Upper Boise River Basin (6,963 km²), which spans four elevation bands: rain dominated, rain–snow transition, and snow dominated below treeline and above treeline. Residuals, defined as differences between values quantified with high resolution (>50 m) models minus the benchmark model (50 m), of simulated snow‐covered area (SCA) and snow water equivalent (SWE) were generally slight in the aggregated scenarios. This was due to transferring the effects of topography on meteorological variables from the 50‐m model to the coarser scales through aggregation. Residuals in SCA and SWE in the distributed 100‐m simulation were greater than those of the aggregated 750 m. Topographic features such as slope and aspect were simplified, and their gradient was reduced due to coarsening the topography from the 50‐ to 100‐m resolution. Therefore, solar radiation was overestimated, and snow drifting was modified and caused substantial SCA and SWE underestimation in the distributed 100‐m model relative to the 50‐m model. Large residuals were observed in the wet year and at the highest elevation band when and where snow mass was large. These results support that model accuracy is substantially reduced with model scales coarser than 50 m.     

Controls of streamflow generation in small catchments across the snow–rain transition in the Southern Sierra Nevada,California

Processes controlling streamflow generation were determined using geochemical tracers for water years 2004–2007 at eight headwater catchments at the Kings River Experimental Watersheds in southern Sierra Nevada. Four catchments are snow‐dominated, and four receive a mix of rain and snow. Results of diagnostic tools of mixing models indicate that Ca²⁺, Mg²⁺, K⁺ and Cl^? behaved conservatively in the streamflow at all catchments, reflecting mixing of three endmembers. Using endmember mixing analysis, the endmembers were determined to be snowmelt runoff (including rain on snow), subsurface flow and fall storm runoff. In seven of the eight catchments, streamflow was dominated by subsurface flow, with an average relative contribution (% of streamflow discharge) greater than 60%. Snowmelt runoff contributed less than 40%, and fall storm runoff less than 7% on average. Streamflow peaked 2–4 weeks earlier at mixed rain–snow than snow‐dominated catchments, but relative endmember contributions were not significantly different between the two groups of catchments. Both soil water in the unsaturated zone and regional groundwater were not significant contributors to streamflow. The contributions of snowmelt runoff and subsurface flow, when expressed as discharge, were linearly correlated with streamflow discharge (R² of 0.85–0.99). These results suggest that subsurface flow is generated from the soil–bedrock interface through preferential pathways and is not very sensitive to snow–rain proportions. Thus, a declining of the snow–rain ratio under a warming climate should not systematically affect the processes controlling the streamflow generation at these catchments. Copyright © 2012 John Wiley & Sons, Ltd.     

Measuring snow accumulation and ablation dynamics during rain-on-snow events: innovative measurement techniques

Rain-on-snow (ROS) is the primary generator of peak flow events in mountainous coastal regions of North America. Uncertainty remains as to the role of forest canopy interception leading up to and during ROS events. Much of this uncertainty can be attributed to a lack of suitable techniques to collect data during ROS, due in part to the dynamic nature of climatic conditions, particularly related to snow accumulation and melt. We supplemented a meteorological network with non-weighing snow melt lysimeters, suspended spring scales to measure snow throughfall and an automated time lapse photography network to monitor state of precipitation (rain vs. snow), snow accumulation/ablation, canopy interception and unloading of snow from the canopy. Image analysis software allowed for the extraction of data from images. Rapid loading and unloading of snow from the canopy, closely linked to changes in temperature, was observed using this approach. We were also able to continuously monitor throughfall snow water equivalent using low cost suspended spring scales. This experimental design allowed us to capture information previously unavailable without direct observation. Copyright © 2008 John Wiley & Sons, Ltd.     

Modified passive capillary samplers for collecting samples of snowmelt infiltration for stable isotope analysis in remote,seasonally inaccessible watersheds 2: field evaluation

Twelve modified passive capillary samplers (M‐PCAPS) were installed in remote locations within a large, alpine watershed located in the southern Rocky Mountains of Colorado to collect samples of infiltration during the snowmelt and summer rainfall seasons. These samples were collected in order to provide better constraints on the isotopic composition of soil‐water endmembers in the watershed. The seasonally integrated stable isotope composition (δ¹⁸O and δ²H) of soil‐meltwater collected with M‐PCAPS installed at shallow soil depths < 10 cm was similar to the seasonally integrated isotopic composition of bulk snow taken at the soil surface. However, meltwater which infiltrated to depths > 20 cm evolved along an isotopic enrichment line similar to the trendline described by the evolution of fresh snow to surface runoff from snowmelt in the watershed. Coincident changes in geochemistry were also observed at depth suggesting that the isotopic and geochemical composition of deep infiltration may be very different from that obtained by surface and/or shallow‐subsurface measurements. The M‐PCAPS design was also used to estimate downward fluxes of meltwater during the snowmelt season. Shallow and deep infiltration averaged 8·4 and 4·7 cm of event water or 54 and 33% of the measured snow water equivalent (SWE), respectively. Finally, dominant shallow‐subsurface runoff processes occurring during snowmelt could be identified using geochemical data obtained with the M‐PCAPS design. One soil regime was dominated by a combination of slow matrix flow in the shallow soil profile and fast preferential flow at depth through a layer of platy, volcanic rocks. The other soil regime lacked the rock layer and was dominated by slow matrix flow. Based on these results, the M‐PCAPS design appears to be a useful, robust methodology to quantify soil‐water fluxes during the snowmelt season and to sample the stable isotopic and geochemical composition of soil‐meltwater endmembers in remote watersheds. Copyright © 2009 John Wiley & Sons, Ltd.