Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 13, Article number: 20619 (2023)
6
Metrics details
The liquid water around the Antarctic Ice Sheet plays a key role in modulating both the vulnerability of ice shelves to hydrofracturing and ice discharge from outlet glaciers. Therefore, it needs to be adequately constrained for precise future projections of ice-mass loss and global sea-level rise. Although glacial lake outburst floods (GLOFs) pose one of the greatest risks in glacierized mountainous regions, any long-term monitoring of Antarctic ice-marginal lakes and their associated potential for GLOFs has been neglected until recently owing to the limited number of such events reported in Antarctica. Here we present direct evidence of repeated GLOFs from Lake Kaminotani-Ike, an ice-sheet-dammed lake in East Antarctica, via an analysis of historical aerial photographs and recent satellite data. Two GLOFs occurred in 1969–1971 and 2017, with discharge volumes of (8.6 ± 1.5) × 107 and (7.1 ± 0.4) × 107 m3, respectively, making them two of the largest GLOFs in Antarctica. A southerly oceanward pathway beneath the ice sheet is the most likely drainage route of these GLOF events based on the available surface- and bed-elevation datasets. Furthermore, the 2017 event occurred during the austral winter, thereby implying the possibility of year-round active subglacial networks in Antarctica. Our results highlight that studies on Antarctic ice-marginal lakes provide an opportunity to better understand Antarctic hydrological processes and emphasize the need for both detailed monitoring of ice-marginal lakes and detailed surveying of the subglacial environments of the Antarctic Ice Sheet.
The Antarctic Ice Sheet (AIS) is the largest potential contributor to sea-level rise1, 2. Recent studies have revealed that the hydrology of Antarctica plays a key role in sea-level rise by modulating the ice dynamics of AIS outlet glaciers and ice shelves3,4,5,6,7. The ongoing increase in supraglacial lakes has the potential to destabilize ice shelves via hydrofracturing3, 7, 8, although supraglacial streams can weaken this process by reducing lake storage via water transport9. Subsequent acceleration of the inland ice may occur if the ice shelf collapses10,11,12 and/or meltwater drainage reaches the subglacial environment3,4,5. Outburst floods from subglacial lakes have accelerated the downglacier flow of outlet glaciers by altering the subglacial hydrological environment5, 13. However, the hydrological processes behind such events are not well constrained because of the remoteness of the subglacial environment. Despite the importance of cryo-hydrology in the evolution of AIS, studies of ice-marginal and/or ice-dammed lakes are extremely limited, thereby limiting our understanding of AIS hydrology.
Glacial lake outburst floods (GLOFs), which abruptly discharge their stored water and lake materials downstream, are a common phenomenon that has been observed at ice-marginal lakes in glacierized regions worldwide, with the exception of Antarctica, and are recognized as one of the greatest risks in mountain societies14,15,16,17. Furthermore, these freshwater pulses can impact ocean circulation, biochemical circulation, and marine ecosystems18,19,20. For example, GLOFs from the largest-scale glacial lakes that were located along the periphery of past ice sheets weakened ocean circulation, thereby modulating climate21, 22. However, GLOFs from ice-marginal and ice-dammed lakes in Antarctica have largely been neglected due to extremely limited reports of GLOFs until recently23,24,25,26. Although outbursts from subglacial, supraglacial, and sub-surface lakes are generally associated with either a collapse of the surrounding ice or ice doline formation27,28,29, the subglacial and englacial hydrological processes behind such events are not well constrained due to the rarity of these events and the inherent difficulty in obtaining direct observations. The lack of ice-marginal lake monitoring in Antarctica limits our understanding of the impacts that GLOFs have on ice-sheet/glacier dynamics and the surrounding environment. Furthermore, constraints on GLOFs in Antarctica may advance our understanding of lake–ice interactions in harsh environments, even though they pose a lower risk compared to other regions.
Lake Kaminotani-Ike (LKI) is an ice-dammed (and perennially ice-covered ice-marginal) lake in East Antarctica. LKI is situated ~ 60 km south of Syowa Station and is located between the Skarvsnes outcrop and the East AIS (EAIS; Fig. 1a). Liquid precipitation rarely occurs across the Lützow-Holm Bay region. For example, a rain event on 24 December 2017 at Syowa Station marked the first recorded since 201230. Therefore, the lake water likely originates from subglacial discharge that is directly fed by EAIS and meltwater of solid precipitation over the catchment. Members of the 59th Japanese Antarctic Research Expedition (JARE59) observed several unusual features at LKI on 22 November 2017, including the exposed lake floor, a vertical ice cliff at the boundary between the lake and EAIS, and snow and ice blocks scattered along the lake shore31 (Fig. 1b). All of these features implied rapid surface lowering. Although information about LKI is limited, previous studies have suggested the possibility of dynamic changes at LKI32, 33 and captured details on the drainage and periodicity of a nearby ice/snow-dammed lake23, 24.
Overview of the LKI GLOFs. (a) Overview of the study region. The background image is a Landsat 8 satellite image that was taken on 22 February 2017. The inset map shows the location of the study region relative to Antarctica. (b) Field photographs taken on 10 January 2018, after the 2017 LKI GLOF. (c, d) Landsat 8 satellite images taken on (c) 22 February and (d) 11 October 2017. The contours were taken from the REMA-DEMs on (c) 29 September 2016 and (d) 1 October 2017. Black arrows highlight the location of the exposed lake floor shown in (b). (e, f) PALSAR-2 satellite backscatter intensity images taken on (e) 13 April and (f) 15 May 2017. (g, h) Orthoimages generated from aerial photographs taken on (g) 21 January 1969 and (h) 24 February 1971. Surface-elevation changes during the (i) 29 September 2016 to 1 October 2017 and (j) 21 January 1969 to 24 February 1971 periods.
Aerial photographs had been acquired across Lützow-Holm Bay, East Antarctica, by JARE since 1957 (Supplementary Fig. 1). The photographs had been taken across the Skarvsnes region during the austral summer over the 1957–2000 period. Development of SfM techniques enable us to trace surface elevation by constructing digital elevation models (DEMs) from historical aerial photographs. This aerial photograph collection, coupled with recent satellite data (imagery and DEMs), provides a ~ 60-year time series of the surface state and surface elevation of LKI. In this paper, we present direct evidence of repeated GLOFs from an Antarctic ice-dammed lake and infer their potential periodicity via an analysis of long-term historical aerial photographs and recent satellite data.
We found two GLOF events in 1969–1971 and 2017 from satellite imagery and aerial photographs (we hereafter refer to the GLOF periods as P1 and P2, respectively). A comparison of Landsat 8 satellite images acquired on 22 February and 11 October 2017 clearly highlighted changes in the surface features of LKI (Fig. 1c and d). The lake surface was flat due to frozen lake ice on 22 February (Fig. 1c), whereas on 11 October the surface was rugged and the lake-floor terrain was exposed (Fig. 1d), with the latter remote-sensing-based features supported by field observations (Fig. 1b). Furthermore, PALSAR-2 synthetic aperture radar (SAR) images acquired on 13 April and 15 May 2017 highlighted clear changes in both the lake-surface features and backscatter intensities between the two images (Fig. 1e and f). Although the lake surface was flat and possessed a uniform spatial distribution of backscatter intensities on 13 April 2017 (Fig. 1e), rugged features and an increase in backscatter intensity were observed at the periphery of the lake on 15 May 2017 (Fig. 1f). The 15 May PALSAR-2 image detected a lake shoreline with a similar shape to that in the 11 October 2017 optical image (Fig. 1d), which indicates that the drainage of P2 occurred sometime during the April–May 2017 period and was at most a one-month-long drainage event. We detected another LKI drainage event during the 1969–1971 period that exhibited similar surface feature changes to those associated with the 2017 GLOF event, whereby a flat surface was observed until 1969 and then became rugged, with the lake-floor terrain being exposed in 1971 (Fig. 1g and h). No other abrupt changes in the surface features of LKI have been observed in the aerial photographs and satellite images (Supplementary Fig. 2).
We traced the LKI surface-elevation changes using satellite-based DEMs (REMA34, 35-, ALOS-, GSI-, and TDX-DEMs; see the Methods section) and DEMs that were constructed by applying the structure-from-motion (SfM) technique to the JARE aerial photographs (SfM-DEMs). DEMs clearly capture a drop in surface elevation across LKI in 1969–1971 and 2017 (Fig. 1c,d, and 2c and Supplementary Fig. 2). The surface elevation dropped by 65.8 ± 9.1 and 55.0 ± 2.1 m, which indicates discharge of (8.6 ± 1.5) × 107 and (7.1 ± 0.4) × 107 m3 from LKI during P1 and P2, respectively (Figs. 1i, j, 2, and 3). A maximum surface lowering of 80.5 ± 2.1 m was observed near EAIS front during P2 (Fig. 1i). This deepest part was not captured in the 1971 DEM, which suggests that either the exact drainage of P1 occurred long enough before the acquisition date (24 February 1971) for refilling to have already begun, or LKI may not have fully drained during P1 (Fig. 2). The total discharge volumes from LKI during P1 and P2 are four times larger than the 2019 GLOF event at Lake Untersee, a glacier-dammed lake in East Antarctica25, thereby making these LKI events two of the largest GLOFs reported at an ice-marginal lake in Antarctica to date.
LKI elevation profiles. Surface elevation along profile AA′ (Fig. 1c), which crosses LKI, during the (a) 1962–1971, (b) 1992–2017, and (c) 2018–2020 periods. (d) Close-up view of (c) to concentrate on the lake-surface uplift after the 2017 GLOF. Note that the color schemes in (b) and (c) are different.
Lake-surface elevation and volume time series. Time series of (a) lake-surface elevation and (b) relative lake volume, which is referenced to the mean of the REMA period. The “DEM” data points are derived from the DEMs (SfM-, REMA-, GSI-, TDX-, and ALOS-DEMs). The “REMAempty + Image” data points are derived from the lake shoreline method, which employed REMAempty in (a). The “ICESat-2” data points are from the ICESat-2 ATL06 product (Supplementary Fig. 7).
We analyzed all of the available elevation datasets (aerial photography, satellite laser altimetry, and a combination of optical satellite imagery and a post-GLOF DEM) to construct a dense time series of the lake-surface elevation and relative lake volume (Fig. 3). The lake volume increased by 1.7 × 106 m3 a−1 between P1 and P2 (the 1971–2017 period) and 2.5 × 106 m3 a−1 after P2 (the 2017–2021 period). The time series clearly highlights two GLOFs that are followed by periods of gradual lake filling (Fig. 3a). The changes in lake-surface elevation and lake volume exhibit a typical saw-tooth pattern (Fig. 3) that is similar to the cyclic patterns of rapid drainage and slow filling that have been observed at ice-dammed lakes in other regions36,37,38.
Outbursts from ice-dammed lakes are generally caused by either a break in the seal (between the glacier bottom and bedrock) or overflow of the ice dam/lake margin. Seal-break mechanisms are classified as (1) flotation36, 39, (2) conduit formation beneath/within the ice dam40, and (3) mechanical collapse41. Evidence of these different outburst mechanisms have been observed in the field at Antarctic glacial lakes, even though Antarctic GLOF studies are limited. For example, the 2019 Untersee GLOF was initiated by water outflow at the lake margin, which was confirmed by field inspections25. The Yatsude Valley GLOF was probably caused via the subglacial hydrological system, with the formation of a near-circular hole at/near the base of the ice dam inferred23. There were no observations of a discernible hole, trace of failure in the ice dam, or surface-water overflow during the 2017 field inspection and in the images that were acquired after the 2017 LKI GLOF event. However, the lake-surface elevations immediately prior to the 1969–1971 and 2017 LKI GLOFs were quite similar. Therefore, we assumed that the triggering mechanism of the LKI GLOFs was related to the flotation criterion of the ice dam, as observed for GLOFs from other ice-dammed lakes36, 42. Even if the ice dam achieved flotation, the glacier uplift may not be detectable in the satellite-based DEMs due to the expected small-scale vertical motion of the ice dam (several centimeters)39.
This flotation hypothesis could be investigated by analyzing the ice-surface, ice-sheet bed, and lake-surface elevations. The elevation difference between the lake level and flotation level, (Delta z), is denoted as follows:
where (rho_{i}) and (rho_{w}) are the ice and water densities, respectively; and (h_{l}), (h_{s}), (h_{b}), and (h_{f}) are the lake-surface, ice-surface, ice-sheet bed, and flotation threshold elevations, respectively. We employed (h_{b}) from the BedMachine dataset43 because it is currently the only available bed-elevation dataset that covers our study region. The regions where (Delta z) ≤ 0 indicate that the potential for flotation exists when the subglacial environment is hydrologically connected to the lake. The spatial distribution of (Delta z) highlights that the lake level is below the flotation threshold around LKI, thereby ensuring that the lake is dammed by the ice sheet (Fig. 4a). The area to the south of the LKI–EAIS boundary is the closest to the flotation criterion, whereas the other areas along the ice-dammed boundary rapidly increase to large positive (Delta z) values, such that these areas would prohibit flotation (Fig. 4b).
Estimation of discharge routing. (a) Spatial (Delta z) distribution, generated using the REMA-DEM on 29 September 2016 and the BedMachine version 3 dataset43. Yellow lines indicate where (Delta z) = 0. The background image is an orthoimage generated from aerial photographs taken on 19 January 2000. The box indicates the extent of Fig. 5. (b) Cross section of the surface and bed elevations along profile ABC, with the position of the profile shown in (a). Black and red lines indicate the REMA-DEM surface elevations before (29 September 2016) and after (1 October 2017) the 2017 GLOF event, respectively. The gray line with dots indicates the subglacial elevation from BedMachine43. The broken line indicates the floating condition of ice-dam, which satisfies the flotation criteria based on the ice-sheet and lake-surface elevations on 29 September 2016.
The southerly route (approximately parallel to profile BC in Fig. 4b) is the most likely drainage route. However, there is a large discrepancy between the bed elevation from BedMachine43 and floating condition of ice-dam that calculated from the surface elevations of the ice sheet and lake, which implies that the available bed-elevation data do not satisfy the flotation threshold. A > 10-m increase in the lake level is required to reach flotation based on the available bed-elevation constraints. Nevertheless, the actual threshold is less than ~ 20 m due to the effect of the gradual opening of the subglacial seal44. This assumption is based on the idea of ice deformation on the lakeward side as a “cantilever” effect, whereby the accumulated lake water induces buoyancy-driven uplift of the glacier-ice. Furthermore, the errors in the bed-elevation dataset introduce uncertainties in the exact location of the drainage route owing to the coarse resolution (500 m) and large uncertainties (100–1000 m) of the employed bed-elevation dataset. We note that these large uncertainties are due to a lack of measurements across the study region and the fact that the BedMachine dataset43 is designed for large-scale studies that incorporate either the entirety of Antarctica or regional catchments. We observed > 100 m of ice-cliff retreat (Fig. 5a and b) and subsidence of the ice-sheet surface and collapsed ice near the lake (Fig. 5c) after the 1969–1971 event. These features are similar to the ice-sheet surface features that have been observed after subglacial drainage at the Antarctic Peninsula27. The observed retreat following the 2017 event was smaller than that following the 1969–1971 event, and there was no clear subsidence/collapse feature in the ice-sheet surface (Fig. 5d–f). Nevertheless, the surface subsidence that was observed following the 1969–1971 event occurred in a region that did not satisfy the flotation criterion (Fig. 4), which implies uncertainties in the bed-elevation dataset43 for this region. Furthermore, other ice-dammed lakes located on the northeastern side of LKI showed neither changes in area nor water overflow. This implies that the northern drainage route was not activated, which supports a southerly drainage route. Therefore, we conclude that the LKI GLOFs occurred owing to flotation of the ice dam, with the southerly route serving as the most likely drainage pathway.
Comparison of the surface features at LKI and the surrounding EAIS. Changes in the surface features and elevation during the (a–c) 1969–1971 GLOF and (d–f) 2017 GLOF events. Orthoimages on (a) 21 January 1969 and (b) 24 February 1971. (c) Surface-elevation change between 21 January 1969 and 24 February 1971. Shaded relief from REMA on (d) 29 September 2016 and (e) 1 October 2017. (f) Surface elevation change between the 29 September 2016 and 1 October 2017 REMAs. Magenta lines denote the LKI margins prior to the GLOF events, and cyan lines denote the edge of the ice sheet after the GLOF events.
The 2017 PALSAR-2 images revealed that the 2017 LKI GLOF event occurred during the austral winter (Fig. 1e and f). Although flotation of the ice dam is the most likely mechanism of the LKI GLOFs, it is unclear what triggered the breach of the water-pressure threshold, which then allowed the formation of a drainage pathway beneath the ice dam. Ice/snow meltwater could contribute to the increase in lake volume during the melt season in other regions42, 45. However, the surface temperature was always below 0 °C across the study region, and no significant snowfall events were observed at Syowa Station during the April–May 2017 period (Supplementary Fig. 3). Additional water inputs from the subglacial environment, such as the drainage from either a nearby subglacial lake or a developed subglacial hydrological system, may assist in reaching the water-pressure threshold for LKI drainage. However, there are no known subglacial lakes near LKI, and the subglacial thermal regime near LKI is unknown.
The volume of LKI increased by 1.7 × 106 m3 a−1 during the 1971–2017 period and 2.5 × 106 m3 a−1 during the 2017–2021 period, despite such a limited supply of liquid water via atmospheric processes. Furthermore, there were no distinct signals from summer meltwater input in the surface-elevation and lake-volume time series (Fig. 3). These results imply a constant water input from the subglacial hydrological network on a year-round, such that the subglacial water pressure at the ice dam can reach the threshold via constant subglacial discharge during the austral winter, even if there are no other water inputs into LKI. Although the water supply is independent of the atmospheric conditions, it does depend on the geothermal flux and subglacial catchment geometry, which are basically invariable. The presence of an active subglacial network has been suggested by direct observations of submarine groundwater discharge in Lützow-Holm Bay46 and an equipotential hydrostatic model47. Nevertheless, the actual trigger of the LKI GLOFs is not constrained at this stage because the subglacial environment is generally frozen around the periphery of EAIS48.
Our long-term observations provide the first evidence of repeated GLOFs from an ice-marginal lake in Antarctica (Figs. 2 and 3). Although the lake-volume time series implies a similar pattern to those of cyclic GLOFs, the cyclicity of the LKI GLOFs is inconclusive due to the small number of the events. If we apply the simple assumption that the time between the two GLOFs defines the periodicity of LKI GLOFs, then LKI possesses a GLOF periodicity of ~ 50 years, which is longer than the GLOF periodicities at ice-marginal lakes in other regions36, 49, 50. Meltwater availability and/or rainfall are key factors that contribute to these differences in GLOF periodicity, with local subglacial hydrological processes and the physical properties (temperature, softness, and permeability) of the ice damming the lake, which are currently unconstrained for LKI and the surrounding area. Recent increases in rainfall and intensive melting along the periphery of EAIS, which have been attributed to changes in atmospheric conditions51, 52, may both accelerate the occurrence of the next LKI GLOF event and increase the frequency of future GLOF events at LKI and other ice-marginal lakes. Nevertheless, there are limited examples of GLOFs associated with ice-dammed lakes in Antarctica, and the time series of these events are too short to distinguish their periodicity23,24,25,26. Thus, continued monitoring of the ice-dammed lakes in Antarctica is needed to further discuss the potential periodicity of GLOFs and how a changing global climate may impact this periodicity.
The large-scale GLOFs that have been observed from ice-marginal lakes in Greenland53 also highlight the need to monitor ice-marginal lakes in Antarctica, as continued warming will alter AIS dynamics and atmospheric conditions in the coming decades and centuries54, 55. Detailed mapping of the ice-dammed and ice-marginal lakes around AIS is recommended to assess the contribution of such lakes as potential GLOF sources, as has been done for the Greenland Ice Sheet56. Furthermore, our results imply a possible subglacial water input into the lake, even during the austral winter. Subglacial water has been observed beneath the Greenland Ice Sheet and is recognized as a key factor in modulating ice movement57, 58. In fact, substantial accelerations in glacier flow have been observed after lake drainage in Greenland59, 60 and Alaska61. Although the subglacial hydrology of these glacierized regions remains unknown, these examples imply that Antarctic GLOFs could potentially trigger accelerations in ice flow. Therefore, more detailed measurements of the subglacial topography around the periphery of AIS are required to both elucidate potential GLOF triggers and map where ice-marginal lakes may form owing to the future retreat and/or thinning of AIS, as has been conducted in other glacierized regions62.
Our results also imply the sporadic influx of freshwater to the ocean. However, our understanding of the freshwater supply from AIS to the surrounding ocean is currently limited. Recent modeling studies have shown that there is a limited contribution to the total basal melt of the ice shelves, with enhanced melt in sub-ice-shelf channels63, 64 due to subglacial freshwater discharge. Although the amount of freshwater that is discharged via LKI GLOF events is much less than the discharge in these modeling studies63, 64, a sporadic freshwater pulse could promote a similar enhanced melt effect if the water discharges toward ice shelves. Nutrients and sediment transport are also expected owing to the freshwater transport, which may have an influence on nearby ocean ecosystems considering changes in the chemical properties of Antarctic ice-marginal lakes after GLOF events13. Although constraints on the subglacial discharge from AIS to the ocean are currently quite limited, such that the impact of this freshwater flux remains unknown, it is important to monitor such lake and freshwater transport to advance our understanding of this phenomenon along the periphery of AIS.
In this study, we revealed the occurrence of repeated GLOFs from LKI, an ice-dammed lake in East Antarctica, during the 1962–2021 period. Several types of data (SfM-/satellite-DEMs, satellite altimetry, and a combination of the shoreline position and lake bathymetry) were utilized to quantify the changes in lake-surface elevation and lake volume. We estimated (8.6 ± 1.5) × 107 and (7.1 ± 0.4) × 107 m3 of water discharge during the 1969–1971 and 2017 GLOFs, respectively, which are two of the largest known GLOFs in Antarctica. Surface- and bed-elevation constraints have indicated that the GLOFs likely occurred via flotation of the ice dam, with water discharge occurring through a southerly drainage pathway from LKI. Furthermore, the 2017 GLOF occurred during the austral winter (April–May), thereby suggesting a year-round active subglacial hydrological network because there were no water inputs into LKI from surface processes.
Our results, coupled with the lack of subglacial constraints in the region, emphasize the importance of elucidating the subglacial drainage system of AIS and highlight that continuous monitoring of ice-marginal lakes may provide important insights into AIS subglacial hydrology. Radar and/or seismic surveys are required to provide more accurate and detailed constraints on the subglacial topography and internal structure of ice dams. Therefore, our results will encourage future monitoring of Antarctic ice-marginal lakes and additional observations of AIS hydrological processes.
We analyzed aerial photographs and optical satellite images to delineate the surface features across LKI. JARE has acquired aerial photographs across Lützow-Holm Bay since 1957, with aerial photographs taken across the Skarvsnes region during the austral summers of 1957, 1962, 1969, 1971, 1983, 1992, 1993, and 2000. We obtained these freely available photographs, which were scanned at 400 dpi, from the Geospatial Information Authority of Japan (GSI; https://www.gsi.go.jp/antarctic/06.html) for our time-series analysis. The spatial resolutions of the aerial photographs are 0.83–2.16 m. We also used optical satellite images acquired by the Landsat 1 Multispectral Scanner (MSS), Landsat 4 MSS/Thematic Mapper (TM), Landsat 5 MSS/TM, Landsat 8 Operational Land Imager (OLI), Sentinel-2 Multi Spectral Instrument (MSI), and Terra Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) (DatalistOpticalSatellite.xlsx). The spatial resolutions of these satellite images are 60 m for MSS, 30 m for TM, 15 m for OLI and ASTER, and 10 m for MSI. False-color composite images were constructed using the red, green, and near-infrared band images that were included in each product. Note that we used pansharpened false-color Landsat 8 images, which were constructed using the Landsat 8 panchromatic images in the QGIS software.
We analyzed both SAR and optical satellite imagery to investigate the intra-annual variations in the 2017 lake-surface features. We used two SAR images that were acquired on 13 April and 15 May 2017 by the Phased Array type L-band Synthetic Aperture Radar 2 (PALSAR-2) mounted on the Advanced Land Observing Satellite 2 (ALOS-2) satellite for this purpose. Radiometric calibrations, multi-looking (two looks and one look in the azimuth and range directions, respectively), speckle reduction via the Refined Lee filter, and ellipsoidal corrections were applied to the Level 1.1 data (single-look complex) using the freely available Sentinel Applications Platform (SNAP) software developed by the European Space Agency. A terrain correction was then applied using the REMA-composite DEM in SNAP.
The Landsat MSS and PALSAR-2 images possessed enormous horizontal geo-location errors (~ 2.4 km for MSS and ~ 300 m for PALSAR-2). We corrected these errors using the Georeferencer Plugin in the QGIS software and manually located ground control points (GCPs). Four GCPs were located within each image to determine their near-edge positions.
We used five types of DEMs for the Skarvsnes region to retrieve the surface-elevation changes across LKI: SfM-DEMs, a GSI-DEM, an ALOS-DEM, REMA-DEMs, and a TDX-DEM.
We applied the SfM tool in Metashape Pro version 1.7.1 software (©Agisoft, 2021) to the JARE aerial photographs to construct the SfM-DEMs. We first applied manual masking to eliminate the camera frames that were included in the periphery of the aerial photographs. We then followed the approach outlined in a previous study33 as follows. We first performed initial image matching, alignment, and three-dimensional (3-D) model reconstruction without any GCPs. We then extracted 53 GCPs from the 1993 DEM and orthoimage of the Skarvsnes region33 and manually located these GCPs on the aerial photographs (Supplementary Fig. 4). We subsequently performed the above-mentioned 3-D model reconstruction using the geo-located aerial photographs to construct a geo-located 3-D model. Finally, we constructed geo-located DEMs and orthoimages from the 3-D models. The spatial resolutions of the DEMs and orthoimages were 2.9–4.2 and 1.4–2.1 m, respectively. Note that the 1983 aerial photographs were not used for the DEM construction because the photographs over the Skarvsnes region were taken along a single flight path. We also did not use the 1957 images because they were oblique images, although they have enabled us to assess the state of LKI in 1957 (Supplementary Fig. 1).
The GSI-DEM was constructed from a pair of 2.5-m resolution stereo images that were acquired on 18 January 2009 by the ALOS/Panchromatic Remote-sensing Instrument for Stereo Mapping (PRISM). GSI performed the DEM construction, which possessed a 10-m spatial resolution. The vertical accuracy of the GSI-DEM was evaluated to be 3.82 m based on a comparison of static geodetic survey constraints and the GSI-DEM elevations at a series of outcrops across the Skarvsnes region33. We also applied the photogrammetric tools in the Ames Stereo Pipeline software65 version 2.7.0 to the nadir and backward stereo-pair ALOS/PRISM images that were acquired on 19 November 2007 to construct a DEM.
We used the REMA-DEMs34, 35 and other DEMs to investigate the recent surface-elevation changes across LKI. The REMA-DEMs were constructed by applying the photogrammetric tools in the SETSM software66 to panchromatic stereo images acquired by the WorldView-1/-2/-3 and GeoEye-1 satellites. We used 22 strip REMA-DEMs (version 4.1) to retrieve the surface-elevation changes across LKI. Each strip REMA-DEM was constructed from a pair of satellite stereo images and, therefore, possessed a precise acquisition date. The spatial resolutions of the original distributed strip REMA-DEMs are 2 m. We used the composite REMA-DEM as the reference DEM in this study because the composite REMA-DEM was constructed by mosaicking the strip REMA-DEMs and possessed a reliable vertical coordinate34. The spatial resolution of the composite REMA-DEM was 10 m, and the vertical accuracy was reported to be < 1 m34.
A TDX-DEM was generated from X-band SAR data that were acquired simultaneously by the TerraSAR-X add-on for Digital Elevation Measurement (TanDEM-X) and TerraSAR-X on 8 June 2017 during the bistatic operation mode. We applied interferometric SAR processing to the SAR data pair and then converted the obtained interferometric phase of each pixel to a surface height with a precisely determined perpendicular baseline length (Bp) between the two satellites. We used GCPs from AW3D30 version 3 (https://www.eorc.jaxa.jp/ALOS/en/dataset/aw3d30/aw3d30_e.htm) to determine the precise orbits during the Bp estimation. The generated TDX-DEM was processed using the GAMMA software version 20220630.
We applied the following corrections to the DEMs to remove the complex errors that are often included in DEMs67. We followed the method proposed by Nuth and Kääb67, whereby we removed the horizontal shifts, elevation-dependent biases, and sensor-dependent biases (along-/across-satellite track biases) from the original DEMs. We only applied horizontal shift and elevation-dependent bias corrections to the SfM-DEMs because the along-/across-satellite tracks were not determined for the SfM-DEMs. We then subtracted the second-order surfaces from the SfM-DEMs to remove any doming/bowl errors that are typically included in SfM-DEMs68, 69. The surfaces were fitted via the least-squares method to the difference between each SfM-DEM and the reference DEM. Note that the TDX-DEM includes long-wavelength noise because the TDX-DEM was constructed via SAR interferometry of the TanDEM-X and TerraSAR-X satellites. We subtracted the noise by applying the least-square method to estimate the second-order surfaces after the elevation-dependent bias-removal correction. We corrected the ALOS-DEM by removing the horizontal shifts and elevation-dependent biases, and the GSI-DEM, SfM93-DEM, and strip REMA-DEMs by removing the vertical offsets. Finally, all of the DEMs were aligned to the grid of the reference DEM to yield DEMs with a consistent 10-m spatial resolution (Supplementary Figs. 4 and 5).
The above corrections were performed on DEMs that were set to the reference DEM elevation over the outcropping bedrock across the Skarvsnes region (yellow region in Supplementary Fig. 6b), where the elevation differences among the various constructed DEMs and reference DEM were assumed to be negligible. We first subtracted the differences between each DEM and the reference DEM over the outcropping bedrock and then applied elevation-dependent, sensor-dependent, and doming/bowl corrections that were estimated via least-squares fitting to the bedrock-corrected DEMs. We excluded the regions across LKI where the slope was steep (> 15°) relative to the reference DEM during the corrections. We also removed any pixels in a given DEM that possessed differences of more than three standard deviations between the corrected DEM and the reference DEM. We then took the standard deviations of the difference between each corrected DEM and the reference DEM to be the uncertainty of that DEM. The elevation uncertainties in the corrected SfM-DEMs, ALOS-DEM, TDX-DEM, and REMA-strip-DEMs were 2.7–5.9, 2.7, 5.0, and 1.0–6.0 m, respectively (Supplementary Table 1).
The lake-surface elevation was derived from the DEMs as the approximately constant elevation value across LKI, such as the 162 m above sea level (a.s.l.) contour in Fig. 1c. We assumed that the elevation uncertainty was equal to the vertical uncertainty in the DEMs. We also used ice-surface elevation measurements from the ICESat-2/ATLAS ATL06 product70 version 5. These consisted of two ICESat-2 tracks that crossed the lake surface on 28 February 2019. The elevation values inside the lake area were averaged to obtain the lake-surface elevation (Supplementary Fig. 7), and the standard deviation was taken to be the measurement uncertainty.
We also estimated the lake-surface elevation using the optical satellite imagery and the 1 October 2017 REMA-DEM (REMAempty). REMAempty was acquired shortly after the lake drainage, when the lake level was low; therefore, it should represent the lake-bed geometry (Supplementary Fig. 8). We assumed that the REMAempty contour lines in the lake area represented the lake extent at an arbitrary elevation, such that the lake-surface elevation was determined by comparing the REMAempty contour lines and surface elevations in the optical satellite images. The elevation of the best-fit contour line was defined as the lake-surface elevation for the date when the satellite image was acquired. We iterated this procedure to determine the upper and lower elevation bounds, which were often different based on the features along the lake margin (Supplementary Fig. 9). The range between the upper and lower bounds was assumed to be the measurement uncertainty. We used the false-color composite images that were constructed from the Landsat 5–8, Sentinel-2, and ASTER imagery (DatalistOpticalSatellite.xlsx) for this procedure.
We calculated the relative lake water volume (dV) from the elevation difference between each DEM and the reference DEM as follows:
where (dhleft( t right)) is the relative height of the lake surface at time (t) relative to the reference DEM, and (A_{lake} left( t right)) is the transient lake area, which is estimated from REMAempty. (overline{dh} left( t right)) is the mean elevation difference relative to the reference DEM over the lake. We also used the lake-surface elevations that were retrieved from the ICESat-2 data to calculate the relative lake volume from pseudo DEMs, whereby the two-dimensional lake-surface elevations from inside the known bounds of LKI were used to infer (dhleft( t right)) and (A_{lake} left( t right)) from REMAempty. The uncertainty in the lake-volume change (sigma_{dV}) is assumed to be:
where ({sigma }_{dhleft(tright)}) is the uncertainty in the elevation difference. Note that we subtracted the lake-volume change owing to the migration of the ice-sheet margin because the LKI–EAIS boundary was recognized.
We investigated the possibility of the LKI ice dam reaching flotation. We used Eq. (1) to test this hypothesis under the assumption that the ice where (Delta z) < 0 can go afloat when the subglacial environment is hydrologically connected to the lake due to the power balance. The local bed elevation was taken from BedMachine43 version 3.
The REMA-DEMs were downloaded using the REMA Strip DEM extent index and REMA Mosaic DEM extent index (https://www.pgc.umn.edu/data/rema/); each URL is summarized in Supplementary Data S2 (DatalistREMA.xlsx). The Landsat images were downloaded from Earth Explorer (https://earthexplorer.usgs.gov/). The Sentinel-2 images were downloaded from https://scihub.copernicus.eu/. The ASTER images were downloaded from https://gbank.gsj.jp/madas/. The SfM-DEMs constructed in this study, lake-surface elevations, and relative volume time series for LKI are available from https://doi.org/10.17592/001.202311202 and https://doi.org/10.17592/001.2023112001, respectively. The TanDEM-X and TerraSAR-X data were downloaded from https://eoweb.dlr.de/egp/. The precipitation data for Syowa Station were downloaded from https://www.data.jma.go.jp/gmd/risk/obsdl/.
The map of the study site (Fig. 1a) was generated using open-source GIS software (QGIS 3.22). The other figures were generated using MATLAB 2023a (https://mathworks.com/products/matlab.html). The MATLAB scripts were written using standard MATLAB functions and functions that are included in M_Map71 version 1.4 m, which is a freely available toolbox.
Bamber, J. L., Westaway, R. M., Marzeion, B. & Wouters, B. The land ice contribution to sea level during the satellite era. Environ. Res. Lett. 13, 063008. https://doi.org/10.1088/1748-9326/aac2f0 (2018).
Article ADS Google Scholar
DeConto, R. M. et al. The Paris climate agreement and future sea-level rise from Antarctica. Nature 593, 83–89. https://doi.org/10.1038/s41586-021-03427-0 (2021).
Article ADS CAS PubMed Google Scholar
Kingslake, J. et al. Widespread movement of meltwater onto and across Antarctic ice shelves. Nature 544, 349–352. https://doi.org/10.1038/nature22049 (2017).
Article ADS CAS PubMed Google Scholar
Trusel, L. D., Pan, Z. & Moussavi, M. Repeated tidally induced hydrofracture of a supraglacial lake at the Amery Ice Shelf grounding zone. Geophys. Res. Lett. 49, e2021GL095661. https://doi.org/10.1029/2021GL095661 (2022).
Article ADS Google Scholar
Willis, M. J., Herried, B. G., Bevis, M. G. & Bell, R. E. Recharge of a subglacial lake by surface meltwater in northeast Greenland. Nature 518, 223–227. https://doi.org/10.1038/nature14116 (2015).
Article ADS CAS PubMed Google Scholar
Stearns, L., Smith, B. & Hamilton, G. Increased flow speed on a large East Antarctic outlet glacier caused by subglacial floods. Nat. Geosci. 1, 827–831. https://doi.org/10.1038/ngeo356 (2008).
Article ADS CAS Google Scholar
Lai, C. Y. et al. Vulnerability of Antarctica’s ice shelves to meltwater-driven fracture. Nature 584, 574–578. https://doi.org/10.1038/s41586-020-2627-8 (2020).
Article ADS CAS PubMed Google Scholar
Bell, R. E. et al. Antarctic surface hydrology and impacts on ice-sheet mass balance. Nat. Clim. Change 8, 1044–1052. https://doi.org/10.1038/s41558-018-0326-3 (2018).
Article ADS Google Scholar
Bell, R. et al. Antarctic ice shelf potentially stabilized by export of meltwater in surface river. Nature 544, 344–348. https://doi.org/10.1038/nature22048 (2017).
Article ADS CAS PubMed Google Scholar
Rignot, E. et al. Accelerated ice discharge from the Antarctic Peninsula following the collapse of Larsen B ice shelf. Geophys. Res. Lett. 31(18), L18401. https://doi.org/10.1029/2004GL020697 (2004).
Article ADS CAS Google Scholar
Scambos, T. A., Bohlander, J. A., Shuman, C. A. & Skvarca, P. Glacier acceleration and thinning after ice shelf collapse in the Larsen B embayment, Antarctica. Geophys. Res. Lett. 31, L18402. https://doi.org/10.1029/2004GL020670 (2004).
Article ADS CAS Google Scholar
Wuite, J. et al. Evolution of surface velocities and ice discharge of Larsen B outlet glaciers from 1995 to 2013. Cryosphere 9(3), 957–969. https://doi.org/10.5194/tc-9-957-2015 (2015).
Article ADS Google Scholar
Wingham, D. et al. Rapid discharge connects Antarctic subglacial lakes. Nature 440, 1033–1036. https://doi.org/10.1038/nature04660 (2006).
Article ADS CAS PubMed Google Scholar
Veh, G. et al. Less extreme and earlier outbursts of ice-dammed lakes since 1900. Nature 614, 701–707. https://doi.org/10.1038/s41586-022-05642-9 (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Lützow, N., Veh, G., & Korup, O. A global database of historic glacier lake outburst floods, Earth Syst. Sci. Data https://doi.org/10.5194/essd-15-2983-2023 (2023).
Article Google Scholar
Veh, G. et al. Trends, breaks, and biases in the frequency of reported glacier lake outburst floods. Earth’s Futur. 10, e2021EF002426. https://doi.org/10.1029/2021EF002426 (2022).
Article ADS Google Scholar
Harrison, S. et al. Climate change and the global pattern of moraine-dammed glacial lake outburst floods. Cryosphere 12, 1195–1209. https://doi.org/10.5194/tc-12-1195-2018 (2018).
Article ADS Google Scholar
Tiberti, R. et al. Food web complexity of high mountain lakes is largely affected by glacial retreat. Ecosystems 23, 1093–1106. https://doi.org/10.1007/s10021-019-00457-8 (2020).
Article Google Scholar
Meerhoff, E., Castro, L. R., Tapia, F. J. & Pérez-Santos, I. Hydrographic and Biological impacts of a glacial lake outburst flood (GLOF) in a Patagonian fjord. Estuaries Coasts. 42, 132–143. https://doi.org/10.1007/s12237-018-0449-9 (2019).
Article CAS Google Scholar
Sommaruga, R. When glaciers and ice sheets melt: Consequences for planktonic organisms. J. Plankton Res. 37(3), 509–518. https://doi.org/10.1093/plankt/fbv027 (2015).
Article PubMed Google Scholar
Glasser, N. et al. Glacial lake drainage in Patagonia (13–8 kyr) and response of the adjacent Pacific Ocean. Sci. Rep. 6, 21064. https://doi.org/10.1038/srep21064 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Thorndycraft, V. R. et al. Glacial lake evolution and Atlantic-Pacific drainage reversals during deglaciation of the Patagonian Ice Sheet. Quat. Sci. Rev. 203, 102–127. https://doi.org/10.1016/j.quascirev.2018.10.036 (2019).
Article ADS Google Scholar
Sawagaki, T., Miura, H. & Iwasaki, S. Discovery of an ice cave in the Yatsude Valley, Langhovde, Dronning Maud Land, East Antarctica. Polar Sci. 2, 287–294. https://doi.org/10.1016/j.po;ar.2008.10.002 (2008).
Article ADS Google Scholar
Kudoh, S., Tanabe, Y., Inoue, T., Imura, S. & Kanda, H. Breaching of a perennial snow dam below Lake Hyoga Ike in the Langhovde region of the Sôya Coast, East Antarctica: Probable effect of disturbance events on the distribution and colonization of flora within/around the lake. Nankyoku Shiryô Antarctic Record. 53(1), 114–122. https://doi.org/10.15094/00009488 (2009).
Article Google Scholar
Faucher, B. et al. Glacial lake outburst floods enhance benthic microbial productivity in perennially ice-covered Lake Untersee (East Antarctica). Commun. Earth Environ. 2, 211. https://doi.org/10.1038/s43247-021-00280-x (2021).
Article ADS Google Scholar
Pryakhina, G. V., Boronina, A. S., Porov, S. V. & Chetverova, A. A. Hydrological studies of lake outbursts in the Antarctic Oases. Russ. Meteorol. Hydrol. Hydrol. Stud. 45(2), 118–123. https://doi.org/10.3103/S1068373920020077 (2020).
Article Google Scholar
Hodgson, D. A., Jordan, T. A., Ross, N., Riley, T. R. & Fretwell, P. T. Drainage and refill of an Antarctic Peninsula subglacial lake reveal an active subglacial hydrological network. Cryosphere 16, 4797–4809. https://doi.org/10.5194/tc-16-4797-2022 (2022).
Article ADS Google Scholar
Pryakhina, G., Chetverova, A., Ryzhova, E. & Grigoreva, S. Formation of a large ice depression on Dålk Glacier (Larsemann Hills, East Antarctica) caused by the rapid drainage of an englacial cavity. J. Glaciol. 67(266), 1121–1136. https://doi.org/10.1017/jog.2021.58 (2021).
Article Google Scholar
Warner, R. C. et al. Rapid formation of an ice doline on Amery Ice Shelf, East Antarctica. Geophys. Res. Lett. 48, e2020GL091095. https://doi.org/10.1029/2020GL091095 (2021).
Article ADS PubMed PubMed Central Google Scholar
Mizuno, T., Shigeoka, H., Ono, F., Mori, Y. & Umezawa, G. Meteorological observations at Syowa Station in 2017 by the 58th Japanese Antarctic Research Expedition. Nankyoku Shiryô (Antarctic Record). 64, 132–177. https://doi.org/10.15094/00015931 (2020).
Article Google Scholar
Suganuma, Y., Tanabe, Y., Katsuki, K., Shibata, D. & Kawamata, M. Report on lake/marine sediment coring from ice surface (JARE-58/59). Nankyoku Shiryô (Antarctic Record). 62, 15–42. https://doi.org/10.15094/00015002 (2018).
Article Google Scholar
Murayama, H. General characteristics of the Antarctic lakes near Syowa station. Nankyoku Shiryô (Antarctic Record). 58, 43–62. https://doi.org/10.15094/00007903 (1977).
Article Google Scholar
Kawamata, M., Doi, K., Sawagaki, T. & Suganuma, Y. Generation and utility of a digital surface model based on archived aerial photographs taken by the Japanese Antarctic research expedition. Geogr. Rev. Jpn. Ser. A. 94(1), 1–16. https://doi.org/10.4157/grj.94.1 (2021).
Article Google Scholar
Howat, I. M., Porter, C., Smith, B. E., Noh, M.-J. & Morin, P. The Reference elevation model of Antarctica. Cryosphere 13, 665–674. https://doi.org/10.5194/tc-13-665-2019 (2019).
Article ADS Google Scholar
Howat, I, et al., The reference elevation model of Antarctica—Strips, Version 4.1. https://doi.org/10.7910/DVN/X7NDNY, Harvard Dataverse, V1, (last access: 10 March 2023). (2022).
Grinsted, A., Hvidberg, C. S., Campos, N. & Dahl-Jensen, D. Periodic outburst floods from an ice-dammed lake in East Greenland. Sci. Rep. 7, 9966. https://doi.org/10.1038/s41598-017-07960-9 (2017).
Article ADS CAS PubMed PubMed Central Google Scholar
Russell, A., Carrivick, J. L., Ingeman-Neilsen, T., Yde, J. C. & Williams, M. A new cycle of jökulhlaups at Russell Glacier, Kangerlussuaq, West Greenland. J. Glaciol. 57(202), 238–246. https://doi.org/10.3189/002214311796405997 (2011).
Article ADS Google Scholar
Meyer, C. & Shuler, T. V. Breaching of an ice dam at Qorlortossup tasia, south Greenland. Ann. Glaciol. 42, 297–302. https://doi.org/10.3189/172756405781812989 (2005).
Article ADS Google Scholar
Sugiyama, S., Bauder, A., Huss, M., Riesen, P. & Funk, M. Triggering and drainage mechanisms of the 2004 glacier-dammed lake outburst in Gornergletscher, Switzerland. J. Geophys. Res. 113, F04019. https://doi.org/10.1029/2007JF000920 (2008).
Article ADS Google Scholar
Bazai, N. Z. et al. Increasing glacial lake outburst flood hazard in response to surge glaciers in the Karakoram. Earth-Sci. Rev. 212, 103432. https://doi.org/10.1016/j.earscirev.2020.103432 (2021).
Article Google Scholar
Clague, J. J. & O’Connor, J. E. Glacier-related outburst floods. Snow and Ice-Related Hazards, Risk, and Disasters (Second edition). Hazards Disasters Ser. https://doi.org/10.1016/B978-0-12-817129-5.00019-6 (2021).
Article Google Scholar
Carrivick, J. L. et al. Ice-dammed lake drainage evolution at Russell Glacier, West Greenland. Front. Earth Sci. 5, 100. https://doi.org/10.3389/feart.2017.00100 (2017).
Article ADS Google Scholar
Morlighem, M. et al. Deep glacial troughs and stabilizing ridges unveiled beneath the margins, of the Antarctic ice sheet. Nat. Geosci. 13, 132–137. https://doi.org/10.1038/s41561-019-0510-8 (2020).
Article ADS CAS Google Scholar
Nye, J. F. Water flow in glaciers: Jökulhlaups, tunnels, and veins. J. Glaciol. 17(76), 181–207. https://doi.org/10.3189/S002214300001354X (1976).
Article ADS Google Scholar
Huss, M., Bauder, A., Werder, M., Funk, M. & Hock, R. Glacier-dammed lake outburst events of Gornersee, Switzerland. J. Glaciol. 53(181), 189–200. https://doi.org/10.3189/172756507782202784 (2007).
Article ADS Google Scholar
Uemura, T., Taniguchi, M. & Shibuya, K. Submarine groundwater discharge in Lützow-Holm Bay, Antarctica. Geophys. Res. Lett. 38(8), L08402. https://doi.org/10.1029/2010GL046394 (2011).
Article ADS Google Scholar
Sawagaki, T. & Hirakawa, K. Hydrostatic investigations on subglacial meltwater: Implications for the formation of streamlined bedforms and subglacial lakes, East Antarctica. Polar Geosci. 15, 123–147. https://doi.org/10.15094/00003114 (2002).
Article Google Scholar
Fujita, S. et al. Radar diagnosis of the subglacial conditions in Dronning Maud Land East Antarctica. Cryosphere 6, 1203–1219. https://doi.org/10.5194/tc-6-1203-2012 (2012).
Article ADS Google Scholar
Rickman, R. L. & Rosenkrans, D. S. Hydrologic conditions and hazards in the Kennicott River Basin, Wrangell-St. Elias national park preserve. Alaska Water-Resour. Investig. Rep. https://doi.org/10.3133/wri964296 (1997).
Article Google Scholar
Rick, B., McGrath, D., McCoy, S. W. & Armstrong, W. H. Unchanged frequency and decreasing magnitude of outbursts from ice-dammed lakes in Alaska. Nat. Commun. 14, 6138. https://doi.org/10.1038/s41467-023-41794-6 (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Vignon, É., Roussel, M.-L., Gorodetskaya, I. V., Genthon, C. & Berne, A. Present and future of rainfall in Antarctica. Geophys. Res. Lett. 48, e2020GL092281. https://doi.org/10.1029/2020GL092281 (2021).
Article ADS Google Scholar
Turner, J. et al. An extreme high temperature event in coastal East Antarctica associated with an atmospheric river and record summer downslope winds. Geophys. Res. Lett. 49, e2021GL097108. https://doi.org/10.1029/2021GL097108 (2022).
Article ADS Google Scholar
Harper, J., Meierbachtol, T., Humphrey, N., Saito, J. & Stansberry, A. Generation and fate of basal meltwater during winter, western Greenland Ice Sheet. Cryosphere 15, 5409–5421. https://doi.org/10.5194/tc-15-5409-2021 (2021).
Article ADS Google Scholar
Pitcher, L. H. et al. Direct observation of winter meltwater drainage from the Greenland Ice Sheet. Geophys. Res. Lett. 47, e2019GL086521. https://doi.org/10.1029/2019GL086521 (2020).
Article ADS Google Scholar
Carrivick, J. L. & Tweed, F. S. A review of glacier outburst floods in Iceland and Greenland with a megaflood perspective. Earth-Sci. Rev. 196, 102876. https://doi.org/10.1016/j.earscirev.2019.102876 (2019).
Article Google Scholar
Tewari, K., Mishra, S. K., Salunke, P. & Dewan, A. Future projections of temperature and precipitation for Antarctica. Environ. Res. Lett. 17, 014029. https://doi.org/10.1088/1748-9326/ac43e2 (2022).
Article ADS Google Scholar
Beaumet, J. et al. Significant additional Antarctic warming in atmospheric bias-corrected ARPEGE projections with respect to control run. Cryosphere 15(8), 3615–3635. https://doi.org/10.5194/tc-15-3615-2021 (2021).
Article ADS Google Scholar
How, P. et al. Greenland-wide inventory of ice marginal lakes using a multi-method approach. Sci. Rep. 11, 4481. https://doi.org/10.1038/s41598-021-83509-1 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Liang, Q. et al. Filling and drainage of a subglacial lake beneath the Flade Isblink ice cap, northeast Greenland. Cryosphere 16, 2671–2681. https://doi.org/10.5194/tc-16-2671-2022 (2022).
Article ADS Google Scholar
Maier, N., Andersen, J. K., Mouginot, J., Gimbert, F. & Gagliardini, O. Wintertime supraglacial lake drainage cascade triggers large-scale ice flow response in Greenland. Geophys. Res. Lett. 50, e2022GL102251. https://doi.org/10.1029/2022GL102251 (2023).
Article ADS Google Scholar
Yang, R. et al. Glacier surface speed variations on the Kenai Peninsula, Alaska, 2014–2019. JGR Earth Surf. 127, e2022JF006599. https://doi.org/10.1029/2022JF006599 (2022).
Article ADS Google Scholar
Furian, W., Maussion, F. & Schneider, C. Projected 21st-century glacial lake evolution in high mountain Asia. Front. Earth Sci. 10, 821798. https://doi.org/10.3389/feart.2022.821798 (2022).
Article Google Scholar
Nakayama, Y., Cai, C. & Seroussi, H. Impact of subglacial freshwater discharge on pine island ice shelf. Geophys. Res. Lett. 48(18), e2021GL093923. https://doi.org/10.1029/2021GL093923 (2021).
Article ADS Google Scholar
Wei, W. et al. Getz Ice Shelf melt enhanced by freshwater discharge from beneath the West Antarctic Ice Sheet. Cryosphere 14(4), 1399–1408. https://doi.org/10.5194/tc-14-1399-2020 (2020).
Article ADS Google Scholar
Shean, D. E. et al. An automated, open-source pipeline for mass production of digital elevation models (DEMs) from very-high-resolution commercial stereo satellite imagery. ISPRS J. Photogramm. Remote Sens. 116, 101–117. https://doi.org/10.1016/j.isprsjprs.2016.03.012 (2016).
Article ADS Google Scholar
Noh, M.-J. & Howat, I. M. Automated stereo-photogrammetric DEM generation at high latitudes: Surface Extraction with TIN-based Search-space Minimization (SETSM) validation and demonstration over glaciated regions. GISci. Remote Sens. 52(2), 198–217. https://doi.org/10.1080/15481603.2015.1008621 (2015).
Article Google Scholar
Nuth, C. & Kääb, A. Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change. Cryosphere 5, 271–290. https://doi.org/10.5194/tc-5-271-2011 (2011).
Article ADS Google Scholar
Sanz-Ablanedo, E., Chandkker, J. H., Ballesteros-Pérez, P. & Rodríguez-Pérez, R. Reducing systematic dome errors in digital elevation models through better UAV flight design. Earth Surf. Process. Landf. 45, 2134–2147. https://doi.org/10.1002/esp.4871 (2020).
Article ADS Google Scholar
Girod, L. & Filhol, S. Absolute coregistration and doming correction in adverse conditions, or how to retrieve snow depth from drone flights. ISPRS Ann Photogram. Remote Sens. Spat. Inf. Sci. 3, 375–379. https://doi.org/10.5194/isprs-annals-V-3-2020-375-2020 (2020).
Article Google Scholar
Smith, B., S. et al. ATLAS/ICESat-2 L3A Land Ice Height, Version 5. Boulder, Colorado USA. NASA national snow and ice data center distributed active archive center. https://doi.org/10.5067/ATLAS/ATL06.005 (last access: 10 March 2023). (2021)
Pawlowicz, R. “M_Map: A mapping package for MATLAB”, version 1.4m, [Computer software], available online at https://www.eoas.ubc.ca/rich/map.html (2020).
Download references
The authors thank Sakae Kudoh for providing insightful information on the lake limnological research in the region and Shin Sugiyama for providing the ALOS/PRISM images and assisting in the photogrammetric analysis of the aerial photographs. The PALSAR-2 SLC data are shared among the ALOS4ice cryospheric science team and are provided under a cooperative research contract (JAXA/EORA3#ER3A2N060); the ownership of the PALSAR-2 data belongs to JAXA. We used TerraSAR-X/TanDEM-X data that were provided by the German Aerospace Centre (DLR) within the framework of an approved proposal (ID XTI_GLAC6706). The authors extend a special thanks to the field members who took the aerial photographs during JARE-1/6/10/12/24/32/34/41 and the field members who conducted field observations in Skarvsnes during JARE59. We thank scientific editor Shiqiao Zhou and five anonymous reviewers for providing constructive comments that greatly improved the manuscript. This work was supported by ROIS-DS-JOINT (023RP2023), and S.H. was supported by JST SPRING, Grant Number JPMJSP2119 (2021–2022).
Institute of Low Temperature Science, Hokkaido University, Sapporo, Japan
Shuntaro Hata
Creative Research Institution, Hokkaido University, Sapporo, Japan
Shuntaro Hata
Civil Engineering Research Institute for Cold Region, Public Works Research Institute, Sapporo, Japan
Moto Kawamata
National Institute of Polar Research, Tokyo, Japan
Koichiro Doi
The Graduate University for Advanced Studies (SOKENDAI), Hayama, Japan
Koichiro Doi
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
S.H., M.K., and K.D. designed the study. S.H. processed the satellite and lake-surface elevation data, analyzed the data, produced the figures, and wrote the initial manuscript. S.H. constructed the SfM-DEMs under the guidance of M.K. K.D. processed the TanDEM-X and TerraSAR-X data, and produced the TDX-DEM. All of the authors contributed to the discussions and final manuscript.
Correspondence to Shuntaro Hata.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Hata, S., Kawamata, M. & Doi, K. Outbursts from an ice-marginal lake in Antarctica in 1969–1971 and 2017, revealed by aerial photographs and satellite data. Sci Rep 13, 20619 (2023). https://doi.org/10.1038/s41598-023-47522-w
Download citation
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-47522-w
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.