Search

GRACE and ICESat Antarctic mass-balance differences are resolved utilizing their dependencies on corrections for changes in mass and volume of the same underlying mantle material forced by ice-loading changes. Modeled gravimetry corrections are 5.22 times altimetry corrections over East Antarctica (EA) and 4.51 times over West Antarctica (WA), with inferred mantle densities 4.75 and 4.11 g cm−3. Derived sensitivities (Sg,Sa) to bedrock motion enable calculation of motion (δB0) needed to equalize GRACE and ICESat mass changes during 2003–08. For EA,δB0is −2.2 mm a−1subsidence with mass matching at 150 Gt a−1, inland WA is −3.5 mm a−1at 66 Gt a−1, and coastal WA is only −0.35 mm a−1at −95 Gt a−1. WA subsidence is attributed to low mantle viscosity with faster responses to post-LGM deglaciation and to ice growth during Holocene grounding-line readvance. EA subsidence is attributed to Holocene dynamic thickening. With Antarctic Peninsula loss of −26 Gt a−1, the Antarctic total gain is 95 ± 25 Gt a−1during 2003–08, compared to 144 ± 61 Gt a−1from ERS1/2 during 1992–2001. Beginning in 2009, large increases in coastal WA dynamic losses overcame long-term EA and inland WA gains bringing Antarctica close to balance at −12 ± 64 Gt a−1by 2012–16.

Conversion of altimetry-derived ice-sheet volume change to mass requires an understanding of the evolution of the combined ice and air content within the firn column. In the absence of suitable techniques to observe the changes to the firn column across the entirety of an ice sheet, the firn column processes are typically modelled. Here, we present new 40-year simulations of firn processes over the Greenland and Antarctic Ice Sheets using the Community Firn Model and atmospheric reanalysis variables. A dataset of more than 250 measured depth-density profiles from both ice sheets provides the basis of the calibration of the dry-snow densification scheme. The resulting scheme results in a reduction in the rate of densification, relative to a commonly used semi-empirical model, through a decreased dependence on the accumulation rate, a proxy for overburden stress. The modelled firn column runoff, when combined with atmospheric variables from MERRA-2, generates realistic mean integrated surface mass balance values for the Greenland (+361 Gt yr−1) and Antarctic (+2623 Gt yr−1) ice sheets when compared to published model-ensemble means. We find that seasonal volume changes associated with firn air content are approximately 3 times larger than those associated with surface mass balance; however, when averaged over multiple years, ice and air-volume fluctuations within the firn column are of comparable magnitudes. Between 1996 and 2019, the Greenland Ice Sheet lost more than 5 % of its firn air content indicating a reduction in the total meltwater retention capability. Nearly all (>98 %) of the meltwater produced over the Antarctic Ice Sheet is retained within the firn column through infiltration and refreezing.

Mass changes of the Antarctic ice sheet impact sea-level rise as climate changes, but recent rates have been uncertain. Ice, Cloud and land Elevation Satellite (ICESat) data (2003–08) show mass gains from snow accumulation exceeded discharge losses by 82 ± 25 Gt a−1, reducing global sea-level rise by 0.23 mm a−1. European Remote-sensing Satellite (ERS) data (1992–2001) give a similar gain of 112 61 Gt a−1. Gains of 136 Gt a−1 in East Antarctica (EA) and 72 Gt a−1 in four drainage systems (WA2) in West Antarctic (WA) exceed losses of 97 Gt a−1 from three coastal drainage systems (WA1) and 29 Gt a−1 from the Antarctic Peninsula (AP). EA dynamic thickening of 147 Gt a−1 is a continuing response to increased accumulation (>50%) since the early Holocene. Recent accumulation loss of 11 Gt a−1 in EA indicates thickening is not from contemporaneous snowfall increases. Similarly, the WA2 gain is mainly (60 Gt a−1) dynamic thickening. In WA1 and the AP, increased losses of 66 ± 16 Gt a−1 from increased dynamic thinning from accelerating glaciers are 50% offset by greater WA snowfall. The decadal increase in dynamic thinning in WA1 and the AP is approximately one-third of the long-term dynamic thickening in EA and WA2, which should buffer additional dynamic thinning for decades.

Variations in accumulation rate As(t) and temperature Ts(t) at the surface of firn cause changes in the rate of firn compaction (FC) and surface height H(t) that do not involve changes in mass, and therefore need to be accounted for in deriving mass changes from measured H(t). As the effects of changes in As(t) and Ts(t) propagate into the firn, the FC rate is affected with a highly variable and complex response time. The H(t) during measurement periods depend on the history of As(t) and Ts(t) prior to the measurements. Consequently, knowledge of firn response times to climate perturbations is important to estimate the required length of the time series of As(t) and Ts(t) used in FC models. We use our numerical FC model, which is time-dependent on both temperature and accumulation rate, to examine the response times of both H(t) and the rates of change dH(t)/dt to variations in As(t) and Ts(t) using sample perturbations and climate data for selected sites in Antarctica. Our results show that the response times for dH(t)/dt, which are of particular interest, are much shorter than the responses of H(t). Typical response times of dH(t)/dt are from several years to Ts(t) from satellite observations since 1982 and As(t) from meteorological data since 1979 are essentially of sufficient length to correct for FC height changes for measurements beginning in 1992.

