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

Glacier facies of Vestfonna (Svalbard) based on SAR images and GPR measurements

The warming climate of the Arctic affects the mass budget of glaciers, and changes in the distribution of glacier facies are indicative of their response to climate change. The glacial mass budget over large land ice masses can be estimated by remote sensing techniques, but selecting an efficient remote sensing method for recognizing and mapping glacier facies in the Arctic remains a challenge. In this study, we compared several methods of distinguishing the facies of the Vestfonna ice cap, Svalbard, based upon Synthetic Aperture Radar (SAR) images and terrestrial high frequency Ground Penetrating Radar (GPR) measurements. Glacier zones as determined using the backscattering coefficient (sigma0) of SAR images were compared against GPR data, and an alternative application of Internal Reflection Energy (IRE) calculated from terrestrial GPR data was also used for differentiating the extent of glacier facies. The IRE coefficient was found to offer a suitable method for distinguishing glacier zones and for validating SAR analysis. Furthermore, results of analysis of fully polarimetric Phased Array type Lband Synthetic Aperture Radar (ALOS PALSAR) and European Remote Sensing Synthetic Aperture Radar (ERS-2 SAR) images were compared with the IRE coefficient classification. Especially promising method is H-α segmentation, where the glacier zone boundaries corresponded very well with both GPR visual interpretation and IRE classification results. The IRE coefficient's simplicity of calculation makes it a good alternative to the subjective GPR visual interpretation method, where results strongly depend on the operator's level of experience. We therefore recommend for GPR profiles to be used for additional validation of SAR image analysis in studies of glacier facies on the High Arctic ice masses

Newly identified climatically and environmentally significant high-latitude dust sources

Dust particles from high latitudes have a potentially large local, regional, and global significance to climate and the environment as short-lived climate forcers, air pollutants, and nutrient sources. Identifying the locations of local dust sources and their emission, transport, and deposition processes is important for understanding the multiple impacts of high-latitude dust (HLD) on the Earth's systems. Here, we identify, describe, and quantify the source intensity (SI) values, which show the potential of soil surfaces for dust emission scaled to values 0 to 1 concerning globally best productive sources, using the Global Sand and Dust Storms Source Base Map (G-SDS-SBM). This includes 64 HLD sources in our collection for the northern (Alaska, Canada, Denmark, Greenland, Iceland, Svalbard, Sweden, and Russia) and southern (Antarctica and Patagonia) high latitudes. Activity from most of these HLD sources shows seasonal character. It is estimated that high-latitude land areas with higher (SI ≥0.5), very high (SI ≥0.7), and the highest potential (SI ≥0.9) for dust emission cover >1 670 000 km2, >560 000 km2, and >240 000 km2, respectively. In the Arctic HLD region (≥60∘ N), land area with SI ≥0.5 is 5.5 % (1 035 059 km2), area with SI ≥0.7 is 2.3 % (440 804 km2), and area with SI ≥0.9 is 1.1 % (208 701 km2). Minimum SI values in the northern HLD region are about 3 orders of magnitude smaller, indicating that the dust sources of this region greatly depend on weather conditions. Our spatial dust source distribution analysis modeling results showed evidence supporting a northern HLD belt, defined as the area north of 50∘ N, with a “transitional HLD-source area” extending at latitudes 50–58∘ N in Eurasia and 50–55∘ N in Canada and a “cold HLD-source area” including areas north of 60∘ N in Eurasia and north of 58∘ N in Canada, with currently “no dust source” area between the HLD and low-latitude dust (LLD) dust belt, except for British Columbia. Using the global atmospheric transport model SILAM, we estimated that 1.0 % of the global dust emission originated from the high-latitude regions. About 57 % of the dust deposition in snow- and ice-covered Arctic regions was from HLD sources. In the southern HLD region, soil surface conditions are favorable for dust emission during the whole year. Climate change can cause a decrease in the duration of snow cover, retreat of glaciers, and an increase in drought, heatwave intensity, and frequency, leading to the increasing frequency of topsoil conditions favorable for dust emission, which increases the probability of dust storms. Our study provides a step forward to improve the representation of HLD in models and to monitor, quantify, and assess the environmental and climate significance of HLD

Newly identified climatically and environmentally significant high-latitude dust sources

