MEASUREMENT OF FROST-INDUCED SNOWMELT RUNOFF IN A FOREST SOIL

Snow interception in a coniferous stand leads to considerable short-range variability in snowcover depth, which in turn affects the water and heat regime of the soil. To study the coupling between snow accumulation, frost penetration, and hydrological response, plot-scale experiments were conducted in a subalpine spruce forest. The stony, sandy–loamy Spodosol was highly permeable and had an organic layer of 5–15 cm thickness. Within two plots, one underneath a tree crown and one in a canopy gap, we measured near-surface runoff, soil temperature, and liquid water content. Snow and frost depths varied more in space than between two winter periods at given locations. Frost penetration was greater near the trunk, where a higher portion of snowmelt water drained downslope close to the surface than in the gap due to frost-induced reduction of infiltration. In both years, the spring snowmelt occurred over two distinct periods. During the first snowmelt, the water percolated primarily through the frozen layer and part of it probably refroze within the frozen layer, thereby raising the total water and ice content. During the second event, near-surface runoff was more pronounced.     

Estimating winter evaporation in boreal forests with operational snow course data

Snow course measurements from 11 sites located in eastern and northern Finland were used to estimate the total interception evaporation of a winter season for different forest categories. We categorized the sites based on forest density and tree species. Results showed that interception loss from gross precipitation increased with forest density and approached 30% for a forest with the highest density class. Interception loss for the most common forest density class was 11%. Interception losses were slightly larger in spruce forests than in pine, deciduous, or mixed forests. We provide suggestions as to how to design snow surveys to estimate wintertime interception evaporation better. Rough terrain and transition zones between forest and open areas should be avoided. Since evaporation fraction was strongly dependent on tree crown characteristics, snow course data should include direct estimates of canopy closure. Qualitative observations made by different observers should be given a reference frame to ensure comparability of records from different sites. Copyright © 2003 John Wiley & Sons, Ltd.     

Representation of canopy snow interception,unloading and melt in a parsimonious snowmelt model

Recent improvements in the Utah Energy Balance (UEB) snowmelt model are focused on snow–vegetation–atmosphere interactions to understand how different types of vegetation affect snow processes in the mountains of Western USA. This work presents field work carried out in the Rocky Mountains of Northern Utah to evaluate new UEB model algorithms that represent the processes of canopy snow interception, sublimation, mass unloading and melt. Four years' continuous field observations showed generally smaller accumulations of snow beneath the forest canopies in comparison with open (sage and grass) areas, a difference that is attributed to interception and subsequent sublimation and redistribution of intercepted snow by wind, much of it into surrounding open areas. Accumulations beneath the denser forest (conifer) canopies were found to be less than the accumulation beneath the less dense forest (deciduous) canopies. The model was able to represent the accumulation of snow water equivalent in the open and beneath the deciduous forest quite well but without accounting for redistribution tended to overestimate the snow water equivalent beneath the conifer forest. Evidence of redistribution of the intercepted snow from the dense forest (i.e. conifer forest) to the adjacent area was inferred from observations. Including a simple representation of redistribution in the model gave satisfactory prediction of snow water equivalent beneath the coniferous forest. The simulated values of interception, sublimation and unloading were also compared with previous studies and found in agreement. Copyright © 2013 John Wiley & Sons, Ltd.     

MID-WINTER STEMFLOW DRAINAGE FROM BIGTOOTH ASPEN (POPULUS GRANDIDENTATA MICHX.) IN CENTRAL MASSACHUSETTS