[1] Sea ice thicknesses derived from NASA's Ice, Cloud, and Land Elevation Satellite (ICESat) altimetry data are examined using two different approaches, buoyancy and empirical equations, and at two spatial scales—ICESat footprint size (70 m diameter spot) and Advanced Microwave Scanning Radiometer (AMSR-E) pixel size (12.5 km by 12.5 km) for the Bellingshausen and Amundsen Seas of west Antarctica. Ice thickness from the empirical equation shows reasonable spatial and temporal distribution of ice thickness from 2003 to 2009. Ice thickness from the buoyancy equation, however, additionally needing snow depth information derived from the AMSR-E, shows an overestimation in terms of maximum, mean (+63% to 75%), and standard deviation while underestimation in modal thickness (−20%) as compared with those from the empirical equation approach. When ICESat snow freeboard is used as the snow depth in the buoyancy equation, i.e., the zero ice freeboard assumption, the derived ice thicknesses match well with those from the empirical equation approach, within 5% overall. The AMSR-E, therefore, may underestimate snow depth and accounts for ~95% of the ice thickness overestimation as compared with the buoyancy approach. The empirical equation derived ice thickness shows a consistent asymmetrical distribution with a long tail to high values, and seasonal median values ranging from 0.8 to 1.4 m over the 2003–2009 period that are always larger than the corresponding modal values (0.6–1.1 m) and lower than the mean values (1.0–1.6 m), with standard deviation of 0.6–1.0 m. An overall increasing trend of 0.03 m/year of mean ice thickness is found from 2003 to 2009, although statistically insignificant (p = 0.11) at the 95% confidence level. Starting from autumn, a general picture of seasonal mean, modal, and median ice thickness increases progressively from autumn to spring and decreases from spring to the following autumn, when new thin ice dominates the ice thickness distribution. The asymmetric shape of the thickness distribution reflects the key role of ice deformation processes in the evolution of the thickness distribution. The statistical properties of the thickness distribution interannually (high range of mean thickness and standard deviation) indicate the variability of deformation processes. However, spring ice volume, the product of ice mean thickness and areal extent computed for the spring maximum, shows variability year to year but is primarily dominated by ice extent variability, with no increasing or decreasing trend over this record length. The dependence of the volume on the ice extent primarily suggests that ice thickness changes have also not covaried with the ice extent losses seen over the satellite record in this region, unlike the Arctic. These properties reflect the interactive processes of ice advection, thermodynamic growth and ice deformation that all substantially influence ice mass balance in the Bellingshausen-Amundsen Seas region.

We calculated Antarctic sea-ice freeboard and thickness for ICESat and ATM campaigns. A Gridded ATM freeboard map of nine IceBridge campaigns in October 2009, 2010, and 2011 is shown in Figure 1.

Ice velocities observed in 2005/06 at three GPS stations along the Sermeq Avannarleq flowline, West Greenland, are used to characterize an observed annual velocity cycle. We attempt to reproduce this annual ice velocity cycle using a 1-D ice-flow model with longitudinal stresses coupled to a 1-D hydrology model that governs an empirical basal sliding rule. Seasonal basal sliding velocity is parameterized as a perturbation of prescribed winter sliding velocity that is proportional to the rate of change of glacier water storage. The coupled model reproduces the broad features of the annual basal sliding cycle observed along this flowline, namely a summer speed-up event followed by a fall slowdown event. We also evaluate the hypothesis that the observed annual velocity cycle is due to the annual calving cycle at the terminus. We demonstrate that the ice acceleration due to a catastrophic calving event takes an order of magnitude longer to reach CU/ETH (‘Swiss’) Camp (46 km upstream of the terminus) than is observed. The seasonal acceleration observed at Swiss Camp is therefore unlikely to be the result of velocity perturbations propagated upstream via longitudinal coupling. Instead we interpret this velocity cycle to reflect the local history of glacier water balance.

Mass-balance analysis of the Greenland ice sheet based on surface elevation changes observed by the European Remote-sensing Satellite (ERS) (1992-2002) and Ice, Cloud and land Elevation Satellite (ICESat) (2003-07) indicates that the strongly increased mass loss at lower elevations (

Changes in ice-sheet surface elevation are caused by a combination of ice-dynamic imbalance, ablation, temporal variations in accumulation rate, firn compaction and underlying bedrock motion. Thus, deriving the rate of ice-sheet mass change from measured surface elevation change requires information on the rate of firn compaction and bedrock motion, which do not involve changes in mass, and requires an appropriate firn density to associate with elevation changes induced by recent accumulation rate variability. We use a 25 year record of surface temperature and a parameterization for accumulation change as a function of temperature to drive a firn compaction model. We apply this formulation to ICESat measurements of surface elevation change at three locations on the Greenland ice sheet in order to separate the accumulation-driven changes from the ice-dynamic/ablation-driven changes, and thus to derive the corresponding mass change. Our calculated densities for the accumulation-driven changes range from 410 to 610 kgm–3, which along with 900 kgm–3 for the dynamic/ablation-driven changes gives average densities ranging from 680 to 790 kgm–3. We show that using an average (or ‘effective’) density to convert elevation change to mass change is not valid where the accumulation and the dynamic elevation changes are of opposite sign.

Sea-ice freeboard heights for 17 ICESat campaign periods from 2003 to 2009 are derived from ICESat data. Freeboard is combined with snow depth from Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) data and nominal densities of snow, water and sea ice, to estimate sea-ice thickness. Sea-ice freeboard and thickness distributions show clear seasonal variations that reflect the yearly cycle of growth and decay of the Weddell Sea (Antarctica) pack ice. During October–November, sea ice grows to its seasonal maximum both in area and thickness; the mean freeboards are 0.33–0.41m and the mean thicknesses are 2.10–2.59 m. During February–March, thinner sea ice melts away and the sea-ice pack is mainly distributed in the west Weddell Sea; the mean freeboards are 0.35–0.46m and the mean thicknesses are 1.48–1.94 m. During May–June, the mean freeboards and thicknesses are 0.26–0.29m and 1.32–1.37 m, respectively. the 6 year trends in sea-ice extent and volume are (0.023±0.051)×106 km2 a–1 (0.45% a–1) and (0.007±0.092)×103 km3 a–1 (0.08% a–1); however, the large standard deviations indicate that these positive trends are not statistically significant.

We derive mass changes of the Greenland ice sheet (GIS) for 2003–07 from ICESat laser altimetry and compare them with results for 1992–2002 from ERS radar and airborne laser altimetry. The GIS continued to grow inland and thin at the margins during 2003–07, but surface melting and accelerated flow significantly increased the marginal thinning compared with the 1990s. The net balance changed from a small loss of 7 ± 3 Gt a−1 in the 1990s to 171 ± 4 Gt a−1 for 2003–07, contributing 0.5 mm a−1 to recent global sea-level rise. We divide the derived mass changes into two components: (1) from changes in melting and ice dynamics and (2) from changes in precipitation and accumulation rate. We use our firn compaction model to calculate the elevation changes driven by changes in both temperature and accumulation rate and to calculate the appropriate density to convert the accumulation-driven changes to mass changes. Increased losses from melting and ice dynamics (17–206 Gt a−1) are over seven times larger than increased gains from precipitation (10–35 Gt a−1) during a warming period of ∼2 K (10 a)−1 over the GIS. Above 2000 m elevation, the rate of gain decreased from 44 to 28 Gt a−1, while below 2000 m the rate of loss increased from 51 to 198 Gt a−1. Enhanced thinning below the equilibrium line on outlet glaciers indicates that increased melting has a significant impact on outlet glaciers, as well as accelerating ice flow. Increased thinning at higher elevations appears to be induced by dynamic coupling to thinning at the margins on decadal timescales.

