Hydrologie Applications of Space Technology (Proceedings of the Cocoa Beach Workshop, Florida, August 1985). IAHS Publ. no. 160, 1986.

Snow cover mapping for runoff simulations based on Landsat-MSS data in an aipine basin

M. F, BAUMGARTNER & K, SEIDEL Institute for Communication Technology, Swiss Federal Institute of Technology, Zurich, Switzerland H, HAEFNER S K. I, ITTEN Remote Sensing Laboratories, Department of Geography, University of Zurich, Switzerland J. MARTINEC Federal Institute for Snow S Avalanche Research, Weissfluhjoch-, Switzerland

Abstract For the Martinec-Rango Snowmelt Runoff Model (SRM) a precise determination of the daily changes of the snow cover is the basic input variable. A method was developed for digital snow- cover classification with Landsat-MSS data including an extra­ polation for partly clouded scenes and areas with not available remote sensing data. The runoff simulations were tested in a large basin of the Swiss Alps (-Felsberg, 3249 km2, 571- 3614 m a.s.l.) regarding the heterogeneous topographic and cli­ matic conditions.

Introduction and Objectives The assessment and forecast of the quantities of surface runoff is a key objective for water management. In mountain areas the runoff in spring and summer is primarily affected by the changes of the seasonal snow cover. Various tests have demonstrated that the Martinec-Rango Snow- melt Runoff Model (SRM; Martinec, 1980) is suited to forecast surface runoff under mountaineous conditions (Martinec, 1970 and 1975). The accuracy of the forecast depends especially on a measurement of the daily changes of the snow cover. An up-to- date, fast and exact recording of these changes is only possible by remote sensing techniques (Haefner, 1980; Itten, 1980). Therefore, it is the purpose of this paper to outline a me­ thod to map and determine the snow cover in a large basin by Landsat-MSS data. The satellite images have been processed digitally using an interactive image processing system (Besse, 1981; Egeli, 1983). Special attention has to be given to the difficult and var­ iable topographical and climatical conditions in the Alps. Ser­ ious problems were the often occuring clouds and not recorded data, which complicate the determination of the snow cover. Even less suited images have to be considered, in order to get enough remote sensing information to monitor the varying snow cover during a melting period. The areas with 'missing data' have to be estimated using an extrapolation technique, based on 191 192 M.F.Baumgartner et al.

a digital terrain model. The different climatological conditions in the basin and the frequent changes of rain- and snowfall during the melting period complicate the application of the SRM and have to be respected for discharge simulation. A solution is presented to overcome these obstacles for a rather large catchment area in the Swiss Alps and results are presented.

Study Site

The test site 'Felsberg' is located in the eastern part of the Swiss Alps (Fig. 1) and includes the four Rhine-tributaries 'Vorderrhein', '', 'Landwasser' and ''. Its size is 3249 km , the basin relief 3043 m (lowest point: stream gauge 'Felsberg', 571 m a.s.l., highest point: 'Toedi', 3614 m a.s.l.). The ice covered area (glaciers) totals 3% at the end of the hydrological year (September 30th).

FIG.l Location of the 'Felsberg' basin: area of Switzerland: 41,293 km2 area of 'Felsberg' basin (shaded area).- 3,249 km

The basin can be divided into three different climatological regions. The western part ('Vorderrhein') has a significant higher amount of precipitation than the eastern part ('Landwas­ ser' and 'Albula'). The southern part ('Hinterrhein') is influ­ enced by the climate of the South Alps, i.e. very high precipi­ tations for a short time during summer and autumn, compared to the more evenly distributed precipitations in the northern part of the Alps. This differences can be explained by diffe­ rently oriented mountain chains and the influence of the wester­ ly winds. The water regime in spring and summer depends mainly on snowmelt, in summer and autumn on rainfall. In many valleys of the Alps the runoff is controlled by reser­ voirs for hydroelectric power generation. This raises a serious problem in the comparison of the measured with the calculated runoff on a daily basis. The recorded (actual) daily minima and maxima and weekly variations reflect the needs for electricity Snow cover mapping for runoff simulations 193 and does not show the natural runoff.

Snow Cover Determination by Satellite Data

For the determination of the changing snow cover during the snowmelt period Landsat-MSS digital recordings were used. Im­ portant procedures for the processing of the satellite images are :

