Search

Arctic lakes are sensitive to climate change, and the timing and duration of ice presence and absence (i.e., ice phenology) on the lake surface can be used as a climate indicator. In this study of Linnévatnet, one of the largest lakes on Svalbard, we compare inferences of lake ice duration from satellite data with continuously monitored lake water temperature and photographs from automatic cameras. Visible surface reflectance data from the moderate-resolution imaging spectroradiometer (MODIS) were used to observe the change in the lake-wide mean surface reflectance of Linnévatnet from 2003–2019, and smoothing splines were applied to the to determine the date of summer ice-off (also called “break-up end”—BUE). Similarly, BUE and fall ice-on (or “freeze-up end”—FUE) were determined from lake-wide mean time series of Sentinel-1 microwave backscatter from 2014–2019. Overall, the ice timing dates identified from the satellite observations agree well with the in-situ observations (RMSE values of approximately 2–7 days for BUE and FUE, depending on the method and in-situ dataset), lending confidence to the accuracy of remote sensing of lake ice phenology in remote Arctic regions. Our observations of Linnévatnet indicate that BUE dates do not have a significant trend, while FUE dates have been occurring approximately 1.5 days later per year during the study period. These results support an overall decrease in annual duration of lake ice cover in this part of Svalbard.

As an important part of the cryosphere, lake ice is a sensitive indicator of climate change. Remote sensing technology can quickly and accurately monitor the process of its formation and decay, among which Moderate Resolution Imaging Spectroradiometer (MODIS) images are the most widely used data in the remote sensing monitoring of lake ice. The reasonable selection of monitoring methods is of great significance to grasp the dynamic process and response to climate change of lake ice. In this study, five commonly used remote sensing monitoring methods of lake ice based on MODIS MOD09GA data, including the single band threshold method (SBT), reflectance difference threshold method (RDT), normalized difference snow index method (NDSI), modified normalized difference snow index method (MNDSI) and lake ice index method (LII), were selected to compare their accuracies in extracting lake ice extent by combining them with four evaluation metrics of accuracy, precision, recall and mean intersection over union (MIoU). In addition, the ability of the high-precision LII method for extracting long time series lake ice phenology and its applicability to multiple types of lakes were verified. The results showed that compared with the NDSI method, the other four methods more easily distinguished between lake ice and lake water by setting thresholds. The SBT method and the RDT method had better extraction effects in the freezing process and the melting process, respectively. Compared with the NDSI and MNDSI methods, the LII method showed a significant improvement in lake ice extraction over the entire freeze–thaw cycle, with the smallest mean monitoring error of 1.53% for the percentage of lake ice area in different periods. Meanwhile, the LII method can be used to determine long term lake ice phenology dates and had good performance in extracting lake ice for different types of lakes on the Qinghai–Tibet Plateau with the optimal threshold interval of 0.05~0.07, which can be used for lake ice monitoring and long-term phenological studies in this region.

Ice surface albedo is an important factor in various optical remote sensing technologies used to determine the distribution of snow or melt water on the ice, and to judge the formation or melting of lake ice in winter, especially in cold and arid areas. In this study, field measurements were conducted at Wuliangsuhai Lake, a typical lake in the semi-arid cold area of China, to investigate the diurnal variation of the ice surface albedo. Observations showed that the diurnal variations of the ice surface albedo exhibit bimodal characteristics with peaks occurring after sunrise and before sunset. The curve of ice surface albedo with time is affected by weather conditions. The first peak occurs later on cloudy days compared with sunny days, whereas the second peak appears earlier on cloudy days. Four probability density distribution functions-Laplace, Gauss, Gumbel, and Cauchy-were combined linearly to model the daily variation of the lake ice albedo on a sunny day. The simulations of diurnal variation in the albedo during the period from sunrise to sunset with a solar altitude angle higher than 5 degrees indicate that the Laplace combination is the optimal statistical model. The Laplace combination can not only describe the bimodal characteristic of the diurnal albedo cycle when the solar altitude angle is higher than 5 degrees, but also reflect the U-shaped distribution of the diurnal albedo as the solar altitude angle exceeds 15 degrees. The scale of the model is about half the length of the day, and the position of the two peaks is closely related to the moment of sunrise, which reflects the asymmetry of the two peaks of the ice surface albedo. This study provides a basis for the development of parameterization schemes of diurnal variation of lake ice albedo in semi-arid cold regions.