The Geoscience Laser Altimeter System (GLAS) aboard ICESat, launched in January 2003, has been designed to detect and monitor changes in the cryosphere. The first objective of this paper is to present high-resolution ice-surface elevation maps derived from GLAS data, using geostatistical analysis. In a regional study of Walgreen Coast and Northern Ellsworth Land, West Antarctica, differences in the representation of geographic and morphologic features in maps based on ERS-1 radar altimeter data and on GLAS data are investigated, with the result that in particular in topographically complex coastal areas and the margin of the ice sheet the improvement in precision and accuracy of the laser altimeter is significant. A second, applied objective is to map elevation changes in Pine Island Glacier, a glacier that plays a key role in the question of stability of the West Antarctic Ice Sheet and has been changing rapidly in recent years. Results of elevation differencing of 2003-GLAS-data and 1995-ERS-1-radar-altimeter-data DEMs (1) show that thinning rates have been increasing and (2) are applied to attribute the observed changes in Pine Island Glacier to internal processes in the glacier, related to dynamic thinning. More generally, this application serves to demonstrate that GLAS data facilitate study of cryospheric change.

Sea ice surviving the summer melt season to become multi-year ice in the Arctic Ocean is of interest because multi-year ice significantly affects the ice-thickness distribution and the dynamics and thermodynamics of the ice pack in subsequent seasons. However, the amount of ice surviving summer melting has not been well determined because the time of the minimum ice area varies from region to region. A concept of local temporal minimum (LTM) accounts for non-simultaneity of the melt–freeze transition by determining the minima ice concentrations (CLTM) on local spatial scales. CLTM are calculated for 25 km gridcells using 24 years (1979–2002) of satellite passive-microwave data. The total area of ice surviving the summer melt (ALTM) is given by spatial integration of CLTM. Over 24 years, the average ALTM is 2.6 × 106 km2 (excluding ∼0.7 × 105 km2 above 84° N). In contrast, the average area (3.8 × 106 km2) of all ice types (ASM), measured when the total (simultaneous) ice cover is a minimum in daily maps in mid-September, is an often-used estimate of ice surviving the summer melting that is ∼45% too large. Over 24 years, the ALTM decreased by 9.5 ± 2.2% (10 a)−1 (0.27 ± 0.06 × 106 km2 (10 a)−1), which is similar to the rate of decline of ASM and about three times the rate of the annual average. The time-of-occurrence of the LTM averaged over the perennial ice pack increased by 8 days from around 11 to 19 August, indicating a later ending of the melt season by about 3 days (10 a)−1 as the summer pack declines. Estimates of multi-year ice in midwinter from passive microwave observations are ∼17% smaller than ALTM, suggesting that the microwave algorithm does not measure all the multi-year ice.

Changes in the surface elevation of the Greenland and Antarctic ice sheets and ice shelves caused by variations in the rate of firn compaction are calculated with a time-dependent firn densification model driven by two decades (1982–2003) of satellite-observed monthly surface temperatures. The model includes the effects of melting and refreezing, both the direct changes in density and the subsequent effects on the densification rate. As previously shown, the temperature-dependent rate of densification is largest in summer, but changes in winter temperatures also have a significant effect. Over the last decade, climate warming has enhanced the rate of compaction and lowered the average surface elevation of Greenland by 1.8 cma-1 and most of West Antarctica by 1.9 cma–1. In East Antarctica, a small cooling raised the average surface elevation by 0.14 cma–1.

Changes in ice mass are estimated from elevation changes derived from 10.5 years (Greenland) and 9 years (Antarctica) of satellite radar altimetry data from the European Remote-sensing Satellites ERS-1 and -2. For the first time, the dH/dt values are adjusted for changes in surface elevation resulting from temperature-driven variations in the rate of firn compaction. The Greenland ice sheet is thinning at the margins (–42 ± 2Gta¯1 below the equilibrium-line altitude (ELA)) and growing inland (+53 ± 2Gta-1 above the ELA) with a small overall mass gain (+11 ± 3Gta–1; –0.03 mma–1 SLE (sea-level equivalent)). The ice sheet in West Antarctica (WA) is losing mass (–47 ± 4Gta–1) and the ice sheet in East Antarctica (EA) shows a small mass gain (+16 ± 11 Gta–1) for a combined net change of –31 ± 12 Gta–1 (+0.08mma–1 SLE). The contribution of the three ice sheets to sea level is +0.05±0.03mma–1. The Antarctic ice shelves show corresponding mass changes of –95 ± 11 Gta–1 in WA and +142 ± 10Gta–1 in EA. Thinning at the margins of the Greenland ice sheet and growth at higher elevations is an expected response to increasing temperatures and precipitation in a warming climate. The marked thinnings in the Pine Island and Thwaites Glacier basins of WA and the Totten Glacier basin in EA are probably ice- dynamic responses to long-term climate change and perhaps past removal of their adjacent ice shelves. The ice growth in the southern Antarctic Peninsula and parts of EA may be due to increasing precipitation during the last century.

The Ice, Cloud and land Elevation Satellite (ICESat) in its 8 day repeat orbit mode provided data not only on the along-track surface slope, but also on the cross-track surface slope from adjacent repeat ground tracks. During the first 36 days of operation, four to five such repeat orbits occurred within 1 km in the cross-track direction. This provided an opportunity to use ICESat data to measure surface slope in the cross-track direction at 1 km scale. An algorithm was developed to calculate the cross-track surface slope. Combining the slopes in the cross-track and along-track directions gives a three-dimensional surface slope at 1 km scale. The along-track surface slope and surface roughness at 10km scale are also calculated. A comparison between ICESat surface elevation and a European Remote-sensing Satellite (ERS-1) 5 km digital elevation model shows a difference of 1–2 m in central Greenland where the surface slope is small, and >20m at the edge of Greenland where the surface slope is large. The large elevation difference at the edge is most likely due to the slope-induced error in radar altimeter measurement. Accurate surface slope data from ICESat will help to correct the slope-induced error of radar altimeter missions such as Geosat, ERS-1 and ERS-2.

Vapor-transfer theory is incorporated into a previous firn-densification model to investigate the effect of vapor-transfer processes on densification in firn within 10 m of the surface. The densification rate in the model is governed by the change of overburden pressure (determined by the accumulation rate), the firn temperature, and the temperature gradient. The time of exposure to temperature gradients at shallow depths is a critical factor determining the importance of vapor-transfer processes. In high-accumulation and high-temperature conditions such as for the Greenland ice sheet, the temperature gradient and vapor transfer are less important due to the shorter exposure times. The high summer temperatures dominate the rate of densification and annual variations in density. In low-accumulation and low-temperature conditions, such as for inland Antarctica, the vapor transfer driven by the temperature gradient has a stronger effect on the densification rate, and temperature-driven processes are less important. These factors determine both the rate of density increase with depth and the amplitudes of annual variations in density with depth.

