Delineation of flood generating processes and their hydrological response

Knowledge about flood generating processes can be beneficial for numerous applications. Especially in the context of climate change impact assessment, daily patterns of meteorological and catchment state conditions leading to flood events (i.e., storylines) may be of value. Here, we propose an approach to identify storylines of flood generation using daily weather and snow cover observations. The approach is tested for and applied to a typical pre‐Alpine catchment in the period between 1961 and 2014. Five precipitation parameters were determined that describe temporal and spatial characteristics of the flood associated precipitation events. The catchment's snow coverage was derived using statistical relationships between a satellite‐derived snow cover climatology and station snow measurements. Moreover, (pre‐) event snow melt sums were estimated using a temperature‐index model. Based on the precipitation and catchment state parameters, 5 storylines were identified with a cluster analysis: These are (a) long duration, low intensity precipitation events with high precipitation depths, (b) long duration precipitation events with high precipitation depths and episodes of high intensities, (c) shorter duration events with high or (d) low precipitation intensity, respectively, and (e) rain‐on‐snow events. The event groups have distinct hydrological characteristics that can largely be explained by the storylines' respective properties. The long duration, high intensity storyline leads to the most adverse hydrological response, namely, a combination of high peak magnitudes, high volumes, and long durations of threshold exceedance. The results show that flood generating processes in mesoscale catchments can be distinguished on the basis of daily meteorological and catchment state parameters and that these process types can explain the hydrological flood properties in a qualitative way. Hydrological simulations of daily resolution can thus be analysed with respect to flood generating processes.     

Surface water input from snowmelt and rain throughfall in western juniper: potential impacts of climate change and shifts in semi‐arid vegetation

In the western USA, shifts from snow to rain precipitation regimes and increases in western juniper cover in shrub‐dominated landscapes can alter surface water input via changes in snowmelt and throughfall. To better understand how shifts in both precipitation and semi‐arid vegetation cover alter above‐ground hydrological processes, we assessed how rain interception differs between snow and rain surface water input; how western juniper alters snowpack dynamics; and how these above‐ground processes differ across western juniper, mountain big sagebrush and low sagebrush plant communities. We collected continuous surface water input with four large lysimeters, interspace and below‐canopy snow depth data and conducted periodic snow surveys for two consecutive water years (2013 and 2014). The ratio of interspace to below‐canopy surface water input was greater for snow relative to rain events, averaging 79.4% and 54.8%, respectively. The greater surface water input ratio for snow is in part due to increased deposition of redistributed snow under the canopy. We simulated above‐ground energy and water fluxes in western juniper, low sagebrush and mountain big sagebrush for two 8‐year periods under current and projected mid‐21st century warmer temperatures with the Simultaneous Heat and Water (SHAW) model. Juniper compared with low and mountain sagebrush reduced surface water input by an average of 138 mm or 24% of the total site water budget. Conversely, warming temperatures reduced surface water input by only an average of 14 mm across the three vegetation types. The future (warmer) simulations resulted in earlier snow disappearance and surface water input by 51 and 45 days, respectively, across juniper, low sagebrush and mountain sagebrush. Information from this study can help land managers in the sagebrush steppe understand how both shifts in climate and semi‐arid vegetation will alter fundamental hydrological processes. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

Gridded climate data products are an alternative to instrumental measurements as inputs to rainfall–runoff models

Rainfall–runoff models are widely used to predict flows using observed (instrumental) time series of air temperature and precipitation as inputs. Poor model performance is often associated with difficulties in estimating catchment‐scale meteorological variables from point observations. Readily available gridded climate products are an underutilized source of temperature and precipitation time series for rainfall–runoff modelling, which may overcome some of the performance issues associated with poor‐quality instrumental data in small headwater monitoring catchments. Here we compare the performance of instrumental measured and E‐OBS gridded temperature and precipitation time series as inputs in the rainfall–runoff models “PERSiST” and “HBV” for flow prediction in six small Swedish catchments. For both models and most catchments, the gridded data produced statistically better simulations than did those obtained using instrumental measurements. Despite the high correspondence between instrumental and gridded temperature, both temperature and precipitation were responsible for the difference. We conclude that (a) gridded climate products such as the E‐OBS dataset could be more widely used as alternative input to rainfall–runoff models, even when instrumental measurements are available, and (b) the processing applied to gridded climate products appears to provide a more realistic approximation of small catchment‐scale temperature and precipitation patterns needed for flow simulations. Further research on this issue is needed and encouraged.     

