Next Article in Journal
A Decadal Global Climatology of Ice Cloud Fraction with Their Microphysical and Optical Properties Inferred from the CALIPSO and Reanalysis Data
Previous Article in Journal
HierarchyNet: Hierarchical CNN-Based Urban Building Classification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dark Glacier Surface of Greenland’s Largest Floating Tongue Governed by High Local Deposition of Dust

1
Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung, D-27560 Bremerhaven, Germany
2
Department of Geosciences, University of Bremen, D-28359 Bremen, Germany
3
Institute of Applied Mechanics, Department of Mechanical and Process Engineering, TU Kaiserslautern, D-67653 Kaiserslautern, Germany
4
Environmental Radioactivity Laboratory, Institute of Nuclear and Radiological Science and Technology, Energy and Safety, P.O. Box 60037, 15310 Agia Paraskevi, Greece
5
CNR-IIA Institute of Atmospheric Pollution Research, National Research Council of Italy, 50019 Sesto Fiorentino, Italy
6
CNR-ISP Institute of Polar Sciences, National Research Council of Italy, 00010 Montelibretti, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(22), 3793; https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223793
Submission received: 2 October 2020 / Revised: 9 November 2020 / Accepted: 14 November 2020 / Published: 19 November 2020
(This article belongs to the Section Environmental Remote Sensing)

Abstract

:
Surface melt, driven by atmospheric temperatures and albedo, is a strong contribution of mass loss of the Greenland Ice Sheet. In the past, black carbon, algae and other light-absorbing impurities were suggested to govern albedo in Greenland’s ablation zone. Here we combine optical (MODIS/Sentinel-2) and radar (Sentinel-1) remote sensing data with airborne radar and laser scanner data, and engage firn modelling to identify the governing factors leading to dark glacier surfaces in Northeast Greenland. After the drainage of supraglacial lakes, the former lake ground is a clean surface represented by a high reflectance in Sentinel-2 data and aerial photography. These bright spots move with the ice flow and darken by more than 20% over only two years. In contrast, sites further inland do not exhibit this effect. This finding suggests that local deposition of dust, rather than black carbon or cryoconite formation, is the governing factor of albedo of fast-moving outlet glaciers. This is in agreement with a previous field study in the area which finds the mineralogical composition and grain size of the dust comparable with that of the surrounding soils.

Graphical Abstract

1. Introduction

Surface melt of the Greenland Ice Sheet contributes to about 60% of the mass loss of the ice sheet in the past decade [1,2]. Strong surface melt appears naturally in the lower elevated zones around the margins of the ice sheet, that are prone to higher summer air temperatures and lower albedo. Albedo itself is driven by light absorbing impurities and grain size. This effect has been important for the paleo-evolution of the Greenland Ice Sheet [3] and will remain a governing factor in future [4].
The surface of the ice sheet may consist of several constituents beside snow, firn or bare ice of different grain size: empty or water filled crevasses and thin fractures [5], surface melt water lakes [6], rivers [7], saturated patches [8] and impurities. Impurities may range from insoluble light absorbing impurities (LIA) as black carbon (BC, soot), dust arising from local sources or remote deserts, algae, organic matter and fine debris, e.g., sedimentary rock particles from nearby mountains. Dust is also used to date ice sheet retreat and advance from ice cores [9]. While some impurities are directly connected to emissions from human activities, others may only be indirectly affected by anthropogenic effects, like changes in wind fields, that affect erosion and transport of debris from nearby mountains.
Ryan et al. [10] conducted a UAV-borne survey over a prominent dark zone of in the western part of the Greenland Ice Sheet. They found that 73% of the variation in albedo can be attributed to three constituents of the surface: dust, black carbon and pigmented algae. The effect of cryoconite in this area is negligible as it only accounts for <1%.
The optical behavior of the ice-sheet components affect strongly the evolution of the melting season since while snow and firn are capable of reflecting the short-wave infrared radiations, water and ice are absorbing terms in such a wavelength domain [11]. While the visible (VIS) albedo ( λ < 0.7 μ m) is above 0.9 and 0.7 for snow/firn and water/ice, respectively, the near infrared (NIR) albedo ( λ > 0.7 μ m) is, in fact, below 0.1 for water and ice [11]. The occurrence of those components on the surface, combined with the meteo-climatological conditions of the season, control the evolution of the melting season, that can be accelerated by the removal of the uppermost snow cover. The importance of this process is governed by the duration of the snow cover and the distance of the snow line from the ice margin [12]. Inter-annual variations are controlled by the screen effect produced by the snow cover that is more reflective than ice cover, which are more sensitive to light absorption and consequently to heating and melting. Furthermore, Shimada et al. [13] focused the attention on the different ice types outcropping during the melting season, evidencing how the presence of bare-ice covers are correlated to the air temperature, and consequently to the capacity of melting on reaching ice layers characterized by higher impurity contents. Several authors [11,14,15] parameterised the snow albedo as a function of snow physical parameters and solar illumination. The snow physical parameters that govern albedo are snow grain size and shape, snow impurity concentrations, liquid water content, and their layer structure in the snowpack, together with snow depth. When the snow depth is optically thin, the underlying surface albedo also affects the snow albedo. While the near-infrared albedo strongly depends on snow grain size [11], the visible albedo is more dependent on the concentrations of snow impurities such as black carbon (BC) and mineral dust [11]. The growth of grains in wet sintering is reducing the albedo dominantly in the NIR, as NIR reflectance is governed by grain size [16]. This is suspected to be a self enhancing effect with more absorption leading to more albedo reduction [4].
Black carbon is a major agent of climate forcing through its role in changing snow/ice albedo [17]. There are several processes that change the optical properties of aerosol BC after deposition in the snow and metamorphism within the snow pack [18]. The long-term transport of BC aerosol over the last decades has been document at several Arctic sites [19,20,21]. However, our understanding of the long term trend of BC transport and subsequent deposition across the Arctic appears fragmented.
Recently, a study has found 10 year long trends to exhibit no statistically significant change for several areas including Svalbard (Ny-Ålesund), Greenland (Summit) and a slight increasing trend at Alert [22]. Despite global BC may have increased from 1990 to 2010, Arctic BC decreased dramatically over that same time span [19,23,24]. At Alert (Canada) the decrease from 1989–2003 reaches about 50% [24]. The decline of Arctic BC has been attributed to the dissolution of the Soviet Union, resulting in the closing of emission sources in Russia and Eastern Europe, and depopulation of the Russian Arctic [25]. Estimated BC deposition in the Arctic varies between 8 and more than 13 TgC yr 1 and is highly dependent on aerosol size containing the BC and the mixing state [26] with dry deposition contributing up to 40% of the total removal [27].
Direct BC deposition measurements and documented behaviour in the snow pack are still limited despite previous efforts to include microphysical properties and chemistry related to aging in order to provide a successful simulation of BC concentrations deposited over the Arctic [25,27]. In addition, the few available measurements of the dry deposition of submicron particles over snow display a large variability (0.02 to 0.33 cm s 1 [28,29]). For areas with a high ratio of dust to BC, the scavenging effect will enforce an increasingly dominance of dust over BC [30]. This is what we presume to take place in our study area.
Nearby our study site, a field experiment aiming at the characterization of ice types, albedo and impurity content was conducted in 2006 [31]. The team of Bøggild et al. [31] found, that the mineralogical composition of most impurities matches that of the surrounding soils, suggesting that most of the material is derived from local sources. If equally distributed, these impurities are found to reduce the albedo of Holocene ice to 0.2–0.4. Interestingly, ice with more impurities per unit area, but aggregated in cryoconite holes, reaches an albedo of 0.6 as the impurities are shaded by the walls of the hole.
Shimada et al. [13] investigated a prominent dark area of the Greenland Ice Sheet by means of optical remote sensing using the MODIS sensors. They introduced a distinction between so called bare ice and dark ice using the reflectance at different wavelengths. Bare ice is defined at λ = 860 nm (band 2) with a reflectance of less than 0.6, while dark is defined at λ = 660 nm (band 1) with a reflectance of less than 0.4. Shimada et al. [13] find a maximum of 1.5% of the ice sheet to consist of dark ice. In the following we are referring to this classification, using Sentinel-2 bands of similar wavelength.
It has also been found that wind-blown sedimentary rock particles that form a layer of debris are the dominant trigger for surface melt ponds on the George VI Ice Shelf in Antarctica [32]. The setting of George VI Ice Shelf with the nearby mountains of Alexander Land is similar to our study area at the margin of the North-East Greenland Ice Stream (NEGIS) and in particular the lower range of 79 Glacier (79 NG).
In this study, we take advantage of changes of drainage of supraglacial lakes leaving a bright surface and measure changes in reflectance over few years and set this into the context of surface properties along a fast flowing outlet glacier and a nearby slow moving area with exposed ice age ice.