Seepage of geological methane through sediments of Arctic lakes might contribute conceivably to the atmospheric methane budget. However, the abundance and precise locations of such seeps are poorly quantified. For Lake Neyto, one of the largest lakes on the Yamal Peninsula in Northwestern Siberia, temporally expanding regions of anomalously low backscatter in C-band SAR imagery acquired in late winter and spring have been suggested to be related to seepage of methane from hydrocarbon reservoirs. However, this hypothesis has not been verified using in-situ observations so far. Similar anomalies have also been identified for other lakes on Yamal, but it is still uncertain whether or how many of them are related to methane seepage. This study aimed to document similar lake ice backscatter anomalies on a regional scale over four study regions (the Yamal Peninsula and Tazovskiy Peninsulas, the Lena Delta in Russia, the National Petroleum Reserve Alaska) during different years using a time series based approach on Google Earth Engine (GEE) that quantifies changes of σ0 from the Sentinel-1 C-band SAR sensor over time. An algorithm for assessing the coverage that takes the number of acquisitions and maximum time between acquisitions into account is presented, and differences between the main operating modes of Sentinel-1 are evaluated. Results show that better coverage can be achieved in extra wide swath (EW) mode, but interferometric wide swath (IW) mode data could be useful for smaller study areas and to substantiate EW results. A classification of anomalies on Lake Neyto from EW Δσ0 images derived from GEE showed good agreement with the classification presented in a previous study. Automatic threshold-based per-lake counting of years where anomalies occurred was tested, but a number of issues related to this approach were identified. For example, effects of late grounding of the ice and anomalies potentially related to methane emissions could not be separated efficiently. Visualizations of Δσ0 images likely reflect the temporal expansions of anomalies and are expected to be particularly useful for identifying target areas for future field-based research. Characteristic anomalies that clearly resemble the ones observed for Lake Neyto could be identified solely visually in the Yamal and Tazovskiy study regions. All data and algorithms produced in the framework of this study are openly provided to the scientific community for future studies and might potentially aid our understanding of geological lake seepage upon the progression of related field-based studies and corresponding evaluations of formation hypotheses.

Lake ice is an important component in understanding the local climate as changes in temperature have an impact on the timing of key ice phenology events. In recent years, there has been a decline in the in-situ monitoring of lake ice events in Canada and microwave remote sensing imagery from synthetic aperture radar (SAR) is more widely used due to the high spatial resolution and response of backscatter to the freezing and melting of the ice surface. RADARSAT-2 imagery was used to develop a threshold-based method for determining lake ice events for mid-latitude lakes in Central Ontario from 2008 to 2017. Estimated lake ice phenology events are validated with ground-based observations and are compared against the Moderate Resolution Imaging Spectroradiometer (MODIS band 2). The threshold-based method was found to accurately identify 12 out of 17 freeze events and 13 out of 17 melt events from 2015–2017 when compared to ground-based observations. Mean absolute errors for freeze events ranged from 2.5 to 10.0 days when compared to MODIS imagery while the mean absolute error for water clear of ice (WCI) ranged from 1.5 to 7.1 days. The method is important for the study of mid-latitude lake ice due to its unique success in detecting multiple freeze and melting events throughout the ice season.

