The design of post‐mining landscapes using geomorphic principles

Nature can provide analogues for post‐mining landscapes in terms of landscape stability and also in terms of the rehabilitated structure ‘blending in’ with the surrounding undisturbed landscape. In soil‐mantled landscapes, hillslopes typically have a characteristic pro?le that has a convex upper hillslope pro?le with a concave pro?le lower down the slope. In this paper hillslope characteristic form is derived using the area–slope relationship from pre‐mining topography at two sites in Western Australia. Using this relationship, concave hillslope pro?les are constructed and compared to linear hillslopes in terms of sediment loss using the SIBERIA erosion model. It is found that concave hillslopes can reduce sediment loss by up to ?ve times that of linear slopes. Concave slopes can therefore provide an alternative method for the construction of post‐mining landscapes. An understanding of landscape geomorphological properties and the use of erosion models can greatly assist in the design of post‐mining landscapes. Copyright © 2003 John Wiley & Sons, Ltd.     

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

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

Coupling the WRF model with a temperature index model based on remote sensing for snowmelt simulations in a river basin in the Altay Mountains,north‐west China

Snowmelt water is an important freshwater resource in the Altay Mountains in north‐west China; however, warming climate and rapid spring snowmelt can cause floods that endanger both public and personal property and safety. This study simulates snowmelt in the Kayiertesi River catchment using a temperature index model based on remote sensing coupled with high‐resolution meteorological data obtained from National Centers for Environmental Prediction (NCEP) reanalysis fields that were downscaled using the Weather Research Forecasting model and then bias corrected using a statistical downscaled model. Validation of the forcing data revealed that the high‐resolution meteorological fields derived from the downscaled NCEP reanalysis were reliable for driving the snowmelt model. Parameters of the temperature index model based on remote sensing were calibrated for spring 2014, and model performance was validated using Moderate Resolution Imaging Spectroradiometer snow cover and snow observations from spring 2012. The results show that the temperature index model based on remote sensing performed well, with a simulation mean relative error of 6.7% and a Nash–Sutcliffe efficiency of 0.98 in spring 2012 in the river of Altay Mountains. Based on the reliable distributed snow water equivalent simulation, daily snowmelt run‐off was calculated for spring 2012 in the basin. In the study catchment, spring snowmelt run‐off accounts for 72% of spring run‐off and 21% of annual run‐off. Snowmelt is the main source of run‐off for the catchment and should be managed and utilized effectively. The results provide a basis for snowmelt run‐off predictions, so as to prevent snowmelt‐induced floods, and also provide a generalizable approach that can be applied to other remote locations where high‐density, long‐term observational data are lacking. Copyright © 2016 John Wiley & Sons, Ltd.     

Restoring earth surface processes through landform design. A 13‐year monitoring of a geomorphic reclamation model for quarries on slopes

The application of geomorphic principles to land reclamation after surface mining has been reported in the literature since the mid‐1990s, mostly from Australia, Canada and the USA. This paper discusses the reclamation problems of contour mining and quarries on slopes, where steep gradients are prone to both mass movement and water erosion. To address these problems simultaneously, a geomorphic model for reclaiming surface mined slopes is described. Called the ‘highwall–trench–concave slope’ model, it was first applied in the 1995 reclamation of a quarry on a slope (La Revilla) in Central Spain. The geomorphic model does not reproduce the original topography, but has two very different sectors and objectives: (i) the highwall–trench sector allows the former quarry face to evolve naturally by erosion, accommodating fallen debris by means of a trench constructed at the toe of the highwall; (ii) the concave‐slope base sector, mimicking the landforms of the surrounding undisturbed landscape, promotes soil formation and the establishment of self‐sustaining, functional ecosystems in the area protected from sedimentation by the trench. The model improves upon simple topographic reconstruction, because it rebuilds the surficial geology architecture and facilitates re‐establishment of equilibrium slopes through the management and control of geomorphic processes. Thirteen years of monitoring of the geomorphic and edaphic evolution of La Revilla reclaimed quarry confirms that the area is functioning as intended: the highwall is backwasting and material is accumulating at the trench, permitting the recovery of soils and vegetation on the concave slope. However, the trench is filling faster than planned, which may lead to run‐off and sedimentation on the concave slope once the trench is full. The lesson learned for other scenarios is that the model works well in a two‐dimensional scheme, but requires a three‐dimensional drainage management, breaking the reclaimed area into several watersheds with stream channels. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluating restoration of man‐made slopes: a threshold approach balancing vegetation and rill erosion