Impacts of climate change in three hydrologic regimes in British Columbia,Canada

Hydrologic modelling has been applied to assess the impacts of projected climate change within three study areas in the Peace, Campbell and Columbia River watersheds of British Columbia, Canada. These study areas include interior nival (two sites) and coastal hybrid nival–pluvial (one site) hydro‐climatic regimes. Projections were based on a suite of eight global climate models driven by three emission scenarios to project potential climate responses for the 2050s period (2041–2070). Climate projections were statistically downscaled and used to drive a macro‐scale hydrology model at high spatial resolution. This methodology covers a large range of potential future climates for British Columbia and explicitly addresses both emissions and global climate model uncertainty in the final hydrologic projections. Snow water equivalent is projected to decline throughout the Peace and Campbell and at low elevations within the Columbia. At high elevations within the Columbia, snow water equivalent is projected to increase with increased winter precipitation. Streamflow projections indicate timing shifts in all three watersheds, predominantly because of changes in the dynamics of snow accumulation and melt. The coastal hybrid site shows the largest sensitivity, shifting to more rainfall‐dominated system by mid‐century. The two interior sites are projected to retain the characteristics of a nival regime by mid‐century, although streamflow‐timing shifts result from increased mid‐winter rainfall and snowmelt, and earlier freshet onset. Copyright © 2012 John Wiley & Sons, Ltd.     

Insights on stream temperature processes through development of a coupled hydrologic and stream temperature model for forested coastal headwater catchments

Stream temperature controls a number of biological, chemical, and physical processes occurring in aquatic environments. Transient snow cover and advection associated with lateral throughflow inputs can have a dominant influence on stream thermal regimes for headwater catchments in the rain‐on‐snow zone. Most existing stream temperature models lack the ability to properly simulate these processes. We developed and evaluated a conceptual‐parametric catchment‐scale stream temperature model that includes the role of transient snow cover and lateral advection associated with throughflow. The model consists of routines for simulating canopy interception, snow accumulation and melt, hillslope throughflow runoff and temperature, and stream channel energy exchange processes. The model was used to predict discharge and stream temperature for a small forested headwater catchment near Vancouver, Canada, using long‐term (1963–2013) weather data to compute model forcing variables. The model was evaluated against 4 years of observed stream temperature. The model generally predicted daily mean stream temperature accurately (annual RMSE between 0.57 and 1.24 °C) although it overpredicted daily summer stream temperatures by up to 3 °C during extended low streamflow conditions. Model development and testing provided insights on the roles of advection associated with lateral throughflow, channel interception of snow, and surface–subsurface water interactions on stream thermal regimes. This study shows that a relatively simple but process‐based model can provide reasonable stream temperature predictions for forested headwater catchments located in the rain‐on‐snow zone.     

A simple method to combine snow height and meteorological observations to estimate winter precipitation at sub-daily resolution

ABSTRACT

Measuring winter solid and liquid precipitation with high temporal resolution in remote or higher elevation regions is a challenging task because of undercatch and power supply issues. However, the number of micro-meteorological stations and ultrasonic height sensors in mountain regions is steadily increasing. To gain more benefit from such stations, a new simple approach for EStimating SOlid and LIquid Precipitation (ESOLIP) is presented. The method consists of three main steps: (1) definition of precipitation events using micro-meteorological data, (2) quantification of solid and liquid precipitation using wet-bulb temperature and filtered snow height and (3) calculation of fresh snow density. ESOLIP performance was validated using data from a heated rain gauge, snow pillow and daily manual observations both for single precipitation events and over three winter seasons. Results proved ESOLIP as an effective approach for precipitation quantification, where snow height observations and basic meteorological measurements (air temperature, solar radiation, wind speed, relative humidity), but no reliable rain gauges are available.     

