Cold Regions Hydrological Modelling Platform Research Articles

Simulation of the impact of future changes in climate on the hydrology of Bow River headwater basins in the Canadian Rockies

This study diagnoses the impact of projected changes in climate and glacier cover on the hydrology of several natural flowing Bow River headwater basins in the Canadian Rockies: the Bow River at Lake Louise (∼420.7 km2), the Pipestone River near Lake Louise (∼304.2 km2), the Bow River at Banff (∼2192.2 km2) all of which drain the high elevation, snowy, partially glaciated Central Range, and the Elbow River at Calgary (∼1191.9 km2), which drains the drier Front Ranges and foothills, using models created using the modular, flexible, physically based Cold Regions Hydrological Modelling platform (CRHM). Hydrological models were constructed and parameterised in CRHM from local research results to include relevant streamflow generation processes for Canadian Rockies headwater basins, such as blowing snow, avalanching, snow interception and sublimation, energy budget snow and glacier melt, infiltration to frozen and unfrozen soils, hillslope sub-surface water redistribution, wetlands, lakes, evapotranspiration, groundwater flow, surface runoff and open channel flow. Surface layer outputs from Weather Research and Forecasting (WRF) model simulations for the current climate and for the late 21st century climate under a “business-as-usual” scenario, Representative Concentration Pathway 8.5 (RCP8.5) at 4-km resolution, were used to force model simulations to examine the climate change impact. A projected glacier cover under a “business-as-usual” scenario (RCP8.5) was incorporated to assess the impact of concomitant glacier cover decline. Uncalibrated model simulations for the current climate and glacier coverage showed useful predictions of snow accumulation, snowmelt, and streamflow when compared to surface observations from 2000 to 2015. Under the RCP8.5 climate change scenario, the basins of the Bow River at Banff and Elbow River at Calgary will warm up by 4.7 and 4.5 °C respectively and receive 12% to 15% more precipitation annually, with both basins experiencing a greater proportion of precipitation as rainfall. Peak snow accumulation in Bow River Basin will slightly rise by 3 mm, whilst it will drop by 20 mm in Elbow River Basin, and annual snowmelt volume will increase by 43 mm in Bow River Basin but decrease by 55 mm in Elbow River Basin. Snowcovered periods will decline by 37 and 46 days in Bow and Elbow river basins respectively due to suppressed snow redistribution by wind and gravity and earlier melt. The shorter snowcovered period and warmer, wetter climate will increase evapotranspiration and glacier melt, if the glaciers were held constant, and decrease sublimation, lake levels, soil moisture and groundwater levels. The hydrological responses of the basins will differ despite similar climate changes because of differing biophysical characteristics, climates and hydrological processes generating runoff. Climate change with concomitant glacier decline is predicted to increase the peak discharge and mean annual water yield by 12.23 m3 s−1 (+11%) and 11% in the higher elevation basins of the Bow River but will decrease the mean annual peak discharge by 3.58 m3 s−1 (-9%) and increase the mean annual water yield by 18% in the lower elevation basin of the Elbow River. This shows complex and compensatory hydrological process responses to climate change with the reduced glacier contribution reducing the impact of higher precipitation in high elevation headwaters and drier soil conditions and lower spring snowpacks reducing peak discharges despite increased precipitation during spring runoff in the Front Range and foothills headwaters under a warmer climate.

Simulating the hydrological impacts of land use conversion from annual crop to perennial forage in the Canadian Prairies using the Cold Regions Hydrological Modelling platform

Abstract. The Red River is one of the largest contributing sources of discharge and nutrients to the world's 10th largest freshwater lake, Lake Winnipeg. Conversion of large areas of annual cropland to perennial forage has been proposed as a strategy to reduce both flooding and nutrient export to Lake Winnipeg. Such reductions could occur either via a reduction in the concentration of nutrients in runoff or through changes in the basin-scale hydrology, resulting in a lower water yield and the concomitant export of nutrients. This study assessed the latter mechanism by using the physically based Cold Regions Hydrological Modelling platform to examine the hydrological impacts of land use conversion from annual crops to perennial forage in a subbasin of the La Salle River basin in Canada. This basin is a typical agricultural subbasin in the Red River Valley, characterised by flat topography, clay soils, and a cold subhumid, continental climate. Long-term simulations (1992–2013) of the major components of water balance were compared between canola and smooth bromegrass, representing a conversion from annual cropping systems to perennial forage. An uncertainty framework was used to represent a range of fall soil saturation status (0 % to 70 %), which governs the infiltration to frozen soil in the subsequent spring. The model simulations indicated that, on average, there was a 36.5 ± 6.6 % (36.5 ± 7.2 mm) reduction in annual cumulative discharge and a 29.9 ± 16.3 % (2.6 ± 1.6 m3 s−1) reduction in annual peak discharge due to forage conversion over the assessed period. These reductions were driven by reduced overland flow 52.9 ± 12.8 % (28.8 ± 10.1 mm), increased peak snowpack (8.1 ± 1.5 %, 7.8 ± 1.6 mm), and enhanced infiltration to frozen soils (66.7 ± 7.7 %, 141.5 ± 15.2 mm). Higher cumulative evapotranspiration (ET) from perennial forage (34.5 ± 0.9 %, 94.1 ± 2.5 mm) was also predicted by the simulations. Overall, daily soil moisture under perennial forage was 18.0 % (57.2 ± 1.2 mm) higher than that of crop simulation, likely due to the higher snow water equivalent (SWE) and enhanced infiltration. However, the impact of forage conversion on daily soil moisture varied interannually. Soil moisture under perennial forage stands could be either higher or lower than that of annual crops, depending on antecedent spring snowmelt infiltration volumes.