2. Study Site

In this study we focus on the outlet region of the North East Greenland Ice Stream (NEGIS, Figure 1). This region is of high importance as 11% of the Greenland Ice Sheet is drained through the two main outlet glaciers Nioghalvfjerdsbræ (79 NG in the following) and Zachariæ Isstrøm [33]. Next to these fast flowing outlet glaciers which reach maximum velocities of >2 km a 1 moderate velocities are prevailing allowing older ice to reach the surface of the ice sheet. The study region is further characterized by many supraglacial lakes which fill and drain frequently [5]. For a detailed analysis, we have selected three flow lines that encompasses the accumulation to the ablation zone as shown in Figure 1. (i) a profile along a slow moving area that exhibits Holocene and Pleistocene ice at the surface allowing for a comparison with Bøggild et al. [31], (ii) a profile along the northern branch of 79 NG, following an airborne survey crossing supraglacial lakes and (iii) a profile along the southern branch of 79 NG including a supraglacial lake that has annually drained for some time, allowing to assess the change in reflectance due to in-situ deposition of impurities. In the following we address these profiles with the names profile “Old Ice”, “Northern Branch” and “Southern Branch” respectively. Figure 1 shows the study area including the locations of all three profiles.

3. Materials and Methods

3.1. Optical Remote Sensing: MODIS

The evolution of the snow cover on the glaciated and not-glaciated areas has been described using high-time resolved optical data based on the Moderate Resolution Imaging Spectroradiometer (MODIS) aboard on the Terra & Aqua platforms. While the M.D10A.006 product can provide the Terra/Aqua Snow Cover Daily Global 500 m [34,35], the M.D09GA.006 can support the estimation of the Terra/Aqua Surface Reflectance Daily Global 1 km and 500 m [36,37]. The first product is characterized by the information concerning the cloud cover for each pixel, useful for performing cloud masking, and the estimation of the Normalized Difference Snow Index (NDSI), suitable for discriminating snow-covered pixels from other types of surfaces (especially soil and rocks). Since the two platforms provide two overpass per day, it is possible to have a daily revisit of the study site even if the spatial resolution is coarser than other platforms. The second product provides the daily spectral reflectance, atmospherically corrected, for 7 bands ranging from 400 nm to 2150 nm. We considered for each pixel a time-series of 246 layers (two layers per day) having 123 days every melting season: from the Julian day 122 (1st May) to 246 (31st August). All of the cloud-covered pixels can be discarded and a daily check between the two available platforms can confirm the final estimation both in terms of NDSI and spectral reflectance at 500 m resolution. The analysis has been focused on characterizing firstly the different melting seasons, from 2016 to 2018, looking at the Snow Cover Area (SCA) on bare soil/rock areas. Similarly the analysis has been concerned on estimating the occurrence of snow cover on glaciated areas, estimating the first day without snow, the “First Snow-Free Day” (FSFD) and the lowest reflectance value at 660 nm (band 1). The first discrimination, between snow/firn and ice surfaces, has been approached following the threshold proposed by Shimada et al. [13], based on reflectance at 860 nm (band 2) below 0.6. The estimation of the lowest value at 660 nm has been based on defining for each pixel the time window where snow cannot be detected and calculating the reflectance median in that time interval. The importance of λ 660 derives from Shimada et al. [13] that proposed a threshold of 0.4 for discriminating between “clean” ice and “dirty” ice.

3.2. Optical Remote Sensing: Sentinel-2

