Search

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

Satellite altimetry has been proven to measure water levels in lakes and rivers effectively. The Sentinel-3A satellite is equipped with a dual-frequency synthetic aperture radar altimeter (SRAL), which allows for inland water levels to be measured with higher precision and improved spatial resolution. However, in regions at middle and high latitudes, where many lakes are covered by ice during the winter, the non-uniformity of the altimeter footprint can substantially impact the accuracy of water level estimates, resulting in abnormal readings when applying standard SRAL synthetic aperture radar (SAR) waveform retracking algorithms (retrackers). In this study, a modified method is proposed to determine the current surface type of lakes, analyzing changes in backscattering coefficients and brightness temperature. This method aligns with ground station observations and ensures consistent surface type classification. Additionally, a dual-threshold algorithm that addresses the limitations of the original bimodal algorithm by identifying multiple peaks without needing elevation correction is introduced. This innovative approach significantly enhances the precision of equivalent water level measurements for ice-covered lakes. The study retrieves and compares the water level data of nine North American lakes covered by ice from 2016–2019 using the dual-threshold and the SAMOSA-3 algorithm with in situ data. For Lake Athabasca, Cedar Lake, Great Slave Lake, Lake Winnipeg, and Lake Erie, the root mean square error (RMSE) of SAMOSA-3 is 39.58 cm, 46.18 cm, 45.75 cm, 42.64 cm, and 6.89 cm, respectively. However, the dual-threshold algorithm achieves an RMSE of 6.75 cm, 9.47 cm, 5.90 cm, 7.67 cm, and 5.01 cm, respectively, representing a decrease of 75%, 79%, 87%, 82%, and 27%, respectively, compared to SAMOSA-3. The dual-threshold algorithm can accurately estimate water levels in ice-covered lakes during winter. It offers a promising prospect for achieving long-term, continuous, and high-precision water level measurements for middle- and high-latitude lakes.

Lake ice phenology (LIP), hiding information about lake energy and material exchange, serves as an important indicator of climate change. Utilizing an efficient technique to swiftly extract lake ice information is crucial in the field of lake ice research. The Bayesian ensemble change detection (BECD) algorithm stands out as a powerful tool, requiring no threshold compared to other algorithms and, instead, utilizing the probability of abrupt changes to detect positions. This method is predominantly employed by automatically extracting change points from time series data, showcasing its efficiency and accuracy, especially in revealing phenological and seasonal characteristics. This paper focuses on Bosten Lake (BL) and employs PMRS data in conjunction with the Bayesian change detection algorithm. It introduces an automated method for extracting LIP information based on the Bayesian change detection algorithm. In this study, the BECD algorithm was employed to extract lake ice phenology information from passive microwave remote sensing data on Bosten Lake. The reliability of the passive microwave remote sensing data was further investigated through cross-validation with MOD10A1 data. Additionally, the Mann–Kendall non-parametric test was applied to analyze the trends in lake ice phenology changes in Bosten Lake. Spatial variations were examined using MOD09GQ data. The results indicate: (1) The Bayesian change detection algorithm (BCDA), in conjunction with PMRS data, offers a high level of accuracy and reliability in extracting the lake ice freezing and thawing processes. It accurately captures the phenological parameters of BL’s ice. (2) The average start date of lake ice freezing is in mid-December, lasting for about three months, and the start date of ice thawing is usually in mid-March. The freezing duration (FD) of lake ice is relatively short, shortening each year, while the thawing speed is faster. The stability of the lake ice complete ice cover duration is poor, averaging 84 days. (3) The dynamic evolution of BL ice is rapid and regionally distinct, with the lake center, southwest, and southeast regions being the earliest areas for ice formation and thawing, while the northwest coastal and Huang Shui Gou areas experience later ice formation. (4) Since 1978, BL’s ice has exhibited noticeable trends: the onset of freezing, the commencement of thawing, complete thawing, and full freezing have progressively advanced in regard to dates. The periods of full ice coverage, ice presence, thawing, and freezing have all shown a tendency toward shorter durations. This study introduces an innovative method for LIP extraction, opening up new prospects for the study of lake ecosystem and strategy formulation, which is worthy of further exploration and application in other lakes and 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.

Landfast ice plays an important role in the nearshore hydrodynamics of large lakes, such as the dampening of surface waves and currents. In this study, previously developed landfast ice basal stress parameterizations were added to an unstructured grid hydrodynamic ice model to represent the effects of grounded ice keels and tensile strength of ice cover. Numerical experiments using this model were conducted to evaluate the development of coastal landfast ice in Lake Superior. A sensitivity study of the free parameters was conducted from December 2018 to May 2021 to cover both high and low ice cover winters in Lake Superior and was compared against observations from the United States National Ice Center. The model reproduces the annual variation in coastal landfast ice in Lake Superior, particularly in shallow nearshore areas and the semi-enclosed bays in the northern regions of the lake. Experiments also show that the growth of landfast ice is mainly controlled by the free parameter that controls the critical ice thickness for the activation of basal stress. Overall, the model tends to underestimate the extent of coastal landfast against observations.

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.

