1. Introduction
A problem in mapping the polar sea-ice concentrations in both hemispheres has been the false indication of sea ice over the open ocean and at the ice edge. These spurious sea-ice concentrations result from the presence of atmospheric water vapor, non-precipitating cloud liquid water, rain and sea-surface roughening by surface winds. While these effects are relatively minor at polar latitudes in winter, they result in serious weather contamination problems at all latitudes in summer (Reference Gloersen and CavalieriGloersen and Cavalieri, 1986; Reference CavalieriCavalieri and others, 1992).
The severity of weather contamination on Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave/Imager (SSM/I) sea-ice concentration in summer is illustrated in Figure 1 for both the Northern (18 August 1988) and Southern (5 February 1989) Hemispheres. Figure 1a shows large spurious sea-ice concentrations as high as 100% in the Pacific Ocean, 35% adjacent to the Kamchatka Peninsula and as high as 50% adjacent to Newfoundland. Over the Southern Ocean (Fig. 1b), spurious concentrations of up to 35% encircle Antarctica, with a maximum concentration of 80% north of the Weddell Sea. Although the following discussion focuses on the Arctic, the technique is equally applicable to the Antarctic.
Fig. 1. MM/I sea-ire concentration maps unfiltered for weather effects for (a) the Northern Hemisphere on 18 August 1987, and (b) the Southern Hemisphere on 5 February 1989.
The problem was addressed for sea-ice concentrations derived from the Nimbtts-7 Scanning Multichannel Microwave Radiometer (SMMR) data through the development of a weather filter (Reference Gloersen and CavalieriGloersen and Cavalieri, 1986). The filter is based on the polarization (PR) and the spectral gradient ratio (GR) distribution of ice-free and ice-covered seas. PR is defined in terms of the 18.0 GHz channels (PR 18) and GR is defined in terms of the 18.0 and 37.0 GHz channels (GR(37/18)). Gloersen and Cavalieri used a semi-empirical radiative-transfer model (RTM) to calculate PR(18) and GR(37/18) values for various atmospheric and sea-surface conditions. The (PR, GR) distributions obtained are illustrated for both the Arctic (Fig. 2a) and Antarctic (Fig. 2b) for typical wintertime conditions. Figure 2 also illustrates the sensitivity in (PR, GR) resulting from changes in sea-surface temperature (SST), near-surface wind (NSW), cloud liquid water (L) and atmospheric water vapor (W). Each plot in Figure 2 also shows the NASA team algorithm triangle defined by the vertices for open water (OW), first-year (FY) and multi-year (MY) sea ice.
Fig. 2. MMR spectral gradient ratio GR(37/l8) vs polarization PR (18) for (a) the Arctic for 3-7 February 1979, and (b) the Antarctic for 27-31 October 1986, showing the effects of ocean temperature, winds and atmospheric cloud liquid water and water vapor. (From Reference Gloersen, Campbell, Cavalieri, Comiso, Parkinson and ZwallyGloersen and others (1992).)
In Figure 2, each value calculated from the RTM is given in brackets using the convention: [SST, NSW, L, W]. The first point plotted corresponds to a specular ocean (no winds, clouds or atmospheric water vapor) at the freezing point (271 K) and is labeled [271, 0, 0, 0]. Next, the RTM was run, varying one parameter at a time. The SST changed from 271 K to 306 K [306. 0, 0, 0]; NSW changed to 36 ms−1; L to 0.08 cm; and W to 8 cm. These points span the region of the OW cluster and thus provide information on the various contributing conditions giving rise to the observed distribution in (PR, GR) space. By changing only SST, the point moves to the lower left of the main OW cluster, an area devoid of data, indicating that there are no warm SSTs with no wind and a dry atmosphere in the region mapped. However, the main cluster is bounded by the other three points, together with the specular point, which indicates that variations in wind, cloud and water vapor result in the observed scatter of points over OW areas.
The weather filter for the Nimbus-7 SMMR sea-ice algorithm consists of setting the sea-ice concentration to zero if GR(37/18) is greater than 0.08. While this threshold technique is effective in eliminating most of the spurious sea-ice concentrations, it also eliminates actual sea-ice concentrations less than 12% for first-year ice and 8% for multi-year ice (Reference Gloersen and CavalieriGloersen and Cavalieri, 1986). This is not considered a serious problem for most applications, since the ice edge has previously been defined as the 15% ice-concentration contour as determined from satellite radiometers.
Application of a GR (37/19) filter for use with sea-ice concentration maps derived with the SSM/I sensor is less successful. Figure 3a shows au SSM/I sea-ice concentration map for 23 June 1988, generated with the GR (37/19) threshold filler. Erroneous ice concentrations up to about 50% are found over open oceans at low latitudes. The difference in effectiveness of this filter results from the closer proximity of the 19.4 GHz SSM/I channels to the center of the 22.2 GHz atmospheric water-vapor line than that of the 18.0 GHz SMMR channels. As shown below in section 2, the proximity of the 19.4 GHz channel to the 22.2 GHz water-vapor line is used to provide an improved weather filter (Fig. 3b).
Fig. 3. SM/I sea-ice concentration map for 23 June 1988. (a) Using only the GR(37/19) filter; (b) Using both the GR(37/19) and the GR(22/19) filers.
A plot of the atmospheric water-vapor attenuation coefficient for a frequency range of 10 30 GHz is show in Figure 4. The attenuation coefficient was calculated using the formulation given by Reference Ulaby, Moure and FungUlaby and others (1981) and, assuming standard sea-level pressure, a surface temperature of 300 K, and a water-vapor density of 1 gm−3. Examination of Figure 4 shows that 18.0 GHz is about one-third up the wing of the water-vapor line, whereas 19.34 GHz is about half-way up the line. Thus, the SSM/I 19.35 GHz channels are more sensitive to changes in atmospheric water vapor. This sensitivity results in greater contamination problems.
Fig. 4. Atmospheric water-vapor absorption coefficient as a function of frequency for standard sea-level pressure, a surface temperature of 300 K and a water-vapor density of 1 g m−3, following the formulation of Reference Ulaby, Moure and FungUlaby and others (1981).
2. Physical Basis of the SSM/I Filter
In this section, we briefly review the microwave polarization and spectral dependency of emission by ocean and atmosphere for the purpose of describing the physical basis of the SSM/I weather filter. The general theory of microwave emission from the ocean, atmosphere and ice-covered surfaces has been reviewed elsewhere (Reference Gloersen and BarathGloersen and Barath, 1977; Reference SwiftSwift, 1980; Reference Wilheit and ChangWilheit and Chang, 1980; Reference Gloersen, Campbell, Cavalieri, Comiso, Parkinson and ZwallyGloersen and others, 1992; St. Reference GermainGermain, 1993).
The radiation received by the Satellite and expressed as a brightness temperature (T B) consists of several components
where T B,λp is measured brightness temperature as a function of wavelength (λ) and polarization (p), ∈λP is total surface emissivity, T s is physical temperature of the surface, k λ is atmospheric opacity, T B,λ,up is total atmospheric upwelling radiation, T B,λ,down is total atmospheric downwelling radiation and T B,c is cosmic background radiation (2.7 K.)
The microwave emission (∈λ,p) from the ocean surface is a function of wavelength, polarization, angle of incidence, surface temperature and salinity. The surface emissivity will increase with increasing wavelength but will decrease with increasing surface temperature. Sea-surface salinity variability affects the emissivity only at wavelengths longer than about 6 cm (Reference Wilheit and ChangWilheit and Chang, 1980) and so is not of concern here, because all SSM/I wavelengths are 1.55 cm or shorter. The change in emissivity with angle of incidence depends on the polarization. Emissivity increases with angle of incidence at vertical polarization and decreases with increasing angle of incidence at horizontal polarization. Because the SSM/I is a conically scanning radiometer, its angle of incidence is constant at 53°.
In addition to an increase in sea-surface temperature, the presence of a surface wind will also increase the surface microwave emission through the roughening of the sea surface. The amount of increase again depends on wavelength and angle of incidence. For example, at a wavelength of 1.55 cm, the radiometric brightness increases by 1.2 K at 0° angle of incidence to about 1.8 K at 70° for each m s−1 increase in wind speed (Reference Nordberg, Conway, Ross and WilheitNordberg and others, 1971).
Atmospheric emission will also contribute to the radiometric brightness temperature measured by the satellite through the atmospheric opacity term (k) in Equation (1). The opacity depends on the total integrated amount of atmospheric oxygen, water vapor, non-precipitating cloud liquid water and rain through the entire atmosphere along the line of sight of the radio-meter. The contribution from oxygen is a constant for a given wavelength as is the cosmic background radiation term (T B,c) in Equation (1). The atmospheric variables of interest then are the water vapor (V), non-precipitating cloud liquid water (L) and rain rate (R).
The total atmospheric opacity (K) may be approximated by
For oxygen, water vapor and cloud liquid water, we use the following linear expressions following St. Reference GermainGermain (1993):
where A and B are wavelength- and temperature-dependent coefficients linearly, relating the total integrated quantity to opacity. C is the wavelength-dependent coefficient for oxygen. The expressions for these coefficients have been derived by St. Reference GermainGermain (1993) and are given here for each of the three SSM/I frequencies of interest:
These expressions are used in the RTM for obtaining the brightness temperatures for various atmospheric conditions.
An approximate expression for the dependency of atmospheric opacity on rain rate has been given by Reference Gloersen and BarathGloersen and Barath (1977):
for k < 0.4 where a = 3.51 × 10−2 + 5.55 × 10−2λ− 6.42 × 10−3λ2 b = 5.14 × 10−2λ−1.85, R is rain rate in mm h−1 and h is height of the rain column (km). The rainfall models used to derive these expressions assume a Marshall-Palmer distribution of rain droplets and utilize a calculation which includes Mie scattering and absorption Reference Gloersen and BarathGloersen and Barath, 1977).
Of direct relevance to the problem al hand is the variation of atmospheric water vapor, precipitating and non-precipitating liquid water, near-surface winds and rainfall rates on PR and GR, the two independent variables used in the SSM/I sea-ice algorithm. PR(19)s and GR(37/19)s were generated for a wide range of near-surface winds, cloud liquid water and atmospheric water vapor using the radiative-transfer model described by St. Reference GermainGermain (1993). Similarly, (PR, GR) pairs were generated for a range of rain rates using the Reference Wilheit and ChangWilheit and Chang (1980) model.
The variation of PR and GR with each of these atmospheric parameters is shown in Figure 5. The RTM was run using ocean-surface temperatures of 299 K (Fig. 5a) and 271 K (Fig. 5b) representative of mid and polar latitudes, respectively. As in Figure 2 for the SMMR, the SSM/I sea-ice algorithm is represented in (PR, GR) space by a triangle with curvilinear sides. The GR(37/19) threshold for the SSM/I is 0.05. For a sea-surface temperature of 299 K, we see thai this GR (37/19) threshold value eliminates cloud liquid-water amounts of up to 0.27 cm, rain rates of up to 15 mm h−1, water vapor of up to almost 2 cm and wind speeds of up to 20 m s−1) (Fig. 5a). For an ocean-surface temperature of 271 K, typical near the ice edge, the GR(37/19) threshold eliminates cloud liquid-water amounts of up to 0.20 cm, rain rates of up to 25 mm h−1, water vapor of less than 4 cm and wind speeds of up to 30 m s−1 (Fig. 5b). Since these atmospheric parameter values are seldom realized at polar latitudes, the GR(37/19) filter generally works well in most cases. It is at lower latitudes or in summer with water-vapor values greater than 2 cm that the GR(37/19) filter fails.
Fig. 5. PR(19) vs GR(37/19) showing the distribution of surface wind speed (0-50 m s−1), cloud liquid water (0-0.3 cm) and atmospheric water vapor (0-6 cm) and rain rate (0-50 mm h−1), calculated using radiative-transfer models with surface temperatures of (a) 299 k and (b) 271 K. The CR(37\19) threshold value (0.05) was used in the first SSM/I weather filter and is indicated by the horizontal solid line.
At this threshold level, ice concentrations tip to 15% are also eliminated. By adjusting the threshold to a lower GR value, more of the water-vapor effects would be eliminated but at the expense of eliminating ice concentrations greater than 15%. The problem, then, is to eliminate the Water-vapor effects, while minimizing the loss of ice-concentration information.
As Figure 4 shows, the proximity of the 19.4 GHz channel to the 22.2 GHz atmospheric water-vapor line is the reason for the greater SSM/I sensitivity to water-vapor variations. However, this sensitivity can be turned to an advantage through the definition of an atmospheric water-vapor filter. The difference between the 22.2 and 19.4 GHz vertical polarization channels is highly sensitive to small increases in atmospheric water vapor in the range 0-4 cm and is used as the basis of the additional weather filler needed to eliminate the water-vapor contamination illustrated in Figure 3a. This difference is normalized by the sum of the brightness temperatures, thereby defining a new spectral gradient ratio, GR(22/19):
The GR approach was chosen to help minimize the effect of ice-temperature variations in cases where we have a field of view with both ice-free and ice-covered seas. The rationale for this approach has been discussed in detail by Reference Gloersen and CavalieriGloersen and Cavalieri (1986). A threshold value of 0.045 was selected for the GR(22/19) filter as a trade-off between eliminating spurious sea-ice concentrations and losing low ice-concentration information at the ice edge. The effectiveness of this additional filter and the loss of sea-ice information is discussed next.
3. Discussion of Results
A plot of GR(37/19) versus GR(22/10) illustrates the potential effectiveness of each filter in eliminating the effects of each of the atmospheric components (Fig. 6). The chosen threshold values for each GR are also indicated in the figure by horizontal and vertical lines. As was noted in the discussion of Figure 5a, a threshold value of 0.05 for GR(37/19) eliminates cloud liquid-water amounts of up to 0.27 cm, rain rates of up to 15 mm h−1, water vapor of less than about 2 cm and wind speeds of up to 20 m s−1 from contributing to spurious sea-ice concentrations at low latitudes. From Figure 6, we see that the GR (22/19) threshold of 0.045 is effective only for water vapor. All other atmospheric variables fall below the 0.045 threshold. Thus, both GR filters are required.
Fig. 6. GR(22/19) vs GR(37/19) showing the distribution of surface wind speed (0-25 m s−1), cloud liquid water (0-0.3 cm) and atmospheric water vapor (0-5 cm) and rain rate (0-100 mm h−1) calculated using radiative-transfer models. The GR(22/19) and GR(37\19) threshold values used in the SSM/I weather filler are indicated by the horizontal and vertical solid lines.
The model results shown in Figure 6 also help determine the conditions under which both filters will fail. The horizontal and vertical lines defined by the two threshold values define an area in the lower left part of the figure. Values of cloud liquid water, water vapor, wind speed and rain rates within this area will still result in spurious sea-ice concentrations. Extrapolation of these atmospheric parameters toward higher values tends toward the origin of the plot (PR = 0, GR = 0), a condition corresponding to a perfect black-body emitter of radiation. Some examples of these extreme conditions are given below.
The sea-ice concentration map for 23 June 1988 is shown in Figure 3b after the application of the additional GR (122/19) filter. All of the spurious ice concentrations are filtered, except for two small areas of high concentration off the coast of Nova Scotia. These, as well as other examples of areas of spurious high ice concentrations, are not of particular concern, because they generally occur at latitudes where sea ice is known not to be present. Of much greater concern is the presence of spurious concentrations at or near the ice edge.
An ice-concentration difference map for 23 June shows the spurious ice concentrations that are filtered (Fig. 7). Except for two small areas, one in the Beaufort Sea near Banks Island and the other in the Barents Sea between Novaya Zemlya and Svalbard, the sea-ice edge is preserved and, as just noted, almost all of the spurious concentrations are filtered. The concern is that some actual ice-concentration information is lost in these two areas. Cloud cover precluded the use of AYHRR imagery to verify the location of the ice edge.
As an alternative test, the sea-ice concentration variability along two transects cutting across each of these two areas was examined to help determine which of the concentrations may be real and which concentrations result from passing weather systems. Presumably weather-related, spurious concentrations will show greater day-to-day variability. The first transect examined cuts across the ice from a point in the open-water area near Banks Island in the Beaufort Sea to a point well within the ice pack north of the area in question (see Fig. 7). The sea-ice concentrations are plotted in Figure 8 for a 4 d period from 21 to 24 June 1988. The four transects are plotted first from ice-concentration maps generated without the use of either weather filter (Fig. 8a and c) and then with both filters (Fig. 8b and d).
Fig. 8. ea-ice concentrations from the two transects shown in Figure 7 for a sequence of 4 d in 1988. The Beaufort Sea ice concentrations calculated (a) without and (b) with weather filters. The Barents, Sea ice concentrations calculated (c) without and (d) with weather fitters.
Without the use of weather filters and for all 4 days, ice concentrations range from 4% to 16% at pixel positions 0-4 (see Fig. 8a) and sharply increase at pixel position 5 until pixel 13 is reached. The ice concentrations then level off and vary from about 90% to 100% to the end of the transect (pixel position 21). Comparison of Figure 8a and b shows that for all 4 days, only those concentrations at pixel positions 0-4 were eliminated by the filters. Since the ice concentrations (90-100%) within the ice pack varied as much as those over the open-water area, one should consider the possibility that weather effects may contribute to some of this variability. A study by Reference MaslanikMaslanik (1992) has suggested that water vapor and cloud liquid water cause overestimates of total ice concentration over areas of high ice concentration, but the effects found were smaller than the observed variations here. The days with the highest (day 174) and lowest (day 172) concentrations over the polynya area (pixels 0-4) have comparable concentrations within the pack ice (Fig. 8a). Thus, weather effects over the pack ice are apparently negligible.
The transect for the Barents Sea region (Fig. 8c and d) shows much greater ice-concentration variability starting from between pixel positions 45 and 50. This also corresponds to a sharp break in the rapid decrease in ice concentration (presumably the actual ice edge). Application of both weather filters in this case results in a sharp ice edge with zero concentrations at pixel positions ranging from 44 to 48.
4. Summary
The Nimbus-7 SMMR technique for filtering spurious sea-ice concentrations over open ocean and near the ice edge is less successful with the DMSP SSM/I because of the closer proximity of the 19.4GHz SSM/I channels to the 22.2GHz water-vapor line than that of the 18.0 GHz SMMR channels. An additional filter, based on a combination of the 19.4 and the 22.2 GHz SSM/I channels has been developed to address this problem. The new filter effectively reduces the spurious ice concentrations, especially in the higher latitudes, where such ambiguity poses the greatest problem.
Analysis based on an RTM indicates that the performance of the filler is optimal for winds less than 30 m s−1, cloud liquid water less than 0.2 m, water vapor less than 24 cm and rain rates less than 12 mm h−1. Weather exceeding these conditions is uncommon in high latitudes, which supports the observed effectiveness of the filter.
While the technique works well in reducing weather effects in most cases, the trade-off in terms of actual ice-concentration information has not been fully quantified. Our analysis indicates that the actual ice-edge location (defined as 15% ice concentration) is not affened by the filter. Further validation studies using SAR, AVHRR or OLS data are required to determine conclusively the effect of the filter on ice-edge-location accuracy and the extent to which actual ice information may be lost.
Acknowledgements
The authors gratefully acknowledge many helpful discussions with A. Chang and P. Gloersen during the course of this work. A. Chang kindly provided the brightness temperatures from his rain-rate model used in this paper. We also thank M. Goodberlet for providing the open-ocean wind-speed model. M. Martillo and J. Sims of Hughes STX Corp. provided the programming and graphics support. D.J. C. was supported partly by the NASA Cryospheric Processes Program (RTOP 971-578-32-20) and partly by the NASA EOS Project 22904-15-05).
