Arctic Sea Ice Thickness Estimation from CryoSat-2 Satellite Data Using Machine Learning-Based Lead Detection

Satellite altimeters have been used to monitor Arctic sea ice thickness since the early 2000s. In order to estimate sea ice thickness from satellite altimeter data, leads (i.e., cracks between ice floes) should first be identified for the calculation of sea ice freeboard. In this study, we proposed novel approaches for lead detection using two machine learning algorithms: decision trees and random forest. CryoSat-2 satellite data collected in March and April of 2011-2014 over the Arctic region were used to extract waveform parameters that show the characteristics of leads, ice floes and ocean, including stack standard deviation, stack skewness, stack kurtosis, pulse peakiness and backscatter sigma-0. The parameters were used to identify leads in the machine learning models. Results show that the proposed approaches, with overall accuracy >90%, produced much better performance than existing lead detection methods based on simple thresholding approaches. Sea ice thickness estimated based on the machine learning-detected leads was compared to the averaged Airborne Electromagnetic (AEM)-bird data collected over two days during the CryoSat Validation experiment (CryoVex) field campaign in April 2011. This comparison showed that the proposed machine learning methods had better performance (up to r = 0.83 and Root Mean Square Error (RMSE) = 0.29 m) compared to thickness estimation based on existing lead detection methods (RMSE = 0.86-0.93 m). Sea ice thickness based on the machine learning approaches showed a consistent decline from 2011-2013 and rebounded in 2014.open0

A 10-year record of Arctic summer sea ice freeboard from CryoSat-2

Satellite observations of pan-Arctic sea ice thickness have so far been constrained to winter months. For radar altimeters, conventional methods cannot differentiate leads from meltwater ponds that accumulate at the ice surface in summer months, which is a critical step in the ice thickness calculation. Here, we use over 350 optical and synthetic aperture radar (SAR) images from the summer months to train a 1D convolution neural network for separating CryoSat-2 radar altimeter returns from sea ice floes and leads with an accuracy >80%. This enables us to generate the first pan-Arctic measurements of sea ice radar freeboard for May–September between 2011 and 2020. Results indicate that the freeboard distributions in May and September compare closely to those from a conventional ‘winter’ processor in April and October, respectively. The freeboards capture expected patterns of sea ice melt over the Arctic summer, matching well to ice draft observations from the Beaufort Gyre Exploration Program (BGEP) moorings. However, compared to airborne laser scanner freeboards from Operation IceBridge and airborne EM ice thickness surveys from the Alfred Wegener Institute (AWI) IceBird program, CryoSat-2 freeboards are underestimated by 0.02–0.2 m, and ice thickness is underestimated by 0.28–1.0 m, with the largest differences being over thicker multi-year sea ice. To create the first pan-Arctic summer sea ice thickness dataset we must address primary sources of uncertainty in the conversion from radar freeboard to ice thickness

Arctic lead detection using a waveform mixture algorithm from CryoSat-2 data

We propose a waveform mixture algorithm to detect leads from CryoSat-2 data, which is novel and different from the existing threshold-based lead detection methods. The waveform mixture algorithm adopts the concept of spectral mixture analysis, which is widely used in the field of hyperspectral image analysis. This lead detection method was evaluated with high-resolution (250 m) MODIS images and showed comparable and promising performance in detecting leads when compared to the previous methods. The robustness of the proposed approach also lies in the fact that it does not require the rescaling of parameters (i.e., stack standard deviation, stack skewness, stack kurtosis, pulse peakiness, and backscatter sigma(0)), as it directly uses L1B waveform data, unlike the existing threshold-based methods. Monthly lead fraction maps were produced by the waveform mixture algorithm, which shows interannual variability of recent sea ice cover during 2011-2016, excluding the summer season (i.e., June to September). We also compared the lead fraction maps to other lead fraction maps generated from previously published data sets, resulting in similar spatiotemporal patterns

Monitoring and Characterization of Arctic Sea Ice using Radar Altimetry