Lake ice thickness is a sensitive indicator of climate change largely through its dependency on near-surface air temperature and on-ice snow mass (depth and density). Monitoring of the seasonal variations and trends in ice thickness is also important for the operation of winter ice roads that northern communities rely on for the movement of goods as well as for cultural and leisure activities (e.g., snowmobiling). Therefore, consistent measurements of ice thickness over lakes is important; however, field measurements tend to be sparse in both space and time in many northern countries. Here, we present an application of L-band frequency Global Navigation Satellite System (GNSS) Interferometric Reflectometry (GNSS-IR) for the estimation of lake ice thickness. The proof of concept is demonstrated through the analysis of Signal-to-Noise Ratio (SNR) time series extracted from Global Positioning System (GPS) constellation L1 band raw data acquired between 8 and 22 March (2017 and 2019) at 14 lake ice sites located in the Northwest Territories, Canada. Dominant frequencies are extracted using Least Squares Harmonic Estimation (LS-HE) for the retrieval of ice thickness. Estimates compare favorably with in-situ measurements (mean absolute error = 0.05 m, mean bias error = −0.01 m, and root mean square error = 0.07 m). These results point to the potential of GPS/GNSS-IR as a complementary tool to traditional field measurements for obtaining consistent ice thickness estimates at many lake locations, given the relatively low cost of GNSS antennas/receivers.

Changes to ice cover on lakes throughout the northern landscape has been established as an indicator of climate change and variability, expected to have implications for both human and environmental systems. Monitoring lake ice cover is also required to enable more reliable weather forecasting across lake-rich northern latitudes. Currently, the Canadian Ice Service (CIS) monitors lakes using synthetic aperture radar (SAR) and optical imagery through visual interpretation, with total lake ice cover reported weekly as a fraction out of ten. An automated method of classification would allow for more detailed records to be delivered operationally. In this research, we present an automatic ice-mapping approach which integrates unsupervised segmentation from the Iterative Region Growing using Semantics (IRGS) algorithm with supervised random forest (RF) labeling. IRGS first locally segments homogeneous regions in an image, then merges similar regions into classes across the entire scene. Recently, these output regions were manually labeled by the user to generate ice maps, or were labeled using a Support Vector Machine (SVM) classifier. Here, three labeling methods (Manual, SVM, and RF) are applied after IRGS segmentation to perform ice-water classification on 36 RADARSAT-2 scenes of Great Bear Lake (Canada). SVM and RF classifiers are also tested without integration with IRGS. An accuracy assessment has been performed on the results, comparing outcomes with author-generated reference data, as well as the reported ice fraction from CIS. The IRGS-RF average classification accuracy for this dataset is 95.8%, demonstrating the potential of this automated method to provide detailed and reliable lake ice cover information operationally.

Satellite remote sensing has transformed our understanding of Earth processes. One component of the Earth system where large uncertainties remain are Arctic and boreal freshwater lakes. With only short periods of open water due to annual ice cover, lake productivity in these regions is extremely sensitive to warming induced changes in ice cover. At the same time, productivity dynamics in these lakes vary enormously, even over short distances, making it difficult to understand these potential changes. A major impediment to an improved understanding of lake dynamics has been sparsely distributed field measurements, in large part due to the complexity and expense of conducting scientific research in remote northern latitudes. This project overcomes that hurdle by using a new set of ‘eyes in the sky’, the Planet Labs CubeSat fleet, to observe 35 lakes across 3 different arctic-boreal ecoregions in western North America. We extract time series of lake reflectance to identify ice-out and green-up across three years (2017–2019). We find that lakes with later ice-out have significantly faster green-ups. Our results also show ice-out varies latitudinally by 38 days from south to north, but only varies across years by ~9 days. In contrast, green-up varied between years by 22 days in addition to showing significant spatial variability. We compare PlanetScope to Sentinel-2 data and independently validate our ice-out estimates, finding an ice-out mean absolute difference (MAD) ~9 days. This study demonstrates the potential of using CubeSat imagery to monitor the timing and magnitude of ice-off and green-up at high spatiotemporal resolution., Remote Sensing, 13 (13), ISSN:2072-4292