Under winter conditions, stemflow drainage in forested ecosystems is often assumed to be a negligible component of the hydrological cycle. This paper reports on mid-winter stemflow drainage from the broadleaved deciduous tree species Populus grandidentata. Stemflow volumes from this species at air temperatures of < 0°C were found to be comparable to rainfall-generated stemflow during summer. Over the three-month period January–March 1993, stemflow ranged from 5.4 to 9.9% of the incident gross precipitation. Expressed as depth equivalents per unit trunk basal area, these stemflow inputs ranged from 1.8 to 4.9 m. These concentrated mid-winter inputs of liquid water to the bases of canopy trees were attributable to: (1) snow interception by the leafless woody frame of each tree; (2) snow retention by glazed ice precipitation associated with the snowfall event; (3) increased temperature at the bark/snow interface caused by the low albedo of the bark tissue; and (4) convergence of snowmelt drainage from steeply inclined upthrust primary branches. The hydrological and ecological significance of liquid water inputs to the forest floor under sub-zero conditions are discussed. © 1997 by John Wiley & Sons Ltd.     

Routing of canopy drip in the snowpack below a spruce crown

Snow that is retained by a forest canopy may either sublimate or evaporate directly from the crown or drop as snow clumps or meltwater to the ground. Redistributed snow and meltwater affect the snow structure and prevent the formation of mechanically weak layers, which is the prerequisite for avalanche formation in forests. In this paper we describe the results of dye tracer experiments conducted in a subalpine forest near Davos, Switzerland. Before a snowfall event we stained snow‐free branches of a spruce with a dye tracer solution. After snowfall the coloured meltwater dripping from the branches down on to the snowpack stained the percolation pathways of the meltwater in the snowpack. Photographs of the snow profiles indicate that the meltwater seeped almost vertically through the isothermal snowpack to the soil surface not exceeding the projected crown edge. Meltwater of different events moves along different preferential flow channels in the snow suggesting that old channels are not non‐conducting and additional meltwater fronts create new channels. Copyright © 1999 John Wiley & Sons, Ltd.     

The simulation of heat and water exchange in the boreal spruce forest by the land-surface model SWAP

All previous versions of a physically based land-surface model SWAP have assumed for simplicity that vegetation is fully covered by snow during the cold season. Such assumption is reasonable only for the regions dominated by short vegetation or for warm climates where snow processes are absent. The major goals of this paper are (i) modification of the latest version of SWAP by incorporation of tall vegetation into the cold-season parameterizations to make the model applicable for simulating heat and water transfer within a boreal forest biome and (ii) validation of the modified version using the data from a forested catchment located in the boreal zone. Modification of SWAP required to parameterize radiative and turbulent exchange between the forest crown and forest floor, partitioning snowfall between interception by the canopy (in doing so, snow interception differs from rain interception) and falling to the ground, formation of snow cover on the forest crown and forest floor including snow accumulation (both in solid and liquid fractions), snow evaporation, and snowmelt. The advanced model was validated using a set of hydrometeorological data measured during 18 years (1966–1983) at the Tayozhniy catchment (covered by boreal spruce forest), Valdai, Russia. Simulations of annual and monthly snow/rain interception, daily runoff at the catchment outlet, snow density, snow depth, snow water equivalent, soil water storage in three layers (0–20, 0–50 and 0–100 cm), and monthly evapotranspiration from the catchment were compared with observations. Analysis of the results of validation shows that the new version of the model SWAP reproduces the heat and water exchange processes occurring in mid-latitude boreal forest quite reasonable.     

Rainfall interception by an isolated evergreen oak tree in a Mediterranean savannah

Redistribution of ground‐level rainfall and interception loss by an isolated Quercus ilex tree were measured over 2 years in a Mediterranean oak savannah. Stemflow, meteorological variables and sap flow were also monitored. Rainfall at ground level was measured by a set of rain‐gauges located in a radial layout centred on the tree trunk and extending beyond the crown limits. Interception loss was computed as the difference between the volume of rainwater that would reach the ground in the absence of the tree and the volume of water that actually fell on the ground sampling area (stemflow included). This procedure provided correct interception loss estimates, irrespective of rainfall inclination. Results have shown a clear non‐random spatial distribution of ground‐level rainfall, with rainwater concentrations upwind beneath the crown and rain‐shadows downwind. Interception loss amounted to 22% of gross rainfall, per unit of crown‐projected area. Stand interception loss, per unit of ground area, was only 8% of gross rainfall and 28% of tree evapotranspiration. These values reflect the low crown cover fraction of the stand (0·39) and the specific features of the Mediterranean rainfall regime (predominantly with few large storms). Nevertheless, it still is an important component of the water balance of these Mediterranean ecosystems. Copyright © 2005 John Wiley & Sons, Ltd.     