Department of Urban and Environmental Engineering (Environmental Science and Engineering)Launching CryoSat-2, which is a current radar altimeter mission for the monitoring of polar region enables to produce monthly based sea ice thickness since April 2010. The Sea ice thickness cannot be measured directly by satellite. Sea ice freeboard that is an elevation above sea level can be converted in to sea ice thickness by assuming hydrostatic equilibrium. Sea ice leads (e.g., linear cracks in sea ices) are regarded as sea surface tie points for the estimation of sea ice freeboard. Identifying the sea ice leads is one of the core factors to retrieve sea ice thickness. The surface elevation is estimated by the use of Threshold First maxima Retracker Algorithm (TFMRA) for a 40% threshold using CryoSat-2 L1b data and the leads are detected by machine learning approaches such as decision trees and random forest. The machine learning produces better accuracy for the sea ice thickness than previous simple thresholding approach, validating EM-31, airborne sea ice thickness observations. A novel method to overcome previous threshold based lead detection methods for identifying leads is developed, which is waveform mixture algorithm that linear mixture analysis is applied in terms of waveforms. The waveform mixture algorithm can distinguish leads without beam behavior parameters and backscatter sigma-0 but just use waveforms, which is less affected by updating baseline for CryoSat-2. In addition to the development of the algorithms, a scientific research is carried out. Causes for sea ice anomaly phenomenon in November 2016 is investigated. Eventually, sea ice the volume derived by thickness is used for the analysis of sea ice extent minimum in November 2016 and suggest a new insight of sea ice minimum phenomenon. Unlike sea ice extent, the sea ice volume is not a minimum in November 2016. However, since the base period for sea ice volume is short, it is hard to mention climatology of sea ice volume.ope

Detection and classification of sea ice from spaceborne multi-frequency synthetic aperture radar imagery and radar altimetry

The sea ice cover in the Arctic is undergoing drastic changes. Since the start of satellite observations by microwave remote sensing in the late 1970\u27s, the maximum summer sea ice extent has been decreasing and thereby causing a generally thinner and younger sea ice cover. Spaceborne radar remote sensing facilitates the determination of sea ice properties in a changing climate with the high spatio-temporal resolution necessary for a better understanding of the ongoing processes as well as safe navigation and operation in ice infested waters.The work presented in this thesis focuses on the one hand on synergies of multi-frequency spaceborne synthetic aperture radar (SAR) imagery for sea ice classification. On the other hand, the fusion of radar altimetry observations with near-coincidental SAR imagery is investigated for its potential to improve 3-dimensional sea ice information retrieval.Investigations of ice/water classification of C- and L-band SAR imagery with a feed-forward neural network demonstrated the capabilities of both frequencies to outline the sea ice edge with good accuracy. Classification results also indicate that a combination of both frequencies can improve the identification of thin ice areas within the ice pack compared to C-band alone. Incidence angle normalisation has proven to increase class separability of different ice types. Analysis of incidence angle dependence between 19-47\ub0 at co- and cross-polarisation from Sentinel-1 C-band images closed a gap in existing slope estimates at cross-polarisation for multiyear sea ice and confirms values obtained in other regions of the Arctic or with different sensors. Furthermore, it demonstrated that insufficient noise correction of the first subswath at cross-polarisation increased the slope estimates by 0.01 dB/1\ub0 for multiyear ice. The incidence angle dependence of the Sentinel-1 noise floor affected smoother first-year sea ice and made the first subswath unusable for reliable incidence angle estimates in those cases.Radar altimetry can complete the 2-dimensional sea ice picture with thickness information. By comparison of SAR imagery with altimeter waveforms from CryoSat-2, it is demonstrated that waveforms respond well to changes of the sea ice surface in the order of a few hundred metres to a few kilometres. Freeboard estimates do however not always correspond to these changes especially when mixtures of different ice types are found within the footprint. Homogeneous ice floes of about 10 km are necessary for robust averaged freeboard estimates. The results demonstrate that multi-frequency and multi-sensor approaches open up for future improvements of sea ice retrievals from radar remote sensing techniques, but access to in-situ data for training and validation will be critical

A long-term record of sea ice thickness in the Canadian Arctic