Comparison of the distribution of seasonal variations in surface elevation derived from a firn-densification–elevation model with observed variations derived from ERS-1/-2 satellite radar altimetry shows close similarity in the patterns of the amplitude of the variations over the North Greenland ice sheet. The amplitudes of the seasonal variations decrease from west to east and from south to north, determined by the accumulation rate and the surface-temperature distribution pattern. Several methods of estimating the amplitude of the seasonal variation in the observations are compared, including the use of a three-frequency sinusoidal function derived from the modeled seasonal variation that is asymmetric. The resulting correlation coefficient between the observed amplitude, estimated with the three-frequency function, and the modeled amplitude is 0.66 and the slope is 0.7. Residual differences may be caused by interannual variability in accumulation and temperature and other approximations in the model.

An ice-sheet flowline model is used to simulate the flow of ice along two particle paths toward the onset to Ice Stream D, West Antarctica. One path is near the centre line of the main tributary to the ice stream, while the second passes by the Byrd Station borehole site. In this paper, we analyze the flow of the moderately fast-flowing tributaries in terms of ice-fabric anisotropy and estimate the steady-state ice-flow regions with the compatible developed crystal orientation fabrics along two particle paths. Comparison between modelled isochrones and internal layers detected from radio-echo sounding surveys in the area is used to suggest that flow upstream of the onset to Ice Stream D appears to have been stable since at least the Last Glacial Maximum.

Ice flow at a location in the equilibrium zone of the west-central Greenland Ice Sheet accelerates above the midwinter average rate during periods of summer melting. The near coincidence of the ice acceleration with the duration of surface melting, followed by deceleration after the melting ceases, indicates that glacial sliding is enhanced by rapid migration of surface meltwater to the ice-bedrock interface. Interannual variations in the ice acceleration are correlated with variations in the intensity of the surface melting, with larger increases accompanying higher amounts of summer melting. The indicated coupling between surface melting and ice-sheet flow provides a mechanism for rapid, large-scale, dynamic responses of ice sheets to climate warming.

Seasonal variations of firn density in ice-sheet firn layers have been attributed to variations in deposition processes or other processes within the upper firn. A recent high-resolution (mm-scale) density profile, measured along a 181 m core from Antarctica, showed small-scale density variations with a clear seasonal cycle that apparently was not related to seasonal variations in deposition or known near-surface processes (Gerland and others, 1999). A recent model of surface elevation changes (Zwally and Li, in press) produced a seasonal variation in firn densification, and explained the seasonal surface elevation changes observed by satellite radar altimeters. In this study, we apply our one-dimensional time-dependent numerical model of firn densification that includes a temperature-dependent formulation of firn densification based on laboratory measurements of grain growth. The model is driven by a steady-state seasonal surface temperature and a constant accumulation rate appropriate for the measured Antarctic ice core. The modeled seasonal variations in firn density show that the layers of snow deposited during spring to mid-summer that have the highest temperature history compress to the highest density, and the layers deposited during later summer to autumn that have the lowest temperature history compress to the lowest density. The initial amplitude of the seasonal difference of about 0.13 reduces to about 0.09 in 5 years and asymptotically to 0.0 at greater depth, which is consistent with the core measurements.

Seasonal and interannual variations in surface elevation at the Greenland summit are modeled using a new temperature-dependent formulation of firn densification and are compared with elevations from European Remote-sensing Satellite (ERS-1/-2) radar altimetry. The rate constant and activation energy, usually set as constants in the Arrhenius-type relation, are strongly temperature-dependent, based on measurements of crystal-growth rates. A multiplicative factor in the densification rate accounts for differences between densification and grain-growth rates and is chosen to match the modeled and measured density profiles from 0 to 40 m. The stronger temperature dependence produces a significant seasonal cycle in the densification rate in the upper firn. Much of the densification and consequent surface lowering occur within 3 months in late spring/early summer, followed by a build-up from accumulation. Modeled elevation changes, using automatic weather station measurements of temperature and accumulation and modeled precipitation, agree well with observations. The respective seasonal amplitudes are 18 and 25 cm peak-to-peak with minima in mid-summer. The minimum elevation in 1995 is driven mainly by a temporary accumulation decrease and secondarily by warmer temperatures. Increased compaction driven by a summer warming trend decreases the modeled elevation (1992–99) by 20 cm, but accumulation increases in latter years dominate the overall 4.2 cm a−1 trend.

Estimates of balance mass flux and the depth-averaged ice velocity through the cross-section aligned with the equilibrium line are produced for each of six drainage systems in Greenland. (The equilibrium line, which lies at approximately 1200 m elevation on the ice sheet, is the boundary between the area of net snow accumulation at higher elevations and the areas of net melting at lower elevations around the ice sheet.) Ice drainage divides and six major drainage systems are delineated using surface topography from ERS (European Remote Sensing) radar altimeter data. The net accumulation rate in the accumulation zone bounded by the equilibrium line is 399 Gt/yr and net ablation rate in the remaining area is 231 Gt/yr. (1 GigaTon of ice is 1090 kM(exp 3). The mean balance mass flux and depth-averaged ice velocity at the cross-section aligned with the modeled equilibrium line are 0.1011 Gt kM(exp -2)/yr and 0.111 km/yr, respectively, with little variation in these values from system to system. The ratio of the ice mass above the equilibrium line to the rate of mass output implies an effective exchange time of approximately 6000 years for total mass exchange. The range of exchange times, from a low of 3 ka in the SE drainage system to 14 ka in the NE, suggests a rank as to which regions of the ice sheet may respond more rapidly to climate fluctuations.

As the planet's principal cold traps, the martian polar regions have accumulated extensive mantles of ice and dust that cover individual areas of approximately 10(6) km2 and total as much as 3-4 km thick. From the scarcity of superposed craters on their surface, these layered deposits are thought to be comparatively young--preserving a record of the seasonal and climatic cycling of atmospheric CO2, H2O, and dust over the past approximately 10(5)-10(8) years. For this reason, the martian polar deposits may serve as a Rosetta Stone for understanding the geologic and climatic history of the planet--documenting variations in insolation (due to quasiperiodic oscillations in the planet's obliquity and orbital elements), volatile mass balance, atmospheric composition, dust storm activity, volcanic eruptions, large impacts, catastrophic floods, solar luminosity, supernovae, and perhaps even a record of microbial life. Beyond their scientific value, the polar regions may soon prove important for another reason--providing a valuable and accessible reservoir of water to support the long-term human exploration of Mars. In this paper we assess the current state of Mars polar research, identify the key questions that motivate the exploration of the polar regions, discuss the extent to which current missions will address these questions, and speculate about what additional capabilities and investigations may be required to address the issues that remain outstanding.

