Hostname: page-component-78c5997874-8bhkd Total loading time: 0 Render date: 2024-11-05T10:19:01.579Z Has data issue: false hasContentIssue false

Optical properties of snow in backscatter

Published online by Cambridge University Press:  08 September 2017

S. Kaasalainen
Affiliation:
Finnish Geodetic Institute, Geodeetinrinne 2, PO Box 15, FIN-02431 Masala, Finland E-mail: [email protected] Finnish Meteorological Institute/Arctic Research Center, Tähteläntie 62, FIN-99600 Sodankylä, Finland
M. Kaasalainen
Affiliation:
Department of Mathematics and Statistics, University of Helsinki, PO Box 68, FIN-00014 Helsinki, Finland
T. Mielonen
Affiliation:
Finnish Geodetic Institute, Geodeetinrinne 2, PO Box 15, FIN-02431 Masala, Finland E-mail: [email protected]
J. Suomalainen
Affiliation:
Finnish Geodetic Institute, Geodeetinrinne 2, PO Box 15, FIN-02431 Masala, Finland E-mail: [email protected]
J.I. Peltoniemi
Affiliation:
Finnish Geodetic Institute, Geodeetinrinne 2, PO Box 15, FIN-02431 Masala, Finland E-mail: [email protected]
J. Näränen
Affiliation:
Finnish Geodetic Institute, Geodeetinrinne 2, PO Box 15, FIN-02431 Masala, Finland E-mail: [email protected] Department of Astronomy, University of Helsinki, PO Box 14, FIN-00014 Helsinki, Finland
Rights & Permissions [Opens in a new window]

Abstract

We present an overview and the first systematic results on the backscatter properties of snow samples with different grain properties. During a 2 year study we investigated the effect of the apparent characteristics of the samples on the intensity enhancement in the direction of exact backscatter (retroreflection), using a specific instrument designed for backscatter measurements. We observed that a sharp peak in intensity (the hot spot) occurs for most types of snow, which is not observable with traditional goniometers. The key factors in the peak properties are the temperature (related to the changes in grain structure) and grain shape and size. These results form the basis of a larger backscatter data archive which is being applied in the systematic study and exploitation of ‘hot spots’ in remote sensing of natural targets, as well as in the development of airborne laser intensity measurement.

Type
Research Article
Copyright
Copyright © International Glaciological Society 2006

1. Introduction

The optical properties of snow play a major role in remote sensing of the cryosphere and mapping of icy regions. Methods for deriving the key factors of snow- and ice-cover variations, such as albedo, metamorphism and grain characteristics (Reference Grenfell, Warren and MullenGrenfell and others, 1994; Reference Fily, Dedieu and DurandFily and others, 1997, 1998, 1999), are developed for more effective detection of changes in snow cover. Mapping of snow- and ice-covered areas provides important indications of climate change, hydrological processes and water resources (Reference Fily, Leroux, Lenoble, Sergent, Schmitt, Bergh and FestouFily and others, 1998; Reference Magagi and BernierMagagi and Bernier, 2003).

Remote sensing of snow and other land surfaces is based on spectral and directional features of light scattered from the surfaces (see Reference Fily, Leroux, Lenoble, Sergent, Schmitt, Bergh and FestouFily and others, 1998; Reference Nolin and LiangNolin and Liang, 2000). The bidirectional reflectance distribution function (BRDF) describes the directional properties of light reflected from a surface, and it is used both in directional corrections of, and as a source of information in, remotesensing data (Reference Verstraete and PintyVerstraete and Pinty, 2001, and references therein). Comparisons with models and in situ measurements are important for more reliable information to be derived from remote instruments. Computational studies of the directional reflectance of snow are found in Reference Leroux, Lenoble, Brogniez, Hovenier and De HaanLeroux and others (1998) and Reference Mishchenko, Dlugach, Yanovitskij and ZakharovaMishchenko and others (1999). Among the most recent goniometric experiments of the bidirectional reflectance of snow are those by Reference Aoki, Aoki, Fukabori, Hachikubo, Tachibana and NishioAoki and others (2000), Reference Painter and DozierPainter and Dozier (2004) and Reference Peltoniemi, Kaasalainen, Näränen, Matikainen and PiironenPeltoniemi and others (2005).

