We introduce the OSI-450, the SICCI-25km and the SICCI-50km climate data records of gridded global sea-ice concentration. These three records are derived from passive microwave satellite data and offer three distinct advantages compared to existing records: first, all three records provide quantitative information on uncertainty and possibly applied filtering at every grid point and every time step. Second, they are based on dynamic tie points, which capture the time evolution of surface characteristics of the ice cover and accommodate potential calibration differences between satellite missions. Third, they are produced in the context of sustained services offering committed extension, documentation, traceability, and user support. The three records differ in the underlying satellite data (SMMR & SSM/I & SSMIS or AMSR-E & AMSR2), in the imaging frequency channels (37 GHz and either 6 or 19 GHz), in their horizontal resolution (25 or 50 km), and in the time period they cover. We introduce the underlying algorithms and provide an evaluation. We find that all three records compare well with independent estimates of sea-ice concentration both in regions with very high sea-ice concentration and in regions with very low sea-ice concentration. We hence trust that these records will prove helpful for a better understanding of the evolution of the Earth's sea-ice cover.
Sea ice concentration has been retrieved in polar regions with satellite microwave radiometers for over 30 years. However, the question remains as to what is an optimal sea ice concentration retrieval method for climate monitoring. This paper presents some of the key results of an extensive algorithm inter-comparison and evaluation experiment. The skills of 30 sea ice algorithms were evaluated systematically over low and high sea ice concentrations. Evaluation criteria included standard deviation relative to independent validation data, performance in the presence of thin ice and melt ponds, and sensitivity to error sources with seasonal to inter-annual variations and potential climatic trends, such as atmospheric water vapour and water-surface roughening by wind. A selection of 13 algorithms is shown in the article to demonstrate the results. Based on the findings, a hybrid approach is suggested to retrieve sea ice concentration globally for climate monitoring purposes. This approach consists of a combination of two algorithms plus dynamic tie points implementation and atmospheric correction of input brightness temperatures. The method minimizes inter-sensor calibration discrepancies and sensitivity to the mentioned error sources.
We report on results of a systematic inter-comparison of 10 global sea-ice concentration (SIC) data products at 12.5 to 50.0 km grid resolution for both the Arctic and the Antarctic. The products are compared with each other with respect to differences in SIC, sea-ice area (SIA), and sea-ice extent (SIE), and they are compared against a global wintertime near-100 % reference SIC data set for closed pack ice conditions and against global year-round ship-based visual observations of the sea-ice cover. We can group the products based on the concept of their SIC retrieval algorithms. Group I consists of data sets using the self-optimizing EUMETSAT OSI SAF and ESA CCI algorithms. Group II includes data using the Comiso bootstrap algorithm and the NOAA NSIDC sea-ice concentration climate data record (CDR). The standard NASA Team and the ARTIST Sea Ice (ASI) algorithms are put into group III, and NASA Team 2 is the only element of group IV. The three CDRs of group I (SICCI-25km, SICCI-50km, and OSI-450) are biased low compared to a 100 % reference SIC data set with biases of −0.4 % to −1.0 % (Arctic) and −0.3 % to −1.1 % (Antarctic). Products of group II appear to be mostly biased high in the Arctic by between +1.0 % and +3.5 %, while their biases in the Antarctic range from −0.2 % to +0.9 %. Group III product biases are different for the Arctic, +0.9 % (NASA Team) and −3.7 % (ASI), but similar for the Antarctic, −5.4 % and −5.6 %, respectively. The standard deviation is smaller in the Arctic for the quoted group I products (1.9 % to 2.9 %) and Antarctic (2.5 % to 3.1 %) than for group II and III products: 3.6 % to 5.0 % for the Arctic and 4.0 % to 6.5 % for the Antarctic. We refer to the paper to understand why we could not give values for group IV here. We discuss the impact of truncating the SIC distribution, as naturally retrieved by the algorithms around the 100 % sea-ice concentration end. We show that evaluation studies of such truncated SIC products can result in misleading statistics and favour data sets that systematically overestimate SIC. We describe a method to reconstruct the non-truncated distribution of SIC before the evaluation is performed. On the basis of this evaluation, we open a discussion about the overestimation of SIC in data products, with far-reaching consequences for surface heat flux estimations in winter. We also document inconsistencies in the behaviour of the weather filters used in products of group II, and we suggest advancing studies about the influence of these weather filters on SIA and SIE time series and their trends.
We assess different methods and input parameters, namely snow depth, snow density and ice density, used in freeboard-to-thickness conversion of Arctic sea ice. This conversion is an important part of sea ice thickness retrieval from spaceborne altimetry. A data base is created comprising sea ice freeboard derived from satellite radar altimetry between 1993 and 2012 and co-locate observations of total (sea ice + snow) and sea ice freeboard from the Operation Ice Bridge (OIB) and CryoSat Validation Experiment (CryoVEx) airborne campaigns, of sea ice draft from moored and submarine upward looking sonar (ULS), and of snow depth from OIB campaigns, Advanced Microwave Scanning Radiometer (AMSR-E) and the Warren climatology (Warren et al., 1999). We compare the different data sets in spatiotemporal scales where satellite radar altimetry yields meaningful results. An inter-comparison of the snow depth data sets emphasizes the limited usefulness of Warren climatology snow depth for freeboard-to-thickness conversion under current Arctic Ocean conditions reported in other studies. We test different freeboard-to-thickness and freeboard-to-draft conversion approaches. The mean observed ULS sea ice draft agrees with the mean sea ice draft derived from radar altimetry within the uncertainty bounds of the data sets involved. However, none of the approaches are able to reproduce the seasonal cycle in sea ice draft observed by moored ULS. A sensitivity analysis of the freeboard-to-thickness conversion suggests that sea ice density is as important as snow depth.
In order to derive long-term changes in sea-ice volume, a multi-decadal sea-ice thickness record is required. CryoSat-2 has showcased the potential of radar altimetry for sea-ice mass-balance estimation over the recent years. However, precursor altimetry missions such as Environmental Satellite (Envisat) have not been exploited to the same extent so far. Combining both missions to acquire a decadal sea-ice volume data set requires a method to overcome the discrepancies due to different footprint sizes from either pulse-limited or beam-sharpened radar echoes. In this study, we implemented an inter-mission-consistent surface-type classification scheme for both hemispheres, based on the waveform pulse peakiness, leading-edge width, and surface backscatter. In order to achieve a consistent retracking procedure, we adapted the threshold first-maximum retracker algorithm, previously used only for CryoSat-2, to develop an adaptive retracker threshold that depends on waveform characteristics. With our method, we produce a global and consistent freeboard data set for CryoSat-2 and Envisat. This novel data set features a maximum monthly difference in the mission-overlap period of 2.2 cm (2.7 cm) for the Arctic (Antarctic) based on all gridded values with spatial resolution of 25 km × 25 km and 50 km × 50 km for the Arctic and Antarctic, respectively.