Kjetil Melvold

Glaciers in Svalbard: mass balance, runoff and freshwater flux

Gain or loss of the freshwater stored in Svalbard glaciers has both global implications for sea level and, on a more local scale, impacts upon the hydrology of rivers and the freshwater flux to fjords. This paper gives an overview of the potential runoff from the Svalbard glaciers. The freshwater flux from basins of different scales is quantified. In small basins (A < 10 km2), the extra runoff due to the negative mass balance of the glaciers is related to the proportion of glacier cover and can at present yield more than 20% higher runoff than if the glaciers were in equilibrium with the present climate. This does not apply generally to the ice masses of Svalbard, which are mostly much closer to being in balance. The total surface runoff from Svalbard glaciers due to melting of snow and ice is roughly 25 ± 5 km3 a?1, which corresponds to a specific runoff of 680 ± 140 mm a?1, only slightly more than the annual snow accumulation. Calving of icebergs from Svalbard glaciers currently contributes significantly to the freshwater flux and is estimated to be 4 ± 1 km3 a?1 or about 110 mm a?1.

On the Net Mass Balance of the Glaciers and Ice Caps in Svalbard, Norwegian Arctic

Abstract The ice masses of Svalbard cover an area of ca. 36 600 km2, and are thus among the largest glaciated areas in the Arctic. Annual mass balance measurements have been carried out on several Svalbard glaciers over up to 30 yr. However, these glaciers extend over only 0.5% of the total ice-covered area. The measured mean net balance has been negative and no changing trend has been observed. On some glaciers and larger ice caps, the mean net balance has also been measured at different altitudes by detecting radioactive reference layers from nuclear fallouts in 1963 and 1986 in shallow ice cores. The net balance/altitude curves have been estimated for thirteen different regions in Svalbard, and combined with digital elevation models of all Svalbard ice masses used to calculate the net balance in each 100-m altitude interval. The net loss of mass through iceberg calving was estimated and appears to be an important component of the net mass loss from Svalbard ice masses. The overall total net balance is slightly negative, −4.5 ± 1 km3 yr−1, giving a specific net balance of ca. −120 ± 30 mm yr−1 over the archipelago. The contribution of ice caps and glaciers on Svalbard to global sea-level change is, therefore, close to 0.01 mm yr−1 as an average value over the last 30 years, which is less negative than former estimates.

Evolution of a surge-type glacier in its quiescent phase : Kongsvegen, Spitsbergen, 1964-95

Kongsyegen is a 102 km 2 sub-polar polythermal) surge-type glacier in northwest Spitsbergen, Svalbard, It surged just before 1948 and is currently in its quiescent phase. Measurements of surface geometry since 1966 show a retreat of the front and strong thinning of up to 75 m in the ablation area, and a build-up of up to 32 m in the accumulation area. Present-day annual velocities along the glacier are low, from 1.4 up to 3.6 m a -1 The measured mean net balance for the period 1987-94 and the balance reconstructed back to 1967 show a weak positive balance of about 0.1 m w.e. The measured actual ice flux is low and the mass transfer down-glacier at the ELA is only about 3 20% of that required for steady state. Thus, the glacier is building up towards a new surge. The total thickening rate on Kongsvegen is somewhat higher than in other cases from Svalbard, but it is small compared with other well-studied surge-type glaciers in Alaska and he Pamirs. This relatively low rate of change is a function of the low accumulation rate and the relatively cold climate compared to other areas and is common for surge-type glaciers in Svalbard.

A mean net accumulation pattern derived from radioactive layers and radar soundings on Austfonna, Nordaustlandet, Svalbard

