Melting Characteristics of Snow Cover on Tidewater Glaciers in Hornsund Fjord, Svalbard

In recent years, the Svalbard area, especially its southern section, has been characterised by an exceptionally thin snow cover, which has a significant impact of the annual mass balance of glaciers. The objective of this study was to determine melting processes of the snow cover deposited on 11 glaciers that terminate into Hornsund Fjord during the melting period of 2014. The study included analyses of snow pits and snow cores, meteorological data collected from automatic weather stations and Polish Polar Station Hornsund, and supervised classification of six Landsat 8 images for assessing the progress of snow cover melting. The calculated Snow-Covered Area (SCA) varied from 98% at the beginning of the melting season to 43% at the end of August. The melting vertical gradient on Hansbreen was −0.34 m 100 m−1, leading to surface melting of −1.4 cm water equivalent (w.e.) day−1 in the ablation zone (c. 200 m a.s.l. (above sea level)) and −0.7 cm w.e. day−1 in the accumulation zone (c. 400 m a.s.l.). Furthermore, the study identified several observed features such as low snow depth in the accumulation zone of the Hornsund glaciers, a large proportion of the snow layers (12–27%) produced by rain-on-snow events, and a frequent occurrence of summer thermal inversions (80% annually), indicating that the area is experiencing intensive climate changes.


Introduction
Snow cover plays a significant role in shaping the global climate system. It forms a highly-insulating layer, reflecting up to 90% of the solar energy [1]. In Polar Regions, intensive surface melting processes, and thus also radical changes in the albedo, may occur over just a few days, during the late-spring melt. Recent decades have seen noticeable changes in the climate worldwide, with the most pronounced transformations of the geographic environment observed in the Arctic, especially in its glacierized areas [2]. Given the high sensitivity of the European sector of the Arctic to meteorological conditions, the area is extremely important for understanding the climate processes that are taking place on a global scale.
At present, most glaciers worldwide experience mass loss [3]. Existing models used to determine the mass balance, retention and meltwater runoff into the ocean clearly lack precise and regular spatial measurements of the snow cover melting [4,5]. A majority of hydro-glaciological studies focus on small and medium-size glaciers [6], which have a greater contribution to feeding the ocean with freshwater than the Antarctic and Greenland ice sheets [7].
Increased ablation results in an intensified supply of meltwater, which has strong implications for the entire Arctic environment, for example, by increasing the basal sliding layer of glaciers, transporting mineral and organic matter beyond the basin boundaries, and locally changing the   The Hornsund region has an oceanic climate, with precipitation prevailing over evaporation. From the west, local climate is influenced by the warm West Spitsbergen Current, increasing the amount of precipitation and reducing the annual amplitudes of air temperature. The east coast remains under the influence of cold sea currents from the north, resulting in local cooling [27]. The multi-annual  mean air temperature at the WMO (World Meteorological Organization) 01003 (the Polish Polar Station Hornsund) was −4.3 • C [28], while the mean precipitation totals in the same period was 434.4 mm, with a prevalence of liquid forms (44%), occurring mainly in summer months-July and August [29]. The increased air temperature in winter months and the frequency of rain-on-snow (RoS) events are an important indicator of ongoing local climate change [24,30]. For the years 1979-2010, the positive air temperature trend was 0.97 • C per decade [31].
Glaciological data: among the glaciers terminating into Hornsund, continuous mass balance monitoring is only conducted on Hansbreen. The data include measurements from six mass-balance stakes corresponding to the elevation zones: 50, 100, 200, 300, 400 and 500 m a.s.l., located along the central axis of the glacier. The measurements were taken c. once a month. In further surface melt calculations, both snow melt (variable density) and ice melt (constant density, ρ = 917 kg m −3 ) were distinguished. Changes of snow density for wet snow conditions have been estimated based on summer snow pits [24,32].
Determination of the general physical properties of the snow cover: A total of seven snow pits were completed in accordance with the International Classification for Seasonal Snow on the Ground [33]. Three were located in the ablation zone: Hansbreen (H4), Storbreen (SG), Flatbreen-the upper section of Hornbreen (FJ); three in the equilibrium-line region: Hansbreen (H6), Storbreen (SG), Flatbreen (FK); and one in the accumulation zone: Hansbreen (H9). The study mainly relied on snow depth, mean density, snow water equivalent (SWE), and a contribution of dense layers: Melt-Freeze Crusts (MFcr) and Ice Formations (IF). Bulk snow density was measured by means of a simple weight gauge (cross-section area: 50 cm 2 , length of tube: 60 cm, weight scale: 0.005 kg).
Additionally, snow cores were sampled in the accumulation zones of Storbreen (SI) and Flatbreen (FL), providing information about seasonal snow depth, mean density of snow cover and SWE. The coring system had a core diameter: 7 ± 0.2 cm. The samples were weighed using a spring weight scale with an accuracy of ±0.3%.
All the snow studies were carried out for dry snow cover in April-May 2014, at the end of the accumulation season.
Satellite data: detailed remote sensing analysis included six Landsat 8 satellite images covering the melting period of 2014. Selected spectral bands (1-5) of the Operational Land Imager sensor, with a wavelength range of 433-885 nm, were used for further analysis. The spatial resolution was 30 m per pixel.
Remote sensing in measuring snow cover surface classification: Satellite data allowed a comparison of the temporary Snow-Covered Area (SCA) at successive melting stages on the selected glaciers. The images were subjected to Top of Atmosphere (ToA) reflectance correction [34], and the glacierized areas were delimited. Detailed analyses were performed using colour compositions: in natural colour (bands 1, 2, 3) and false colour (bands 2, 3, 4).
In order to map the glacier surface, supervised classification was performed using the Maximum Likelihood algorithm [35]. For a better visualisation and interpretation purposes, a pan-sharpening function was applied, which increased the spatial resolution of the multispectral image bands from 30 m to c. 15 m. The output images showed six classes: snow, ice, shadowed snow, shadowed ice, debris and water. A cloud mask was applied for images with clouds impeding proper interpretation [36]. Mask was obtained from band 6 and any artefacts caused by similar intensity of the cloud pixels and the individual classes were corrected manually in the last stage of the post-processing.

