Search

Information about sea ice surface topography and related deformation is crucial for studies of sea ice mass balance, sea ice modeling, and ship navigation through the ice pack. The Ice, Cloud, and land Elevation Satellite-2 (ICESat-2), part of the National Aeronautics and Space Administration (NASA) Earth Observing System, has been on orbit for over 4 years, sensing the sea ice surface topography with six laser beams capable of capturing individual features such as pressure ridges. To assess the capabilities and uncertainties of ICESat-2 products, coincident high-resolution measurements of sea ice surface topography are required. During the yearlong Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the Arctic Ocean, we successfully carried out a coincident underflight of ICESat-2 with a helicopter-based airborne laser scanner (ALS), achieving an overlap of more than 100 km. Despite the comparably short data set, the high-resolution centimeter-scale measurements of the ALS can be used to evaluate the performance of ICESat-2 products. Our goal is to investigate how the sea ice surface roughness and topography are represented in different ICESat-2 products as well as how sensitive ICESat-2 products are to leads and small cracks in the ice cover. Here, we compare the ALS measurements with ICESat-2's primary sea ice height product, ATL07, and the high-fidelity surface elevation product developed by the University of Maryland (UMD). By applying a ridge-detection algorithm, we find that 16 % (4 %) of the number of obstacles in the ALS data set are found using the strong (weak) center beam in ATL07. Significantly higher detection rates of 42 % (30 %) are achieved when using the UMD product. While only one lead is indicated in ATL07 for the underflight, the ALS reveals many small, narrow, and only partly open cracks that appear to be overlooked by ATL07.

In the Arctic, multi-year sea ice is being rapidly replaced by seasonal sea ice. Baffin Bay, situated between Greenland and Canada, is part of the seasonal ice zone. In this study, we present a long-term multi-mission assessment (2003–2020) of spring sea ice thickness in Baffin Bay from satellite altimetry and sea ice charts. Sea ice thickness within Baffin Bay is calculated from Envisat, ICESat, CryoSat-2, and ICESat-2 freeboard estimates, alongside a proxy from the ice chart stage of development that closely matches the altimetry data. We study the sensitivity of sea ice thickness results estimated from an array of different snow depth and snow density products and methods for redistributing low-resolution snow data onto along-track altimetry freeboards. The snow depth products that are applied include a reference estimated from the Warren climatology, a passive microwave snow depth product, and the dynamic snow scheme SnowModel-LG. We find that applying snow depth redistribution to represent small-scale snow variability has a considerable impact on ice thickness calculations from laser freeboards but was unnecessary for radar freeboards. Decisions on which snow loading product to use and whether to apply snow redistribution can lead to different conclusions on trends and physical mechanisms. For instance, we find an uncertainty envelope around the March mean sea ice thickness of 13 % for different snow depth/density products and redistribution methods. Consequently, trends in March sea ice thickness from 2003–2020 range from −23 to 17 cm per decade, depending on which snow depth/density product and redistribution method is applied. Over a longer timescale, since 1996, the proxy ice chart thickness product has demonstrated statistically significant thinning within Baffin Bay of 7 cm per decade. Our study provides further evidence for long-term asymmetrical trends in Baffin Bay sea ice thickness (with −17.6 cm per decade thinning in the west and 10.8 cm per decade thickening in the east of the bay) since 2003. This asymmetrical thinning is consistent for all combinations of snow product and processing method, but it is unclear what may have driven these changes.

In Release 001 and 002 of the ICESat-2 sea ice products, candidate height segments used to estimate the reference sea surface height for freeboard calculations included two surface types: specular and smooth dark leads. We found that the uncorrected photon rates, used as proxies of surface reflectance, are attenuated due to clouds resulting in the potential misclassification of sea ice as dark leads, biasing the reference sea surface height relative to those derived from the more reliable specular returns. This results in higher reference sea surface heights and lower estimated ice freeboards. The resolution of available cloud flags from the ICESat-2 atmosphere data product is too coarse to provide useful filtering at the lead segment scale. In Release 003, we have modified the surface-reference-finding algorithm so that only specular leads are used. The consequence of this change can be seen in the composites of mean freeboard of the Arctic and Southern oceans. Broadly, coverages have decreased by ∼10–20 % because there are fewer leads (by excluding the dark leads), and the composite means have increased by 0–4 cm because of the use of more consistent specular leads.