The increased reflectance at small phase angles (light-source-target-detector angle; see Fig. 1) is observed as a ‘hot spot’ in remote images for targets such as forest and vegetation (Reference Hapke, DiMucci, Nelson and SmytheHapke and others, 1996; Reference Comacho-de-Coca, Bréon, Leroy and Garcia-HaroComacho-de-Coca and others, 2004). The increase is particularly pronounced at zero phase and its immediate vicinity, which is often referred to as ‘exact backscattering’ to distinguish it from backward scattering in the same hemisphere in which the light source is placed (also termed retrosolar scattering in remote sensing). Backscattering has recently gained more importance in terrestrial remote sensing, since there are new instruments that observe at, or close to, exact backscatter (e.g. POLDER and HyMap; see Reference Comacho-de-Coca, Bréon, Leroy and Garcia-HaroComacho-de-Coca and others, 2004). The most important observational property of exact backscatter is that it cannot be reached with goniometers, as in these the detector necessarily blocks the light source. However, the effect at zero-angle vicinity undoubtedly contains a significant amount of information on the target material: as Van Reference Van Albada and LagendijkAlbada and Lagendijk (1985) note, the effect is simply much too sharp to be observed without a beam-splitter.

Fig. 1. The small-angle laser goniometer. Retroreflection from the sample is observed through the beam-splitter by the chargecoupled device (CCD) camera (left). Non-zero phase angles (source-target-detector) are obtained by rotating the beam-splitter (right). The distance from beam-splitter to sample was 12cm, and from the CCD to the sample 120 cm.

There are few ground reference data available on natural hot spots, even though the backscatter effect has been measured systematically for various substances in laboratory experiments. Laboratory experiments on backscattering patterns have been carried out mostly for powders and particles, and applied in the fields of optics, condensed matter physics, astronomy and medicine (for more references see Reference Yoon, Ghosh Roy and StraightYoon and others, 1993; Reference ShkuratovShkuratov and others, 2002; Reference Kaasalainen and RautiainenKaasalainen and others, 2005a). Small-scale laboratory measurements of natural targets have been made, for example for clumps of grass or needles (Reference Hapke, DiMucci, Nelson and SmytheHapke and others, 1996). However, a large number of measurements, particularly field measurements made more or less in situ, would be required to understand and exploit better the hot spot signatures and establish a remote-sensing ground reference databank.

Snow (back)scattering, exact or general, is still a relatively unexplored field. It has been studied mainly for snow crystals in the atmosphere (e.g. Reference Borovoi, Grishin, Naats and OppelBorovoi and others, 2000) and mapping of snow cover at zero phase angle by radar or microwave, with connections to physical parameters (Reference Fily, Leroux, Lenoble, Sergent, Schmitt, Bergh and FestouFily and others, 1998; Reference Strozzi and MätzlerStrozzi and Mätzler, 1998; Reference Magagi and BernierMagagi and Bernier, 2003). This is essentially the objective in the optical regime as well. In astronomy, strong backscatter peaks (called opposition effect in this context) have been observed for numerous ice-covered surfaces (including those with water ice), such as rings and satellites (Reference Domingue and VerbiscerDomingue and Verbiscer, 1997; Reference Poulet, Cuzzi, French and DonesPoulet and others, 2002). Optical backscattering of snow on the ground has been addressed in only a few studies thus far (Reference Piironen, Muinonen, Keränen, Verstraete, Meneti and PeltoniemiPiironen and others, 2000; Reference ShkuratovShkuratov and others, 2002; Reference Kaasalainen, Karttunen, Piironen, Virtanen, Näränen and LiljeströmKaasalainen and others, 2003), with a few peaks observed for some samples. These studies have focused mainly on the effects of impurities. Systematic backscattering experiments are lacking, and no description of snow backscatter, or the key parameters related to it, is available. This lack of information is one of the main incentives for undertaking the present study. Another key factor is that this study makes use of a standard experimental set-up that has been used widely in backscatter studies of materials other than snow. Therefore its use here is a natural extension of these applications.

On the theoretical side, there are a number of studies on exact backscattering in both optics and light scattering, so the need to access the zero phase-angle vicinity properly is also important for developing the theoretical methods. Indeed, much of current backscatter theory has remained tentative due to a lack of the systematic experiments needed for building any hypotheses into a theory. In fact, current theoretical results and recent experiments have several mutual inconsistencies or uncertainties (Reference Nelson, Smythe, Hapke and HaleNelson and others, 2002), so there is an obvious need for new measurements. Although the backscatter effects are extremely strong, even their basic physical interpretation requires both experimental and theoretical studies. The current physical interpretation of hot spots is based on shadow hiding in rough surfaces and coherent backscatter, which is a microscopic-scale interference effect (Reference Hapke, DiMucci, Nelson and SmytheHapke and others, 1996; Reference ShkuratovShkuratov and others, 2002).

This paper presents the results of a 2 year campaign of snow backscattering measurements, the first extensive study with systematic variation of snow parameters. The existing, well-established techniques for measurement of backscattering were applied to snow surfaces to obtain information on the backscattering properties of snow and to search for the main contributors to snow hot spot. The potential applications to glaciology and snow studies via remote sensing are obvious as airborne laser intensity measurements are under development (Reference Geist, Lutz and StötterGeist and others, 2003; Reference Kaasalainen and RautiainenKaasalainen and others, 2005b). The laser altimetry measurement deals with the retroreflected laser beam, but the intensity information has not been discussed thus far. The search for a tool for change detection is of primary importance to glaciology in this context. It is essential to realize that this approach is quite feasible even when exact quantitative inference of the physical properties of the target is uncertain due to insufficient theoretical modelling, noisy data, or even fundamental ambiguities in the inverse problem: change detection is possible without any parametric modelling of the physical properties of the target.