In the Tibetan Plateau (TP), the changes of lake ice phenology not only reflect regional climate change, but also impose substantial ecohydrological impacts on the local environment. Due to the limitation of ground observation, remote sensing has been used as an alternative tool to investigate recent changes of lake ice phenology. However, uncertainties exist in the remotely sensed lake ice phenology owing to both the data and methods used. In this paper, three different remotely sensed datasets are used to investigate the lake ice phenology variation in the past decade across the Tibetan Plateau, with the consideration of the underlying uncertainties. The remotely sensed data used include reflectance data, snow product, and land surface temperature (LST) data of moderate resolution imaging spectroradiometer (MODIS). The uncertainties of the three methods based on the corresponding data are assessed using the triple collocation approach. Comparatively, it is found that the method based on reflectance data outperforms the other two methods. The three methods are more consistent in determining the thawing dates rather than the freezing dates of lake ice. It is consistently shown by the three methods that the ice-covering duration in the northern part of the TP lasts longer than that in the south. Though there is no general trend of lake ice phenology across the TP for the period of 2000–2015, the warmer climate and stronger wind have led to the earlier break-up of lake ice.

As one of the most sensitive indicators of global climate change, seasonal ice-covered lakes are attracting gaining attention worldwide. As a large seasonal ice-covered lake located in Northern Mongolia, Lake Khovsgol not only provides important freshwater resources for the local population but also serves as a means of water transportation in summer and an important land-based activity for residents in winter. In this study, we used the sub-pixel offset technique with multi-temporal Sentinel-2 optical images to estimate the time-series displacement of lake ice in Lake Khovsgol from 7 December 2020 to 17 June 2021. With the processing of 112 Sentinel-2 images, we obtained 27 pairs of displacement results at intervals of 5, 10, and 15 days. These lake ice movement results covered three stages from ice-on to ice-off. The first stage was the lake ice growth period, which lasted 26 days from 7 December 2020 to 3 January 2021. Ice formation started from the south and extended northward, with a displacement of up to 10 m in 5 days. The second stage was the active phase of the ice cover, which took place from 3 January 2021 to 18 April 2021. Maximum displacement values reached 12 m in the east and 11 m in the north among all observations. The value of the lake ice movement in the north–south direction (NS) was found to be larger than in the east–west direction (EW). The third stage was the melting period, which closed on 17 June 2021. In comparison to the freezing date of November in past years, our results demonstrate the ice-on date of Lake Khovsgol has been delayed to December, suggesting a possible reason that the seasonal ice-covered lake located at the middle latitude has been affected by global warming. In addition, the lake ice movement of our results can reveal the regional climate characteristic. This study is one of the few cases to reveal the distribution characteristics and changing trends of lake ice on the Mongolia Plateau, providing a rare reference for lake ice research in this region.

An automated ice-mapping algorithm has been developed and evaluated using data from the GOES-13 imager. The approach includes cloud-free image compositing as well as image classification using spectral criteria. The algorithm uses an alternative snow index to the Normalized Difference Snow Index (NDSI). The GOES-13 imager does not have a 1.6 µm band, a requirement for NDSI; however, the newly proposed Mid-Infrared Sea and Lake Ice Index (MISI) incorporates the reflective component of the 3.9 µm or mid-infrared (MIR) band, which the GOES-13 imager does operate. Incorporating MISI into a sea or lake ice mapping algorithm allows for mapping of thin or broken ice with no snow cover (nilas, frazil ice) and thicker ice with snow cover to a degree of confidence that is comparable to other ice mapping products. The proposed index has been applied over the Great Lakes region and qualitatively compared to the Interactive Multi-sensor Snow and Ice Mapping System (IMS), the National Ice Center ice concentration maps and MODIS snow cover products. The application of MISI may open additional possibilities in climate research using historical GOES imagery. Furthermore, MISI may be used in addition to the current NDSI in ice identification to build more robust ice-mapping algorithms for the next generation GOES satellites.