The management of reclaimed slopes derived from industrial and civil activities (e.g. surface mining and road construction) requires the development of practical stability analysis approaches that integrate the processes and mechanisms that rule the dynamics of these ubiquitous emerging ecosystems. This work describes a new modelling approach focused on stability analysis of water‐limited reclaimed slopes, where interactive relationships between rill erosion and vegetation regulate ecosystem stability. Our framework reproduces two main groups of possible trends along the temporal evolution of reclaimed slopes: successful trends, characterized by widespread vegetation development and the effective control of rill erosion processes; and gullying trends, characterized by the progressive loss of vegetation and a sharp logistic increase in erosion rates. Furthermore, this analytical approach allows the determination of threshold values for the state variables (i.e. vegetation cover and rill erosion) that drive the system's stability, facilitating the identification of critical situations that require specific human intervention (e.g. revegetation or, in very problematic cases, revegetation combined with rill network destruction) to ensure the long‐term sustainability of the restored ecosystem. The application of our threshold analysis framework in Mediterranean‐dry reclaimed slopes derived from surface coal mining (the Teruel coalfield in central‐eastern Spain) showed a good field‐based performance. Therefore, we believe that this model is a valuable contribution for the management of water‐limited reclaimed systems, including those associated with rill erosion, as it provides a tool for the evaluation of restoration success and can play an important role in decision‐making during ecosystem restoration in severely disturbed landscapes. Copyright © 2011 John Wiley & Sons, Ltd.     

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

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

Relating snow dynamics and balsam fir stand characteristics,Montmorency Forest,Quebec

Snow water equivalent was measured during three springs on north‐ and south‐exposed sites representing a range of stand structure and development stages of Quebec's balsam fir forest. Maximum snow water equivalent of the season, mean seasonal snowmelt rate, snowmelt season duration and total snowmelt season degree‐day factor were related to canopy height, canopy density, light interception fraction and basal area of the stands using random coefficient models. Seasonal mean snowmelt rate was better explained by stand characteristics (R² from 0·41 to 0·61) than was maximum snow water equivalent (R² from 0·08 to 0·23). The best relationship was found with light interception, which explained 61% of snowmelt rate variability between stands. These relationships were not significantly affected by stand aspect (Pr ≥ S = 0·14 or higher), as snow dynamics seemed less dependent on aspect than on stand characteristics. Snowmelt recovery rates could be used by forest planners to establish an acceptable time step for the harvesting of different parts of a watershed in order to prevent peak flow augmentations. Copyright © 2005 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.     

Quantifying the role of the snowpack in generating water available for run‐off during rain‐on‐snow events from snow pillow records

Rain‐on‐snow events have generated major floods around the world, particularly in coastal, mountainous regions. Most previous studies focused on a limited number of major rain‐on‐snow events or were based primarily on model results, largely due to a lack of long‐term records from lysimeters or other instrumentation for quantifying event water balances. In this analysis, we used records from five automated snow pillow sites in south coastal British Columbia, Canada, to reconstruct event water balances for 286 rain‐on‐snow events over a 10‐year period. For large rain‐on‐snow events (event rainfall >40 mm), snowmelt enhanced the production of water available for run‐off (WAR) by approximately 25% over rainfall alone. For smaller events, a range of antecedent and meteorological factors influenced WAR generation, particularly the antecedent liquid water content of the snowpack. Most large events were associated with atmospheric rivers. Rainfall dominated WAR generation during autumn and winter events, whereas snowmelt dominated during spring and summer events. In the majority of events, the sensible heat of rain contributed less than 10% of the total energy consumed by snowmelt. This analysis illustrated the importance of understanding the amount of rainfall occurring at high elevations during rain‐on‐snow events in mountainous regions.     

Observations of urban snow properties in Calgary,Canada