Evaluation of snow water equivalent datasets over the Saint‐Maurice river basin region of southern Québec

A 10‐km gridded snow water equivalent (SWE) dataset is developed over the Saint‐Maurice River basin region in southern Québec from kriging of observed snow survey data for evaluation of SWE products. The gridded SWE dataset covers 1980–2014 and is based on manual gravimetric snow surveys carried out on February 1, March 1, March 15, April 1, and April 15 of each snow season, which captures the annual maximum SWE (SWEM) with a mean interpolation error of ±19%. The dataset is used to evaluate SWEM from a range of sources including satellite retrievals, reanalyses, Canadian regional climate models, and the Canadian Meteorological Centre operational snow depth analysis. We also evaluate a number of solid precipitation datasets to determine their contribution to systematic errors in estimated SWEM. None of the evaluated datasets is able to provide estimates of SWEM that are within operational requirements of ±15% error, and insufficient solid precipitation is determined to be one of the main reasons. The Climate System Forecast Reanalysis is the only dataset where snowfall is sufficiently large to generate SWEM values comparable to observations. Inconsistencies in precipitation are also found to have a strong impact on year‐to‐year variability in SWEM dataset performance and spread. Version 3.6.1 of the Canadian Land Surface Scheme land surface scheme driven with ERA‐Interim output downscaled by Version 5.0.1 of the Canadian Regional Climate Model was the best physically based model at explaining the observed spatial and temporal variability in SWEM (root‐mean‐square error [RMSE] = 33%) and has potential for lower error with adjusted precipitation. Operational snow products relying on the real‐time snow depth observing network performed poorly due to a lack of real‐time data and the strong local scale variability of point snow depth observations. The results underscore the need for more effort to be invested in improving solid precipitation estimates for use in snow hydrology applications.     

An approach to using snow areal depletion curves inferred from MODIS and its application to land surface modelling in Alaska

Snowcover areal depletion curves inferred from the moderate resolution imaging spectroradiometer (MODIS) are validated and then applied in NASA's catchment‐based land surface model (CLSM) for numerical simulations of hydrometeorological processes in the Kuparuk River basin (KRB) of Alaska. The results demonstrate that the MODIS snowcover fraction f derived from a simple relationship in terms of the normalized difference snow index compares well with Landsat values over the range 20 ≤ f ≤ 100%. For f < 20%, however, MODIS 500 m subpixel data underestimate the amount of snow by up to 13% compared with Landsat at spatial resolutions of 30 m binned to equivalent 500 m pixels. After a bias correction, MODIS snow areal depletion curves during the spring transition period of 2002 for the KRB exhibit similar features to those derived from surface‐based observations. These results are applied in the CLSM subgrid‐scale snow parameterization that includes a deep and a shallow snowcover fraction. Simulations of the evolution of the snowpack and of freshwater discharge rates for the KRB over a period of 11 years are then analysed with the inclusion of this feature. It is shown that persistent snowdrifts on the arctic landscape, associated with a secondary plateau in the snow areal depletion curves, are hydrologically important. An automated method is developed to generate the shallow and deep snowcover fractions from MODIS snow areal depletion curves. This provides the means to apply the CLSM subgrid‐scale snow parameterization in all watersheds subject to seasonal snowcovers. Improved simulations and predictions of the global surface energy and water budgets are expected with the incorporation of the MODIS snow data into the CLSM. Copyright © 2005 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.     

A temperature‐precipitation‐based model of thirty‐year mean snowpack accumulation and melt in Oregon,USA

