Overview of Ice-Sheet Mass Balance and Dynamics from ICESat Measurements

The primary purpose of the ICESat mission was to determine the present-day mass balance of the Greenland and Antarctic ice sheets, identify changes that may be occurring in the surface-mass flux and ice dynamics, and estimate their contributions to global sea-level rise. Although ICESat's three lasers were planned to make continuous measurements for 3 to 5 years, the mission was re-planned to operate in 33-day campaigns 2 to 3 times each year following failure of the first laser after 36 days. Seventeen campaigns were conducted with the last one in the Fall of 2009. Mass balance maps derived from measured ice-sheet elevation changes show that the mass loss from Greenland has increased significantly to about 170 Gt/yr for 2003 to 2007 from a state of near balance in the 1990's. Increased losses (189 Gt/yr) from melting and dynamic thinning are over seven times larger'than increased gains (25 gt/yr) from precipitation. Parts of the West Antarctic ice sheet and the Antarctic Peninsula are losing mass at an increasing rate, but other parts of West Antarctica and the East Antarctic ice sheet are gaining mass at an increasing rate. Increased losses of 35 Gt/yr in Pine Island, Thwaites-Smith, and Marie-Bryd.Coast are more than balanced by gains in base of Peninsula and ice stream C, D, & E systems. From the 1992-2002 to 2003-2007 period, the overall mass balance for Antarctica changed from a loss of about 60 Gt/yr to near balance or slightly positive

Satellite-derived ice data sets no. 1: Antarctic monthly average microwave brightness temperatures and sea-ice concentrations, 1973 - 1976

A summary data set concerning 4 years of Antarctic sea-ice conditions was created and is available on magnetic tape. The data were derived from electrically scanning microwave radiometer brightness temperatures and were mapped into a polar stereographic grid enclosing the 50 deg S latitude circle. The grid size varies from about 32 by 32 sq km at the poles to about 28 by 28 sq km at 50 deg S. The microwave brightness temperatures of Antarctic sea ice are predominantly characteristic of first-year ice with an emissivity of 0.92 at 19 GHz frequency. Sea ice concentrations were calculated from the brightness temperature data for each grid element with an algorithm that uses an emissivity value of 0.92 and an ice physical temperature estimate from climatological surface air temperatures. Monthly, multiyear monthly, and yearly maps of brightness temperatures and sea ice concentrations were created for the 4 years, except for 7 months for which useable data were insufficient

Satellite-derived ice data sets no. 2: Arctic monthly average microwave brightness temperatures and sea ice concentrations, 1973-1976