Lake ice is a significant component of the cryosphere due to its large spatial coverage in high-latitude regions during the winter months. The Laurentian Great Lakes are the world’s largest supply of freshwater and their ice cover has a major impact on regional weather and climate, ship navigation, and public safety. Ice experts at the Canadian Ice Service (CIS) have been manually producing operational Great Lakes image analysis charts based on visual interpretation of the synthetic aperture radar (SAR) images. In that regard, we have investigated the performance of the semi-automated segmentation algorithm “glocal„ Iterative Region Growing with Semantics (IRGS) for lake ice classification using dual polarized RADARSAT-2 imagery acquired over Lake Erie. Analysis of various case studies indicated that the “glocal„ IRGS algorithm could provide a reliable ice-water classification using dual polarized images with a high overall accuracy of 90.4%. However, lake ice types that are based on stage of development were not effectively identified due to the ambiguous relation between backscatter and ice types. The slight improvement of using dual-pol as opposed to single-pol images for ice-water discrimination was also demonstrated.

With the increasing availability of SAR imagery in recent years, more research is being conducted using deep learning (DL) for the classification of ice and open water; however, ice and open water classification using conventional DL methods such as convolutional neural networks (CNNs) is not yet accurate enough to replace manual analysis for operational ice chart mapping. Understanding the uncertainties associated with CNN model predictions can help to quantify errors and, therefore, guide efforts on potential enhancements using more–advanced DL models and/or synergistic approaches. This paper evaluates an approach for estimating the aleatoric uncertainty [a measure used to identify the noise inherent in data] of CNN probabilities to map ice and open water with a custom loss function applied to RADARSAT–2 HH and HV observations. The images were acquired during the 2014 ice season of Lake Erie and Lake Ontario, two of the five Laurentian Great Lakes of North America. Operational image analysis charts from the Canadian Ice Service (CIS), which are based on visual interpretation of SAR imagery, are used to provide training and testing labels for the CNN model and to evaluate the accuracy of the model predictions. Bathymetry, as a variable that has an impact on the ice regime of lakes, was also incorporated during model training in supplementary experiments. Adding aleatoric loss and bathymetry information improved the accuracy of mapping water and ice. Results are evaluated quantitatively (accuracy metrics) and qualitatively (visual comparisons). Ice and open water scores were improved in some sections of the lakes by using aleatoric loss and including bathymetry. In Lake Erie, the ice score was improved by ∼2 on average in the shallow near–shore zone as a result of better mapping of dark ice (low backscatter) in the western basin. As for Lake Ontario, the open water score was improved by ∼6 on average in the deepest profundal off–shore zone.

The timing of lake ice-off regulates biotic and abiotic processes in Arctic ecosystems. Due to the coarse spatial and temporal resolution of available satellite data, previous studies mainly focused on lake-scale investigations of melting/freezing, hindering the detection of subtle patterns within heterogeneous landscapes. To fill this knowledge gap, we developed a new approach for fine-resolution mapping of Pan-Arctic lake ice-off phenology. Using the Scene Classification Layer data derived from dense Sentinel-2 time series images, we estimated the pixel-by-pixel ice break-up end date information by seeking the transition time point when the pixel is completely free of ice. Applying this approach on the Google Earth Engine platform, we mapped the spatial distribution of the break-up end date for 45,532 lakes across the entire Arctic (except for Greenland) for the year 2019. The evaluation results suggested that our estimations matched well with both in situ measurements and an existing lake ice phenology product. Based on the generated map, we estimated that the average break-up end time of Pan-Arctic lakes is 172 ± 13.4 (measured in day of year) for the year 2019. The mapped lake ice-off phenology exhibits a latitudinal gradient, with a linear slope of 1.02 days per degree from 55°N onward. We also demonstrated the importance of lake and landscape characteristics in affecting spring lake ice melting. The proposed approach offers new possibilities for monitoring the seasonal Arctic lake ice freeze–thaw cycle, benefiting the ongoing efforts of combating and adapting to climate change.

