Introduction
An ablation climate study was made at the margin of the Greenland ice sheet in Kronprins Christian Land, eastern north Greenland (Fig. 1). The study was part of a 2 year programme on world sea-level changes supported by the European Community. The specific objective of the study was to collect data on climate, ablation, radiation and ice temperature to estimate the sources and sinks of the energy balance during the period of measurements (8-27 July 1993) and to compare conditions with those found in West Greenland (Reference Braithwaite, Olesen and OerlemansBraithwaite and Olesen, 1989, Reference Braithwaite and Olesen1990).
The fieldwork was carried out at two study sites. One was located on the tundra in front of the ice sheet (base camp). The other was established on the ice sheet (glacier station). The glacier station was located about 300 m inland of the ice margin at an elevation about 50 m above the ground-level at the edge. The site was well-exposed without marked slope, i.e. it represents a large area close to the ice-sheet margin. The surrounding area was covered by hummocks about 0.5 m in height and 5 m in wavelength, i.e. the surface topography is rougher than the site in West Greenland (Camp IV; Fig. 1), where Reference AmbachAmbach (1963) made energy-balance studies. The glacier station was a walk of only 10 min from the base camp and was visited several times daily during the experiment. The position of the glacier station was determined by the global positioning system (GPS) to be latitude 79 °54 ’43” N and longitude 24 °04 ’25” W at an altitude of 380 m a.s.l. Ablation was determined as the average of readings at ten stakes drilled within an area of about 100 m2. Two thermistor strings were also installed to measure temperature gradients in the top 3 m of the glacier surface. At the same time as the ablation measurements, meteorological data air temperature, relative humidity and wind speed) were obtained once-hourly and recorded by data loggers at base camp and at the glacier station. Radiation conditions were studied by continuous logging of global radiation and all-wave (short-and longwave) incoming radiation at the base camp, while albedo was measured at the glacier station. Weather conditions were rather stable throughout the whole field period with nearly continuous sunshine, little cloud amount, nearly constant temperatures (daily means of 3-6 °C) and strong winds from the ice sheet (daily means of 3-8 m s−1). Cloud observations (amount and type) were made six times a day using World Meteorological Organization classification schemes. The interior of Kronprins Christian Land is obviously drier and has more sunshine than the coast in general, e.g. Station Nord is notorious for fog and late-lying snow. The tine weather conditions during the period of observation were probably not completely typical as there was a strong anticyclone over Greenland for most of the field period (Reference Geb and NaujokatGeb and Naujokat. 1993).
Reference AmbachAmbach (1963) showed for die upper ablation area (Camp IV) that the energy gain at the surface in summer is caused mainly by the radiative fluxes as the sensible-and latent-heat fluxes almost cancel each other. At ΕΤΗ Camp (Reference OhmuraOhmura and others, 1991, Reference Ohmura1992; Fig. 1), which is located close to the mean equilibrium-line altitude (ELA), net radiation provides almost all of the ablation energy (Reference Ohmura, Konzelmann, Forrer, Abe-Ouchi and ToritaniOhmura and others, in press), Radiation was also expected to be the main source of ablation energy at the present study site. In the lower ablation area of a glacier or an ice sheet, large variations of the surface reflectivity (albedo) can be expected. This is due mainly to different ice conditions, e.g. age. surface undulation, contamination, slope and orientation. These variations in albedo affect the energy balance and the ablation rate, and the relation between small-scale variations of ablation and albedo is treated as the main point of the present paper.
The smaller, less important energy-balance components, such as turbulent fluxes and conductive-heat flux in the ice. were evaluated by methods requiring only light-weight equipment. Reference MorrisMorris (1989) suggested that accuracy in estimating energy-balance components should be commensurate with their relative importance in the overall energy balance. We extend this concept by replacing “accuracy” with “logistic cost”. For example, the turbulent components were estimated from simple climatic data using profile-aerodynamic formulae Reference AmbachAmbach, 1986), thus saving the weight of heavy instruments. By analogy with preliminary geological surveys, such an approach may be termed a reconnaissance energy-balance study and it is hoped that the present results legitimate the approach.
Ablation
Previous experience has shown that ablation measurements involve considerable error and that ablation itself varies greatly even on a scale of metres (Reference Olesen, Braithwaite and OerlemansOlesen and Braithwaite, 1989), suggesting that mam measurements are needed to obtain a representative value for a site.
The ten stakes were read daily close to 1900 h local time (about 1715 h solar time). As stakes are usually surrounded In an “ablation hollow” of 0.1-0.3 m diameter, stake readings were made by measuring from the top of each stake to a straight edge laid on the upstream side of the stake which defines an average level of the surrounding surface. Differences in successive daily readings are converted into ablation values assuming an ice density of 900 kg m−3. We find it convenient to treat ablation as a positive rather than negative quantity as recommended by Reference Anon.Anon. (1969) to avoid the clumsiness of having to say that (negative) ablation decreases with increasing availability of energy.
The stakes were placed so that five were located in “dark” areas and five in “light” areas. These were sections of the ice surface which had “dark” and “light” color, respectively, based on visual inspection at the beginning of the field work. The difference between the two kinds of surface appears to be whether debris is spread evenly over die surface (“dark”) or concentrated into a few deep cryconite holes surrounded by very white ice (“light”).
Ablation variations From day-to-day and between stakes were analysed by the simple linear model of Reference LliboutryLliboutry (1974):
where yjt is the ablation at stake j on day t. αj is the mean ablation at stake j. βt is the ablation deviation which is the same for all stakes on day t and is the error in the model which is assumed to be random. Measurement errors contribute to the error term as long as they are random but if they are systematic, i.e. shared by all stakes, they contribute to the ablation deviation. According to Reference Braithwaite, Olesen and OerlemansBraithwaite and Olesen (1989), readings at the individual slakes should be approximately equally correlated with the ablation deviation βjt. At ablation stake B, however, the correlation coefficient between the ablation deviation and data is only 0.36 while correlations for other stakes are in the range 0.83-0.95. Unfortunately, stake Β was located on a steep slope which made it difficult to use the straight-edge method. Because of inconsistency with the other stakes, data from that stake are dropped from any further analysis. The ablation variations at the remaining nine stakes are shown in Table 1. The error in the linear model for the nine-stake, 20 d matrix is surprisingly low with a standard deviation of only ±5 kg m−2 d−1 which is certainly small compared with the range ±13 to ±19 kg m−2 d−1 quoted by Reference BraithwaiteBraithwaite (1985, p. 21-22).
The hypothesis that the four “dark” stakes (A, C, D and G) have the same mean ablation as the live “light” stakes (E, F. H. I and J) is tested with the Student t test (Reference KreyszigKreyszig, 1968, p. 219). The mean and standard deviations for the two data sets are 44 ± 9 kg m −2 d−1 and 39 ± 10 kg m−2d−1, respectively. The mean values are different (at 5% significance level) and it can be concluded that “dark” stakes have significantly greater ablation than “light” stakes. The time series of ablation energy for “light” and “dark” stakes are presented in Figure 2.
The glacier surface consisted almost always of a 20-30 mm thick “weathering crust” (Reference Müller and KeelerMüller and Keeler, 1969) with a relatively low density. The effect of this crust on the evaluation of ablation deserves discussion. The ablation in a period is the loss of material in the surface layer of thickness h relative to a stake fixed in the ice. The total surface-layer thickness h consists of a layer of ice of thickness hi, and a layer of crust of thickness hc. The ablation is given by:
where Δh is the change in surface-layer thickness, Δhc is the change in crust thickness, and ρv and ρi are the densities of crust and ice, respectively. As the depth of crust is not observed and its density is unknown, the ablation is evaluated using only the first term with an ice density of 900 kg m−3 which refers to slightly bubbly glacier ice. Some attempts were made in 1986 at Qamanârssûp sermia to measure the density of crust with a miniature corer (unpublished data from R.J.Braithwaite). Although these data cannot be very accurate, a crust density of 500-700 kg m−3 was indicated. Over the whole 20 d period of the field experiment with an average Δh. of 0.9 m, the second term, involving change in crust thickness, is entirely negligible. For the daily ablation readings, however, the second term could be a significant source of error although the weathering crust never completely disappeared in the present experiment as observed by Mâler and Keeler:1969) in their data sets. It should be noted that an error of this kind is not detected by the Lliboutry model because the error is more-or-less common to all stakes.
Radiative Fluxes
Due to the importance of radiative (fluxes in the energy balance, priority was given to accurate radiation measurements. At the base camp, global radiation and all-wave incoming radiation were measured directly with a Swissteco SS-25 pyranometer and a Swissteco ST-25 pyrradiometer, respectively. The longwave incoming radiation was then calculated as the difference between the all-wave incoming radiation and the global radiation, and compensated for the emission loss of the instrument (σTi 4, where σ is the Stefan-Boltzmann constant (= 5.67 × 10−8 Wm−2 K−4) and Ti the temperature of the pyrradiometer in kelvin. The cosine error of the upfacing pyranometers was corrected for zenith angle > 70 using a polynomial function, whereby only the part of direct solar radiation was taken into account. Based on measurements made under similar meteorological conditions at ΕΤΗ Camp in June 1990, it was assumed that 70% of global radiation is caused by direct solar radiation (Reference KonzelmannKonzelmann. 1994). Due to a possible underestimation of the all-wave incoming radiation due to thermal convection as proposed by Reference Ohmura and GilgenOhmura and Gilgen 1993), the longwave incoming radiation was corrected accordingly. Shortwave radiative fluxes at the glacier station were measured with a Swissteco SW-2 two-component pyranometer (albedometer). The instrument was mounted on a tripod at a height of about 1 m and was relocated every second day to sample as many types of surface as possible. In general, global radiation at the glacier station and at base camp agreed closely. However, to avoid possible systematic errors due to any tilt of the radiometer at the glacier station, global radiation measured at base camp was used to calculate albedo together with shortwave reflected radiation at the glacier station. Assuming melting conditions all day at the glacier station, longwave outgoing radiation was set at 316 Wm−2 according to the Stefan Boltzmann law.
The response time of the radiometers is about 5 s and the signals of the instruments were continuously recorded every 15 s. Every 30 min the average was computed and stored. The operation of the instruments and data loggers was checked several times a day. The uncertainty in the longwave incoming radiation was set at ±10 Wm−2 according to the correction method mentioned above and an instrument intercomparison performed by Reference DeLuisi, Dehne, Vogt, Konzelmann, Ohmura, Keevallik and KärnerDeLuisi and others (1993).
Turbulent Fluxes
The sensible-heat flux (SHF) is estimated by the method of Reference AmbachAmbach (1986) from measured wind speed and temperature:
where Ks is the exchange coefficient for turbulent-heat flux, Ρ is the atmospheric pressure, u2 is the wind speed at 2 m above the glacier surface and T2 is the air temperature at 2 m above the glacier surface. The glacier surface is assumed to be at the melting point. For a Prandtl-type neutral boundary layer with logarithmic profiles for wind speed, temperature and vapour pressure, the exchange coefficient is given by:
where cp is the specific heat of air with constant pressure (1005 J kg−1 Κ−1), k is von Karman’s constant (0.41), ρ0, is the density of air in the standard condition (1.29 kg m−3), b0 is the standard atmospheric pressure (1.013 × 105 Pa), z is the instrument height (2 m) and Z0w and Z0T are the roughness lengths for logarithmic profiles of wind and temperature, respectively.
The latent-heat flux (LHF) is similarly estimated from measured data for wind speed and vapour pressure:
where KL is the exchange coefficient for latent-heat flux and Δe2 is the difference between vapour pressure at 2 m above the glacier surface and that at the melting glacier surface. For the same assumptions as before, the exchange coefficient is given by:
where L is the latent heal of evaporation or sublimation as appropriate ((2.514 or 2.849) × 106J kg−1) and Z0e , is the roughness length for the logarithmic profile of water vapour.
Following Reference AmbachAmbach (1986), the roughness length for wind over an ice surface is assumed to be 2.0 × 10−1 m while the roughness lengths for temperature and water vapour are both assumed to be 6.0 × 10−6 m. As formulated by Reference AmbachAmbach (1986), the exchange coefficients KS ; and KL are valid for a neutral boundary layer while the air just over a glacier surface may be rather stable.
The turbulent-heat fluxes were calculated every hour according to Equations (1) and (6) using the air temperature, relative humidity and wind speed measured at 2 m. Some radiation heating of the temperature/humidity screen was noted whereby temperatures Suddenly rose when wind speeds were low. and the data were corrected by hand. Hourly values of the calculated turbulent fluxes and half-hourly values of the radiative fluxes were then summed up to daily values ending at 1900 h local time, so they could be compared directly with ablation measurements made at that time.
Ice Temperatures and Conductive-Heat Flux in the Ice
It was expected that the conductive-heat flux in the ice would be a significant energy sink in contrast to the ice-sheet margin in West Greenland where it was assumed by Reference Braithwaite and OlesenBraithwaite and Olesen (1990) to be negligibly small. Englacial temperatures were therefore measured every day to calculate the conductive-heat flux in the ice. The measurements were made at two strings with thermistors which were drilled at 0.5 m intervals down to a maximum of 3 m depth.
Measured temperatures were even lower than expected with the 0 °C isotherm only 0.2-0.3 m below the surface (even with the possible help of radiation warming of the thermistor cables). This probably accounts for the rather difficult drilling conditions encountered when placing the stakes, i.e. a wet drill barrel penetrating very cold ice and frequently jamming. During the period of measurements, the ice temperatures at the various depths rose by 0.2- 0.3 Κ d−1.
The rate of change of ice temperature with time is given by:
where ρi is density of ice (900 kg m), c is specific heat of ice (2009 J kg −1 Κ−1), Η(z) is the englacial-heat flux and z is the depth below the surface. Assuming that equals a constant C during the period of measurements, then H(z) is a linear function of depth:
where H0 is the heat flux through the glacier surface at z = 0 m. Equation (8) suggests that H0 can be calculated as the intercept in a regression equation of heat flux versus depth where the heat flux H(z) is calculated from englacial-temperature data according to the equation
where Κ is the thermal conductivity of ice, 2.1 Wm −1 Κ −1 according to Reference PatersonPaterson (1981, p. 186). The heat fluxes for each day are strongly correlated with depth with generally similar equations and, assuming that any differences are due to statistical fluctuations, the mean heat flux H0 for the whole period was calculated by combining all the available data into a single regression equation (Fig. 3). The intercept is 17.6 Wm−1 and the slope of the regression line corresponds to a constant temperature change of 0.21 Κ d −1 during the field period. Extrapolation of the regression line (Fig. 3) to greater depth suggests that the heat flux becomes zero at about 4 m depth but it is more likely that the approach to zero flux is asymptotic at greater depths than 4 m.
Energy Balance
The energy used for ablation ABL is given by:
where ABL is calculated from measured ablation using the latent heat of fusion (3.34 × 105 J kg−1), CHF is conductive-heat flux in the ice, SWR is shortwave net radiation, LWR is longwave net radiation, SHF is sensible-heat flux, LHF is latent-heat flux, and ERR is the total error due to measurement errors, simplifying assumptions and disregarded terms in the energy-balance equation. Defined as above, ERR can be regarded as an unknown, extra source of energy which is estimated from the residual of the other terms in Equation (10). The magnitude of ERR is a useful check on the accuracy of the other terms. For calculation of SWR an albedo of 0.48 was used. The condutive-heat flux in the ice was taken as 17.6 Wm−2. The time series of all components of the energy balance is presented in Figure 4.
The calculated energy balance for the 20 d is shown in Table 2 and summarised in Table 3. As discussed below, the main error is probably due to the turbulent fluxes so that they probably constitute more than the 26% of ablation energy shown in the table.
Error in the Calculated Energy Balance
On average, the error term ERR represents nearly 11% of the total energy balance. In view of the uncertainty of the various assumptions, it is good that the mean error is not even larger. For example, Reference Braithwaite and OlesenBraithwaite and Olesen (1990) found errors of −40 to + 29% of monthly ablation at Nordbogletscher and Qamanârssûp sermia. In terms of variance, the error corresponds to about one-third of the ablation variance. There is therefore a substantial source of both systematic and random error in the energy-balance calculation.
The error ERR has almost no correlation with the radiation terms SWR and LWR. On the other hand, the error has moderate correlation with the turbulent fluxes, i.e. r =+0.63 for SHF and r =-0.53 for LHF, which suggests that they may be the main sources of systematic error in the calculated energy balance. For example, if both SHF and LHF are multiplied by a factor of 1.5 before calculating the energy balance, the mean error is essentially reduced to zero (the error standard deviation is little changed). In order to assess the plausibility of a 50% error in the turbulent fluxes it is necessary to examine in more detail how they are calculated, i.e. with regard to stability and to the assumed surface roughness.
The exchange coefficient for sensible-heat flux (Reference AmbachAmbach, 1986) given by Equation (4) does not take account of the stability of the boundary layer over the melting glacier. Reference Price and DunnePrice and Dunne (1976) propose stability corrections in terms of the bulk Richardson number, and application of this reduces the sensible-heat flux in the present study by 5-18%, thus increasing the error in the calculated energy balance rather than reducing it.
The Reference AmbachAmbach (1986) formulation uses different roughness lengths for wind, temperature and water vapour (Over ice) while Reference Price and DunnePrice and Dunne (1976) and Reference MooreMoore (1983) assume the same roughness lengths (over snow). Putting Z0W = Z0T = Z(0e = 10−3 m in Equations (4) and (6) would exactly account for the mean error in the calculated energy balance. Such a roughness length falls well within the wide range of values quoted in the literature for ice, e.g. by Reference Grainger and ListerGrainger and Lister (1966), Reference Strelen and WendlerStrelen and Wendler (1968), Reference Wendler and WellerWendler and Weiler (1974), Reference PoggiPoggi (1977) and Reference Hay and FitzharrisHay and Fitzharris (1988), although Reference MorrisMorris (1989) suggests that some of the larger roughness values may be caused by slope errors. It appears, therefore, that the calculation of the energy balance is extremely sensitive to the assumptions about surface roughness in agreement with Reference MunroMunro: 1989). This is an important point for future attempts to couple ice-sheet models lo general circulation models (GCMs).
Effect of Albedo on Ablation
Measurements of albedo are normally made at one point and therefore depend strongly on the local surface conditions. For an area of 50 m × 100 m in the upper ablation zone (Camp IV), Reference AmbachAmbach (1963) gives a variation of the albedo from 0.34 to 0.58 for a particular day. Preliminary results from an altitudinal profile of albedo on bare ice close to the present study site show albedo variations of between 0.30 and 0.62 (personal communication from H. Oerter).
On a total of 16 d, albedo was measured at different places inside the area of ablation measurements. Parts with “dark” surface (8 d) and “light” surface (8 d) were selected subjectively. For the “dark” surface a mean albedo of 0.43 was found and the “light” surface showed a mean value of 0.53, while the mean albedo of all sites was 0.48. The albedo difference between “light” and “dark” surfaces is here 0.10. For a daily mean global radiation of 320 Wm−2, this gives a difference in calculated melt of 8.3 kg m−2 d−1 which is about 20% of the mean ablation over the whole data set (Table 1). The difference in mean ablation between the four “dark” and the live “light” stakes is 5 kg m−2 d−1 while the greatest inter-stake difference is 11 kg m−2 d−1 between stakes A and I (Table 1).
The above discussion refers to separate ablation and albedo variations within an area, and it would be interesting to extend the discussion to the relation between ablation and albedo for individual stakes. However, it is impossible to measure albedo exactly at a stake because the stake and the albedometer interfere with each other, and a more subjective procedure for estimating albedo at each stake is used. In this method, the surface conditions under the albedometer were compared on each day to the situations at nearby stakes, and the measured albedo value was assigned to one or more stakes that appeared to have a similar surface. The albedo for each stake was then determined as the average of the assigned albedo values. The estimated albedos of the four “dark” stakes (A, C, D and G) agree closely but there is quite a large range for the five “light” stakes. Average albedos for the groups are higher than previously stated, but this simply reflects the different sampling. The relation between mean daily ablation and the estimated albedo for the individual stakes is illustrated by Figure 5. The correlation shows a value of r = −0.82, which again suggests that albedo is a main factor in ablation variations within a small area.
The above indicates that a point measurement of albedo is not sufficient for accurate energy-balance calculations in ablation areas. Therefore, an albedo value for a larger area should be determined, either by measurements from high towers (Langleben, 1968) or by enough point measurements at ground level to sample small-scale albedo variations as attempted here.
Comparison with West Greenland
In West Greenland, energy-balance studies were carried out by Reference Braithwaite and OlesenBraithwaite and Olesen (1990) at Nordbogletscber (NBG) and Qamanârssûp sermia (QAM), and by Van de Reference Wal and OerlemansWal and Oerlemans (1994) near Søndre Strømfjord in environments similar to those found at the present study site (KCL). At all of these stations the net radiation was measured at one point and no albedo variation was taken into account. Reference AmbachAmbach (1963) performed a more sophisticated energy-balance study at Camp IV, but the station was located in the upper ablation area, and the albedo variations and their effect on ablation were unfortunately not analysed in detail. The energy balance from different experiments is summarized in Table 3. Net radiation is the main source of ablation energy at all locations. Turbulent fluxes are an important energy source in the ablation zone compared to the situation close to the ELA (Reference AmbachAmbach, 1963; Reference Ohmura, Konzelmann, Forrer, Abe-Ouchi and ToritaniOhmura and others, in press), corresponding to an ablation rate of about 10-20 kg m−2 d−1. In contrast to West Greenland where heat conduction into the ice is assumed to be small, it is an important heat sink in the present study. The possible ablation is reduced by about 4.5 kg m−2 d−1.
Ablation rates and degree-day factors for KCL, QAM and NBG are summarized in Table 4. Despite the low elevation of the station compared with the two West Greenland sites, the ablation rate is quite low. This is partly due to lower air temperature, but that is offset by a higher degree-day factor. At KCL a value of 9.8 mm d−1 °C−1 was found compared with 7.2 – 8.1 mm d−1 °C−1 at the West Greenland sites (Reference BrahlrwaiteBraithwaite, 1992), which may be the result of higher wind speeds. Reference Boggild, Reeh and OerterBoggild and others (1994) also got a higher degree-day factor (9.6 mm d−1 °C−1) at Storstrammen in northeast Greenland (Fig. 1).
Conclusions
Net radiation is the major source of ablation energy at the margin of the Greenland ice sheet in Kronprins Christian Land as in other parts of Greenland. Albedo varies greatly within a small area and, because of the high income of global radiation, causes relatively large differences in ablation at stakes close to each other. Small-scale albedo variations must therefore be carefully sampled to obtain representative albedo values for large-scale energy-balance calculation.
The calculated energy balance in the present study also has a substantial error which may be caused by underestimation of the turbulent fluxes using the aerodynamic formulae of Reference AmbachAmbach (1986). The reason is probably that the glacier surface is rougher than assumed. The conductive-heat flux in the ice is a substantial heat sink and reduces the energy available for ablation. The average ablation rate in Kronprins Christian Land in July is low compared with values in West Greenland, reflecting lower temperature although this is partly offset by a higher degree-day factor due to high wind speed. The degree-day factor is significantly higher in north Greenland than the ones found in West Greenland.
Acknowledgements
This paper is published by permission of The Geological Survey of Greenland (GGU). K. Schroff (Department of Geography, Swiss Federal Institute of Technology) and O. Olesen (GGU) prepared instruments. Logistic support was supplied by GGU’s camp at Centrum So. Kronprins Christian Land, led by GGU Statsgeolog N. Henriksen. The research is supported by the European Community under contract number EV5V-CT91-0051 which is coordinated by the Climate Research Unit, University of East Anglia. The first draft of this paper was prepared in November 1993 when R.J. Braithwaite was a guest of Professor A. Ohmura. Department of Geography. Swiss Federal Institute of Technology (grant No. 0-04-509-93). Professor Ohmura also made useful continents on the manuscript.