The paper is organized as follows: in section 2 we present the experimental and sampling methods; in section 3 we discuss the results, with particular emphasis on melting and refreezing, specular reflections and grain properties; and in section 4 we present our conclusions and discuss plans for future work.

2. Methods

2.1. The experiment

From the first experimental demonstration of coherent backscatter (Van Reference Van Albada and LagendijkAlbada and Lagendijk, 1985), experimental studies of backscatter have been developed in many fields, and the techniques based on the use of a beamsplitter have become well established, as well as the use of charge-coupled device (CCD) detectors in the measurement of (backscattered) laser intensity (e.g. Reference Yoon, Ghosh Roy and StraightYoon and others, 1993; Reference Wiersma, van Albada, B.A. and LagendijkWiersma and others, 1995; Reference Nelson, Hapke, Smythe and SpilkerNelson and others, 2000; Reference ShkuratovShkuratov and others, 2002). The measurement presented here is based on these techniques, and has been used widely in laser-based backscatter measurement in optics and related fields for laser intensity measurement and standard CCD photometry used in the laboratory and, for example, in astronomical studies.

The measurements were made using the small-angle goniometer constructed particularly for backscattering measurements, as shown in Figure 1. A 632.8 nm He-Ne laser was used as a light source. The backscatter geometry is worked out using a beam-splitter. The sample was placed on a rotator to smooth out the laser speckle and obtain a larger sample area as the spot (about 4mm in diameter) moved on the surface. Transmittance filters of 10% and 3% were used in front of the light source to keep the signal in the linear range in the CCD frame. More detailed descriptions of the instrument are provided by Reference Näränen, Kaasalainen, Peltoniemi, Heikkilä, Granvik and SaarinenNäränen and others (2004) and Reference Kaasalainen and RautiainenKaasalainen and others (2005a, b). We modified the original laboratory set-up to be robust enough for field conditions.

Two separate measurement campaigns were carried out: the first in February/March 2004, producing a set of preliminary results and mainly aimed at instrument development; and the second in February/March 2005 with an improved instrument and measurement facilities. Both campaigns took place at the Sodankyla Observatory (67.368° N, 26.633° E, 179ma.s.l.) of the Arctic Research Center (Finnish Meteorological Institute). The site is located in a riverside area covered with coniferous forest and surrounded by typical Lappish marshland. The 2004 measurements were made outside or in a covered entrance room with no heating, and the 2005 measurements were made inside a small hut with no heating. The hut had no windows, which allowed measurements to be made during the daytime (a dark workspace was needed to minimize the background noise). The space was occasionally heated (e.g. up to –5°C) during periods as low as –30°C, which prohibited optical measurements (because of precipitation of humidity into the surfaces of the optics).

Typically 25 data points were measured for each phase curve (averaged from five images at the same angle). Fast data acquisition was important to minimize the effects of grain sintering and metamorphism on the sample. On the other hand, camera exposures of sufficient duration were needed to allow a large fraction of the (rotating) sample to be exposed. The range in phase angle varied from –3° to +10–15°; the positive and negative signs denote the viewing geometry around zero phase, from where the spot moved towards and away from the light source, respectively (see Fig. 1). The nominal angular resolution (determined from the source and detector apertures) of the instrument varied between 1.3° (2004) and 0.4° (2005). This causes some rounding of the peaks near zero, especially in the 2004 data, but the differences in the peak properties are well observable because of the extreme sharpness of the effect itself (for more details of the effects of resolution see Reference Kaasalainen and RautiainenKaasalainen and others, 2005a).

2.2. Sampling

To achieve systematic sampling, we made several measurements over a long timescale; there is little opportunity for a priori sample variation with natural snow targets, especially in field conditions. The device has to operate on a solid foundation to achieve the proper optical alignment (cf. Reference Kaasalainen and RautiainenKaasalainen and others, 2005a), so the samples from the snowpack were transferred into a small plastic container, which was placed underneath the beam-splitter. This results in small sample sizes (the sample thicknesses varied between 3 and 5 cm; see Tables 1 and 2), but this does not affect the comparable properties of the backscatter peak shape, since the directional signal from the platform adds to the diffuse (rather than the directional) reflection as it is scattered on the sample layer above. The sample surfaces were smoothed with a spade (essentially just wiping out the excess snow), except for hard (refrozen) layers which were extracted as a chunk (see the layers photographed in Fig. 2). The approximate sample densities were measured by weighing the sample in its container (of known volume), and the grain properties were analyzed visually using a 0.1mm scaled field magnifier loupe. In the 2005 campaign, we also took images of the grains through the loupe (see Figs 310).

