Hostname: page-component-78c5997874-v9fdk Total loading time: 0 Render date: 2024-11-19T02:16:06.371Z Has data issue: false hasContentIssue false

Toward the satellite monitoring of glacier changes on Mount Kenya

Published online by Cambridge University Press:  20 January 2017

Stefan Hastenrath*
Affiliation:
Department of Meteorology, University of Wisconsin-Madison, 1225 West Dayton Street, Madison, WI 53706, U.S.A.
Rights & Permissions [Opens in a new window]

Abstract

For a pilot study of change detection of small tropical mountain glaciers from SPOT satellite, Mount Kenya was chosen because a map at scale 1 : 5000 in a local coordinate system, based on aerial photography in September 1987, provided a control for mapping of horizontal ice extent, and because the precipitous relief offered a particular challenge for technique development. A terrain model in compatible map projection was constructed from the topography of the 1987 map. A SPOT subscene of February 1988 was geocoded and terrain corrected to conform with the terrain model and identifiable terrain features. Another 26 control points established on the ground were entered in the digital array. The pixels were re-sampled to a nominal resolution of 10 m. This digital map was output at scale 1 : 5000. The glacier boundaries and the location of the terrain control points agreed with the September 1987 map within the mapping accuracy commensurate with the 1 : 5000 scale. This satellite-based mapping provides an internally consistent reference for determining the changes of Mount Kenya’s glaciers from February 1988 to a SPOT sampling intended for the mid-1990s. This pilot study further indicates the feasibility of glacier change detection on various other high mountains of the tropics.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1993

Introduction

In the context of global climatic change, the variations of tropical mountain glaciers merit particular attention. The glacierized mountain regions of the tropics are not easily accessible and systematic mappings from terrestrial or airborne photogrammetry have been rare. Monitoring from satellite appears an obvious option, but the spatial resolution of LANDSAT was far too coarse to assess changes in the small tropical glaciers. The considerably finer resolution of SPOT opens the prospect for applications to glacier change detection in the tropics. With this motivation Mount Kenya was chosen for a pilot study, because our extensive earlier work (Reference Hastenrath and KrussHastenrath, 1984; Hastenrath and others, 1989; Reference HastenrathHastenrath and Kruss, 1992a, b) has examined the sensitivity of its glaciers to climatic variations and has created an unusual topographic basis, and because the extreme local relief posed particular challenges to technique development. The present paper reports the results of a satellite-based digital mapping of Mount Kenya’s glaciers, and discusses the potential of this approach for the small glaciers on Mount Kenya and other high mountains of the tropics.

Background

Mount Kenya has been an object of our research on tropical glacier and climate variations over two decades (Reference Hastenrath and KrussHastenrath, 1984). Field work since the mid 1980s included the gradual establishment of a network of ground-control points all around the glacierized peak region. This used electromagnetic distance measuring equipment (EDM) and optical theodolite (Reference Hastenrath and KrussHastenrath and others, 1989) and linked up with topographic control points created by the IGY Mount Kenya Expedition in 1957–58 (Reference CharnleyCharnley, 1959). We marked the control points in the terrain with white paint in preparation for aerial photography flown on 3 September 1987. Of the total of 26 ground control points 14 were used for aerial triangulation and stereo-plotting. A map was produced at scale 1:5000 (Reference Hastenrath and KrussHastenrath and others, 1989; World Glacier Monitoring Service, 1992), conforming in scale and area with an earlier map based on terrestrial photogrammetry in January 1963 (Forschungsunternehmen Nepal-Himalaya, 1967). Both maps used a local coordinate system stemming from the IGY Mount Kenya Expedition. A simplified version of our 1987 map reduced to scale 1 : 15 000 is shown in Figure 1. A central frame from the 3 September 1987 aerial survey is reproduced in Figure 2 at approximately the same scale. Noteworthy with regard to Figures 1 and 2 is the extreme local relief of about 1 km within a 2 × 2 km horizontal area.

From a comparison of the 1963 and 1987 maps it was possible to determine quantitatively not only the area decrease but also the thickness changes of the glaciers on Mount Kenya over the past quarter century (Reference Hastenrath and KrussHastenrath and others, 1989). The ice thinning amounted to about 15m for all glaciers, regardless of topographic location and thus exposure to solar radiation. The ice thickness change, in turn, was evaluated in terms of the glacier heat budget, and on these grounds we estimated the sensitivity of the glaciers to changes of the atmospheric environment that may be steered by the greenhouse effect through various lines of causality (Reference HastenrathHastenrath and Kruss, 1992a, b). In particular, we concluded that a slight increase of astmospheric humidity — itself largely controlled by conditions over the Indian Ocean — was instrumental in reducing sublimation and thus making energy available for enhanced melting.