A summary data set for four years (mid 70's) of Arctic sea ice conditions is available on magnetic tape. The data include monthly and yearly averaged Nimbus 5 electrically scanning microwave radiometer (ESMR) brightness temperatures, an ice concentration parameter derived from the brightness temperatures, monthly climatological surface air temperatures, and monthly climatological sea level pressures. All data matrices are applied to 293 by 293 grids that cover a polar stereographic map enclosing the 50 deg N latitude circle. The grid size varies from about 32 X 32 km at the poles to about 28 X 28 km at 50 deg N. The ice concentration parameter is calculated assuming that the field of view contains only open water and first-year ice with an ice emissivity of 0.92. To account for the presence of multiyear ice, a nomogram is provided relating the ice concentration parameter, the total ice concentration, and the fraction of the ice cover which is multiyear ice

Changes in the Mass Balance of the Greenland Ice Sheet in a Warming Climate During 2003-2009

Mass changes of the Greenland ice sheet (GIS) derived from ICESat and GRACE data both show that the net mass loss from GIS during 2003-2009 is about 175 Gt/year, which contributes 0.5mm/yr global sea-level rise. The rate of mass loss has increased significantly since the 1990's when the GIS was close to mass balance. Even though the GIS was close to mass balance during the 1990's, it was already showing characteristics of responding to8 warmer climate, specifically thinning at the margins and thickening inland at higher elevations. During 2003-2009, increased ice thinning due to increases in melting and acceleration of outlet glaciers began to strongly exceed the inland thickening from increases in accumulation. Over the entire GIS, the mass loss between the two periods, from increased melting and ice dynamics, increased by about 190 Gt/year while the mass gain, from increased precipitation and accumulation, increased by only about 15Gt/year. These ice changes occurred during a time when the temperature on GIS changed at rate of about 2K/decade. The distribution of elevation and mass changes derived from ICESat have high spatial resolution showing details over outlet glaciers, by drainage systems, and by elevation. However, information on the seasonal cycle of changes from ICESat data is limited, because the ICESat lasers were only operated during two to three campaigns per year of about 35 days duration each. In contrast, the temporal resolution of GRACE data, provided by the continuous data collection, is much better showing details of the seasonal cycle and the inter-annual variability. The differing sensitivity of the ICESat altimetry and the GRACE gravity methods to motion of the underlying bedrock from glacial isostatic adjustment (GIA) is used to evaluate the GIA corrections provided by models. The two data types are also combined to make estimates of the partitioning of the mass gains and losses among accumulation, melting, and ice discharge from outlet glaciers

Matrix partitioning and EOF/principal component analysis of Antarctic Sea ice brightness temperatures

A field of measured anomalies of some physical variable relative to their time averages, is partitioned in either the space domain or the time domain. Eigenvectors and corresponding principal components of the smaller dimensioned covariance matrices associated with the partitioned data sets are calculated independently, then joined to approximate the eigenstructure of the larger covariance matrix associated with the unpartitioned data set. The accuracy of the approximation (fraction of the total variance in the field) and the magnitudes of the largest eigenvalues from the partitioned covariance matrices together determine the number of local EOF's and principal components to be joined by any particular level. The space-time distribution of Nimbus-5 ESMR sea ice measurement is analyzed

Modeling of Firn Compaction for Estimating Ice-Sheet Mass Change from Observed Ice-Sheet Elevation Change

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 kg/cu m, which along with 900 kg/cu m for the dynamic/ablation-driven changes gives average densities ranging from 680 to 790 kg/cu m. 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

Surface topography of the Greenland Ice Sheet from satellite radar altimetry

Surface elevation maps of the southern half of the Greenland subcontinent are produced from radar altimeter data acquired by the Seasat satellite. A summary of the processing procedure and examples of return waveform data are given. The elevation data are used to generate a regular grid which is then computer contoured to provide an elevation contour map. Ancillary maps show the statistical quality of the elevation data and various characteristics of the surface. The elevation map is used to define ice flow directions and delineate the major drainage basins. Regular maps of the Jakobshavns Glacier drainage basin and the ice divide in the vicinity of Crete Station are presented. Altimeter derived elevations are compared with elevations measured both by satellite geoceivers and optical surveying

Persistent surface snowmelt over Antarctica (1987–2006) from 19.35 GHz brightness temperatures

[1] Persistent melting (e.g., continuing for more than three days or for one consecutive day and night) is mapped in Antarctica (1987–2006) using night- and day-time Special Sensor Microwave Imager brightness temperatures (Tb) at 19.35 GHz, horizontal polarization. Snowmelt is indicated when Tb and relative daily difference exceed threshold values, respectively Tc and ΔT, computed for each pixel and year, or when both daytime and nighttime Tb exceed Tc. Results from an electromagnetic model suggest that the minimum detectable liquid water content ranges between 0.2 and 0.5%, in volume. We find that melting areas have been moving inland since 1987. A first-time extensive melting (1987–2006) is detected over the Transantarctic Mountains on January 2005, 875 Km inland and 2000 m above sea level. Melting extent and index have been decreasing over Antarctica, since 1987, although either positive and negative trends are observed from a sub-continental scale analysis

Microwave emission from snow and glacier ice

The microwave brightness temperature for snow fields was studied assuming that the snow cover consists of closely packed scattering spheres which do not interact coherently. The Mie scattering theory was used to compute the volume scattering albedo. It is shown that in the wavelength range from 0.8 to 2.8 cm, most of the micro-radiation emanates from a layer 10 meters or less in thickness. It is concluded that it is possible to determine snow accumulation rates as well as near-surface temperature

Polar microwave brightness temperatures from Nimbus-7 SMMR: Time series of daily and monthly maps from 1978 to 1987

A time series of daily brightness temperature gridded maps (October 25, 1978 through August 15, 1987) were generated from all ten channels of the Nimbus-7 Scanning Multichannel Microwave Radiometer orbital data. This unique data set can be utilized in a wide range of applications including heat flux, ocean circulation, ice edge productivity, and climate studies. Two sets of data in polar stereographic format are created for the Arctic region: one with a grid size of about 30 km on a 293 by 293 array similar to that previously utilized for the Nimbus-5 Electrically Scanning Microwave Radiometer, while the other has a grid size of about 25 km on a 448 by 304 array identical to what is now being used for the DMSP Scanning Multichannel Microwave Imager. Data generated for the Antaractic region are mapped using the 293 by 293 grid only. The general technique for mapping, and a quality assessment of the data set are presented. Monthly and yearly averages are also generated from the daily data and sample geophysical ice images and products derived from the data are given. Contour plots of monthly ice concentrations derived from the data for October 1978 through August 1987 are presented to demonstrate spatial and temporal detail which this data set can offer, and to show potential research applications