The cold regions hydrological modelling platform for hydrological diagnosis and prediction based on process understanding

Cold regions involve hydrological processes that are not often addressed appropriately in hydrological models. The Cold Regions Hydrological Modelling platform (CRHM) was initially developed in 1998 to assemble and explore the hydrological understanding developed from a series of research basins spanning Canada and international cold regions. Hydrological processes and basin response in cold regions are simulated in a flexible, modular, object-oriented, multiphysics platform. The CRHM platform allows for multiple representations of forcing data interpolation and extrapolation, hydrological model spatial and physical process structures, and parameter values. It is well suited for model falsification, algorithm intercomparison and benchmarking, and has been deployed for basin hydrology diagnosis, prediction, land use change and water quality analysis, climate impact analysis and flood forecasting around the world. This paper describes CRHM’s capabilities, and the insights derived by applying the model in concert with process hydrology research and using the combined information and understanding from research basins to predict hydrological variables, diagnose hydrological change and determine the appropriateness of model structure and parameterisations.

Recent hydrological response of glaciers in the Canadian Rockies to changing climate and glacier configuration

Abstract. Mountain snow and ice greatly influence the hydrological cycle of alpine regions by regulating both the quantity of and seasonal variations in water availability downstream. This study considers the combined impacts of climate and glacier changes due to recession on the hydrology and water balance of two high-elevation basins in the Canadian Rockies. A distributed, physically based, uncalibrated glacier hydrology model developed in the Cold Regions Hydrological Modelling platform (CRHM) was used to simulate the glacier mass balance and basin hydrology of the Peyto and Athabasca glacier basins in Alberta, Canada. Bias-corrected reanalysis data were used to drive the model. The model calculates the water balance of glacierized basins, influenced by the surface energy and mass balance, and considers the redistribution of snow by wind and avalanches. It was set up using hydrological response units based on elevation bands, surface slope, and aspect, as well as changing land cover. Aerial photos, satellite images and digital elevation models (DEMs) were assimilated to represent the changing configurations of glacier area and the exposure of ice and firn. Observations of glacier mass balance, snow, and glacier ice surface elevation changes at glacier and alpine tundra meteorological stations and streamflow discharge at the glacier outlets were used to evaluate the model performance. Basin hydrology was simulated over two periods, 1965–1975 and 2008–2018, using the observed glacier configurations for those time periods. Both basins have undergone continuous glacier loss over the last 3 to 5 decades, leading to a 6 %–31 % reduction in glacierized area, a 78 %–109 % increase in ice exposure, and changes to the elevation and slope of the glacier surfaces. Air temperatures are increasing, mainly due to increasing winter maximum and summer minimum daily temperatures. Annual precipitation has increased by less than 11 %, but rainfall ratios have increased by 29 %–44 %. The results show that changes in both climate and glacier configuration have influenced the melt rates and runoff and a shift of peak flows in the Peyto Glacier basin from August to July. Glacier melt contributions increased/decreased from 27 %–61 % to 43 %–59 % of the annual discharges. Recent discharges were 3 %–19 % higher than in the 1960s and 1970s. The results suggest that increased exposure of glacier ice and lower surface elevation due to glacier thinning were less influential than climate warming in increasing streamflow. Streamflow from these glaciers continues to increase.

Diagnosing changes in glacier hydrology from physical principles using a hydrological model with snow redistribution, sublimation, firnification and energy balance ablation algorithms