High‐resolution, spatially extensive climate grids can be useful in regional hydrologic applications. However, in regions where precipitation is dominated by snow, snowmelt models are often used to account for timing and magnitude of water delivery. We developed an empirical, nonlinear model to estimate 30‐year means of monthly snowpack and snowmelt throughout Oregon. Precipitation and temperature for the period 1971–2000, derived from 400‐m resolution PRISM data, and potential evapotranspiration (estimated from temperature and day length) drive the model. The model was calibrated using mean monthly data from 45 SNOTEL sites and accurately estimated snowpack at 25 validation sites: R² = 0·76, Nash‐Sutcliffe Efficiency (NSE) = 0·80. Calibrating it with data from all 70 SNOTEL sites gave somewhat better results (R² = 0·84, NSE = 0·85). We separately applied the model to SNOTEL stations located < 200 and ≥ 200 km from the Oregon coast, since they have different climatic conditions. The model performed equally well for both areas. We used the model to modify moisture surplus (precipitation minus potential evapotranspiration) to account for snowpack accumulation and snowmelt. The resulting values accurately reflect the shape and magnitude of runoff at a snow‐dominated basin, with low winter values and a June peak. Our findings suggest that the model is robust with respect to different climatic conditions, and that it can be used to estimate potential runoff in snow‐dominated basins. The model may allow high‐resolution, regional hydrologic comparisons to be made across basins that are differentially affected by snowpack, and may prove useful for investigating regional hydrologic response to climate change. Published in 2011 by John Wiley & Sons, Ltd.     

Köppen Climate Types in Observed and Simulated Climates

The Köppen climate classification was applied to the observed gridded climatological sets and the outputs of four general circulation models (GCMs) over the continents of the Earth. All data had been acquired via the Data Distribution Centre established by the Intergovernmental Panel on Climate Change. The ability of the GCMs to simulate the Köppen climate zones identified in the real data was explored and possible future (global warming) changes in the climate types' distribution for each GCM were assessed. Differences in the area distributions derived from the GCMs' recent climate simulations give evidence about uncertainties generally involved in climate models. As to the global warming simulations, all GCM projections of warming climate (horizon 2050) show that the zones representing tropical rain climates and dry climates become larger, and the zones identified with boreal forest and snow climates together with the polar climates are smaller.     

Investigation of the snow‐cover dynamics in the Upper Euphrates Basin of Turkey using remotely sensed snow‐cover products and hydrometeorological data

The Euphrates and Tigris rivers serve as the most important water resources in the Middle East. Precipitation in this region falls mostly in the form of snow over the higher elevations of the Euphrates Basin and remains on the ground for nearly half of the year. This snow‐covered area (SCA) is a key element of the hydrological cycle, and monitoring the SCA is crucial for making accurate forecasts of snowmelt discharge, especially for energy production, flood control, irrigation, and reservoir‐operation optimization in the Upper Euphrates (Karasu) Basin. Remote sensing allows the detection of the spatio‐temporal patterns of snow cover across large areas in inaccessible terrain, such as the eastern part of Turkey, which is highly mountainous. In this study, a seasonal evaluation of the snow cover from 2000 to 2009 was performed using 8‐day snow‐cover products (MOD10C2) and the daily snow‐water equivalent (SWE) product. The values of SWE products were obtained using an assimilation process based on the Helsinki University of Technology model using equal area Special Sensor Microwave Imager (SSM/I) Earth‐gridded advanced microwave scanning radiometer—EOS daily brightness‐temperature values. In the Karasu Basin, the SCA percentage for the winter period is 80–90%. The relationship between the SCA and the runoff during the spring period is analysed for the period from 2004 to 2009. An inverse linear relationship between the normalized SCA and the normalized runoff values was obtained (r = 0·74). On the basis of the monthly mean temperature, total precipitation and snow depth observed at meteorological stations in the basin, the decrease in the peak discharges, and early occurrences of the peak discharges in 2008 and 2009 are due to the increase in the mean temperature and the decrease in the precipitation in April. Copyright © 2011 John Wiley & Sons, Ltd.     

Projected 21st century climate change on snow conditions over Shasta Dam watershed by means of dynamical downscaling