Dust particles from high latitudes have a potentially large local, regional, and global significance to climate and the environment as short-lived climate forcers, air pollutants, and nutrient sources. Identifying the locations of local dust sources and their emission, transport, and deposition processes is important for understanding the multiple impacts of high-latitude dust (HLD) on the Earth\u27s systems. Here, we identify, describe, and quantify the source intensity (SI) values, which show the potential of soil surfaces for dust emission scaled to values 0 to 1 concerning globally best productive sources, using the Global Sand and Dust Storms Source Base Map (G-SDS-SBM). This includes 64 HLD sources in our collection for the northern (Alaska, Canada, Denmark, Greenland, Iceland, Svalbard, Sweden, and Russia) and southern (Antarctica and Patagonia) high latitudes. Activity from most of these HLD sources shows seasonal character. It is estimated that high-latitude land areas with higher (SI ≥0.5), very high (SI ≥0.7), and the highest potential (SI ≥0.9) for dust emission cover >1 670 000 km$^{2}$, >560 000 km$^{2}$, and >240 000 km$^{2}$, respectively. In the Arctic HLD region (≥60$^{∘}$ N), land area with SI ≥0.5 is 5.5 % (1 035 059 km$^{2}$), area with SI ≥0.7 is 2.3 % (440 804 km$^{2}$), and area with SI ≥0.9 is 1.1 % (208 701 km$^{2}$). Minimum SI values in the northern HLD region are about 3 orders of magnitude smaller, indicating that the dust sources of this region greatly depend on weather conditions. Our spatial dust source distribution analysis modeling results showed evidence supporting a northern HLD belt, defined as the area north of 50$^{∘}$ N, with a “transitional HLD-source area” extending at latitudes 50–58∘ N in Eurasia and 50–55$^{∘}$ N in Canada and a “cold HLD-source area” including areas north of 60$^{∘}$ N in Eurasia and north of 58$^{∘}$ N in Canada, with currently “no dust source” area between the HLD and low-latitude dust (LLD) dust belt, except for British Columbia. Using the global atmospheric transport model SILAM, we estimated that 1.0 % of the global dust emission originated from the high-latitude regions. About 57 % of the dust deposition in snow- and ice-covered Arctic regions was from HLD sources. In the southern HLD region, soil surface conditions are favorable for dust emission during the whole year. Climate change can cause a decrease in the duration of snow cover, retreat of glaciers, and an increase in drought, heatwave intensity, and frequency, leading to the increasing frequency of topsoil conditions favorable for dust emission, which increases the probability of dust storms. Our study provides a step forward to improve the representation of HLD in models and to monitor, quantify, and assess the environmental and climate significance of HLD

Newly identified climatically and environmentally significant high-latitude dust sources

Peer reviewe

Comparison of Three Methods for Distinguishing Glacier Zones Using Satellite SAR Data

Changes in glacier zones (e.g., firn, superimposed ice, ice) are good indicators of glacier response to climate change. There are few studies of glacier zone detection by SAR that are focused on more than one ice body and validated by terrestrial data. This study is unique in terms of the dataset collected—four C- and L-band quad-pol satellite SAR images, Ground Penetrating Radar data, shallow glacier cores—and the number of land ice bodies analyzed, namely, three tidewater glaciers in Svalbard and one ice cap in Iceland. The main aim is to assess how well popular methods of SAR analysis perform in distinguishing glacier zones, regardless of factors such as the morphologic differences of the ice bodies, or differences in SAR data. We test and validate three methods of glacier zone detection: (1) Gaussian Mixture Model–Expectation Maximization (GMM-EM) clustering of dual-pol backscattering coefficient (sigma0); (2) GMM-EM of quad-pol Pauli decomposition; and (3) quad-pol H/α Wishart segmentation. The main findings are that the unsupervised classification of both sigma0 and Pauli decomposition are promising methods for distinguishing glacier zones. The former performs better at detecting the firn zone on SAR images, and the latter in the superimposed ice zone. Additionally, C-band SAR data perform better than L-band at detecting firn, but the latter can potentially separate crevasses via the classification of sigma0 or Pauli decomposition. H/α Wishart segmentation resulted in inconsistent results across the tested cases and did not detect crevasses on L-band SAR data

Seeking the sources of dust: geochemical and magnetic studies on “cryodust” in glacial cores from Southern Spitsbergen (Svalbard, Norway)