Table 1. Data from 2004. Snow sample characterization (approximate) and backscattering properties from the fits of the exponential function (see section 3.1)

Table 2. Data from 2005. Snow sample characterization (approximate) and backscattering properties from the fits of the exponential function (see section 3.1). R = sample taken from the snowpack on top of the river ice; T = surface tilted. The approximate sample thickness was 4 cm, except for the melted/refrozen samples for which it was 3cm for samples taken between 1 and 3 March and 3.5cm for both samples taken on 13 March

Fig. 2. Refrozen layers from the 2004 campaign measured (a) below (11 March 2004) and (b) above 0°C (19 March 2004).

Fig. 3. (a) Melting and (b) frozen surface grains measured on 2 March 2005. The grain images were taken with a millimeter scale (more clearly visible in the bottom left image in Fig. 4).

Fig. 4. (a) Melting grains (same sample from 2 March 2005 as shown in Fig. 3a) and (b) same sample after refreezing. The millimeter scale visible in the bottom left image is the same in all the grain images.

Fig. 5. (a) Melting new snow produces a stronger backscatter effect than does (b) the refrozen sample (both from 3 March 2005). Both surfaces were tilted to elimintate the effect of specular reflection.

Fig. 6. (a) Refrozen (rounded) grains (same sample as shown in Fig. 5b) and (b) grains of a coarser structure taken from 40cm depth in the snowpack (also measured on 3 March 2005). The difference in the smoothness of the phase curve is visible (cf. also Fig. 9).

Fig. 7. Effect of compression on surface (needle) grains: (a) slightly compressed (7 March 2005) and (b) compressed overnight (8 March 2005).

Fig. 8. Comparison of the effect of grain size: (a) 0.7–1.0mm rounded grains from 5–10cm depth in the snowpack (12 March 2005) and (b) same sample slightly melted and refrozen (13 March 2005).

Fig. 9. Effect of grain-size: (a) 0.2–0.5mm surface grains from previous day’s snowfall (8 March 2005) and (b) 0.1–0.3mm surface grains from previous day’s snowfall (8 March 2005; same sample as shown in Fig. 7b).

Fig. 10. Larger grain sizes (than those shown in Figs 8 and 9): (a) depth hoar with grains of 2mm from about 80cm depth in the snowpack (4 March 2005) and (b) a similar sample with grains of 2–3mm (10 March 2005).

The snow properties were classified according to Reference ColbeckColbeck and others (1990). The characterization presented here is sufficient to enable us to highlight the differences in backscattering peaks of different types, and we concentrated on the characteristics that were shown to have an effect on backscatter in previous experiments for other samples (Reference Nelson, Hapke, Smythe and SpilkerNelson and others, 2000; Reference ShkuratovShkuratov and others, 2002; Reference Näränen, Kaasalainen, Peltoniemi, Heikkilä, Granvik and SaarinenNäränen and others, 2004; Reference Kaasalainen and RautiainenKaasalainen and others, 2005a).

3. Results and Discussion

We employed the standard practice of relative CCD photometry, normalizing all intensities to the Spectralon spectroscopic standard (Labsphere Inc.). For this we used its phase curve average, excluding the smallest phase angles where its (linearly polarized) reflectance is not constant. To compare more easily the intensity levels at 0° and 5° (i.e. outside the backscatter peak) phase angles, a fourparameter empirical function was fitted to most data to achieve a better approximation of the intensity level at 0° and 5° and compute the peak amplitude, I(0°)/I(5°), and the half-width-half-maximum (HWHM, which is a measure of the apparent width of the backscatter peak) phase angle (see, e.g., Reference Nelson, Hapke, Smythe and SpilkerNelson and others, 2000; Reference Kaasalainen and RautiainenKaasalainen and others, 2005b). The exponential function gives the relative intensity, I(a), as a function of phase angle: I(a) = a exp(—α/d) + b + kα, where a, d, k and b are empirical parameters (see also Reference Piironen, Muinonen, Keränen, Verstraete, Meneti and PeltoniemiPiironen and others, 2000; Reference Kaasalainen and RautiainenKaasalainen and others, 2005a). The intensities, I(0°), peak amplitudes, I(0°)/I(5°), and HWHMs are tabulated with the sample characterizations in Tables 1 and 2.