A comprehensive glacier hydrology model was developed within the Cold Regions Hydrological Modelling platform (CRHM) to include modules representing wind flow over complex terrain, blowing snow redistribution and sublimation by wind, snow redistribution by avalanches, solar irradiance to sloping surfaces, surface sublimation, glacier mass balance and runoff, meltwater and streamflow routing. The physically based glacier hydrology model created from these modules in CRHM was applied to simulate the hydrology of the instrumented, glacierized and rapidly deglaciating Peyto and Athabasca glacier research basins in the Canadian Rockies without calibration of parameters from streamflow. It was tested against observed albedo, point and aggregated glacier mass balance, and streamflow and found to successfully simulate surface albedo, snow redistribution, snow and glacier accumulation and ablation, mass balance and streamflow discharge, both when driven by in-situ observations and reanalysis forcing data. Long term modelling results indicate that the increases in discharge from the 1960s to the present are due to increased glacier ice melt contributions, despite declining precipitation and snow melt.

The sensitivity of snow hydrology to changes in air temperature and precipitation in three North American headwater basins

Whether or not the impact of warming on mountain snow and runoff can be offset by precipitation increases has not been well examined, but it is crucially important for future downstream water supply. Using the physically based Cold Regions Hydrological Modelling Platform (CRHM), elasticity (percent change in runoff divided by change in a climate forcing) and the sensitivity of snow regimes to perturbations were investigated in three well-instrumented mountain research basins spanning the northern North American Cordillera. Hourly meteorological observations were perturbed using air temperature and precipitation changes and were then used to force hydrological models for each basin. In all three basins, lower temperature sensitivities of annual runoff volume (≤ 6% °C−1) and higher sensitivities of peak snowpack (−17% °C−1) showed that annual runoff was far less sensitive to temperature than the snow regime. Higher and lower precipitation elasticities of annual runoff (1.5 – 2.1) and peak snowpack (0.7 – 1.1) indicated that the runoff change is primarily attributed to precipitation change and, secondarily, to warming. A low discrepancy between observed and simulated precipitation elasticities showed that the model results are reliable, and one can conduct sensitivity analysis. The air temperature elasticities, however, must be interpreted with care as the projected warmings range beyond the observed temperatures and, hence, it is not possible to test their reliability. Simulations using multiple elevations showed that the timing of peak snowpack was most sensitive to temperature. For the range of warming expected from North American climate model simulations, the impacts of warming on annual runoff, but not on peak snowpack, can be offset by the size of precipitation increases projected for the near-future period 2041–2070. To offset the impact of 2 °C warming on annual runoff, precipitation would need to increase by less than 5% in all three basins. To offset the impact of 2 °C warming on peak snowpack, however, precipitation would need to increase by 12% in Wolf Creek in Yukon Territory, 18% in Marmot Creek in the Canadian Rockies, and an amount greater than the maximum projected at Reynolds Mountain in Idaho. The role of increased precipitation as a compensator for the impact of warming on snowpack is more effective at the highest elevations and higher latitudes. Increased precipitation leads to resilient and strongly coupled snow and runoff regimes, contrasting sharply with the sensitive and weakly coupled regimes at low elevations and in temperate climate zones.

Advances in the simulation of nutrient dynamics in cold climate agricultural basins: Developing new nitrogen and phosphorus modules for the Cold Regions Hydrological Modelling Platform

Excess nutrients in aquatic ecosystems is a major water quality problem globally. Worsening eutrophication issues are notable in cold temperate areas, with pervasive problems in many agriculturally dominated catchments. Predicting nutrient export to rivers and lakes is particularly difficult in cold agricultural environments because of challenges in modelling snow, soil, frozen ground, climate, and anthropogenic controls. Previous research has shown that the use of many popular small basin nutrient models can be problematic in cold regions due to poor representation of cold region hydrology. In this study, the Cold Regions Hydrological Modelling Platform (CRHM), a modular modelling system, which has been widely deployed across Canada and cold regions worldwide, was used to address this problem. CRHM was extended to simulate biogeochemical and transport processes for nitrogen and phosphorus through a complex of new process-based modules that represent physicochemical processes in snow, soil and freshwater. Agricultural practices such as tillage and fertilizer application, which strongly impact the availability and release of soil nutrients, can be explicitly represented in the model. A test case in an agricultural basin draining towards Lake Winnipeg shows that the model can capture the extreme hydrology and nutrient load variability of small agricultural basins at hourly time steps. It was demonstrated that fine temporal resolutions are an essential modelling requisite to capture strong concentration changes in agricultural tributaries in cold agricultural environments. Within these ephemeral and intermittent streams, on average, 30%, 31%, 20%, and 16% of the total annual load of nitrate (NO3), ammonium (NH4), soluble reactive phosphorus (SRP), and particulate phosphorous (partP)NO3, NH4, SRP and partP occurred during the episodic snowmelt freshet (∼9 days, accounting for 21% of the annual flow), but shows extreme temporal variation. The new nutrient modules are critical tools for predicting nutrient export from small agricultural drainage basins in cold climates via better representation of key hydrological processes, and a temporal resolution more suited to capture dynamics of ephemeral and intermittent streams.