The detailed analysis with high-resolution data has been approached using Sentinel-2 images. This satellite constellation is composed by two platforms that are active since 2015 and 2017 in the framework of the European Copernicus EO Programme [38]. The Copernicus Open Access Hub (https://scihub.copernicus.eu/) provides L1C data that must be converted from top-of-atmosphere to bottom-of-atmosphere reflectance. The analysis has been developed selecting cloud-free images and processing those data with the Sen2Cor toolbox [39]. The atmospheric correction process required some input data that we retrieved from the AERONET network Giles et al. [40]. Considering the available values of the 500 nm aerosol optical depth (<0.1) at Thule and Ittoqqortoormiit (since there are no available station in the study area), we used as input data a visibility of 40 km. The total columnar ozone content was selected referencing to Bahramvash Shams et al. [41] and we used as input values 400 and 300 DU, respectively for April 2017 and July 2017, The aerosol type was selected to be maritime and the atmosphere type as mid latitude winter, following the indications of König M and N [42]. The preparation of a L2A final product, with 10 m resolution and surface reflectance, has provided information about the spectral optical behaviour (12 bands with wavelengths ranging from 440 nm to 2200 nm) of the surface at selected dates characterized by a low cloud coverage. The Sentinel-2 bands are distributed as follows: B1 at 443 nm, B2 at 490 nm, B3 at 560 nm, B4 at 1665 nm, B5 at 705 nm, B6 at 740 nm, B7 at 783 nm, B8 at 842 nm, B9 at 945 nm, B11 at 1614 nm, B12 at 2202 nm. We defined for the next analysis two normalized difference indexes useful for enhancing the different contribution of visible, near-infrared and short-wave infrared wavelength domains.
I n d e x 1 = ( B 3 B 8 ) / ( B 3 + B 8 ) I n d e x 2 = ( B 8 B 11 ) / ( B 8 + B 11 )

3.3. Radar Remote Sensing: Sentinel-1 Polarimetry

Dual-polarisation SAR radar imagery can be used to distinguish between different types of reflection. A flat surface (as, e.g., a water body with no significant waves) causes a specular reflection, which means that the side-looking radar system hardly receives any return power at all. On a rough surface the reflection is dominated by Bragg-like scattering [43] and the incident horizontal polarisation is mainly reflected back as a horizontal polarised wave (HH). In contrast, volume scattering (e.g., in the firn pack) causes many random reflections within the volume and a part of the incident horizontally polarised wave is reflected in vertical polarisation (HV). The change in polarisation can therefore be used to distinguish between the scattering mechanism, that can subsequently be used for inferring material properties. For this study we are only using Sentinel-1 scenes in HH and HV polarisation for two different time periods, one representing frozen conditions and one representing summer, melt conditions.
The Single Look Complex (SLC) S1 scenes (IW mode) are processed using GAMMA Software [44] with the following steps: (1) calibration for instrumental noise, (2) mosaicing, (3) multilooking for speckle reduction and (4) georeferencing the scenes to EPSG:3413. The resulting products has a ground pixel resolution of 50 × 50 m. Topographic effects in the radar backscatter are corrected with a line-of-sight area normalization that transfers σ 0 to γ 0 using the ArcticDEM [45]. The details of polarimetric processing is presented in [5].
We selected two periods for S1 acquisition, one close to the acquisition of our airborne data in April 2018 and one during the melt season.

3.4. Roughness

To estimate surface roughness along the profiles we analyzed laser scanner data conducted with a RIEGL-LMSQ560 during the airborne survey in April 2018. The geo-located laser scanner point cloud (Latitude, Longitude, Elevation referenced to WGS84) was obtained using laser range and angle measurements, post processed GNSS trajectory (PPP solution using Waypoint 8.4) and attitude information (Inertial Measuring Unit: Honeywell Laseref V) within our processing scheme, following Hofton et al. [46].
Surface roughness, here defined as standard deviation of a high pass filtered laser scanner elevation model (ALS-DEM), was estimated at 1 km along track spacing. A series of 20% overlapping ALS-DEMs, each with a pixel resolution of 1 m and dimension of 1200 m × 1200 m, were conducted from the geo-located laser scanner point cloud data using an inverse distance weighting interpolation scheme. The final ALS-DEMs were referenced to the EPSG:3413 coordinate system.
High pass filtering (subtracting a smoothed version of the ALS-DEM using a 300 m kernel) was applied to remove large scale topographic features, as we were only interested in small scale roughness in the order of meters (e.g., sastrugis).

3.5. Firn Modelling

To simulate the density and age of the firn at the ice sheet surface, we relate spatial points along various flow lines, including the three selected sites, to points in time. This is done using the present time surface velocity field. For every flowline a one dimensional firn densification model is forced over time with corresponding mean annual values of the surface temperature and accumulation rate. The data is obtained from the regional climate model RACMO [47].
The model follows an approach of continuum mechanics solving conservation of mass, linear momentum and energy. The needed constitutive equations are derived from sintering mechanics, as the densification of firn can be compared to the process of hot isostatic pressing, a common process in sintering. This physics based approach of firn densification modelling provides the possibility to simulate not only the density of the firn profile but also its age, temperature, grain size and various other properties [48].
The implementation of the model allows to simulate varying accumulation and ablation rates at high temporal resolution. For accumulation a constant density of ρ 0 = 360 kg m 3 is chosen to represent recent deposited firn. In case of a negative accumulation rate the ablated mass is removed from the top of the simulated firn profile, uncovering older and more dense firn.

4. Results

4.1. Annual Deposition of Impurities

At many locations of 79 NG we find the glacier surface to be clean after the drainage of supraglacial lakes. Some examples are shown in Figure 2. For quantifying the change in reflectance over time, we found unfortunately only one appropriate lake, due to constraints such as S2 launch, cloud free coverage and snow free coverage with frozen surface. This lake is shown in Figure 3. We named this lake ‘Lake iCUPE’. Lake iCUPE has drained annually in the past years on: 21 August 2015, 17 July 2016, 16 July 2017, 1 August 2018, 17 July 2019 and 22 June 2020. After the drainage the surface appears bright, while over time this bright spots darkens, as can be seen in the right panel of Figure 3. As in 2019 major cracks further south appeared, which reached into our lake ground, the undisturbed time period ends 2018. In 2015 there has been no S2 coverage, thus we end up with the time period 2016–2018 to quantify the change in reflectance.
From a Sentinel-2 scene acquired on 28 July 2017 three brighter spots were digitized manually and the timing of their center positions was quantified from a flowline analysis (Figure 3). By projecting the locations of the bright spots back in time we found the centers of the two former lake grounds to be at the actual lake position on 11 August 2015 and 24 July 2016. As the glacier moves with a speed of about 1000 m/a, the lake ground (∼800–900 m length) of year 2 is separated from the year 1’s lake ground surface by about 100–150 m. A comparison of the reflectance of the lake ground over four subsequent patches of former lake ground indicates the in-situ darkening by deposition: the reflectance drops within the first year after drainage from 0.8 to 0.68 and further to 0.62 two years after drainage (Figure 3). After that, the former lake ground cannot be distinguished from surrounding areas. Here we assume that snow albedo is connected to reflectance as analogous but we are only making qualitative use of the terms.
We infer that the change in reflectance is due to local deposition of dust, as we can rule out lag-deposit to govern this effect: a comparison with a further upstream situated lake, Lake CE (Figure 1), shows that in that area, changes of reflectance with distance (time) is not significant (see Figure S1). Ablation rates and flow speed together amounts to roughly the same amount of melt at the surface, leaving to similar amounts of lag-deposit. This is in conflict to the strong change in reflectance at the floating tongue (Lake iCUPE), that can only be achieved by deposition in addition to lag-deposits. We hypothesize that this deposition is caused by surrounding soil and mountains (despite that we have not analysed any samples), as Bøggild et al. [31] have shown dust to come indeed from nearby mountains. Although Lake iCUPE is somewhat away from the former field site of the study of [31], it is very likely that at about the same elevation the source of the dust is similar. Fresh wet and dry deposition from non-local sources is anticipated to have orders of magnitude lower impact and is thus unlikely to cause such severe change in reflectance. Ref. [18] have concluded that annual deposited BC was causing a marginal reduction of the snow albedo, while post-depositional processes may cause a stronger direct and indirect impact on the snow albedo via enhanced metamorphism processes but in any case it has several order of magnitude lower impact than Holocene ice.
In addition to that, grain-size increase due to fluid sintering is similarly unlikely to be the cause, as the time scale for grain-growth is by far larger, consequently such a change in reflectance due to change in specific surface area would take considerably more time than just one year. Having ruled out these two causes only a local source is plausible. In order to check if this is changing with altitude, we analyse a second lake location further upstream in an area where lakes drain, too.
In order to broaden this, we are analysing profiles along flow lines of the glacier, that are shown in Figure 1. One profile is along the southern branch, crossing Lake iCUPE, another profile in the northern branch of 79 NG, and a third one is ending at the field site of the study of Bøggild et al. [31]. We start with the first, which we name ?Southern Branch?, picking up the findings of Bøggild et al. [31], who found Holocene and Pleistocene ice at the glacier surface.

4.2. Profile “Southern Branch”

The last along flow profile addressed as “Southern Branch” and situated furthest south (see Figure 1) is dedicated to cross “Lake iCUPE”. This supraglacial lake is used to study the effect of in-situ deposition of impurities in detail as described in the following section. Figure 4 is showing a similar figure than for profiles “Old Ice” and “Northern Branch”. Also there the transition into an impermeable surface layer (pore-close off density is ∼830 kg/m 3 ) is at an elevation of about 900 to 1000 m a.s.l. The backscatter shows more variation in the upstream area in this profile. Again the crevassed area is coincident with higher backscatter, both in summer and winter conditions. Similar to profile “Northern Branch”, the downstream part shows significant lower reflectance values, indicating a surface with high impurity content (Figure 3).
The average for the areas with elevation above 1000 m a.s.l. (distance from the ice margin > 20 km) showed on April about 0.99, 1.00 and 0.066 for bands 3, 8 and 11, respectively. The values in areas closer to the ice margin (<1000 m a.s.l.) were 0.95, 0.96 and 0.10 for the above mentioned bands and same season. The indexes were also in this case close to each other highlighting the occurrence of similar surface types (Index1: −0.017 and −0.004: Index2: 0.88 and 0.81). The same analysis evidences on July reflectance values significantly different: 0.94, 0.78 and 0.034 far from the ice margin; 0.78, 0.60 and 0.031 close to the ice margin. The Vis-NIR index (Index 1) highlighted limited differences (0.10 against 0.13) as well as the NIR-SWIR index (Index 2) (0.92 against 0.91).
In addition to these profiles we were using AWI’s snow radar (UWB-M) to detect reflections in the upper layers. Normally, summer layers are forming a stratigraphy that is also often used to determine accumulation rates over past decades to centuries [49]. Figure S4 displays radargrams along some short transects. Our analysis shows a decline in penetration depth with the transition into the impurity covered area. In the accumulation zone the stratigraphy can be resolved down to 13 m, with a decline to 1.5 m at the onset of the ablation zone. Further downstream no reflections of the last summer were found at all.

4.3. Profile “Northern Branch”

The data along profile “Northern Branch”, representing the fast moving part of 79 NG (shown in Figure 1), is presented in Figure 5. Firn modelling indicates a transition of a permeable to impermeable surface layer at around 120 km in an upstream located topographic plateau at about 1000 m a.s.l. (see also Figures S2 and S3). The profile crosses supraglacial lakes (grey bars in Figure 5), over which HH and HV backscatter drops, both in summer and winter conditions [5]. In winter the difference of HH and HV backscatter remains nearly constant down to the crevassed area, while in summer the difference increases towards lower elevations, indicating a decreasing volume scattering due to wet surface conditions. Starting from around 140 km surface roughness increases and becomes the dominating effect on radar backscatter and optical reflectance. With the onset of high roughness more of the radar signal is reflected back to the sensor. Moreover, a part of this increase in backscatter can also be explained by side reflections from sharp crevasses that are dominant in this area.
Most striking is the difference in backscatter between upstream and downstream areas indicating differences in surface characteristics: In the upstream (i.e., higher elevated) region, the backscatter is dominated by the ratio firn pack penetration while further downstream the crevassed surface roughness dominates the signal. Upstream in winter, the backscatter is nearly constant for elevations above 1000 m and decreases quite slowly in lower elevations. In contrast in summer, the backscatter drops very fast around 1500 m a.s.l., indicating a significantly wetter surface below.
Looking at reflectance values the average for the areas with elevation above 1000 m a.s.l. (distance from the ice margin > 20 km) showed on April about 0.97, 1.00 and 0.073 for bands 3, 8 and 11, respectively. The values in areas closer to the ice margin (<1000 m a.s.l.) were 0.93, 0.92 and 0.09 for the above mentioned bands and same season. The indexes were also in this case close to each other highlighting the occurrence of similar surface types (Index1: −0.016 and 0.010: Index2: 0.86 and 0.82). The same analysis evidences on July reflectance values significantly different: 0.89, 0.75 and 0.022 far from the ice margin; 0.71, 0.51 and 0.009 close to the ice margin. The Vis-NIR index (Index 1) highlighted a less significant difference (0.15 versus 0.09) as well as the NIR-SWIR index (Index 2) (0.94 versus 0.95).

4.4. Profile “Old Ice”

Figure 6 presents various data sets along the flow line of profile “Old Ice”, located in a slow moving area outside the main trunk of 79 NG (see also Figure 1). As the flowline is heading towards the area studied by Bøggild et al. [31] we mark in Figure 1) their study site, allowing for a direct comparison. Visual inspection of neither optical nor radar satellite imagery indicate any heavily crevassed region in this profile, nor lakes, but numerous rivers. Our firn modelling results indicate that the ablation zone is consisting predominantly of an ice surface with a density above pore-close off (see also Figures S2 and S3), blocking surface melt water to penetrate into the glacier, fostering surface run-off. The steep slope in this area leads is further intensifying run-off through rivers. In-situ deposition of impurities from nearby soil, as observed by [31] is likely to either remain on the surface or to be washed off with the run-off. Sentinel-2 data indicates a change from clean late Holocene ice to early Holocene ice and also the transition to Pleistocene ice, as also reported by Bøggild et al. [31]. With the onset of dirty Holocene ice, both, winter and summer periods show a smooth decline in HH and HV backscatter. Also the transition in reflectance of various bands in S2 is smooth, with a drop of about 0.1 in the visible range over the last 60 km towards the ice margin. In details, referring to areas in Figure 6, the average reflectance value for the areas with elevation above 1000 m a.s.l. (distance from the ice margin > 20 km) showed on April about 0.98, 1.00 and 0.075 for bands 3, 8 and 11, respectively. The values in areas closer to the ice margin (<1000 m a.s.l.) were 0.94, 0.95 and 0.08 for the above mentioned bands and same date. The indexes were both close to each other highlighting the occurrence of similar surface types (Index1: −0.012 and −0.007: Index2: 0.86 and 0.85). The same analysis evidences on July reflectance values significantly different: 0.89, 0.77 and 0.026 far from the ice margin; 0.79, 0.55 and 0.009 close to the ice margin. While the Vis-NIR index (Index 1) highlighted an important difference (0.18 against 0.07), the NIR-SWIR index (Index 2) showed limited increase (0.93 against 0.97).