Large deviations in the data are typical for snow (cf. earlier measurements in backscatter and other direction, (e.g. Reference Kaasalainen, Karttunen, Piironen, Virtanen, Näränen and LiljeströmKaasalainen and others, 2003; Reference Peltoniemi, Kaasalainen, Näränen, Matikainen and PiironenPeltoniemi and others, 2005)). 2005)). Another feature in snow measurements is that the calibrated intensity levels are lower than reported snow albedo values at 633 nm (see Reference Grenfell, Warren and MullenGrenfell and others, 1994). This results from the laser beam entering deep into the snow layer due to strong multiple scattering inside the layer, which, unlike for other types of samples or standards, is not observable with all photometric techniques. Since the diffuse radiation only adds to the background intensity, and the strongest directional signal comes from the surface, the comparison of backscatter peak characteristics is still possible in a relative sense. These results must be considered preliminary, but they give a qualitative overview of the most typical effects present in the backscatter direction. Further investigation of laser-based snow measurement techniques remains an important objective of future study. It is also apparent that, for snow data, many features are visible from a large bulk of repeated measurements rather than a single curve. Representative datasets are plotted in Figures 212, and photographs of the grains taken with the aid of a millimeter-scaled viewing loupe are presented. No small-scale photographs are available for the 2004 samples, but large-scale images are presented to give an overall view of the sample structure.

Fig. 11. Comparison with data from 2004: (a) aged snow from about 40 cm depth in the snowpack and (b) depth hoar from 60 cm depth in the snowpack. The samples (both measured on 23 March 2004) are photographed in a sample cup of about 11cm diameter.

Fig. 12. Data from 2004 measurements: (a) new snow from the surface and (b) aged snow from about 20 cm depth in the snowpack (both samples measured on 23 March 2004). The effect of rounded grain structure (compared with fresh fallen grains of a finer structure) on the smoothness of the phase curve is also clear in this case. The sample grains are photographed in a sample cup of about 11cm diameter.

3.1. Melting and refreezing

As seen from Figures 212, the graphs are smoothest and the results more accurate for melted and refrozen samples. These samples also produced the strongest and most distinguishable peaks, whereas much greater deviation in data is observed for the samples measured below 0°C, which may also have masked the peak. This is observed most clearly in Figure 3 for the two similar samples, the other of which (Fig. 3a) was measured after heating the hut up to 3°C. There is a notable difference in the relative peak amplitude (Table 2). A strong peak was still observed a few hours later after refreezing the same sample (Fig. 4), which suggests that the role of grain structure (which changes greatly in the process of melting) is stronger than that of temperature. This is also supported by comparison with the 2004 data: there is no substantial difference in the peak amplitudes of a refrozen surface measured below and above the freezing point (Fig. 2). In earlier studies of the effect of liquid water on snow reflectance (see Reference WarrenWarren, 1982, and references therein), the spectral reflectance of refrozen snow was found to be very similar to that of melting snow. Temperature seems to have an indirect effect by affecting the grain structure and hence the optical properties.

3.2. Specular reflection from the surface

The effect of specular reflection from a wet surface can be investigated by comparing the measurement in Figure 4a with that in Figure 5a, and Figure 4b with Figure 5b. The samples are otherwise similar, but those in Figure 5 have been tilted (by tilting the other edge of the sample cup about 1cm higher) to avoid the normal incidence into the sample surface, and hence the possible specular reflection back to the camera at 0° phase angle. As strong backscatter peaks are still observed (see also Table 2), it is apparent that they are not caused by specular reflection from the surface. More studies on the role of these reflections on bidirectional reflectance of snow would still be worthwhile.

3.3. Grain proprties

The mean grain size of the samples was approximated as the diameter of the longest dimension of the grains. As discussed above, factors other than grain size seem to have a stronger effect on the backscatter peak. However, comparing the samples measured in freezing temperatures, the curve shapes appear sharper for those with larger grains, whereas the smallest peaks occur for those with a grain size of < 0.5mm (Table 2). The samples with large grains mostly have smaller reflectances, which emphasizes the relative increase in brightness towards 0°. This is in agreement with the 2004 data, for which the samples with grain sizes > 1mm generally produced the largest peaks (Table 1). Also, depth hoars have large enhancements from 5° towards 0°, but there is deviation in the data which is likely to be caused by the grain structure (see the images of the typical depthhoar grains in Fig. 10).

The effect of grain rounding is presented in Figures 6 and 9. The effect is also seen for 2004 data in Figure 12, where the phase curve of the rounded grains from 20 cm depth in the snowpack is also smoother than that of the finer-structured column/hexagon-type grains, and the errors are smaller. Overall, the winter of 2004 was warmer, with several periods of temperatures above zero, which seems to have had an effect on all samples: all the grain structures were observed to be smoother than those observed in 2005, which were mostly coarser in structure. Therefore the results agree with observations for melted and refrozen samples (Figs 3 and 4): the smoother structure causes a stronger hot spot.

We also attempted to determine the effects of packing density by mechanically compressing some samples or allowing the sample to flatten out overnight (Fig. 7). The only observed effects were an increase in surface brightness and greater deviation in data. Thus factors other than packing density had a stronger effect on the backscatter peak properties. We found no strong links between surface brightness and the amplitude or HWHM of the peak (Table 2). The reflectance also increased for the samples taken from the dense surface layer on top of the river ice, but here too the other properties (such as grain size or the decomposed and faceted structure) played a stronger role.