Diagnosis of future changes in hydrology for a Canadian Rockies headwater basin

Abstract. Climate change is anticipated to impact the hydrology of the Saskatchewan River, which originates in the Canadian Rockies mountain range. To better understand the climate change impacts in the mountain headwaters of this basin, a physically based hydrological model was developed for this basin using the Cold Regions Hydrological Modelling platform (CRHM) for Marmot Creek Research Basin (∼9.4 km2), located in the Front Ranges of the Canadian Rockies. Marmot Creek is composed of ecozones ranging from montane forests to alpine tundra and alpine exposed rock and includes both large and small clearcuts. The model included blowing and intercepted snow redistribution, sublimation, energy-balance snowmelt, slope and canopy effects on melt, Penman–Monteith evapotranspiration, infiltration to frozen and unfrozen soils, hillslope hydrology, streamflow routing, and groundwater components and was parameterised without calibration from streamflow. Near-surface outputs from the 4 km Weather Research and Forecasting (WRF) model were bias-corrected using the quantile delta mapping method with respect to meteorological data from five stations located from low-elevation montane forests to alpine ridgetops and running over October 2005–September 2013. The bias-corrected WRF outputs during a current period (2005–2013) and a future pseudo global warming period (PGW, 2091–2099) were used to drive model simulations to assess changes in Marmot Creek's hydrology. Under a “business-as-usual” forcing scenario, Representative Concentration Pathway 8.5 (RCP8.5) in PGW, the basin will warm up by 4.7 ∘C and receive 16 % more precipitation, which will lead to a 40 mm decline in seasonal peak snowpack, 84 mm decrease in snowmelt volume, 0.2 mm d−1 slower melt rate, and 49 d shorter snow-cover duration. The alpine snow season will be shortened by almost 1.5 months, but at some lower elevations there will be large decreases in peak snowpack (∼45 %) in addition to a shorter snow season. Declines in the peak snowpack will be much greater in clearcuts than under mature forest canopies. In alpine and treeline ecozones, blowing snow transport and sublimation will be suppressed by higher-threshold wind speeds for transport, in forest ecozones, sublimation losses from intercepted snow will decrease due to faster unloading and drip, and throughout the basin, evapotranspiration will increase due to a longer snow-free season and more rainfall. Runoff will begin earlier in all ecozones, but, as a result of variability in surface and subsurface hydrology, forested and alpine ecozones will generate the greatest runoff volumetric increases, ranging from 12 % to 25 %, whereas the treeline ecozone will have a small (2 %) decrease in runoff volume due to decreased melt volumes from smaller snowdrifts. The shift in timing in streamflow will be notable, with 236 % higher flows in spring months and 12 % lower flows in summer and 13 % higher flows in early fall. Overall, Marmot Creek's annual streamflow discharge will increase by 18 % with PGW, without a change in its streamflow generation efficiency, despite its basin shifting from primarily snowmelt runoff towards rainfall-dominated runoff generation.

Assimilating snow observations to snow interception process simulations

AbstractSnow interception is a crucial hydrological process in cold regions needleleaf forests, but is rarely measured directly. Indirect estimates of snow interception can be made by measuring the difference in the increase in snow accumulation between the forest floor and a nearby clearing over the course of a storm. Pairs of automatic weather stations with acoustic snow depth sensors provide an opportunity to estimate this, if snow density can be estimated reliably. Three approaches for estimating fresh snow density were investigated: weighted post‐storm density increments from the physically based Snobal model, fresh snow density estimated empirically from air temperature (Hedstrom, N. R., et al. [1998]. Hydrological Processes, 12, 1611–1625), and fresh snow density estimated empirically from air temperature and wind speed (Jordan, R. E., et al. [1999]. Journal of Geophysical Research, 104, 7785–7806). Automated snow depth observations from adjacent forest and clearing sites and estimated snow densities were used to determine snowstorm snow interception in a subalpine forest in the Canadian Rockies, Alberta, Canada. Then the estimated snow interception and measured interception information from a weighed, suspended tree and a time‐lapse camera were assimilated into a model, which was created using the Cold Regions Hydrological Modelling platform (CRHM), using Ensemble Kalman Filter or a simple rule‐based direct insertion method. Interception determined using density estimates from the Hedstrom‐Pomeroy fresh snow density equation agreed best with observations. Assimilating snow interception information from automatic snow depth measurements improved modelled snow interception timing by 7% and magnitude by 13%, compared to an open loop simulation driven by a numerical weather model; its accuracy was close to that simulated using locally observed meteorological data. Assimilation of tree‐measured snow interception improved the snow interception simulation timing and magnitude by 18 and 19%, respectively. Time‐lapse camera snow interception information assimilation improved the snow interception simulation timing by 32% and magnitude by 7%. The benefits of assimilation were greatly influenced by assimilation frequency and quality of the forcing data.