- Preprocessing, as a geometric transformation of the images to a basic coordinate system and resampling to a unique pixelsize - Determination of an adequate classification scheme and training samples for a supervised classification - Snow cover classification - Superimposition of all data into a multivariate data set - Extrapolation to estimate the snow cover for areas with not available data

- Problem Spots of Snow Cover Determination

Due to the increase of the solar irradiance in spring and summer clouds develop in the late morning (Rott', 1978). The Landsat-MSS sensors (Tab. 1) are recording the radiation reflected from the clouds and therefore the earth surface is invisible. In addition the spectral bands do not allow an automated separation between snow and clouds. In a false color representation a different-

Tab. 1: Landsat-MSS characteristics Landsat 1-3 Landsat 4,5

Altitude 920 km 705 km Period 102 min 9 9 min Frame area 185*185 km2 185*185 km2 Resolution 56*79 m2 56*79 m2 Repetition rate 18 days 16 days Band 4 0.5-0.6 um 0.5-0.6 um Band 5 0.6-0.7 um 0.6-0.7 um Band 6 0.7-0.8 um 0.7-0.8 um Band 7 0.8-1.1 um 0.8-1.1 um

iation becomes possible, based on the cloud forms and the cloud shadows. For every image a digital cloud mask was produced by visual inspections on a interactive image processing system. Together with missing frame segments (not available remote sen-• sing data) a special procedure was established: an extrapolation technique, based on elevation, aspect and slope for each pixel, allows an estimation of the snow coverage in the unknown parts (Seidel, 1983). This and the request for the snow coverage in different elevation zones lead to the use of a Digital Terrain Model (DTM) as a suitable auxiliary information source for the processing of satellite data. Because of the repetition rate of Landsat (Tab. 1) and the often by heavyly cloud coverage obscured scenes, a sufficient, 194 M.F.Baumgartner et al.

number of images during a snowmelt season is very seldom re­ corded. Considering even less suitable scenes, for the snowmelt period 1982 sufficient data were available. The following dates with frame numbers 209/027 and 209/028 have been used:

(1) 2 5-MAR82 1.4% cloudcover/18.8% not available (2) 18-MAY82 15.1% cloudcover/13.1% not available (3) 5-JUN82 20.9% cloudcover/ 2.1% not available (4) 11-JUL82 not available (5) 29-JUL82 10.9/0.20.9% cloudcover/ 3.4% not available It was necessary to compensate in (1) for missing lines by doub­ ling the neighboring lines and in (3) for misregistered lines and channels. The frame (209/027) of scenes (4) and (5) did only cover a part of the basin. They were digitally mosaiced with the southwards adjacent frames (209/028). The cloud coverage on the record of July 29th has been re­ duced by a pixelwise comparison with the record of July 11th, i.e. in each doubtful case the picture element in (5) has been assigned to the category of the corresponding element in (4). This procedure is not possible as soon as new snowfall has changed the situation in the mean time.

Digital Terrain Model and Basin Perimeter

The DTM is used for the extrapolation process mentioned above For the test site 'Felsberg' the DTM was interpolated from the (250*250) rn grid into a raster size of (100*100) m applying the Akima method (Akima, 1974), which is based on a bicubic in terpolation. In order to determine the contribution to the total snow coverage from 'obscured' and 'not available' areas we used a method published earlier (Seidel, 1983). We devided the conti­ nuum of elevation-aspect-slope values into 205 classes:

- 5 elevation zones (560-1100, 1100-1600, 1600-2100, 2100-2600, 2600-3600 m a.s.l.) - 8 aspect classes (N, NE, E, SE, S, SW, W, NW, N) - 5 slope classes (<0-15, 16-22, 23-28, 29-36, 37-90c - 1 class for flat terrain The basin boundary has been digitized from the National Topo­ graphic Map 1:50'000.

Preprocessing The Landsat-MSS digital products were received in UTM coordi­ nates. They were transformed to the coordinate system of the National Topographic Map (based on a conform, Special-UTM pro­ jection) to facilitate the comparison of the different data sets (satellite images, DTM, masks). With the standard reference point approach, the Landsat images and the map were registered by an affine backward transformation (Steiner, 1977). At the same time a nearest neighbor resampling from 56*79 m MSS reso­ lution to the common raster size of (100*100) m was carried out. Snow cover mapping for runoff simulations 195 - Classification

The different Landsat-MSS scenes have been analysed with re­ spect to snow cover. During a supervised classification process each picture element was assigned to one of the following cate­ gories: 'snow free1, 'transition zone' and 'snow covered1. For the different scenes we had to use different classification schemes due to recording conditions: illumination, atmospheric conditions, areal extent of snow, reflectivity of the snow sur­ face. The accuracy of the classification increases with the number of classes. For the most complex situation we had to consider a detailed list of classes:

- Water - Forest (conifereous) - Crop fields - Meadow - Pasture - Wet grassland (wet and long grass from the previous year, beneath the transition zone, watered by the melting snow) - Transition zone (snow patches, approx. 50 % snow covered) - Dark snow (less reflecting snow due to soiling or debris) - Snow (in shadow) - Snow (in sun) For the supervised classification, training samples for each of these classes have been selected. Experience and familiarity with the terrain is necessary to delineate representative samples. Special attention is recommended during the time of maximum snowmelt. Uncarefully selected training samples cause serious misclassifications. In parallel to the complexity of the set of training samples an appropriate expensive classification procedure was found to be adequate ranging from Parallelepiped Discriminant (PPD), Minimum Distance, Mahalanobis Distance and Maximum Likelihood analyses. In addition, a reduction of the variables using Prin­ cipal Components as well as Fisher's Linear Discriminant Vec­ tors (Fisher, 1936) was carried out. Best classification re­ sults were always obtained using two MSS channels (band 5/band 7). Finally in the classified images, the detailed classes were condensed into the above mentioned main categories: 'snow free', 'transition zone' and 'snow covered'. A pictorial representation is given in Fig. 2. Until now, the areas marked as 'clouds' or 'not available' were dis­ regarded and have to be related to the three main categories with a method outlined in the following.

- Snow Cover Extrapolation The classified Landsat-MSS images and the DTM classes are super­ imposed into a multivariate data set. As mentioned above, the region can be subdivided into 205 classes using a DTM. Each class is characterised by defined ranges of elevation, aspect and slope. Under the assumption that for each class a unique relative snow coverage is given we assigned to all cloud covered picture elements the same relative snow coverage as detected for the same class in the cloudfree portion. Finally, the snow 196 M.F.Baumgartner et al.

LftNOSAT-MSS' 18. MftV 1982 SNOWCOVER CLASSIFICftTIOH

FIG.2 Classification result of a partly cloud covered scene.

coverage is summed up for each elevation zone separately, as demanded by the SRM. The contribution of the transition zone was weighted with 50 %. From these values the depletion curves have been constructed (Fig. 3), representing the relative snow coverage for each elevation zone during the snowmelt period in 1982.

RPRIL - SEPTEMBER 1982 FIG.3 Depletion curves for five elevation zones A to E.

The Martinec-Rango Snowmelt Runoff Model

The SRM-Model is designed to simulate or forecast the daily discharge in mountain basins, resulting mainly from snowmelt but also from precipitation (Martinec, 1983). Climatological variables, such as daily air temperature, pre­ cipitation and snow coverage are required. The values were taken from the records of one or more base stations. The discharge was measured at the stream gauge 'Felsberg' and has been used to Snow cover mapping for runoff simulations 197 estimate the overall accuracy of the simulated runoff. The dis­ charge is heavily influenced by the hydroelectric power gene­ ration so that the daily, weekly or seasonal runoff does not show the natural discharge. In addition, basin dependent parameters are required to specify the global hydrological and topographical properties. With reference to the 'Dischma' Valley (Martinec, 1970) the above mentioned five elevation zones were selected. For the test site 'Felsberg' the recession coefficient was derived from the known coefficient of the 'Dischma' Valley (Martinec, 1963). Based on the snow cover measurements a snowmelt runoff simu­ lation with the SRM was carried out for the period April 1st to September 30th, 1982.

Simulations and Results

Daily air temperature and precipitation are necessary input variables for the SRM. In order to get representative values for the rather large 'Felsberg'-basin we used the arithmetic means of 7 base stations located within the catchment area: Alveneu (1175 m a.s.l.), Arosa (1847 m a.s.l.), Chur (586 m a.s.l.), Davos (1590 m a.s.l.), Disentis (1180 m a.s.l.), Hinterrhein (1619 m a.s.l.), Weissfluhjoch (2667 m a.s.l.). We found the resulting discharge curve from this 'synthetic base station' very similar to the curve based on the 'Disentis' temperature and precipitation records. That will say for fur­ ther computations the 'Disentis' values are representative for the 'Felsberg'-basin. In Fig. 4 the simulated and measured discharge curves are presented for comparison.

o-| 1 , 1 1 ! 1 , , , 1 , 1 1 , 1 , 1 1 H 0 50 100 150 APRIL - SEPTEMBER 1982 FIG.4 Simulated (C) and measured (M) discharge: runoff M: 2778 x 10^ m* runoff C: 2937 x JO6 m3 relative difference: 5.7%

In view of the operation of reservoirs the daily simulated flows do not always agree with the measured discharge. The simulated seasonal runoff volume, however, does not deviate very much from the measured volume. It has to be emphasized that during the simulation with the 198 M.F.Baumgartner et al.

SRM no 'updating' was carried out, i.e. the measured discharge was at no time used to correct the deviation.

Conclusions With the applied procedures it is possible to determine the changing snow cover with Landsat-MSS data even in a large al­ pine basin. Furthermore methods have been developed to esti­ mate the snow coverage for partly clouded areas. The runoff can be simulated with the SRM if enough suitable satellite images during the snowmelt period are available. To this effect it is necessary to test a possible data complementa­ tion by NOAA-AVHRR imagery.

Acknowledgement Many thanks to Dr. A. Rango and E. Major for making available the SRM computer program for this study.

References Akima, H., 1974, A Method of Bivariate Interpolation and Smooth Surface Fitting Based on Local Procedures, Commun. Ass. Comput. Mach., Vol. 17, No. 1, pp. 18-20, 26-31. Besse, L., Seidel, K., Kuebler, 0., 1981, Large Scale Multipur­ pose Interactive Image Processing Facility at ETH Zuerich, Design and Digital Image Processing Systems, James L. Mannos, Editor, Proc. SPIE 301, pp. 154-161. Egeli, E., Gerig, G., Klein, F., Kuebler,0., 1983, A Hardware and Software Optimized Program System for Interactive Image Processing, Architecture and Algorithms for Digital Image Processing, Andre Oosterlinck, Per-Eric Danielsson, Editors, Proc. SPIE 435, pp. 134-138. Fisher, R.A., 1936, The Use of Multiple Measurements in Taxo- nomic Problems, Annals of Eugenics 7, Part II, pp. 179-188. Haefner, H., 1980, Snow surveys from Earth Resources Satellites in the Swiss Alps - A Review on Six Years' Research, Remote Sensing Series, Department of Geography, University of Zurich, pp. 1-65. Itten, K.I., 1980, Snow Mapping from Spaceplatforms, in: SAL0- M0NS0N V.V. + BHAVSAR P.D. (Edit.): The Contribution of Space Observations to Water Resources Management, COSPAR Symp., Bangalore, India. Martinec, J., 1963, Seasonal Forecasts of Discharge for Operation of Reservoirs (In Czech, English Abstract), Prace a Studie VUV Praha-Podbaba, Vol. 110, p. 27. Martinec, J., 1970, Study of Snowmelt Runoff Process in two Representative Watersheds with Different Elevation Range, Symp. on the Results of Research on Representative and Experimental Basins, Wellington, IAHS Publ. No. 96, pp. 29-39. Martinec, J., 1975, New Methods in Snowmelt Runoff Studies in Representative Basins, In: The Hydrological Characteristics of River Basins, Proc. of the Tokyo Symp., IAHS Publ. No. 117, pp. 99-107. Snow cover mapping for runoff simulations 199

Martinec, J., 1980, Snowmelt-Runoff Forecasts Based on Auto­ matic Temperature Measurements, Proceedings at the Oxford Symposium, IAHS Publ. No. 129, pp. 239-246. Martinec, J., Rango, A., Major, E., 1983, The Snowmelt Runoff Model (SRM) User's Manual, NASA Reference Publication 1100, Washington, D.C., 110 pp. Rott, H., 1978, Zur Schneekartierung in alpinen Einzugsgebie- ten aus Satellitenbildern, Zeitschrift fuer Gletscherkunde und Glazioloqie, Bd. 14, H. 1, pp. 81-93."" Seidel, K., Ade, F., Lichtenegger, J., 1983, Augmenting LAND- SAT MSS Data with Topographic Information for Enhanced Registration and Classification, IEEE Transactions on Geo- science and Remote Sensing, Vol. GE-21, No. 3, pp. 252-258. Steiner, D., Kirby, M.E., 1977, Geometrical Referencing of Landsat Images by Affine Transformation and Overlaying of Map Data, Photogrammetria, Vol. 33, pp. 41-7 5.