Natural mineral particulate matter deposited from aerosols and trapped in glaciers— herein defined as “cryodust”—may be an excellent indicator of atmospheric circulation, if terrestrial sources of dust can be identified. In this study, we analyzed the composition of cryodust in shallow ice cores taken from five glaciers in Southern Spitsbergen (Svalbard Archipelago, Northern Norway). The chemical composition, magnetic properties and radiogenic ages of individual grains were measured, where possible, to provide indicators of source areas. To identify mineral and rock fragments, solid particulates were examined by Scanning Electron Microscope fitted with a backscattered electron and Energy Dispersive Spectroscopic detectors. An Electron MicroProbe was employed for the U-Th-Pb chemical dating of monazite grains. Magnetic measurements comprised analyses of magnetic susceptibility ( ) vs. temperature (T) variations and determination of magnetic hysteresis parameters. Monazite ages span 445–423 Ma, consistent with mineral growth during the Caledonian orogeny. Caledonian rocks are exposed in the Nordaustlandet area of North-Eastern Svalbard, and this is the most probable source for monazite grains. Magnetic analyses show a predominance of ferrous (FeII) over ferric (FeIII) phases, consistent with a lack of input from subtropical sources. The results from both methods are consistent with local sources of dust from exposures in the Svalbard archipelago

Changes of glacier facies on Hornsund glaciers (Svalbard) during the decade 2007–2017

Changes in glacier facies (glacier zones), such as firn or superimposed ice (SI), are good indicators of glacier response to climate change. They are especially important for fast-warming Svalbard, where only a few glaciers are under glaciological mass balance monitoring. This paper presents a first study of changes of glacier facies extent for three tidewater glaciers located in southern Spitsbergen (Svalbard) and it is based on both satellite remote sensing and terrestrial data analysis, covering two time spans: 2007–2017 for Hansbreen and 2012–2017 for Storbreen and Hornbreen. Satellite remote sensing analysis include unsupervised classification of Synthetic Aperture Radar (SAR) data from both decommissioned (ENVISAT ASAR) and modern satellite missions (RADARSAT-2, Sentinel-1). The results of the SAR classification are compared to the information on glacier zones retrieved from terrestrial data, i.e. shallow cores and visual interpretation of 800 MHz Ground Penetrating Radar (GPR) profiles. In addition, a novel application of the Internal Reflection Power (IRP) coefficient as an objective method of distinguishing glacier zones based on GPR data is discussed. Changes in glacier facies areas over time are analysed, as well as their correlation to Hansbreen's mass balance. The main finding of the study is that firn and SI of Hansbreen, Storbreen and Hornbreen significantly decreased over the study period. For example, due to continuous negative mass balance between 2010 and 2017, the contribution of firn area to Hansbreen's total area decreased ca. 14% (cumulative firn area loss during that time: ~45%) whereas since 2012 SI has not been distinguished as a vast area on this glacier. In addition, an east–west gradient of firn area loss was observed as a result of differences in local climate conditions. Therefore, for the common time span (i.e. 2012–2017) Hansbreen recorded a ca. 12% loss of firn contribution to glacier area whereas Hornbreen recorded ca. 9%. Finally, application of the IRP coefficient as an objective method of glacier zones discrimination by GPR data gave very good results, so the method is recommended for future analysis of glacier zones instead, or as a support, to popular visual interpretation of the GPR profiles

Quality Assessment and Glaciological Applications of Digital Elevation Models Derived from Space-Borne and Aerial Images over Two Tidewater Glaciers of Southern Spitsbergen

In this study, we assess the accuracy and precision of digital elevation models (DEM) retrieved from aerial photographs taken in 2011 and from Very High Resolution satellite images (WorldView-2 and Pléiades) from the period 2012−2017. Additionally, the accuracy of the freely available Strip product of ArcticDEM was verified. We use the DEMs to characterize geometry changes over Hansbreen and Hornbreen, two tidewater glaciers in southern Spitsbergen, Svalbard. The satellite-based DEMs from WorldView-2 and Pléiades stereo pairs were processed using the Rational Function Model (RFM) without and with one ground control point. The elevation quality of the DEMs over glacierized areas was validated with in situ data: static differential GPS survey of mass balance stakes and GPS kinematic data acquired during ground penetrating radar survey. Results demonstrate the usefulness of the analyzed sources of DEMs for estimation of the total geodetic mass balance of the Svalbard glaciers. DEM accuracy is sufficient to investigate glacier surface elevation changes above 1 m. Strips from the ArcticDEM are generally precise, but some of them showed gross errors and need to be handled with caution. The surface of Hansbreen and Hornbreen has been lowering in recent years. The average annual elevation changes for Hansbreen were more negative in the period 2015−2017 (−2.4 m a−1) than in the period 2011−2015 (−1.7 m a−1). The average annual elevation changes over the studied area of Hornbreen for the period 2012−2017 amounted to −1.6 m a−1. The geodetic mass balance for Hansbreen was more negative than the climatic mass balance estimated using the mass budget method, probably due to underestimation of the ice discharge. From 2011 to 2017, Hansbreen lost on average over 1% of its volume each year. Such a high rate of relative loss illustrates how fast these glaciers are responding to climate change