Shifting Hydrological Processes in a Canadian Agroforested Catchment due to a Warmer and Wetter Climate

This study examines the hydrological sensitivity of an agroforested catchment to changes in temperature and precipitation. A physically based hydrological model was created using the Cold Regions Hydrological Modelling platform to simulate the hydrological processes over 23 years in the Acadie River Catchment in southern Québec. The observed air temperature and precipitation were perturbed linearly based on existing climate change projections, with warming of up to 8 °C and an increase in total precipitation up to 20%. The results show that warming causes a decrease in blowing snow transport and sublimation losses from blowing snow, canopy-intercepted snowfall and the snowpack. Decreasing blowing snow transport leads to reduced spatial variability in peak snow water equivalent (SWE) and a more synchronized snow cover depletion across the catchment. A 20% increase in precipitation is not sufficient to counteract the decline in annual peak SWE caused by a 1 °C warming. On the other hand, peak spring streamflow increases by 7% and occurs 20 days earlier with a 1 °C warming and a 20% increase in precipitation. However, when warming exceeds 1.5 °C, the catchment becomes more rainfall dominated and the peak flow and its timing follows the rainfall rather than snowmelt regime. Results from this study can be used for sustainable farming development and planning in regions with hydroclimatic characteristics similar to the Acadie River Catchment, where climate change may have a significant impact on the dominating hydrological processes.

Modelling the effects of permafrost loss on discharge from a wetland‐dominated, discontinuous permafrost basin

AbstractPermafrost degradation in the peat‐rich southern fringe of the discontinuous permafrost zone is catalysing substantial changes to land cover with expansion of permafrost‐free wetlands (bogs and fens) and shrinkage of forest‐dominated permafrost peat plateaux. Predicting discharge from headwater basins in this region depends upon understanding and numerically representing the interactions between storage and discharge within and between the major land cover types and how these interactions are changing. To better understand the implications of advanced permafrost thaw‐induced land cover change on wetland discharge, with all landscape features capable of contributing to drainage networks, the hydrological behaviour of a channel fen sub‐basin in the headwaters of Scotty Creek, Northwest Territories, Canada, dominated by peat plateau–bog complexes, was modelled using the Cold Regions Hydrological Modelling platform for the period of 2009 to 2015. The model construction was based on field water balance observations, and performance was deemed adequate when evaluated against measured water balance components. A sensitivity analysis was conducted to assess the impact of progressive permafrost loss on discharge from the sub‐basin, in which all units of the sub‐basin have the potential to contribute to the drainage network, by incrementally reducing the ratio of wetland to plateau in the modelled sub‐basin. Simulated reductions in permafrost extent decreased total annual discharge from the channel fen by 2.5% for every 10% decrease in permafrost area due to increased surface storage capacity, reduced run‐off efficiency, and increased landscape evapotranspiration. Runoff ratios for the fen hydrological response unit dropped from 0.54 to 0.48 after the simulated 50% permafrost area loss with a substantial reduction of 0.47 to 0.31 during the snowmelt season. The reduction in peat plateau area resulted in decreased seasonal variability in discharge due to changes in the flow path routing, with amplified low flows associated with small increases in subsurface discharge, and decreased peak discharge with large reductions in surface run‐off.

Detecting intercepted snow on mountain needleleaf forest canopies using satellite remote sensing

