Search

Lake ice phenology serves as a sensitive indicator of climate change in the lake-rich Northeast China. In this study, the freeze-up date (FUD), break-up date (BUD), and ice cover duration (ICD) of 31 lakes were extracted from a time series of the land water surface temperature (LWST) derived from the combined MOD11A1 and MYD11A1 products for the hydrological years 2001 to 2021. Our analysis showed a high correlation between the ice phenology measures derived by our study and those provided by hydrological records (R2 of 0.89) and public datasets (R2 > 0.7). There was a notable coherence in lake ice phenology in Northeast China, with a trend in later freeze-up (0.21 days/year) and earlier break-up (0.19 days/year) dates, resulting in shorter ice cover duration (0.50 days/year). The lake ice phenology of freshwater lakes exhibited a faster rate of change compared to saltwater lakes during the period from HY2001 to HY2020. We used redundancy analysis and correlation analysis to study the relationships between the LWST and lake ice phenology with various influencing factors, including lake properties, local climate factors, and atmospheric circulation. Solar radiation, latitude, and air temperature were found to be the primary factors. The FUD was more closely related to lake characteristics, while the BUD was linked to local climate factors. The large-scale oscillations were found to influence the changes in lake ice phenology via the coupled influence of air temperature and precipitation. The Antarctic Oscillation and North Atlantic Oscillation correlate more with LWST in winter, and the Arctic Oscillation correlates more with the ICD.

Lakes on the Qinghai Tibet Plateau (QTP) are widely distributed spatially, and they are mostly seasonally frozen. Due to global warming, the thickness and phenology of the lake ice has been changing, which profoundly affects the regional climate evolution. There are a few studies about lake ice in alpine regions, but the understanding of climatological characteristics of lake ice on the QTP is still limited. Based on a field experiment in the winter of 2022, the thermal conductivity of Qinghai Lake ice was determined as 1.64 W·m−1·°C−1. Airborne radar ice thickness data, meteorological observations, and remote sensing images were used to evaluate a quasi-steady ice model (Leppäranta model) performance of the lake. This is an analytic model of lake ice thickness and phenology. The long-term (1979–2017) ice history of the lake was simulated. The results showed that the modeled mean ice thickness was 0.35 m with a trend of −0.002 m·a−1, and the average freeze-up start (FUS) and break-up end (BUE) were 30 December and 5 April, respectively, which are close to the field and satellite observations. The simulated trend of the maximum ice thickness from 1979 to 2017 (0.004 m·a−1) was slightly higher than the observed result (0.003 m·a−1). The simulated trend was 0.20 d·a−1 for the FUS, −0.34 d·a−1 for the BUE, and −0.54 d·a−1 for the ice duration (ID). Correlation and detrending analysis were adopted for the contribution of meteorological factors. In the winters of 1979–2017, downward longwave radiation and air temperature were the two main factors that had the best correlation with lake ice thickness. In a detrending analysis, air temperature, downward longwave radiation, and solar radiation contributed the most to the average thickness variability, with contributions of 42%, 49%, and −48%, respectively, and to the maximum thickness variability, with contributions of 41%, 45%, and −48%, respectively. If the six meteorological factors (air temperature, downward longwave radiation, solar radiation, wind speed, pressure, and specific humidity) are detrending, ice thickness variability will increase 83% on average and 87% at maximum. Specific humidity, wind, and air pressure had a poor correlation with ice thickness. The findings in this study give insights into the long-term evolutionary trajectory of Qinghai Lake ice cover and serve as a point of reference for investigating other lakes in the QTP during cold seasons.