We use repeat high-resolution airborne geophysical data consisting of laser altimetry, snow, and Ku-band radar and optical imagery acquired in 2014, 2016, and 2017 to analyze the spatial and temporal variability in surface roughness, slope, wind deposition, and snow accumulation at 88∘ S, an elevation bias validation site for ICESat-2 and potential validation site for CryoSat-2. We find significant small-scale variability ( km) in snow accumulation based on the snow radar subsurface stratigraphy, indicating areas of strong wind redistribution are prevalent at 88∘ S. In general, highs in snow accumulation rate correspond with topographic lows, resulting in a negative correlation coefficient of r2=-0.32 between accumulation rate and MSWD (mean slope in the mean wind direction). This relationship is strongest in areas where the dominant wind direction is parallel to the survey profile, which is expected as the geophysical surveys only capture a two-dimensional cross section of snow redistribution. Variability in snow accumulation appears to correlate with variability in MSWD. The correlation coefficient between the standard deviations of accumulation rate and MSWD is r2=0.48, indicating a stronger link between the standard deviations than the actual parameters. Our analysis shows that there is no simple relationship between surface slope, wind direction, and snow accumulation rates for the overall survey area. We find high variability in surface roughness derived from laser altimetry measurements on length scales smaller than 10 km, sometimes with very distinct and sharp transitions. Some areas also show significant temporal variability over the course of the 3 survey years. Ultimately, there is no statistically significant slope-independent relationship between surface roughness and accumulation rates within our survey area. The observed correspondence between the small-scale temporal and spatial variability in surface roughness and backscatter, as evidenced by Ku-band radar signal strength retrievals, will make it difficult to develop elevation bias corrections for radar altimeter retrieval algorithms.

Icebergs in polar regions affect water salinity, alter marine habitats, and impose serious hazards on maritime operations and navigation. These impacts mainly depend on the iceberg volume, which remains an elusive parameter to measure. We investigate the capability of TanDEM-X bistatic single-pass synthetic aperture radar interferometry (InSAR) to derive iceberg subaerial morphology and infer total volume. We cross-verify InSAR results with Operation IceBridge (OIB) data acquired near Wordie Bay, Antarctica, as part of the OIB/TanDEM-X Antarctic Science Campaign (OTASC). While icebergs are typically classified according to size based on length or maximum height, we develop a new volumetric classification approach for applications where iceberg volume is relevant. For icebergs with heights exceeding 5 m, we find iceberg volumes derived from TanDEM-X and OIB data match within 7 %. We also derive a range of possible iceberg keel depths relevant to grounding and potential impacts on subsea installations. These results suggest that TanDEM-X could pave the way for future single-pass interferometric systems for scientific and operational iceberg mapping and classification based on iceberg volume and keel depth.

In this paper we develop a CryoSat-2 algorithm to retrieve the surface elevation of the air–snow interface over Antarctic sea ice. This algorithm utilizes a two-layer physical model that accounts for scattering from a snow layer atop sea ice as well as scattering from below the snow surface. The model produces waveforms that are fit to CryoSat-2 level 1B data through a bounded trust region least-squares fitting process. These fit waveforms are then used to track the air–snow interface and retrieve the surface elevation at each point along the CryoSat-2 ground track, from which the snow freeboard is computed. To validate this algorithm, we compare retrieved surface elevation measurements and snow surface radar return power levels with those from Operation IceBridge, which flew along a contemporaneous CryoSat-2 orbit in October 2011 and November 2012. Average elevation differences (standard deviations) along the flight lines (IceBridge Airborne Topographic Mapper, ATM – CryoSat-2) are found to be 0.016 cm (29.24 cm) in 2011 and 2.58 cm (26.65 cm) in 2012. The spatial distribution of monthly average pan-Antarctic snow freeboard found using this method is similar to what was observed from NASA's Ice, Cloud, and land Elevation Satellite (ICESat), where the difference (standard deviation) between October 2011–2017 CryoSat-2 mean snow freeboard and spring 2003–2007 mean freeboard from ICESat is 1.92 cm (9.23 cm). While our results suggest that this physical model and waveform fitting method can be used to retrieve snow freeboard from CryoSat-2, allowing for the potential to join laser and radar altimetry data records in the Antarctic, larger (∼30 cm) regional differences from ICESat and along-track differences from ATM do exist, suggesting the need for future improvements to the method. Snow–ice interface elevation retrieval is also explored as a potential to obtain snow depth measurements. However, it is found that this retrieval method often tracks a strong scattering layer within the snow layer instead of the actual snow–ice interface, leading to an overestimation of ice freeboard and an underestimation of snow depth in much of the Southern Ocean but with promising results in areas such as the East Antarctic sector.