Snow interception in cold regions needleleaf forest canopies is a crucial process that controls local snow accumulation and redistribution over >20% of the Earth's land surface. Various ground-based methods exist to measure intercepted snow load, however all are based on single-tree measurements and are difficult to implement. No research has focussed on detecting large areal intercepted snow loads and no studies have assessed the use of satellite observations. In this study, four remote sensing indices (NDSI, NDVI, albedo, and land surface temperature (LST)) were retrieved from Landsat images to study their sensitivity to canopy intercepted snow and the possibility of using them to detect the presence of intercepted snow. The results indicate that presence of intercepted snow on canopy increased NDSI and albedo, but decreased NDVI. Intercepted snow presence also decreased the areal variability of NDSI and NDVI while increasing that of albedo. For these three indices, the differences between snow-free and snowcovered canopies were correlated to topography and forest canopy cover. Of these indices, NDSI changed the greatest. Intercepted snow noticeably decreased the LST difference between forest and open areas in springtime while the influence in wintertime was relatively smaller. An intercepted snow detection approach that uses both NDSI and NDVI to classify pixels into either snowcovered canopy or other (snow-free canopy and non-forest areas) is proposed here. A case study applying this approach compared remote sensing detection to simulations by the snow interception and sublimation model implemented in the Cold Regions Hydrological Modelling platform (CRHM). This used local meteorological observations from the pine, spruce and fir forest covered Marmot Creek Research Basin in the Canadian Rockies. The remote sensing detection of intercepted snow agreed well with CRHM simulations for continuous forests (83%) and less well for sparse forests (72%) and clearings with small trees (70%). Therefore, the approach is suitable for intercepted snow detection over continuous evergreen canopies. This technique provides a new capability for large-scale snow interception model validation and data assimilation to cold regions hydrological forecasting models.

Simulating cold-region hydrology in an intensively drained agricultural watershed in Manitoba, Canada, using the Cold Regions Hydrological Model

Abstract. Etrophication and flooding are perennial problems in agricultural watersheds of the northern Great Plains. A high proportion of annual runoff and nutrient transport occurs with snowmelt in this region. Extensive surface drainage modification, frozen soils, and frequent backwater or ice-damming impacts on flow measurement represent unique challenges to accurately modelling watershed-scale hydrological processes. A physically based, non-calibrated model created using the Cold Regions Hydrological Modelling platform (CRHM) was parameterized to simulate hydrological processes within a low slope, clay soil, and intensively surface drained agricultural watershed. These characteristics are common to most tributaries of the Red River of the north. Analysis of the observed water level records for the study watershed (La Salle River) indicates that ice cover and backwater issues at time of peak flow may impact the accuracy of both modelled and measured streamflows, highlighting the value of evaluating a non-calibrated model in this environment. Simulations best matched the streamflow record in years when peak and annual discharges were equal to or above the medians of 6.7 m3 s−1 and 1.25 × 107 m3, respectively, with an average Nash–Sutcliffe efficiency (NSE) of 0.76. Simulation of low-flow years (below the medians) was more challenging (average NSE < 0), with simulated discharge overestimated by 90 % on average. This result indicates the need for improved understanding of hydrological response in the watershed under drier conditions. Simulation during dry years was improved when infiltration was allowed prior to soil thaw, indicating the potential importance of preferential flow. Representation of in-channel dynamics and travel time under the flooded or ice-jam conditions should also receive attention in further model development efforts. Despite the complexities of the study watershed, simulations of flow for average to high-flow years and other components of the water balance were robust (snow water equivalency (SWE) and soil moisture). A sensitivity analysis of the flow routing model suggests a need for improved understanding of watershed functions under both dry and flooded conditions due to dynamic routing conditions, but overall CRHM is appropriate for simulation of hydrological processes in agricultural watersheds of the Red River. Falsifications of snow sublimation, snow transport, and infiltration to frozen soil processes in the validated base model indicate that these processes were very influential in stream discharge generation.

Diagnosis of the hydrology of a small Arctic basin at the tundra-taiga transition using a physically based hydrological model

A better understanding of cold regions hydrological processes and regimes in transitional environments is critical for predicting future Arctic freshwater fluxes under climate and vegetation change. A physically based hydrological model using the Cold Regions Hydrological Model platform was created for a small Arctic basin in the tundra-taiga transition region. The model represents snow redistribution and sublimation by wind and vegetation, snowmelt energy budget, evapotranspiration, subsurface flow through organic terrain, infiltration to frozen soils, freezing and thawing of soils, permafrost and streamflow routing. The model was used to reconstruct the basin water cycle over 28years to understand and quantify the mass fluxes controlling its hydrological regime. Model structure and parameters were set from the current understanding of Arctic hydrology, remote sensing, field research in the basin and region, and calibration against streamflow observations. Calibration was restricted to subsurface hydraulic and storage parameters. Multi-objective evaluation of the model using observed streamflow, snow accumulation and ground freeze/thaw state showed adequate simulation. Significant spatial variability in the winter mass fluxes was found between tundra, shrubs and forested sites, particularly due to the substantial blowing snow redistribution and sublimation from the wind-swept upper basin, as well as sublimation of canopy intercepted snow from the forest (about 17% of snowfall). At the basin scale, the model showed that evapotranspiration is the largest loss of water (47%), followed by streamflow (39%) and sublimation (14%). The models streamflow performance sensitivity to a set of parameter was analysed, as well as the mean annual mass balance uncertainty associated with these parameters.