Lake ice, one of the most direct lake physical characteristics affected by climate change, can reflect small-scale environmental changes caused by the atmosphere and hydrology, as well as large-scale climate changes such as global warming. This study uses National Oceanic and Atmospheric Administration, Advanced Very High Resolution Radiometer (NOAA AVHRR), MOD09GQ surface reflectance products, and Landsat surface reflectance Tier 1 products, which comprehensively used RS and GIS technology to study lake ice phenology (LIP) and changes in Qinghai Lake. Over the past 38 years, freeze-up start and freeze-up end dates were gradually delayed by a rate of 0.16 d/a and 0.19 d/a, respectively, with a total delay by 6.08 d and 7.22 d. The dates of break-up start and break-up end showed advancing trends by −0.36 d/a and −0.42 d/a, respectively, which shifted them earlier by 13.68 d and 15.96 d. Overall, ice coverage duration, freeze duration, and complete freeze duration showed decreasing trends of −0.58 d/a, −0.60 d/a, and −0.52 d/a, respectively, and overall decreased by 22.04 d, 22.81 d, and 9.76 d between 1980 and 2018. The spatial pattern in the freeze–thaw of Qinghai Lake can be divided into two areas; the west of the lake area has similar spatial patterns of freezing and ablation, while, in the east of the lake area, freezing and ablation patterns are opposite. Climate factors were closely related to LIP, especially the accumulated freezing degree-day (AFDD) from October to April of the following year. Furthermore, freeze-up start time was more sensitive to changes in wind speed and precipitation.

The lake ice phenology variations are vital for the land–surface–water cycle. Qinghai Lake is experiencing amplified warming under climate change. Based on the MODIS imagery, the spatio-temporal dynamics of the ice phenology of Qinghai Lake were analyzed using machine learning during the 2000/2001 to 2019/2020 ice season, and cloud gap-filling procedures were applied to reconstruct the result. The results showed that the overall accuracy of the water–ice classification by random forest and cloud gap-filling procedures was 98.36% and 92.56%, respectively. The annual spatial distribution of the freeze-up and break-up dates ranged primarily from DOY 330 to 397 and from DOY 70 to 116. Meanwhile, the decrease rates of freeze-up duration (DFU), full ice cover duration (DFI), and ice cover duration (DI) were 0.37, 0.34, and 0.13 days/yr., respectively, and the duration was shortened by 7.4, 6.8, and 2.6 days over the past 20 years. The increased rate of break-up duration (DBU) was 0.58 days/yr. and the duration was lengthened by 11.6 days. Furthermore, the increase in temperature resulted in an increase in precipitation after two years, the increase in precipitation resulted in the increase in DBU and decrease in DFU in corresponding years, and decreased DI and DFI after one year.

Ice phenology data of 22 large lakes of the Northern Hemisphere for 40 years (1979–2018) have been retrieved from passive microwave remote sensing brightness temperature (Tb). The results were compared with site-observation data and visual interpretation from Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectivity products images (MOD09GA). The mean absolute errors of four lake ice phenology parameters, including freeze-up start date (FUS), freeze-up end date (FUE), break-up start date (BUS), and break-up end date (BUE) against MODIS-derived ice phenology were 2.50, 2.33, 1.98, and 3.27 days, respectively. The long-term variation in lake ice phenology indicates that FUS and FUE are delayed; BUS and BUE are earlier; ice duration (ID) and complete ice duration (CID) have a general decreasing trend. The average change rates of FUS, FUE, BUS, BUE, ID, and CID of lakes in this study from 1979 to 2018 were 0.23, 0.23, −0.17, −0.33, −0.67, and −0.48 days/year, respectively. Air temperature and latitude are two dominant driving factors of lake ice phenology. Lake ice phenology for the period 2021–2100 was predicted by the relationship between ice phenology and air temperature for each lake. Compared with lake ice phenology changes from 1990 to 2010, FUS is projected to be delayed by 3.1 days and 11.8 days under Representative Concentration Pathways (RCPs) 2.6 and 8.5 scenarios, respectively; BUS is projected to be earlier by 3.3 days and 10.7 days, respectively; and ice duration from 2080 to 2100 will decrease by 6.5 days and 21.9 days, respectively.