Fig .1. Map of Mount Kenya’s glaciers on 3 September 1987. The key below lists for the north and south sides of the mountain and from west to east the names of the glaciers indicated in the map by number. Triangles indicate terrain control points used in the 1987 mapping, and open circles other control points surveyed and marked in the terrain, with letter code referring to Reference Hastenrath and KrussHastenrath and others (1989). Crosses represent terrain features identified for geocoding of digital SPOT scene. Square dots mark the highest peaks (Batian 5199 m, and Nelion 5188 m). Height contours are at 200 m spacing and scale is 1:15 000 (source: Reference Hastenrath and KrussHastenrath and others, 1989). Dot raster denotes glaciers and dash raster lakes. North side: 13 Cesar, 14 Joseph, 16 Northey, 1 Krapf, 2 Gregory; south side: 10 Tyndall, 9 Heim, 8 Forel, 7 Diamond, 6 Darwin, 4 Lewis.

With such climatic objectives, remote sensing by satellite is of interest for determining area and area changes, if not elevations and thickness changes. The potential and limitations of LANDSAT imagery have been explored by Reference Allison and PetersonAllison and Peterson (1989) for New Guinea and by Young and Hastenrath (1991) for East Africa. The 80 m resolution of LANDSAT was found not suitable for area change detection of the small tropical glaciers.

Evaluation of Spot Scene

The pixel size of SPOT is 20 m for multi-channel, and 10 m for panchromatic scenes. Resampling allows sharpening to nominal resolutions of 10 and 5 m, respectively. These values should be seen in perspective with the 1963–87 terminus retreats of Mount Kenya’s glaciers of mostly between 50 and 150 m, and the scale of our 1987 map (1 mm on the map corresponding to 5 m in the terrain). From these magnitude considerations it was found desirable to map the ice extent on Mount Kenya in scale and coordinate system consistent with our map based on aerial photogrammetry (Reference HastenrathHastenrath and others, 1989), from a SPOT pass as close as possible to the 3 September 1987 date of our map. From a search of the SPOT files, a cloudless panchromatic shot of 3 March 1987 appeared most desirable, but this was found damaged in the SPOT archives. The second choice then was a cloudless multi-channel subscene of 4 February 1988, and this proved to be of good quality.

Fig .2. Air photograph of Mount Kenya’s glaciers on 3 September 1987, of flight (PHOTOMAP Kenya Ltd.) from which map was produced by aerial photogrammetry (Reference HastenrathHastenrath and others, 1989). Approximate scale 1: 15000.

In preparation for the processing of the SPOT subscene, a digital terrain model at 5 m cell size was constructed from the topography in our September 1987 map (Reference HastenrathHastenrath and others, 1989), and a compatible map projection was created (polyconic projection, Clarke 1880 Ellipsoid, consistent with ARC 1960 datum used in Kenya; central meridian 37°18’30ʺ;E, latitude origin 0°09’00ʺS, false easting 2700.00, false northing 2165.00). Prominent terrain features were then identified on the September 1987 map (Fig. 1) and in the February 1988 satellite image. On this basis, the digital subscene was geocoded and terrain corrected, using the procedure described by Reference Hastenrath, Rostom and CaukwellMarvin and others (1987), to conform with the terrain model and identifiable terrain features. The 26 control points established on the mountain from a precision ground survey (see above), but not visible on SPOT, were entered in the digital array. The pixels were resampled to a nominal resolution of 10 m.

Results and Outlook

The digital map with date 4 February 1988 created as described above was produced in the form of a photo-map at scale 1 : 5000, so as to be compared against our 3 September 1987 map constructed from aerial photogrammetry (Reference HastenrathHastenrath and others, 1989). A digital computer print at scale 1 : 15 000 is reproduced in Figure 3 for comparison with Figures 1 and 2. The glacier boundaries and the location of the terrain control points on the SPOT photo-map (4 February 1988) and on the air-photogrammetric map (3 September 1988) agreed within the mapping accuracy commensurate with the 1:5000 scale. Various digital enhancements were found possible on the screen of a Macintosh computer. These may bring out the fine structure of crevasse patterns and moraine systems. A distinctly seasonal feature is apparent on the 4 February 1988 SPOT subscene (Fig. 3), namely the brightness of the glaciers on the north side of the mountain, at variance with the 3 September 1987 air photograph (Fig. 2), but consistent with the experience that the glaciers in those quadrants are snow-covered in boreal winter but become bare in summer.