Elevations measured by the Mars Orbiter Laser Altimeter have yielded a high-accuracy global map of the topography of Mars. Dominant features include the low northern hemisphere, the Tharsis province, and the Hellas impact basin. The northern hemisphere depression is primarily a long-wavelength effect that has been shaped by an internal mechanism. The topography of Tharsis consists of two broad rises. Material excavated from Hellas contributes to the high elevation of the southern hemisphere and to the scarp along the hemispheric boundary. The present topography has three major drainage centers, with the northern lowlands being the largest. The two polar cap volumes yield an upper limit of the present surface water inventory of 3.2 to 4.7 million cubic kilometers.

Elevations from the Mars Orbiter Laser Altimeter (MOLA) have been used to construct a precise topographic map of the martian north polar region. The northern ice cap has a maximum elevation of 3 kilometers above its surroundings but lies within a 5-kilometer-deep hemispheric depression that is contiguous with the area into which most outflow channels emptied. Polar cap topography displays evidence of modification by ablation, flow, and wind and is consistent with a primarily H 2 O composition. Correlation of topography with images suggests that the cap was more spatially extensive in the past. The cap volume of 1.2 × 10 6 to 1.7 × 10 6 cubic kilometers is about half that of the Greenland ice cap. Clouds observed over the polar cap are likely composed of CO 2 that condensed out of the atmosphere during northern hemisphere winter. Many clouds exhibit dynamical structure likely caused by the interaction of propagating wave fronts with surface topography.

An updated compilation of oxygen-isotope ratio data for 562 sites in Antarctica shows a significant increase in the number of sites and an improvement in the representation of the coastal zone oyer previous versions. The data hase consists of ratio values (δ18O; multi-year mean 18O/16O relative to Standard Mean Ocean Water, in ‰) compiled as the dependent variable, together with data for the so-called independent variables: latitude, surface elevation, mean annual surface temperature and mean annual shortest distance to open ocean denoted by the 20% sea-ice concentration boundary. The problem of covariation between so-called independent variables is minimized using stepwise regression analyses. A general model is described using all the field data, and the regional variation at drainage-system scale is assessed by contrasting models for two physiographically distinct regions. in addition, entity-specific models are determined using data subsets for the conterminous grounded-ice and ice-shelf areas. Inversions of the specific models are applied to a 100 km grid data base to produce two contoured distributions of the ratio, one based on field data, and the other on remotely sensed data. The difference between these produces residuals that, relative to the summation of standard errors of the models, are small in most of the interior area of the ice sheet, and large in several areas of mountain and coastal regions, where interpolation and extrapolation of field data are particularly unreliable. Remotely sensed data generally produce ratio values which are isotopically cooler.

A compilation of mean values of the oxygen-isotope ratio relative 10 Standard mean ocean waterfor 22 siIes representative of conditions in northeastern Canada is complemented with data on mean annual surface temperature, latitude, surface elevation, and mean annual shortest distance to open ocean denoted by the 10% sea-ice concentration boundary. Stepwise regression analysis is used to develop a multivariate model suitable to inter the distribution ofin an area of complex topography and possibly mixed sourcing of advected water vapor. The best model is produced by a run in the backward mode at the 95% confidence level in which only temperature, latitude and distance to the open ocean remain in the model (the correlation coefficient is 0.915, the adjusted coefficient of determination is 0.809, the root mean square residual is 1.62). This model is similar to the bestpredictive model derived elsewhere for Greenland, suggesting a common principal source of advected moisture.

Mean oxygen-isotope ratio values relative to standard mean ocean water (δ18 O, in ‰) reported for 406 sites in Antarctica are compiled together with data on mean annual surface temperature, latitude, surface elevation, and mean annual shortest distance to open ocean denoted by the 20% sea-ice concentration boundary. Stepwise regression analyses with δ18O as the dependent variable are used as a model-building procedure based on statistical, rather than physical, criteria. Multivariate models sensitive to covariation between independent variables are defined using the whole dataset (N406 where N denotes the number of sites), as well as sub-sets for areas of conterminous grounded ice (N206) and ice shell (N110). The models show improvement over bivariate regression models. Distance to the open ocean enters all models at the second step. Inversions of the set and sub-set models applied to a database for 1351 gridpoint locations 100 km apart (it excludes the regions of Graham Land and eastern Palmer Land) are used to produce contoured distributions of δ18O. These may be used to assess the effects of atmospheric advection, as well as derive ice-flow adjustments for δ18O series obtained from deep-core or ablation-zone samples. Suggestions are made to improve model reliability.

The annual rate of net mass accumulation at the surface in the Antarctic and Greenland ice sheets is determined from firn emissivity based on Nimbus-5 ESMR and Nimbus-7 THIR data. In this study the determinations are limited to the areas of dry-snow facies and are based on a hyperbolic function of emissivity. Two coefficients of the function are selected for particular regions of each ice sheet after a comparison with field data selected for their reliability (82 stations in East Antarctica, 69 stations in West Antarctica and 89 stations in Greenland). Derived accumulation values are produced for grid-point locations 100 km apart which cover 56–94% of the dry-snow areas and 32–58% of the accumulation areas of each ice sheet. These values are compared with interpolated values obtained from the latest contoured compilations of field data. The means of derived values for East and West Antarctica are 12% and 39% larger, respectively, than the mean obtained from interpolated values, suggesting that the isopleth patterns as drawn in the compilation of field data lead to underestimates. The mean of derived values for Greenland is 5% smaller than the mean obtained from interpolated values, suggesting that the compilation of field data may lead to small overestimates that are within the error of determination and the variability of accumulation. Improving facies zonation and the determination of coefficients for the areas of upper percolation facies should improve these preliminary assessments.