Lake ice phenology serves as a sensitive indicator of climate change in the lake-rich Northeast China. In this study, the freeze-up date (FUD), break-up date (BUD), and ice cover duration (ICD) of 31 lakes were extracted from a time series of the land water surface temperature (LWST) derived from the combined MOD11A1 and MYD11A1 products for the hydrological years 2001 to 2021. Our analysis showed a high correlation between the ice phenology measures derived by our study and those provided by hydrological records (R2 of 0.89) and public datasets (R2 > 0.7). There was a notable coherence in lake ice phenology in Northeast China, with a trend in later freeze-up (0.21 days/year) and earlier break-up (0.19 days/year) dates, resulting in shorter ice cover duration (0.50 days/year). The lake ice phenology of freshwater lakes exhibited a faster rate of change compared to saltwater lakes during the period from HY2001 to HY2020. We used redundancy analysis and correlation analysis to study the relationships between the LWST and lake ice phenology with various influencing factors, including lake properties, local climate factors, and atmospheric circulation. Solar radiation, latitude, and air temperature were found to be the primary factors. The FUD was more closely related to lake characteristics, while the BUD was linked to local climate factors. The large-scale oscillations were found to influence the changes in lake ice phenology via the coupled influence of air temperature and precipitation. The Antarctic Oscillation and North Atlantic Oscillation correlate more with LWST in winter, and the Arctic Oscillation correlates more with the ICD. [ABSTRACT FROM AUTHOR]

Lake ice thickness (LIT) is one of the key climate variables in the lake ice domain, but there are currently large uncertainties in the retrieval of LIT. We present and validate a new LIT retrieval method that utilizes ICESat-2 data to assist CryoSat-2 echo peak selection, aiming to improve the accuracy of LIT retrieval and enable data acquisition without on-site measurements. The method involves screening out similar ICESat-2 and CryoSat-2 tracks based on time and space constraints. It also involves dynamically adjusting the range constraint window of CryoSat-2 waveforms based on the high-precision lake ice surface ellipsoid height obtained from ICESat-2/ATL06 data. Within this range constraint window, the peak selection strategy is used to determine the scattering interfaces between snow-ice and ice-water. By utilizing the distance between the scattering horizons, the thickness of the lake ice can be determined. We performed the ice thickness retrieval experiment for Baker Lake in winter and verified it against the on-site measurement data. The results showed that the accuracy was about 0.143 m. At the same time, we performed the ice thickness retrieval experiment for Great Bear Lake (GBL), which does not have on-site measurement data, and compared it with the climate change trend of GBL. The results showed that the retrieval results were consistent with the climate change trend of GBL, confirming the validity of the proposed method.

Lake ice thickness (LIT) is a sensitive indicator of climate change, identified as a thematic variable of Lakes as an Essential Climate Variable (ECV) by the Global Climate Observing System (GCOS). Here, we present a novel and efficient analytically based retracking approach for estimating LIT from high-resolution Ku-band (13.6 GHz) synthetic-aperture radar (SAR) altimetry data. The retracker method is based on the analytical modeling of the SAR radar echoes over ice-covered lakes that show a characteristic double-peak feature attributed to the reflection of the Ku-band radar waves at the snow–ice and ice–water interfaces. The method is applied to Sentinel-6 Unfocused SAR (UFSAR) and Fully Focused SAR (FFSAR) data, with their corresponding tailored waveform model, referred to as the SAR_LIT and FFSAR_LIT retracker, respectively. We found that LIT retrievals from Sentinel-6 high-resolution SAR data at different posting rates are fully consistent with the LIT estimations obtained from thermodynamic lake ice model simulations and from low-resolution mode (LRM) Sentinel-6 and Jason-3 data over two ice seasons during the tandem phase of the two satellites, demonstrating the continuity between LRM and SAR LIT retrievals. By comparing the Sentinel-6 SAR LIT estimates to optical/radar images, we found that the Sentinel-6 LIT measurements are fully consistent with the evolution of the lake surface conditions, accurately capturing the seasonal transitions of ice formation and melt. The uncertainty in the LIT estimates obtained with Sentinel-6 UFSAR data at 20 Hz is in the order of 5 cm, meeting the GCOS requirements for LIT measurements. This uncertainty is significantly smaller, by a factor of 2 to 3 times, than the uncertainty obtained with LRM data. The FFSAR processing at 140 Hz provides even better LIT estimates, with 20% smaller uncertainties. The LIT retracker analysis performed on data at the higher posting rate (140 Hz) shows increased performance in comparison to the 20 Hz data, especially during the melt transition period, due to the increased statistics. The LIT analysis has been performed over two representative lakes, Great Slave Lake and Baker Lake (Canada), demonstrating that the results are robust and hold for lake targets that differ in terms of size, bathymetry, snow/ice properties, and seasonal evolution of LIT. The SAR LIT retrackers presented are promising tools for monitoring the inter-annual variability and trends in LIT from current and future SAR altimetry missions.