We present the snow-accumulation distribution over Austfonna, Nord-austlandet, Svalbard, based on 29 shallow ice cores that were retrieved from this ice cap during 1998 and 1999. Mean annual net accumulation is deduced from radioactive layers resulting from the 1954-74 atmospheric nuclear tests (maximum in 1963) and the Chernobyl accident (1986). The Chernobyl layer was located in 19 ice cores in the accumulation area, and the nuclear test layer was located in two deeper ice cores. In addition, the spatial variation of the depth of winter 1998/99 snowpack was mapped using snow probing, ground-penetrating radar methods and pit studies. The altitudinal gradient of the mean annual net mass balance and the altitude of the mean equilibrium line are determined along five transects ending at the top of the ice cap. The mean annual net mass balance and the equilibrium-line altitudes show a high degree of asymmetry between the western and eastern parts of Austfonna, in accordance with the distribution of winter accumulation. Large interannual variations of the accumulation exist. However, the study of the mean annual net mass balance shows no trend for two different time periods, 1963-86 and 1986 to the date of the drillings (1998/99).

Flow field of Kronebreen, Svalbard, using repeated Landsat 7 and ASTER data

Abstract Knowledge about the spatio-temporal distribution of fast-flowing Arctic glaciers is still limited. Kronebreen, Svalbard, in particular, includes the confluence − and the dynamic interplay − of the fast-flowing Kronebreen and the currently slow-flowing Kongsvegen. In this study, image-matching techniques on the basis of repeated Landsat 7 Enhanced Thematic Mapper Plus (ETM+) pan and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite data are applied in order to derive surface velocity fields of the lowermost 10 km of Kronebreen for the annual periods 1999/2000, 2000/01, 2001/02 and a 40 day period around July 2001. This work perfectly complements differential synthetic aperture radar interferometry (DInSAR) studies available for Kronebreen. A complete surface velocity field is now available from combining the DInSAR studies for the upper part of the glacier and the optical image-matching study presented here. The data obtained within this study are also compared to velocity data of 1964, 1986, 1990 and 1996. As also suggested by previous studies, a significant spatio-temporal variability of the spring/summer and annual ice speeds becomes evident.

Geometry changes on Svalbard glaciers: mass-balance or dynamic response?

Abstract The geometry of glaciers is affected by both the mass balance and the dynamics. We present repeated GPS measurements of longitudinal altitude profiles on three glaciers in Svalbard and show that surface altitude changes alone cannot be used to assess the mass balance. The three measured glaciers are in different dynamic modes, and the observed changes in geometry are strongly affected by the dynamics. Nordenskiöldbreen shows no significant change in the geometry, indicating that the mass balance is in steady state with the dynamics. On Amundsenisen the surface has lowered by 1.5–2.0 ma–1 in the lower part of the accumulation area at 520–550m a.s.l., indicating that the ice flux is higher than the mass-balance input, probably due to a surge advance of the glacier further downstream affecting the higher part of the drainage area. On Kongsvegen the opposite situation was found. Here the geometry of the profile showed a clear build-up of 0.5 ma–1 in the accumulation area and a lowering of 1 ma–1 in the lower part of the ablation area. The ice velocity is very low, giving a negligible vertical velocity component and an ice flux that is far smaller than the mass-balance flux, indicating that the glacier is building up towards a surge advance. Our results show that if mapping of height changes is to be used to monitor the response of the glaciers to climate change, both surface net mass-balance data and dynamic data are needed.

The distribution of snow accumulation across the Austfonna ice cap, Svalbard: direct measurements and modelling

We present an analysis of the spatial variability in the snow accumulation on the Austfonna ice cap in Svalbard, Norway, based on the results of field investigations conducted in the spring of 1999, 2004 and 2005. During the campaigns ground penetrating radar measurements at 500 and 800 MHz were collected along profiles, along with additional manual snow sounding and pit stratigraphy work. The analysis of the data reveals a consistent pattern in the spatial distribution of the snow accumulation over the ice-cap, and therefore enables the compilation of a multiple regression model of the snow distribution across the ice cap. Once validated, the results of the model complement the information derived from direct measurements, with an accumulation index for every point on the ice cap, thus enhancing the accuracy of future mass balance studies and dynamic modelling of Austfonna. As very few direct meteorological measurements are performed in the eastern part of Svalbard, the accumulation measurements on the Austfonna ice cap provide valuable integrated information about winter precipitation in this region.

