1 Introduction
Several sea ice concentration (SIC) algorithms have been developed for passive microwave data (PMW). They differ in the usage of, e.g., frequencies and polarizations of the PMW data () or the usage of static or dynamic tie points (; ; ). From the resulting datasets, it is apparent that the Arctic sea ice extent (SIE) in September has been decreasing at a rate of about 12 % per decade since the launch of modern satellite multi-frequency microwave radiometers in 1978 (; ). This negative sea ice trend in the Northern Hemisphere started in the 1970s (; ), though regional trends can differ, as seen for example within the Barents Sea (). In the Antarctic there are large regional differences in SIE trends, but, until recently, the overall trend was positive due to sea ice dynamics (; ). This has changed, however, in the last decade as a result of several record lows, and, as such, overall trends have shifted to a more homogeneous pattern () and are now slightly negative in summer, November to February (NDJF). Until recently, the slightly positive trend was believed to be part of the long-term natural variability that overshadowed the effects of global warming starting in the 1960s (; ; ; ; ; ). In order to fully understand the drivers of sea ice variability, extending the sea ice data record backwards in time is essential.
Globally, SIE information prior to the satellite data record came largely from ice charts and ship observations. While there have been efforts to include these data in long-term assessments of sea ice change, the data are typically provided at relatively coarse spatial and temporal resolutions (1° grid) (), interpolating in both time and space (). Only satellite-based datasets offer the ability to cover both hemispheres at improved spatial and temporal resolutions and generally show consistency in processing methods (). The SIC derived from NIMBUS 5 Electrically Scanning Microwave Radiometer (ESMR) data was previously processed by . Here, we apply a new processing method that is comparable to the EUMETSAT/ESA CCI SIC record from 1978 and onwards (see ; ). Compared to , this method reduces atmospheric noise regionally over both ice and water surfaces and uses the pre-processed data to develop a SIC algorithm calibration that is effective in removing both instrument drift and offsets. Seasonal sea ice signature variations are removed by using dynamical tie points. Lastly, the algorithm calculates temporally and spatially varying uncertainty estimates. The ESMR SIC data are presented on the same grid and with the same masking as the EUMETSAT/ESA CCI record, which makes these two records directly comparable. This and the modern processing chain mentioned above warrant the reprocessing presented in this article.
In Sect. 2, the satellite and reanalysis data are described, including the formatting and initial filtering of the data. Section 3 introduces the radiative transfer model (RTM) used for the atmospheric correction, while Sect. 4 describes the dynamical tie points, the SIC algorithm with uncertainty estimations, the land-spillover method, and data flags assigned during post-processing. In Sect. 5, the resulting SIC dataset is presented and compared to other datasets. Finally, Sect. 6 consists of a discussion, and Sect. 7 provides the conclusions of this work.
2 The NIMBUS 5 ESMR instrument and data
The NIMBUS 5 ESMR instrument was a cross-track scanner that measured at 78 scan positions perpendicular to the flight track with a maximum incidence angle of about 64° on both sides (). No direct observations at nadir were made; the closest positions were at °. Its near-circular orbit had a height of about 1112 km and an inclination of 81°. The phased-array antenna dimensions were cm, and the spatial resolution was about 25 km near nadir, increasing to about km at the edges of the swath (). The full swath was about 3100 km, with the varying incidence angle and the spatial resolution giving very good (unprecedented) daily coverage in polar regions with no gaps, i.e., no pole holes. The ESMR onboard the NIMBUS 5 satellite was a one-channel 19.35 GHz horizontally polarized microwave radiometer operating from 11 December 1972 until 16 May 1977 (1617 d) with some interruptions (see the list of days with missing files in Appendix A2). Due to a hot-load anomaly, there are major data gaps between March to May and again in August 1973. Another major data gap occurs from 3 June until 14 September 1975 because the ground segment was used for receiving NIMBUS 6 data instead (). When operation was resumed in September 1975, the instrument was only operated approximately every other day. From late 1976 to the end of the mission, operation was highly irregular. The last file in the dataset is from 16 May 1977. The data have recently been made available online by NASA in the original tape archive format (TAP files).
2.1 Formatting and co-location of brightness temperatures and ECMWF ERA5 data
The ESMR data were retrieved from the online data archive of the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) (). This dataset contains, along with a number of instrument and geographical parameters, calibrated 19.35 GHz brightness temperatures expressed in units of kelvin. The raw data were recovered by NASA from the magnetic tapes, called calibrated brightness temperature tapes (CBTTs), where they were stored in the original binary TAP file format, each file corresponding to a particular orbit ().
All variables in the TAP files were read using online NASA software and converted to NetCDF format without changing the original data structure, creating raw ESMR NetCDF files. Each data point in the TAP file was matched with the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 re-analysis data (; ) that were nearest in time and space and appended to the raw ESMR NetCDF file serving as input to the processing chain. The resulting data are structured in arrays line by line (across track). Appendix summarizes the variables included in the NetCDF files.
2.2 Initial filtering and correction of brightness temperatures
NASA provides a correction of the brightness temperature data to account for the lobe structure, antenna loss, and angular brightness temperature () variation (). According to NASA, the correction was needed because “The cause of the gross variations in antenna properties which were observed soon after launch has been determined to be a cross-polarized grating lobe […] The problem does not exist for the near-nadir beam positions so these positions are unaffected. […] An empirical calibration has been developed which removes the effect of the lobe structure and antenna loss, which vary with position, and roughly corrects for angular variations in viewing geometry.” (, p. 400).
Originally, it was planned to use only lobe-corrected values with their natural angular dependency, but we did not find a way to extract this from the NASA-provided dataset. Essentially, only the combined lobe and angular correction, which is a function of brightness temperature, can be removed altogether from the data NASA provides. Thus, the values do not vary as a function of incidence angle, as would be expected for values from the sea surface and sea ice.
Despite the corrections done by NASA, the ESMR data still contain erroneous values, scan lines, sudden jumps in the calibration, and other obvious artifacts.
Since the ESMR data contain corrupted data and erroneous scan lines, filtering is needed before the data can be used for sea ice mapping. The filters that we apply are described in Eqs. (1)–(4). They are applied in the same order as described here. If only a single data point or a couple of scan lines are affected, only those data points and scan lines are removed from the file. If the whole file is corrupted then it is deleted.
An initial analog filter is used for filtering erroneous values and scan lines. The filter is based on the 16 analog voltage entries in the data. The NIMBUS 5 ESMR user's guide () does not explain very well what the 16 entries really are, or what unit the voltages are stored in, but jumps in these analog signals correspond to anomalous values. Our analog filter computes the absolute gradient in the analog signals, and anything over a value of 10 is removed. This threshold was estimated experimentally.
Next, several other filters are employed using the processed (in K) from the previous step. The filters are applied in the following order.
Data that are non-physical and outside the expected range for sea and ice surfaces are removed. Only data points that lie inside the range specified in Eq. () are kept:
1 The next filter removes erroneous scan lines (across track rows). Consecutive scan lines should not differ by more than 50 K, as shown in Eq. (). The threshold of 50 K was estimated experimentally. 2 where in Eq. () is an across-track row of values and is an index along the track, while is an index across the track. is the maximum across-track index of a row, i.e., for a complete row with valid data points for all 78 incidence angles, .
Afterwards, single outliers are removed if they do not satisfy Eq. (): 3 where is a single-pixel and is an along-track index. The threshold of 150 K was selected manually after identifying erroneous single-pixel outliers in the data. The last filter in Eq. () removes neighboring values which are locked to the same value, i.e., values for which the following equation equals zero: 4 where is again a single-pixel with an along-track index . The decision to compare seven consecutive values was made based on qualitative experiments. Since the filter is used universally for all incidence angles, the search window varies but covers a minimum distance of 175 km.
The outer data points of the swath edges show significantly higher noise levels and a coarser resolution than the near-nadir data points (). Therefore, after the filtering, we additionally remove the four outermost data points of the swath, corresponding to incidence angles between 57° and 64° on both sides. The new outer edge of the swath data is then at °, similar to that of modern microwave radiometers used for sea ice retrieval: 50–55° (AMSR, ; SMMR, ; SSM/I, ).
Before the filtering, the dataset contained 13 496 orbital data files; after filtering, there are 10 649 () good files left. A complete list of days where data are missing after filtering and no SIC could be retrieved is given in Appendix .
3 The radiative transfer modelThe radiative transfer model (RTM) was developed specifically for atmospheric noise reduction and is comparable to the RTMs used in and but with the addition that this ESMR RTM (Eq. ) can be applied for different incidence angles over both ocean and ice. The RTM takes as input the atmospheric columnar water vapor [mm or kg m], the 10 m wind speed [m s], the atmospheric columnar cloud liquid water [mm or kg m], the sea surface temperature [K], the ice emitting layer temperature [K], the sea ice concentration [0–1], and the incidence angle [°]. In return, it simulates the top-of-the-atmosphere 19.35 GHz at horizontal polarization:
5
The RTM uses the atmospheric part of the model described in to compute the atmospheric emission, transmissivity, and reflectivity at the sea surface (open water) together with added modules for simulating the sea ice emissivity (Fig. ) and open water reflectivity as a function of incidence angle.
Figure 1
The first-year ice emissivity used in the RTM.
[Figure omitted. See PDF]
For the sea ice emissivity, a look-up table is produced from a simulation using a combined sea ice thermodynamic and emission model during the Arctic winter on first-year ice. The thermodynamical and emission model setup and the simulations are described in . The emissivities as a function of incidence angle are shown in Fig. , and the look-up table is given in Table .
Table 1Sea ice emissivity look-up table ( is the incidence angle and is the first-year ice emissivity).
0 | 5 | 10 | 15 | 20 | 25 | 30 | 35 | 40 | 45 | 50 | 55 | 60 | 65 | |
0.91 | 0.91 | 0.91 | 0.91 | 0.90 | 0.90 | 0.89 | 0.89 | 0.88 | 0.86 | 0.85 | 0.83 | 0.80 | 0.77 |
Sea water permittivity, which is used to estimate the calm sea reflectivity (Eq. ) as a function of temperature, is computed using equation E64 (p. 2046) in . The permittivity is almost invariant with respect to the water salinity at 19 GHz, and a constant value of 34 ppt is used here.
The calm sea (Fresnel) power reflection coefficient, , as a function of the relative permittivity, , and the incidence angle, , for a lossy medium is computed using Eq. (1.52) in , i.e., 6 where and are abbreviations for: 7 and 8 The relative permittivity of the water surface is a complex number. The calm sea emissivity, , is then 9
The rough water surface emissivity component, , is added to the calm sea emissivity, , to produce the total sea water emissivity, . Between ESMR incidence angles of 0 and 63°, the sensitivity of the 19.35 GHz horizontal polarization, , to wind speed is an almost linear function of the incidence angle (), i.e., 10 and 11 where is the incidence angle in degrees, is the wind speed, and is the sea surface temperature [K].
This combination of and follows the procedure described in .
The resulting brightness temperature is a linear combination of the sea water and sea ice emission weighted by the SIC following : 12 where is the upwelling brightness temperature from the atmosphere, is the atmospheric transmissivity, is the water surface emissivity, is the sea ice emissivity, is the reflection reduction factor due to water surface roughness, is the downwelling brightness temperature, and is the cosmic background brightness temperature (2.7 K).
EMSR-simulated values and emissivities have been compared with other values simulated using other RTMs for AMSR (), SMMR (), and SSM/I () for a constant incidence angle of 55°, which is close to the incidence angles of other instruments and RTMs (53–55°). The comparison showed that the values of the ESMR RTM are within the range (within approx. 2 K) of values obtained with the other models and therefore seem to be reasonable given the differences in instrument center frequencies and measurement geometries. Note that in the correction procedure, the difference between two simulated values is used to minimize model biases. Even if the absolute values of the RTM-simulated values were biased, this bias would be removed by taking the difference between two simulated values, which is the only part used in the correction.
4 Derivation of the sea ice concentrationThe RTM is an essential part of the SIC derivation for applying an atmospheric noise reduction. The following section presents the calculations of dynamical tie points and the SIC algorithm along with uncertainty estimations. Lastly, the post-processing, including the land-spillover method and data flag assignment, is described. A flow chart illustrating this processing chain can be found in the ESA CCI ESMR product user guide (PUG) ().
4.1 Tie points and geophysical noise reduction
Tie points are typical signatures of 100 % ice and open water (0 % ice) and are used in SIC algorithms as a reference for estimating the total ice fraction per satellite pixel . The aim of using dynamical tie points is to reduce SIC biases that may result from seasonal and interannual variations in (), instrument drift, and RTM and ERA5 biases. For example, argue that the variations in the average open-water values near the ice edge are affected primarily by variations in instrument calibration, and they describe the drop followed by a sharp increase seen in 1975 (see Fig. ) as an instrument drift issue.
Our ESMR tie points are derived on a daily basis from the swath files. Regions of open water and high SICs are selected for each hemisphere, resulting in two regions of sea ice and two of open water for both hemispheres. The ERA5 SIC prior to 1979 is based on the HadISST2.0.0.0 dataset (), which mainly utilizes digitized sea ice charts for this period (). The two main data sources are the Walsh dataset (; ; ) and National Ice Center (NIC) charts (). The datasets also consist of several data types besides ice charts, e.g., ship observations and satellite data (both infrared and microwave observations, including data from ESMR). The selection of the four tie points is based on a criteria set for SIC from ERA5 – the distance from the ice edge, the observed brightness temperature, the latitude, and the sea surface temperature – as shown in Table . The distance from the ice edge criterion is imposed by putting a threshold on the mean SIC of a grid point box; this should be larger than 80 % for ice tie points or less than 1 % for open water points. While computed daily, these are subsequently combined into 15 d running-mean tie points that are 7 d ahead and 7 d behind the processed date shown in Fig. . The 15 d averaging period is maintained even at the beginning and end of the dataset and when there are data gaps; i.e., if there are gaps, the averaging is done over the days available in the 15 d period (this is also the case in the beginning, where the first 7 d are missing). Figure depicts the 15 d averaged tie points through time. It shows that the ice tie points follow a seasonal pattern, while the water tie points are relatively constant. The tie-point criteria from Table ensure that each daily tie point is based on many observations, which results in stable tie points. The 15 d interval was chosen experimentally, so the variations seem reasonable and one is still able to identify calibration issues such as jumps (e.g., that seen in 1976).
Table 2
Criteria for the four different tie points.
Ice | Ocean | |||
---|---|---|---|---|
Arctic | latitude | latitude | ||
sea ice concentration (ERA5) | sea ice concentration (ERA5) | |||
mean sea ice concentration (ERA5) of | mean sea ice concentration (ERA5) of | |||
a grid point box | a grid point box | |||
100 K brightness temperature K | sea surface temperature (ERA5) K | |||
90 K brightness temperature K | ||||
Antarctic | latitude | latitude | ||
sea ice concentration (ERA5) | sea ice concentration (ERA5) | |||
mean sea ice concentration (ERA5) of | mean sea ice concentration (ERA5) of | |||
a grid point box | a grid point box | |||
100 K brightness temperature K | sea surface temperature (ERA5) K | |||
90 K brightness temperature K |
Figure 2
The 2-weekly tie points for Arctic and Antarctic ice and water after correction. The boxes show the period May–July 1976, when there were obvious instrument calibration issues.
[Figure omitted. See PDF]
Figure 3
Mean atmospheric water vapor for all grid points included in the four tie points. Water vapor data are from ERA5 ().
[Figure omitted. See PDF]
The per-grid-point correction term is the difference between a simulated reference obtained using mean values of the total column water vapor [kg m] in the atmosphere (), 10 m wind speed [m s] (), total column cloud liquid water [kg m] in the atmosphere (), sea surface temperature (), and ice emitting layer temperature () as input to the RTM and a simulated obtained using the actual ERA5 values (, , , , and ) for the grid point. is not corrected for the cloud liquid water, , so the mean is input to both the reference and the actual simulation. can both be negative and positive, and the values have a reduced sensitivity to the geophysical noise sources (, , , ) after correction. The fact that the correction term is the difference between two RTM simulations minimizes the impact of biases in the model and the ERA5 data.
The correction term is added to the measured , i.e., 13 Here, 14 where is the open water tie point and the ice tie point. Following , is computed as 15 where is the 2 m air temperature, which is taken from the ERA5 data. Horizontal bars above the variables indicate that they are daily mean values for the cluster of points selected for the tie point. The mean water vapor, , at the tie point is shown in Fig. .
Figure shows the correction term, , for 1 January 1974 over open water in the Northern and Southern hemispheres, respectively. The path length through the atmosphere is longest at high incidence angles and shortest near nadir, and thus the absolute value of the correction is largest at high incidence angles. For example, when the atmosphere is driest in the reference compared to the actual simulation, the ends of the corrected turn negative, while they turn strongly positive when the reverse is true.
Figure 4
Difference between the values before and after correction with a mean reference for open water tie points in the (a) Northern Hemisphere and (b) Southern Hemisphere.
[Figure omitted. See PDF]
The correction works best over open water areas, where it acts as only an atmospheric correction. The RTM appears to better simulate the relevant emission processes in the atmosphere, and the ERA5 data more accurately quantify the atmospheric noise sources. Over sea ice, geophysical noise sources are related to processes in the snow and ice profile () which are not characterized by the RTM, except for the emitting layer temperature . , which is used as input to the RTM, is estimated from the 2 m air temperature in the ERA5 data using Eq. (). This is important because the ESMR is a single-channel instrument and thus the and also the derived are sensitive to the emitting layer temperature.
The standard deviations of the brightness temperature for water points in both hemispheres before and after the correction are shown in Fig. .
Figure 5
Standard deviations of values before and after correction in January 1974 for the (a) Northern Hemisphere and (b) Southern Hemisphere. The plots only show ocean points which fulfill the ERA5 SIC and sea surface temperature (SST) criteria described in Table 2.
[Figure omitted. See PDF]
4.2 The sea ice concentration (SIC) and its uncertaintySIC () is estimated using the measured brightness temperature () and the open water () and ice () tie points, i.e.,
16
Because the RTM requires as input, is processed iteratively in two steps:
-
The is first estimated using uncorrected values and tie points derived from uncorrected data. The estimate is truncated to the interval between 0 and 1 and an open water filter is applied, forcing all values of less than 0.15 to 0.
-
The estimate from step (1) is used in the RTM calculation (Eq. ) together with ERA5 data for the geophysical noise reduction of the values, and is then estimated again in a second iteration, this time using the corrected values and corrected tie points. The mean values of used in the reference simulation are averages weighted by the obtained using the mean water and ice tie-point values respectively, i.e., is used as a ratio to mix the two tie-point values to create mean values of the numerical weather prediction (NWP) data for any sea ice concentration.
The total SIC uncertainty is the combination of two components: (1) the algorithm uncertainty, which includes instrument noise and tie-point variability (geophysical noise), and (2) the resampling uncertainty, which is the uncertainty due to data resampling.
The algorithm uncertainty is the squared sum of three independent components following : 17 where the first term in Eq. () represents variations due to instrument noise, estimated to a brightness temperature error of 3 K ().
Without the instrument noise term, which is already included in the two tie-point uncertainties, the second and third terms in Eq. () are used to compute the algorithm uncertainty, : 18 where is the water tie-point error (1 standard deviation of the daily tie point here) and is the ice tie-point error (e.g., 1 standard deviation of the daily tie point). The water and ice tie-point errors are weighted by the SIC, and all three errors are normalized with the ice–water brightness temperature contrast and the 2-weekly tie points. The algorithm uncertainty is computed based on swath data.
The resampling uncertainty is the maximum difference - minimum in a pixel window. The resampling uncertainty is computed based on resampled data (e.g., ).
The total uncertainty is the squared sum of the algorithm and resampling uncertainties, i.e., 19
The two uncertainty components and the total uncertainty are included in the data file.
4.3 Land-spillover correction and post-processingThe land-spillover correction follows the procedure described in . A 5 by 5 pixel neighborhood of the land mask (EASE2 version 2, by OSI-SAF) is analyzed to determine which coastal points should be corrected. The land mask is divided into two classes: land points, which are given a value of 90 % SIC, and open ocean points. If the difference between the original land mask and the mean mask calculated for the 5 by 5 window is larger than the previously estimated SIC (the RTM-corrected and resampled SIC), i.e., the SIC is smaller than the theoretical value of the land spillover, the SIC value is set to 0 % and the status flag variable of the dataset is raised to 8.
Additionally, a monthly climatology (also from OSI-SAF; the same version as for the land mask) is used to set SIC to 0 % and to mark open water points by a climatology boundary, which is indicated by a status flag value of 64. Afterwards, the land mask is also used to mark lakes and coastal areas with status flags 2 and 32, respectively. An overview of all status flag values is shown in Table .
Table 3
Description of values of the status flag variable of the dataset.
no flag/flag 0 | Nominal retrieval by the SIC algorithm |
flag 1 | Land |
flag 2 | Lake |
flag 4 | SIC is set to zero by the open water filter |
flag 8 | SIC value is changed to correct for land spillover effects |
flag 16 | Handle with caution: the 2 m air temperature is high at this position, and this might be false ice |
flag 32 | Coast |
flag 64 | SIC is set to zero since the position is outside the maximum sea ice climatology |
flag 128 | Point is not accepted but no other flags are raised |
The results of the post-processing are included in the daily NetCDF files for the Northern and Southern hemispheres, respectively.
5 ResultsA list of all output variables in the daily SIC files and a short description of them can be seen in Appendix . Examples of monthly means of the SIC and the mean uncertainty can be seen in Figs. and .
It is worth noting that the coverage in Fig. is complete and, because of the ESMR's wide swath width of 3100 km and its inclination, the North Pole is covered, in contrast to the satellite microwave radiometers that have followed the NIMBUS 5 ESMR, which have a “pole hole”. The area covered by multi-year ice in the central Arctic has a lower SIC than the first-year ice regions. This is a consequence of the one-channel SIC algorithm, which has inherent ambiguity in SIC, ice type, and emitting layer temperature variations. In Figs. b and b, it can be seen that the uncertainties are largest near the ice edges, as expected. This is due to the resampling uncertainty, which dominates near the ice edge, where the spatial variability of is high. Coastal regions also show higher uncertainties for this reason, since the land-spillover correction is first applied to the SIC after the uncertainty estimations.
Figure 6
Monthly mean (a) SIC and (b) uncertainty for February 1973 in the Northern Hemisphere. Water areas with no ice and uncertainty due to the ocean climatology are displayed in dark blue.
[Figure omitted. See PDF]
Figure 7
Monthly mean (a) SIC and (b) uncertainty for October 1974 in the Southern Hemisphere. Water areas with no ice and uncertainty due to the ocean climatology are displayed in dark blue.
[Figure omitted. See PDF]
The SIC shows interesting sea ice features in the years 1972–1977. One such feature is the Odden ice tongue () extending eastward from the East Greenland Current, which is visible in Fig. (around 73° N, 0° E), while another feature is the Maud Rise Polynya (), an open water area encircled by sea ice in the Southern Hemisphere, which can be seen in Fig. (around 65° S, 0° E). Both examples were much larger in extent in the 1970s and occurred more frequently then than they do today (; ; ).
The daily coverage of valid data points that passed all filtering varied a lot through the ESMR's operating period. While there was nearly full coverage for the first few months, it got much worse after the summer of 1975, when the instrument only recorded data every second day. An example of the poor coverage is shown for May 1976 in Fig. .
Figure 8
Example of poor monthly coverage in May 1976 in the Northern Hemisphere for (a) the mean SIC and (b) the number of days with valid values. Water areas that showed full daily coverage due to the ocean climatology are displayed in dark blue to avoid misleading comparisons. The number of days with valid data is indicated by the colorbar.
[Figure omitted. See PDF]
Monthly averaged SICs are derived to compare our results against other datasets. Only months with 99 % coverage were used in the comparison, i.e., 99 % of all grid points are covered at least once per month. Monthly mean SIEs are calculated from the monthly SICs using the threshold . In Fig. , the ESMR dataset (orange line) is shown together with the OSI-SAF climate data record (CDR; blue line) for 1979–2022 () and the sea ice extent derived from the NSIDC's NIMBUS 5 ESMR ice concentration product (green line) (), using the same threshold for all products ().
The comparison shows comparable SIE levels around 1980. In general, our ESMR dataset has slightly higher monthly SIE values than the NSIDC's ESMR product, even though the seasonal pattern is the same.
The mean difference between our sea ice extent and that from NSIDC is 0.24 million km in the Arctic and 0.59 million km in the Antarctic for the whole dataset.
For the Northern Hemisphere, the SIE during the operational period of NIMBUS 5 ESMR (1972 to 1977) seems to have been similar in magnitude to the SIE during the operational period of NIMBUS 7 SMMR (from 1978 to 1987), with the ESMR's minimum extents being slightly higher than the SMMR's ones. In the Southern Hemisphere, the values for the second half of the 1970s seem to have been around the same magnitude as the largest SIE during the 2014/2015 season.
Figure 9
Monthly sea ice extent time series for the (a) Arctic and (b) Antarctic based on the threshold. The orange curve shows values in the ESMR dataset which have 99 % coverage of the hemisphere monthly, while the blue curve is based on the SIC products from the OSI-SAF (), where the threshold was applied as in . The green line represents the NSIDC's ESMR SIC product ().
[Figure omitted. See PDF]
6 DiscussionComparisons between different sea ice products and the new ESMR dataset proved to be more difficult than initially expected, since not only the processing algorithms but also the land masks, map projections, and dataset grids differ. We were not able to find two independent SIC datasets for 1978 onwards and 1972–1977 which share exactly the same land mask.
Thus, it was decided at the beginning of the processing to use the same land mask as the OSI-420 product (1978 onwards) () for our ESMR dataset, i.e., a 25 km equal-area grid (EASE-2 version 2) land mask, to at least ensure a fair comparison between these two datasets. The NSIDC ESMR dataset (green line in Fig. ) used a different land mask with a polar stereographic projection (), but it was still compared to our ESMR dataset.
The difference in SIE is therefore also influenced by the different projections of the data; however, the area difference between the projections is relatively small (only a few thousand km), so even a re-projection is expected to yield minimal differences compared to the differences caused by the use of different land masks. The comparison of different land masks is complicated by the varying sea ice extent, which exposes more or less land throughout its annual cycle and thus changes the number of grid points affected by the land mask.
The land mask land area differs between the OSI-SAF and NSIDC ESMR land masks. A comparison of land and ocean points between the NSIDC land mask and the OSI-SAF land mask showed a difference for the Northern Hemisphere of 460 000 km (north of 60° N), where the NSIDC land mask has more land, while the difference is the opposite and much smaller in the Southern Hemisphere – only 79 000 km (south of 60° S) – with the OSI-SAF's land mask having slightly more land. More land points in the land mask result in fewer available grid points for potential sea ice. The difference in Southern Hemisphere sea ice extent is significantly larger and opposite to the expected contrast from the land mask difference. In the Northern Hemisphere, it is not very clear how much of the SIE difference can be accounted for by the land mask or by algorithm differences. However, since the SIE differences vary a lot, and the differences in the Southern Hemisphere clearly cannot be explained by the land mask difference alone, it is likely that most of the SIE differences come from the algorithms and processing methods, such as the atmospheric correction and tie point calculation.
The blue curve in Fig. is based on OSI-SAF's SIC products OSI-450 and OSI-430-b (). The shown extend corresponds to the SIE of the OSI-420 product () but with a 30 % instead of a 15 % sea ice threshold applied, matching the OSI-402-d sea ice extent product ().
Compared to the more common 15 %, the 30 % threshold was better suited to a comparison between different ESMR SIE datasets due to the relatively high noise level, which can be seen from the total uncertainty in Fig. . The uncertainty algorithm has been applied for easier data assessment and comparison to other datasets.
A large number of ESMR data are currently filtered out, and the 99 % threshold for inclusion in the monthly timeline in particular filters out the second half of the ESMR data, where large data gaps occurred, as seen in Fig. .
The filters worked as expected and removed erroneous values from the raw data. A reprocessing of the data to rescue more of the currently filtered out data points from the 20 % of ESMR data files is planned. Two approaches will be used to improve the data selection during the next processing. One approach consists of making changes to the set of filters, i.e., performing an adjustment of the filter thresholds, testing different filters, and the possible incorporation of an incidence angle dependency into the data selection. The other approach is to recalibrate some of the erroneous data files that have shown some systematic offsets to rescue whole swaths.
To reduce the uncertainty caused by atmospheric noise, the brightness temperatures were corrected with an RTM that used several atmospheric parameters from NWP (ERA5) data, such as water vapor and wind, as input. This correction showed a consistent reduction of the standard deviation of the brightness temperatures for water points in both hemispheres, as can be seen from Fig. . The correction was less steady over ice surfaces, since the RTM did not describe all relevant processes related to the snow and ice, which are the main noise source over sea ice. Such processes include sea ice deformation, the creation of leads or ridges, and changes in the snow layer, e.g., changes in snow depth, snow density, and grain size but also melting and refreezing, which influences scattering and emission processes inside the layer. Atmospheric noise caused by, e.g., water vapor and cloud liquid water influences the over sea ice, but its influence is smaller: only around one-third of the total noise (). By correcting for atmospheric effects with ERA5 data, we might have introduced some noise into the angular dependency of the SIC due to the use of an incident-angle-dependent emissivity in the RTM (Fig. and Table ).
To avoid biases from the RTM and the NWP data, dynamical tie points have been used, which also calibrate the algorithm for seasonal variations and instrument drift. However, we currently use mean tie points that are independent of the incident angle. Therefore, a possible improvement to a future version of the dataset might be accomplished by using angle-dependent tie points instead.
Even after filtering the data for obvious errors, it is clear that there are still issues with the absolute calibration of the instrument (). For example, after the hot-load anomaly in 1973, the ocean in the Southern Hemisphere is several kelvin below the level before the anomaly, and in 1976 there is a dip in May and June followed by a sharp increase in (). Low-frequency (with timescales days) variations and regional variations at hemispherical scales are compensated for by the dynamical SIC algorithm tie points ().
In spite of data gaps and calibration issues, the experimental NIMBUS satellite program was very successful. Applying modern processing methodologies, including dynamical tie points and atmospheric noise reduction of the values, reduces the noise over both ice and open water consistently. This newly processed ESMR sea ice dataset extends the existing sea ice climate data record (CDR) with an important period from the 1970s. This extension of the SIE record contributes to the United Nations Sustainable Development Goals (SDGs) related to climate change by providing more observations for longer-term assessments of Arctic and Antarctic sea ice changes.
7 Data availability
The newly processed ESMR data have been released through the ESA CCI Open Data Portal:
8 Conclusions
In this paper, we have presented a new SIC dataset covering the period 1972–1977, obtained using the ESMR data from the NIMBUS 5 satellite. The dataset consists of resampled daily NetCDF files for the Northern and Southern hemispheres. The SIC, associated uncertainties, and processing flags are included in the dataset. The uncertainties follow the same principles as those of the EUMETSAT SIC CDR, including both algorithm and resampling uncertainties. The choice to use the same land mask, spatial grid, and projection as for EUMETSAT's SIC CDR make comparisons between the time periods easier.
A comparison to NSIDC's ESMR SIC product and the OSI-SAF CDR showed that the seasonal pattern is very similar to NSIDC's ESMR SIC product, but our product shows systematically larger SIE values which can not be explained by differences between land masks alone. For the Northern Hemisphere, our SIE values match the levels in the 1980s from the OSI-SAF CDR with the same land mask, while values in the Southern Hemisphere in the 1970s are larger than those in the 1980s.
Compatibility with the EUMETSAT's SIC CDR was achieved by using a similar processing chain. The processing included atmospheric noise reduction with the use of an RTM and the ERA5 atmospheric data, which lowered the standard deviation of the s consistently. Additionally, dynamical tie points were used to avoid biases from the RTM and NWP data as well as to adjust for seasonal variability and instrument biases. To ensure better data assessment in analysis and models and easier comparison to other datasets, temporally and regionally varying uncertainty estimates have been included in our ESMR dataset.
Appendix A Table A1The data variables in the NetCDF processing input file and a description of each variable.
Satellite variables | |
---|---|
Time | Time of data [year, month, d, h, min, s] |
Brightness_temperature | Brightness temperature of the 78 scan spots [K ] |
Latitude | Latitude of the 78 scan spots [° ] |
Longitude | Longitude of the 78 scan spots [° ] |
Pitch_fine_error | Pitch fine error [° ] |
Roll_fine_error | Roll fine error [° ] |
RMP_rate | Rate Measuring Package (RMP) indicated rate high [] |
NADIR_LAT | Sub-satellite latitude [° ] |
NADIR_LON | Sub-satellite longitude [° ] |
Height | Satellite height [km] |
Digital_b | A set of 1 bit status words to indicate the positions of each of the command relays (users guide p. 83) |
Status_indicator_1 | A bit status word |
Status_indicator_2 | A bit status word |
Data_source | A bit status word |
Beam_position | A bit status word |
PGM_id | Unique identification number assigned to program that prepared tapes |
HOT_MEAN | Hot load mean [] |
HOT_RMS | RMS of hot load [] |
COLD_MEAN | Cold load mean [] |
COLD_RMS | RMS of cold load [] |
MUX_1 | Average antenna temperature |
MUX_2 | Average phase shifter temperature |
MUX_3 | Ferrite switch temperature |
MUX_4 | Ambient load temperature |
MUX_5 | Reference load temperature |
MUX_6 | Automatic gain control |
Analog_0 | Analog signals (voltages) |
Analog_1 | Analog signals (voltages) |
... | ... |
Analog_15 | Analog signals (voltages) |
ERA5 variables | |
u10 | component of the wind speed at 10 m (parallel to latitude) [m s] |
v10 | component of the wind speed at 10 m (parallel to longitude) [m s] |
t2m | 2 m air temperature [K] |
istl1 | Internal temperature of ice [K] |
... | ... |
istl4 | Internal temperature of ice [K] |
lsm | Land–sea mask |
msl | Mean sea level pressure [hPa] |
siconc | Sea ice concentration [0–1] |
sst | Sea surface temperature [K] |
skt | Skin temperature [K] |
tcw | Total column water [kg m] |
tcwv | Total column water vapor [kg m] |
era_time | Valid time for analysis |
The output variables stored in the daily NetCDF files and a description of each variable.
ice_conc | Filtered sea ice concentration obtained using atmospheric correction of brightness |
temperatures and open water filters [%] | |
raw_ice_conc_values | Raw sea ice concentration estimates as retrieved by the algorithm [%] |
total_standard_error | Total uncertainty (one standard deviation) of sea ice concentration [%] |
smearing_standard_error | Smearing uncertainty of sea ice concentration [%] |
algorithm_standard_error | Algorithm uncertainty of sea ice concentration [%] |
status_flag | Status flag bits for the sea ice concentration, as described in Table |
Tb_corr | Corrected brightness temperatures [K] |
Tb | Uncorrected brightness temperatures [K] |
time | Time of data [year, month, day] |
xc | coordinate of projection [km] |
xy | coordinate of projection [km] |
lat | Latitude [°] |
lon | Longitude [°] |
Missing dates.
Year | Jan | Feb | Mar | Apr | May | Jun |
---|---|---|---|---|---|---|
– | 28 | 1–31 | 1–30 | 1–27 | – | |
12–15 | 10–11 | 24–31 | 1–5 | 14–15 | – | |
– | 21 | 24,31 | 1–8,16–30 | 1–2 | 3–30 | |
1, 3, 5, 7, | 2, 6, 12, | 1, 3, 5, 7, 9, | 2, 4, 6, 8, 10, | 2, 6, 8, | 13, 15, 17, | |
9, 13, 15, 17, | 14, 16, 18, | 11, 13, 15, 19, | 12, 14, 16, 18, 20, | 10, 12, 20, | 19, 21, 23, | |
19, 21, 23, 27 | 22, 24, 26, 28 | 21, 23, 25, 27, 31 | 22, 24, 28, 30 | 24, 30 | 125, 27, 29 | |
9–18, 23, | 4, 6–8, 10, | 4, 6, 8, 10, | 1–29 | 1, 3, 5, 7, 9, 11–16 | ||
25, 27, 29, 31 | 12, 14, 16, 18, | 12, 14, 16, 20, | ||||
20-22, 24, 26, 28 | 26, 28, 30-31 | |||||
Year | Jul | Aug | Sep | Oct | Nov | Dec |
31 | ||||||
27–28 | 1–31 | 1–4, 13–16, 22–23 | 14–15 | 9–10, 29–30 | 1, 5, 14-17 | |
17–19 | 1–8, 13–14 | – | 22–23, 25–27 | 1 | – | |
1–31 | 1-17, 19, 21-25, | 5-6, 8, 10, | 2, 4, 6-8, | 1, 3, 5, 7, 9, 11, | 1, 3, 5, 7, 9, 11, | |
27, 29, 31 | 12, 14, 16, 18, | 10, 12, 14, 18, | 13, 15, 17, 19, | 13, 15, 17, 19, 21, | ||
20, 22, 24, | 20, 22, 24, 26 | 21, 23, 25, 27, 29 | 23, 25, 27, 29, 31 | |||
26, 28, 30 | ||||||
1, 3, 5, 7, 9, | 2, 4, 6, 8, 10, | 1, 5, 11, | 1, 3, 5, 7, 9, | 2, 4, 6, 8, | 2, 4, 6, 8, | |
11, 13, 15, 17, 19, 21, | 12, 14, 16, 18, 20, | 13, 15, 17, | 13, 15, 17, 19, 21, | 10, 12, 14, 16, | 10, 12, 14-18, | |
23, 25, 27, 29, 31 | 22, 24, 26, 28, 30 | 23, 25, 27 | 23, 25, 27, 29, 31 | 18, 20, 22, 24, 26, 28, 30 | 20, 22, 24, 31 | |
26, 28, 30 | ||||||
Author contributions
WK performed the experiments, investigated the data, developed software, and wrote the manuscript, with contributions from all authors. RT contributed to the conceptualization, methodology, software, investigation, original draft, review, and editing of the manuscript and provided supervision throughout the work. JS contributed with critical feedback that shaped the interpretation and presentation of the data and improved the manuscript through reviewing and editing.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Acknowledgements
This work is part of the ESA Climate Change Initiative (ESA CCI) program's Sea Ice CCI (Sea_Ice_cci) project and the Danish National Centre for Climate Research (NCKF) at the Danish Meteorological Institute (DMI). The authors would like to thank Tadea Veng, Roberto Saldo, Thomas Lavergne, Atle Sørensen, and Leif Toudal Pedersen for their inputs and valuable feedback. We would also like to thank Imke Sievers, John Andrew Dawson, and Andreas Svejgaard Jensen for developing and testing the data filters described in Sect. 2.
Review statement
This paper was edited by Alexander Fraser and reviewed by two anonymous referees.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
The Electrically Scanning Microwave Radiometer (ESMR) instrument onboard the NIMBUS 5 satellite was a one-channel microwave radiometer that measured the 19.35 GHz horizontally polarized brightness temperature (
Here, we have reprocessed the entire dataset using a modern processing methodology that includes the implementation of pre-processing filtering, dynamical tie points, and a radiative transfer model (RTM) together with numerical weather prediction (NWP) for atmospheric correction. We present the one-channel sea ice concentration (SIC) algorithm and the model for computing temporally and spatially varying SIC uncertainty estimates. Post-processing steps include resampling to daily grids, land-spillover correction, the application of climatological masks, the setting of processing flags, and the estimation of sea ice extent, monthly means, and trends. This sea ice dataset derived from the NIMBUS 5 ESMR extends the sea ice record with an important reference from the mid-1970s. To make it easier to perform a consistent analysis of sea ice development over time, the same grid and land mask as used for EUMETSAT's OSI-SAF SMMR-based sea-ice climate data record (CDR) were used for our ESMR dataset. SIC uncertainties were included to further ease comparison to other datasets and time periods.
We find that our sea ice extent in the Arctic and Antarctic in the 1970s is generally higher than those available from the National Snow and Ice Data Center (NSIDC) Distributed Active Archive Center (DAAC), which were derived from the same ESMR dataset, with mean differences of 240 000 and 590 000 km
The new ESMR SIC dataset has been released as part of the ESA Climate Change Initiative (ESA CCI) program and is publicly available at 10.5285/34a15b96f1134d9e95b9e486d74e49cf
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details

1 National Space Institute, Technical University of Denmark (DTU Space), 2800 Lyngby, Denmark; Danish Meteorological Institute (DMI), National Centre for Climate Research (NCKF), Copenhagen, Denmark
2 National Space Institute, Technical University of Denmark (DTU Space), 2800 Lyngby, Denmark
3 Centre for Earth Observation Science (CEOS), University of Manitoba, Winnipeg, Canada; Department of Earth Sciences, University College London (UCL), London, UK; National Snow and Ice Data Center (NSIDC), University of Colorado, Boulder, Colorado, USA