Assessing the effects of post‐pine beetle forest litter on snow albedo

The effect of forest litter on snow surface albedo has been subject to limited study, mainly in the hardwood‐dominated forests of the northeastern United States. Given the recent pine beetle infestation in Western North America and associated increases in litter production, this study examines the effects of forest litter on snow surface albedo in the coniferous forests of south‐central British Columbia. Measured changes in canopy transmittance provide an indication of canopy loss or total litterfall over the winter of 2007–2008. Relationships between percent litter cover, an index of albedo, snow depth, and snow ablation during the 2008 melt season are compared between a mature, young, and clearcut coniferous stand. Results indicate a strong feedback effect between canopy loss and subsequent enhanced shortwave transmittance, and litter accumulation on the snow surface from that canopy loss. However, this relationship is confounded by other variables concurrently affecting albedo. While results suggest that a relatively small percent litter cover can have a significant effect on albedo and ablation, further research is underway to extract the litter signal from that of other factors affecting albedo, particularly snow depth. Copyright © 2010 John Wiley & Sons, Ltd.     

Snow Mass Quantification and Avalanche Victim Search by Ground Penetrating Radar

Ground penetrating radar (GPR) systems can be used in many applications of snow and ice research. The information from the GPR is used to identify and interpret layers, objects and different structures in the snow. A commercially available GPR system was further developed to work in the rough environment of snow and ice. The applied GPR is a 900 MHz system that easily reaches snow depths of up to 10 meters. The system was calibrated in the course of several manual snow depth measurements during each survey. The depth resolution depends on the snow type and is around ±0.1 m. The GPR system is carried alongside a line of interest and is triggered by an odometer wheel at regular adjustable steps. All equipment is mounted in a sledge and is pulled by a snowmobile over the snow surface. This setup allows for an efficient coverage of several kilometers of terrain profiles. The radar profiles give a real time two-dimensional impression of structures and objects and the interface between snow and the underlying ground. The actual radar profile is shown on a screen on the sledge allowing the immediate marking of objects and structures. During the past three years the instrument was successfully used for the study of snow distributions, for the detection of glacier crevasses under the snow cover, and for the search of avalanche victims in avalanche debris. The results show the capability of the instrument to detect persons and objects in the snow cover. In the future, this device may be a new tool for avalanche rescue operations. Today, the size and weight of the system prevents the accessing of very steep slopes and areas not accessible to snowmobiles. Further developments will decrease the size of the system and make it a valuable tool to quantify snow masses in avalanche release zones and run-out areas.     

Quantifying snow water equivalent using terrestrial ground penetrating radar and unmanned aerial vehicle photogrammetry

This study demonstrates the potential value of a combined unmanned aerial vehicle (UAV) Photogrammetry and ground penetrating radar (GPR) approach to map snow water equivalent (SWE) over large scales. SWE estimation requires two different physical parameters (snow depth and density), which are currently difficult to measure with the spatial and temporal resolution desired for basin-wide studies. UAV photogrammetry can provide very high-resolution spatially continuous snow depths (SD) at the basin scale, but does not measure snow densities. GPR allows nondestructive quantitative snow investigation if the radar velocity is known. Using photogrammetric snow depths and GPR two-way travel times (TWT) of reflections at the snow-ground interface, radar velocities in snowpack can be determined. Snow density (RSN) is then estimated from the radar propagation velocity (which is related to electrical permittivity of snow) via empirical formulas. A Phantom-4 Pro UAV and a MALA GX450 HDR model GPR mounted on a ski mobile were used to determine snow parameters. A snow-free digital surface model (DSM) was obtained from the photogrammetric survey conducted in September 2017. Then, another survey in synchronization with a GPR survey was conducted in February 2019 whilst the snowpack was approximately at its maximum thickness. Spatially continuous snow depths were calculated by subtracting the snow-free DSM from the snow-covered DSM. Radar velocities in the snowpack along GPR survey lines were computed by using UAV-based snow depths and GPR reflections to obtain snow densities and SWEs. The root mean square error of the obtained SWEs (384 mm average) is 63 mm, indicating good agreement with independent SWE observations and the error lies within acceptable uncertainty limits.     