Warming and Melting Mass loss from the ice sheets of Greenland and Antarctica account for a large fraction of global sea-level rise. Part of this loss is because of the effects of warmer air temperatures, and another because of the rising ocean temperatures to which they are being exposed. Joughin et al. (p. 1172 ) review how ocean-ice interactions are impacting ice sheets and discuss the possible ways that exposure of floating ice shelves and grounded ice margins are subject to the influences of warming ocean currents. Estimates of the mass balance of the ice sheets of Greenland and Antarctica have differed greatly—in some cases, not even agreeing about whether there is a net loss or a net gain—making it more difficult to project accurately future sea-level change. Shepherd et al. (p. 1183 ) combined data sets produced by satellite altimetry, interferometry, and gravimetry to construct a more robust ice-sheet mass balance for the period between 1992 and 2011. All major regions of the two ice sheets appear to be losing mass, except for East Antarctica. All told, mass loss from the polar ice sheets is contributing about 0.6 millimeters per year (roughly 20% of the total) to the current rate of global sea-level rise.

ICESat and Envisat altimetry data provide measurements of the instantaneous sea surface height (SSH) across the Arctic Ocean, using lead and open water elevation within the sea ice pack. First, these data were used to derive two independent mean sea surface (MSS) models by stacking and averaging along-track SSH profiles gathered between 2003 and 2009. The ICESat and Envisat MSS data were combined to construct the high-resolution ICEn MSS. Second, we estimate the 5.5-year mean dynamic topography (MDT) of the Arctic Ocean by differencing the ICEn MSS with the new GOCO02S geoid model, derived from GRACE and GOCE gravity. Using these satellite-only data we map the major features of Arctic Ocean dynamical height that are consistent with in situ observations, including the topographical highs and lows of the Beaufort and Greenland Gyres, respectively. Smaller-scale MDT structures remain largely unresolved due to uncertainties in the geoid at short wavelengths.

Firn-temperature profiles are calculated in a thermal model using continuous surface temperatures derived from automatic weather station data and passive-microwave data in the Greenland summit region during the period 1987–99. the results show that significant interannual variations of mean summer (June–August) and annual temperatures occur in the top 15 m, in addition to the normal seasonal cycle of firn temperature. At 5 m depth, the seasonal cycle is damped to 13% of the surface seasonal range, but even at 15m about 1% or 0.6˚C of the seasonal cycle persists. Both summer and mean annual temperatures decrease from 1987 to 1992, followed by a general increasing trend. Interannual variability is 5˚C at the surface, but is dampened to 3.2˚C at 5 m depth and 0.7˚C at 15 m depth. Dampening of the interannual variability with depth is slower than dampening of the seasonal cycle, because of the longer time constant of the interannual variation. the warmer spring and summer temperatures experienced in the top 5 m, due to both the seasonal cycle and interannual variations, affect the rate of firn densification, which is non-linearly dependent on temperature. During the 12 year period 1987–99, the annual mean surface temperature is –29.2˚C, and the annual mean 15 m temperature is –30.1˚C, which is >1˚C warmer than a 15 mborehole temperature representing the period around 1959 and warmer than the best-fit temperature history by Alley and Koci (1990) back to AD 1500.

Geosat-altimeter wave forms from the Greenland and Antarctic ice sheets are analyzed using an algorithm based upon a combined surface-and volume-scattering model. The results demonstrate that sub-surface volume-scattering occurs over major parts of the ice sheets. Quantitative estimates of geographic variations in the near-surface ice-sheet properties are derived by retracking individual altimeter wave forms. The derived surface properties correlate with elevation, latitude and microwave brightness-temperature data. Specifically, the extinction coefficient of snow obtained by this method varies from 0.48 to 0.13 m−1 over the latitudes from 65° to 72°N on the central part of the Greenland ice sheet and from 0.20 to 0.10 m−1 over a section of Wilkes Land in East Antarctica where the elevation increases from 2550 to 3150 m.Analysis of passive-microwave data over East Antarctica shows that the brightness temperature increases with elevation as the extinction coefficient decreases. Larger snow grain-sizes occur at lower elevations of the ice sheet because of higher mean annual temperatures. The larger grain-sizes increase the extinction coefficient of snow and decrease the emitted energy (brightness temperature) from greater snow depths. The passive-microwave data are also used to determine the average number of melt d year−1 (1979–87) for the central part of the Greenland ice sheet. For latitudes from 65° to 68.5° N, the average number of melt days decreases from 3.5 to 0.25 d year, whereas no melt events are observed for latitudes above 69°N over the 8 year period. Snow subjected to alternate melting and freezing has enhanced grain-sizes compared to that of dry snow. This accounts for the larger values and larger spatial variations of ke on the Greenland ice sheet compared to East Antarctica, where surface temperatures are never high enough to cause surface melting.

We have recently completed an analysis that examines in detail the spatial and temporal variations in global sea-ice coverage from 26 October 1978, through 20 August 1987. The sea-icemeasurements we analyzed are derived from data collected by a multispectral, dual-polarized, constant incidence-angle microwave imager, the Scanning Multichannel Microwave Radiometer (SMMR) on board the NASA Nimbus 7 satellite. The characteristics of the SMMR have permitted a more accurate calculation of total sea-ice concentrations (fraction of ocean area covered by sea ice) than earlier single-channel instruments and, for the first time, a determination of both multiyear sea-ice concentrations and physical temperatures of the sea-ice pack. An estimate of the SMMR wintertime total ice concentration accuracy of ± 7% in both hemispheres has been obtained. As this is an improvement over the estimated accuracies of previous microwave sensors, we are able to present improved calculations of the sea-ice extents (areas enclosed by the 15% ice concentration boundaries), sea-ice concentrations, and open-water areas within the ice margins. This analysis will be published in a book, Arctic and Antarctic sea ice, 1978–1987: satellite passive microwave observations and analysis, due for publication in1992. Some highlights from the analysis are presented in this paper.

A new, detailed marine gravity field for the persistently ice-covered Arctic Ocean, derived entirely from satellite data, reveals important new tectonic features in both the Amerasian and Eurasian basins. Reprocessed Geoscience Laser Altimeter System (GLAS) data collected by NASA's Ice Cloud and land Elevation Satellite (ICESat) between 2003 and 2005 have been combined with 8 years worth of retracked radar altimeter data from ESA's ERS-2 satellite to produce the highest available resolution gravity mapping of the entire Arctic Ocean complete to 86 degrees N. This ARCtic Satellite-only (ARCS) marine gravity field uniformly and confidently resolves marine gravity to wavelengths as short as 35 km. ARCS relies on a Gravity Recovery and Climate Experiment (GRACE)-only satellite gravity model at long (> 580 km) wavelengths and plainly shows tectonic fabric and numerous details imprinted in the Arctic seafloor, in particular, in the enigmatic Amerasian Basin (AB). For example, in the Makarov Basin portion of the AB, two north-south trending lineations are likely clues to the highly uncertain seafloor spreading history which formed the AB.

