1. Introduction
Antarctic sea ice plays a key role in the ocean-atmosphere heat, impulse, and matter exchange and is of profound importance for species as living and hatching ground. The Antarctic sea-ice cover is highly variable and—in contrast to its Arctic counterpart—does not yet show a decrease in its overall extent over the satellite era [
1,
2]. One of the key shortcoming in future climate scenario investigations as well as in the understanding of the inter-annual regional Antarctic sea-ice extent variability is given by the still widely unknown Antarctic sea-ice thickness distribution. The lack of this information has led to efforts to, for example, re-construct Antarctic sea-ice thickness distribution by means of numerical modeling (see e.g., [
3]).
Sea-ice thickness information can be obtained most accurately by in situ measurements carried out during expeditions. However, only a limited amount of such measurements is available [
4] and the vast majority of the sea-ice cover remains unsampled. What alternatives do exist? Ship-based visual observations carried out according to the Antarctic Sea Ice Processes and Climate (ASPeCt) protocol [
5,
6] are of limited spatio-temporal coverage as well, as illustrated by the track maps shown in [
7]. Upward looking sonar (ULS) data are—in contrast to the Arctic—only available from moorings [
8,
9] and provide very useful temporal but quite limited spatial information about sea-ice thickness. Only recently, autonomous underwater vehicles (AUV) equipped with ULS devices allowed a first glimpse of Antarctic sea-ice thickness from below also over kilometer scales [
10]. Satellite remote sensing is required to obtain a circum-Antarctic sea-ice thickness distribution.
Synergistic supervised classification and analysis of observations from various remote sensing sensors is used for so-called ice charts which provide, among sea-ice concentration, type and state of development, also information about the sea-ice thickness range. This information has been exploited by [
11,
12] who derived seasonal and annual estimates of the Antarctic sea-ice thickness and volume for 1995–1998. Satellite infrared temperature measurements can be used to obtain sea-ice thickness but only for thin ice regions [
13,
14]. Imaging satellites such as active and passive microwave sensors have also been used to retrieve Antarctic sea-ice thickness. The sea-ice thickness range within which these sensors provide useful information is, however, under debate and is a function of frequency, region, and season [
14,
15,
16,
17,
18,
19,
20,
21]. Passive microwave sensors measure the brightness temperature, which can be used to retrieve primarily the thickness of thin sea ice [
14,
15], i.e., sea ice thinner than ~0.2 m. This is sufficient to estimate ice production in Antarctic coastal polynyas [
16]. The retrieval is based on the relationship between the microwave emissivity of the sea ice and its dielectric properties. These are primarily influenced by the sea-ice salinity which is a function of temperature and thickness of the sea ice. The salinity of thin and first-year sea ice is larger at the surface than in the deeper layers of the sea ice; the salinity profile has a C-shape. The salinity therefore influences the penetration depth of the microwave radiation into the sea ice. This limits sea-ice thickness retrieval using most passive microwave sensors to maximum ~0.2 m. The method of Aulicino et al. [
18], which combines the microwave frequencies used by [
14,
15,
16] in a different way, seems to succeed in retrieving larger sea-ice thickness values. However, this method is limited to values below 1 m and has only been tested for the southern Ross and Weddell Seas. Another alternative is given by the Soil Moisture and Ocean Salinity (SMOS) sensor. The SMOS sensor operates at L-band and its brightness temperatures are routinely used to retrieve sea-ice thickness up to ~0.7 m in the Arctic [
19]; SMOS data are currently also being tested for Antarctic sea-ice thickness retrieval. Active microwave sensors measure the amount of microwave radiation emitted to and scattered back by the sea ice. Like is the case for brightness temperatures, the amount of radiation scattered back is a function of the salinity of the sea ice and therefore of its thickness. For Antarctic sea ice, only a few studies exist which have used active microwave data to estimate sea-ice thickness directly from radar backscatter e.g., [
17]. The limited coverage of high-resolution synthetic aperture radar (SAR) images, which are required for this approach, is one of the main impediments. Other limitations are given by ambiguities in the radar backscatter signal. These two limitations are also valid for another, indirect method to retrieve sea-ice thickness from SAR imagery. This method exploits the fact that the wave spectrum of ocean swell entering the sea ice changes as a function of sea-ice thickness [
20,
21]. The main application area for this method is the pancake ice of the marginal ice zone. The thickness of sea ice into which swell does not penetrate cannot be retrieved with this method.
What remains is satellite altimetry. Two types of sensors exist: radar altimetry to estimate sea-ice freeboard—provided that the radar signal is reflected at the ice-snow interface—and laser altimetry to estimate total (sea ice + snow) freeboard. While radar altimetry has been successfully employed in the Arctic, e.g., Laxon et al. [
22,
23], difficulties arise in the Antarctic [
24]. An unknown fraction of the radar signal is not reflected at the ice-snow interface but somewhere inside the snow layer or even at the snow surface [
25]. Work to improve Antarctic sea-ice freeboard retrieval with radar altimetry is ongoing [
26] but not part of the present study which focuses on laser altimetry.
For laser altimetry, snow properties are less of an issue for the retrieval of the total freeboard and therefore sea-ice thickness retrieval. By using data of the Geoscience Laser Altimeter System (GLAS) aboard the Ice, Cloud and land Elevation Satellite (ICESat) Antarctic sea-ice thickness has been estimated in various regions. Markus et al. [
27] focused on East Antarctic sea ice, Xie et al. [
28] focused on the Bellingshausen and Amundsen Seas, and Zwally et al. [
29], Yi et al. [
30] and Kern and Spreen [
31] focused on the Weddell Sea. Kurtz and Markus [
32] were the first to provide circum-Antarctic sea-ice thickness and volume estimates based on ICESat data. However, there is evidence that the main assumption made in [
32]: sea-ice freeboard is zero everywhere and the elevation measured by ICESat equals snow depth, could have resulted in a considerable sea-ice thickness underestimation [
31,
33]. The main driver for that assumption was to avoid inclusion of snow-depth information into the sea-ice thickness retrieval. Snow depth on Antarctic sea ice is complex, difficult to obtain, and seems to be unbiased only when retrieved from satellite observations over undeformed seasonal sea ice (see e.g., [
34]).
Motivated by these results we carry out a qualitative inter-comparison of sea-ice thickness retrieved from ICESat total freeboard using a number of different approaches—including the data of [
32]. Except this data, all approaches presented in this paper are based on total freeboard retrieved from ICESat data following [
31]. None of these assume that sea-ice freeboard is zero. Instead we apply different ways to treat snow depth on sea ice, e.g., a snow depth climatology, a modified sea-ice density, or different empirical approaches to directly convert total freeboard into sea-ice thickness. At first glance, we find that sea-ice thickness obtained with the empirical approaches seem to agree better with in situ and ship-based observations during winter and spring. However, if we take into account that the observations are potentially biased low, then the “classical” approach for Antarctic sea ice [
29,
30], named SICCI approach [
31] here, and the approach utilizing a modified sea-ice density seem to perform best in comparison to the above-mentioned observations and other independent observations detailed in the paper.
We lay out the approaches in the following section after the description of the data. We present our results in
Section 3 and we discuss these in
Section 4, before we give some concluding remarks in
Section 5.
5. Conclusions
We compute Antarctic sea-ice thickness from ICESat laser altimeter observations which provide total (sea ice + snow) freeboard. We investigate four different approaches, which are based on Archimedes’ principle; only one of these (SICCI) requires actual snow depth information. The other three either use climatological snow depth information (MandC), assume that ICESat freeboard equals snow depth (KandM), or assume a one-layer system in which the impact of snow on sea ice is represented by a modified sea-ice density (Worby). In addition, we investigate empirical approaches (OC2013) which allow computation of sea-ice thickness directly from total freeboard and which neither require snow depth nor densities [
4]. We compute sea-ice thickness for ICESat measurement periods between 2004 and 2008 at 100 km grid resolution. We focus on winter and spring periods. The results are inter-compared with each other and discussed qualitatively in view of the limited amount of independent sea-ice thickness information. Data sets of the sea-ice thickness retrieved with the SICCI and the Worby 1-layer approach are provided as
supplementary material to this article and can also be accessed from the Integrated Climate Data Center (
http://icdc.zmaw.de/1/projekte/eas-sea-ice-ecv0.html). Our main conclusions are as follows:
- (i)
The main limitation of approaches like SICCI and MandC is the lack of accurate reliable snow depths. Available snow-depth data sets are based on satellite microwave radiometry. For these, snow depth is underestimated for deformed sea ice and potentially also for flooded sea ice, i.e., regions of negative sea-ice freeboard. For the snow-depth data set used in this paper, the winter-to-spring snow-depth evolution also does not seem to be correct [
34]. A snow depth which is biased low explains the observed overestimation of sea-ice thickness, particularly for spring, and also the observed too large winter-to-spring increase in modal sea-ice thickness. However, the main limitation of the SICCI approach is at the same time its main advantage. Only by incorporating information about the snow depth which is contemporary to the freeboard observations one can properly resolve the smaller-scale sea-ice thickness variations.
- (ii)
An approach like KandM, where the observed total freeboard is assumed to equal the snow depth (=sea-ice freeboard is assumed to be zero) has the potential to provide more accurate sea-ice thickness estimates than SICCI in regions where (1) the above-mentioned assumption is valid and where (2) the snow depth used by the SICCI approach is biased. We find, however, that the average KandM sea-ice thickness is too small in comparison to the other approaches (
Section 3.1 and
Section 3.2) and to available independent sea-ice thickness observations (
Section 4.5)—in agreement to [
33]. The same applies to the empirical OC2013 approaches (
Section 3.2), except the one denoted OCEA which is based on measurements obtained in the East Antarctic sector (see
Figure 1, right).
- (iii)
Our conceptually new attempt to derive sea-ice thickness, which takes the sea ice-snow system as one layer (Worby), provides sea-ice thickness values which seem to be consistent with independent observations and with regard to the winter-to-spring sea-ice thickness evolution (
Section 3.2 and
Section 4.5). For this approach we assume that the density of the one-layer system with a thickness equaling snow depth plus sea-ice thickness can be described as a linear combination of the densities of snow and sea ice (see
Section 2.2.5).
Currently used sea ice and snow densities need to be reviewed and potentially revised. Recent results from expeditions point to sea-ice densities which might be considerably lower and more variable than thought before [
57]. Sea-ice density can cause considerable systematic uncertainties in the obtained sea-ice thickness as illustrated in
Section 4.4.
The above-mentioned limitations might also apply for sea-ice thickness retrieval planned with ICESat-2 [
52]. While the design of the new sensor aboard ICESat-2 will substantially improve the accuracy of freeboard retrievals and provide improved spatial-temporal coverage [
52], the limitations caused by inaccurate snow depth and other input parameters like snow and ice density could counterbalance the improvements in freeboard retrieval to be achieved with ICESat-2. In this context, we recommend to (1) Improve snow depth on sea ice estimates and provide information about precision and accuracy. We refer to [
59,
60] for first attempts into this direction; (2) Further investigate the potential of using a one-layer system for sea-ice thickness retrieval. For this paper, we utilize climatological values of sea-ice thickness and snow depth from visual ship-based observations [
6] to derive the one-layer ice density. It would be desirable to take this information from more recent observations. Continuation of programs supporting expeditions to carry out co-ordinated ground-based, underwater equipment-based, and airborne measurements of snow and ice properties is mandatory for any further development and improvement of algorithms and the validation of the obtained results.