4.5. Snow, Bare and Dark Ice Surface

The evolution of the snow cover can only be assessed with remote sensing data with a short repetition time, as provided by the MODIS platform. While this platform provides lower spatial resolution, it ensures higher time resolution that can be used to obtain a time series of spectral reflectance during each melting season. This approach was followed by Shimada et al. [13] with the aim of detecting bare ice and dark ice on the surface. We performed such an analysis in the 79 NG area for three consecutive years 2016, 2017 and 2018, focusing particularly on the relationship with the mountains surrounding the floating tongue of the glacier. The reduction of the fractional snow cover occurring on the outcropping rocks showed (Figure 7) that the snow cover started to melt by the end of May – beginning of June in 2016 and 2017 (respectively: 6 June 2016, 30 May 2017) and by the end of June in 2018. The areas closer to the outcropping rocks, ROI1 (120 pixels) and ROI4 (700 pixels) in Figure 1a, showed a change in spectral behavior. In the middle of June in 2016 and 2017 and at the beginning of July in 2018 reflectance values of up to 0.8 at λ 860 nm (characteristic for snow/firn surface) changed to the ones associated with a bare ice surface (<0.6 at λ 860 nm [13]). Both areas, ROI1 and ROI4, reached the lowest values of reflectance (close to 0.4 at λ 860 nm ) in the beginning of August in 2016 and 2017 and in the middle of August in 2018. In contrast, areas ROI2 (678 pixels) and ROI3 (920 pixels), that are situated further upstream and farther away from the mountains, showed consistently high reflectance values throughout the melting period, including a limited period of bare ice emergence in 2016 and 2017. The delayed melting season of 2018 preserved the snow cover further upstream. Looking at the spatial distribution of the different reflectance values, the above mentioned near-infrared band showed higher values in the inner part of the ice sheet with a distance from the outcropping rocks up to 40 km in 2016 and 2017, and 20 km in 2018. In order to identify areas in which dark ice recursively occurred during the various melting seasons, the red band ( λ 660 nm ) was examined with regard to its spatial distribution (Figure 8), since it represents the indicator for the detection of dark ice, as defined by Shimada et al. [13].

