Introduction
Avalanche- and flood-warning authorities need precise snow-property data. Hydropower stations in alpine and Nordic countries depend on the snow parameters of their catchments to predict the filling stage of their water-storage reservoirs. The distribution of liquid-water flow paths in a snowpack is also important to know for solute transport and wet-snow avalanche formation.
Radar remote sensing has the potential to provide this information if accurate ground reference measurements are available for calibration. For this purpose, we developed and tested a new in situ sensor system for snow moisture and density determination in a snowpack. The system consists of up to 25 m long radio-frequency transmission lines cast in a flat PVC band. The dielectric coefficient of the snow along the flat-band cable is measured with a time-domain reflect-ometer (TDR) at high frequencies and with a low-frequency impedance analyzer. Ice shows significant differences in dielectric properties at high and low frequencies, which enables the determination of the volumetric ice, air and water content, as well as of the snowpack density using measurements at different frequencies. The cable may either be mounted sloping from a mast to the ground to measure vertical snow-profile characteristics or be placed on the snow surface prior to a snowfall to measure the horizontal variations of the dielectric coefficient in the adjacent snow layer.
The objective of this paper is to present the basic principles of the sensor and to show its performance at a high-alpine test site in Switzerland during two winter seasons.
Materials and Methods
The snow cable sensor
The sensor consists of a three-wire flat-band cable with a plastic insulation as shown in Figure 1. The total capacitance C of the equivalent circuit is a non-linear function of the dielectric properties εm of the surrounding material (Reference HuebnerHuebner and others, 1997):
The parameters C1, C2 and C3 (Fig. 2) are measured (travel time and wave impedance) and calculated with electromagnetic field simulations. The values for symmetric excitement are compared inTable 1 together with the transmission line inductance.
There is good agreement between the measurement and the calculation. The deviation for C2 is due to the manufacturing process with its tolerance of ±0.1mm of insulation thickness.
Air-gap correction
One of the challenging problems for continuous measurements in snow is the disturbing interaction of the sensor with the surrounding material. An air gap may form around the sensor due to mechanical influences (induced by wind) or partial melting. In order to detect and correct the effect of an air gap, the following procedure is applied. The three-wire sensor cable is measured in two modes, with small and with large spacing as shown in Figure 3.
In the large spacing mode, the air gap has less influence on the measurement of the dielectric coefficient than in the small spacing mode. This influence can be calculated and a correction algorithm can be derived (Reference Brandelik, Huebner, Doepke and WunderleBrandelik and others, 1998).
In addition to the air-gap correction algorithm, a snow profile reconstruction code has been developed (Reference Becker, Brandelik, Hübner, Schädel, Scheuermann and SchlägerBecker and others, 2002) which allows the spatial determination of the liquid-water content or the snow-layer density (for situations without liquid water) along the cable sensor.
Field test
The performance of the sensor system was tested during two winter seasons (2001/02 and 2002/03) at the experimental site Weissfluhjoch, Davos, Switzerland, located at an altitude of 2550 m (Fig. 4). The site is snow-covered from late October to June, and the winter climate is harsh, with minimum air temperatures of down to –25°C and high wind speeds. In both winters, four flat-band cables were set up: one of them sloping from a mast to the soil surface to measure vertical properties of the snowpack (snow depth, density profile, snow water equivalent (SWE)), and three of them horizontally at different depths of the snowpack to capture the spatial variability of liquid-water content and snow density. These horizontal cables were placed on the snow surface at different dates of the winter (Table 2) and followed the natural settling of the snowpack. In addition, in the second winter a dummy sloping cable was mounted at Weissfluh-joch, as well as at a sub-alpine site in central Switzerland (Alptal (1200m)), to examine the influence of the cable on the snowpack during the snowmelt. The electronic measurement devices of the system (TDR cable tester Tektronix 1502B; impedance analyzer HP-VNWA 8712; personal computer (PC)) were placed in a shelter and connected with 50 m (first winter) and 35 m (second winter) coaxial cables to the flat-band sensor.
To validate the snow-cable measurements, snow profiles were excavated every 15 days close to the sensor, where the snow density was determined layer-wise. Also, an ultrasonic snow-depth sensor recorded the snow depth continuously, and meltwater outflow from the snowpack was measured in a 2 62 mlysimeter next to the sloping-cable sensor. Standard meteorological variables were monitored continuously throughout both winter seasons.
Results
Overall behaviour of the new snow-sensing system
During the test at Weissfluhjoch we closely inspected the behaviour of the snow-band cable and the electronic system under harsh alpine winter conditions. In particular, the following critical components of the system were examined:
Suspension and set-up of the sloping cable
The original suspension, consisting of a wooden fastener attached with ropes to two iron springs, failed twice during the first winter: after a heavy wind storm in the second half of January, and towards the end of the winter season, when the snow load (maximum SWE: 740mm) became too heavy and a spring broke. Consequently, an improved suspension was designed for the second winter, consisting of a cylinder with a mechanical brake (Fig. 5), which kept a constant pressure on the cable. The suspension adjusted continuously when the snow load increased, by releasing the cable.
The slope angle of the cable sensor is one of the key factors for the formation of an air gap around the flat-band cable, as well as for the bending of the cable inside the snow-pack. In the first winter we set up the cable at an angle of 30° to the ground surface and tried to keep it in its original position by fixing it on three vertical poles. Such poles, however, proved to have negative impacts during the snowmelt by preferentially melting the snow and generating air gaps around the pole. Also a sharp bend of the cable formed at the head of the mid-pole after snowpack settling during the melting period (Fig. 6). In addition, the poles were not stable enough to withstand the heavy snow load and bent considerably towards the end of the season. Consequently, for the second winter we decided to set up the cable at a steeper angle (45°) and to omit the supporting poles. The steeper installation reduced the snow load on the cable, and the cable was kept in its original position during the whole winter. However, we noted increased flickering of the cable in moderate to high wind speeds in the absence of supporting poles. These movements of the cable most likely generated the large air gaps observed at the surface, and may have caused a cable failure during a snowstorm in January 2003. A further complication was a 90° twist of the cable around its own axes, which may have increased the water flow along the cable (Fig. 6).
Apart from these problems, we believe that the new suspension at a steeper angle has improved the installation, and could be further improved by the combination with springs to reduce the lateral movements of the cable.
Air gap
Air gaps along the flat-band cable are highly undesirable since they not only falsify the electronic signal, but may also transfer melting water preferentially along the cable. The influence of an air gap on the measurement results can be partially corrected with a correction method (Reference HuebnerHuebner and others, 1997), but with respect to a potential preferential water transport along the cable it was crucial to thoroughly inspect the size of the air gap and evaluate its importance as a water-flow path.
At both sites the air gap showed up in different sizes and shapes, depending on wind speed, wind direction and melting intensity (Fig. 7). The assumption of a symmetric air gap, which is made in the correction algorithm, was not appropriate for many situations. There were dates (27 February 2003 at Alptal) when an up to 5 cm opening formed on one side of the band cable, whereas the other side was in close contact with the snow. At the Weissfluhjoch site the air gap was mostly generated by wind. Hence, the compass orientation of the flat-band cable was decisive for the shape of the air gap. We also observed that air gaps were able to fill up again when the wind conditions had changed. Our overall conclusion was that the shape of the air gap at the snow surface was very dynamic and unpredictable. Inside the snowpack, we found no evidence for significant voids around the cable when we excavated the flat-band cable at the sub-alpine site, Alptal, on 1 March 2003. This was at a stage when the snow-pack had settled to a density of 290 kgm–3 and liquid melt-water had penetrated to the bottom of the snowpack (Fig. 8).
Close inspection of the flat-band cable was also carried out to examine whether the snow band had conducted melt-water preferentially during the snowmelt. A dye tracer had been applied 15 days earlier (prior to the snowmelt) on the snow surface (area indicated in Fig. 8). Avertical profile was excavated along the cable and perpendicularly at the lower end of the cable. No stained water was detected at the lower end of the cable, whereas to the left and right of the cable end considerable stained areas were detected originating from lateral water transport along snow layers. There is evidence that some dye did follow the cable for some distance before finding alternative flow paths to the side of and below the cable, but we did not find indications of extensive lateral water flow expressed in terms of intensely stained areas following the direction of the cable. We think that we can conclude in general that water flow along the cable is not greater than the natural lateral water flow in the snowpack.
Concerning thermal effects of the sensor on the surrounding snow, we found some thermal radiation effects of a horizontal cable only once, shortly after it was laid out on the snow surface and when the coverage by fallen snow was very shallow. But as soon as the cable was covered with 410 cm of snow, this effect was no longer found.
The electronic system
Except for some weeks in December 2001, when an erroneous setting prevented the impedance analyzer from running, and in April 2002, when the breakdown of an internal electronic component forced us to replace the impedance analyzer, the electronic system ran more or less failure-free throughout the first winter season. The system uptime was also high in the second winter (80%). The automatic start of the measurements occasionally failed at the interface between the controlling PC unit and the Tektronix. Since this occurred mainly at very low temperatures, we think this is a temperature problem of the electronic components of the Tektronix device, which was not built to operate under such low temperatures.
Comparison of snow sensor measurements with manual measurements
Dielectric constant, bulk snow density and liquid-water content of snowpack (mean value along the cable) were calculated from the raw signals. And the moisture distribution along one of the horizontal cables was reconstructed for different stages of the winter. The natural settling of the snow cover was well reflected in the horizontal cable measurements (Fig. 9). The snow density increased from initially 100–200 kgm–3 to approximately 500 kgm–3 at the end of winter, which was in accordance with manual snow-density measurements of the corresponding snow layers. The two cables located at two different depths of the snowpack also reproduced correctly a higher density in the lower part of the profile in February (cable 1) and a faster compaction in the upper part during the snowmelt (cable 2). However, cable 2 seems to have systematically underestimated density, which we think may be explained by spatial differences in the measurement spots of cable and manual density measurements.
No liquid water was detected with the horizontal cables until the end of April; after that the snowpack had reached its maximum depth (Fig. 10). This is supported by the automated snow temperature measurements, which indicate dry snow conditions and temperatures well below zero before this period. Once the snowmelt set in, we measured the steadily increasing liquid-water content, with the horizontal cables indicating the downward penetration of the wetting front. The upper cable (No. 2) not only reflected an earlier and faster increase in liquid-water content than the lower cable, it also reacted more clearly to the diurnal variation caused by melting during the day and refreezing during the night.
With regard to the spatial variation of liquid water, the horizontal flat-band cable nicely demonstrated the formation of preferential water-flow paths adjacent to the cable (Fig. 11), which is the natural water-transport process in a melting snowpack (Reference Marsh and WooMarsh and Woo, 1984). The horizontal cable at about 1 m below the snow surface indicated in mid-May an emerging water-conducting zone at 14 m from the cable end, as well as a newly developing flow finger at 8 m from the cable end.
Discussion and Conclusions
The measurement results yielded good correspondence of the snowpack density with manual reference measurements taken twice a month at the measurement site. Also, the determination of liquid-water content from measurements of the dielectric constants according to Reference Tiuri, Sihvola, Nyfors and HallikainenTiuri and others (1984) gave plausible results, both compared to lysimeter data taken on the test field and with regard to the spatial variation of flow fingers that we normally experience in a natural snowpack. The device’s suitability for these measurements has been demonstrated, but the instrumentation set-up needs to be improved to withstand the harsh conditions of alpine winter seasons. This will be a task for future improvements.
We learned that the suspension of the cable sensors is a crucial point in the measurement. For successful measurements it is indispensable to reduce the following sources of error:
Air gaps between cable sensor and surrounding snow must be minimized, although it is possible to correct for minor effects.
The suspension of the cable sensors must allow the cable to settle with the snowpack, but at the same time must be robust enough to withstand heavy winds and snowl oads.
Acknowledgements
The presented work was carried out within the European Union project SNOWPOWER (fifth framework, NNE5/2000/251).