Microclimate and Topographic Determinants of the Structure of the Snow Cover on Hornsund Glaciers and Its Melting Characteristics
The year 2014 was the second warmest year, after 2012, in the entire observation history in Hornsund, with a mean annual air temperature of −1.3 • C. The monthly average temperature for all months of the year was significantly higher than multi-annual records. This was particularly clear in winter, from January to March, with the deviation reaching 9.5 • C in February. It was also a dry year, with annual precipitation of 410 mm, which represented 90.6% of the multi-annual 1979-2013. The lowest precipitation was recorded in winter months, which are particularly important for snow cover formation. The monthly precipitation in February (7.6 mm) represented only 26.7% of the mean long-term values for this month [37]. This resulted in lower than average snow depth, notably in the accumulation zones of investigated glaciers.
The local thermal conditions are strongly determined by the character of the surface. This is particularly evident in the summer, when tundra, which is not insulated by snow, is strongly heated and undergoes intensive evaporation. This causes the number of positive degree days (PDD) to vary significantly from site to site ( Table 2). The accumulated total of PDD ranged from 164.3 (H9) to 585.4 (HOR), and positive mean diurnal air temperature was recorded nearly twice as frequently in Hornsund (47% of the year) as in the accumulation zone of Hansbreen (25% of the year). On account of continuing low air temperature in April caused by northerly and north-easterly cyclonic and anticyclonic situations, April was the only month without a PDD at any of the measuring stations.  The melting period in the glacierized area started relatively late, i.e., only towards the end of the second half of May, and lasted until the beginning of October. An unusual meteorological situation was observed in the summer months, especially in July (on 16 days), when strong thermal inversions occurred, with the mean diurnal air temperature in the accumulation zone 3.6 • C higher than the temperature in the ablation zone ( Figure 2).
In recent years, warm spells during winter months have been occurring more frequently, accompanied by liquid precipitation, so-called rain-on-snow events [29]. They are an important indicator of on-going climate change in the polar regions [38,39]. Occasionally, they even occur in highly elevated glacierized zones, and producing very compact layers in the snowpack-Melt-Freeze crusts (MFcr) and Ice Formations (IF), which slows down the meltwater percolation processes. This was the case on Hansbreen, where the share of these layers in the total snow depth exceeded 27% in the ablation zone, 20% at the equilibrium-line, and 12% in the accumulation zone ( Table 3). The mean snow density was similar across the sites and amounted to c. 400 kg m −3 . A noteworthy finding was the very low snow depth in the accumulation zones of the glaciers, on which the physical characteristics of the snow cover were measured in more detail. crusts (MFcr) and Ice Formations (IF), which slows down the meltwater percolation processes. This was the case on Hansbreen, where the share of these layers in the total snow depth exceeded 27% in the ablation zone, 20% at the equilibrium-line, and 12% in the accumulation zone ( Table 3). The mean snow density was similar across the sites and amounted to c. 400 kg m −3 . A noteworthy finding was the very low snow depth in the accumulation zones of the glaciers, on which the physical characteristics of the snow cover were measured in more detail.  The measurements taken at the mass-balance stakes in 2014 show that the melting vertical gradient along the Hansbreen axis depends on the elevation, reaching c. −0.34 m w.e. 100 m −1 . Figure  3 presents mean values of the surface melt rate per day, calculated for periods between observations. The melt rate at the lowermost ablation zone was three times higher than in the upper areas of the glacier, resulting in the lowering of the glacier by c. 2.44 m at the terminus. The highest melt rate was recorded in July, when most of the ablation zone lost its snow cover, which protects the glacier surface against increased absorption of shortwave solar radiation, due to higher albedo. In July and August, the snow cover melt rate in the uppermost part of Hansbreen (500 m a.s.l.) was higher than in the 400 m a.s.l. zone, which can be attributed to frequent thermal inversions and wind activity. The measuring point is located at the ice divide, whereby the air masses moving southwards cause snow transfer and reduce the snow cover.  The measurements taken at the mass-balance stakes in 2014 show that the melting vertical gradient along the Hansbreen axis depends on the elevation, reaching c. −0.34 m w.e. 100 m −1 . Figure 3 presents mean values of the surface melt rate per day, calculated for periods between observations. The melt rate at the lowermost ablation zone was three times higher than in the upper areas of the glacier, resulting in the lowering of the glacier by c. 2.44 m at the terminus. The highest melt rate was recorded in July, when most of the ablation zone lost its snow cover, which protects the glacier surface against increased absorption of shortwave solar radiation, due to higher albedo. In July and August, the snow cover melt rate in the uppermost part of Hansbreen (500 m a.s.l.) was higher than in the 400 m a.s.l. zone, which can be attributed to frequent thermal inversions and wind activity. The measuring point is located at the ice divide, whereby the air masses moving southwards cause snow transfer and reduce the snow cover.

