Supraglacial meltwater accumulation on ice shelves may have important implications for future sea level rise. Despite recent progress in the understanding of Antarctic surface hydrology, potential influences on ice shelf stability as well as links to environmental drivers remain poorly constrained. In this study, we employ state-of-the-art machine learning on Sentinel-1 synthetic aperture radar (SAR) and optical Sentinel-2 satellite imagery to provide new insight into the inter-annual and intra-annual evolution of surface hydrological features across six major Antarctic Peninsula and East Antarctic ice shelves. For the first time, we produce a high-resolution record of supraglacial lake extent dynamics for the period 2015–2021 at unprecedented 10 m spatial resolution and bi-weekly temporal scale. Through synergetic use of optical and SAR data, we obtain a more complete mapping record also enabling the delineation of buried lakes. Our results for Antarctic Peninsula ice shelves reveal below-average meltwater ponding during most of melting seasons 2015–2018 and above-average meltwater ponding throughout summer 2019–2020 and early 2020–2021 considering years 2015–2021 as a reference period. Meltwater ponding on investigated East Antarctic ice shelves was far more variable, with above-average lake extents during most 2016–2019 melting seasons and below-average lake extents during 2020–2021, considering the reference interval 2016–2021. This study is the first to investigate relationships with climate drivers both spatially and temporally including time lag analysis. The results indicate that supraglacial lake formation in 2015–2021 is coupled to the complex interplay of local, regional and large-scale environmental drivers with similar driving factors over both ice sheet regions. In particular, varying air temperature, solar radiation and wind conditions influenced supraglacial lake formation over all six ice shelves despite strong local to regional discrepancies, as revealed through pixel-based correlation analysis. Furthermore, regional climatic conditions were shown to be influenced by Southern Hemisphere atmospheric modes showing large-scale impacts on the spatio-temporal evolution of supraglacial lakes as well as on above- or below-average meltwater ponding with respect to the period 2015–2021. Finally, the local glaciological setting, including melt–albedo feedbacks and the firn air content, was revealed to strongly influence supraglacial lake distribution. Recent increases in Antarctic Peninsula surface ponding point towards a further reduction in the firn air content, implying an increased risk for ponding and hydrofracture. In addition, lateral meltwater transport was observed over both Antarctic regions with similar implications for future ice shelf stability.
With accelerating global climate change, the Greenland and Antarctic ice
sheets are exposed to increasing environmental pressure (Meredith et al., 2019). Recent disintegration of Antarctic Peninsula (API) ice shelves was
associated with atmospheric warming causing extensive meltwater ponding and
hydrofracture to initiate ice shelf collapse (Banwell et al., 2013; Bell et al., 2018; Cook and Vaughan, 2010; Leeson et al., 2020; Rignot et al., 2004; Scambos et al., 2004). Ice shelf collapse leads to a reduction in buttressing, initiating glacier acceleration and causing API ice mass loss to increase from
Despite recent advances in the understanding of Antarctic surface hydrology, detailed analyses of the spatio-temporal evolution of Antarctic supraglacial lakes as well as implications for ice shelf stability and links to environmental controls remain poorly constrained. For improved evaluation of seasonal fluctuations in Antarctic surface hydrology as well as associated drivers and impacts, automated mapping methods with spatio-temporal transferability are required. In this context, recent progress in the mapping of Antarctic surface hydrology results from large-scale mapping efforts using optical (Dell et al., 2020; Dirscherl et al., 2020; Halberstadt et al., 2020; Kingslake et al., 2017; Moussavi et al., 2020; Stokes et al., 2019) and synthetic aperture radar (SAR) (Dirscherl et al., 2021) data from spaceborne remote sensing and has revealed widespread occurrence of surface ponding particularly on the API and East Antarctic Ice Sheet (EAIS) (Arthur et al., 2020b). Employing state-of-the-art machine learning method developments using Sentinel-1 SAR (Dirscherl et al., 2021) and optical Sentinel-2 data (Dirscherl et al., 2020) enables the generation of fused supraglacial lake extent classification products at high spatial and temporal resolution. This allows sensor-specific advantages and disadvantages to be exploited and overcome, including the retrieval of maps during polar darkness or cloud cover with Sentinel-1 SAR. Furthermore, the penetration capability of Sentinel-1 enables the mapping of buried and slightly frozen lakes. On the other hand, Sentinel-2 allows the retrieval of mappings for lakes that cannot be identified in SAR data due to effects such as severe wind roughening. Providing more complete mapping records than single-sensor classifications, combined Sentinel-1 and Sentinel-2 supraglacial lake extent mappings are the primary data source within the present study.
Exploiting the archive of Sentinel-1 and Sentinel-2, this study aims to (i) investigate seasonal and inter-annual supraglacial lake extent dynamics in 2015–2021 over six major Antarctic ice shelves on the API and EAIS, (ii) analyse anomalous supraglacial lake occurrence above or below average with respect to the observational period 2015–2021, (iii) establish a link between spatio-temporal variability in supraglacial lake occurrence and associated environmental controls including both the near-surface climate and the local glaciological setting, and (iv) evaluate the influence of teleconnections and Southern Hemisphere atmospheric modes on climate and lake extent variability. While it is beyond the scope of this study to identify specific short-term drivers of lake occurrence or absence, we mostly focus on the analysis of climate variability during the entire investigation period 2015–2021.
Overview map of Antarctica with the study regions on the Antarctic Peninsula (API)
The six major ice shelves investigated in this study are shown in Fig. 1. Of
these, three are situated on the south-western API, and three are distributed around the margin of the EAIS. The choice of study regions was made in agreement with known supraglacial lake occurrences across different
glaciological and climatic settings (Banwell et al., 2021; Dell et al., 2020; Dirscherl et al., 2021; Kingslake et al., 2017; Spergel et al., 2021) as well as satellite data availability. The study regions on the API are located in the Bellingshausen Sea sector and cover the neighbouring George VI, Bach and Wilkins ice shelves (Fig. 1a–c). Separated by Alexander Island (71
Sentinel-2 optical image extracts (RGB) showing supraglacial lake occurrence on
The three investigated ice shelves on the EAIS are Riiser-Larsen, Nivlisen
and Amery ice shelves, each stretching over an area of
In this study, the area of interest (AOI) chosen for statistical analyses over each ice shelf region was delineated in agreement with confirmed
supraglacial lake occurrences during at least one time step within our
investigation period of 2015–2021 as well as on the basis of resampled ERA5-Land
pixel footprints (turquoise area in Fig. 1a–f). For George VI, Bach, Wilkins, Riiser-Larsen, Nivlisen and Amery ice shelves, this resulted in AOI
subsets of 27 324, 9504, 15 660, 1512, 5292 and 41 760 km
For automated mapping of Antarctic supraglacial lake extents, we jointly use Sentinel-1 and Sentinel-2 satellite imagery. The Sentinel-1 SAR data used in this study were acquired in ground range detected (GRD) format and interferometric wide (IW) swath mode, delivering data products at 10 m pixel spacing and in horizontal transmit-and-receive polarization (HH). While the general repeat interval of Sentinel-1 corresponds to 6 d at the Equator, up to daily revisit frequencies can be achieved in high latitudes. For each study region, we acquired all available Sentinel-1 data covering the period 1 November to 31 March in 2015–2021 within pre-defined boundaries, including the six ice shelves as well as surrounding grounded area. Over Riiser-Larsen, Nivlisen and Amery ice shelves, data availability during the 2015–2016 melting season was not sufficient, restricting the Sentinel-1 data selection to the period 2016–2021. In contrast, Sentinel-1 data over George VI, Bach and Wilkins ice shelves were available during all melting seasons 2015–2021. Overall, this resulted in a total number of 3075 Sentinel-1 acquisitions covering all six study regions.
Regarding optical Sentinel-2 data, we selected all imagery with a scene cloud cover
For analysis of climatological controls on Antarctic supraglacial lake
occurrence in 2015–2021, we use hourly ERA5-Land atmospheric reanalysis data
(Munoz Sabater, 2019). ERA5-Land replaces the former ERA-Interim product of the European Centre for Medium-Range Weather Forecasts (ECMWF) and provides reanalysis data for a variety of land variables at a spatial resolution of 0.1
The ERA5-Land variables used in this study include 2 m air temperature, surface net solar radiation, wind, total precipitation and snowmelt. Data were retrieved for the period 1 October to 1 March of each melting season in 2015–2021. October and November data were included for evaluation of potential time lag effects of climate variables. For precipitation data, we retrieved hourly estimates for all remaining months in 2015–2020 in order to allow the evaluation of links between previous-year precipitation and supraglacial lake formation.
To start with, ERA5-Land temperature data provide the best estimate of the
seasonal cycle of near-surface temperature over coastal Antarctica and
particularly the API (Gossart et al., 2019; Tetzner et al., 2019). ERA5-Land temperature data outperform other reanalysis datasets considering its low bias and low mean absolute error of 2
The satellite data were pre-processed using sensor-specific processing tools. For Sentinel-1 SAR data, we used the open-source Sentinel Application Platform (SNAP). Pre-processing of Sentinel-1 SAR imagery includes orbit
correction, thermal noise removal, radiometric calibration, speckle filtering and terrain correction. For the latter, the TanDEM-X PolarDEM 90 m of Antarctica (Wessel et al., 2021) was used. In order to prepare the data for prediction of supraglacial lakes (Sect. 3.2.2), each pre-processed Sentinel-1
image was normalized and tiled into
Pre-processing of ERA5-Land data included a range of individual processing
steps for each variable. Hourly estimates of ERA5-Land air temperature were
resampled to daily averages, and temperature measured in kelvin was converted
to degrees Celsius (
For prediction of supraglacial lake extents in normalized Sentinel-1 image
patches, we use a convolutional neural network (CNN) based on residual U-Net
(Dirscherl et al., 2021). In detail, a standard U-Net architecture was complemented with residual connections as well as an Atrous Spatial Pyramid Pooling (ASPP) module for improved performance. In this context, the choice of a deep learning network over traditional image processing techniques was mainly made due to the similar appearance of Antarctic supraglacial lakes and surface features such as wet snow, blue ice or topographic shadow in single-polarized SAR imagery as well as due to effects such as wind roughening and speckle noise (Dirscherl et al., 2021). Model training was performed on a graphics processing unit (GPU) cluster using 21 200 Sentinel-1 image patches extracted from 58 manually labelled Sentinel-1 image scenes covering melting seasons 2018–2019 and 2019–2020 over 13 different Antarctic regions. Following semantic segmentation of the image patches, the resultant prediction probability maps were reshaped into their original image size and thresholded. In particular, all pixels above a
probability threshold of 0.5 were accepted as supraglacial lake pixels. As
part of post-processing, remaining errors in the binary classification maps
were reduced by applying morphological filtering as well as topographic and
coastline masking using a resampled version of the TanDEM-X PolarDEM as well
as a dilated Sentinel-1 coastline product of 2018 (Baumhoer et al., 2021). The performance of the Sentinel-1 classification approach was evaluated by means of 23 000 random point samples collected over 10 independent test sites that were not presented to the model during training and that contained surface features that are difficult to discriminate from lakes as well as lakes of varying appearances. Their evaluation as part of a confusion matrix suggests the good functionality of the algorithm with an average
Supraglacial lake classification in optical Sentinel-2 imagery was performed
using an advanced version of the machine learning algorithm published in
Dirscherl et al. (2020). In detail, the algorithm is based upon a random forest (RF) classifier trained on a high-performance computing (HPC) cluster with selected Sentinel-2 bands and indices covering 14 different Antarctic regions during three austral summers as well as manually created labels and topographic variables derived from the Antarctic TanDEM-X digital elevation model (DEM). Similarly to Sentinel-1, post-classification involved bands, indices, topographic and coastline masking, and morphological erosion. The main modifications in the advanced version of the RF algorithm include (i) the integration of more training data over shadow; (ii) the improvement of post-classification, e.g. through fine-tuning of variable thresholds; and (iii) the use of prediction probability maps for efficient slush elimination. Overall, this classification strategy returned an average
Following supraglacial lake extent mapping, we aggregated all data to bi-weekly time series. To start with, Sentinel-1 and Sentinel-2 classification products were combined through decision-level fusion, thus the mosaicking of all classifications to maximum lake extent classifications for each bi-weekly time interval between November and March of a melting season. To interpolate data gaps during specific months over Riiser-Larsen and Nivlisen ice shelves (Fig. S1d and e in the Supplement), we used the respective long-term mean of the period 2016–2021 at bi-weekly scale. Since data gaps were mostly present for months of low lake coverage, this interpolation strategy was considered acceptable. For analysis of environmental controls on Antarctic supraglacial lake occurrence, we converted the bi-weekly lake extent time series to fractional lake extent considering ERA5-Land grids as AOIs (see Fig. 1). Similarly, the daily ERA5-Land data were aggregated to average bi-weekly time series, and the grid values over the AOIs were extracted. For temperature, we also calculated bi-weekly maxima to capture temperature extreme events and positive degree days, resulting in seven climate variables used for statistical analysis. In order to investigate drivers of anomalous meltwater ponding and to analyse where above- or below-normal climatic conditions result in above- or below-normal lake formation with respect to the investigation period 2015–2021, we additionally converted the bi-weekly fractional lake extent and climate time series to time series of anomalies. In this study, anomalies are calculated with respect to the long-term mean 2015–2021 of each bi-weekly time interval. Finally, for joint analysis of January lake extents and total precipitation in a year preceding lake coverage, we extracted the corresponding aggregates and grid values over each study region. In this study, previous-year precipitation was used as an indirect indicator for the FAC of snow, thus for analysis of relationships between supraglacial lake formation and the state of the snowpack. Additionally, January lake extents were considered for analysis of links with annual Southern Annular Mode (SAM) (Marshall, 2003, 2018) as well as for evaluation of the inter-annual recurrence of lakes, both spatially and temporally supporting the investigation of local controls.
To establish a link between environmental variables and fractional lake extents, we use the Pearson correlation. Pearson's correlation coefficient (
As shown in Fig. 3a–c, supraglacial lake extents on the API varied strongly
both inter-annually and intra-annually. For George VI, Bach and Wilkins ice
shelves, supraglacial lake extents were highest during melting seasons 2019–2020 and 2020–2021, reaching peak values of
Supraglacial lake extent dynamics over
Intra-annual supraglacial lake extent anomalies over
Figure 5a–c illustrate the annual recurrence of lakes during January 2016–2021. As can be seen, supraglacial lakes tend to form at low elevations and low surface slopes downstream of the grounding line. On Bach and Wilkins ice shelves, lakes are clustered closer to the grounding line than on George VI Ice Shelf, where lakes accumulated in flow stripes over vast parts of the flat ice shelf during the last two melting seasons. Near Wilkins Ice Shelf, lakes extended furthest onto the adjacent grounded ice, where rock outcrop is present (Fig. 5c). Further, supraglacial lakes on George VI Ice Shelf are clustered closer to the grounding line and rock outcrop in years with low supraglacial lake coverage (see Fig. 5a and b), as reflected in higher recurrence times of lakes up to 5–6 years.
Inter-annual January recurrence of supraglacial lakes over
George VI
Supraglacial lake extent dynamics on the EAIS (Fig. 3d–f) differ from the
pattern observed on the API. Over all three ice shelves, comparatively high
lake coverage was observed during the 2016–2019 melting seasons, while
particularly low lake coverage was found for the 2020–2021 melting season. This pattern is also reflected in the calculation of supraglacial lake extent
anomalies for the overall period 2016–2021 (Fig. 4d–f). The highest
supraglacial lake extents over Riiser-Larsen (
The annual recurrence of lakes during January shows clustered supraglacial lake occurrence near the grounding line of Riiser-Larsen Ice Shelf, where a blue ice region is present (Figs. 2d and 5d); lakes spreading across Nivlisen Ice Shelf, with a progressive advance in a north-easterly direction (Figs. 5e and S3); and widespread lake occurrence on the southern and eastern Amery Ice Shelf (Fig. 5f). During years of peak lake occurrence, lakes expanded further onto the flat ice shelves, while years with lower lake coverage were characterized by lakes clustering closer to the grounding line, rock outcrop or blue ice (Figs. 1d–f, 2d, 5d–f and S3). Over all three ice shelves, lakes formed at low elevations and low surface slopes. At the same time, the January recurrence of lakes reached up to 4–5 years in local depressions upstream of the grounding line (Fig. 5d–f).
Figure 6 presents the results of the statistical correlation analysis between climate variables and fractional lake extents (see Tables S1 and S2 in the Supplement for individual correlation values). As can be seen, bi-weekly mean and maximum air temperature correlates positively at lag 0–1 over all three ice shelves considering both spatial averages (Fig. 6a) and anomalies of spatial averages (Fig. 6b). Moreover, lake extents over George VI Ice Shelf correlate positively with spatial averages of lag 2 solar radiation, and lake extents over Wilkins Ice Shelf correlate negatively with anomalies of lag 1 solar radiation and lag 1 wind direction. Lastly, surface melt correlates positively over Bach and Wilkins ice shelves considering spatial averages and anomalies.
Linear correlation analysis over EAIS ice shelves (Fig. 6a and b) returned significant positive lag 0–3 relationships with bi-weekly mean and maximum air temperature considering both spatial averages and anomalies. Likewise, solar radiation correlates positively at varying lags over all three ice shelves. Over Riiser-Larsen Ice Shelf, significant negative lag 0 correlation was found for precipitation and wind magnitude, and significant positive lag 0 correlation was found for wind direction considering spatial averages and anomalies. Further, statistical analysis over Nivlisen Ice Shelf returned negative lag 1 correlation for wind direction. Finally, surface melt correlates positively at lag 1 over Riiser-Larsen and Amery ice shelves.
Correlation with annual SAM (Table 1) indicates a strong negative linear
relationship over all three API ice shelves. Over EAIS ice shelves,
Results of the statistical correlation analysis between fractional
lake extents and previous-year precipitation as well as annual Southern
Annular Mode (SAM), with
Figures 7 and 8 illustrate selected results for pixel-based correlations, chosen on the basis of Fig. 6. The figures show only significant pixels with the highest correlations considering all performed cross-correlations as well as the corresponding lags. Remaining pixel-based correlations can be found in the Supplement (Figs. S4 and S5). Pixel-based correlation over George VI Ice Shelf (Fig. 7a–h) revealed positive relationships for temperature and solar radiation correlating at higher lags in the central ice shelf section compared to locations near the grounding line. Likewise, statistical analysis over Bach Ice Shelf (Fig. 7i–p) implies predominant positive correlations for temperature and radiation. Over both ice shelves, snowmelt correlates positively, and wind direction correlates positively in the east and negatively in the west. Spatial correlations over Wilkins Ice Shelf were rather heterogeneous (see Figs. 5c and 7q–x).
Pixel-based Pearson correlations (
Pixel-based Pearson correlations (
Over EAIS ice shelves (Fig. 8a–x), positive correlation was found for bi-weekly average air temperature and solar radiation. While lags were varying for each ice shelf, they were generally higher on the ice shelf area than on adjacent inland ice and near the grounding line. Moreover, correlation of lake extents with wind magnitude and precipitation over Riiser-Larsen Ice Shelf revealed a negative relationship at lag 0, and correlation with wind magnitude over Amery Ice Shelf revealed negative and positive relationships at varying lags along the south-eastern and eastern grounding zone.
Finally, Fig. 9 outlines the spatial correlation pattern between previous-year precipitation and January fractional lake extents. As correlations result from only few observations, potentially leading to unreliable values of
Pixel-based correlations (
Starting with local control factors, correlation analysis with precipitation revealed API lake coverage to be peaking in years with previous-year low precipitation (Table 1, Fig. 9a–c). This highlights that supraglacial lake formation is promoted over more compact, air-depleted firn (Arthur et al., 2020a; Leeson et al., 2020), while fresh snow dampens lake evolution allowing the percolation of meltwater into the low-density firn layer before it can accumulate on the surface (Kuipers Munneke et al., 2014; Lenaerts et al., 2017; Stokes et al., 2019). In fact, low precipitation in 2017, 2019 and 2020 might have led to increased firn air depletion, potentially contributing to facilitated meltwater ponding in 2017–2018, 2019–2020 and 2020–2021. Apart from the local glaciological setting, also the local topography and ice shelf geometry control lake distribution. Over George VI Ice Shelf, meltwater accumulated in flow stripes controlled by the prevailing compressive flow regime (Fig. 5a) (LaBarbera and MacAyeal, 2011; Hambrey et al., 2015). Further, lakes on all three ice shelves appeared at low elevations and low surface slopes (Fig. 5a–c), being favourable to meltwater ponding (Arthur et al., 2020b; Stokes et al., 2019). Additionally, lake ponding was dictated by the ice shelf geometry and grounding line location where lakes tended to cluster (Fig. 5a–c). Other local factors that could have influenced lake distribution are the refreezing of meltwater in firn pore space and channelized basal melting, causing a modified ice surface topography (Spergel et al., 2021; Dow et al., 2018; McGrath et al., 2012).
In addition, cross-correlation analysis demonstrated a positive relationship
between API supraglacial lake formation and air temperature, also considering temperature extremes, with lag 1 prevailing in regions where lakes form during high-melt years and lag 0 in regions where lakes form during low-melt years closer to the grounding line and rock outcrop (Figs. 5a–c, 6, and 7b, j, r). Moreover, the observed correlation pattern reveals that anomalously high or low lake extents result from anomalously high or low air temperatures (Fig. 6b). On the one hand, this indicates that meltwater ponding is initiated directly after temperature rises, particularly in regions where meltwater accumulation is promoted through melt–albedo feedbacks near rock outcrop, as present upstream of the grounding line of Wilkins and George VI ice shelves (Fig. 5a and c), while lake formation farthest from these regions occurs at slower timescales. On the other hand, meltwater transport across the ice shelves (Kingslake et al., 2017) later in the melting season can likely be attributed to higher lags downstream of the grounding lines (see connected drainage systems in Fig. 2a–c). The detected link between supraglacial lake formation and air temperature was previously
reported for other Antarctic regions (Kingslake et al., 2015; Langley et al., 2016; Dell et al., 2020) and is well in line with the analysis in Banwell et al. (2021), suggesting that high meltwater ponding in 2019–2020 over George VI Ice
Shelf results from sustained air temperatures
Over George VI and Bach ice shelves, solar radiation exhibited a positive relationship with supraglacial lake occurrence, particularly in regions where lakes cluster during high-melt years (Figs. 5a, b and 7c, d, k, l). This points towards the amount of incoming solar radiation being an equally important driver for supraglacial lake formation on API ice shelves. Similar conclusions were drawn by Laffin et al. (2021), who proposed that enhanced solar radiation is a main driver for foehn-induced melt on Bach and George VI ice shelves. In contrast, correlation with anomalies of solar radiation over Wilkins Ice Shelf returned a negative link at lag 1, which could indicate that lake formation results from enhanced longwave radiation, e.g. due to cloud cover and thus the higher importance of effects such as sensible heat transport, as discussed in Laffin et al. (2021) and Banwell et al. (2021). As correlations over Wilkins Ice Shelf were rather heterogeneous overall, the detected pattern could also reflect spurious correlation. However, weak negative lag 0–1 correlation was also found for George VI and Bach ice shelves (Table S2), reflecting an anomalously low availability of incoming solar radiation during the last two melting seasons (Fig. S6c and d).
Furthermore, cross-correlation of supraglacial lake extents over Wilkins Ice Shelf with anomalies of wind direction (Fig. 6b) point towards the existence of a negative link over regions of frequent lake coverage (Figs. 5c and 7u). Analysing the corresponding time series in more detail, negative correlations seem to be caused by predominant occurrence of north(-easterly) winds (i.e. north-easterly winds that are inclined more towards northerly direction) preceding recent increased lake ponding (Fig. S6). At the same time, positive correlations were visible, reflecting short-lived shifts to north-westerly winds preceding periods of high lake formation. Similar shifts in wind direction are also visible along the western and eastern grounding lines of George VI and Bach ice shelves (Fig. 7e and m). Wind conditions are a known driver for supraglacial lake occurrence on the API, with previous studies reporting that warm air temperatures over George VI Ice Shelf are promoted by sensible heat transport with warm north-westerly and north-easterly winds as well as a shift towards more frequent north-easterlies (Banwell et al., 2021). Moreover, foehn winds are known to influence surface melting on the north-eastern API (Cape et al., 2015; Datta et al., 2019; Luckman et al., 2014; Turton et al., 2020) and can equally be observed on the western API with (north-)easterly prevailing wind direction (Laffin et al., 2021). This suggests the occurrence of both foehn conditions and sensible heat transport on the API. Foehn winds could also explain the clustering of lakes near the grounding lines, where foehn wind influence is assumed to be highest (Lenaerts et al., 2017), as well as the positive link with lag 0–1 snowmelt (Fig. 7) (Cape et al., 2015; Turton et al., 2020). Further, the detected relationship is in good agreement with our observations of air temperature and solar radiation, where radiation potentially influences foehn conditions, and foehn conditions influence air temperature.
For investigation of large-scale climatic drivers, annual SAM was correlated with January lake extents. Our results imply a negative relationship for all investigated API ice shelves (Table 1), mostly due to occurrence of a negative SAM year in 2019 (Fig. 10) coincident with enhanced melting in 2019–2020 (Fig. 4a–c). Although the detected relationship seems counterintuitive because negative SAM years usually lead to cold conditions on the API, Bevan et al. (2020) suggest that teleconnections between the tropics and high latitudes in 2019 are responsible for the observed 2019–2020 melt and negative SAM. In detail, these teleconnections caused unusual perturbations of Southern Hemisphere atmospheric flow as well as sudden stratospheric warming (SSW), thus causing SAM to be negative and the ongoing positive Indian Ocean Dipole (IOD) in 2019 (BOM, 2021), represented via the dipole mode index (DMI) (Fig. 10), to bring warm air to the API via northerly winds. This overarching driver is well in line with correlations between lake extents and air temperature as well as wind and could potentially explain increased 2019–2020 meltwater ponding. As DMI returned to neutral values in 2020, other overarching drivers, including SAM returning to positive values (Fig. 10) as well as modified pressure conditions in the Amundsen Sea Low (ASL) (NSIDC, 2021), are believed to have contributed to high and low lake extents during the early and late 2020–2021 melting season, respectively.
Annual observations of Southern Annular Mode (SAM) and dipole mode
index (DMI) between 1981 and 2020. DMI captures the east-to-west sea surface
temperature gradient across the tropical Indian Ocean. Years 2015–2020 are
shaded in turquoise, representing our study period. SAM data are from Marshall (2018), and DMI data are available at
Overall, our results indicate that the complex interplay of local, regional and large-scale environmental controls dictates supraglacial lake formation on the API. In particular, the local glaciological and microclimatic setting including melt–albedo feedbacks, sensible heat exchange or the ice shelf topography was found to play an important role for supraglacial lake distribution. Likewise, long-term precipitation conditions influence the FAC and were found to be a primary control on the location and strength of meltwater ponding over all ice shelves. On the other hand, regional climatic conditions and particularly air temperature, solar radiation and (foehn) wind conditions are among the main controlling factors for temporal dynamics of supraglacial lakes in 2015–2021 and are mainly driven by large-scale atmospheric modes including SAM, IOD or ASL. In consequence of the good agreement between lake evolution and detected driving factors on all ice shelves (Table 2), a similar response of other API ice shelves resulting in high meltwater coverage during 2019–2020 and 2020–2021 is suggested.
Summary of suggested controls on 2015–2021 meltwater ponding across
six Antarctic ice shelves. The results present a summary of pixel-based and
spatially averaged correlations, also considering weak correlations and
neglecting spurious results where appropriate. Positive and negative signs
indicate positive and negative relationships, respectively. The signs “
For EAIS ice shelves, the observed negative relationship with short-term
precipitation over Riiser-Larsen Ice Shelf again highlights that supraglacial lake formation is controlled by the local glaciological setting and a reduced FAC (Figs. 6 and 8g, h). At the same time, weak negative correlation (
Furthermore, cross-correlation analysis with air temperature indicates that supraglacial lake formation is driven by temperature at lags 0–2 (Figs. 6 and 8a, b, i, j, q, r). Similarly to what is found for API ice shelves, the spatial pattern of temperature time lags (Fig. 8b, j and r) revealed larger lags for regions where lakes are spreading onto the ice shelves during years of high meltwater ponding and lower lags close to and upstream of the grounding lines. This likely reflects the temporal evolution of drainage systems through lateral meltwater transport across the ice shelves (Dell et al., 2020; Spergel et al., 2021; Kingslake et al., 2017) (Figs. 2d–f, 5d–f and S3) or a longer response time to air temperature, particularly in regions where the local glaciological setting does not promote supraglacial lake formation through melt–albedo feedbacks.
Cross-correlation with solar radiation returned positive relationships at lags 0–2 over all three ice shelves (Figs. 6 and 8c, d, k, l, s, t). As above, time lags were highest where lakes expanded farthest onto the ice shelves during years of high lake coverage, likely reflecting similar effects as described for air temperature also given the widespread distribution of blue ice or rock near all three EAIS ice shelves (Figs. 1d–f, 2d, and 5e, f). Similarly, the availability of incoming solar radiation was shown to contribute to enhanced melting over Roi Baudouin Ice Shelf, EAIS, where katabatic-wind-driven blue ice and firn exposure initiated albedo lowering and enhanced absorption of solar radiation (Lenaerts et al., 2017).
Moreover, we note that supraglacial lake formation on Riiser-Larsen Ice Shelf is dictated by prevailing low wind speeds from a (south-)easterly direction following periods of anomalously high wind speeds from a north-easterly direction (Figs. 6, 8e, S5a, and S7a, b). Considering the detected pattern, we suggest that meltwater ponding on Riiser-Larsen is driven by katabatic wind forcing. A similar pattern can be observed over Nivlisen Ice Shelf (Fig. 8m). Here, shifts in wind direction from low-speed south-easterly winds to strong (south-)easterly katabatic winds preceding high supraglacial lake coverage can be attributed to the observed correlations (Fig. S7c and d, Table S2) and meltwater ponding, e.g. by preserving the nearby blue ice region near the Wohlthat Mountains (Fig. 1) (Horwath et al., 2006). The spatial correlation pattern over Amery Ice Shelf is twofold (Figs. 8u, S5i and S7e, f). Along the north-east grounding line, the correlation pattern revealed that supraglacial lakes peak shortly after occurrence of high-speed (south-)easterly winds, likely reflecting the influence of katabatic winds flowing onto the ice shelf via the adjacent eastern plateau (Stokes et al., 2019). In contrast, supraglacial lakes in the south of Amery Ice Shelf peak during low-speed winds from the south(-east) or south-(west) following high-speed southerly katabatic winds. This indicates that katabatic winds drive melting in agreement with the local topography, determining the direction of wind inflow. Easterly katabatic winds warm adiabatically through vertical mixing and can lead to the exposure of albedo-lowering blue ice and firn through wind scouring, which further enhances surface melt, as observed over other EAIS ice shelves (Arthur et al., 2020a; Lenaerts et al., 2017). In fact, lag 1 correlation with modelled snowmelt events over Riiser-Larsen and Amery ice shelves (Figs. 6 and S5c, d, k, l) could reflect localized katabatic wind-induced peaks in melting.
In addition, the negative correlation with SAM over Riiser-Larsen Ice Shelf (Table 1) could point towards Southern Hemisphere atmospheric modes similarly affecting western East Antarctica, e.g. due to the spatial proximity to the Weddell Sea. In fact, January 2020 lake extents were above average with respect to the overall period 2016–2021, in agreement with a positive DMI and negative SAM (Fig. 10). Correlations with SAM over Nivlisen and Amery ice shelves were weakly positive, potentially reflecting the opposite climate effects SAM can have on the API and EAIS (Kwok and Comiso, 2002; Wachter et al., 2020). As mentioned, also the influence of an atypical ASL driving below-average temperature in 2020–2021 (NSIDC, 2021) has to be considered, most likely explaining the observed low lake extents over EAIS ice shelves in 2020–2021.
To summarize, EAIS supraglacial lake formation in 2016–2021 was mainly controlled by the regional near-surface climate, with air temperature, solar radiation and (katabatic) wind conditions showing the strongest influence. In this context, some of the observed climate variability can be explained by large-scale climatic drivers including SAM, IOD or ASL. Further, local controls including melt–albedo feedbacks, the firn air content and the ice surface topography showed strong influence on EAIS supraglacial lake formation and distribution.
Table 2 provides a summary of all detected drivers for the API and EAIS. Even though each ice shelf showed slightly different driving factors for supraglacial lake formation, local controls as well as regional climate variability with respect to air temperature, solar radiation and wind conditions were found to influence all six investigated regions in a similar manner. This points towards their transferability to other ice sheet regions as well as their role as primary driving factors. Furthermore, the potential influence of large-scale atmospheric modes at continent-wide scale is highlighted through the influence of the intense ASL in 2020–2021 on air temperature and lake formation across all six Antarctic ice shelves. Considering this, similar trends in lake evolution can be expected, e.g. in terms of low lake extents over other Antarctic ice shelves during the 2020–2021 summer melting season. Likewise, also the influence of SAM and IOD was shown to influence adjacent ice sheet sections (Table 2) over the API and EAIS, making it likely that similar lake evolution can be observed in surrounding regions, e.g. in the form of high lake extents during the 2019–2020 melting season over the API. At the same time, it has to be noted that the complex interplay of climate variables might have led to non-linear relations not captured in our correlations, and the comparatively low number of observations might have led to correlations being flagged as insignificant despite the existence of causal links. Further, also the coarse resolution of climate reanalysis data might have affected correlation values particularly at a small scale. Finally, more detailed investigation of teleconnections would require additional atmospheric indices.
Future atmospheric warming is expected to enhance meltwater production and ponding on Antarctic ice shelves (Bell et al., 2018; Kingslake et al., 2017; Trusel et al., 2015; Gilbert and Kittel, 2021). Extensive surface ponding threatens ice shelf stability due to variations in stress, causing flexure (Banwell et al., 2019) and ultimately vertical fracturing and ice shelf collapse, effectively reducing ice shelf buttressing (Banwell and Macayeal, 2015; Dunmire et al., 2020; Leeson et al., 2020; Rott et al., 2018; Scambos et al., 2004). Ice shelves that are vulnerable to hydrofracturing are widespread around Antarctica (Alley et al., 2018; Lai et al., 2020), particularly when they are pre-weakened by fractures and crevasses (Lhermitte et al., 2020) or firn air depletion (Kuipers Munneke et al., 2014; Lenaerts et al., 2017; Spergel et al., 2021). Lateral meltwater transport across ice shelves delivers meltwater to regions with only little or no surface melt, effectively reducing the FAC and making them more vulnerable to hydrofracture (Dell et al., 2020). This is particularly prevalent in the absence of an efficient drainage system transporting meltwater to the ocean (Bell et al., 2017). Given repeated meltwater ponding and drainage across the investigated ice shelves during the last 5–6 melting seasons (Fig. 5), future ice shelf destabilization could be imminent considering the classification of respective ice shelf areas as potentially vulnerable to hydrofracture (Lai et al., 2020). This particularly applies to ice shelves on the API, where exceptional meltwater ponding during the last two melting seasons likely contributed to further firn air depletion, as suggested for Larsen C (Bevan et al., 2020). Future enhanced precipitation due to an enhanced moisture holding capacity of warmer air may partly compensate this (Bell et al., 2018; Bengtsson et al., 2011; Krinner et al., 2007; Lenaerts et al., 2016) even though additional research is required for a better evaluation of all factors, e.g. considering an increasing presence of firn aquifers under high-precipitation scenarios with similar effects on ice shelf stability (Bell et al., 2018; van Wessem et al., 2021). Similarly, increasingly positive IODs under high-emission scenarios may counteract API warming in the absence of perturbations such as SSW (Bevan et al., 2020).
This study presents the first intra-annual and inter-annual investigation of
Antarctic surface hydrology across six major Antarctic ice shelves, employing
state-of-the-art machine learning on spaceborne observations from Sentinel-1
SAR and optical Sentinel-2. We derive fused classification products at
unprecedented 10 m spatial resolution and bi-weekly temporal scale to map
Antarctic supraglacial lake extent in 2015–2021. In this context, the
integration of SAR data enabled the delineation of buried and partly frozen
lakes, providing a more complete mapping record than single-sensor products.
For the investigated Antarctic Peninsula ice shelves, our results reveal
anomalously low supraglacial lake coverage during most 2015–2018 melting seasons and anomalously high lake coverage during 2019–2020 and 2020–2021 summers considering the years 2015–2021 as a reference period. Peak lake
coverage reached 805 km
To establish a link between intra-annual supraglacial lake extent dynamics in 2015–2021 and potential environmental drivers, we employ multi-temporal linear correlation analysis at different time lags using spatio-temporal observations of ERA5-Land climate reanalysis data. Moreover, inter-annual observations of supraglacial lake extents are used to investigate the influence of large-scale atmospheric modes as well as of the local glaciological setting, including the state of the firn layer. Results of the statistical analysis with spatial averages and anomalies as well as with pixel-based values return a strong coupling between temporal fluctuations in Antarctic Peninsula supraglacial lake extents and the regional near-surface climate, including air temperature, the amount of incoming solar radiation and (foehn) wind conditions, with lake extents above or below average being linked to anomalies of climate variables. In addition, correlations with large-scale atmospheric drivers including annual Southern Annular Mode revealed their influence on the regional near-surface climate and thus supraglacial lake formation. Further, investigation of previous-year precipitation implied the strong influence of reduced-snow conditions, promoting supraglacial lake formation through a reduction in the firn air content. Also, pixel-based analysis revealed a faster response time of lake evolution to climate drivers in regions where surface features such as blue ice or rock initiate melt–albedo feedbacks. Similarly to Antarctic Peninsula ice shelves, air temperature, solar radiation and (katabatic) wind conditions were found to be primary drivers for supraglacial meltwater ponding on East Antarctic ice shelves. Likewise, the influence of Southern Hemisphere atmospheric modes on the regional near-surface climate as well as the local glaciological setting, including the ice shelf geometry, elevation and melt–albedo feedbacks, is found to exert strong control on supraglacial lake formation. Moreover, the influence of a reduced firn air content over Riiser-Larsen and particularly Amery Ice Shelf was related to above-normal supraglacial lake formation.
To summarize, our findings provide important insight into present-day Antarctic surface hydrology and associated drivers and could be particularly useful for an improved representation of meltwater transport across Antarctic ice shelves in modelling efforts (Buzzard et al., 2018). Specifically, the created dataset on bi-weekly supraglacial lake extent dynamics over six Antarctic ice shelves is particularly useful for future analyses on supraglacial lake characteristics or temporal variations over individual lakes and drainage structures. In addition, the results of this study underline the complex interplay of climate variables at different time lags and the local glaciological setting, which influences supraglacial lake formation both spatially and temporally, with similar implications for other ice sheet regions. Furthermore, the value of ERA5-Land and Copernicus Sentinel data was highlighted, with the growing archive of satellite data enabling more detailed analyses on supraglacial lake dynamics in the future. In particular, future work will focus on the generation of circum-Antarctic supraglacial lake extent mapping products, to be made available via the GeoService of the Earth Observation Center (EOC) at the German Aerospace Center; the expansion of statistical analyses to continental scale; and the application of advanced causal analysis techniques. Once available, the integration of higher-resolution model outputs from RACMO (Regional Atmospheric Climate Model) (van Wessem et al., 2018) as well as more detailed investigation of teleconnections with additional atmospheric indices will be considered to better constrain climate drivers.
The Supplement provides results from the statistical correlation of fractional supraglacial lake extents and bi-weekly spatial averages (Table S1) and anomalies (Table S2) over all six investigated Antarctic ice shelves. Moreover, Fig. S1 provides information on the share of data originating from Sentinel-1 and Sentinel-2 in bi-weekly classifications. Figure S2 shows larger map plots of Fig. 5, and Fig. S3 shows the temporal development of January supraglacial lake extents in 2017–2021 over Nivlisen Ice Shelf. Figures S4 and S5 provide pixel-based correlation results for all variables not shown in Figs. 7 and 8. Figures S6 and S7 show auxiliary plots on spatially averaged wind direction and lake extents for Wilkins, Riiser-Larsen, Nivlisen and Amery ice shelves to complement the conclusions made in Sect. 5. Figure S6 also shows results for lag 1 solar radiation.
The bi-weekly maximum lake extent mapping products for all six ice shelves
are available via the EOC GeoService of the German Aerospace Center (DLR)
(
The supplement related to this article is available online at:
MD designed the study, conducted the analysis and wrote the original draft of the manuscript. CK and AD assisted in the study design. MD, CK and AD contributed to the discussion of the results and were involved in editing the manuscript.
The contact author has declared that neither they nor their co-authors have any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We thank the Copernicus Climate Change Service (C3S) Climate Data Store for providing ERA5-Land reanalysis data. Moreover, we thank the European Union Copernicus programme for providing Sentinel-1 and Sentinel-2 data. Matthias Hofmann and Andreas Jahnke, DLR, are acknowledged for their support with Sentinel data access, and Julian Zeidler and Soner Uereyen, DLR, are acknowledged for IT support and proofreading, respectively. The authors also acknowledge Gareth Marshall for granting open access to the Southern Annular Mode index and NOAA/ESRL for dipole mode index data. Finally, we thank the editor (Johannes Fürst) and two anonymous reviewers for their constructive and helpful remarks on the manuscript.
The article processing charges for this open-access publication were covered by the German Aerospace Center (DLR).
This paper was edited by Johannes J. Fürst and reviewed by two anonymous referees.