Lake ice thickness (LIT) is a sensitive indicator of climate change, identified as a thematic variable of Lakes as an Essential Climate Variable (ECV) by the Global Climate Observing System (GCOS). Here, we present a novel and efficient analytically based retracking approach for estimating LIT from high-resolution Ku-band (13.6 GHz) synthetic-aperture radar (SAR) altimetry data. The retracker method is based on the analytical modeling of the SAR radar echoes over ice-covered lakes that show a characteristic double-peak feature attributed to the reflection of the Ku-band radar waves at the snow–ice and ice–water interfaces. The method is applied to Sentinel-6 Unfocused SAR (UFSAR) and Fully Focused SAR (FFSAR) data, with their corresponding tailored waveform model, referred to as the SAR_LIT and FFSAR_LIT retracker, respectively. We found that LIT retrievals from Sentinel-6 high-resolution SAR data at different posting rates are fully consistent with the LIT estimations obtained from thermodynamic lake ice model simulations and from low-resolution mode (LRM) Sentinel-6 and Jason-3 data over two ice seasons during the tandem phase of the two satellites, demonstrating the continuity between LRM and SAR LIT retrievals. By comparing the Sentinel-6 SAR LIT estimates to optical/radar images, we found that the Sentinel-6 LIT measurements are fully consistent with the evolution of the lake surface conditions, accurately capturing the seasonal transitions of ice formation and melt. The uncertainty in the LIT estimates obtained with Sentinel-6 UFSAR data at 20 Hz is in the order of 5 cm, meeting the GCOS requirements for LIT measurements. This uncertainty is significantly smaller, by a factor of 2 to 3 times, than the uncertainty obtained with LRM data. The FFSAR processing at 140 Hz provides even better LIT estimates, with 20% smaller uncertainties. The LIT retracker analysis performed on data at the higher posting rate (140 Hz) shows increased performance in comparison to the 20 Hz data, especially during the melt transition period, due to the increased statistics. The LIT analysis has been performed over two representative lakes, Great Slave Lake and Baker Lake (Canada), demonstrating that the results are robust and hold for lake targets that differ in terms of size, bathymetry, snow/ice properties, and seasonal evolution of LIT. The SAR LIT retrackers presented are promising tools for monitoring the inter-annual variability and trends in LIT from current and future SAR altimetry missions. [ABSTRACT FROM AUTHOR]