Assessment of interannual variations in the surface mass balance of 18 Svalbard glaciers from the Moderate Resolution Imaging Spectroradiometer/Terra albedo product

[1] We estimate annual anomalies of the surface mass balance of glaciers on Svalbard for the period 2000–2005 (six years), by calculating the so-called ‘‘satellite-derived mass balance’’ (Bsat) from time series of satellite-derived surface albedos. The method needs no other input variables. Surface albedos are extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS)/Terra albedo product. We validate the MODIS albedos by comparing them with in situ measurements on Kongsvegen, and we find a low root-meansquare error of 0.04 for higher-quality MODIS data. Confidence in the MODIS product is also provided by realistic profiles of albedo along glacier centerlines. We apply the method to 18 glaciers that are evenly distributed over the archipelago. Correlation coefficients of time series of Bsat and direct measurements of the annual mass balance on Kongsvegen and Hansbreen are highly significant (0.94 and 0.82, respectively). Moreover, spatial distributions of the anomalies for individual years are coherent. Disadvantages of the method are that absolute values of the mass balance cannot be determined and that the interannual variability is underestimated. The latter might be corrected by equations to be established with mass balance models.

Mass Balance Methods on Kongsvegen, Svalbard

On the glacier Kongsvegen (102 km2) in northwest Spitsbergen, Svalbard, traditional mass balance measurements by stake readings and snow surveying have been conducted annually since 1987. In addition, repeated global positioning system (GPS) profiling, shallow core analysis and ground-penetrating radar (GPR) surveying have been applied. The purpose of this paper is to evaluate the input from the different methods, especially the GPS profiling, using the results from the traditional direct method as a reference. The annual flow rate on Kongsvegen is low (2 − 3 m a−1), and the emergence velocity is almost negligible. Thus the geometry changes of the glacier, i.e. the change in altitude per distance from the head of the glacier, should reflect the change in net balance of the glacier. The mean annual altitude change from the longitudinal, centreline GPS profiles was compared to the direct stake readings and showed a very good agreement. On Kongsvegen the measured actual ice flux is so low that the mass transfer down-glacier at the mean equlibrium line altitude is less than 10% of what is needed to maintain steady-state geometry. This is clearly shown in the changing altitude profiles. GPS profiling can be used on large glaciers in remote areas to monitor geometry changes, ice flow and net mass balance changes. However, it requires that the centreline profile changes are representative for the area/altitude intervals, i.e. that the accumulation and ablation pattern is evenly distributed. For this purpose the GPR surveying quickly gave the snow distribution variability over long distances. Shallow cores drilled in different altitudes in the accumulation area were analysed to detect radioactive reference layers from the fallout after the Chernobyl accident in 1986, and showed very good agreement to the direct measured net balance. Thus older reference horizons from bomb tests in 1962 could be used to extend the net balance series backwards.

Multiscale spatial variability of lidar-derived and modeled snow depth on Hardangervidda, Norway

Abstract This study presents results from an Airborne Laser Scanning (ALS) mapping survey of snow depth on the mountain plateau Hardangervidda, Norway, in 2008 and 2009 at the approximate time of maximum snow accumulation during the winter. The spatial extent of the survey area is >240 km2. Large variability is found for snow depth at a local scale (2 m2), and similar spatial patterns in accumulation are found between 2008 and 2009. The local snow-depth measurements were aggregated by averaging to produce new datasets at 10, 50, 100, 250 and 500 m2 and 1 km2 resolution. The measured values at 1 km2 were compared with simulated snow depth from the seNorge snow model (www.senorge.no), which is run on a 1 km2 grid resolution. Results show that the spatial variability decreases as the scale increases. At a scale of about 500 m2 to 1 km2 the variability of snow depth is somewhat larger than that modeled by seNorge. This analysis shows that (1) the regional-scale spatial pattern of snow distribution is well captured by the seNorge model and (2) relatively large differences in snow depth between the measured and modeled values are present.