Impact of climate warming on snow processes in Ny-Ålesund, a polar maritime site at Svalbard

The impact of observed changes in air temperature and precipitation from 1969 to 2013 and climate projections for 2050 and 2080 at Ny-Ålesund, an arctic research station on Spitzbergen Island in the Svalbard Archipelago on snow hydrological processes, were analyzed using snow accumulation and ablation algorithms in the physically based Cold Regions Hydrological Modelling platform (CRHM). The climate projections were obtained from phase 5 of the Coupled Model Intercomparison Project (CMIP5), with a focus on the snow-dominated period (October to June). To identify the potential effects of increasing temperature and precipitation, a model sensitivity analysis (1°C to 5°C), with and without a 25% increase in precipitation, was run on CRHM snow processes. The results indicated that the greatest observed warming was during the early snow season (October–February), with increases of 0.8 and 0.9°C decade−1 for maximum (Tmax) and minimum (Tmin) temperatures, respectively. There was also a significant increase in annual and winter precipitation (24mm decade−1). The late snow season (March–June) also had a marked increase in temperature (0.5 and 0.69°C decade−1 for Tmax and Tmin respectively), but no significant change in precipitation. These changes lead to a significant increase in the number of days with rainfall rather than snowfall. The sensitivity analysis indicated that mean snow water equivalent snowpack will decrease by 10.2% (early snow season) and 11.1% (late snow season) per degree of increased air temperature. For each degree of temperature increase, the modelled peak snow-water-equivalent (SWE) declined by 6.9%, duration of snowpack declined 11days, and the number of days with rain increased 43% for the early snow season and 12.8% for the late snow season. A warmer climate also leads to markedly decreased surface snow sublimation and the fraction of snowfall eroded and transported by blowing snow. For most snowpack parameters analyzed, the response to warming accelerates with increased warming, especially above 3°C. A 25% increase in precipitation partially counteracted the response to warming, with the greatest effect on peak SWE.

Description of current and future snow processes in a small basin in the Bavarian Alps

Snow cover dynamics in alpine regions play a crucial role in view of the water balance of head water catchments. The temporal storage of water in form of snow and ice leads to a decoupling of precipitation and runoff. Changes in the volume and the temporal dynamics of the snow storage lead to modified runoff regimes and can influence the frequency of low flow events and floods. For a better estimation of the possible range and direction of future changes, projection runs can be realized by using process-based models. In this study, the Cold Regions Hydrological Modelling platform (CRHM) is used to compile such a model for simulating the snow cover development within research catchment Zugspitze (RCZ; 11.4 km2/Germany). Therefore, the catchment is divided into four hydrological response units (HRUs), able to cover the physiographic characteristics in four elevation zones. The model is evaluated over snow depth measurements. The range of variability within and differences between the HRUs are analyzed, and future projections (2001–2100) are performed on the basis of three different WETTREG realizations. It could be shown that CRHM is able to reproduce the snow cover dynamics very well and that the ongoing climate change does have an identifiable influence on the average extent and size of the snow storage. Furthermore, it could be shown that variations in snow cover dynamics within the RCZ are strongly connected to NAO.

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

AbstractThe 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.

Impact of antecedent conditions on simulations of a flood in a mountain headwater basin

A devastating flood struck Southern Alberta in late June 2013, with much of its streamflow generation in the Front Ranges of the Rocky Mountains, west of Calgary. To better understand streamflow generation processes and their sensitivity to initial conditions, a physically based hydrological model was developed using the Cold Regions Hydrological Modelling platform (CRHM) to simulate the flood for the Marmot Creek Research Basin (~9.4 km2). The modular model includes major cold and warm season hydrological processes including snow redistribution, sublimation, melt, runoff over frozen and unfrozen soils, evapotranspiration, subsurface runoff on hillslopes, groundwater recharge and discharge and streamflow routing. Uncalibrated simulations were conducted for eight hydrological years and generally matched streamflow observations well, with a NRMSD of 52%, small model bias (−3%) and a Nash–Sutcliffe efficiency (NSE) of 0.71. The model was then used to diagnose the responses of hydrological processes in 2013 flood from different ecozones in Marmot Creek: alpine, treeline, montane forest and large and small forest clearings to better understand spatial variations in the flood runoff generation mechanisms. To examine the sensitivity to antecedent conditions, ‘virtual’ flood simulations were conducted using a week (17 to 24 June 2013) of flood meteorology imposed on the meteorology of the same period in other years (2005 to 2012), or switched with the meteorology of one week in different months (May to July) of 2013. Sensitivity to changing precipitation and land cover was assessed by varying the precipitation amount during the flood and forest cover and soil storage capacity in forest ecozone. The results show that runoff efficiency increases rapidly with antecedent snowpack and soil moisture storage with the highest runoff response to rainfall from locations in the basin where there are recently melted or actively melting snowpacks and resulting high soil moisture or frozen soils. The impact of forest canopy on flooding is negligible, but flood peak doubles if forest canopy removal is accompanied by 50% reduction in water storage capacity in the basin. Copyright © 2016 John Wiley & Sons, Ltd.

