Introduction
It is estimated that the mean surface air temperature of the Earth will rise at a rate in the order of 0.3° C per decade during the next decades due to increasing concentrations of CO2 and other trace gases in the atmosphere (Reference Houghton, Jenkins and EphraumsHoughton and others, 1990, Reference Houghton, Callander and Varney1992). This rate of warming will have pronounced effects on glaciers and ice caps and lead to major run-off changes in glacierized areas. The Nordic research project “Climate change and energy production” (Reference SælthunSælthun, 1992) was started in 1991 with the aim of estimating the hydrological effects of global warming in the Nordic countries with special emphasis on the possible consequences for the operation and planning of hydroelectric power plants.
As a part of this project, a degree-day glacier mass-balance model for temperate glaciers has been developed in order to estimate the effects of global warming on glacier mass balance and run-off from glacierized areas. This model has been tested on the glaciers Sätujökull (part of Hofsjökull), central Iceland, Nigardsbreen part of Jostedalsbreen. southern Norway, and Qamauârssûp sermia (outlet glacier from the Greenland ice sheet), West Greenland. In this paper, we describe the application of the mass-balance model to these glaciers and compute the predicted effect of climatic warming on the mass balance and on the run-off from these glaciers. The stability of the model parameters over the period of measurements is analyzed in order to assess whether these parameters are likely to change as a consequence of the climate changes which are considered.
Mass-balance data from glaciers and ice sheets contain implicit information about the dependence of glacier mass balance on climate. The measured mass balance varies with elevation and from year to year, mainly as a consequence of variations in temperature and precipitation. It should, therefore, be possible to use measured variations in the mass balance together with meteorological data to parameterize the relation between glacier mass balance and climate using a suitable glacier mass-balance model. The mass-balance model may then be used to estimate likely glacier mass-balance changes resulting from hypothetical climatic changes.
The mass-balance measurements from the three glaciers come from different climate regimes and should therefore be a good test of the applicability of a degree-day glacier mass-balance model over a wide range of conditions. Sálujökull is an outlet glacier from an ice cap in the middle of the Icelandic highland; Nigardsbreen is a narrow valley glacier with a wide elevation range; and the measurements from Qamanârssûp sermia come from the ablation area of the Greenland ice sheet. The measurements from Nigardsbreen cover 24 years, which can be divided into sub-sets with significantly different climate Conditions, and are therefore especially suitable for an analysis of the stability of model parameters. Detailed energy-balance and degree-day glacier mass-balance modelling has previously been applied to mass balance and meteorological data from Qamanärssșp sermia, and our results can therefore be compared to previous results from tins glacier (cf. Reference Braithwaite, Olesen and OerlemansBraithwaite and Olesen, 1989, Reference Braithwaite and Olesen1990a, Reference Braithwaite and Olesen1993).
The mass-balance modelling presented here is a continuation of the modelling described by Reference Laumann and NLaumann and Reeh (1993) which applied a similar degree-day glacier mass-balance model to three glaciers in Norway. The model used here has been improved to accommodate daily and monthly input data instead of yearly data, and the calibration of the model is based on yearly mass-balance measurements rather than long-term mass-balance averages. We place special emphasis on the stability of the model parameters by analyzing how well the model explains yearly mass-balance variations and we test the applicability of the degree-day model by applying it to three glaciers in very different climate regimes.
Model Description
The glacier mass-balance model, which is termed MBT (an acronym for Mass Balance of Temperate glaciers), is based on existing degree-day hydrological and glacier mass-balance models which have been used extensively in the Nordic countries Reference Jóhannesson and LaumannJóhannesson and Laumann, 1993): the HBV model (Reference BergströmBergström, 1976; Reference SælthunSælthun, 1990), the NAM2 model (Reference GottliebGottlieb, 1980; Reference Holm and EinarssonHolm and Einarsson, 1992), the MB1 model Braithwaite. 1984) and a glacier mass-balance model described by Reference ReehReeh (1991) and Reference Laumann and NLaumann and Reeh 1993).
A degree-day model was chosen for this analysis in preference to an energy-balance model for the following reasons:
-
1. The formulation of the mass balance in terms of temperature and precipitation is advantageous for simulations of the response of glacier mass balance to climatic changes, since climate scenarios are typically given in terms of temperature and precipitation changes. Scenarios for future changes in cloudiness, cloud height, wind speed, relative humidity and other variables which are often required in energy-balance models (Reference Braithwaite and OlesenBraithwaite and Olesen, 1990b; Reference OerlemansOerlemans, 1992) are in general not available.
-
2. The use of energy-balance models to compute glacier mass balance over whole ice caps or glaciers raises difficult questions about the variation of input data and model parameters over the area covered by the ice cap or glacier Reference BraithwaiteBraithwaite 1992). Except for temperature, the input data are in general not readily extrapolated away from measurement locations. Similar problems are, of course, also encountered for degree-day models but they are not as serious as for energy-balance models.
-
3. Energy-balance models do not seem to perform better than degree-day models in practical hydrological applications when modelling daily run-off on catchment scale (Reference BergströmBergström and others, 1992).
-
4. Degree-day coefficients show a high degree of stability with time and across a wide range of different climatic regimes Reference Braithwaite and OlesenBraithwaite and Olesen. 1984, Reference Braithwaite, Olesen and Oerlemans1989; Reference BergströmBergström and others, 1992; Braith-Waite, 1992).
-
5. A degree-day glacier mass-balance model can benefit from the considerable practical experience witli similar glacier (MB1) and hydrological models (HBV, NAM2) which has been built up in Denmark, Greenland, Iceland, Norway and Sweden.
-
6. Experience from a degree-day mass-balance model can be used to improve existing hydrological models in partly glacierized catchments (Reference Einarsson, óhannesson, Kern-Hansen, D and REinarsson and Jöhannesson, 1994).
The mass-balance model determines the precipitation, snow accumulation, melting of snow and ice, and re-freezing as a function of altitude based on observed temperature and precipitation at a meteorological station. The total ablation, mass balance and run-off can then be computed. Precipitation and temperature on the glacier are computed assuming a linear precipitation variation with altitude and a constant temperature lapse rate with altitude.
The model computations are based on daily, monthly or yearly precipitation values. The temperature is Specified as a series of daily or monthly mean values or as a sinusoidal function with a period of 1 year according to the following equation
where t is time since the beginning of the mass-balance year in days), A is the length of the year in (days), T a is the average temperature of the year, and T co and T si are coefficients. A statistical approach is used in the determination of the number of positive degree-days and snow accumulation (see below) when the temperature is given as monthly mean values or as a sinusoidal function. Temperature deviations from the monthly mean or from the sinusoidal function are assumed to be normally distributed with a standard deviation σ (cf. Reference BraithwaiteBraithwaite, 1984; Reference ReehReeh, 1991). Precipitation is assumed to be evenly distributed in time in the statistical computation of snow accumulation from monthly and yearly input values.
Precipitation on the glacier is assumed to fall as snow, if the temperature at the altitude in question is below a specified threshold. This threshold is typically about 1° C. A wind correction with separate correction factors for rain and snow may be applied to the precipitation data before they are used in the mass-balance computations. When the temperature is specified as monthly mean values, the fraction f of the precipitation that falls as snow in a given month is given by
where T is temperature, T s, is the rain/snow threshold temperature and T m is the monthly mean temperature. When the temperature is specified as a sinusoidal function, the fraction of the precipitation that falls as snow in a given period from day D 1 to day D 2 is given by
where D = D 2 − D 1 + 1 is the period length in days and T d(t) given by Equation (1). The fraction f determined from Equations (2) and (3) is equal to the fraction of the period with T < T s.
Melting of snow and ice is computed from the number of positive degree-days (PDD), using different degree-day factors (amount of melting per PDD) for snow and ice (cf. Reference Braithwaite, Olesen and OerlemansBraithwaite and Olesen, 1989; Reference ReehReeh, 1991). When the temperature is specified as monthly mean values, PDD is given by
where the number of days in a year has been distributed evenly among the months for simplicity. When the temperature is specified as a sinusoidal function, PDD is given by
When the snow thickness becomes less than a specified threshold, the degree-day factor is found as a weighted average of the degree-day factors for snow and ice. The reason for this is that the snow-line is not a sharp well-defined line at a certain altitude. Rather, it represents a transition from a surface of clean ice to a surface completely covered by snow, where patches of clean ice and snow will be mixed.
The model computes the amount of stored meltwater or rain (which may be frozen) by assuming that the storage is a certain fraction of the amount of remaining snow. This fraction is specified as a model parameter and is typically between 5% and 10%. This storage leads to a delay in the onset of glacier run-off with respect to the start of melting on the glacier.
The ablation is defined as the negative of the melting (of both snow and ice) plus the refrozen and/or stored liquid water. The mass balance is the sum of the accumulation and the ablation. The run-off from each elevation on the glacier is given as the precipitation minus the mass balance.
The model computes mass balance over a specified time interval which is usually whole or part of a mass-balance year. A more detailed description of the model has been given by Jöhannesson and others (1993).
Mass-Balance Measurements and Model Calibration
The mass-balance model calculates total accumulation, total ablation and therefore the mass balance at specified elevations on the glacier surface between given dates. These quantities are not always directly comparable to traditional glacier mass-balance measurements. In order to calibrate the model, using measured mass balance, it is important to take account of the measurement dates and to consider the exact meaning of the various terms in published mass-balance data (cf. Reference Østrem and BrugmanØstiem and Brugman, 1991).
Until 1975, mass balance on Nigardsbreen was measured frequently during the summer by field assistants, in order to obtain the balance at the end of the stratigraphic balance year, both in the accumulation area and on the glacier tongue. The published figures are therefore derived from measurements made on different dates and at different locations on the glacier. In later years, the same date was used over the whole glacier, although new snow King on the summer surface was not included in the summer balance. In addition, winter balance in the ablation area was usually measured by probing rather than by stake measurements, such that melting of ice during the previous autumn (after the summer-balance measurements) was counted as part of the summer balance instead of the winter balance. Thus, comparison of modelled and measured mass balance for Nigardsbreen is complicated.
The mass balance of Hofsjökull and Qamanärssûp sermia is measured between known calendar dates and it is therefore straightforward to compare the model results with the measured summer, winter or yearly balance.
The purpose of model calibration is to find the set of parameters which leads to the smallest prediction error when the mass balance found by the model is compared with the actual mass-balance measurements. When measurements of both winter and summer balance are used, the model is run in two steps, first from the autumn to the spring and then from the spring to the autumn, in order to compute winter- and summer-balance values which are directly comparable with the measurements. For Nigardsbreen, one should ideally use different dates for different years and different altitudes and adjust the computed balance for late summer snowfall and melting of ice during the autumn after the summer-balance measurements, but there are practical difficulties which prevent us from making this correction (see below).
The model uses a total of 13 parameters when the mass balance is computed from daily temperature and precipitation data (cf. Reference Jóhannesson, Sigurôsson, T and KennettJóhannesson and others. 1993). Some of the parameters are known (i.e. the altitudes of the meteorological stations). Others are given fixed values which are not changed in the calibration. These are the snow/rain threshold temperature, the rain and snow correction factors, the precipitation correction and the starting elevation for the precipitation gradient, the threshold snow thickness for blending of the degree-day factors for snow and ice, and the refreezing ratio. This haves four parameters, i.e. the degree-day factors for snow and ice, and the temperature and precipitation gradients, which have to be determined by calibration.
The snow/rain threshold has a similar value to that used in hydrological models in the Nordic countries. The value of the refreezing ratio is chosen on the basis of unpublished measurements of the liquid-water content and temperature of snow in Norway. It is similar to the hydrologically determined ratio which is used in the HBV model in Sweden and Norway. The values of the fixed model parameters for the three glacers are given in the tables in the following section.
The model parameters which are considered in the calibration are determined by running the model for all years within the calibration period, and finding the parameter values which minimize the sum of squares of the residuals, i.e. the differences between the model output and the measurements for every year. The minimization is carried out iteratively using a non-linear least-squares routine (nls) which is a part of the statistical package S (Reference Chambers and HastieChambers and Hastie, 1992). This calibration procedure is somewhat different from the method used by Reference Laumann and NLaumann and Reeh (1993) for Nigardsbreen. The latter is based on a comparison of the average computed accumulation and ablation with average measured winter and summer balances. However, in the method described here, the mass balance which corresponds to the mass-balance measurements (winter, summer or year) is computed directly. Furthermore, the year-to-year variability of the measured mass balance is explicitly considered when the best model parameters are determined.
The computed summer balance is dependent on the degree-day factors and on the temperature gradient, but they are relatively insensitive to the precipitation gradient. The computed Winter balance, on the other hand, is most sensitive to the value of the precipitation gradient. The degree-day factors and the temperature gradient for Nigardsbreen and Sátujükull are found by minimizing the residuals in the summer balance after a suitable value of the precipitation gradient has been found. The precipitation gradient is determined by trial and error, primarily on the basis of the winter-balance measurements. This is done in order to prevent undesirable effects from errors in the predicted winter balance on the degree-day factors and the temperature gradient. For Qamanärssûp sermia, yearly mass-balance measurements are used for simplicity and the precipitation gradient is set to zero based on scattered winter-balance measurements personal communication from R.J. Braithwaite ).
Results
Sätujökull, central Iceland
Mass-balance measurements have been carried out on Sätujökull by Orkustofnun since 1988 (Reference SigurðssonSiguròsson, 1989, Reference Sigurðsson1991, Reference Sigurðsson1993). The model compulations are based on temperature and precipitation measurements at Hveravellir (64°52′N, 19°34°W: 641 m a.s.l., about 33km from the glacier). Daily meteorological measurements at Hveravellir are available from the Icelandic Meteorological Institute from 1965. Data from the five mass-balance years 1987-88 to 1991-92 are used in the analysis below. The spring and autumn mass-balance measurements on Sätujökull are carried out between 10 and 20 May and between 10 and 30 September each year. For simplicity, the starting points of the winter and summer seasons in the compulations are put equal to 1 October and 1 May for all years.
The fixed model parameters for Sätujökull are given in Table 1.
Fixed model parameters for Sätujökull
The rain- and snow-correction factors rko and sko are recommended values for the meteorological station at Hveravellir (Sigurósson, 1990). The starting elevation for the precipitation gradient computations elq, is chosen equal to the lowest elevation on the glacier. The precipitation correction pko is used to compute the precipitation at this elevation from the corrected precipitation at Hveravellir. The value of pko was chosen by comparison of precipitation data from Hveravellir and measured winter balance at the lowest elevations on the glacier.
The calibration of the degree-day factors and the temperature gradient was carried out first for each of the 5 years, individually. The parameter estimates showed considerable scatter. A special case was the estimated degree-day factor for snow for the mass-balance year 1990-91, which was 0.01 m w.e. °−C−1 larger than the degree-day factor for ice, 0.007 m w.e. °−C−1 and almost twice as large as the estimated degree-day factor for snow for the other years. The glacier was covered by black ash from an eruption from Mount Hekla during this year. The ash greatly increased the melting of snow, because of the reduced albedo of the snow surface. The melting of the ice, on the other hand, was not affected to the same extent, because the ash was quickly washed off the smooth ice surface. The data from the mass-balance year 1990-91 were therefore not used in the remainder of the calibration.
The precipitation gradient was determined as 55% per 100 m. The degree-day factors and the temperature gradient were determined by minimizing the residuals from the summer-balance measurements. The best-fit parameters are given in Table 2.
Best-fit parameter values determined from the measured summer and winter balance of Sátujokull
These parameter values explain 95% of the variance in the measured summer balance of Sátujökull at the individual elevations, but only 42% of the variance in the winter-balance measurements. The estimated degree-day factor for snow is consistent with previous estimates of degree-day factors for snow in Iceland (0.004-0.010 m w.e. ° C−1 d−1 cf. Reference BergströmBergstrom and others, 1992), and somewhat higher than in Scandinavia (0.0018-0.005 m w.e. ° C−1 d−1. cf. Reference BergströmBergström and others, 1992). The estimated degreeday factor for ice is similar to values determined for glaciers in West and northeast Greenland (0.0077 and 0.0094 m w.e. ° C−1 d−1, cf. Reference Braithwaite, Olesen and OerlemansBraithwaite and Olesen, 1989; øeggild and others, 1994).
Figure 1a - e shows measured and computed mass balance for the five mass-balance years. The parameters in Tables 1 and 2 are used for all the years. Figure 1f) shows the average of the mass balance measurements for the four mass-balance years 1987-88, 1988-89, 1989-90 and 1991-92, together with the computed mass balance using the mean yearly precipitation and a sinusoidal representation of the temperature record at Hveravellir for the same years. The standard deviation of the temperature deviation from the sinusoidal temperature approximation was chosen as σ = 3,.32° C, based on the standard deviation of measured temperature deviations from the sinusoidal function during, the melting season (May-September). The figures show that the variation of the summer balance with elevation and between years is in general well explained (except for 1990-91). Year-to-year fluctuations in the winter balance are, on the other hand, not well explained by the model and the variation of the winter balance with elevation is not well predicted. This is probably due to the large spatial variability in precipitation. The minimum in rhe winter balance around 1000 m a.s.l. is caused by wind transport of snow which cannot be explained by the linear increase of precipitation with elevation used in the model. The parameter set of Table 2 leads to an underestimate in the summer melting for the year 1990-91 when the glacier was covered by ash. A better fit was obtained using an increased degree-day factor for snow (dds = 0.009 m w.e.° C−1 d−1) and a higher value of the precipitation gradient (grp= 99% per 100 m) for this year and the results are shown as dashed curves in Figure 1d). The average results for the 4 years 1987-88. 1988-89, 1989-90 and 1991-92 Fig. 1f) are better than for the individual years, for both the summer and the winter balance.
Measured (symbols) and computed (solid curues) mass balance (winter (Δ), summer (Ο) and year (+ )) at Sátujökul for the mass-balance years 1987-88 to 1991 92, and for the average of the years 1987-88, 1988-89, 1989-90 and 1991-92. Model parameters are the same for all the years (cf. Tables 1 and 2). The dashed curves for 1990-91 are produced with a different parameter set (ddi = 0.0077 in w.e. ° C d−1 dds = 0.009 m w.e. ° C d−1, grt = 0.53 ° C per 100m and grp = 90% per 100m).
Figure 2a) shows a scatter diagram of the measured and computed mass balance of Sátujökull for the mass-balance years 1987-88, 1988-89, 1989-90 and 1991-92. Figure 2)a shows the variation of the mass balance with elevation and from year to year on the glacier. The summer-balance values show a large spreading and are well predicted by the model (model explains 95% of variance). The winter-balance values show much less spreading, although similar deviations from the diagonal, and the model performance is therefore worse (model explains 42% of variance) than for the summer balance. Figure 2b) shows the mass balance averaged over all elevations, in order to assess the ability of the model to explain year-to-year mass-balance variations on the basis of year-to-year temperature and precipitation fluctuations. It is seen that the summer-balance variations are rather well predicted by the model (model explains 79% of variance). Year-to-year variations in the winter balance are much smaller and are not well predicted by the model. It should be borne in mind that the above results for Sátujökull are based on measurements from only four mass-balance years.
Measured and computed mass balance of Sátujökull for the mass-balanceyears 1987-88, 1988-89, 1989-90 and 1991-92. (a) Winter (Δ), summer (Ο) and yearly (+) mass balance of individual years at each elevation. (b) Elevation-averaged values for each year.
Nigardsbreen, southern Norway
Mass-balance measurements have been carried out on Nigardsbreen by NVE since 1962 (Norges Vassdrags- og Energivork (NVE), 1964-91). The model computations are based on temperature and precipitation measurements at Bjørkehaug (61°40′N, 07°18′E, 324m a.s.l.) about 5 km from the glacier. Daily meteorological measurements at Bjørkehaug are available from the Norwegian Meteorological Institute (DNMI) for all years since the mass-balance measurements were started, with the exception of parts of the mass-balance years 1981-82 and 1982-83. Data from the mass-balance years 1964-65 to 1980-81 and 1983-84 to 1989-90 (a total of 24 years) are used in the analysis below. Previously published values from below 550 m a.s.l. were not used, as they are based on an extrapolation below the lowest measurements.
The mass-balance measurements were made on different dates on different parts of the glacier in many years. Accurate information about the measurement dates is lacking in many of the mass-balance reports from NVE. It has not been possible to account for this problem with the measurement dates. We have therefore defined the starting points of the winter and summer seasons as 10 September and 15 May for all years. These dates are less than 2-3 weeks in error for most years, but can be up to 2 months in error in some years. We have experimented with different dates for the years where accurate measurement dates are available and have found that the model output is significantly improved in some cases.
The fixed-model parameters for Nigardsbreen are given in Table 3.
Scattered temperature measurements carried out on the tongue of Nigardsbreen at the same elevation as the Bjørkehaug meteorological station, give between 1° and 2° C lower temperatures than at Bøerkehaug. This is taken into account by lowering the elevation of the temperature station in the computations by 200 m. The rain- and snow-correction factors rko and sko are similar to the default values in the Norwegian version of the HBV model.
Fixed-model parameters for Nigardsbreen
The precipitation gradient was found by using the average winter balance for all 24 years and determined to he 9% per 100 m. The degree-day factors and the temperature gradient were determined by minimizing the residuals from the summer-balance measurements. The best-fit parameters are given in Table 4.
Best-fit parameter values determined from the measured summer and winter balance of Nigardsbreen
These parameter values explain 94% of the variance in the measured summer balance and 67% of the variance in the winter balance at the individual elevations. The estimated degree-day factors are consistent with the results of previous studies (Norges Vassdrags- og Energivork (NVE), 1964-91; Reference Furmyr and TollanFurmyr and Tollan, 1975). The degree-day factors used in Reference Laumann and NLaumann and Reeh (1993) (0.0055 m w.e. ° C−1 d−1 for ice and 0.004 m w.e. ° C−1 d−1 for snow) are slightly lower than found here, because they have not lowered the elevation of Bjørke-haug. The estimated temperature gradient of 0.58° C per 100 m is in good agreement with the gradient estimated temperature measurements on the glacier tongue and at Steinmannen (61°41′N, 07°08′E, 1634 m a.s.l.)
Figure 3)a and b shows the measured and the computed mass balance for two of the 24 mass-balance years using the parameters given in Tables 3 and 4 for both years. Figure 3a) shows one of the six best years (1965-66) and Figure 3b) shows one of the six worst years (1980-81). Figure 3c) shows the average of the massbalance measurements for all 24 mass-balance years together with the computed mass balance using the mean yearly precipitation and a sinusoidal representation of the temperature record at Bjorkehaug for the same years. The standard deviation of the temperature deviations from the sinusoidal temperature approximation was chosen to be σ = 3.12° C, based on the standard deviation of measured temperature deviations from the sinusoidal function during the melting season (May-September). This figure shows that the variation of the summer and winter balance with elevation is in general well explained.
Measured (symbols) and computed (solid curves) mass balance (winter (Δ) summer (Ο) end year (+ )) at Nigardsbreen for the mass-balnnce years 1965-66 and 1980-81 and for the average of the years 1964-65 to 1980-81 and 1983-84 to 1989-90. Model parameters are the same for all the years (cf. Tables 3 and 4).
Figure 4) shows scatter diagrams of the measured and computed mass balance of Nigardsbreen similar to Figure 2) for Sátujökull. The quality of the model prediction of the summer balance at the individual elevations in Figure 4)a (model explains 94% of variance) is similar to that for and the predicted winter balance (model explains 67% of variance) is better than for Sátujökull. The elevation-averaged values in Figure 4b) show that the model is able to explain year-to-year summer-balance on the basis of year-to-year temperature and precipitation fluctuations reasonably well (the model explains 62% of variance) but year-to-year winterbalance variations are not so well predicted (the model explains 48% of variance).
Measured and computed mass balance d Nigardsbreen for the mass-balanreyrars 1964-65 to 1980-81 and 1983-84 to 1989-90. (a) winter (Δ) summer (Ο) end year (+) mass balance of individual years at each elevation, (b) Eleuation-averaged values for each year.
Qamanärssûp sermia, West Greenland
The mass-balance data from Qamanärssûp sermia are “centre-line” stakes from the mass-balance years 1979-80 to 1986-87 and span the elevation range 110-1410 m a.s.l. Missing mass-balance data are estimated by a statistical procedure using data from the other stakes (personal communication from R.J. Braithwaite). The mass-balance year extends from 1 September to 31 August. The model computations are based on monthly temperature averages computed from temperature meas-urements at Qamanärssûp sermia Base Camp (64°29′N, 49°29′W, 760 m a.s.l.). Missing temperature values are estimated from temperatures at Nuuk/Godthåb, about 150 km away from Qamanärssûp sermia (Reference Braithwaite and OlesenBraithwaite and Olesen, 1993). The mass-balance and temperature measurements were carried out by the Geological Survey of Greenland (Grønlands geologiske undersøgelse, GGU) (Reference Braithwaite, Olesen and OerlemansBraithwaite and Olesen, 1989; Reference Olesen, Braithwaite and OerlemansOlesen and Braithwaite, 1989).
The mass-balance model is derived for temperate glaciers and ice caps. The Greenland ice sheet is, of course, not temperate. However, the mass-balance data from Qamanärssûp sermia come from the ablation area of the ice sheet, where internal accumulation is not particularly important due to the high net ablation.
The degree-day computations are based on monthly temperature averages assuming temperature deviations from the monthly mean are normally distributed with a standard deviation σ. The standard deviation of the temperature deviations was chosen as σ = 3.5° C. The degree-day computations above for Sátujökull and Nigardsbreen were, however, based on daily temperature data. Reference Jóhannesson, Sigurôsson, T and KennettJóhannesson and others (1993) showed that degree-day sums computed from daily temperature data with the MBT model are consistent with degree-day sums computed from monthly temperature data (cf. Reference BraithwaiteBraithwaite, 1984) and that the degree-day sums are not sensitive to the exact value of σ used.
Accumulation or winter balance data at the stakes on Qamanärssûp sermia is not available. Winter-balance measurements for stakes at 1200 and 1410 m a.s.l., together with measured precipitation at Nuuk/Godthåb and measured summer precipitation at the Qamanärssûp sermia base camp, indicate that the yearly precipitation on the glacier is between 0.5 and 1.0 m w.e., with little gradient with elevation. A fixed value of 0.7 m w.e. for the yearly precipitation with uniform distribution between the months of the year is used in the model computations.
Other values of the precipitation did not significantly affect the calibrated values of the temperature gradient and the degree-day factor for ice. Varying the precipitaction from year to year, based on variation in the measured precipitation at Nuuk/Godthåb, turned out to have little effect on the model calibration. Wind correction for snow or rain is not used and zero precipitation gradient with elevation is specified. The fixed model parameters for Qamanärssûp sermia are given in Table 5.
Fixrd-modelparameter rfijr Qamanärssûp sermia
The mass-balance measurements (Fig. 5) exhibit considerable local deviations from a smooth variation with elevation, mainly due to systematic spatial variations in the accumulation over the ice sheet which can in some cases be explained by the local glacier landscape (personal communication from R.J. Braithwaite). Below about 370 m a.s.l., the measurements show zero or even slightly negative mass-balance gradient with elevation over an elevation range of about 300 m which may to some extent be caused by local advection of cold air from areas adjacent to the lower parts of the glacier (personal communication from R.J. Braithwaite). This feature of the mass-balance profile cannot he adequately modelled by linear precipitation and temperature gradients. Including the mass-balance measurements below 370 m a.s l. in the model calibration leads to a set of unrealistic model parameters which greatly underestimate the massbalance gradient above 370 m a.s.l. (Reference JóhannessonJóhannesson, 1993). Since the overall shape of the mass-balance profile above 370 m a.s.l. is characterized by a relatively monotonic increase of the net balance with elevation (cf. Fig. 5), the calibration procedure was performed for data in the elevation range 370-1410 m a.s.l. The results of this calibration are given in Table 6.
Best-fit parameters determined from the measured yearly net balance in the elevation range 370-1410 m a.s.l. on Qamanärssûp sermia
Measured yearly mass balance (+) and computed yearly accumulation, mass balance and ablation (solid and dashed curves) at Qamanärssûp sermia for the mass-balance years 1980-81 and 1983-84 and for the average of the years 1979-80 lo 1986-87. Model parameters are the same for all the years (cf. Tables 5 and 6). The solid curses correspond to the parameter set in Table 6 derived from mass-balance data in the elevation range 370-1410 m a.s.l. The dashed curves correspondto the best-fit parameter set when points below 370 m a.s.l. are included.
These parameter values explain 84% of the variance in the measured yearly net balance at the individual elevations above 370 m a.s.l. The parameter values in Table 6 are in fairly good agreement with the average degree-day factors estimated by Reference Braithwaite, Olesen and OerlemansBraithwaite and Olesen (1989, Reference Braithwaite and Olesen1993) from mass-balance measurements at stake 751 close to Base Camp, which are 0.0077 m w.e. ° C−1 d−1 and 0.0079 m w.e. ° C−1 d−1, respectively, for ice and 0.0025 m w.e. ° C−1 d−1 for snow.
Figure 5) a and b show the measured yearly mass balance and the computed yearly accumulation, ablation and mass balance for two of the eight mass-balance years using the parameters given in Tables 5 and 6. The model results for 1980-81 (Fig. 5a) are the best for 8 years and the model results for 1983-84 (Fig. 5b) are the worst. Figure 5)c shows the average of the mass-balance measurements for the eight mass-balance years, together with the computed yearly accum-ulation, ablation and mass balance based on the average monthly temperature and precipitation for the same years. (Figure 5c) also shows the model results corresponding to best-fit model parameters when data below 370 m a.s.l are included.)
Figure 6) shows scatter diagrams of the measured and computed yearly mass balance of Qamanärssûp sermia similar to Figures 2 and 4 for Sátujökull and Nigardsbreen, respectively. The quality of the model prediction of the yearly balance in Figure 6a) (the model explains 84% of variance) is slightly worse than shown in Figures 2a and 4a for the summer balance or yearly balance of Sátujökull and Nigardsbreen, respectively. The elevation-averaged values in Figure 6b) show that the model is able to explain year-to-year variations in the mass balance on the basis of year-to-year temperature fluctuations reasonably well (the model explains 69% of variance).
Parameter Stability
An important question with regard to modelling of climate changes is whether model parameters determined from the current climate can be used to predict changes associated with a different climate. The glacier mass-balance model used here reproduces large variations in the mass balance with elevation and from year to year using the same parameter set. The temperature change over the altitude range of the glaciers is much greater than the expected CO2-induced temperature change during the next century, and year-to-year variations in the temperature are of a magnitude similar to the expected climatic warming during the next 50-100 years. Furthermore, degree-day factors for ice determined from measurements in Iceland, Norway and Greenland show relatively little spreading across a wide range of different climatic regimes. A possible climatically induced change in the degree-day factors may therefore he expected to have only a small effect on the ablation compared to the large change in the ablation caused directly by the temperature change.
Measured and computed yearly mass balance of Qamanärssûp sermia above 370 m a.s.l. for the mass-balance years 1979-80 to 1986-87. (a) Yearly mass balance (+) of individual years at each elevation, (b) Elevalion-averaged values for each year.
Another related question is whether degree-day models are suitable for estimating mass-balance changes due to climatic warming at all, because they relate all snow and ice melting to temperature through the degree-day sum. The largest energy-balance component, however, is typically net radiation, which is often not well correlated with temperature. Thus, if the ablation is modelled as linearly dependent on the degree-day sum, one may expect the model to predict too high sensitivity to temperature changes, because in reality only a part of the energy available for melting is directly related to temperature. However, experience shows that ablation and temperature (or rather the degree-day sum) are well correlated, whereas there is a rather poor correlation between ablation and net radiation (Reference BraithwaiteBraithwaite, 1981). This question may be addressed on the basis of Figures 2b, 4b and 6b which show that the model reproduces 60-80% of the measured year-to-year variance in the elevation-averaged summer balance or net balance. In addition, the slope through the summer-balance or net-balance points in the figures is close to 1, which indicates that sensitivity of the model to temperature variations is not in fact too high. If the sensitivity of the model had been significantly overestimated, then the summer-balance and net-balance points in Figures 2b, 4b and 6b would have had a trend with a slope greater than 1.
The question of model-parameter stability can he studied in some detail for Nigardsbreen, because the measurements there cover a long time period which can be divided into sub-sets with significantly different temperature and precipitation. In order to study the stability of the degree-day factors and the temperature gradient found by non-linear minimization, the calibration was performed for the 12 first years and the 12 last years separately. The results are shown in Table 7.
The value of the degree-day factor for ice is evidently quite stable (cf. Table 4). The variations in the degree-day factor for snow and the temperature gradient are greater. The reason for this difference is that there is a strong correlation (r = 0.94) between the parameters dds and grt both for Nigardsbreen and for Sátujökull (cf. Reference SigurðssonJóhannesson and others, 1993). As a consequence of the high correlation between dds and grt, data errors and simplifying assumptions in the model will lead to greater uncertainty in the estimates of the degree-day factor for snow and in the temperature gradient than in the degree-day factor for ice.
The effect of temperature on the estimated parameter values was further investigated by sorting the massbalance years according to the average temperature of each year. The model parameters were then estimated for the 12 coldest and the 12 warmest years separately. The results are shown in Table 8.
Again, variations in the degree-day factor for snow are greater than variations in the degree-day factor for ice (cf. Table 7). There is no significant reduction in the estimated degree-day factors with temperature which would he expected if the sensitivity of glacier mass balance to temperature was highly overestimated by degree-day mass-balance models. The changes in the parameter values could theoretically arise from random mass-balance fluctuations with variance similar to the variance of the residuals in Tables 7 and 8 (i.e. the changes in the parameter values are consistent with the computed standard error of the parameter estimates).
Best-fit degree-day factors and temperature gradient for Nigardsbreen determined fir two 12 year periods. Taw is the average temperature of the mass-balance years in each period at Bøerkehaug
The variance of the residuals corresponding to all 24 years is reduced from 0.39 m2w.e. (cf. Table 4) to 0.37 m2 w.e. by using separate parameters for two 12 year periods (cf. Table 7) and to 0.36 m2 w.e. by using separate parameters for the 12 coldest and 12 warmest years (cf. Table 8). This is a relatively small reduction in the variance for an effective doubling in the number of free model parameters, and supports our assertion that model parameters are stable over the temperature variations of these 24 years. Similar results have been obtained when mass-balance years were divided into four groups of 6 years sorted according to the average temperature of each year (Reference Jóhannesson, Sigurôsson, T and KennettJóhannesson and others, 1993).
The temperature differences considered above (1.1° C difference between 12 coldest/warmest years, 1.8° C difference between six coldest/warmest years) are of a similar magnitude to the climatic warming which is considered in the following section (2° C). We assume therefore these parameters will not change significantly for the climate scenarios considered here.
Climate Changes and Glacier Mass Balance
For the North Atlantic area, the climate scenarios which are used in the project “Climate change and energy production” specify a warming rate of 0.25° C per decade in mid-summer and 0.35° C per decade in mid-winter with a sinusoidal variation through the year. Precipitation is predicted to increase by 5% per degree of warming. The effect of a warmer climate on glacier mass balance is briefly analyzed here by computing the mass balance of the three glaciers for a temperature increase of 2° C, with and without an accompanying 10% precipitation increase. This is similar to the average warming that may occur over a time period of approximately 70 years to the scenario.
Figure 7) shows the average mass balance of Sátujökull, Nigardsbreen and Qamanärssûp sermia (solid curves, cf. Figures 1, 3 and 5), together with the mass balance predicted for 2° C warmer climate with the same precipitation and for 2° C warmer climate with 10% more precipitation. This figure shows that the accumulation of Qamanärssûp sermia is rather insensitive to the climate change but the predicted winter balance of and Nigardsbreen at the lower elevations is reduced even when the precipitation is increased by 10%, because a larger fraction of the precipitation falls as rain in a warmer climate. The summer balance (yearly ablation for Qamanärssûp sermia) and the net balance of all glaciers are predicted to decrease (become more negative) by between 2.4 and 2.9 m w.e. year−1 at the lowest elevations, and by between 1.1 and 1.3 m w.e. year−1 at the highest elevations. A warming of 2° C is predicted to lead to a 220 m or 180 m rise in the equilibrium-line altitude of Sátujökull, for no precipitaction change or 10% precipitation increase, respectively. The corresponding rise in the equilibrium-line altitude for Nigardsbreen is 280 m or 230 m, respectively. All the points on Qamanärssûp sermia which are considered here are below the equilibrium line and the rise in the equilibrium-line altitude is therefore not determined for Qamanärssûp sermia.
At an elevation of 790 m a.s.l., the yearly ablation at Qamanärssûp sermia is predicted to increase by 1.9 m w.e. This is in good agreement with the previous results of Reference Braithwaite and OlesenBraithwaite and Olesen (1993) who predicted that the yearly ablation at stake 751 at 790 m a.s.l. is increased by about l m w.e. ° C−1 of warming using a degree-day model. Reference Braithwaite and OlesenBraithwaite and Olesen (1990a) calculated a somewhat lower increase in the melting during the months June-August using an energy-balance model.
Best-fit degree-day factors and temperature gradient for Nigardsbreen determined for the 12 coldest and the 12 warmest mass-balance years
Modelled winter, summer and net balance of Sátujökull (a) and Nigardsbreen (b) and modelled yearly accumulation, ablation and mass balance of Qamanärssûp sermia (c) under the present climate (solid curves) and a warmer climate (dashed curves). For each set of three curves of the same line type, the lefthand curve is the summer balance (yearly ablation for Qamanärssûp sermia), the middle curue is the net balance and the righthand curve is the winter balance (yearly accumulation for Qamanärssûp sermia). The curves for the present climate are identical to the curves labelled “Average” in Figure 1, 3 and 5. Model parameters are given in Tables 1-6.
The above results were computed with a uniform distribution of the warming within the year. Moderately higher warming during the winter months according to the scenario described above gives essentially the same results as uniform warming if the temperature increase during the summer is the same in both cases. Recent coupled GCM model computations (Reference Manabe, Spehnan and SloufferManabe and others, 1992) indicate that there is little seasonal difference in the warming in the North Atlantic area, but winter warming over the Arctic Ocean and north Greenland is predicted to be much higher than warming during the summer. Reference Bøggild, Reeh and OerterBøggild and others (1994) showed that the large seasonal difference in the warming, which is predicted for north Greenland, may lead to a significant increase in snow accumulation which may partly compensate for the effect increased ablation on the net balance of the ice sheet.
Climate Changes and Glacier Run-Off
It is possible to estimate the hydrological consequences of changes in glacier mass balance due to climate change from the mass-balance computations which are described in the previous section. Figure 8) shows glacier run-off at three elevations from Sátujökull (a) and Nigardshreen (b) for the present climate, together with the computed runoff for 2° C warmer climate with the same precipitation and for 2° C warmer climate with 10% more precipitation. The results are obtained by integrating the model for the calendar months starting on 1 October for Sátujökull and for 30 or 31 d periods starting on 10 September for Nigardsbreen. Figure 8a) and b shows the total run-off as a function of time and Figure 8c) and d shows the computed increase in run-off due to the warming. The curves for Sátujökull (Fig. 8a) are difficult to distinguish from each other because the run-off curves at 1400 and 1800 m a.s.l. for the warmer climates almost coincide with the run-off curves at 1000 and 1400 m a.s.l. for the present climate.
Figure 8) shows that the run-off during the ablation season increases by 50-150 1 s−l km−2 for Sátujökull and by 50-200 1 s−1 km−2 for Nigardsbreen. As expected, the change is greater at the lower elevations. The onset of glacier run-off is about 1 month earlier and the ablation season is about 2 months longer in the warmer climates than under present climatic conditions. The 10% precipitation increase has only a marginal effect on run-off.
Discussion and Conclusions
The degree-day mass-balance model described above works well for three glaciers in significantly different climate regimes. With the long time series of mass-balance data from Nigardsbreen, it was shown that the model parameters are relatively stable m tune, and that climate changes within the measurement period do not lead to a modification of the parameters. This indicates that the model is suitable for a study of the effect of climate changes on glacier mass balance. The model explains over 80% of the variance in measured summer or yearly net balance at the individual elevations on the glaciers and it furthermore explains between 60% and 80% of the measured year-to-year variance in the yearly summer or net balance averaged over all the elevations on each glacier. Model performance is worse for winter balance, particularly for Sátujökull where the model explains only 42% of the variance. It can be seen in Figure 1) that the winter-balance curves are not well modelled using a linear precipitation gradient but that the form of the curves is reasonably constant from year to year. Thus, prediction of winter balance could probably be improved using an alternative (non-linear) precipit-ation distribunon with elevation. This pattern could be defined by elevation-dependent correction factors determined by least-squares fitting as described by Reference Bøggild, Reeh and OerterBøggild and others (1994).
Glacier run-off as a function of time from the beginning of the mass-balance year for Sátujökull (a) and Nigardsbreen (b). The time axis for Sátujökull is labelledwilh the month, but the time axis for Nigardsbreen is labelled with the number of 30-31 d periods starting on 10 September. Solid curves are for the present climate, short-dashed curves are for a° C warmer climate with the same precipitation and long-dashed curves ore for a 2° C warmer climate wilh 10% mori precipitation, (c) and (d) show the increase in run-off due to the warmer climates. Results from three elevations are plotted for each glarier. The curves are labelled with the elevations in m a.s.l. The results, with and without the precipitation increase, are often so close that the short- and long-dashed curves cannot be distinguished.
The degree-day factors found here are in relatively good agreement with parameter values which have previously been used for hydrological and glaciological models in Iceland, Norway and Greenland. The degree-day factors for ice for Sátujökull and Qamanärssûp sermia are higher than for Nigardshreen, which may be partly explained by the windiness of the climate in Iceland and Greenland. High wind speeds tend to increase sensible-heat flux from the air to the glacier surface. The calibration indicates a temperature lapse rate between 0.5° and 0.7° C per 100 m for the glaciers. This is close to the saturated adiabatic lapse rate near the surface of the Earth at temperatures near 0° C and similar to values which are most often used in hydrological models (0.6-0.7° C per 100 m).
Model calibration is perhaps the most critical part of and glaciological modelling studies. If the modelsare too unconstrained, it may be possible to fit the measurements quite well for a wide range of model parameters, resulting possibly in an unrealistic model response to changes in input data. It is therefore important to make maximum use of available measurements to constrain the model parameters as much as possible. This is the reason that the model calibration performed here is based on yearly mass-balance data instead of an average over the measurement period. The year-to-year variations in the measurements constrain the model parameters and the degree to which the model reproduces these variations is a check of the validity of the calibration.
Calibration of degree-day glacier mass-balance models is somewhat different from calibration of traditional hydrological degree-day snowmelt models which are calibrated using river discharge measurements. Glacier mass-balance measurements have very coarse time resolution, with typically only two measurements per year at each altitude. However, glacier mass-balance measurements tend to have a good altitude coverage, often covering the entire altitude range of the glacier. Glacier mass-balance models are therefore suitable For estimating the variation of the run-off with elevation, but suitable for estimating the time-dependence of the melting. River-discharge measurements, on the other hand, usually have good tine resolution but poor altitude and spatial coverage. Therefore, hydrological models are suitable for analyzing time-dependent processes which affect the melting, i.e. storage of liquid water in the snow-pack, the effect of ageing of the snow on albedo, refreezing of meltwater, etc., but they are not particularly suitable for analyzing altitude-dependence of parameters. Combined mass-balance/run-off models calibrated using mass-balance data and high-resolution discharge data can lead to improved estimates of model parameters, together with estimates of their time-dependence.
Acknowledgements
This study was carried out as a part of the project “Climate change and energy production”, a joint project between Norway, Sweden, Finland, Denmark and Iceland sponsored by the Nordic Council of Ministers. Meteorological data were provided by the Icelandic, Norwegian and Danish meteorological institutes. Massbalance data from Qamanärssûp sermia were provided by R. Braithwaite (GGU), who also gave useful advice during the course of this study.