This effort of mapping the glaciers of Mount Kenya from satellite faced as major difficulties the extreme local relief, the reference to a local rather than absolute coordinate system, and the demanded fine scale of 1:5000. The results are encouraging. Given the prospect of an intensifying greenhouse effect, tropical glaciers may continue to retreat at an accelerated rate (Reference HastenrathHastenrath and Kruss, 1992a, b). The work completed will provide an internally consistent reference for determining the changes of Mount Kenya’s glaciers from February 1988 to a SPOT sampling intended for the mid-1990s, which may further be complemented by a remapping from purpose-flown aerial photography. The present pilot study for Mount Kenya also indicates the feasibility of detecting, with adequate accuracy, changes of the relatively small glaciers on other tropical mountains.

Acknowledgements

Fig .3. Digital map of Mount Kenya’s glaciers produced from multi-channel SPOT subscene of 4 February 1988. Points identified in the map Figure 1 are shown by small dots. Scale is 1:15000.

This study was supported by U.S. National Science Foundation Grants EAR-9016657 and EAR-9217211. Processing of the digital satellite information was performed at STX Remote Sensing Services, Lanham, MD. The glacier research is authorized by the Office of the President, Republic of Kenya.

References

Allison, I. Peterson, J.A.. 1989 Glaciers of Irian Jaya, Indonesia. U.S. Geol. Sun. Prof. Pap. 1386H, 1120.Google Scholar
Charnley, F.E.. 1959 Some observations on the glaciers of Mt. Kenya. J. Glaciol., 3(26), 483492.CrossRefGoogle Scholar
Forschungsunternehmen Nepal-Himalaya.1967 Mount Kenya 1 : 5000 Wien, Kartographische Anstalt Freytag–Berndt und Artaria. Google Scholar
Hastenrath, S. 1984 The glaciers of equatorial East Africa. Dordrecht, etc., D.Reidel Publishing Co. CrossRefGoogle Scholar
Hastenrath, S. Kruss, P.D.. 1992a. Greenhouse indicators in Kenya. Nature, 355(6360), 503504.CrossRefGoogle Scholar
Hastenrath, S. Kruss, P.D.. 1992b. The dramatic retreat of Mount Kenya’s glaciers between 1963 and 1987: greenhouse forcing. Ann. Glaciol., 16, 127133.CrossRefGoogle Scholar
Hastenrath, S. Rostom, R. Caukwell, R.A.. 1989 Variations of Mount Kenya’s glaciers, 1963–87. Erdkunde, 43(3), 202210.Google Scholar
Marvin, J.W. Labovitz, M.L. Wolfe, R.. 1987 Derivation of a fast algorithm to account for distortions due to terrain in earth–viewing satellite sensor images. IEEE Trans. Geosci. Remote Sensing, GE– 25(2), 244251.CrossRefGoogle Scholar
World Glacier Monitoring Service In press. Fluctuations of glaciers 198690 (Vol.VI). Zurich, IAHS–UNEP–UNESCO.Google Scholar
Young, J.A.T., Hastenrath, S.. 1991 Glaciers of Africa. U.S. Geol. Surv. Prof. Pap. 1386G, 4970.Google Scholar
Figure 0

Fig .1. Map of Mount Kenya’s glaciers on 3 September 1987. The key below lists for the north and south sides of the mountain and from west to east the names of the glaciers indicated in the map by number. Triangles indicate terrain control points used in the 1987 mapping, and open circles other control points surveyed and marked in the terrain, with letter code referring to Hastenrath and others (1989). Crosses represent terrain features identified for geocoding of digital SPOT scene. Square dots mark the highest peaks (Batian 5199 m, and Nelion 5188 m). Height contours are at 200 m spacing and scale is 1:15 000 (source: Hastenrath and others, 1989). Dot raster denotes glaciers and dash raster lakes. North side: 13 Cesar, 14 Joseph, 16 Northey, 1 Krapf, 2 Gregory; south side: 10 Tyndall, 9 Heim, 8 Forel, 7 Diamond, 6 Darwin, 4 Lewis.

Figure 1

Fig .2. Air photograph of Mount Kenya’s glaciers on 3 September 1987, of flight (PHOTOMAP Kenya Ltd.) from which map was produced by aerial photogrammetry (Hastenrath and others, 1989). Approximate scale 1: 15000.

Figure 2

Fig .3. Digital map of Mount Kenya’s glaciers produced from multi-channel SPOT subscene of 4 February 1988. Points identified in the map Figure 1 are shown by small dots. Scale is 1:15000.