Simulating cold regions hydrological processes using a modular model in the west of China

The Cold Regions Hydrological Model platform (CRHM), a flexible object-oriented modeling system, was devised to simulate cold regions hydrological processes and predict streamflow by its capability to compile cold regions process modules into purpose-built models. In this study, the cold regions hydrological processes of two basins in western China were evaluated using CRHM models: Binggou basin, a high alpine basin where runoff is mainly caused by snowmelt, and Zuomaokong basin, a steppe basin where the runoff is strongly affected by soil freezing/thawing. The flexibility and modular structure of CRHM permitted model structural intercomparison and process falsification within the same model framework to evaluate the importance of snow energy balance, blowing snow and frozen soil infiltration processes to successful modeling in the cold regions of western China. Snow accumulation and ablation processes were evaluated at Binggou basin by testing and comparing similar models that contained different levels of complexity of snow redistribution and ablation modules. The comparison of simulated snow water equivalent with observations shows that the snow accumulation/ablation processes were simulated much better using an uncalibrated, physically based energy balance snowmelt model rather than with a calibrated temperature index snowmelt model. Simulated seasonal snow sublimation loss was 138 mm water equivalent in the alpine region of Binggou basin, which accounts for 47 % of 291 mm water equivalent of snowfall, and half of this sublimation loss is attributed to 70 mm water equivalent of sublimation from blowing snow particles. Further comparison of simulated results through falsification of different snow processes reveals that estimating blowing snow transport processes and sublimation loss is vital for accurate snowmelt runoff calculations in this region. The model structure with the energy balance snowmelt and blowing snow components performed well in reproducing the measured streamflow using minimal calibration, with R2 of 0.83 and NSE of 0.76. The influence of frozen soil and its thaw on runoff generation was investigated at Zuomaokong basin by comparing streamflow simulated by similar CRHM models with and without an infiltration to frozen soil algorithm. The comparison of simulated streamflow with observation shows that the model which included an algorithm describing frozen soil infiltration simulated the main runoff events for the spring thawing period better than that which used an unfrozen infiltration routine, with R2 of 0.87 and NSE of 0.79. Overall, the test results for the two basins show that hydrological models that use appropriate cold regions algorithms and a flexible spatial structure can predict cold regions hydrological processes and streamflow with minimal calibration and can even perform better than more heavily calibrated models in this region. Given that CRHM and most of its algorithms were developed in western Canada, this is encouraging for predicting hydrology in ungauged cold region basins around the world.

The effect of slope aspect on the response of snowpack to climate warming in the Pyrenees

The aim of this study was to analyse the effect of slope aspect on the response of snowpack to climate warming in the Pyrenees. For this purpose, data available from five automatic weather stations were used to simulate the energy and mass balance of snowpack, assuming different magnitudes of an idealized climate warming (upward shifting of 1, 2 and 3 °C the temperature series). Snow energy and mass balance were simulated using the Cold Regions Hydrological Modelling platform (CRHM). CRHM was used to create a model that enabled correction of the all-wave incoming radiation fluxes from the observation sites for various slope aspects (N, NE, E, SE, S, SW,W,NW and flat areas), which enabled assessment of the differential impact of climate warming on snow processes on mountain slopes. The results showed that slope aspect was responsible for substantial variability in snow accumulation and the duration of the snowpack. Simulated variability markedly increased with warmer temperature conditions. Annual maximum snow accumulation (MSA) and annual snowpack duration (ASD) showed marked sensitivity to a warming of 1 °C. Thus, the sensitivity of the MSA in flat areas ranged from 11 to 17 % per degree C amongst the weather stations, and the ASD ranged from 11 to 20 days per degree C. There was a clear increase in the sensitivity of the snowpack to climate warming on those slopes that received intense solar radiation (S, SE and SW slopes) compared with those slopes where the incident radiation was more limited (N, NE and NW slopes). The sensitivity of the MSA and the ASD increased as the temperature increased, particularly on the most irradiated slopes. Large interannual variability was also observed. Thus, with more snow accumulation and longer duration the sensitivity of the snowpack to temperature decreased, especially on south-facing slopes.