We found no correlation of the HWHM values with any observable parameter in the samples, except for the 2004 data, for which the HWHM appears to increase with increasing grain size; however, more measurements are required.

4. Conclusion

This study confirms the earlier assumption that snow in general exhibits a backscatter peak (Reference Piironen, Muinonen, Keränen, Verstraete, Meneti and PeltoniemiPiironen and others, 2000; Reference Kaasalainen, Karttunen, Piironen, Virtanen, Näränen and LiljeströmKaasalainen and others, 2003), i.e. most snow types present a sharp intensity surge at phase angles approaching zero. This is also the first and thus far most extensive overview of backscatter (hot spot) effects from different snow types. The key factors in the existence and properties of the peak are temperature (the effect of which is indirect and related to grain rounding) and the grain structure and size, as the (refrozen) samples with the largest grain/agglomerate size produced the most prominent peaks. The results observed here are preliminary, but they represent the first systematic approach to the snow parameter variation in the investigation of snow backscatter.

As far as the physical interpretation of the backscatter peak is concerned, it seems that this effect is dominated by coherent backscatter since the peaks are narrow and other experiments indicate that shadowing starts to take effect at larger angles. This is one explanation for the fact that snow does not show backward scattering in large-angle-scale bidirectional reflectance measurements, whereas rough substances show substantial shadowing effects (cf. lichens in Reference Kaasalainen and RautiainenKaasalainen and Rautiainen 2005).

There is evidence from other fields that backscattering peaks can be applied in target characterization, but the development of methods to be used in remote sensing of snow calls for a large data archive, for which these measurements are a starting point. More accurate laser measurement techniques for snow are also a subject of further investigation. Because of the combined effects of snow properties, backscatter measurements alone may not be suitable for accurate retrieval of surface parameters, but combined with simultaneous data from other sources they could provide an important source of information about the ongoing changes in snow and ice surfaces. Laser intensity data from glacier altimetry may in some cases be a major source of information on a particular surface, which makes the further investigation of laser backscatter intensity worthwhile. Since change detection is possible from repeated (relative) experiments, our results support the application of backscatter measurements in the study and calibration of intensity information from laser scanning, which has recently become a common method in glacier altimetry. Exploiting the intensity information would enable an airborne laser-based change detection method for snow and ice surfaces, which would have a strong impact on climate study.

Acknowledgements

We thank the scientific and technical staff at Arctic Research Center for their support during the experiments, J. Torppa at Helsinki University Observatory for help in the measurements, and J. Piironen and H. Karttunen for comments and ideas during the course of the snow study. This work was financially supported by the The Finnish Society of Sciences and Letters (2004 campaign) and the Academy of Finland.

References