Lakes on the Qinghai Tibet Plateau (QTP) are widely distributed spatially, and they are mostly seasonally frozen. Due to global warming, the thickness and phenology of the lake ice has been changing, which profoundly affects the regional climate evolution. There are a few studies about lake ice in alpine regions, but the understanding of climatological characteristics of lake ice on the QTP is still limited. Based on a field experiment in the winter of 2022, the thermal conductivity of Qinghai Lake ice was determined as 1.64 W·m−1·°C−1. Airborne radar ice thickness data, meteorological observations, and remote sensing images were used to evaluate a quasi-steady ice model (Leppäranta model) performance of the lake. This is an analytic model of lake ice thickness and phenology. The long-term (1979–2017) ice history of the lake was simulated. The results showed that the modeled mean ice thickness was 0.35 m with a trend of −0.002 m·a−1, and the average freeze-up start (FUS) and break-up end (BUE) were 30 December and 5 April, respectively, which are close to the field and satellite observations. The simulated trend of the maximum ice thickness from 1979 to 2017 (0.004 m·a−1) was slightly higher than the observed result (0.003 m·a−1). The simulated trend was 0.20 d·a−1 for the FUS, −0.34 d·a−1 for the BUE, and −0.54 d·a−1 for the ice duration (ID). Correlation and detrending analysis were adopted for the contribution of meteorological factors. In the winters of 1979–2017, downward longwave radiation and air temperature were the two main factors that had the best correlation with lake ice thickness. In a detrending analysis, air temperature, downward longwave radiation, and solar radiation contributed the most to the average thickness variability, with contributions of 42%, 49%, and −48%, respectively, and to the maximum thickness variability, with contributions of 41%, 45%, and −48%, respectively. If the six meteorological factors (air temperature, downward longwave radiation, solar radiation, wind speed, pressure, and specific humidity) are detrending, ice thickness variability will increase 83% on average and 87% at maximum. Specific humidity, wind, and air pressure had a poor correlation with ice thickness. The findings in this study give insights into the long-term evolutionary trajectory of Qinghai Lake ice cover and serve as a point of reference for investigating other lakes in the QTP during cold seasons. [ABSTRACT FROM AUTHOR]

In this work, backscattering signatures of snow-covered lake ice and sea ice from X- and Ku-band synthetic aperture radar (SAR) data are investigated. The SAR data were acquired with the ESA airborne SnowSAR sensor in winter 2012 over Lake Orajärvi in northern Finland and over landfast ice in the Bay of Bothnia of the Baltic Sea. Co-incident with the SnowSAR acquisitions, in situ snow and ice data were measured. In addition, time series of TerraSAR-X images and ice mass balance buoy data were acquired for Lake Orajärvi in 2011–2012. The main objective of our study was to investigate relationships between SAR backscattering signatures and snow depth over lake and sea ice, with the ultimate objective of assessing the feasibility of retrieval of snow characteristics using X- and Ku-band dual-polarization (VV and VH) SAR over freshwater or sea ice. This study constitutes the first comprehensive survey of snow backscattering signatures at these two combined frequencies over both lake and sea ice. For lake ice, we show that X-band VH-polarized backscattering coefficient (σo) and the Ku-band VV/VH-ratio exhibited the highest sensitivity to the snow depth. For sea ice, the highest sensitivity to the snow depth was found from the Ku-band VV-polarized σo and the Ku-band VV/VH-ratio. However, the observed relations were relatively weak, indicating that at least for the prevailing snow conditions, obtaining reliable estimates of snow depth over lake and sea ice would be challenging using only X- and Ku-band backscattering information.

Lake ice phenology is an indicator of past and present climate, it is sensitive to regional and global climate change. In the past few decades, the climate of Central Asia has changed significantly due to global warming and anthropogenic activities. However, there are few studies on the lake ice phenology in Central Asia. In this study, the lake ice phenology of 53 lakes in Central Asia were extracted using MODIS daily LST products from 2002 to 2020. The results show that MODIS-extracted lake ice phenology is generally consistent with Landsat-extracted and AVHRR-extracted lake ice phenology. Generally, lakes in Central Asia start to freeze from October to December. The trends in the lake ice phenology show strong regional differences. Lakes distributed along the Kunlun Mountains show overall delayed trends in all lake ice phenology variables, while lakes located in southwestern Central Asia show clear advancing trends in the freeze-up start dates (7.06 days) and breakup end dates (6.81 days). Correlations between the phenology of lake ice and local and climatic factors suggest that the ice breakup process and the duration of its complete coverage depend more on heat, while precipitation mainly affects the freezing time of the ice. Wind speed mainly affects the time of completely frozen of ice. In general, the breakup process is more susceptible to climatic factors, while local factors have strong influences on the freeze-up process.