Influence of canopy density on snow distribution in a temperate mountain range

We analyse spatial variability and different evolution patterns of snowpack in a mixed beech–fir stand in the central Pyrenees. Snow depth and density were surveyed weekly along six transects of contrasting forest cover during a complete accumulation and melting season; we also surveyed a sector unaffected by canopy cover. Forest density was measured using the sky view factor (SVF) obtained from digital hemispherical photographs. During periods of snow accumulation and melting, noticeable differences in snow depth and density were found between the open site and those areas covered by forest canopy. Principal component analysis provided valuable information in explaining these observations. The results indicate a high variability in snow accumulation within forest areas related to differences in canopy density. Maximum snow water equivalent (SWE) was reduced by more than 50% beneath dense canopies compared with clearings, and this difference increased during the melting period. We also found significant temporal variations: when melting began in sectors with low SVF, most of the snow had already thawed in areas with high SVF. However, specific conditions occasionally produced a different response of SWE to forest cover, with lower melting rates observed beneath dense canopies. The high values of correlation coefficients for SWE and SVF (r > 0·9) indicate the reliability of predicting the spatial distribution of SWE in forests when only a moderate number of observations are available. Digital hemispherical photographs provide an appropriate tool for this type of analysis, especially for zenith angles in the range 35–55 . Copyright © 2007 John Wiley & Sons, Ltd.     

Infiltration-excess caused by Stemflow in a cyclone-prone tropical rainforest

Tropical rainforest canopy trees that have large projected areas of upwardly inclined branches are capable of funnelling large volumes of rainwater down their trunks. During periods of prolonged heavy rainfall on Mount Bellenden Ker in northeast Queensland, Australia, stemflow volumes were found to be as much as two orders of magnitude greater than the volume of incident rainfall expected in a rain gauge occupying an area equal to the trunk basal area. Stemflow totals ranging from 6000 to 70000 litres were generated by individual trees from 7800 mm of rainfall over two successive wet seasons. The combination of high intensity rainfall and the funnelling effect results in significant quantities of infiltration-excess at the ground surface. Stemflow fluxes as high as 31.4 cm³ min^?1 per cm² of basal area (i.e. the equivalent of 314 mm min^?1) were recorded when rainfall intensity was only 2 mm min^?1. The mean infiltration capacity of the topsoil was determined to be 6.2 mm min^?1. The areas over which the stemflow would have had to spread in order to infiltrate were computed to be as much as 3 m² around the bases of individual canopy trees. Approximations of the distances that the infiltration-excess would have travelled away from the tree bases were calculated by assuming that the infiltration area either expands radially outward in the form of an annulus or extends straight downslope from the tree base.     

The effect of single vegetation elements on wind speed and sediment transport in the Sahelian zone of Burkina Faso