Aoki, T., Aoki, T., Fukabori, M., Hachikubo, A., Tachibana, Y. and Nishio, F.. 2000 Effects of snow physical parameters on spectral albedo and bidirectional reflectance of snow surface. J. Geophys. Res., 105(D8), 10,21910,236. Google Scholar
Borovoi, A., Grishin, I., Naats, E. and Oppel, U.. 2000 Backscattering peak of hexagonal ice columns and plates. Opt. Lett., 25(18), 13881390. Google Scholar
Colbeck, S.C. and 7 others. 1990. The international classification for seasonal snow on the ground. Wallingford, Oxon., International Association of Hydrological Sciences. International Commission on Snow and Ice. Google Scholar
Comacho-de-Coca, F., Bréon, F.M. Leroy, M. and Garcia-Haro, F.J.. 2004 Airborne measurement of hot spot reflectance signatures. Remote Sens. Environ., 90(1), 6375. CrossRefGoogle Scholar
Domingue, D. and Verbiscer, A.. 1997 Re-analysis of the solar phase curves of the icy Galilean satellites. Icarus, 128(1), 4974. Google Scholar
Fily, M., Bourdelles, B., Dedieu, J.P. and Sergent, C.. 1997 Comparison of in situ and Landsat thematic mapper derived snow grain characteristics in the Alps. Remote Sens. Environ., 59(3), 452460. Google Scholar
Fily, M., Leroux, C., Lenoble, J. and Sergent, C.. 1998 Terrestrial snow studies from remote sensing in the Solar spectrum and the thermal infrared. In Schmitt, B., Bergh, C.D. and Festou, M., eds. Solar system ices. Dordrecht, etc., Kluwer Academic Publishers, 421441. CrossRefGoogle Scholar
Fily, M., Dedieu, J.P. and Durand, Y.. 1999 Comparison between the results of a snow metamorphism model and remote sensing derived snow parameters in the Alps. Remote Sens. Environ., 68(3), 254263. Google Scholar
Geist, T., Lutz, E. and Stötter, J.. 2003 Airborne laser scanning technology and its potential for applications in glaciology. In Proceedings of the ISPRS Workshop on 3-D Reconstruction from Airborne Laserscanner and INSAR data. Dresden, International Society for Photogrammetry and Remote Sensing, 101106. Google Scholar
Grenfell, T.C., Warren, S.G. and Mullen, P.C.. 1994 Reflection of solar radiation by the Antarctic snow surface at ultraviolet, visible, and near-infrared wavelengths. J. Geophys. Res., 99(D9), 18,66918,684. Google Scholar
Hapke, B., DiMucci, D., Nelson, R. and Smythe, W.. 1996 The cause of the hot spot in vegetation canopies and soils: shadow-hiding versus coherent backscatter. Remote Sens. Environ., 58(1), 6368. CrossRefGoogle Scholar
Kaasalainen, S. and Rautiainen, M.. 2005 Hot spot reflectance signatures of common boreal lichens. J. Geophys. Res., 110(D20), D20102. (10.1029/2005JD005834.) Google Scholar
Kaasalainen, M., Karttunen, H., Piironen, J., Virtanen, J., Näränen, J. and Liljeström, A.. 2003 Backscattering from snow and ice: laboratory and field measurements. Can. J. Phys., 81(1-2), 135143. Google Scholar
Kaasalainen, S., Peltoniemi, J., Näränen, J., Suomalainen, J.I. Kaasalainen, M. and Stenman, F.. 2005a. Small-angle goniometry for backscattering measurements in the broadband spectrum. Appl. Opt., 44(8), 14851490. Google Scholar
Kaasalainen, S., Ahokas, E., Hyyppa, J. and Suomalainen, J.. 2005b. Study of surface brightness from backscattered laser intensity: calibration of laser data. IEEE Trans. Geosci. Remote Sens., 2(3), 255259. Google Scholar
Leroux, C., Lenoble, J., Brogniez, G., Hovenier, J.W. and De Haan, J.F.. 1998 A model for the bidirectional polarized reflectance of snow. J. Quant. Spectrosc. Radiat. Transfer, 61(3), 273285. Google Scholar
Magagi, R. and Bernier, M.. 2003 Optimal conditions for wet snow detection using RADARSAT SAR data. Remote Sens. Environ., 84(2), 221233. Google Scholar
Mishchenko, M.I., Dlugach, J.M. Yanovitskij, E.G. and Zakharova, N.T.. 1999 Bidirectional reflectance of flat, optically thick particulate layers: an efficient radiative transfer solution and applications to snow and soil surfaces. J. Quant. Spectrosc. Radiat. Transfer, 63(2-6), 409432. Google Scholar
Näränen, J., Kaasalainen, S., Peltoniemi, J., Heikkilä, S., Granvik, M. and Saarinen, V.. 2004 Laboratory photometry of planetary regolith analogs II. Surface roughness and extremes of packing density. Astron. Astrophys., 426, 11031109. CrossRefGoogle Scholar
Nelson, R.M., Hapke, B.W. Smythe, W.D. and Spilker, L.J.. 2000 The opposition effect in simulated planetary regoliths: reflectance and circular polarization ratio change at small phase angle. Icarus, 147(2), 545558. CrossRefGoogle Scholar
Nelson, R.M., Smythe, W.D. Hapke, B.W. and Hale, A.S.. 2002 Low phase angle laboratory studies of the opposition effect: search for wavelength dependence. Planet. Space Sci., 50(9), 849856. Google Scholar
Nolin, A.W. and Liang, S.. 2000 Progress in bidirectional reflectance modeling and applications for surface particulate media: snow and soils. Remote Sens. Rev., 18, 307342. Google Scholar
Painter, T.H. and Dozier, J.. 2004 Measurements of the hemispherical-directional reflectance of snow at fine spectral and angular resolution. J. Geophys. Res., 109(D18), D18115. (10.1029/2003JD004458.) Google Scholar
Peltoniemi, J.I., Kaasalainen, S., Näränen, J., Matikainen, L. and Piironen, J.. 2005 Measurement of directional and spectral signatures of light reflectance by snow. IEEE Trans. Geosci. Remote Sens., 43(10), 22942304. Google Scholar
Piironen, J., Muinonen, K. and Keränen, S.. 2000 Backscattering of light by snow: field measurements. In Verstraete, M., Meneti, M.M. and Peltoniemi, J., eds. Observing land from space: science, customers, and technology. Dordrecht, Kluwer Academic Publishers, 219228. Google Scholar
Poulet, F., Cuzzi, J.N. French, R.G. and Dones, L.. 2002 A study of Saturn’s ring phase curves from HSTobservations. Icarus, 158(1), 224228. CrossRefGoogle Scholar
Shkuratov, Yu. and 9 others. 2002 The opposition effect and negative polarization of structural analogs for planetary regoliths. Icarus, 159(2), 396416. CrossRefGoogle Scholar
Strozzi, T. and Mätzler, C.. 1998 Backscattering measurements of alpine snowcovers at 5.3 and 35 GHz. IEEE Trans. Geosci. Remote Sens., 36(3), 838848. Google Scholar
Van Albada, M.P. and Lagendijk, A.. 1985 Observation of weak localization of light in a random medium. Phys. Rev. Lett., 55(24), 26922695. Google Scholar
Verstraete, M. and Pinty, B.. 2001 Introduction to special section: modeling, measurement, and exploitation of anisotropy in the radiation field. J. Geophys. Res., 106(D11), 11,90311,908. Google Scholar
Warren, S.G. 1982 Optical properties of snow. Rev. Geophys. Space Phys., 20(1), 6789. Google Scholar
Wiersma, D.S., van Albada, M.P. B.A., van Tiggelen and Lagendijk, A.. 1995 Experimental evidence for recurrent multiple scattering events of light in disordered media. Phys. Rev. Lett., 74(21), 41934196. CrossRefGoogle ScholarPubMed
Yoon, G., Ghosh Roy, D.N. and Straight, R.C.. 1993 Coherent backscattering in biological media: measurement and estimation of optical properties. Appl. Opt., 32(4), 580585. Google Scholar
Figure 0

