Introduction
Sea-ice concentration has already been successfully used as a primary variable in certain sea-1ce models for input and verification (Reference Parkinson and WashingtonParkinson and Washington 1979, Reference Hibler and AckleyHibler and Ackley 1982), whereas sea-ice thickness is a primary variable in other models (e.g. Reference CoonCoon 1980) that have been limited by a lack of suitable data. The small but potentially-observable changes in ice-sheet elevation will also require models specifically designed to study these changes. Traditionally, after ice models are constructed, data requirements have been considered by remote-sensing specialists. However, the physics of the emission and propagation of electromagnetic radiation and the available technology combine to dictate the observable parameters with very little flexibility. Fortunately, these factors have united in recent years to begin to yield significant information about ice variability on a global scale, and from daily to decadal time scales. Because ice models are inherently more flexible than sources of long-term data, models need to be adapted or developed to complement the few parameters that can be obtained only from satellite remote sensing.
Determination of the variability of polar ice over climatically relevant time scales requires repetitive synoptic data on all high-latitude regions. Satellite data provide information on present climatic conditions in polar regions as well as information useful in studies of the time-dependent behavior of sea ice, ice sheets, ice shelves, and continental snow cover. Significant progress has been made in recent years on the observation of snow and ice during all weathers and in all seasons by passive microwave imaging, in particular. The seasonal cycle of sea-ice extent is now comprehensively described on a regional and hemispheric basis. Observations of ice-sheet surface climate parameters and changes in surface elevation have been initiated, as have satellite observations of changes in ice-shelf fronts. Recently, Reference Matson and WiesnetMatson and Wiesnet (1981) and Reference Dewey and HeimDewey and Heim (1982) described visual satellite observations of snow cover in the northern hemisphere since 1966. Advances have also been made in the microwave observation of the seasonal cycle of continental snow cover (Reference Kunz, Patil and RottKunzi and others 1982). This paper provides an overview of some recent results of observing the variability of polar ice, specifically sea ice, ice sheets, and ice shelves.
Sea-Ice Variability
The primary variation of sea ice is the seasonal cycle, and the interannual variations of the seasonal cycle are of particular climatic interest. The interannual variability of sea ice includes significant year-to-year variations, short-term variations (over several years), and possible long-term trends (extending over a decade and longer). Primary variables describing the sea-ice cover are geographical extent, concentration or actual ice area, and the amount of multiyear ice pack. Other parameters such as sea-ice thickness and drift velocities are important for studies of ice dynamics and ice processes, but here we concentrate on parameters that can be obtained over the required spatial and temporal scales for climate studies.
Sea-ice extent
Sea-ice extent, defined as the area of ocean covered by at least 10 to 1535 of sea ice, is the most widely used parameter describing sea-ice coverage. Before the advent of satellite imaging in the 1960s, information on ice extent was very limited, particularly in the Antarctic, and complete spatial and temporal coverage did not begin until the launch of the first passive microwave imager known as the electrically scanning microwave radiometer (ESMR) on the Nimbus-5 satellite in December 1972 (see, for example, Reference Gloersen, Wilheit, Chang, Nordberg and CampbellGloersen and others 1974, Reference ZwallyZwally and Gloersen 1977).
The longest data set on monthly sea-ice cover, with nearly-complete coverage of the Arctic from 1953 to 1977, was compiled by Reference Walsh and JohnsonWalsh and Johnson (1979) in order to examine the relationship between atmospheric variability and sea-ice variability; statistically significant correlations between the dominant modes in the atmosphere and in the sea ice were deter-mined. In the Antarctic, characteristics of the seasonal cycle, including regional contrasts and interannual variations, were examined in detail using a four-year (1973–76) ESMR data set (Reference Zwally, Comiso, Parkinson, Campbell, Carsey and GloersenZwally and others 1983[a]). The most notable regional difference in the Antarctic is the more rapid growth of the ice extent to a near maximum in the Ross Sea sector (160°E to 130°W) between March and July, in contrast to a slower growth to a maximum in the other regions in September. In all regions of the Antarctic, the seasonal decay is very rapid from October to January. The asymmetry in the Antarctic cycle, caused by the faster decay in spring and summer than the growth in fall and winter, is in contrast to a nearly symmetrical seasonal cycle in the Arctic around a March maximum and a September minimum (Reference Walsh and JohnsonWalsh and Johnson (1979) noted that the Arctic growth rate is slightly more rapid than the decay).
The seasonal cycle of sea ice in the Southern Ocean for ten years (1973–82), shown in Figure 1 which is adapted from Reference Zwally, Parkinson and ComisoZwally and others (1983[c]), exhibits an interesting and unexplained variability. Over periods of three to five years ice extents with greater winter maxima are associated with ice extents with lesser summer minima (Fig.1). However, over the full nine-year period, the trends exhibited by linear-fits to the data are nearly the same in each of the four seasons (Reference Zwally, Parkinson and ComisoZwally and others 1983[c]). Regional interannual variability is large, with sea ice decreasing in some regions while increasing in others. A significant decrease in overall annual mean ice extent on the Southern Ocean during the mid-1970s, suggested by Reference Kukla and GavinKukla and Gavin (1981) to be reflecting a C02-induced warming, has not been maintained. In particular, ice extent in the Weddell Sea sector (60°W to 20°E) has rebounded after a large decrease concurrent with a major oceanographic anomaly, i.e. the Weddell polynya. A decrease of similar magnitude in the Ross Sea sector followed the Meddel1 decrease by three years with a minimum in 1980. The mid-1970s decrease was also preceded by an increase in ice extent from 1966 to 1972, again indicating the presence of short-term components of variation that obscure any long-term trends that might be caused by a CO2 warming.
Modeling of the seasonal cycle over several years by Reference Parkinson and BindschadlerParkinson and Bindschadler (1984, in press) has provided estimates of the sensitivity of the amplitude of the Antarctic seasonal cycle to changes in atmospheric temperature. The aforementioned modulation of the seasonal cycle, shown by the association of greater winter maxima with lesser summer minima and conversely, presents a challenging problem for simulation by models. The amplitude modulations (of ~30%) occurred in several individual sectors, but with an apparent phase shift of about 2 to 3 a successively from the Weddell, to the Ross, and to the Bel1ingshausen-Amundsen seas sector (130°W to 60°W). Since it is unlikely that the cause of this modulation is limited to the sea-ice system, models that account for relevant oceanic and atmospheric interactions will probably be required to simulate and explain the observed interannual variations.
Sea-ice concentration
Sea-ice concentration, i.e., the percentage of a given ocean area covered by sea ice, has been derived from the ESMR data for the Southern Ocean to an estimated accuracy of 15% at each location (Reference Comiso and ZwallyComiso and Zwally 1982). The actual ice area is obtained by spatial integration of the concentration. The open water within the ice pack is the difference between the ice extent and the actual ice area, and the mean ice concentration is the actual ice area divided by the total ice extent. Any two of’these four parameters are sufficient to describe the area of ice coverage and its compactness.
In the Antarctic, the seasonal cycle of open water clearly differs from the seasonal cycles of ice extent and actual ice area, as shown in Figure 2.
Typically, the area of open water increases linearly for about nine months from February to November and then decreases sharply, although there is considerable monthly variation about the mean behavior. During the growth season, the total open-water area within the ice pack naturally increases as the increase in total ice extent expands the region in which open water can occur. The amount of open water is also related to the mean ice concentration, which reflects the rates of ice divergence and refreezing within the ice pack and the resulting areas of open leads and polynyas. The area of open water continues to increase for one to two months after the time of maximum ice extent, a feature which is also shown by the faster decrease in the actual ice area than in the ice extent just after their maxima (Fig.2). This continued increase of open water during the early part of the decay season is consistent with a rapid early decrease in areas of highly concentrated ice and a corresponding decrease in mean concentration (Reference Zwally, Comiso, Parkinson, Campbell, Carsey and GloersenZwally and others 1983[a]).
The mean ice concentration in the Southern Ocean ranges from approximately 60 to 80% in summer and winter, respectively. Although the summer ESMR values are in good agreement with the Navy-National Oceanic and Atmospheric Association (NOAA) sea-ice maps (Reference Godin and BarnettGodin and Barnett 1979) and reported, for example, by Reference Kukla and GavinKukla and Gavin (1981), the Navy-NOAA maps show less than 5% open water in winter in contrast to the 20% values derived from ESMR data. The ESMR-derived concentrations are considered to be more accurate (e.g. Reference Comiso and ZwallyComiso and Zwally 1982, and Reference Zwally, Comiso, Parkinson, Campbell, Carsey and GloersenZwally and others 1983[a]), because the microwave techniques account for all the open water in leads and polynyas, no matter how small the individual areas of open water are. From February to March, the Navy-NOAA estimates of mean ice concentration increase abruptly from 70 to approximately 95% and in April to 99%; these values would be typical of the highly constrained winter ice in the central Arctic Ocean. The ESMR-derived values show a much smoother increase in mean concentration throughout the growth season as the ice cover expands toward the equator. In the Southern Ocean, the unconstrained ice is much more divergent than in the Arctic Ocean, producing large areas of reduced ice concentration, particularly near the ice edge as shown by ESMR. Within the center of the ice pack, much of the 1ce cover has concentrations within the range of 88 to 100% (12% to zero open water). Near the Antarctic coast, reduced concentrations are also observed, due to the opening and closing of near-shore leads and polynyas.
Over the four years (1973–76), the trends in ice extent and actual ice area are very similar in the Antarctic, but differences of short duration are probably the result of differences in synoptic weather patterns. An example of different variability in ice extent, compared to actual ice area, is given by the winter ice covers of 1973 and 1974, for which the year-to-year decrease in actual ice area was greater than the decrease in ice extent (Fig.2). During July, September, and October 1973, less open water occurred within the ice pack than in 1974. The existence of a greater ice extent and a lesser amount of open water suggests that the ice pack was generally less divergent during the winter of 1973. Since there is the possibility of longer-term trends in the differences between ice extent and actual ice area due, for example, to variations in atmospheric circulation or storm frequency, longer-term observation and analysis of both actual ice area and ice extent are needed to describe and understand the observed changes in the sea-ice cover.
The sea-ice concentrations shown in Figure 3 are derived (Reference Cavalieri, Gloersen and CampbellCavalieri and others 1981 and in press) from the multifrequency and dual-polarization passive microwave imager termed the scanning multichannel microwave radiometer (SMMR) on Nimbus-7 which was launched in October 1978. During winter months, the amount of open water within the central Arctic 1ce pack may be only several percent, which is less than the amount generally Indicated by the concentrations in Figure 3 (upper). The variation of ice emissivity with ice type, a characteristic which enables the mapping of multiyear ice as discussed below, also introduces errors in the determination of ice concentration. Additional research is required on the methods for extraction of ice parameters from multi-frequency data over regions having small percentages of open water and mixtures of ice types.
Multiyear ice distribution
The multiyear ice fraction is defined as the percentage of the ice cover consisting of ice that has survived at least one summer’s melt season. This definition of multiyear ice is often used rather than the term “old ice”, which, in WM0 sea-1ce terminology (World Meteorological Organization 1970), is defined to include second-year ice and “multiyear” ice as sub-categories, with “multiyear” defined as ice that is older than second-year. Due to brine drainage processes in the freeboard, Arctic multiyear ice has a structure that produces a lower microwave emissivity than first-year ice (see, for example, Reference Wilheit, Nordberg, Blinn, Campbell and EdgertonWilheit and others 1972, Reference Gloersen, Zwally, Chang, Hall, Campbell and RamseierGloersen and others 1978). The lower microwave emissivity allows determination of the multiyear ice fraction (Fig.3 (lower)) from the multifrequency SMMR data (Reference Cavalieri, Gloersen and CampbellCavalieri and others 1981 and in press). However, in the Weddell Sea region of Antarctica, where extensive brine drainage does not occur (Reference Gow, Ackley, Weeks and GovonlGow and others 1982), the multiyear ice is radiometrically similar to first-year ice and is not distinguished from first-year ice by passive microwave techniques (Reference Zwally, Comiso, Parkinson, Campbell, Carsey and GloersenZwally and others 1983[a]).
As more is learned about the microwave properties of multiyear ice, in particular, various distinctions of ice types, or clarifications of terminology appropriate to remote sensing, may be needed. For example, as discussed in Reference Campbell, Gloersen, Zwally, Dyer and ChryssostomidisCampbell and others (1983), refrozen melt ponds appear to have microwave emissivities similar to first-year ice. Consequently, the area of refrozen melt ponds on multiyear ice floes in the Arctic is probably classified by the passive microwave techniques as first-year ice, even though the substrata of the melt ponds is multiyear ice. Therefore, the present determination of the multiyear ice fraction, as shown in Figure 3 (lower), would be the fraction of the ice cover consisting of multiyear ice floes multiplied by the fraction of the multiyear ice floes not covered by melt ponds having the emissivity of first-year ice. Confirmation is needed of this interpretation, which provides an explanation of the discrepancy between the low values (~ 40%) of the multiyear ice fraction in the region several hundred kilometers from the edge of the multiyear ice pack and other higher estimates (N Untersteiner private communication). The largest values of the multiyear fraction (>90%) indicated by the SMMR data are correctly located in the region of the known highest multiyear ice fraction, and lower coverage of melt ponds, north of the Canadian Arctic.
Interannual variability of the extent of the multiyear ice pack can be studied by observation of the extent of sea ice at the summer minimum (see, for example, Reference CarseyCarsey 1982), or by observation of the multiyear ice fraction throughout the winter months. Changes in the actual area of multiyear ice during winter, due to ridging and other processes, should be small, but the distribution of multiyear ice can be expected to change due to overall drift of the pack and continual movement within the pack. Studies of the seasonal and interannual changes in the extent and distribution of multiyear ice should be possible with sequential data such as those displayed in Figure 3 (lower).
North of Alaska, in February 1979, a band of predominantly first-year ice several hundred kilometers wide extended from the Beaufort Sea to the East Siberian Sea (Fig.3). The edge of the pack around the time of the previous summer minimum (23–26 September 1978, shown on Fig.3 (lower)), shows that this band of first-year ice was consistent with the unusually large extent of open water there in the summer. Therefore, the SMMR data do in fact outline the approximate position of the multiyear pack in winter. Differences between the boundary at the summer minimum and the boundary of the multiyear pack in February may be due in part to the influence of new ice growth on the ice edge near the summer minimum and in part to drift of the multiyear pack during winter.
In some regions during winter months the boundary of the multiyear ice pack can also be discerned in the single-wavelength ESMR data as a strong spatial gradient in the brightness temperature. In the Laptev and East Siberian seas during November this boundary is exhibited by a decrease in brightness temperature from 240 K, over areas covered totally by first-year ice, to 225 K over multiyear ice, in a distance of about 100 km from one to the other. Analysis of the 232 K contour, the approximate midpoint of this gradient, in the Laptev Sea around 130°E and 80°N, for example, shows that the contour lies inside the position of the edge at summer minimum. From December 1973 to April 1974, the 232 K contour moved about 200 km northward, consistent with a general movement of the ice pack with the transpolar drift stream. A similar displacement of the multiyear pack during the winter of 1978–79 is evident in Figure 3 (lower). Preliminary analysis of the ESMR data for 1973–76 indicates that the drift of the 232 K contour during winter tends to be less in the East Siberian Sea than in the Laptev Sea, and the drift is quite variable from year to year. These preliminary results illustrate the potential for studies of the extent, distribution, and dispersion of the multiyear ice pack in the Arctic, its seasonal cycle, and its interannual variability, similar to studies of the total ice cover.
Ice-Sheet and Ice-Shelf Variability
The surfaces of ice sheets and ice shelves are subject to the seasonal cycle of temperature, radiation, precipitation, and ablation, as well as longer-term variations of the atmosphere, ocean, and solar insolation. Most of the ice parameters of interest should be averaged over at least one year, consistent with the slow response of glaciers to changes in boundary conditions. One exception is the extent and duration of surface melting in summer, which requires daily observations for studies of interannual variations. Two important parameters which are best observed by satellite in the study of the long-term behavior of ice sheets and ice shelves are changes in surface elevation and in the position of ice-shelf fronts.
Surface elevation and ice margins
At present it is” not known whether the existing ice sheets of Greenland and Antarctica are growing or shrinking. Present estimates of mass balance, which are based on estimates of the mass input (accumulation rate) and mass output (ablation, bottom melting, and iceberg discharge), are probably only accurate to ±50%. A review of results from the few areas which have been studied in detail by surface measurements indicate that either thickening or thinning of the local ice is the rule rather than the exception and steady state appears to be rare (Reference Zwally, Thomas and BindschadlerZwally and others 1981).
Observation of land-based ice margins in isolation is known not to be a very useful indicator of ice-sheet behavior, because the response of the ice sheet margins to changes in other parts of the ice sheet is very slow. Therefore, on ice sheets, observation of the entire surface elevation is needed. In contrast, because ice shelves are floating, changes at the ice-shelf front are more rapidly coupled with other parts of the ice shelf. Compared to ice sheets, the response of ice shelves to climate change can be fast (Reference Thomas, Sanderson and RoseThomas and others 1979) with significant changes possible within a century. For several reasons, such as the controlling effect of ice shelves on ice-sheet discharge upstream, long-term observation of ice shelf behavior from space is especially important.
Satellite laser altimetry (Reference Zwally, Thomas and BindschadlerZwally and others 1981, Reference Bufton, Robinson, Femiano and FlatowBufton and others 1982) provides a potential method of determining net mass balance, and hence changes in ice volume, directly from repeated measurements of surface elevation. Mean elevation changes of about 10 cm over areas of 100 to 5 000 kmz could be detected, allowing studies of regional and overall mass balance. Satellite radar altimetry, while less accurate than laser altimetry, has already been used for mapping the surface topography of portions of the Greenland and Antarctic ice sheets (Reference Zwally, Bindschadler, Brenner, Martin and ThomasZwally and others 1983[b]). The surface topography of southern Greenland, uncorrected for slope-induced errors (Reference Brenner, Bindschadler, Thomas and ZwallyBrenner and others 1983), is shown in Figure 4 together with a preliminary delineation of the major drainage basins (R A Bindschadler private communication).
Ultimately, repetitive satellite measurement of surface elevations along the major flow lines will provide unique information for studies of the present state of ice flow. For example, by measuring changes in surface elevation it should be possible to deter-mine whether the central portions of the ice sheets are in balance with the margins. Other information obtainable from accurate measurements of surface elevation includes the thickness of Ice shelves and estimates of bottom melting rates.
Although changes in the position of ice-shelf fronts can be monitored by Landsat-type imagery with a spatial resolution of 30 to 100 m, the absolute accuracy of geographical location is not so well determined. Recently, Reference Brenner, Bindschadler, Thomas and ZwallyThomas and others (1983) have shown that satellite radar altimetry can be used to map the seaward ice margins to an absolute horizontal accuracy of 100 to 1 000 m (Fig.5). From continual observation by radar altimetry, the short-term variations of the ice front due to calving of major icebergs could be determined, and possibly secular changes in the position of the front also.
Surface climate parameters
Several parameters describing the surface climate of the ice sheets are potentially observable from space, such as the extent and duration of surface melting, accumulation rates (in the dry firn zones), and surface temperature. Some of the progress on observing these parameters is described briefly.
The brightness temperature TB of snow increases markedly with wetness (e.g. Reference Stiles and UlabyStiles and Ulaby 1980) to values as high as 260 K at the ESMR wavelength of 1.55 cm, corresponding to an emissivity of 0.95 at the melting point. The magnitude of the change in TB due to surface melting depends on the value of TB of the dry firn, which ranges form 145 to 220 K. In regions of significant summer melting, such as the Larsen Ice Shelf in Antarctica and the percolation zone of southwest Greenland, TB values for dry firn are about 160 K, corresponding to an emissivity of about 0.7 (e.g. Reference ZwallyZwally and Gloersen 1977). The TB of dry firn also varies with the seasonal temperature cycle, but with a smaller amplitude than that of the surface temperature because the emissivity is less than unity, and also the radiation emanates from below the surface where the seasonal temperature fluctuations are damped (Reference ZwallyZwally 1977, Reference Comiso and ZwallyComiso and others 1982). The rise in TB due to the onset of surface wetness is therefore relatively large and easy to detect. Examples of increases in TB from about 160 K for dry firn to 245 K for wet firn on ice shelves of the Antarctic peninsula are shown in Reference ZwallyZwally and Gloersen (1977). Somewhat smaller increases of TB with wetness on the percolation zone of Greenland are shown in Figure 6. Mapping of the extent and duration of summer melting will require systematic analysis of passive microwave data to delineate the times and regions where TB departs significantly from the smooth seasonal cycle typical of continuously dry firn.
The relationship between accumulation rates and microwave emissivity has been investigated (e.g. Reference Chang, Gloersen, Schmugge, Wilheit and ZwallyChang and others 1976, Reference ZwallyZwally 1977) showing that the microwave emission is sensitive to the grain size distribution with depth. An empirical relationship based on those studies was used by Reference Rotman, Fisher and StaelinRotman and others (1982) to make ice-sheet accumulation maps. Although the microwave methods are promising, additional research is needed because factors other than the accumulation rate (such as the temperature history of the firn) affect the grain-size distribution with depth.
An independent measurement of the physical temperature would greatly assist the mapping of accumulation rates from the microwave measurements, as well as being of direct interest to studies of ice-sheet climate. Thermal infrared measurements from a satellite have been used recently (Reference ComisoComiso 1983) for measurement of sea-ice surface temperatures, but substantial difficulties are encountered in eliminating the effects of clouds in the intervening atmosphere. Further developments with microwave and infrared methods should lead to useful mapping of accumulation rates and surface temperatures of the ice sheets.
Summary
Significant progress has been made in the development of methods for observing the variability of polar ice, especially the seasonal cycle of sea-ice extent. In particular, the interannual variability of the Antarctic sea ice has been shown to exhibit regional and seasonal characteristics over periods of three to five years that differ significantly when compared to decadal time scales. The extent, drift, and dispersion of the multiyear ice pack can be derived from new satellite data. Observation of ice-sheet elevations and the positions of ice-shelf fronts by satellite altimetry has been initiated and should lead to a determination of the long-term behavior of the ice.
The importance of continual synoptic observations of polar ice is illustrated by the large regional and interannual variability of sea Ice obtained from space observations, and by the characteristic regional behavior of ice sheets and ice shelves as suggested by field studies. In some cases, additional development 1n order to measure particular parameters is needed, and in others the principal need is for the systematic processing and analysis of existing satellite data. In general, however, the continued acquisition of satellite data on polar ice cannot be assumed. New sensors such as laser altimeters are needed to measure certain parameters, and for others continued operation of proven sensors, such as passive microwave imagers, is required on polar orbiting satellites.
This research was supported by NASA’s Oceans Program. I thank Dr Don Cavalieri for the SMMR image used in Figure 3.