Snow Cover Melting Characteristics of Hornsund Glaciers Based on Satellite Data
According to the field measurements and meteorological data, the accumulation period lasted until mid-May (see: Section 4.1). The first satellite image used in the analysis was from the end of May, when the surface melting processes had begun. At that time, the contribution of snow on the overall glacier surface reached 98%. The areas assigned to the ice class (2%) mainly include the highly

Snow Cover Melting Characteristics of Hornsund Glaciers Based on Satellite Data
According to the field measurements and meteorological data, the accumulation period lasted until mid-May (see: Section 4.1). The first satellite image used in the analysis was from the end of May, when the surface melting processes had begun. At that time, the contribution of snow on the overall glacier surface reached 98%. The areas assigned to the ice class (2%) mainly include the highly crevassed termini (Figure 4).
A noticeable reduction in the Snow-Covered Area (SCA) was seen at end of June, when the SCA was changing at a rate of −0.86% day −1 . Consequently, the ice-covered area was gradually expanding, and by the end of August, the percentage of the ice class was 11% higher than that of the snow class ( Figure 5). The total SCA reduction during the studied period equalled −0.56% day −1 . When comparing the individual glaciers, the reduction of SCA varied considerably, ranging from −18% (Mühlbacherbreen) to −74% (Körberbreen) with an average of −52%. Mühlbacherbreen was the only glacier with insignificant SCA reduction. Previous papers focusing on snow distribution on Svalbard glaciers highlighted that interior areas are characterized by the local continental microclimate, where lower air temperature persists throughout the year [30,40,41]. This dependency is visible both on the northern and southern coasts of the Hornsund Fjord ( Figure 6). The melting of the snow cover on Hornsund glaciers continued for several more weeks. However, the lack of later satellite data made further analysis impossible. The latest satellite image that was available, dated 21 September, had already been taken after the first snowfall episode.     From the beginning of August, a superimposed ice zone can be distinguished on the surface of some glaciers. During the classification, it was usually assigned to the accumulation zone (the snow class), with certain exceptions. However, it must be pointed out that separating superimposed ice exclusively by means of indirect methods is a highly subjective operation, requiring great interpretative experience and suitable in situ verification.
The setting that significantly improved the separation of snow from ice was the rejection of shortwave infrared bands (SWIR, band 6 and 7), even though they are widely used for the separation of cloud and snow [9]. In addition, band 9 (Cirrus), which is mainly used for cloud detection [34,42,43], was also disregarded in the classification since it generated an increased number of artefacts in areas characterised by low amount of clouds. However, rejecting the SWIR bands had a negative effect on automatic detection of the debris class-fortunately, the class was relatively infrequent, and the number of pixels which represented it was strictly controlled by decreasing the bias for the class and its threshold.
All the material available was characterised by a low value of sun elevation (Table 4), which highlighted the local topographic features of the glaciers, causing significant differences in the lighting of the surface and producing a high proportion of shadowed areas, reaching 12.5% in mid-July. Nevertheless, the classification allowed us to identify the type of the surface within those areas. Although, in general, shadowed snow was classified correctly, shadowed ice often included the debris class, and differentiating them was a major challenge in the classification process. The earliest discernible loose moraine material on the surface of glaciers (debris) appeared in early June. Throughout the melting period, the contribution of debris class was only several percent, but it contributed to melting in a significant way. The class consists of gravel and rocks that heat easily and have the lowest albedo in the study area. This generates positive feedback which accelerates the melting of the adjacent snow and ice cover. The largest accumulation of debris was observed in early August. However, debris is difficult to classify at times of advanced melting because areas of glacier ice with low image tonality, like wide glacier crevasses, were often included in the class.
During the melting period, temporary supraglacial water bodies may appear on the surface of the glaciers. Most of them are small melt ponds (<0.005 km 2 ), only occasionally reaching larger sizes. The largest supraglacial pond was located on the surface of Paierlbreen (0.21 km 2 ) and was formed after the snow cover started to melt intensively. The lifetime of such melt ponds mainly depends on the timing at which the glacial drainage system opens, followed by feeding with meltwater from the glacier surface. Small ponds were seen from the beginning of June until the end of August. Of all the glaciers, the greatest number of such water bodies was observed on Storbreen, where the unique relief of the glacier, i.e., a vast flat area with a small inclination, a low crevasse level, and numerous depressions, is favourable to their seasonal formation. The extra water class was added to the images dated 22 June and 15 July, when the supraglacial ponds were distinct and filled with meltwater, even though their share of the overall area was negligible (0.09%).
With a suitable image quality (no clouds, higher pixel intensity) and length of the glacier terminus, Landsat imagery can also be used for observing the evolution of subglacial drainage systems and for identifying the meltwater outflows, where the water, rich in mineral material, leaves the glacial basins. Such points are present in all Hornsund glaciers, most often at the side parts of their termini. Large amounts of suspended sediment are discernible through their much lighter image tone on the surface of the fjord. The earliest signs of the sediment in its initial form were observed at the beginning of June in front of the Mühlbacherbreen, Storbreen, and Hornbreen. As the melting advanced, the drainage channels were unsealed, and the migration of meltwater became even more noticeable. The greatest turbidity of Hornsund waters was observed in mid-July, especially in the bays facing the termini of the largest glaciers (Storbreen and Hornbreen). This could be seen in all images until the end of August. The largest amounts of sediment are transported from the glaciers in Brepollen, which can be explained by the presence of the long capes of Treskelodden and Meranpynten (in the north and south, respectively), which divide the fjord, preventing the Hornsund waters from mixing. Towards the end of September, the bays facing the termini of the glaciers slowly started to freeze over.
Clouds were a natural factor hindering correct identification of classes. In the summer, considerable cloudiness is typical for the Hornsund area [27]. This is the factor that most significantly affects the quality of satellite data used in remote sensing analyses. For some glaciers, a cloud mask was applied: at the beginning of July (Samarinbreen) and in mid-July (Mendelejevbreen and Kvalfangarbreen). Due to extensive cloudiness, the remote sensing analyses in 2014 had to be discontinued in late August.
Of all the glaciers under study, the greatest problems with accurate classification were posed by those located on the southern side of the Hornsund Fjord. We assume that morphometric differences have an important effect on the quality of classification. The southern glaciers cut strongly into narrow valleys and landforms that are surrounded by higher mountain ranges, which generate more surface shadowing, and produce local cloudiness. In addition, they have a much smaller surface area. All this caused major interpretation difficulties, decreasing the accuracy of the classification results of those glaciers. However, the final results of each Landsat image have a very high kappa coefficient [44] and overall accuracy of over 95% [35], see: Table 4. With a suitable image quality (no clouds, higher pixel intensity) and length of the glacier terminus, Landsat imagery can also be used for observing the evolution of subglacial drainage systems and for identifying the meltwater outflows, where the water, rich in mineral material, leaves the glacial basins. Such points are present in all Hornsund glaciers, most often at the side parts of their termini. Large amounts of suspended sediment are discernible through their much lighter image tone on the surface of the fjord. The earliest signs of the sediment in its initial form were observed at the beginning of June in front of the Mühlbacherbreen, Storbreen, and Hornbreen. As the melting advanced, the drainage channels were unsealed, and the migration of meltwater became even more