Soil loss caused by wind erosion is a widespread phenomenon in the Sahelian zone of West Africa. According to Sahelian farmers, scattered vegetation standing in amongst the crop has the potential for a wind erosion control strategy. This study was conducted to study the effect of single vegetation elements on the pattern of average wind speed and sediment transport. This was done by two experiments that were carried out during the rainy seasons of 2002 and 2003 in north Burkina Faso, West Africa. Wind speeds were measured using three sonic anemometers, at a sampling frequency of 16 Hz. Sediment transport was determined by calculating the mass fluxes from 17 MWAC catchers. In this study, a shrub was defined as a vegetation element with branches until ground and a tree as a vegetation element with a distinctive trunk below a canopy. Behind shrubs wind speed near the soil surface was reduced up to approximately seven times the height of the shrub. The observed reduction in wind speed in the area where wind speed was reduced was 15 per cent on average. At the sides of the shrub, wind speed was increased, by on average 6 per cent. As the area of increase in wind speed is one‐third of the area of decrease in wind speed, the net effect of a shrub is a reduction in wind speed. A similar pattern was visible for the pattern of sediment transport around a shrub. Downwind of a shrub, sediment transport was diminished up to seven times the height of the shrub. Probably most of this material was trapped by the shrub. Trees showed a local increase of wind around the trunk, which is expected to relate to an increase in sediment transport around the trunk. Mass flux measurements of sediment transport were not made, but visual observations in the field substantiate this. Behind the canopy of a tree, a tree acts similarly to a shrub regarding its effects on average wind speed, but as a tree is generally a larger obstacle than a shrub the extent of this effect is larger than for shrubs. Thus, whereas shrubs are more effective than trees regarding their direct effect on soil loss by trapping sand particles near the soil surface, trees are more effective in affecting soil loss indirectly by reducing the wind speed downwind more effectively than shrubs. Therefore, to reduce soil loss in an area, the presence of both trees and shrubs is crucial. Copyright © 2007 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.     

A conceptual model of water yield effects from beetle‐induced tree death in snow‐dominated lodgepole pine forests

In regions of western North America with snow‐dominated hydrology, the presence of forested watersheds can significantly influence streamflow compared to areas with other vegetation cover types. Widespread tree death in these watersheds can thus dramatically alter many ecohydrologic processes including transpiration, canopy solar transmission and snow interception, subcanopy wind regimes, soil infiltration, forest energy storage and snow surface albedo. One of the more important causes of conifer tree death is bark beetle infestation, which in some instances will kill nearly all of the canopy trees within forest stands. Since 1996, an ongoing outbreak of bark beetles (Coleoptera: Scolytidae) has caused widespread mortality across more than 600,000 km² of coniferous forests in western North America, including numerous Rocky Mountain headwaters catchments with high rates of lodgepole pine (Pinus contorta) mortality from mountain pin beetle (Dendroctonous ponderosae) infestations. Few empirical studies have documented the effects of MPB infestations on hydrologic processes, and little is known about the direction and magnitude of changes in water yield and timing of runoff due to insect‐induced tree death. Here, we review and synthesize existing research and provide new results quantifying the effects of beetle infestations on canopy structure, snow interception and transmission to create a conceptual model of the hydrologic effects of MPB‐induced lodgepole pine death during different stages of mortality. We identify the primary hydrologic processes operating in living forest stands, stands in multiple stages of death and long‐dead stands undergoing regeneration and estimate the direction of change in new water yield. This conceptual model is intended to identify avenues for future research efforts. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

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

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

Snow disappearance timing is dominated by forest effects on snow accumulation in warm winter climates of the Pacific Northwest,United States

