Arctic sea ice variations are sensitive to Arctic environmental changes and global changes. Freeboard and thickness are two important parameters in sea ice change research. Satellite altimetry can provide long-time and large-scale sea ice monitoring. We estimated the Arctic sea ice freeboard and its variations for the period from 2002 to 2012 from Envisat satellite altimetry data. To remove geoid undulations, we reprocessed the Envisat data using a newly developed mean sea surface (MSS) model, named DTU18. Residuals in the static geoid were removed by using the moving average technique. We then determined the local sea surface height and sea ice freeboard from the Envisat elevation profiles. We validated our freeboard estimates using two radar freeboard products from the European Space Agency (ESA) Climate Change Initiative (CCI) and the Alfred Wegener Institute (AWI), as well as the Operation IceBridge (OIB) sea ice freeboard product. The overall differences between our estimates and the CCI and AWI data were 0.11 ± 0.14 m and 0.12 ± 0.14 m, respectively. Our estimates show good agreement with the three products for areas of freeboard larger than 0.2 m and smaller than 0.3 m. For areas of freeboard larger than 0.3 m, our estimates correlate better with OIB freeboard than with CCI and AWI. The variations in the Arctic sea ice thickness are discussed. The ice freeboard reached its minimum in 2008 during the research period. Sharp decreases were found in the winters of 2005 and 2007.

Satellite altimeters can be used to derive long-term and large-scale sea ice thickness changes. Sea ice thickness retrieval is based on measurements of freeboard, and the conversion of freeboard to thickness requires knowledge of the snow depth and snow, sea ice, and sea water densities. However, these parameters are difficult to be observed concurrently with altimeter measurements. The uncertainties in these parameters inevitably cause uncertainties in sea ice thickness estimations. This paper introduces a new method based on least squares adjustment (LSA) to estimate Arctic sea ice thickness with CryoSat-2 measurements. A model between the sea ice freeboard and thickness is established within a 5 km × 5 km grid, and the model coefficients and sea ice thickness are calculated using the LSA method. Based on the newly developed method, we are able to derive estimates of the Arctic sea ice thickness for 2010 through 2019 using CryoSat-2 altimetry data. Spatial and temporal variations of the Arctic sea ice thickness are analyzed, and comparisons between sea ice thickness estimates using the LSA method and three CryoSat-2 sea ice thickness products (Alfred Wegener Institute (AWI), Centre for Polar Observation and Modelling (CPOM), and NASA Goddard Space Flight Centre (GSFC)) are performed for the 2018–2019 Arctic sea ice growth season. The overall differences of sea ice thickness estimated in this study between AWI, CPOM, and GSFC are 0.025 ± 0.640 m, 0.143 ± 0.640 m, and −0.274 ± 0.628 m, respectively. Large differences between the LSA and three products tend to appear in areas covered with thin ice due to the limited accuracy of CryoSat-2 over thin ice. Spatiotemporally coincident Operation IceBridge (OIB) thickness values are also used for validation. Good agreement with a difference of 0.065 ± 0.187 m is found between our estimates and the OIB results.

The ice shelf is an important component of the Antarctic system, and the interaction between the ice sheet and the ocean often proceeds through mass variations of the ice shelf. The digital elevation model (DEM) of the ice shelf is particularly important for ice shelf elevation change and mass balance estimation. With the development of satellite altimetry technology, it became an important data source for DEM research of Antarctica. The National Aeronautics and Space Administration (NASA) Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) launched in 2018 is a significant improvement in along-track sampling rate and measurement accuracy compared with previous altimetry satellites. This study uses ordinary kriging interpolation to present new DEMs (ICESat-2 DEM hereinafter) for the three ice shelves (Ross, Filchner–Ronne and Amery) in Antarctica with ICESat-2 altimetry data. Two variogram models (linear and spherical) of ordinary kriging interpolation are compared in this paper. The result shows that the spherical model generally shows better performance and lower standard deviation (STD) than the linear models. The precision of the ultimate DEM was evaluated by NASA Operation IceBridge (OIB) data and compared with five previously published Antarctic DEM products (REMA, TanDEM-X PolarDEM, Slater DEM, Helm DEM, and Bamber DEM). The comparison reveals that the mean difference between ICESat-2 DEM of the Ross ice shelf and OIB is −0.016 m with a STD of 0.918 m, and the mean difference between ICESat-2 DEM of the Filchner–Ronne ice shelf and OIB is −0.533 m with a STD of 0.718 m. The three ICESat-2 DEMs show higher spatial resolution and elevation accuracy than five previously published Antarctic DEMs.