Discussion
The natural Arctic environment is undergoing radical transformation, which is mainly caused by increasing air temperature. The highest rate of warming in Europe, amounting to 3.9 • C per century, has been recorded on Svalbard [45], and out of the 35 years of meteorological observations in Hornsund, the investigated period has been the second warmest. Changes in temperature lead to increased frequency and intensity of warm spells and extension of the melting period [46,47]. The highest standard deviation from monthly mean values was seen in the winter months, which are crucial for the evolution of snow cover, corresponding with earlier studies [24,30,31]. Consequently, the snow depth in the accumulation zones of the glaciers in southern Spitsbergen, was lower in 2014, but not as low as in 2013 [30].
Winter warm spells are often accompanied by intensive rainfall, which freezes on the snow surface or just below it, creating layers of increased density, slowing down the percolation processes and further runoff of meltwater. Increased frequency of the rain-on-snow in recent years across the southern Spitsbergen is an important indicator of on-going climate change, also in other sectors of the Arctic [38,39,48]. RoS events started to be the indirectly investigated in the Hornsund area as early as the 1980s [49]. On Hansbreen, the contribution of the layers formed by those events to the snowpack is high, remaining above 25% in the ablation zone [24,30]. This translates into increased density of snow cover, which is c. 400 kg m −3 compared to the average for Svalbard of 370 kg m −3 [40].
Surface-based inversions, during which air temperature in the troposphere rises with height from the surface, are a frequent atmospheric phenomenon in high latitude climate, but are more characteristic for winters than summers [50,51]. Also, the inversion strength was much greater at that time [52]. In 2014, Hansbreen saw a reverse situation, with summer inversions accounting for 80% of the annual cases. In July, i.e., the warmest month, the phenomenon displayed the highest diurnal frequency, accelerating intensive melting processes in the upper areas of the glacial zones. A similar dependency was observed, e.g., in the Canadian Arctic [53].
The melting rate not only depends on air temperature, the energy balance of the surface, rain-on-snow events, or surface-based inversions, but also on surface reflectance. As a result of the differences between the albedo of snow and ice, the melt rate is lower during spells with the presence of continuous snow cover than in its absence, which means that the snowpack also plays the role of a natural protective layer against the intensification of melting processes [54]. Of all the glaciers in southern Spitsbergen, the longest series of glacial environment monitoring data is available for Hansbreen (since 1989), where several attempts have been already made to determine the snow cover melt rate in different seasons, using both statistical models based on positive-degree days [26], SR50 ultrasonic rangers [24,25], and mass-balance stake measurements [24]. In the above studies, the melt rate in the ablation zone (c. 200 m a.s.l.) ranged between 0.6 cm w.e. d −1 and 0.8 cm w.e. d −1 , which is half the value calculated in this paper. As a result, one can assume that in 2014, the other glaciers of the Hornsund area were also characterised by above-average melting compared to previous years.
Snow cover is also a key element determining the glacial mass balance [55]. In recent years, the declining accumulation of snow cover on the Hornsund glaciers has not been able to compensate for the losses caused by its melting, which results in negative mass balance [30,31]. Similar trends are also observed in other Arctic areas [3]. Apart from surface ablation, ice calving significantly contributes to overall mass loss and represents on average 21% of the summer balance of the Svalbard glaciers [19].
The boundaries of the glacial basins were compared to the most recent measurements [8]. This was used to calculate the changes in the geometry of the glaciers terminating in the Hornsund Fjord, which were caused by calving and terminus retreat. In the years 2004-2014, the surface area of the glaciers decreased by a total of 19.53 km 2 (−2.55% of the surface). The greatest changes were observed on Mendelejevbreen (−5.97%), and the least pronounced ones on Paierlbreen (−1.63%). Earlier studies in Sørkappland demonstrated a decrease in the areal extent by c. 18% between 1936 and 1991 [56]. In the scale of entire Svalbard, glaciers of Wedel Jarlsberg Land and Sørkappland are characterised by the lowest values of geodetic mass balance, −0.6-0.8 m w.e. y −1 [16], and are retreating faster, at c. 70 m y −1 [8], which implies that this area is among the fastest changing ones in the Arctic.
Another good indicator of the condition of glaciers in the context of recent climate change is the fluctuation of Equilibrium-Line Altitude [18]. The results obtained in this study were compared to the glaciological measurements from 1980 [18]. In 1980-2014, all the glaciers under study saw an increase in ELA, which averaged 78 m. Considering that the last satellite image available for 2014 was from the end of August, i.e., before the end of the melting period, the value is likely to be underestimated, even though the melt rate in September was rather low (see: Figure 3). Based on the final results of snow balance modelling, Uszczyk et al. [57] calculated a 71.5% SCA reduction on Hansbreen at the real end of melting season in 2014. The mentioned underestimation equalled c. 6.1% and might adequately represent other glaciers under study.
Remote sensing data supported by direct glaciological measurements provide valuable material, giving a better insight into the general spatial patterns of the melting processes on Svalbard glaciers, e.g., [15][16][17]58,59]. Attempts were made to use common automated snow identification methods, such as the Normalised Difference Snow Index, NDSI [36,60] and the Normalised Difference Snow and Ice Index, NDSII [61,62] as part of the work in order to compare the data with the results of the supervised classification. However, the attempts failed. In accordance with the studies cited above, the threshold value for snow cover detection is considered to be 0.4 (low probability of snow). Below this value, many non-snow pixels are identified as snow [60]. The indices calculated in the present study showed much higher values, also for ice-covered areas (>0.95). In our opinion, both snow indices display serious limitations when it comes to estimation of the melting characteristics of the snow cover present on glaciers because they fail to identify the transition between snow and glacier ice. A similar problem was raised by [63]. Mapping the extent of the temporary snow line and its position relative to the advance of melting was also abandoned. The highly non-continuous character of the snow line requires far-fetched generalisation and could lead to a significant error.