Lake ice thickness (LIT) is one of the key climate variables in the lake ice domain, but there are currently large uncertainties in the retrieval of LIT. We present and validate a new LIT retrieval method that utilizes ICESat-2 data to assist CryoSat-2 echo peak selection, aiming to improve the accuracy of LIT retrieval and enable data acquisition without on-site measurements. The method involves screening out similar ICESat-2 and CryoSat-2 tracks based on time and space constraints. It also involves dynamically adjusting the range constraint window of CryoSat-2 waveforms based on the high-precision lake ice surface ellipsoid height obtained from ICESat-2/ATL06 data. Within this range constraint window, the peak selection strategy is used to determine the scattering interfaces between snow-ice and ice-water. By utilizing the distance between the scattering horizons, the thickness of the lake ice can be determined. We performed the ice thickness retrieval experiment for Baker Lake in winter and verified it against the on-site measurement data. The results showed that the accuracy was about 0.143 m. At the same time, we performed the ice thickness retrieval experiment for Great Bear Lake (GBL), which does not have on-site measurement data, and compared it with the climate change trend of GBL. The results showed that the retrieval results were consistent with the climate change trend of GBL, confirming the validity of the proposed method. [ABSTRACT FROM AUTHOR]

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.

In this work, backscattering signatures of snow-covered lake ice and sea ice from X- and Ku-band synthetic aperture radar (SAR) data are investigated. The SAR data were acquired with the ESA airborne SnowSAR sensor in winter 2012 over Lake Orajärvi in northern Finland and over landfast ice in the Bay of Bothnia of the Baltic Sea. Co-incident with the SnowSAR acquisitions, in situ snow and ice data were measured. In addition, time series of TerraSAR-X images and ice mass balance buoy data were acquired for Lake Orajärvi in 2011–2012. The main objective of our study was to investigate relationships between SAR backscattering signatures and snow depth over lake and sea ice, with the ultimate objective of assessing the feasibility of retrieval of snow characteristics using X- and Ku-band dual-polarization (VV and VH) SAR over freshwater or sea ice. This study constitutes the first comprehensive survey of snow backscattering signatures at these two combined frequencies over both lake and sea ice. For lake ice, we show that X-band VH-polarized backscattering coefficient ( σ o ) and the Ku-band VV/VH-ratio exhibited the highest sensitivity to the snow depth. For sea ice, the highest sensitivity to the snow depth was found from the Ku-band VV-polarized σ o and the Ku-band VV/VH-ratio. However, the observed relations were relatively weak, indicating that at least for the prevailing snow conditions, obtaining reliable estimates of snow depth over lake and sea ice would be challenging using only X- and Ku-band backscattering information. [ABSTRACT FROM AUTHOR]

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.

A thorough understanding of the freshwater ice process received considerable critical attention due to increasing winter recreations and ice engineering. The development of the lake ice process of Chagan Lake was monitored using MODIS and Landsat images over eight consecutive snow seasons from October 2013 to April 2021. We derived the lake ice phenology from an eight-day time series of lake water skin temperature (LWST) provided by MODIS, including freeze-up date, break-up date, and ice cover duration. We discovered a large-scale fracture extending from northwest to southeast that repeatedly appeared on Landsat images since 1986. A novel fractal-based auto-extraction is proposed to extract the length and angle of these fractures. We also carried out a field campaign and an ice ridge was found at the southernmost part of what we observed from the images. Moreover, we explained the fracturing development by thermal changes, wind in lake, and underlying flow. Results show that the lake ice fracture is nearly perpendicular to the dominant wind direction in the cold season, which indicates the crucial role of wind on lake ice fracturing.