Urban winter hydrology has garnered very little attention owing to the general notion that high‐intensity rainfalls are the major flood‐generating events in urban areas. As a result, few efforts have been made to research urban snow and its melt characteristics. This study investigates the characteristics of urban snow that differentiate it from rural snow, and makes recommendations for incorporating these characteristics into an urban snowmelt model. A field study was conducted from the fall of 2001 to the spring of 2002 in the city of Calgary, Canada. Snow depths and densities, soil moisture, soil temperature, snow albedo, net radiation, snow evaporation, and surface temperature were measured at several locations throughout the winter period. The combination of urban snow removal practices and the physical elements that exist in urban areas were found to influence the energy balance of the snowpack profoundly. Shortwave radiation was found to be the main source of energy for urban snow; as a consequence, the albedo of urban snow is a very important factor in urban snowmelt modelling. General observations lead to the classification of snow as one of four types: snow piles, snow on road shoulders, snow on sidewalk edges, and snow in open areas. This resulted in the development of four separate functions for the changing snow albedo values. A study of the frozen ground conditions revealed that antecedent soil moisture conditions had very little impact on frozen ground, and thus frozen ground very nearly always acts as a near impervious area. Improved flood forecasting for urban catchments in cold regions can only be achieved with accurate modelling of urban winter runoff that involves the energy balance method, incorporating snow redistribution and urban snow‐cover characteristics, and using small time steps. Copyright © 2004 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.     

Watershed‐scale controls on snow accumulation in a small montane watershed,southwestern Alberta,Canada

Snow accumulation in mountain headwater basins is a major water source, particularly in semi‐arid environments such as southern Alberta where water resources are stressed and snowmelt supplies more than 80% of downstream runoff. Relationships between landscape predictor variables and snow water equivalent (SWE) were quantified by combining field and LiDar measurements with classification and regression tree analysis over two winter seasons (2010 and 2011) in a small, montane watershed. 2010 was a below average snow accumulation year, while 2011 was well above normal. In both the field and regression tree data, elevation was the dominant control on snow distribution in both years, although snow distribution was driven by melt processes in 2010 and accumulation processes in 2011. The importance of solar radiation and wind exposure was represented in the regression trees in both years. The regression trees also noted the lower importance of canopy closure, slope, and aspect, which was not observed in the field data. This technique could provide an additional method of forecasting annual water supply from melting snow. However, further research is required to address the lack of data collected above treeline, to provide a full‐basin estimate of SWE. Copyright © 2012 John Wiley & Sons, Ltd.     

Factors and interactions controlling infiltration,runoff, and soil loss at the microscale in a patchy Mediterranean semiarid landscape

In semiarid ecosystems, the transfer of water, sediments, and nutrients from bare to vegetated areas is known to be crucial to ecosystem functioning. Rainfall simulation experiments were performed on bare‐soil and vegetated surfaces, on both wet and dry soils, in semiarid shrub‐steppe landscapes of SE Spain to investigate the spatial and temporal factors and interactions that control the fine‐scale variation in water infiltration, runoff and soil loss, and hence the water and sediment flows in these areas. Three types of shrub‐steppe landscapes varying in plant community and physiography, and four types of plant patches (oak shrub, subshrub, tussock grass, and short grass mixed with chamaephytes) were studied. Higher infiltration and lower runoff and soil loss were measured on vegetation patches than on bare soils, for both dry and wet conditions. The oak‐shrub patches produced no runoff, while the subshrub patches showed the highest runoff and soil loss. Despite these differences among patch types, the influence of vegetation patch type on the variables analysed was not significant. The response of bare soil surfaces clearly varied between landscape types, yet the differences were only relevant under dry soil conditions. Stone cover, particularly the cover of embedded stones, and crust cover, were the key explanatory variables for the hydrological behaviour of bare soils. The study documents quantitatively how bare soils and vegetation patches function as runoff sources and runoff sinks, respectively, for a wide range of soil moisture conditions, and illustrates that landscape‐type effects on bare‐soil runoff sources may also exert an important control on the site hydrology, while the role of the vegetation patch type is less important. The effects of the control factors are modulated by antecedent soil moisture, with dry soils showing the most contrasting soil water infiltration between landscapes and surface types. Copyright © 2009 John Wiley & Sons, Ltd.     

A comparison of MODIS and NOHRSC snow‐cover products for simulating streamflow using the Snowmelt Runoff Model