Since 2009, the ultra-wideband snow radar on Operation IceBridge (OIB; a NASA airborne mission to survey the polar ice covers) has acquired data in annual campaigns conducted during the Arctic and Antarctic springs. Progressive improvements in radar hardware and data processing methodologies have led to improved data quality for subsequent retrieval of snow depth. Existing retrieval algorithms differ in the way the air–snow (a–s) and snow–ice (s–i) interfaces are detected and localized in the radar returns and in how the system limitations are addressed (e.g., noise, resolution). In 2014, the Snow Thickness On Sea Ice Working Group (STOSIWG) was formed and tasked with investigating how radar data quality affects snow depth retrievals and how retrievals from the various algorithms differ. The goal is to understand the limitations of the estimates and to produce a well-documented, long-term record that can be used for understanding broader changes in the Arctic climate system. Here, we assess five retrieval algorithms by comparisons with field measurements from two ground-based campaigns, including the BRomine, Ozone, and Mercury EXperiment (BROMEX) at Barrow, Alaska; a field program by Environment and Climate Change Canada at Eureka, Nunavut; and available climatology and snowfall from ERA-Interim reanalysis. The aim is to examine available algorithms and to use the assessment results to inform the development of future approaches. We present results from these assessments and highlight key considerations for the production of a long-term, calibrated geophysical record of springtime snow thickness over Arctic sea ice.

We present an analysis of Arctic sea ice topography using high-resolution, three-dimensional surface elevation data from the Airborne Topographic Mapper, flown as part of NASA's Operation IceBridge mission. Surface features in the sea ice cover are detected using a newly developed surface feature picking algorithm. We derive information regarding the height, volume and geometry of surface features from 2009 to 2014 within the Beaufort/Chukchi and Central Arctic regions. The results are delineated by ice type to estimate the topographic variability across first-year and multi-year ice regimes. The results demonstrate that Arctic sea ice topography exhibits significant spatial variability, mainly driven by the increased surface feature height and volume (per unit area) of the multi-year ice that dominates the Central Arctic region. The multi-year ice topography exhibits greater interannual variability compared to the first-year ice regimes, which dominates the total ice topography variability across both regions. The ice topography also shows a clear coastal dependency, with the feature height and volume increasing as a function of proximity to the nearest coastline, especially north of Greenland and the Canadian Archipelago. A strong correlation between ice topography and ice thickness (from the IceBridge sea ice product) is found, using a square-root relationship. The results allude to the importance of ice deformation variability in the total sea ice mass balance, and provide crucial information regarding the tail of the ice thickness distribution across the western Arctic. Future research priorities associated with this new data set are presented and discussed, especially in relation to calculations of atmospheric form drag.

We develop a physical model capable of simulating the mean echo power of CryoSat-2 SAR- and SARIn-mode waveforms over sea-ice-covered regions. The model simulations are used to show the importance of variations in the radar backscatter coefficient with incidence angle and surface roughness for the retrieval of surface elevation of both sea ice floes and leads. The physical model is used to fit CryoSat-2 waveforms to enable retrieval of surface elevation through the use of lookup tables and a bounded trust region Newton least-squares fitting approach. The use of a model to fit returns from sea ice regions offers advantages over currently used threshold retracking methods, which are here shown to be sensitive to the combined effect of bandwidth-limited range resolution and surface roughness variations. Laxon et al. (2013) have compared ice thickness results from CryoSat-2 and IceBridge, and found good agreement; however consistent assumptions about the snow depth and density of sea ice were not used in the comparisons. To address this issue, we directly compare ice freeboard and thickness retrievals from the waveform-fitting and threshold tracker methods of CryoSat-2 to Operation IceBridge data using a consistent set of parameterizations. The purpose of the comparison is to highlight the physical basis between differences in the retracking methods. For three IceBridge campaign periods from March 2011 to March 2013, mean differences (CryoSat-2 – IceBridge) of 0.144 and 1.351 m are found between the freeboard and thickness retrievals, respectively, using a 50% sea ice floe threshold retracker, while mean differences of 0.019 and 0.182 m are found when using the waveform-fitting method. This suggests the waveform-fitting technique is capable of better reconciling the sea ice thickness data record from laser and radar altimetry data sets through the usage of consistent physical assumptions.

The study of sea ice using airborne remote sensing platforms provides unique capabilities to measure a wide variety of sea ice properties. These measurements are useful for a variety of topics including model evaluation and improvement, assessment of satellite retrievals, and incorporation into climate data records for analysis of interannual variability and long-term trends in sea ice properties. In this paper we describe methods for the retrieval of sea ice thickness, freeboard, and snow depth using data from a multi-sensor suite of instruments on NASA's Operation IceBridge airborne campaign. We assess the consistency of the results through comparison with independent data sets that demonstrate that the IceBridge products are capable of providing a reliable record of snow depth and sea ice thickness. We explore the impact of inter-campaign instrument changes and associated algorithm adaptations as well as the applicability of the adapted algorithms to the ongoing IceBridge mission. The uncertainties associated with the retrieval methods are determined and placed in the context of their impact on the retrieved sea ice thickness. Lastly, we present results for the 2009 and 2010 IceBridge campaigns, which are currently available in product form via the National Snow and Ice Data Center.