Fig. 1. The small-angle laser goniometer. Retroreflection from the sample is observed through the beam-splitter by the chargecoupled device (CCD) camera (left). Non-zero phase angles (source-target-detector) are obtained by rotating the beam-splitter (right). The distance from beam-splitter to sample was 12cm, and from the CCD to the sample 120 cm.

Figure 1

Table 1. Data from 2004. Snow sample characterization (approximate) and backscattering properties from the fits of the exponential function (see section 3.1)

Figure 2

Table 2. Data from 2005. Snow sample characterization (approximate) and backscattering properties from the fits of the exponential function (see section 3.1). R = sample taken from the snowpack on top of the river ice; T = surface tilted. The approximate sample thickness was 4 cm, except for the melted/refrozen samples for which it was 3cm for samples taken between 1 and 3 March and 3.5cm for both samples taken on 13 March

Figure 3

Fig. 2. Refrozen layers from the 2004 campaign measured (a) below (11 March 2004) and (b) above 0°C (19 March 2004).

Figure 4

Fig. 3. (a) Melting and (b) frozen surface grains measured on 2 March 2005. The grain images were taken with a millimeter scale (more clearly visible in the bottom left image in Fig. 4).

Figure 5

Fig. 4. (a) Melting grains (same sample from 2 March 2005 as shown in Fig. 3a) and (b) same sample after refreezing. The millimeter scale visible in the bottom left image is the same in all the grain images.

Figure 6

Fig. 5. (a) Melting new snow produces a stronger backscatter effect than does (b) the refrozen sample (both from 3 March 2005). Both surfaces were tilted to elimintate the effect of specular reflection.

Figure 7

Fig. 6. (a) Refrozen (rounded) grains (same sample as shown in Fig. 5b) and (b) grains of a coarser structure taken from 40cm depth in the snowpack (also measured on 3 March 2005). The difference in the smoothness of the phase curve is visible (cf. also Fig. 9).

Figure 8

Fig. 7. Effect of compression on surface (needle) grains: (a) slightly compressed (7 March 2005) and (b) compressed overnight (8 March 2005).

Figure 9

Fig. 8. Comparison of the effect of grain size: (a) 0.7–1.0mm rounded grains from 5–10cm depth in the snowpack (12 March 2005) and (b) same sample slightly melted and refrozen (13 March 2005).

Figure 10

Fig. 9. Effect of grain-size: (a) 0.2–0.5mm surface grains from previous day’s snowfall (8 March 2005) and (b) 0.1–0.3mm surface grains from previous day’s snowfall (8 March 2005; same sample as shown in Fig. 7b).

Figure 11

Fig. 10. Larger grain sizes (than those shown in Figs 8 and 9): (a) depth hoar with grains of 2mm from about 80cm depth in the snowpack (4 March 2005) and (b) a similar sample with grains of 2–3mm (10 March 2005).

Figure 12

Fig. 11. Comparison with data from 2004: (a) aged snow from about 40 cm depth in the snowpack and (b) depth hoar from 60 cm depth in the snowpack. The samples (both measured on 23 March 2004) are photographed in a sample cup of about 11cm diameter.

Figure 13

Fig. 12. Data from 2004 measurements: (a) new snow from the surface and (b) aged snow from about 20 cm depth in the snowpack (both samples measured on 23 March 2004). The effect of rounded grain structure (compared with fresh fallen grains of a finer structure) on the smoothness of the phase curve is also clear in this case. The sample grains are photographed in a sample cup of about 11cm diameter.