[1] Sea ice freeboard heights in the Weddell Sea of Antarctica are derived from the Ice, Cloud, and Land Elevation Satellite (ICESat) laser altimeter measurements, which have a unique range precision to flat surfaces of 2 cm within 70 m footprints spaced at 172 m along track. Although elevations of flat surfaces can be obtained to an accuracy of � 10 cm (1s) per footprint, direct determination of freeboard heights is precluded by errors in knowledge of the geoid and temporal variability of the ocean surface. Therefore freeboards are determined relative to an ocean reference level detected over areas of open water and very thin ice within the sea ice pack using an along-track filtering method. The open water/thin ice selections show good agreement in the combined analysis of ICESat segments and Envisat Synthetic Aperture Radar (SAR) imagery. The average residual between the ICESat-measured ocean level and the EGM96 geoid is 1.4 m. Estimates of snow depth on the sea ice from AMSR-E passive microwave data along with nominal densities of snow, water, and sea ice are used to estimate sea ice thickness. Four periods of ICESat data in May–June (fall) and October–November (late winter) of 2004 and 2005 between longitudes 298E and 360E are analyzed. In the fall the mean freeboards are 0.28 m in 2004 and 0.29 m in 2005, and the mean thicknesses are 1.33 m in 2004 and 1.52 m in 2005. In late winter the freeboards grew to 0.37 m in 2004 and 0.35 in 2005, and the thicknesses grew to 2.23 m in 2004 and 2.31 m in 2005. The interannual differences in freeboard are small, and the larger interannual change in estimated thickness mainly represents differences in the snow depth estimates. Seasonal changes in the spatial patterns of freeboard and thickness over the 4 months correlate with the expected circulation of sea ice in the Weddell Sea, as indicated by sea ice velocity fields.

Seasonal mean changes in the surface elevation of the ablation zone of West Greenland to 72°N between spring 1985 and summer 1986 are measured using radar altimeter data from the 18-month Geosat Geodetic Mission. Semi-variograms are used to estimate the noise in the data as a function of position on the ice sheet. Mean elevation changes are computed by averaging the elevation differences measured at points where orbits ascending in latitude are later crossed by orbits descending in latitude (or the reverse), with each cross-over difference weighted in proportion to the inverse square of the noise level in the neighborhood of the cross-over point. The mean surface elevation of the ablation zone, relative to spring 1985, ranged from 1.5 ± 0.6 m lower during summer 1985 to 1.7 ± 0.4 m higher during spring 1986.

Analysis of near-coincident ICESat and RADARSAT imagery shows that the retrieved elevations from the laser altimeter are sensitive to new openings (containing thin ice or open water) in the sea ice cover as well as to surface relief of old and first-year ice. The precision of the elevation estimates, measured over relatively flat sea ice, is approx. 2 cm Using the thickness of thin-ice in recent openings to estimate sea level references, we obtain the sea-ice free-board along the altimeter tracks. This step is necessitated by the large uncertainties in the time-varying sea surface topography compared to that required for accurate determination of free-board. Unknown snow depth introduces the largest uncertainty in the conversion of free-board to ice thickness. Surface roughness is also derived, for the first time, from the variability of successive elevation estimates along the altimeter track Overall, these ICESat measurements provide an unprecedented view of the Arctic Ocean ice cover at length scales at and above the spatial dimension of the altimeter footprint.

Publisher Summary This chapter discusses the dynamics of ice sheet and mass balance. Ice-sheet mass balance is the difference between the mass input and the mass output. The total balance includes both the surface mass balance processes and the ice-flow components. Mass is added to the surface from snowfall, condensation, and occasional rainfall. Mass is removed by evaporation, surface and bottom melting, water runoff, and iceberg discharge. The chapter discusses the radar altimeter measurement of ice-sheet surface elevations. Several methods have been used to create the digital elevation models (DEMs) of ice sheet topography from satellite altimetry. All methods involve mapping the data onto evenly spaced grids and correcting for the slope error, either before or after the mapping. Full DEMs of Greenland and Antarctica have been produced using triangularization or gridding procedures. A three-step inversion technique has also been demonstrated in small regions by Remy et al: first a large-scale reference surface is estimated, the residuals are mapped related to the undulations, and finally iteratively corrected for the slope error.

Over the past decade, the Arctic Ocean and Beaufort Sea ice pack has been less extensive and thinner than has been observed during the previous 35 years [e.g., Wadhams and Davis, 2000; Tucker et al., 2001; Rothrock et al., 1999; Parkinson and Cavalieri, 2002; Comiso, 2002]. During the summers of 2007 and 2008, the ice extents for both the Beaufort Sea and the Northern Hemisphere were the lowest on record. Mechanisms causing recent sea ice change in the Pacific Arctic and the Beaufort Sea are under investigation on many fronts [e.g., Drobot and Maslanik, 2003; Shimada et al., 2006]; the mechanisms include increased ocean surface warming due to Pacific Ocean water inflow to the region and variability in meteorological and surface conditions. However, in most studies addressing these events, the impact of sea ice dynamics, specifically deformation, has not been measured in detail.

The principal conclusion of C. H. Davis et al . ([1][1]), that the elevation of the Southern Greenland ice sheet increased at a smaller rate than the 23 cm/year we previously reported ([2][2], [3][3]), is in qualitative agreement with our new analysis. Their statement, however, that “surface