Lake ice phenology is a climate-sensitive indicator. However, ground-based monitoring suffers from the limitations of human vision and the difficulty of its implementation in harsh environments. Remote sensing provides great potential to detect lake ice phenology. In this study, a new automated method was developed to extract lake ice phenology parameters by capturing the temporal pattern of the transitional water/ice phase using a parameterized time function. The method is based on Moderate-Resolution Imaging Spectroradiometer (MODIS) daily temperature products, which have unique potential for monitoring lake ice cover as a result of providing four observations per day at 1 km spatial resolution from 2002 to 2016. Three seasonally ice-covered lakes with different characteristics in different climate regions were selected to test the method during the period of 2002–2016. The temporal pattern of water/ice transition phase was determined on the basis of unfrozen water cover fraction extracted from the MODIS daily temperature data, and was compared with the MODIS snow and reflectance products and Landsat images. A good agreement with an R2 of above 0.8 was found when compared with the MODIS snow product. The annual variation of extracted ice phenology dates showed good consistency with the MODIS reflectance and AMSR-E/2 products. The approach was then applied to nine seasonally ice-covered lakes in northern China from 2002 to 2016. The strongest tendency towards a later freeze-up start date was revealed in Lake Qinghai (6.31 days/10 yr) among the lakes in Tibetan plateau, and the break-up start and end dates rapidly shifted towards earlier dates in Lake Hulun (−3.73 days/10 yr; −5.02 days/10 yr). The method is suitable for estimating and monitoring ice phenology on different types of lakes over large scales and has a strong potential to provide valuable information on the responses of ice processes to climate change.

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° 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°, but also reflect the U-shaped distribution of the diurnal albedo as the solar altitude angle exceeds 15°. 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.

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. [ABSTRACT FROM AUTHOR]

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 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.

Remote sensing of ice phenology for small lakes is hindered by a lack of satellite observations with both high temporal and spatial resolutions. By merging multi-source satellite data over individual lakes, we present a new algorithm that successfully estimates ice freeze and thaw timing for lakes with surface areas as small as 0.13 km2 and obtains consistent results across a range of lake sizes. We have developed an approach for classifying ice pixels based on the red reflectance band of Moderate Resolution Imaging Spectroradiometer (MODIS) imagery, with a threshold calibrated against ice fraction from Landsat Fmask over each lake. Using a filter derived from the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) surface air temperature product, we removed outliers in the time series of lake ice fraction. The time series of lake ice fraction was then applied to identify lake ice breakup and freezeup dates. Validation results from over 296 lakes in Maine indicate that the satellite-based lake ice timing detection algorithm perform well, with mean absolute error (MAE) of 5.54 days for breakup dates and 7.31 days for freezeup dates. This algorithm can be applied to lakes worldwide, including the nearly two million lakes with surface area between 0.1 and 1 km2.

Thermokarst lakes in the Arctic and Subarctic release carbon from thawing permafrost in the form of methane and carbon dioxide with important implications for regional and global carbon cycles. Lake ice impedes the release of gas during the winter. For instance, bubbles released from lake sediments become trapped in downward growing lake ice, resulting in vertically-oriented bubble columns in the ice that are visible on the lake surface. We here describe a classification technique using an object-based image analysis (OBIA) framework to successfully map ebullition bubbles in airborne imagery of early winter ice on an interior Alaska thermokarst lake. Ebullition bubbles appear as white patches in high-resolution optical remote sensing images of snow-free lake ice acquired in early winter and, thus, can be mapped across whole lake areas. We used high-resolution (9–11 cm) aerial images acquired two and four days following freeze-up in the years 2011 and 2012, respectively. The design of multiresolution segmentation and region-specific classification rulesets allowed the identification of bubble features and separation from other confounding factors such as snow, submerged and floating vegetation, shadows, and open water. The OBIA technique had an accuracy of >95% for mapping ebullition bubble patches in early winter lake ice. Overall, we mapped 1195 and 1860 ebullition bubble patches in the 2011 and 2012 images, respectively. The percent surface area of lake ice covered with ebullition bubble patches for 2011 was 2.14% and for 2012 was 2.67%, representing a conservative whole lake estimate of bubble patches compared to ground surveys usually conducted on thicker ice 10 or more days after freeze-up. Our findings suggest that the information derived from high-resolution optical images of lake ice can supplement spatially limited field sampling methods to better estimate methane flux from individual lakes. The method can also be used to improve estimates of methane ebullition from numerous lakes within larger regions.