5. Discussion

We find that the backscatter in the two profiles in the main trunk of the glacier has its lowest values in the low elevated part in which also the optical reflectance is low. In addition to that, the drop of backscatter in wet, summer, conditions is in the order of 10 dB, whereas the drop in the areas with low reflectance is only about 5 dB (Figure 4 and Figure 5). This means that the water content does not affect the absorption in the same way in areas of high or low impurity content at the surface. This needs to be interpreted in relation to the change in backscatter from low to high impurity content of about 10 dB. Once the impurities at the surface reach a significant level, the backscatter is more dominated by impurity content than by moisture. This is further supported by the 2 m air temperature (see Figure S5), which has the same range at the lower elevated floating tongue than in the upper areas (km 70–140 in Figure 4 and Figure 5).
Non local sources annual rate of change is currently stabilized or only slightly increasing compared to past decades. For long term trends in BC, global emissions were estimated to grow by about 15% by 2010 compared to 1990 [50]. The aerosol light absorption levels do not show a significant trend in the Arctic as recently displayed in the analysis by Collaud Coen et al. [22] for long term data from Arctic stations Alert, Summit and Zeppelin-Ny Ålesund.
How will this develop in future? With rising temperatures the area that may undergo melt in future is likely to increase. Given that the floating tongue is experiencing now already strong melt during summer, we consider two scenarios: (i) with increasing melt the number of ponds on the floating tongue is increasing and/or the size of the ponds is increasing. In both cases the albedo becomes closer to that of open ocean. (ii) Increasing melt will lead to more surface water streams that flush the surface around the melt and sweep more dirt into the lakes. The consequence may be that the ice surface is coming to an increased albedo that is predominantly driven by coarsening of the upper layer. Also the area upstream the floating tongue that is currently prone to in-situ deposition from nearby mountains may undergo change. With increasing melt the cycle of lake filling/drainage can increase and lake size may increase, as well as rivers may increase in size and number. All that would lead to a more regular cleaning of the surface and consequently increase in albedo.
Taking the findings of Willeit and Ganopolski [3] into account, we recognise that there is a need of representation of the scavenging effect in firn hydrology modelling, which itself requires an evolution equation of grain size in wet conditions. Similar to that, using an ice sheet model tracer for dust concentration for accounting for albedo of early Holocene and Pleistocene ice in the ablations zones may become an important factor for projections of Greenland’s contribution to sea level rise.
Satellite remote sensing approaches to detect albedo in ablation zones may provide valuable information for assessing the influence over time span of decades and to validate simulated albedo. Approaches like the one by Kokhanovsky et al. [51] that employ a sophisticated atmospheric correction to Sentinel-3’s OLCI instrument obtaining an albedo estimation of high accuracy can potentially be combined with a debris cover ratio, as proposed by Azzoni et al. [52] allowing for separating in future better between debris/dust influence and grain size.

6. Conclusions

Melt in the ablation zone of fast moving outlet glaciers in Greenland that are situated in vicinity to surrounding mountains is an important factor in ice sheet mass loss. We find, that within two years time, the reflectance drops by about 25%, turning a clean ice surface after supraglacial lake drainage into a dirty surface. A change of this order of magnitude in such short time scale is unlikely to be facilitated by grain growth, black carbon or cryoconite formation, but points towards high local dust deposition. Given the magnitude of this effect we suggest that dust deposition and scavenging should be taken into account for simulation of albedo in regional climate models, despite being challenge to model adequately.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-4292/12/22/3793/s1, Figure S1 (Profile across Lake CE), Figure S2 (Simulated age at the surface), Figure S3 (Simulated density at the surface), Figure S4 (Stratigraphy from radargrams), Figure S5 (air temperature).

Author Contributions

Conceptualization, A.H., R.S. (Rosamaria Salvatori), R.S. (Roberto Salzano); PolSAR, L.S., N.N.; Roughness, V.H.; Optical Remote Sensing, R.S. (Rosamaria Salvatori), R.S. (Roberto Salzano); formal analysis; UWB: R.Z., V.H.; BC: K.E.; Firn modelling: T.S., R.M., A.H. All authors contributed to writing, discussed the results and prepared figures jointly. All authors have read and agreed to the published version of the manuscript.

Funding

This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 689443 via project iCUPE (Integrative and Comprehensive Understanding on Polar Environments). The airborne campaign RESURV79 was conducted with AWI’s Polar 6.

Acknowledgments

We want to thank the crew and technical team of Polar 6 for support for RESURV79. We acknowledge the support of Tobias Binder in UWB-M data acquisition.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
79 NG79 N Glacier
MODISModerate Resolution Imaging Spectroradiometer
S1, S2Sentinel-1, Sentinel-2
UWB-MUltra-Wide-Band-Microwave radar
BCblack carbon