This article is an attempt to analyse and compare several selected parameters regarding ice phenomena using the correlation analysis of two lakes, which are benchmark lakes located in the central part of the Kashubian Lakeland. These lakes are: Raduńskie Górne, a larger one (387.2 hm2) and Ostrzyckie, a smaller one covering an area of 308.0 hm2. The analysis covered measurement sequences for the period 1971–2010. The material for analysis regarding data on ice phenomena on Lake Raduńskie Górne was obtained from the University of Gdańsk Limnological Station in Borucino while the data for Lake Ostrzyckie was obtained from the Institute of Meteorology and Water Management (IMGW). Relations between the following parameters of the ice regime were analysed: duration of the ice phenomena, duration of the ice cover, average and maximum thickness of the ice cover. The analysis and data comparisons have revealed that there are strong and very strong relations between the analysed parameters, which made it possible to develop proper statistical models. Hence, should there be a lack of data on the ice-cover pattern for one of these lakes it is possible to recreate them using the elaborated empirical models and data for the other.

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.

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.

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 models are a vital tool for studying the response of ice-covered lakes to changing climates throughout the world. The Canadian Lake Ice Model (CLIMo) is a one-dimensional freshwater ice cover model that simulates Arctic and sub-Arctic lake ice cover well. Modelling ice cover in temperate regions has presented challenges due to the differences in ice composition between northern and temperate region lake ice. This study presents a comparison of measured and modelled ice regimes, with a focus on refining CLIMo for temperate regions. The study sites include two temperate region lakes (MacDonald Lake and Clear Lake, Central Ontario) and two High Arctic lakes (Resolute Lake and Small Lake, Nunavut) where climate and ice cover information have been recorded over three seasons. The ice cover simulations were validated with a combination of time lapse imagery, field measurements of snow depth, snow density, ice thickness and albedo data, and historical ice records from the Canadian Ice Database (for Resolute Lake). Simulations of High Arctic lake ice cover show good agreement with previous studies for ice-on and ice-off dates (MAE 6 to 8 days). Unadjusted simulations for the temperate region lakes show good ice-on timing, but an under-representation of ice thickness, and earlier complete ice-off timing (~3 to 5 weeks). Field measurements were used to adjust the albedo values used in CLIMo, which resulted in improvements to both simulated ice thickness (~3 cm MAE compared to manual measurements), and ice-off timing, within 0 to 7 days (2 days MAE) of observations. These findings suggest regionally specific measurements of albedo can improve the accuracy of lake ice simulations, which further our knowledge of the response of temperate and High Arctic lake ice regimes to climate conditions.

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.

Lake ice is a dominant component of Canada’s landscape and can act as an indicator for how freshwater aquatic ecosystems are changing with warming climates. While lake ice monitoring through government networks has decreased in the last three decades, the increased availability of remote sensing images can help to provide consistent spatial and temporal coverage for areas with annual ice cover. Synthetic aperture radar (SAR) data are commonly used for lake ice monitoring, due to the acquisition of images in any condition (time of day or weather). Using Sentinel-1 A/B images, a high-density time series of SAR images was developed for Lake Hazen in Nunavut, Canada, from 2015−2018. These images were used to test two different methods of monitoring lake ice phenology: one method using the first difference between SAR images and another that applies the Otsu segmentation method. Ice phenology dates determined from the two methods were compared with visual interpretation of the Sentinel-1 images. Mean errors for the pixel comparison of the first difference method ranged 3−10 days for ice-on and ice-off, while average error values for the Otsu method ranged 2−10 days. Mean errors for comparisons of different sections of the lake ranged 0−15 days for the first difference method and 2−17 days for the Otsu method. This research demonstrates the value of temporally consistent image acquisition for improving the accuracy of lake ice monitoring.

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

Dual-polarized airborne passive microwave (PM) brightness temperatures (Tb) at 6.9 GHz H/V, 19 GHz H/V and 37 GHz H/V and spaceborne active microwave (AM) X-band (9.65 GHz VV, VH) backscatter (σ0) are observed coincident to in situ snow and lake-ice measurements collected over two lakes near Inuvik, Canada. Lake-ice thickness is found to be positively correlated with 19 GHz V emission (R = 0.67) and negatively with 19 GHz H emission (R = −0.79), indicating surface ice conditions influence microwave interaction. Lake ice types are delineated from TerraSAR-X synthetic aperture radar (SAR) images using the iterative region growing with semantics (IRGS) segmentation algorithm implemented in the MAGIC (MAp Guided Ice Classification) system. The spatial extent of derived ice type classes correspond well to in situ observations. The overall magnitude of emission at 19 GHz H and X-band VH σ0 increase with the scattering potential of associated ice types (grey/rafted ice). Transects of 6.9 GHz PM and 19 GHz PM exhibit positive relationships with VH σ0 over freshwater lake ice, with the greatest R coefficients at H-pol (R = 0.64, 0.46). Conversely, 6.9 GHz Tb and 19 GHz Tb exhibit negative R coefficients in regions of brackish water due to tubular bubble and brine inclusions in the ice. This study identifies congruency between PM and AM scattering mechanisms over lake ice for the purpose of identifying the influence of ice types on overall microwave interaction within the lake-ice system.

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.