Lake ice phenology is an indicator of past and present climate, it is sensitive to regional and global climate change. In the past few decades, the climate of Central Asia has changed significantly due to global warming and anthropogenic activities. However, there are few studies on the lake ice phenology in Central Asia. In this study, the lake ice phenology of 53 lakes in Central Asia were extracted using MODIS daily LST products from 2002 to 2020. The results show that MODIS-extracted lake ice phenology is generally consistent with Landsat-extracted and AVHRR-extracted lake ice phenology. Generally, lakes in Central Asia start to freeze from October to December. The trends in the lake ice phenology show strong regional differences. Lakes distributed along the Kunlun Mountains show overall delayed trends in all lake ice phenology variables, while lakes located in southwestern Central Asia show clear advancing trends in the freeze-up start dates (7.06 days) and breakup end dates (6.81 days). Correlations between the phenology of lake ice and local and climatic factors suggest that the ice breakup process and the duration of its complete coverage depend more on heat, while precipitation mainly affects the freezing time of the ice. Wind speed mainly affects the time of completely frozen of ice. In general, the breakup process is more susceptible to climatic factors, while local factors have strong influences on the freeze-up process. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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.

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.

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.

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.

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.

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.

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.

Over the past decades, the cryosphere has changed significantly in High Mountain Asia (HMA), leading to multiple natural hazards such as rock–ice avalanches, glacier collapse, debris flows, landslides, and glacial lake outburst floods (GLOFs). Monitoring cryosphere change and evaluating its hydrological effects are essential for studying climate change, the hydrological cycle, water resource management, and natural disaster mitigation and prevention. However, knowledge gaps, data uncertainties, and other substantial challenges limit comprehensive research in climate–cryosphere–hydrology–hazard systems. To address this, we provide an up-to-date, comprehensive, multidisciplinary review of remote sensing techniques in cryosphere studies, demonstrating primary methodologies for delineating glaciers and measuring geodetic glacier mass balance change, glacier thickness, glacier motion or ice velocity, snow extent and water equivalent, frozen ground or frozen soil, lake ice, and glacier-related hazards. The principal results and data achievements are summarized, including URL links for available products and related data platforms. We then describe the main challenges for cryosphere monitoring using satellite-based datasets. Among these challenges, the most significant limitations in accurate data inversion from remotely sensed data are attributed to the high uncertainties and inconsistent estimations due to rough terrain, the various techniques employed, data variability across the same regions (e.g., glacier mass balance change, snow depth retrieval, and the active layer thickness of frozen ground), and poor-quality optical images due to cloudy weather. The paucity of ground observations and validations with few long-term, continuous datasets also limits the utilization of satellite-based cryosphere studies and large-scale hydrological models. Lastly, we address potential breakthroughs in future studies, i.e., (1) outlining debris-covered glacier margins explicitly involving glacier areas in rough mountain shadows, (2) developing highly accurate snow depth retrieval methods by establishing a microwave emission model of snowpack in mountainous regions, (3) advancing techniques for subsurface complex freeze–thaw process observations from space, (4) filling knowledge gaps on scattering mechanisms varying with surface features (e.g., lake ice thickness and varying snow features on lake ice), and (5) improving and cross-verifying the data retrieval accuracy by combining different remote sensing techniques and physical models using machine learning methods and assimilation of multiple high-temporal-resolution datasets from multiple platforms. This comprehensive, multidisciplinary review highlights cryospheric studies incorporating spaceborne observations and hydrological models from diversified techniques/methodologies (e.g., multi-spectral optical data with thermal bands, SAR, InSAR, passive microwave, and altimetry), providing a valuable reference for what scientists have achieved in cryosphere change research and its hydrological effects on the Third Pole.