Snow is an important component of the Earth's climate system and is particularly vulnerable to global warming. It has been suggested that warmer temperatures may cause significant declines in snow water content and snow cover duration. In this study, snowfall and snowmelt were projected by means of a regional climate model that was coupled to a physically based snow model over Shasta Dam watershed to assess changes in snow water content and snow cover duration during the 21st century. This physically based snow model requires both physical data and future climate projections. These physical data include topography, soils, vegetation, and land use/land cover, which were collected from associated organizations. The future climate projections were dynamically downscaled by means of the regional climate model under 4 emission scenarios simulated by 2 general circulation models (fifth‐generation of the ECHAM general circulation model and the third‐generation atmospheric general circulation model). The downscaled future projections were bias corrected before projecting snowfall and snowmelt processes over Shasta Dam watershed during 2010–2099. This study's results agree with those of previous studies that projected snow water equivalent is decreasing by 50–80% whereas the fraction of precipitation falling as snowfall is decreasing by 15% to 20%. The obtained projection results show that future snow water content will change in both time and space. Furthermore, the results confirm that physical data such as topography, land cover, and atmospheric–hydrologic data are instrumental in the studies on the impact of climate change on the water resources of a region.     

Effects of mid‐twenty‐first century climate and land cover change on the hydrology of the Puget Sound basin,Washington

The distributed hydrology–soil–vegetation model (DHSVM) was used to study the potential impacts of projected future land cover and climate change on the hydrology of the Puget Sound basin, Washington, in the mid‐twenty‐first century. A 60‐year climate model output, archived for the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4), was statistically downscaled and used as input to DHSVM. From the DHSVM output, we extracted multi‐decadal averages of seasonal streamflow, annual maximum flow, snow water equivalent (SWE), and evapotranspiration centred around 2030 and 2050. Future land cover was represented by a 2027 projection, which was extended to 2050, and DHSVM was run (with current climate) for these future land cover projections. In general, the climate change signal alone on sub‐basin streamflow was evidenced primarily through changes in the timing of winter and spring runoff, and slight increases in the annual runoff. Runoff changes in the uplands were attributable both to climate (increased winter precipitation, less snow) and land cover change (mostly reduced vegetation maturity). The most climatically sensitive parts of the uplands were in areas where the current winter precipitation is in the rain–snow transition zone. Changes in land cover were generally more important than climate change in the lowlands, where a substantial change to more urbanized land use and increased runoff was predicted. Both the annual total and seasonal distribution of freshwater flux to Puget Sound are more sensitive to climate change impacts than to land cover change, primarily because most of the runoff originates in the uplands. Both climate and land cover change slightly increase the annual freshwater flux to Puget Sound. Changes in the seasonal distribution of freshwater flux are mostly related to climate change, and consist of double‐digit increases in winter flows and decreases in summer and fall flows. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluation of the MODIS snow cover fraction product

Eleven years of daily 500 m gridded Terra Moderate Resolution Imaging Spectroradiometer (MODIS) (MOD10A1) snow cover fraction (SCF) data are evaluated in terms of snow presence detection in Colorado and Washington states. The SCF detection validation study is performed using in‐situ measurements and expressed in terms of snow and land detection and misclassification frequencies. A major aspect addressed in this study involves the shifting of pixel values in time due to sensor viewing angles and gridding artifacts of MODIS sensor products. To account for this error, 500 m gridded pixels are grouped and aggregated to different‐sized areas to incorporate neighboring pixel information. With pixel aggregation, both the probability of detection (POD) and the false alarm ratios increase for almost all cases. Of the false negative (FN) and false positive values (referred to as the total error when combined), FN estimates dominate most of the total error and are greatly reduced with aggregation. The greatest POD increases and total error reductions occur with going from a single 500 m pixel to 3×3‐pixel averaged areas. Since the MODIS SCF algorithm was developed under ideal conditions, SCF detection is also evaluated for varying conditions of vegetation, elevation, cloud cover and air temperature. Finally, using a direct insertion data assimilation approach, pixel averaged MODIS SCF observations are shown to improve modeled snowpack conditions over the single pixel observations due to the smoothing of more error‐prone observations and more accurately snow‐classified pixels. Copyright © 2012 John Wiley & Sons, Ltd.     

Differing controls on river‐ and lake‐water hydrogen and oxygen isotopic values in the western United States

Surface water oxygen and hydrogen isotopic values are commonly used as proxies of precipitation isotopic values to track modern hydrologic processes while proxies of water isotopic values preserved in lake and river sediments are used for paleoclimate and paleoaltimetry studies. Previous work has been able to explain variability in USA river‐water and meteoric‐precipitation oxygen isotope variability with geographic variables. These studies show that in the western United States, river‐water isotopic values are depleted relative to precipitation values. In comparison, the controls on lake‐water isotopic values are not well constrained. It has been documented that western United States lake‐water input values, unlike river water, reflect the monthly weighted mean isotopic value of precipitation. To understand the differing controls on lake‐ and river‐water isotopic values in the western United States, we examine the seasonal distribution of precipitation, evaporation and snowmelt across a range of seasonality regimes. We generate new predictive equations based on easily measured factors for western United States lake‐water, which are able to explain 69–63% of the variability in lake‐water hydrogen and oxygen isotopic values. In addition to the geographic factors that can explain river and precipitation values, lake‐water isotopic values need factors related to local hydrologic and climatic characteristics to explain variability. Study results suggest that the spring snowmelt runs off the landscape via rivers and streams, depleting river and stream‐water isotopic values. By contrast, lakes receive seasonal contributions of precipitation in proportion to the seasonal fraction of total annual precipitation within their watershed. Climate change may alter the ratio of snow to rain fall, affecting water resource partitioning between rivers and lakes and by implication of groundwater. Paleolimnological studies must account for the multiple drivers of water isotopic values; likewise, studies based on the isotopic composition of fossil material need to distinguish between species that are associated with rivers versus lakes. Copyright © 2010 John Wiley & Sons, Ltd.     

Diagnosing snow accumulation errors in a rain‐snow transitional environment with snow board observations

Diagnosing the source of errors in snow models requires intensive observations, a flexible model framework to test competing hypotheses, and a methodology to systematically test the dominant snow processes. We present a novel process‐based approach to diagnose model errors through an example that focuses on snow accumulation processes (precipitation partitioning, new snow density, and snow compaction). Twelve years of meteorological and snow board measurements were used to identify the main source of model error on each snow accumulation day. Results show that modeled values of new snow density were outside observational uncertainties in 52% of days available for evaluation, while precipitation partitioning and compaction were in error 45% and 16% of the time, respectively. Precipitation partitioning errors mattered more for total winter accumulation during the anomalously warm winter of 2014–2015, when a higher fraction of precipitation fell within the temperature range where partition methods had the largest error. These results demonstrate how isolating individual model processes can identify the primary source(s) of model error, which helps prioritize future research.     

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.     

Estimation of winter leaf area index and sky view fraction for snow modelling in boreal coniferous forests: consequences on snow mass and energy balance

Leaf area index (LAI) and canopy coverage are important parameters when modelling snow process in coniferous forests, controlling interception and transmitting radiation. Estimates of LAI and sky view factor show large variability depending on the estimation method used, and it is not clear how this is reflected in the calculated snow processes beneath the canopy. In this study, the winter LAI and sky view fraction were estimated using different optical and biomass‐based approximations in several boreal coniferous forest stands in Fennoscandia with different stand density, age and site latitude. The biomass‐based estimate of LAI derived from forest inventory data was close to the values derived from the optical measurements at most sites, suggesting that forest inventory data can be used as input to snow hydrological modelling. Heterogeneity of tree species and site fertility, as well as edge effects between different forest compartments, caused differences in the LAI estimates at some sites. A snow energy and mass balance model (SNOWPACK) was applied to detect how the differences in the estimated values of the winter LAI and sky view fraction were reflected in simulated snow processes. In the simulations, an increase in LAI and a decrease in sky view fraction changed the snow surface energy balance by decreasing shortwave radiation input and increasing longwave radiation input. Changes in LAI and sky view fraction affected directly snow accumulation through altered throughfall fraction and indirectly snowmelt through the changed surface energy balance. Changes in LAI and sky view fraction had a greater impact on mean incoming radiation beneath the canopy than on other energy fluxes. Snowmelt was affected more than snow accumulation. The effect of canopy parameters on evaporation loss from intercepted snow was comparable with the effect of variation in governing meteorological variables such as precipitation intensity and air temperature. Copyright © 2012 John Wiley & Sons, Ltd.