Over century time scales, the primary effect of ice-sheet/climate-change interactions is vertical growth or shrinkage of the ice in response to changes in precipitation and surface heat flux. Because the dynamic response of large ice masses is generally slower than climate variations experienced during the last few centuries, the ice sheets are unlikely to be in equilibrium with today’s climate. Uncertainty in the current mass balance has been large, at least ±30% or ±2 mm yr−1 in sea-level equivalent. Estimates of annual snow accumulation, iceberg discharge, and peripheral melting of the Antarctic ice sheet (personal communication from S. Jacobs) would suggest a negative mass-balance equivalent to +2 mm yr−1 of sea-level rise, in contrast to other estimates of a small positive balance for both Antarctica (−0.6 ±0.6 mm yr−1 sea level) and Greenland (−0.1 ±0.4 mm yr−1 sea level) (Meier and others, 1985). For some years, measurement of changes in ice-sheet surface elevation by satellite altimetry has been noted as a potential means of determining the overall ice-sheet mass balance and investigating regional variations. Difficulties in deriving elevation change from a set of sequential measurements from several satellites have been primarily a result of the limited precision (about 1–2 m) of satellite radar altimetry, residual orbit errors, and relative uncertainties in the gravity fields and geoid reference levels used for different satellites. However, recent radar altimeter measurements by the U.S. Navy Geosat to 72°N and 72°S provide a sufficient density of repeated measurements of ice elevation for analysis of elevation change during the life of the satellite. The average elevation change at 224 267 orbital crossovers over southern Greenland is +28.3 ±0.4 cm yr−1. The largest values are observed near the summit and over the southern dome, with smaller values in the saddle region and toward the margins. Local gridded values of thickening and thinning rates agree with estimates from surface studies in the vicinity of the EGIG traverse (Steckel, 1977), Dye 3 (Reeh, 1985), and the OSU survey (Kostecka and Whillans, 1988) within the error limits of the respective measurements. Comparisons with elevations measured by the Seasat altimeter in 1978 and continuing Geosat measurements provide information on the temporal continuity of the thickening. The spatial distribution of the elevation changes is used to estimate an average thickening rate and the current rate of oceanic depletion.

Many of the major advances in glaciology during the past 50 years have followed the development and application of new technology for viewing and measuring various characteristics of ice. Microscopes to study ice crystals, radars to probe the internal structure of large ice masses, mass spectrometers to analyze the atomic composition of ice cores, and satellite sensors to measure the global distribution of ice are some of the tools readily adapted by glaciologists. Today, new tools include microcomputers for automatic data logging, large-memory computers for data processing and numerical modeling, sensitive instruments for ice analysis, and satellite sensors for large-scale ice observations. In the future, continued advances in key technologies will help guide the evolution of science questions considered by glaciologists, expanding our view of ice, its fundamental properties, its interactions within the ice–ocean–land–atmosphere system, and its role in the evolution of our global environment.

The Seasat radar altimeter, which was designed to measure ranges to the sea surface, has also provided the most accurate available maps of ice-sheet eleva­ tion. Seasat operated for three months during the austral winter of 1978, when Antarctica was girdled by sea ice. As the satellite approached the continent from the ocean, the altimeter obtained very strong reflections from the sea ice, and it continued to measure ranges to the sea ice even for a short time after passing over the ice front. The measured ranges are oblique distances to the nearest portion of sea ice, and the sequence of oblique ranges gives the position of the sea ice along the ice front. After the satellite crossed the ice front travelling sea­ ward, oblique ranges were measured to the nearest portion of ice shelf. Examples are given from Seasat orbits crossing the Amery and Fimbul ice shelves. The entire Seasat data set provides an opportunity for mapping most of the East Antarctic coastline to an absolute horizontal accuracy of ± (0.1 + 1) km. The resulting map would effectively be an instantan­ eous view of the ice margin, requiring no adjustment for ice movement during the survey period. INTRODUCTION The radar altimeter aboard Seasat was designed to measure ranges to the ocean surface. It transmitted short radar pulses with a beam width (to 3 dB attenu­ ation) of 1. 6 degrees of arc. Pulse length was 3.125 ns (equivalent to -1 m), and pulse rate was 1 020 s-l. Ranges were obtained from the time delay between pulse transmission and receipt of the reflected pulse. Return energy was measured during a data-acquisition window of about 188 ns duration. High resolution was achieved by dividing this window into 63 time gates. In operation, the altimeter initially undertook a low-resolution search for the ocean surface, and then located the center of the high resolution acquisition window to coincide with the time of the return-pulse leading edge. There­ after, the altimeter maintained track of the ocean *Present address: California Institute of Technology, Pasadena, California 91109, USA. surface by continually adjusting the timing of the acquisition window as dictated by the rate of change of previous measured ranges. The tracking worked well over the oceans, where measured ranges changed very slowly. However, over sloping or undulating terrain, the servo-tracking circuit was not sufficiently agile to monitor the rapidly changing ranges, and the alti­ meter frequently lost track of the return pulse. This resulted in short periods (usually a few seconds) when no useful data were obtained. Nevertheless, more than 600 000 surface elevations were measured on the portions of Greenland and Antarctica that lie between ±72' latitude. These are undergoing systematic correction for tracking errors and slope-induced errors as part of the research program of the Oceans and Ice Branch at the Goddard Space Flight Center (Zwally and others in press) . In this paper, we show how radar-altimetry data can be used to map ice-sheet margins. Seaward ice fronts of the ice shelves, and ice margins generally, are associated with a step change in surface eleva­ tion. Mapping the ice front would be straightforward if the altimeter responded immediately to the abrupt change in surface elevation. However, the Seasat radar altimeter responded slowly and usually lost track for a few seconds. Nevertheless, the ice-front position can be inferred accurately from the rela­ tionship between satellite position and measured range prior to losing track. As the satellite passed over the ice front from seaward, it continued to record the return signal from the sea, or sea ice, that lay within the seaward portion of the beam­ limited radar footprint (radius approximately 11 km) . Consequently, the measured range was the oblique distance from the satellite to the sea surface next to the ice front. The variation of measured distance with time gives the ice-front position to within 100 to 1000 m, depending on data quality and the length of time during which oblique ranges were measured. In addition, the data give an indication of the shape of the nearby ice front as related to the orbit track. When the satellite crossed the ice front travel­ ling seaward, oblique ranges were similarly obtained Jet Propulsion Laboratory, 4800 Oak Grove Drive

Passive microwave images of the polar regions, first produced after the launch of the Nimbus-5 Electrically Scanning Microwave Radiometer (ESMR)in December 1972, have become a valuable new source of polar information. Some of the potential applications of this new capability were anticipated. Of these, the sensing of sea ice through clouds and the polar night is probably the most important application for polar research and for operations on the polar seas. Other applications, such as the measurement of certain near-surfaceice sheet parameters, have been formulated more recently. Measurement of various ocean surface parameters is expected from the forthcoming multifrequency microwave observations. Undoubtedly additional uses of passive microwave datawill be conceived and developed. Two remarkable aspects of satellite-borne microwave radiometers are the complete spatial detail obtained by the scanning sensors and the temporal detail provided by continual coverage. For example, the observations of detailed microwave emission patterns over the Antarctic ice sheet should yield information that could not be obtained by surface or even aircraft measurements. Sequences of images produced at three-day intervalsreveal short-term ice sheet and sea ice phenomena that would otherwise be missed.