Sea ice plays a vital role in the Arctic region and affects numerous processes: it influences the radiative budget by reflecting sunlight and acts as a barrier for heat transport between atmosphere and ocean; it influences Arctic ecosystems as a habitat for different species; it is important for hunting and travel for local communities; and it acts as a hazard for marine shipping. Monitoring sea ice, specifically its thickness, is essential in understanding how it is changing with ongoing global warming.This thesis presents a novel method to create a long-term record (1996-2020) for sea ice thickness in the Canadian Arctic and assesses how sea ice thickness changed and what the impacts of these changes are.This thesis initially aimed to extract a long-term sea ice thickness record for the Canadian Arctic from satellite altimetry. However, it revealed that assumptions regarding the snowpack, sea ice density, and processing algorithms highly influence conclusions on sea ice thickness state and trends, and this approach was rejected. Instead, this thesis presents a proxy sea ice thickness product for the Canadian Arctic using ice charts, which for the first time consistently covers the Canadian Arctic Archipelago. In the final research chapter, this sea ice thickness proxy product and ice charts are used to assess sea ice changes in the Canadian Arctic Archipelago and their impact on accessibility.Sea ice has thinned across most of the Canadian Arctic region, with a mean change over the full area of 38.5 cm for November and 20.5 cm for April over the period 1996-2020. Moreover, the marine navigability is shown to increase in the access channels to the Canadian Arctic Archipelago, which enhances the possibilities for resupply for local communities. However, with continuing dynamic influx of old and thick sea ice, there is no change in full navigability of the Northwest Passage connecting the Atlantic and Pacific Oceans

Retrieving Sea Level and Freeboard in the Arctic: A Review of Current Radar Altimetry Methodologies and Future Perspectives

Spaceborne radar altimeters record echo waveforms over all Earth surfaces, but their interpretation and quantitative exploitation over the Arctic Ocean is particularly challenging. Radar returns may be from all ocean, all sea ice, or a mixture of the two, so the first task is the determination of which surface and then an interpretation of the signal to give range. Subsequently, corrections have to be applied for various surface and atmospheric effects before making a comparison with a reference level. This paper discusses the drivers for improved altimetry in the Arctic and then reviews the various approaches that have been used to achieve the initial classification and subsequent retracking over these diverse surfaces, showing examples from both LRM (low resolution mode) and SAR (synthetic aperture radar) altimeters. The review then discusses the issues concerning corrections, including the choices between using other remote-sensing measurements and using those from models or climatology. The paper finishes with some perspectives on future developments, incorporating secondary frequency, interferometric SAR and opportunities for fusion with measurements from laser altimetry or from the SMOS salinity sensor, and provides a full list of relevant abbreviations

Evaluation of Machine Learning Algorithms for the Classification of Lake Ice and Open Water from Sentinel-3 SAR Altimetry Waveforms

