Introduction
Snow plays an important role in hydrological and climatic applications in mountainous regions. Field measurements of snow properties provide point observation but are restricted to limited locations due to rough terrain and harsh weather conditions. Remote sensing provides information about snow properties, which conventionally is not available for large areas due to typical terrain conditions. Snow shows high reflectance in the visible wavelength region and low reflectance in the shortwave infrared (SWIR) region. This spectral property of snow was used for snow-cover monitoring on a regular basis using the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor and the Advanced Wide Field Sensor (AWiFS) (Reference Hall, Riggs and SalomonsonHall and others, 1995, Reference Hall, Riggs, Salomonson, DiGirolamo and Bayr2002; Reference Kulkarni, Singh, Mathur and MishraKulkarni and others, 2006).
The reflectance of fresh snow is approximately 90% in the visible region and decreases at longer wavelengths. Freshly fallen snow almost immediately begins to compact and metamorphose, changing the snowpack characteristics. The albedo of fresh snow in the visible region of the spectrum remains high and decreases slowly with age of snow, but in the near-infrared region the albedo of aging snow decreases considerably as compared to fresh snow (Reference O’Brien and MunisO’Brien and Munis, 1975; Reference Warren and WiscombeWarren and Wiscombe, 1980; Reference Wiscombe and WarrenWiscombe and Warren, 1980). Reference Wiscombe and WarrenWiscombe and Warren (1980) showed that minute concentrations of small highly absorbing particles can lower snow albedo in the visible wavelengths by 5–15%.
Broadband albedo decreases when grain size increases as the snow ages. Melting of snow causes grains to grow clusters and behave as a large single grain (Reference Dozier, Schneider and McGinnisDozier and others, 1981; Warren, 1982). Reference GerlandGerland and others (1999) measured a maximum albedo more than 90% on Svalbard, Norway, before melt onset, and 60% after melt progressed in the spring when the snow was considered old but still clean. Many studies have related the snow reflectance in the near-infrared (NIR) region to grain-size variability using the 1.03 μm absorption feature by its continuum (Reference Dozier, Schneider and McGinnisDozier and others, 1981; Reference Nolin and DozierNolin and Dozier, 1993, Reference Nolin and Dozier2000). Snow-cover area, albedo, grain size, snow temperature, wetness and meteorological parameters are required for any snowmelt runoff model. Reference Painter, Dozier, Roberts, Davis and GreenPainter and others (2003) described the multiple end-member snow-covered area and grain-size (MEMSCAG) model which couples spectral mixture analysis with a radiative transfer to map sub-pixel snow-covered area and the grain size of the snow cover simultaneously, without using topographic data.
In undulating terrain like the Himalaya, and especially in the lower Himalayan region, atmospheric conditions change very fast. Temperature varies significantly between day and night, generating melt–freeze cycles. Very few previous studies have used a field spectroradiometer to understand snow characteristics in the Himalayan region. Spectroradiometer data of Beas basin are used in the present study. Contaminated snow shows a decrease in reflectance in the visible region, whereas grain size within the snowpack has more influence in the NIR region (Reference Singh, Negi, Babu, Kulkarni and SharmaSingh and others, 2005, Reference Singh, Negi and Kulkarni2008). Contamination and the effect of mixed objects like vegetation and soil on snow reflectance for the Himalayan region was reported by Reference Negi, Kulkarni and SemwalNegi and others (2009) where different indices were attempted to provide an understanding of snow pixel identification under different scenarios for broadband sensors. But the increasing number of airborne/satellite hyperspectral sensors like the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) provides an opportunity to measure the absorption spectrum and its utilities to improve applications (Reference Gao, Heidebrecht and GoetzGao and others, 1993; Reference PricePrice, 1998). Reference Green, Painter, Roberts and DozierGreen and others (2006) estimated the abundance of the three phases of water, i.e. water vapour, liquid water and ice, using AVIRIS over Mount Rainier, Washington, USA. Contamination and grain-size variability strongly influence the reflected albedo from snowpack and need to be characterized based on narrowband channel information. Selection and identification of proper bands is essential to reduce the redundancy in acquired information. Snowpack albedo, which is an essential input to a snowmelt model, is affected by contamination in the visible region and by grain size in the NIR region. In the present study, the possibility of identifying the contamination- and grain-size-based metamorphism was explored using field-based hyperspectral data available from the satellite platform. This study attempts to provide a set of information which can be used to identify the contamination and grain size of a snowpack simultaneously.
Study Area
Beas basin was selected as the study area for this experiment; it is in the Pir Panjal range, at 2000–4000 m a.s.l. It is characterized by moderate temperature and high precipitation. Snowfall takes place even in the lower ridges of the basin, and the approach of summer sees snowmelt take place. Due to moderate temperature, the snow cover changes into isothermal snowpack at 0°C and the winter period remains short. A location map of Beas basin is shown in Figure 1. Snow and Avalanche Study Establishment (SASE) observatories, which were used for snow reflectance and other collateral data collection in the present study, are established at Bahang (32.27° N, 77.158E; 2192ma.s.l.), Solang (32.32° N, 77.16° E; 4480 m a.s.l.) and Dhundi (32.36° N, 77.13° E; 3050 m a.s.l.) (Reference Mohile, Pethkar and DeshpandeMohile and others, 2004). These observatories receive a significant amount of seasonal snowfall on average between November and April (Bahang 53 cm, Solang 127cm and Dhundi 219 cm), and temperature remains high in this area compared to the other Himalayan ranges (at Bahang, T max = 21°C and T min = –0.44°C, at Solang, T max = 15°C and T min = –3.9° C and at Dhundi T max=11°C and T min=–4.7°C, from November to April, 1971–97) (Reference Mohile, Pethkar and DeshpandeMohile and others, 2004).
Field Set-Up and Data Collection
We used an ASD spectroradiometer, which collects radiometric quantities like radiance, reflectance and irradiance in the spectral range 350–2500 nm. It collects the data at a spectral resolution of 3 nm in the visible and NIR region, and 10–12 nm in the SWIR region. The bare-fibre 25° field of view was used. Figure 2 depicts the experimental set-up of the ASD instrument during field investigations. A crystal gauge with a magnifying glass was used to estimate the snow grain size. A thermometer was used for snow surface temperature and air-temperature measurements.
Methodology
A varying amount of contamination and grain size was considered. To observe the effects of contamination, contaminant-like soil was used in varying quantities. During the experiment the soil was dried, grinded and filtered. It was weighed into 1, 5, 20 and 30g packets which were equally distributed over the study area. The snow was moist, with a density of 240 kgm−3 and a grain size of 0.1–0.6 mm. The instrument sensor was kept at a height so that the view area on the surface had a diameter of 20 cm. Measurements were taken for clean snow and after spraying 1, 5, 20 and 30g soil equally on the surface. Reflectance data were collected using a bare-fibre sensor (25° field of view), and the average of 15 measurements was determined. Soil contamination was 0.0, 3.18, 15.92, 63.69 and 95.54 mgcm−2 on the snow-covered surface (Fig. 3). The duration of measurements was 10–15 min to avoid any other ambient influences due to change in sun-sensor viewing geometry and physical/meteorological conditions. Albedo was computed for each category of added soil contamination during the experiment. Snow grain-size data were collected in natural conditions and divided into fine (<0.5 mm), medium (<1.0mm) and coarse (1.0–2.0 mm) categories using a crystal gauge (Fig. 4). There is minimal human intervention in the study area, so it is assumed that no other disturbing effects exist.
Field-based hyperspectral data were binned at 10nm intervals using simple averaging which is compatible with the available hyperspectral satellite sensor’s spectral resolution and also with the spectral resolution of the field spectroradiometer at SWIR wavelengths with the assumption of square-wave response function. Principal components analysis (PCA) was carried out using the statistical software SPSS V-10 (SPSS Inc., 1999). Only those principal components with eigenvector values >1 were considered.
Different wavelength regions were chosen based on the influence of physical change on the snowpack characteristics. Position, strength and shape of a spectral curve can provide the information on smoothly varying spectral properties. The absorption peak at 1025 nm was chosen, as this is influenced by the grain-size variability. The relative strength, width and asymmetry of the absorption feature were computed after continuum removal for different wavelength regions. The objective of continuum removal is to isolate the properties of the absorption features from the overall reflectance properties. The relative strength of each band is computed using the formula (Reference Clark and RoushClark and Roush, 1984)
where S is the band strength, Rc is the reflectance of the continuum at the wavelength of Rb, and Rb the reflectance at the band minimum.
Band shape can be parameterized for the purposes of this systematic analysis by two simple parameters such as full width at half-maximum and asymmetry. The width of a feature is the absolute difference in microns between the right and left sides of absorption, where the reflectances are half the strength of the band. Asymmetry is defined as the base 10 logarithm of the sum of the reflectance over the number of channels to the right of the band minimum divided by the sum of the reflectance over the numbers of channels to the left of the band minimum. A zero value represents a perfectly symmetric band, a positive value right asymmetry and a negative value left asymmetry.
where R c→right represents the reflectance from centre to right, and R left→c the reflectance from left to centre, of the wavelength of the absorption peak. This parameterization, like band depth and asymmetry, allows rapid reduction in the dimensionality of a spectroradiometer dataset into discrete key absorption features of the obtained spectrum. The first and second derivatives were computed, which indicated any change in shape of the curve at particular wavelengths. Changes in percentage of reflectance were computed with contaminants (soil) and for varying grain size.
Results and Discussions
The influence of soil contamination and grain size on snow reflectance was considered, with a view to using this information in hyperspectral remote sensing. Figure 3 shows the reflectance behaviour of snow with a changing amount of contamination, and Figure 4 the grain-size variability of snow reflectance. A significant decrease in reflectance due to soil contamination is seen in the visible and NIR region, with a lesser decrease due to change in grain size. It is also observed that the absorption peak at 1025 nm almost diminishes for 63.69 and 95.54 mgcm−2 amount of soil contamination, whereas it is still significant for varying grain size. The relative strength at 1025 nm, after continuum removal, shows a decreasing trend in soil contamination, whereas a reverse pattern of increasing relative strength is observed for grain size (Fig. 5).
Tables 1 and 2 show the relative strength and asymmetric nature of the absorption peak in the NIR region. Relative strength decreases except for the 15.92 mg cm−2 soil contamination level at 1025 nm which can be explained with the help of Table 2. The absorption peak is right asymmetric at 1025 nm and was systematically asymmetric towards the left. Hence, the relative strength at 1000 nm shows a clearer picture of decreasing relative strength due to increase in soil contamination level. Asymmetry at 1085 nm is not very significant.
The relative strength at 1025 nm due to grain size shows an increasing pattern (Table 3). The absorption peak at 1025 nm shows a right asymmetric nature (Table 4). The shape is inclined towards the left side but only from fine grain to medium grain, and no difference was observed in the shape of the peak from medium to coarse grain. However, asymmetry at 1085 nm is systematically asymmetric towards the right. The peak shift is also seen very clearly in the first derivative of soil contamination and grain size (Fig. 6). Figure 6 shows a peak shift in the visible region due to the addition of contamination which is not present for varying grain size. As the contamination level increases from clean snow to highly contaminated snow, the shift moves from 575 nm to 735 nm. The width of the absorption feature was found to be within 6 nm range at 1025 nm from clean to highly contaminated snow. Albedo decreased as the contamination level increased. However, the albedo was drastically reduced (∼25%) for the already small amount of soil contamination, and further reduction in albedo was not significant (∼5%) for the higher amount of contamination. Albedo was reduced by ∼12% from fine to medium grain size, but there was no significant reduction from medium to coarse grain size (Fig. 7). Tables 5 and 6 show reflectance and percentage change in reflectance for a few optimum important bands for contamination and grain-size variability study. SWIR 1 and SWIR 2 show opposite patterns for contamination and grain-size metamorphism, especially for higher levels of contamination. They show a significant reduction in reflectance from fine to medium and coarse grain size.
Relative strength, asymmetry and peak shift show a squeezing of snow spectrum from the visible to the NIR region for contamination which is not present for varying grain size. This is an important observation which shows the possibility of identifying snow metamorphism due to contamination and grain-size variability on a temporal basis using hyperspectral remote sensing.
Conclusions
Radiometric data for different soil contamination and varying grain size were collected in controlled conditions in the Himalayan region. Relative strength, asymmetry, width, first derivative, albedo and percentage change in reflectance were used to analyse the field data. Relative strength, asymmetry and first derivative of reflectance show that as the soil contamination increased, a squeezing of snow spectrum was observed in the visible to NIR region of spectrum. Peak width was not significant for soil contamination and grain size (within 5–6 nm) at 1025 nm. Different levels of contamination resulted in a decrease in relative strength at 1025 nm and 1000 nm. Increasing grain size caused an increase in relative strength at 1025 nm in the Himalayan region, as has also been reported by Reference Nolin and DozierNolin and Dozier (2000) using AVIRIS data. A combination of these two observations is helpful to differentiate contamination- or grain-size-based snow metamorphism using temporal information. Both soil contamination and grain size have a right asymmetric absorption peak at 1025 nm, but soil contamination shows a systematic tendency towards left asymmetry of curve shape with increase in soil contamination level. Albedo significantly decreased due to the already low level of contamination (∼25%) and did not decrease much further for higher soil contamination levels. This observation suggests that contamination can initially drastically reduce the albedo but almost saturates if the snow is further contaminated.
This study shows that the 1025 nm absorption peak is useful not only for detecting grain-size variations but also for studying soil contamination. In the future, relative strength, asymmetry, peak shift, albedo and percentage change in reflectance may be used as a classification scheme to study/identify the cause of snow metamorphism, i.e. contamination/grain size, using hyperspectral remote sensing.
Acknowledgements
We thank R.R. Navalgund (Director, Space Applications Centre) and J.S. Parihar (Deputy Director, Remote Sensing Applications Area) for their kind support during this study. We also thank Dr Ajai (Group Director, MESG) and T.J. Majumdar (Head, Earth Sciences and Hydrology Division (ESHD)/MESG) for their encouragement. We gratefully acknowledge the help of Snow and Avalanche Study Establishment (SASE) Team Members during the field investigations, and of anonymous reviewers for their critical evaluation of and suggestions on the manuscript.