Most high-altitude lakes are more sensitive to global warming than the regional atmosphere. However, most existing climate models produce unrealistic surface temperatures on the Tibetan Plateau (TP) lakes, and few studies have focused on the influence of ice surface albedo on high-altitude lakes. Based on field albedo measurements, moderate resolution imaging spectrometer (MODIS) albedo products and numerical simulation, this study evaluates the ice albedo parameterization schemes in existing lake models and investigates the characteristics of the ice surface albedo in six typical TP lakes, as well as the influence of ice albedo error in the FLake model. Compared with observations, several ice albedo schemes all clearly overestimate the lake ice albedo by 0.26 to 0.66, while the average bias of MODIS albedo products is only 0.07. The MODIS-observed albedo of a snow-covered lake varies with the snow proportion, and the lake surface albedo in a snow-free state is approximately 0.15 during the frozen period. The MODIS-observed ice surface (snow-free) albedos are concentrated within the ranges of 0.14–0.16, 0.08–0.10 and 0.10–0.12 in Aksai Chin Lake, Nam Co Lake and Ngoring Lake, respectively. The simulated lake surface temperature is sensitive to variations in lake ice albedo especially in the spring and winter.

Thermokarst lakes and ponds are major elements of permafrost landscapes, occupying up to 40% of the land area in some Arctic regions. Shallow lakes freeze to the bed, thus preventing permafrost thaw underneath them and limiting the length of the period with greenhouse gas production in the unfrozen lake sediments. Radar remote sensing permits to distinguish lakes with bedfast ice due to the difference in backscatter intensities from bedfast and floating ice. This study investigates the potential of a unique time series of three-year repeat-pass TerraSAR-X (TSX) imagery with high temporal (11 days) and spatial (10 m) resolution for monitoring bedfast ice as well as ice phenology of lakes in the zone of continuous permafrost in the Lena River Delta, Siberia. TSX backscatter intensity is shown to be an excellent tool for monitoring floating versus bedfast lake ice as well as ice phenology. TSX-derived timing of ice grounding and the ice growth model CLIMo are used to retrieve the ice thicknesses of the bedfast ice at points where in situ ice thickness measurements were available. Comparison shows good agreement in the year of field measurements. Additionally, for the first time, an 11-day sequential interferometric coherence time series is analyzed as a supplementary approach for the bedfast ice monitoring. The coherence time series detects most of the ice grounding as well as spring snow/ice melt onset. Overall, the results show the great value of TSX time series for monitoring Arctic lake ice and provide a basis for various applications: for instance, derivation of shallow lakes bathymetry, evaluation of winter water resources and locating fish winter habitat as well as estimation of taliks extent in permafrost.

Remote sensing, and unmanned aerial vehicles (UAVs) in particular, can be a valid tool for assessing the dynamics of cryotic features as frost blisters and to monitor the surface changes and the sublimation rates on perennially frozen lakes that host important ecosystems. In this paper, through the use of these remote sensing techniques, we aim to understand the type of groundwater supply of an Antarctic perennial frozen lake that encompasses two frost blisters (M1 and M2) through the temporal analysis of the features’ elevation changes (frost blisters and lake ice level). The frozen lake is located at Boulder Clay (northern Victoria Land, Antarctica). We relied on several photogrammetric models, past satellite images and ground pictures to conduct differencing of digital elevation models, areal variations and pixel counting. In addition, in situ measurements of the ice sublimation or snow accumulation were carried out. The two frost blisters showed different elevation trends with M1 higher in the past (1996–2004) than recently (2014–2019), while M2 showed an opposite trend, similarly to the ice level. Indeed, the linear regression between M2 elevation changes and the ice level variation was statistically significant, as well as with the annual thawing degree days, while M1 did not show significant results. From these results we can infer that the groundwater supply of M1 can be related to a sublake open talik (hydraulic system) as confirmed also by pressurized brines found below M1, during a drilling in summer 2019. For M2 the groundwater flow is still not completely clear although the hydrostatic system seems the easiest explanation as well as for the uplift of the lake ice.