Lakes cover a significant fraction of the landscape in many northern countries. They play a key role in regulating weather and climate and also have a significant impact on northern communities since the presence (or absence), extent and thickness of lake ice affect transportation (ice roads), food availability, recreational activities, and tourism in wintertime. The drastic decline in in-situ observations of lake ice phenology (i.e., freeze-up and break-up dates and ice cover duration) and lake ice thickness globally over the last three decades make remote sensing technology a viable means for monitoring lake ice conditions. Although satellite radar altimetry has been used in various cryospheric and hydrological studies, little work has been conducted on lake ice compared to, for example, sea ice and the estimation of lake water levels. This study was carried out using Sentinel-3A/B SAR altimetry data acquired over three ice seasons (2018-2019, 2019-2020 and 2020-2021) at 11 large lakes across the Northern Hemisphere. We explored the information provided by radar waveforms to discriminate between open water, first (young) ice, growing ice and melting ice using machine learning models. To characterize the waveforms, seven waveform parameters were derived: Leading Edge Width (LEW), Offset Center of Gravity (OCOG) Width, Pulse Peakiness (PP), backscatter coefficient (Sigma0), late tail to peak power (LTTP), early tail to peak power (ETTP) and the maximum value of the echo power. Four machine learning algorithms including Random Forest (RF), Gradient Boosting Trees (GBT), K Nearest Neighbour (KNN) and Support Vector Machine (SVM) classifiers were tested to assess their capability in classifying the lake surfaces across all years. Manual class labelling based on Sentinel-3 Synthetic Aperture Radar Altimeter (SRAL) waveforms and complementary satellite data (Sentinel-1 imaging SAR data, Sentinel-2 Multispectral Instrument (MSI) Level 1C data, and MODIS Aqua/Terra data) was performed to create training and test samples for the classifiers. Accuracies greater than 95% were achieved across all classifiers using a 4-parameter combination (Sigma0, PP, OCOG Width, and LEW). Amongst all waveform parameters, Sigma0, OCOG width and PP were found to be the most important parameters for discriminating between lake ice and open water. Despite showing comparable classification performances in the overall classification, RF and KNN are found to be a better fit for global lake ice mapping as both are less sensitive to their internal hyperparameters and have faster processing speeds. Additionally, consistent results (>93.7% accuracy in all classifiers) achieved on the accuracy assessment carried out for each lake revealed the strength of the classifiers for spatial transferability. Implementation of RF and KNN could be valuable in a pre-or post-processing step for identifying lake surface conditions under which the retrieval of water level and ice thickness may be limited or not possible and, therefore, inform algorithms currently used for the generation of operational or research products. While the research focused on 11 of the largest lakes of the Northern Hemisphere, the classification approach has potential for application on smaller lakes too since SAR mode data (~300 m along-track resolution) is used in the study

Machine learning tools for pattern recognition in polar climate science

This thesis explores the application of two novel machine learning approaches to the study of polar climate, with particular focus on Arctic sea ice. The first technique, complex networks, is based on an unsupervised learning approach which is able to exploit spatio-temporal patterns of variability within geospatial time series data sets. The second, Gaussian Process Regression (GPR), is a supervised learning Bayesian inference approach which establishes a principled framework for learning functional relationships between pairs of observation points, through updating prior uncertainty in the presence of new information. These methods are applied to a variety of problems facing the polar climate community at present, although each problem can be considered as an individual component of the wider problem relating to Arctic sea ice predictability. In the first instance, the complex networks methodology is combined with GPR in order to produce skilful seasonal forecasts of pan-Arctic and regional September sea ice extents, with up to 3 months lead time. De-trended forecast skills of 0.53, 0.62, and 0.81 are achieved at 3-, 2- and 1-month lead time respectively, as well as generally highest regional predictive skill ($> 0.30$) in the Pacific sectors of the Arctic, although the ability to skilfully predict many of these regions may be changing over time. Subsequently, the GPR approach is used to combine observations from CryoSat-2, Sentinel-3A and Sentinel-3B satellite radar altimeters, in order to produce daily pan-Arctic estimates of radar freeboard, as well as uncertainty, across the 2018--2019 winter season. The empirical Bayes numerical optimisation technique is also used to derive auxiliary properties relating to the radar freeboard, including its spatial and temporal (de-)correlation length scales, allowing daily pan-Arctic maps of these fields to be generated as well. The estimated daily freeboards are consistent to CryoSat-2 and Sentinel-3 to within $< 1$ mm (standard deviations $< 6$ cm) across the 2018--2019 season, and furthermore, cross-validation experiments show that prediction errors are generally $\leq 4$ mm across the same period. Finally, the complex networks approach is used to evaluate the presence of the winter Arctic Oscillation (AO) to summer sea ice teleconnection within 31 coupled climate models participating in phase 6 of the World Climate Research Programme Coupled Model Intercomparison Project (CMIP6). Two global metrics are used to compare patterns of variability between observations and models: the Adjusted Rand Index and a network distance metric. CMIP6 models generally over-estimate the magnitude of sea-level pressure variability over the north-western Pacific Ocean, and under-estimate the variability over the north Africa and southern Europe, while they also under-estimate the importance of regions such as the Beaufort, East Siberian and Laptev seas in explaining pan-Arctic summer sea ice area variability. They also under-estimate the degree of covariance between the winter AO and summer sea ice in key regions such as the East Siberian Sea and Canada basin, which may hinder their ability to make skilful seasonal to inter-annual predictions of summer sea ice