Conclusions
The year 2014 was the second warmest year since 1978. The low precipitation recorded in winter months and the high standard deviation from monthly mean air temperature resulted in reduced snow depth, which was observable in particular in the accumulation zones of the Hornsund glaciers. A noteworthy phenomenon was the presence of deep thermal inversions in the summer. Their share of the overall number of inversions during the year was 80%. July, which was the warmest month, saw the highest frequency of the inversions, adding to the intensive melting processes in the upper areas of the glaciers.
One of the factors slowing down the meltwater percolation and runoff processes is the internal structure of the snow cover, in particular the presence of very compact layers, such as Melt-Freeze Crusts and Ice Formations. In recent years, the winter months have seen a growing frequency of warm spells, with rain-on-snow events, which reach even highly elevated glacial zones, resulting in the presence of the above layers in the snowpack. Their share of the overall snow depth has been high in recent seasons, increasing the mean snow density.
The melting vertical gradient on Hansbreen amounted to −0.34 m w.e. 100 m −1 , and the melt rate of the lowermost ablation zone was three times higher than in the upper parts of the glacier, causing significant surface lowering. The highest melt rate was observed in July, when most of the ablation zone had lost the snow cover which protects the surface of the glacier against increased absorption of short-wave solar radiation.
Landsat 8 data allowed for an effective analysis of the snow cover melting characteristics on Hornsund glaciers and the changes in their geometry, the mapping of the Equilibrium-Line Altitude, the evolution of supraglacial ponds, and the intensity of the outflow of meltwater with suspended organic matter into the fjord. The Snow-Covered Area, calculated for the entire study area, varied from 98% at the beginning of the melting season to 43% at the end of August, and was reduced considerably at a rate of −0.56% day −1 during the period under study. When comparing the individual glaciers, the total SCA reduction showed noticeable spatial diversity, which is related to local microclimatic and topographic conditions. Very intensive surface melting continuing over longer periods leads to gradual transformation of entire glacial systems. In the years 1980-2014, all the glaciers terminating in the Hornsund Fjord saw an increase in ELA, while as a result of glacier calving and terminus retreat, their geometry also changed. In 2004-2014, the surface area of the glaciers reduced overall, i.e., by −2.55%. The greatest changes in the period were observed on Mendelejevbreen, and the least pronounced on Paierlbreen.
The results indicate the future evolution of the natural environment of the Hornsund region.