Forests modify snow processes and affect snow water storage as well as snow disappearance timing. However, forest influences on snow accumulation and ablation vary with climate and topography and are therefore subject to temporal and spatial variability. We utilize multiple years of snow observations from across the Pacific Northwest, United States, to assess forest–snow interactions in the relatively warm winter conditions characteristic of maritime and transitional maritime–continental climates. We (a) quantify the difference in snow magnitude and disappearance timing between forests and open areas and (b) assess how forest modifications of snow accumulation and ablation combine to determine whether snow disappears later in the forest or in the open. We find that snow disappearance timing at 12 (out of 14) sites ranges from synchronous in the forest and open to snow persisting up to 13 weeks longer in the open relative to a forested area. By analyzing accumulation and ablation rates up to the day when snow first disappears from the forest, we find that the difference between accumulation rates in the open and forest is larger than the difference between ablation rates. Thus, canopy snow interception and subsequent loss, rather than ablation, set up longer snow duration in the open. However, at two relatively windy sites (hourly average wind speeds up to 8 and 17 m/s), differential snow disappearance timing is reversed: Snow persists 2–5 weeks longer in the forest. At the windiest sites, accumulation rates in the forest and open are similar. Ablation rates are higher in the open, but the difference between ablation rates in the forest and open at these sites is approximately equivalent to the difference at less windy sites. Thus, longer snow retention in the forest at the windiest sites is controlled by depositional differences rather than by reduced ablation rates. These findings suggest that improved quantification of forest effects on snow accumulation processes is needed to accurately predict the effect of forest management or natural disturbance on snow water resources.     

Snow measurement by GPS interferometric reflectometry: an evaluation at Niwot Ridge,Colorado

Snow is a critical storage component in the hydrologic cycle, but current measurement networks are sparse. In addition, the heterogeneity of snow requires surveying larger areas to measure the areal average. We presented snow measurements using GPS interferometric reflectometry (GPS‐IR). GPS‐IR measures a large area (~100 m²), and existing GPS installations around the world have the potential to expand existing snow measurement networks. GPS‐IR uses a standard, geodetic GPS installation to measure the snow surface via the reflected component of the signal. We reported GPS‐IR snow depth measurements made at Niwot Ridge, Colorado, from October 2009 through June 2010. This site is in a topographic saddle at 3500 m elevation with a peak snow depth of 1.7 m near the GPS antenna. GPS‐IR measurements are compared with biweekly snow surveys, a continuously operating scanning laser system and an airborne light detection and ranging (LIDAR) measurement. The GPS‐IR measurement of peak snowpack (1.36–1.76 m) matches manual measurements (0.95–1.7 m) and the scanning laser (1.16 m). GPS‐IR has RMS error of 13 cm (bias = 10 cm) compared with the laser, although differences between the measurement locations make comparison imprecise. Over the melt season, when the snowpack is more homogenous, the difference between the GPS‐IR and the laser is reduced (RMS = 9 cm, bias = 6 cm). In other locations, the GPS and the LIDAR agree on which areas have more or less snow, but the GPS estimates more snow on the ground on tracks to the west (1.58 m) than the LIDAR (1.14 m). Copyright © 2011 John Wiley & Sons, Ltd.     

Validation of spatial variability of snowpack thickness and density obtained with GPR and TDR methods

We evaluate the reliability of the joint use of Ground Penetrating Radar (GPR) and Time Domain Reflectometry (TDR) to map dry snow depth, layering, and density where the snowpack thickness is highly irregular and the use of classical survey methods (i.e., hand probes and snow sampling) is unsustainable.We choose a test site characterised by irregular ground morphology, slope, and intense wind action (about 3000 m a.s.l., Western Alps, northern Italy) in dry snow conditions and with a snow-depth ranging from 0.3 m to 3 m over a few tens of metres over the course of a season.The combined use of TDR and high-frequency GPR (at a nominal frequency of 900 MHz) allows for rapid high-resolution imaging of the snowpack. While the GPR data show the interface between the snowpack and the ground, the snow layering, and the presence of snow crusts, the TDR survey allows the local calibration of wave speed based on GPR measurements and the estimation of layer densities. From January to April, there was a slight increase in the average wave speed from 0.22 to 0.24 m/ns from the accumulation zone to the eroded zone. The values are consistent with density values in the range of 350–450 kg/m³, with peaks of 600 kg/m³, as gravimetrically measured from samples from snow pits at different times. The conversion of the electromagnetic wave speed into density agrees with the core samples, with an estimated uncertainty of about 10%.