References

  1. van den Broeke, M.R.; Enderlin, E.M.; Howat, I.M.; Kuipers Munneke, P.; Noël, B.P.Y.; van de Berg, W.J.; van Meijgaard, E.; Wouters, B. On the recent contribution of the Greenland ice sheet to sea level change. Cryosphere 2016, 10, 1933–1946. [Google Scholar] [CrossRef] [Green Version]
  2. Mouginot, J.; Rignot, E.; Bjørk, A.A.; van den Broeke, M.; Millan, R.; Morlighem, M.; Noël, B.; Scheuchl, B.; Wood, M. Forty-six years of Greenland Ice Sheet mass balance from 1972 to 2018. Proc. Natl. Acad. Sci. USA 2019, 116, 9239–9244. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Willeit, M.; Ganopolski, A. The importance of snow albedo for ice sheet evolution over the last glacial cycle. Clim. Past 2018, 14, 697–707. [Google Scholar] [CrossRef] [Green Version]
  4. Tedesco, M.; Doherty, S.; Fettweis, X.; Alexander, P.; Jeyaratnam, J.; Stroeve, J. The darkening of the Greenland ice sheet: Trends, drivers, and projections (1981–2100). Cryosphere 2016, 10, 477–496. [Google Scholar] [CrossRef] [Green Version]
  5. Schröder, L.; Neckel, N.; Zindler, R.; Humbert, A. Perennial supraglacial lakes in Northeast Greenland observed by polarimetric SAR. Remote Sens. 2020, 12, 2798. [Google Scholar] [CrossRef]
  6. Selmes, N.; Murray, T.; James, T.D. Fast draining lakes on the Greenland Ice Sheet. Geophys. Res. Lett. 2011, 38. [Google Scholar] [CrossRef]
  7. Smith, L.C.; Chu, V.W.; Yang, K.; Gleason, C.J.; Pitcher, L.H.; Rennermalm, A.K.; Legleiter, C.J.; Behar, A.E.; Overstreet, B.T.; Moustafa, S.E.; et al. Efficient meltwater drainage through supraglacial streams and rivers on the southwest Greenland ice sheet. Proc. Natl. Acad. Sci. USA 2015, 112, 1001–1006. [Google Scholar] [CrossRef] [Green Version]
  8. Yang, K.; Smith, L.C. Supraglacial Streams on the Greenland Ice Sheet Delineated From Combined Spectral–Shape Information in High-Resolution Satellite Imagery. IEEE Geosci. Remote Sens. Lett. 2013, 10, 801–805. [Google Scholar] [CrossRef]
  9. Simonsen, M.F.; Baccolo, G.; Blunier, T. East Greenland ice core dust record reveals timing of Greenland ice sheet advance and retreat. Nat. Commun. 2019, 10. [Google Scholar] [CrossRef] [Green Version]
  10. Ryan, J.; Hubbard, A.; Stibal, M.; Irvine-Fynn, T.; Cook, J.; Smith, L.C.; Cameron, K.; Box, J. Dark zone of the Greenland Ice Sheet controlled by distributed biologically-active impurities. Nat. Commun. 2018, 9. [Google Scholar] [CrossRef] [Green Version]
  11. Gardner, A.S.; Sharp, M.J. A review of snow and ice albedo and the development of a new physically based broadband albedo parameterization. J. Geophys. Res. Earth Surf. 2010, 115. [Google Scholar] [CrossRef]
  12. Ryan, J.; Smith, L.; van As, D.; Cooley, D.; Cooper, M.; Pitcher, L.; Hubbard, A. Greenland Ice Sheet surface melt amplified by snowline migration and bare ice exposure. Sci. Adv. 2019, 5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Shimada, R.; Takeuchi, N.; Aoki, T. Inter-Annual and Geographical Variations in the Extent of Bare Ice and Dark Ice on the Greenland Ice Sheet Derived from MODIS Satellite Images. Front. Earth Sci. 2016, 4, 1–10. [Google Scholar] [CrossRef] [Green Version]
  14. Aoki, T.; Kuchiki, K.; Niwano, M.; Kodama, Y.; Hosaka, M.; Tanaka, T. Physically based snow albedo model for calculating broadband albedos and the solar heating profile in snowpack for general circulation models. J. Geophys. Res. Atmos. 2011, 116. [Google Scholar] [CrossRef]
  15. Domine, F.; Salvatori, R.; Legagneux, L.; Salzano, R.; Fily, M.; Casacchia, R. Correlation between the specific surface area and the short wave infrared (SWIR) reflectance of snow. Cold Reg. Sci. Technol. 2006, 46, 60–68. [Google Scholar] [CrossRef]
  16. Warren, S.G. Optical properties of snow. Rev. Geophys. 1982, 20, 67–89. [Google Scholar] [CrossRef]
  17. He, C.; Flanner, M.G.; Chen, F.; Barlage, M.; Liou, K.N.; Kang, S.; Ming, J.; Qian, Y. Black carbon-induced snow albedo reduction over the Tibetan Plateau: Uncertainties from snow grain shape and aerosol–snow mixing state based on an updated SNICAR model. Atmos. Chem. Phys. 2018, 18, 11507–11527. [Google Scholar] [CrossRef] [Green Version]
  18. Jacobi, H.W.; Obleitner, F.; Da Costa, S.; Ginot, P.; Eleftheriadis, K.; Aas, W.; Zanatta, M. Deposition of ionic species and black carbon to the Arctic snowpack: Combining snow pit observations with modeling. Atmos. Chem. Phys. 2019, 19, 10361–10377. [Google Scholar] [CrossRef] [Green Version]
  19. Eleftheriadis, K.; Vratolis, S.; Nyeki, S. Aerosol black carbon in the European Arctic: Measurements at Zeppelin station, Ny-ålesund, Svalbard from 1998–2007. Geophys. Res. Lett. 2009, 36. [Google Scholar] [CrossRef] [Green Version]
  20. Stone, R.; Sharma, S.; Herber, A.; Eleftheriadis, K.; Nelson, D. A characterization of Arctic aerosols on the basis of aerosol optical depth and black carbon measurements. Elem. Sci. Anthr. 2014, 2. [Google Scholar] [CrossRef] [Green Version]
  21. Petäjä, T.; Duplissy, E.M.; Tabakova, K.; Schmale, J.; Altstädter, B.; Ancellet, G.; Arshinov, M.; Balin, Y.; Baltensperger, U.; Bange, J.; et al. Integrative and comprehensive Understanding on Polar Environments (iCUPE): The concept and initial results. Atmos. Chem. Phys. Discuss. 2020, 2020, 1–62. [Google Scholar] [CrossRef]
  22. Collaud Coen, M.; Andrews, E.; Alastuey, A.; Arsov, T.P.; Backman, J.; Brem, B.T.; Bukowiecki, N.; Couret, C.; Eleftheriadis, K.; Flentje, H.; et al. Multidecadal trend analysis of aerosol radiative properties at a global scale. Atmos. Chem. Phys. Discuss. 2020, 2020, 1–54. [Google Scholar] [CrossRef] [Green Version]
  23. Gong, S.L.; Zhao, T.L.; Sharma, S.; Toom-Sauntry, D.; Lavoué, D.; Zhang, X.B.; Leaitch, W.R.; Barrie, L.A. Identification of trends and interannual variability of sulfate and black carbon in the Canadian High Arctic: 1981–2007. J. Geophys. Res. Atmos. 2010, 115. [Google Scholar] [CrossRef]
  24. Sharma, S.; Andrews, E.; Barrie, L.A.; Ogren, J.A.; Lavoué, D. Variations and sources of the equivalent black carbon in the high Arctic revealed by long-term observations at Alert and Barrow: 1989–2003. J. Geophys. Res. Atmos. 2006, 111. [Google Scholar] [CrossRef]
  25. Sharma, S.; Ishizawa, M.; Chan, D.; Lavoué, D.; Andrews, E.; Eleftheriadis, K.; Maksyutov, S. 16-year simulation of Arctic black carbon: Transport, source contribution, and sensitivity analysis on deposition. J. Geophys. Res. Atmos. 2013, 118, 943–964. [Google Scholar] [CrossRef]
  26. Bond, T.C.; Doherty, S.J.; Fahey, D.W.; Forster, P.M.; Berntsen, T.; DeAngelo, B.J.; Flanner, M.G.; Ghan, S.; Kärcher, B.; Koch, D.; et al. Bounding the role of black carbon in the climate system: A scientific assessment. J. Geophys. Res. Atmos. 2013, 118, 5380–5552. [Google Scholar] [CrossRef]
  27. Liu, J.; Fan, S.; Horowitz, L.W.; Levy, H., II. Evaluation of factors controlling long-range transport of black carbon to the Arctic. J. Geophys. Res. Atmos. 2011, 116. [Google Scholar] [CrossRef]
  28. Contini, D.; Donateo, A.; Belosi, F.; Grasso, F.M.; Santachiara, G.; Prodi, F. Deposition velocity of ultrafine particles measured with the Eddy-Correlation Method over the Nansen Ice Sheet (Antarctica). J. Geophys. Res. Atmos. 2010, 115. [Google Scholar] [CrossRef] [Green Version]
  29. Grönlund, A.; Nilsson, D.; Koponen, I.K.; Virkkula, A.; Hansson, M.E. Aerosol dry deposition measured with eddy-covariance technique at Wasa and Aboa, DronningMaud Land, Antarctica. Ann. Glaciol. 2002, 35, 355–361. [Google Scholar] [CrossRef] [Green Version]
  30. Doherty, S.J.; Grenfell, T.C.; Forsström, S.; Hegg, D.L.; Brandt, R.E.; Warren, S.G. Observed vertical redistribution of black carbon and other insoluble light-absorbing particles in melting snow. J. Geophys. Res. Atmos. 2013, 118, 5553–5569. [Google Scholar] [CrossRef]
  31. Bøggild, C.E.; Brandt, R.E.; Brown, K.J.; Warren, S.G. The ablation zone in northeast Greenland: Ice types, albedos and impurities. J. Glaciol. 2010, 56, 101–113. [Google Scholar] [CrossRef] [Green Version]
  32. Wager, A.C. Flooding of the ice shelf in the GeorgVI sound. Br. Antarct. Surv. Bull. 1972, 28, 71–74. [Google Scholar]
  33. Krieger, L.; Floricioiu, D.; Neckel, N. Drainage basin delineation for outlet glaciers of Northeast Greenland based on Sentinel-1 ice velocities and TanDEM-X elevations. Remote Sens. Environ. 2020, 237, 111483. [Google Scholar] [CrossRef]
  34. Hall, D.; Riggs, G. MODIS/Terra Snow Cover Daily L3 Global 500m SIN Grid, Version 6. 2016. Available online: https://nsidc.org/data/MOD10A1/versions/6 (accessed on 2 October 2020). [CrossRef]
  35. Hall, D.; Riggs, G. MODIS/Aqua Snow Cover Daily L3 Global 500m SIN Grid, Version 6. 2016. Available online: https://ladsweb.modaps.eosdis.nasa.gov/missions-and-measurements/products/MOD09A1/ (accessed on 2 October 2020). [CrossRef]
  36. Vermote, E. MOD09A1 MODIS/Terra Surface Reflectance 8-Day L3 Global 500m SIN Grid V006. 2015. Available online: https://nsidc.org/data/MOD10A1/versions/6 (accessed on 2 October 2020). [CrossRef]
  37. Vermote, E. MYD09A1 MODIS/Aqua Surface Reflectance 8-Day L3 Global 500m SIN Grid V006. 2015. Available online: https://ladsweb.modaps.eosdis.nasa.gov/missions-and-measurements/products/MOD09A1/ (accessed on 2 October 2020). [CrossRef]
  38. Drusch, M.; Bello, U.; Carlier, S.; Colin, O.; Fernandez, V.; Gascon, F.; Hoersch, B. Sentinel-2: ESA’s Optical High-Resolution Mission for GMES Operational Services. Remote Sens. Environ. 2012, 120, 25–36. [Google Scholar] [CrossRef]
  39. Ranghetti, L.; Boschetti, M.; Nutini, F.; Busetto, L. sen2r: An R toolbox for automatically downloading and preprocessing Sentinel-2 satellite data. Comput. Geosci. 2020, 139, 104473. [Google Scholar] [CrossRef]
  40. Giles, D.M.; Sinyuk, A.; Sorokin, M.G.; Schafer, J.S.; Smirnov, A.; Slutsker, I.; Eck, T.F.; Holben, B.N.; Lewis, J.R.; Campbell, J.R.; et al. Advancements in the Aerosol Robotic Network (AERONET) Version 3 database—Automated near-real-time quality control algorithm with improved cloud screening for Sun photometer aerosol optical depth (AOD) measurements. Atmos. Meas. Tech. 2019, 12, 169–209. [Google Scholar] [CrossRef] [Green Version]
  41. Bahramvash Shams, S.; Walden, V.P.; Petropavlovskikh, I.; Tarasick, D.; Kivi, R.; Oltmans, S.; Johnson, B.; Cullis, P.; Sterling, C.W.; Thölix, L.; et al. Variations in the vertical profile of ozone at four high-latitude Arctic sites from 2005 to 2017. Atmos. Chem. Phys. 2019, 19, 9733–9751. [Google Scholar] [CrossRef] [Green Version]
  42. König, M.; Hieronymi, M.; Natascha, O. Application of Sentinel-2 MSI in Arctic Research: Evaluating the Performance of Atmospheric Correction Approaches Over Arctic Sea Ice. Front. Earth Sci. 2019, 7. [Google Scholar] [CrossRef] [Green Version]
  43. Cloude, S. The Dual Polarization Entropy/Alpha Decomposition: A PALSAR Case Study. In Proceedings of the 3rd International Workshop on Science and Applications of SAR Polarimetry and Polarimetric Interferometry; Lacoste, H., Ouwehand, L., Eds.; European Space Agency: Noordwijk, The Netherlands, 2007; ISBN 92-9291-208-7. [Google Scholar]
  44. Wegmüller, U.; Werner, C.; Strozzi, T.; Wiesmann, A.; Frey, O.; Santoro, M. Sentinel-1 Support in the GAMMA Software. Procedia Comput. Sci. 2016, 100, 1305–1312. [Google Scholar] [CrossRef] [Green Version]
  45. Porter, C.; Morin, P.; Howat, I.; Noh, M.J.; Bates, B.; Peterman, K.; Keesey, S.; Schlenk, M.; Gardiner, J.; Tomko, K.; et al. ArcticDEM. 2018. Available online: https://dataverse.harvard.edu/dataset.xhtml?persistentId=doi:10.7910/DVN/OHHUKH (accessed on 2 October 2020). [CrossRef]
  46. Hofton, M.A.; Blair, J.B.; Minster, J.B.; Ridgway, J.R.; Williams, N.P.; Bufton, J.L.; Rabine, D.L. An airborne scanning laser altimetry survey of Long Valley, California. Int. J. Remote Sens. 2000, 21, 2413–2437. [Google Scholar] [CrossRef]
  47. Noël, B.; Van De Berg, J.; Van Meijgaard, E.; Kuipers Munneke, P.; Van De Waal, R.S.W.; Van den Broeke, M.R. Evaluation of the updated regional comate model RACMO2.3: Summer snowfall impact on the Greenland Ice Sheet. Cryosphere 2015, 9, 1831–1844. [Google Scholar] [CrossRef] [Green Version]
  48. Schultz, T.; Müller, R.; Gross, D.; Humbert, A. Modelling the Transformation from Snow to Ice Based on the Underlying Sintering Process. Proc. Appl. Math. Mech. 2020, 20. accepted. [Google Scholar] [CrossRef]
  49. Karlsson, N.B.; Colgan, W.T.; Binder, D.; Machguth, H.; Abermann, J.; Hansen, K.; Pedersen, A.O. Ice-penetrating radar survey of the subsurface debris field at Camp Century, Greenland. Cold Reg. Sci. Technol. 2019, 165, 102788. [Google Scholar] [CrossRef]
  50. Klimont, Z.; Kupiainen, K.; Heyes, C.; Purohit, P.; Cofala, J.; Rafaj, P.; Borken-Kleefeld, J.; Schöpp, W. Global anthropogenic emissions of particulate matter including black carbon. Atmos. Chem. Phys. 2017, 17, 8681–8723. [Google Scholar] [CrossRef] [Green Version]
  51. Kokhanovsky, A.; Box, J.E.; Vandecrux, B.; Mankoff, K.D.; Lamare, M.; Smirnov, A.; Kern, M. The Determination of Snow Albedo from Satellite Measurements Using Fast Atmospheric Correction Technique. Remote Sens. 2020, 12, 234. [Google Scholar] [CrossRef] [Green Version]
  52. Azzoni, R.S.; Senese, A.; Zerboni, A.; Maugeri, M.; Smiraglia, C.; Diolaiuti, G.A. Estimating ice albedo from fine debris cover quantified by a semi-automatic method: The case study of Forni Glacier, Italian Alps. Cryosphere 2016, 10, 665–679. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Overview of the study area of 79 NG, Greenland. Background is a Sentinel-2 image from 28 July 2017. The red line denotes the equilibrium line altitude (ELA) taken from RACMO2 surface mass balance, the blue dotted line the simulated line of ice density at pore-close off. The orange boxes are the region of interest used in Figure 7 and the tent symbol near region 1 marks the location of the study of Bøggild et al. [31].