Remote sensing is an important source of snow‐cover extent for input into the Snowmelt Runoff Model (SRM) and other snowmelt models. Since February 2000, daily global snow‐cover maps have been produced from data collected by the Moderate Resolution Imaging Spectroradiometer (MODIS). The usefulness of this snow‐cover product for streamflow prediction is assessed by comparing SRM simulated streamflow using the MODIS snow‐cover product with streamflow simulated using snow maps from the National Operational Hydrologic Remote Sensing Center (NOHRSC). Simulations were conducted for two tributary watersheds of the Upper Rio Grande basin during the 2001 snowmelt season using representative SRM parameter values. Snow depletion curves developed from MODIS and NOHRSC snow maps were generally comparable in both watersheds: satisfactory streamflow simulations were obtained using both snow‐cover products in larger watershed (volume difference: MODIS, 2·6%; NOHRSC, 14·0%) and less satisfactory streamflow simulations in smaller watershed (volume difference: MODIS, −33·1%; NOHRSC, −18·6%). The snow water equivalent (SWE) on 1 April in the third zone of each basin was computed using the modified depletion curve produced by the SRM and was compared with in situ SWE measured at Snowpack Telemetry sites located in the third zone of each basin. The SRM‐calculated SWEs using both snow products agree with the measured SWEs in both watersheds. Based on these results, the MODIS snow‐cover product appears to be of sufficient quality for streamflow prediction using the SRM in the snowmelt‐dominated basins. Copyright © 2005 John Wiley & Sons, Ltd.     

Snowmelt runoff and soil moisture dynamics on steep subalpine hillslopes

Snowmelt water supplies streamflow and growing season soil moisture in mountain regions, yet pathways of snowmelt water and their effects on moisture patterns are still largely unknown. This study examined how flow processes during snowmelt runoff affected spatial patterns of soil moisture on two steep sub‐alpine hillslope transects in Rocky Mountain National Park, CO, USA. The transects have northeast‐facing and east‐facing aspects, and both extend from high‐elevation bedrock outcrops down to streams in valley bottoms. Spatial patterns of both snow depth and near‐surface soil moisture were surveyed along these transects in the snowmelt and summer seasons of 2008–2010. To link these patterns to flow processes, soil moisture was measured continuously on both transects and compared with the timing of discharge in nearby streams. Results indicate that both slopes generated shallow lateral subsurface flow during snowmelt through near‐surface soil, colluvium and bedrock fractures. On the northeast‐facing transect, this shallow subsurface flow emerged through mid‐slope seepage zones, in some cases producing saturation overland flow, whereas the east‐facing slope had no seepage zones or overland flow. At the hillslope scale, earlier snowmelt timing on the east‐facing slope led to drier average soil moisture conditions than on the northeast‐facing slope, but within hillslopes, snow patterns had little relation to soil moisture patterns except in areas with persistent snow drifts. Results suggest that lateral flow and exfiltration processes are key controls on soil moisture spatial patterns in this steep sub‐alpine location. Copyright © 2014 John Wiley & Sons, Ltd.     

Influences of vegetation disturbance on hydrogeomorphic response following wildfire

Quantifying the linkages between vegetation disturbance by fire and the changes in hydrologic processes leading to post‐fire erosional response remains a challenge. We measured the influence of fire severity, defined as vegetation disturbance (using a satellite‐derived vegetation disturbance index, VDI), landscape features that organize hydrologic flow pathways (relief and elongation ratios), and pre‐fire vegetation type on the probability of the occurrence of post‐fire gully rejuvenation (GR). We combined field surveys across 270 burned low‐order catchments (112 occurrences of GR) and geospatial analysis to generate a probabilistic model through logistic regression. VDI alone discriminated well between catchments where GR did and did not occur (area under the curve = 0.78, model accuracy = 0.72). The strong effect of vegetation disturbance on GR suggests that vegetation exerts a primary influence on the occurrence of infiltration excess run‐off and post‐fire erosion and that major gully erosion will not occur until fire consumes aboveground biomass. Other topographic and local factors also influenced GR response, including catchment elongation, per cent pre‐fire shrub, mid‐slope riparian vegetation, armoured headwaters, firehose effects, and concentration of severe burn in source areas. These factors highlight the need to consider vegetation effects in concert with local topography and site conditions to understand the propensity for flow accumulation leading to GR. We present a process‐based conceptual hydrologic model where vegetation loss from fire decreases rainfall attenuation and surface roughness, leading to accelerated flow accumulation and erosion; these effects are also influenced by interactions between fire severity and landscape structure. The VDI metric provides a continuous measure of vegetation disturbance and, when placed in a hydrologic context, may improve quantitative analysis of burned‐area susceptibility to erosive rainfall, hazard prediction, ecological effects of fire, landform evolution, and sensitivity to climate change. Copyright © 2015 John Wiley & Sons, Ltd.     

Model complexity and data requirements in snow hydrology: seeking a balance in practical applications

We investigate the problem of balancing model complexity and input data requirements in snow hydrology. For this purpose, we analyze the performance of two models of different complexity in estimating variables of interest in snow hydrology applications. These are snow depth, bulk snow density, snow water equivalent and snowmelt run‐off. We quantify the differences between data and model prediction using 18 years of measurements from an experimental site in the French Alps (Col de Porte, 1325 m AMSL). The models involved in this comparison are a one‐layer temperature‐index model (HyS) and a multilayer model (Crocus). Results show that the expected loss in performance in the one‐layer temperature‐index model with respect to the multilayer model is low when considering snow depth, snow water equivalent and bulk snow density. As for run‐off, the comparison returns less clear indications for identification of a balance. In particular, differences between the models' prediction and data with an hourly resolution are higher when considering the Crocus model than the HyS model. However, Crocus is better at reproducing sub‐daily cycles in this variable. In terms of daily run‐off, the multilayer physically based model seems to be a better choice, while results in terms of cumulative run‐off are comparable. The better reproduction of daily and sub‐daily variability of run‐off suggests that use of the multilayer model may be preferable for this purpose. Variation in performance is discussed as a function of both the calibration solution chosen and the time of year. Copyright © 2016 John Wiley & Sons, Ltd.     

The effect of landform variation on vegetation patterning and related sediment dynamics

Semi‐arid ecosystems are often spatially self‐organized in typical patterns of vegetation bands with high plant cover interspersed with bare soil areas, also known as ‘tiger bush’. In modelling studies, most often, straight planar slopes were used to analyse vegetation patterning. The effect of slope steepness has been investigated widely, and some studies investigated the effects of microtopography and hillslope orientation. However, at the larger catchment scale, the overall form of the landscape may affect vegetation patterning and these more complex landscapes are much more prevalent than straight slopes. Hence, our objective was to determine the effect of landform variation on vegetation patterning and sediment dynamics. We linked two well‐established models that simulate (a) plant growth, death and dispersal of vegetation, and (b) erosion and sedimentation dynamics. The model was tested on a straight planar hillslope and then applied to (i) a set of simple synthetic topographies with varying curvature and (ii) three more complex, real‐world landscapes of distinct morphology. Results show banded vegetation patterning on all synthetic topographies, always perpendicular to the slope gradient. Interestingly, we also found that movement of bands – a debated phenomenon – seems to be dependent on curvature. Vegetation banding was simulated on the slopes of the alluvial fan and along the valley slopes of the dissected and rolling landscapes. In all landscapes, local valleys developed a full vegetation cover induced by water concentration, which is consistent with observations worldwide. Finally, banded vegetation patterns were found to reduce erosion significantly as compared to other vegetation configurations. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Hydrologic flowpaths during snowmelt in forested headwater catchments under differing winter climatic and soil frost regimes

Changes in hydrologic flowpaths have important impacts on the timing, magnitude and hydrochemistry of run‐off during snowmelt in forested catchments, but how flowpaths are affected by variation in winter climate and the irregular presence of soil frost remains poorly understood. The depth and extent of soil frost may be expected to increase as snowpack decreases or develops later because of climate change. In this study, we used end‐member mixing analysis to determine daily contributions of snow, forest floor soil water and groundwater to stream run‐off during snowmelt under different soil frost regimes resulting from interannual and elevational variation at the Hubbard Brook Experimental Forest in New Hampshire, USA. We observed greater routing of run‐off through forest floor flowpaths during early snowmelt in 2011, when the snowpack was deep and soil frost was minimal, compared with the early snowmelt in 2012 under conditions of deep and extensive soil frost. The results indicate that widespread soil frost that penetrated the depth of the forest floor decreased the flow signal through the shallowest subsurface flowpaths, but did not reduce overall infiltration of melt waters, as the contribution from the snow‐precipitation end‐member was similar under both conditions. These results are consistent with development of granular soil frost which permits vertical infiltration of melt waters, but either reduces lateral flow in the forest floor or prevents the solute exchange that would produce the typical chemical signature of shallow subsurface flowpaths in streamwater. Copyright © 2016 John Wiley & Sons, Ltd.