This study applied remote sensing techniques to the study of water color in Himalayan glacial lakes as a proxy of suspended solid load. In situ measurements gathered in 5 lakes in October 2014 during satellite data acquisition enabled the characterization of water reflectance and clarity and supported image processing. Field data analysis led to a distinction between 3 water colors and a consequent lake water color classification on a regional scale from Landsat-8 data previously corrected for atmospheric and adjacency effects. Several morphometric parameters (lake size and shape, distance between lake and glacier) were also computed for the lakes thus classified. The results showed spatial and temporal variations in lake water color, suggestive of relationships between glacier shrinkage and the presence of brighter and more turbid water. A finer-scale analysis of the spatial variability of water reflectance on Chola Lake (based on GeoEye-1 data captured on 18 October 2014) showed the contribution of water component absorption from the inflow. Overall, the findings support further research to monitor Himalayan lakes using both Landsat-8 and Sentinel-2 (with its improved resolutions).
Introduction
Climate change is known to affect glacier mass balances and lead to a net loss of ice, due mainly to higher temperatures and declining precipitation (Salerno et al 2008; Wu et al 2012; Thakuri et al 2014). The rising melting rates induced by climate change can have various effects, including an increase in the number and size of glacial lakes (Thompson et al 2012) and the creation of potentially dangerous lakes with a risk of outburst floods (Bolch et al 2008). Meltwater from glaciers is the main factor determining lake volume in the region, but it also carries glacial sediment into the lakes, increasing their turbidity (eg Irwin 1974; Robinson and Matthaei 2007). Glacial meltwater can alter chemical processes in glacial lakes, for instance increasing water conductivity and ionic concentrations (Thies et al 2007; Lami et al 2010). Sediment-rich water becomes less transparent, inhibiting primary production and making it hard for higher organisms in the food web to survive in the lake (Koenings et al 1990). All these changes have been underway in mountainous regions where there are glaciers, and it is important to monitor the close interactions between glaciers and glacial lakes.
For glacial lakes in remote and relatively inaccessible areas, it is impractical to characterize these interactions by means of conventional fieldwork. Remote sensing has become a powerful and convenient option for ongoing synoptic measurements and their retrospective analysis. Satellite data have been widely used to study lakes at high altitudes, in particular in the Himalayas. Recent examples of the use of satellite data on the Himalayan region can be found in Raj and Kumar (2016), who describe the evolution of glacial lakes in the Uttarakhand Himalaya; Raj et al (2013), who produced an inventory of the glacial lakes in the Sikkim Himalaya; Giardino et al (2010), who mapped the suspended particulate matter in lakes in the Mount Everest region; Quincey et al (2007), who identified glacial lake hazards in the Himalayas; and Kargel et al (2005), who found a supraglacial lake on Rongbuk Glacier.
In the present study, satellite data were used to examine the optical characteristics of Himalayan lakes in terms of their radiometric properties, such as reflectance or water color, identifiable by remote sensing. The size of the particles in suspension in the lake and their numerical density change, depending on the influx of water from the glacier; this affects light interaction and lake color. Higher turbidity levels make these lakes look gray, whereas lakes with a lower input of particles in suspension from the glacier appear dark blue (Kargel et al 2005). Lakes with large amounts of suspended sediment generally have a higher visible reflectance than clear waters (Bukata et al 1995). Once satellite data have been corrected for atmospheric effects, changes in water reflectance signals from visible to near-infrared wavelengths depend on differences in the water components in the upper layer of the water body. These variations can be retrieved by applying specific algorithms (eg Odermatt et al 2010; Gholizadeh et al 2016) built on empirical or semiempirical methods that use reflectance at appropriate wavebands as correlates, or by means of such semianalytical and quasi-analytical approaches as spectral inversion, which relies on matching spectral data with bio-optical forward models. Whatever the approach, the calibration and validation of such models take advantage of the availability of site-specific data gathered in the field. For instance, in situ radiometric data have been used to tune suitable regional algorithms for measuring water quality parameters (eg Hommersom et al 2011) and to assess the performance of atmospheric correction methods adopted to obtain water reflectance data from satellite images (eg Brando et al 2016).
The aim of the present work was to characterize the water reflectance in the visible near-infrared (VIS-NIR) region of glacial lakes in the Mount Everest region (Nepal) as a proxy of their color, using in situ radiometric data and satellite acquisitions. During fieldwork in October 2014, the water reflectance of 5 Himalayan lakes was measured at the same time as data were acquired by 2 satellites, GeoEye-1 and Landsat-8. The data from GeoEye-1, with its higher spatial resolution (2 m), enabled water color variability to be assessed within a single lake, whereas the Landsat-8 image at 30 m resolution documented the color variability between different lakes in the region. The classification of the lakes based on the Landsat-8 data is analyzed here in relation to their morphometric parameters (lake size and shape, distance of lake from glacier), and considering the lake color classification proposed by Giardino et al (2010), in an effort to explain recent lake color variations in the region.
Study area
The study area is between 27.6–28.3°N and 86.3–87.1°E, and includes both Nepali and Chinese territory (Figure 1). In the center of this area is the Sagarmatha National Park (SNP), the highest protected area in the world, a UNESCO World Heritage Site since 1979. The SNP is situated in the Solu-Khumbu district and covers 1141 km2 in northeastern Nepal. The area is dominated by Mount Everest and 7 other peaks over 7000 m high, and it is the main source of freshwater for the Nepali people downstream.
FIGURE 1
Map of the study site. Lakes inside SNP are coded with cadastral numbers (1–200) based on Tartari et al (1997, 2008), and other lakes are coded with randomly assigned numbers (beginning at 1001). (Map by the authors)
![i0276-4741-37-1-16-f01.tif](ContentImages/Journals/mred/37/1/MRD-JOURNAL-D-15-00052.1/graphic/WebImages/i0276-4741-37-1-16-f01.jpg)
The glaciers are largely covered in debris, and this alters the energy exchanges between ice and atmosphere (Mattson et al 1993). The area's glacial lakes have different origins and are differently connected with the glaciers: some just formed as melting ponds above the glacier's tongue (supraglacial); some are moraine dammed and emerge directly from the glacier's front (proglacial); and some are not connected directly to a glacier but share the same basin (cirque). Proglacial lakes are considered the most dangerous in the event of major melting of the glaciers behind them because of the risk of outburst floods—as already seen in some lakes in the SNP, such as Nare Drangka in 1977 (Buchroithner et al 1982) and Dig Tsho in 1985 (Vuichard and Zimmermann 1987). The lakes in the study area lie between approximately 4000 and 5700 m above sea level, well above the tree limit, and they are fed mainly by perennial snow and ice.
The climate in the SNP is driven by the Indian monsoon, with successions of dry and wet seasons. Most of the annual precipitation (80%) falls between June and September. The winter is normally dry, though occasional cyclones can cause heavy snowfall (Ueno et al 2001). One such exceptional event occurred during fieldwork in October 2014, when Cyclone Hudhud, originating in the Andaman Sea (Bay of Bengal), intensified on its way toward India, developing into a severe storm and causing several deaths in India (Andhra Pradesh) and Nepal (Annapurna region).
Five lakes in the Khumbu and Imja valleys with different basin morphologies were selected for fieldwork (Table 1 and Figure 1). These lakes are here identified by a univocal sequential code, which is the lake cadastral number (LCN), as assigned by Tartari et al (1997, 2008) in their inventory and surface change detection studies. Lakes LCN9 and LCN10 are located above the CNR-EvK2 Pyramid observatory and constitute a cascade system, with the upper lake feeding the lower one. Lake LCN24 (Chola) is dammed by a lateral moraine separate from the lake basin that conceals the outlet, so the water emerges some meters below the lake. Lake LCN31 is a semicircular lake formed at the end of the Duwo glacier moraine. Lake LCN161 (Imja) began as a small meltwater pond in the 1960s, but has developed into a proglacial lake, raising concern about its stability and the risk of outburst flooding (Bolch et al 2008; Byers et al 2013). The Imja Glacier's ice tongue protrudes into the water, directly feeding the lake. It is a typical proglacial lake, gradually expanding over the years as the glacier melts and retreats. It now poses a threat to the local population (Bajracharya et al 2007) because its basin consists of low sloping glacier surfaces and it collects the meltwater from 3 surrounding glacier tongues, so a moraine break or snow/ice avalanches could have catastrophic consequences. For the other lakes within the study area (see Figure 1), most of those within the SNP and characterized by previously described basin morphologies, the overall assessment concerning their risk of outburst is then conducted by using satellite observations.
TABLE 1
Morphological characteristics of sampled lakes. Morphology is based on Tartari et al (2008).
![i0276-4741-37-1-16-t01.eps](ContentImages/Journals/mred/37/1/MRD-JOURNAL-D-15-00052.1/graphic/WebImages/i0276-4741-37-1-16-t01.jpg)
Material and methods
Field data
Fieldwork was conducted from 14 to 20 October 2014 to collect radiometric and transparency data on 5 lakes between 4532 and 5067 m above sea level (Table 1). Radiometric measurements were taken with a field spectroradiometer (WISP-3, WaterInsight), and water transparency was measured using a Secchi disk (SD) to determine the depth to which visibility remained clear. WISP-3 simultaneously records radiance (40° sky and 40° water) and irradiance spectra in the 400–800 nm range (Hommersom et al 2012). To avoid any influence of the lake bottom on the radiometric measurements, sampling stations were located in areas of deep open water reached with an inflatable boat. For each station, a water sample was collected and analyzed with an iQwtr instrument (Ghezehegn et al 2014), which uses an innovative, low-cost, user-friendly field application to estimate surface water turbidity from photos taken with smartphones (BlueLeg Monitor and WaterInsight).
Satellite images
Two images were acquired by the GeoEye-1 and Landsat-8 Operation Land Imager (OLI) satellite sensors on 18 and 29 October 2014, respectively. The GeoEye-1 acquisition was concurrent with in situ measurements obtained on Lake LCN24. Both sensors are multispectral, with typical blue, green, red, and infrared channels.
For both images, radiometric calibration and correction for atmospheric effects were performed to obtain surface reflectance data. In particular, radiometric correction of the OLI image was based on the specific coefficients suggested by Pahlevan et al (2014) for water applications. The vector version of the 6S code (Second Simulation of a Satellite Signal in the Solar Spectrum) (Vermote et al 1997) was applied to the satellite images to correct for atmospheric effects and measure water reflectance. The atmospheric profile was established from user-defined ozone and water vapor values (Ozone Over Your House, NASA and AERONET station EVK2_CNR), and the continental aerosol model was adopted as the most representative of the study site. The aerosol optical depth was derived from AERONET data (EVK2_CNR station) and set to 0.05. The 6S-derived reflectance values were further corrected for adjacency effects based on the residual between the 6S-derived reflectance from the GeoEye-1 image and the synchronous WISP-3 in situ data obtained for Lake Chola (LCN24).
Three thresholds were chosen to classify the lakes by water color on the corrected Landsat-8 data, based on intensity of water reflectance (rw) and on the observed desaturation of lake colors from dark blue to turquoise to gray (according to Giardino et al 2010): rwVIS_NIR < 10% for clear blue waters (as in lakes LCN10 and LCN31), 10% < rwVIS_NIR < 20% for turquoise lakes (like LCN24); and rwVIS_NIR > 20% for gray, sediment-rich waters (as in LCN161). The value of rwVIS_NIR was calculated from the sum of the rw values (Bukata et al 1995) in the typical VIS-NIR band settings of multispectral sensors like the GeoEye-1 and OLI. For the classified lakes, the shape index proposed by Frohn (2006)—which has been judged to perform well for raster-derived data (Bogaert et al 2000)—was used to test the relationship between a lake's water color and its shape. The shape index defined by Frohn ranges between 0 (for a square lake) and 1 (for a lake more complex in shape). The distances between lake and glacier were also computed on the basis of the Randolph glacier inventory (Arendt et al 2014) using instruments commonly available in geographic information systems.
On a finer scale, the rw values measured in the 545 nm GeoEye-1 band were used to analyze spatial variations in the LCN24 lake, because this spectral band is known to be affected by organic and inorganic sediment in suspension (Kirk 2011).
Results
The WISP-3 rw spectra and SD measurements for the 5 surveyed lakes are shown in Figure 2. Overall, 3 groups of spectra with increasing degrees of brightness were clearly distinguishable by water color. LCN9 and LCN10 are spectrally identical highly transparent lakes with dark blue waters and similar SD depths (about 8–9 m). LCN31 is also a clear lake with a bluish color and reflectance values decreasing from blue toward longer wavelengths because of water absorption. LCN24 and LCN161 have rw spectra typical of higher turbidity levels, as confirmed by the SD depths lower than 1 m. The highest reflectance was measured in LCN161, which looked gray during fieldwork observations, with clearly visible patterns of suspended sediment. Field measurements of total suspended matter concentrations obtained in October 2008 (Giardino et al 2010) support the observed radiometric differences: the total suspended-matter values ranged from 0.63 g m−3 for lake LCN10 to 102 g m−3 for lake LCN161, with intermediate values of 1.7 g m−3 and 6.69 g m−3 for lakes LCN31 and LCN24, respectively. The WISP-3 rw spectra for 4 lakes (all except LCN9) were comparable with the ALOS-derived reflectance of Figure 2 in Giardino et al (2010). The average rw in situ data in the VIS-NIR were then compared with the SD values, and a correlation was identified (r2 = 0.85). The transparency values obtained with the SD (cf. Figure 2) and the iQwtr also revealed a good correlation (r2 = 0.94).
FIGURE 2
Reflectance (rw) spectra measured with the WISP-3 field spectroradiometer during fieldwork. Note that LCN9 and LCN10 are spectrally identical. Measured SD depths are indicated close to each spectrum. (-) means dimensionless.
![i0276-4741-37-1-16-f02.tif](ContentImages/Journals/mred/37/1/MRD-JOURNAL-D-15-00052.1/graphic/WebImages/i0276-4741-37-1-16-f02.jpg)
When compared with in situ radiometric measurements, the 6S-derived reflectance showed the same trend, but with a higher magnitude. This difference in spectral magnitude was attributed to adjacency effects due to fresh snow falling on 14 October 2014, a few days before the Landsat-8 and GeoEye-1 data were acquired. This effect increased the signal by 39% (more in the blue band, and less toward the NIR), a finding generally consistent with the report from Bélanger et al (2007) of an adjacency effect due to sea ice resulting in a similar overestimation of the water-leaving reflectance. To correct for adjacency effects, the 6S-derived reflectance obtained with the Landsat-8 and GeoEye-1 data was adjusted to the WISP-3 in situ data for Lake Chola (LCN24), taken for reference because the field and satellite data collection coincided.
The 3 rwVIS_NIR thresholds defined in the previous section were applied to OLI data and corrected for atmospheric and adjacency effects to classify the lakes. Figure 3 shows the outcome of the classification procedure for 119 lakes. Of the lakes, 52% have mainly (>50% of pixels) gray waters, 24% have mainly blue waters, and 21% have mainly turquoise waters; 3% have mixed colors (half blue and half turquoise). The lakes in Nepal (96 of 119) have a greater diversity of water colors, largely gray (45%) with fewer blue (27%) and turquoise (25%) (the remaining 3% were classified with mixed colors), whereas the Chinese lakes (23 of 119) have a preponderance of gray waters (83%) and even fewer blue (13%) or turquoise (4%). When the classification obtained on the OLI image was compared with one previously obtained on an ALOS image (Giardino et al 2010) for 116 lakes in the area (which is the maximum number of lakes commonly detectable—i.e. included in both scene swaths—and free from cloud cover in both images), 71% of the lakes revealed differences in the percentage composition of their water color over time, the waters tending to become clearer in 38% of the lakes and more turbid in 62%.
FIGURE 3
Water color classification of the 119 investigated lakes. Each bar shows the percentage of pixels classified as blue, turquoise, or gray water for 1 lake. Lake codes are on the left and right; asterisks indicate lakes located in China.
![i0276-4741-37-1-16-f03.tif](ContentImages/Journals/mred/37/1/MRD-JOURNAL-D-15-00052.1/graphic/WebImages/i0276-4741-37-1-16-f03.jpg)
Figure 4 compares water color and lake size and shape. Lakes with turquoise and gray water were generally larger than the blue lakes. The lakes with gray water could have a surface area in excess of 1 km2, indicating relatively large water volumes capable of causing damage in the event of an outburst flood. The shift from blue to gray water coincided with a declining likelihood of a lake being square (ie shape index values approaching 1). Lakes with gray and turquoise water were closer to their glaciers than those with blue water (median <0.25 km for gray-water lakes and >0.5 km for blue-water lakes). Lakes containing gray water also tended to be more elongated and often lay at the end of a glacier tongue (eg LCN161, 1040, 1056).
FIGURE 4
Distribution of water color by lake size and shape index (Frohn 2006). The shape index ranges between 0 (circular shape) and 1 (irregular shape). Box plots show 25th and 75th percentiles; error bars show 10th and 90th percentiles; diamonds show outlier values.
![i0276-4741-37-1-16-f04.tif](ContentImages/Journals/mred/37/1/MRD-JOURNAL-D-15-00052.1/graphic/WebImages/i0276-4741-37-1-16-f04.jpg)
Water color variability was investigated on a finer scale for Lake Chola (LCN24), which has predominantly turquoise waters, captured by the GeoEye-1 acquisition concurrently with our in situ measurements. Figure 5A shows a map of the rw values in the 545-nm band, and Figure 5B shows the correspondence between the rw spectra measured in situ and derived from the satellite data. The rw at 545 nm is lower near the inflow, higher along the southwestern shore, and intermediate along the northeastern shore—constituting 3 zones that can be roughly identified on the image. The average spectral signatures of these zones, measured in situ and derived from satellite data, are plotted in Figure 5B. The southwestern and northeastern parts of the lake have similar spectra, with a slight variability in magnitude (probably due to differences in their suspended sediment loads). The waters near the inflow show slightly lower rw values for the blue-green wavelengths, suggesting the presence of light-absorbing substances entering the lake.
FIGURE 5
Lake Chola (LCN24): (A) map showing rw values in the 545-nm band, captured by GeoEye-1; (B) rw values in the 400–800-nm range of 3 lake zones, with different reflectance characteristics, captured by GeoEye-1, and rw values in the 400–800-nm range measured with WISP-3 in the location marked with a star. (-) means dimensionless.
![i0276-4741-37-1-16-f05.tif](ContentImages/Journals/mred/37/1/MRD-JOURNAL-D-15-00052.1/graphic/WebImages/i0276-4741-37-1-16-f05.jpg)
Discussion
The analysis of lake color in relation to size and shape showed differences that can be explained by the different origins of the lakes. Those with a more circular shape are usually cirque lakes, with a well-defined basin modeled by ancient glacier action and currently fed by river water or rain; lakes with an ellipsoid shape are more typically supraglacial and proglacial, formed directly by the melting of glacier tongues and adapting to their shape. Lakes with gray and turquoise water were found closer to glaciers than blue-water lakes, and their proximity to a glacier points to a possible link between lake water color and glacier melting. The elongated shape and increased surface area of a lake may be indicators of a glacial lake outburst flood risk. ICIMOD (2011) identified 10 lakes in Nepal with a high level of surveillance priority due to their physical parameters and socioeconomic factors. Our study found most (7 of 10) of these potentially dangerous lakes (LCN161, 1021, 1028, 1033, 1048, 1053, and 1059) to be gray in color and located 0 to 75 m from their respective glaciers, with a shape index higher than 0.4.
Some 116 lakes included in our study were also classified by water color in an earlier study (Giardino et al 2010); in general, these became more gray over the intervening 6 years. One possible reason for this trend is glacier shrinkage, which has been amply discussed as a consequence of global climate change. Glacier surface area in the SNP shrank by 261 km2 from 1980 to 2010 (ICIMOD 2014, 2015). The lowest elevation of glacier tongues has also shifted upward for 54% of the glaciers in the region, by an average of 40 m, and 97% of the glaciers have become thinner, on average by 7 m. All these data point to a retreating trend of the glaciers in the study area from 1980 to 2010 and support our hypothesis that glacier shrinkage influences water color in the region's glacial lakes.
With its high spatial resolution, the GeoEye-1 generated data capable of showing variations in water reflectance within a single lake (Figure 5). The mountain slopes on the southwestern and northeastern shores of this lake are quite different, the former mainly consisting of rocks and debris, the latter colonized by sparse vegetation (alpine pastures and shrubs). When the snow melts, the surface flow of water tends to carry some of the debris into the lake, and runoff erosion affects mainly less consolidated slopes, which are usually not covered by vegetation (Fattet et al 2011; Huang et al 2013), as in the case of the southwestern slope. This could explain why the reflectance values were higher on the southwestern side of the lake.
Conclusions
The present study combined in situ and satellite data to characterize the water color of Himalayan lakes in the SNP and surrounding area. In situ radiometric measures (WISP-3) and SD measurements enabled a distinction between 3 water color classes—blue, turquoise, and gray—coinciding with decreasing levels of water clarity and increasing levels of brightness. The comparison between the transparency measurements obtained with the SD and the iQwtr instrument seems promising (even though only 5 samples were available), and would support the future use of innovative app-based instrumentation in fieldwork.
The different spatial resolutions of satellite sensors proved adaptable to different aims, from more detailed analysis of the spatial distribution of the water characteristics in a given lake to mapping of the lakes on a regional scale. Most of the lakes in the study area were classified as having predominantly gray waters, and these lakes also seemed to have larger surface areas and elongated shapes, probably extending directly from the glacier tongue that feeds them—a combination that raises the risk of outburst flooding. When direct surface melting gives rise to small ponds on a glacier's surface that evolve into proglacial lakes, the water of the lake is dammed by the glacier's moraine, which cannot retain it securely in the event of an avalanche or damage to the moraine. Our comparison of data obtained in 2008 and 2014 showed a tendency for lake water to become less clear (more turquoise or gray than blue). The factors behind this change need to be investigated and discussed, but we already know that the past 40 years or so have seen the glaciers in the area (and specifically within the SNP) shrink in size and thickness, and the lowest edge of more than half of them has shifted to higher altitudes.
Satellite remote sensing is a good way to track glacier–lake interactions in scarcely accessible areas such as the Himalayas, but it has some limitations. One concerns the spatial resolution, which restricts the investigation of lakes smaller than the minimum mappable unit (image pixel). Another concerns the orography and adjacent bright lands, which hamper the extraction of pure signals reflected by water because of areas of shadow. With the Sentinel-2 now in orbit (capable of 10-m resolution and a 10-day revisiting time, which will drop to 5 days once the second Sentinel-2 is in orbit), better data are available for capturing the spatial and temporal variability of water color in the Himalayan region. Its 10- and 20-m spatial resolution bands are an important innovation, because spatial resolution remains one of the main factors limiting the application of satellite remote sensing to freshwater ecosystems (Ozesmi and Bauer 2002; Hestir et al 2015).
ACKNOWLEDGMENTS
This study was carried out within the framework of the Ev-K2-CNR Project in collaboration with the Nepal Academy of Science and Technology, as foreseen by the Memorandum of Understanding between Nepal and Italy and supported by contributions from the Italian National Research Council and the Italian Ministry of Foreign Affairs. The study was also part of the GLaSS project (EC FP7, project number 313256). We thank G. P. Gobbi, Principal Investigator of the AERONET CNR-EVK2 station, and M. Laanen, S. Ghebrehiwot, and A. Hommersom for providing the iQwtr instrument. A special thanks to A. Lami and to G. Verza for supporting the fieldwork and to all the Pyramid Observatory local staff for their invaluable help during measurements and trekking. We thank 2 anonymous reviewers for their insightful comments, which helped strengthen this manuscript.