Figure 1. Overview of the study area of 79 NG, Greenland. Background is a Sentinel-2 image from 28 July 2017. The red line denotes the equilibrium line altitude (ELA) taken from RACMO2 surface mass balance, the blue dotted line the simulated line of ice density at pore-close off. The orange boxes are the region of interest used in Figure 7 and the tent symbol near region 1 marks the location of the study of Bøggild et al. [31].
Remotesensing 12 03793 g001
Figure 2. Former lake grounds (bright areas) in the study area as seen by Sentinel-2 on 28 July 2017 (a), 31 July 2016 (b), 5 August 2018 (c) and by aerial photography on 8 August 2013 (d). The white box marks the location of former Lake iCUPE from 2015 to 2017, for more details see Figure 3.
Figure 2. Former lake grounds (bright areas) in the study area as seen by Sentinel-2 on 28 July 2017 (a), 31 July 2016 (b), 5 August 2018 (c) and by aerial photography on 8 August 2013 (d). The white box marks the location of former Lake iCUPE from 2015 to 2017, for more details see Figure 3.
Remotesensing 12 03793 g002
Figure 3. (a) Spectral behaviour of three different lake floors formed in three different years. Coloured markers have been averaged considering the respective region of interests on Sentinel-2 BOA reflectance image (panel b, image acquired on 28 July 2017). Black lines are the spectral Vis-SWIR acquisitions performed at the surface by Bøggild et al. [31]: A—surface cryoconite; B—uniform dark; C—intermediate dark; D—light dark; E—Pleistocene ice; F—Holocene ice; G—snow. (c) sketch explaining the process from lake drainage to subsequent dirt deposition.
Figure 3. (a) Spectral behaviour of three different lake floors formed in three different years. Coloured markers have been averaged considering the respective region of interests on Sentinel-2 BOA reflectance image (panel b, image acquired on 28 July 2017). Black lines are the spectral Vis-SWIR acquisitions performed at the surface by Bøggild et al. [31]: A—surface cryoconite; B—uniform dark; C—intermediate dark; D—light dark; E—Pleistocene ice; F—Holocene ice; G—snow. (c) sketch explaining the process from lake drainage to subsequent dirt deposition.
Remotesensing 12 03793 g003
Figure 4. Profile “Southern Branch” displaying in the upper panel topography and density (left color bar) along the flow line (shown in Figure 1), as well as the small (thin line) and large scale roughness (thick line, coloured with amplitude HH-HV (16 August 2018) of the right color bar). HH and HV backscatter for winter and summer are displayed in the middle panel. Grey bards denote location of supraglacial lakes. The lower panel represents the reflectance of channel B3 (green), Index 1 (red) and Index 2 (grey) of S2 in summer (pale) and winter (dark). The glacier flow direction is from left (upstream) to right (downstream).
Figure 4. Profile “Southern Branch” displaying in the upper panel topography and density (left color bar) along the flow line (shown in Figure 1), as well as the small (thin line) and large scale roughness (thick line, coloured with amplitude HH-HV (16 August 2018) of the right color bar). HH and HV backscatter for winter and summer are displayed in the middle panel. Grey bards denote location of supraglacial lakes. The lower panel represents the reflectance of channel B3 (green), Index 1 (red) and Index 2 (grey) of S2 in summer (pale) and winter (dark). The glacier flow direction is from left (upstream) to right (downstream).
Remotesensing 12 03793 g004
Figure 5. Profile “Northern Branch” with the same panels as Figure 4.
Figure 5. Profile “Northern Branch” with the same panels as Figure 4.
Remotesensing 12 03793 g005
Figure 6. Profile “Old Ice” with nearly same panels as Figure 5, except that no laser scanner based roughness is available here. The glacier flow direction is from left (upstream) to right (downstream). The grey bars indicate here Holocene ice (wide bar) and Pleistocene ice (narrow bar), that is detected in S2 imagery following [31].
Figure 6. Profile “Old Ice” with nearly same panels as Figure 5, except that no laser scanner based roughness is available here. The glacier flow direction is from left (upstream) to right (downstream). The grey bars indicate here Holocene ice (wide bar) and Pleistocene ice (narrow bar), that is detected in S2 imagery following [31].
Remotesensing 12 03793 g006
Figure 7. Evolution of the snow cover area (SCA) over mountainous areas (black line and left y-axis) in relation to the reflectance variations at 860 nm in selected areas (coloured markers, right y-axis) numbered as in Figure 1a. The dashed line represents the reflectance threshold defined by [13] for discriminating between snow cover and ice. All the considered data were obtained analysing MODIS data as described in Section 3.1.
Figure 7. Evolution of the snow cover area (SCA) over mountainous areas (black line and left y-axis) in relation to the reflectance variations at 860 nm in selected areas (coloured markers, right y-axis) numbered as in Figure 1a. The dashed line represents the reflectance threshold defined by [13] for discriminating between snow cover and ice. All the considered data were obtained analysing MODIS data as described in Section 3.1.
Remotesensing 12 03793 g007
Figure 8. Evolution of the snowline (limit between persistent snow cover and ice) and of the extension of dark ice, defined by Shimada et al. [13] from 2016 to 2018.
Figure 8. Evolution of the snowline (limit between persistent snow cover and ice) and of the extension of dark ice, defined by Shimada et al. [13] from 2016 to 2018.
Remotesensing 12 03793 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Humbert, A.; Schröder, L.; Schultz, T.; Müller, R.; Neckel, N.; Helm, V.; Zindler, R.; Eleftheriadis, K.; Salzano, R.; Salvatori, R. Dark Glacier Surface of Greenland’s Largest Floating Tongue Governed by High Local Deposition of Dust. Remote Sens. 2020, 12, 3793. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223793

AMA Style

Humbert A, Schröder L, Schultz T, Müller R, Neckel N, Helm V, Zindler R, Eleftheriadis K, Salzano R, Salvatori R. Dark Glacier Surface of Greenland’s Largest Floating Tongue Governed by High Local Deposition of Dust. Remote Sensing. 2020; 12(22):3793. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223793

Chicago/Turabian Style

Humbert, Angelika, Ludwig Schröder, Timm Schultz, Ralf Müller, Niklas Neckel, Veit Helm, Robin Zindler, Konstantinos Eleftheriadis, Roberto Salzano, and Rosamaria Salvatori. 2020. "Dark Glacier Surface of Greenland’s Largest Floating Tongue Governed by High Local Deposition of Dust" Remote Sensing 12, no. 22: 3793. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223793

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop