Search NASA⌕ Search

SEARCH · Search NASA

Results for “ice sheets”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

Airborne Resurveys of the Southern Greenland Ice Sheet

Airborne laser-altimeter flight lines from 1993 over southern parts of the ice sheet were resurveyed with almost complete repeat coverage. In 1993 and 1994, NASA surveyed the entire Greenland ice sheet by airborne laser altimeter, obtaining surface-elevation profiles with root mean square (rms) accuracies of 10 cm or better (Krabill 1995) along flight lines that crossed all the major catchment basins. In 1998, the ten flight lines flown in 1993 in the south of Greenland were resurveyed with about 99% repeat coverage; flight lines in the north will be resurveyed in 1999. Additional flights in 1998 were over glaciers, identified by E. Rignot, where existing SAR data give information on ice motion.

Krabill, William B.↗

Seasonal Monitoring of Melt and Accumulation Within the Deep Percolation Zone of the Greenland Ice Sheet and Comparison with Simulations of Regional Climate Modeling

Increasing melt over the Greenland Ice Sheet (GrIS) recorded over the past several years has resulted in significant changes of the percolation regime of the ice sheet. It remains unclear whether Greenland's percolation zone will act as a meltwater buffer in the near future through gradually filling all pore space or if near-surface refreezing causes the formation of impermeable layers, which provoke lateral runoff. Homogeneous ice layers within perennial firn, as well as near-surface ice layers of several meter thickness have been observed in firn cores. Because firn coring is a destructive method, deriving stratigraphic changes in firn and allocation of summer melt events is challenging. To overcome this deficit and provide continuous data for model evaluations on snow and firn density, temporal changes in liquid water content and depths of water infiltration, we installed an upward-looking radar system (upGPR) 3.4 m below the snow surface in May 2016 close to Camp Raven (66.4779 deg N, 46.2856 deg W) at 2120 m a.s.l. The radar is capable of quasi-continuously monitoring changes in snow and firn stratigraphy, which occur above the antennas. For summer 2016, we observed four major melt events, which routed liquid water into various depths beneath the surface. The last event in mid-August resulted in the deepest percolation down to about 2.3 m beneath the surface. Comparisons with simulations from the regional climate model MAR are in very good agreement in terms of seasonal changes in accumulation and timing of onset of melt. However, neither bulk density of near-surface layers nor the amounts of liquid water and percolation depths predicted by MAR correspond with upGPR data. Radar data and records of a nearby thermistor string, in contrast, matched very well for both timing and depth of temperature changes and observed water percolations. All four melt events transferred a cumulative mass of 56 kg/(sq. m) into firn beneath the summer surface of 2015. We find that continuous observations of liquid water content, percolation depths and rates for the seasonal mass fluxes are sufficiently accurate to provide valuable information for validation of model approaches and help to develop a better understanding of liquid water retention and percolation in perennial firn.

Heilig, Achim↗

Assessment of Antarctic Ice-Sheet Mass Balance Estimates: 1992 - 2009

Published mass balance estimates for the Antarctic Ice Sheet (AIS) lie between approximately +50 to -250 Gt/year for 1992 to 2009, which span a range equivalent to 15% of the annual mass input and 0.8 mm/year Sea Level Equivalent (SLE). Two estimates from radar-altimeter measurements of elevation change by European Remote-sensing Satellites (ERS) (+28 and -31 Gt/year) lie in the upper part, whereas estimates from the Input-minus-Output Method (IOM) and the Gravity Recovery and Climate Experiment (GRACE) lie in the lower part (-40 to -246 Gt/year). We compare the various estimates, discuss the methodology used, and critically assess the results. Although recent reports of large and accelerating rates of mass loss from GRACE=based studies cite agreement with IOM results, our evaluation does not support that conclusion. We find that the extrapolation used in the published IOM estimates for the 15 % of the periphery for which discharge velocities are not observed gives twice the rate of discharge per unit of associated ice-sheet area than the 85% faster-moving parts. Our calculations show that the published extrapolation overestimates the ice discharge by 282 Gt/yr compared to our assumption that the slower moving areas have 70% as much discharge per area as the faster moving parts. Also, published data on the time-series of discharge velocities and accumulation/precipitation do not support mass output increases or input decreases with time, respectively. Our modified IOM estimate, using the 70% discharge assumption and substituting input from a field-data compilation for input from an atmospheric model over 6% of area, gives a loss of only 13 Gt/year (versus 136 Gt/year) for the period around 2000. Two ERS-based estimates, our modified IOM, and a GRACE-based estimate for observations within 1992 to 2005 lie in a narrowed range of +27 to - 40 Gt/year, which is about 3% of the annual mass input and only 0.2 mm/year SLE. Our preferred estimate for 1992-2001 is - 47 Gt/year for West Antarctica, + 16 Gt/year for East Antarctica, and -31 Gt/year overall (+0.1 mm/year SLE), not including part of the Antarctic Peninsula (1.07 % of the AIS area)

Zwally, H. Jay↗

Evolving Understanding of Antarctic Ice-Sheet Physics and Ambiguity in Probabilistic Sea-Level Projections

Mechanisms such as ice-shelf hydrofracturing and ice-cliff collapse may rapidly increase discharge from marine-based ice sheets. Here, we link a probabilistic framework for sea-level projections to a small ensemble of Antarctic ice-sheet (AIS) simulations incorporating these physical processes to explore their influence on global-mean sea-level (GMSL) and relative sea-level (RSL). We compare the new projections to past results using expert assessment and structured expert elicitation about AIS changes. Under high greenhouse gas emissions (Representative Concentration Pathway [RCP] 8.5), median projected 21st century GMSL rise increases from 79 to 146 cm. Without protective measures, revised median RSL projections would by 2100 submerge land currently home to 153 million people, an increase of 44 million. The use of a physical model, rather than simple parameterizations assuming constant acceleration of ice loss, increases forcing sensitivity: overlap between the central 90% of simulations for 2100 for RCP 8.5 (93-243 cm) and RCP 2.6 (26-98 cm) is minimal. By 2300, the gap between median GMSL estimates for RCP 8.5 and RCP 2.6 reaches >10 m, with median RSL projections for RCP 8.5 jeopardizing land now occupied by 950 million people (versus 167 million for RCP 2.6). The minimal correlation between the contribution of AIS to GMSL by 2050 and that in 2100 and beyond implies current sea-level observations cannot exclude future extreme outcomes. The sensitivity of post-2050 projections to deeply uncertain physics highlights the need for robust decision and adaptive management frameworks.

Kopp, Robert E.↗

A New Glacial Varve Chronology Along the Southern Laurentide Ice Sheet That Spans the Younger Dryas–Holocene Boundary

Glacial varves can detail ice-margin positions and provide a proxy for meltwater discharge at resolutions comparable to those of the Greenland ice core archives, and thus they can be critical paleorecords for assessing the response of both ancient and modern ice sheets to climate change. Here we provide an ~1500 yr varve chronology straddling the Younger Dryas (YD)–Holocene boundary (11.65 cal. kyr B.P.), the first such chronology in North America. The varves are from glacial Lake Agassiz (central North America). The chronology is pinned on accelerator mass spectrometry radiocarbon-dated terrestrial macrofossils at the base of a widespread red-clay bed deposited during flooding from the Lake Superior basin. We illustrate the utility of this record by examining ice-margin retreat and melting through the late Younger Dryas and across the Holocene boundary. The ice margin receded at a constant rate, not only during the late YD, but for at least 300 yr after the onset of the Holocene. In contrast, varve thicknesses increased at the boundary, and a moraine formed over a 50 yr period, perhaps in response to the warming climate. Our expectation is that this time series will continue to be developed, expanded, and refined because it promises to be the longest and most geographically extensive glacial varve data set in North America.

Glacial varves↗

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric↗

Comparisons of Satellite and Airborne Altimetry with Ground-Based Data from the Interior of the Antarctic Ice Sheet

A series of traverses has been conducted for validation of the NASA Ice, Cloud, and land Elevation Satellite 2 on the flat interior of the Antarctic ice sheet. GNSS data collected on three separate 88S Traverses intersect 20% of the ICESat-2 reference ground tracks and have precisions of better than ±7 cm and biases of less than ~4 cm. Data from these traverses were used to assess heights from ICESat-2, CryoSat-2, and Airborne Topographic Mapper (ATM). ICESat-2 heights have better than ±3.3 cm bias and better than ±7.2 cm precision. ATM heights have better than 9.3 cm bias and better than ±9.6 cm precision. CryoSat-2 heights have -38.9 cm of bias and ± 47.3 cm precision. These best-case results are from the flat ice-sheet interior but provide a characterization of the quality of satellite and airborne altimetry.

K. M. Brunt↗

Anomalous snow accumulation over the southeast region of the Greenland ice sheet during 2002-2003 snow season

Our objective is to determine seasonal snow accumulation in the percolation zone of the Greenland ice sheet on the daily-weekly basis over the large scale. Our approach utilizes data from the Greenland Climate Network (GC-Net) and from the SeaWinds Scatterometer on the QuikSCAT satellite (QSCAT) to measure snow accumulation (SA) in the percolation zone of the Greenland ice sheet. GC-Net measurements provide crucial in-situ data to facilitate the interpretation of QSCAT backscatter signature for the development of an algorithm to map SA.

Greenland↗

Determination of Local Slope on the Greenland Ice Sheet Using a Multibeam Photon-Counting Lidar in Preparation for the ICESat-2 Mission

The greatest changes in elevation in Greenland and Antarctica are happening along the margins of the ice sheets where the surface frequently has significant slopes. For this reason, the upcoming Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) mission utilizes pairs of laser altimeter beams that are perpendicular to the flight direction in order to extract slope information in addition to elevation. The Multiple Altimeter Beam Experimental Lidar (MABEL) is a high-altitude airborne laser altimeter designed as a simulator for ICESat-2. The MABEL design uses multiple beams at fixed angles and allows for local slope determination. Here, we present local slopes as determined by MABEL and compare them to those determined by the Airborne Topographic Mapper (ATM) over the same flight lines in Greenland. We make these comparisons with consideration for the planned ICESat-2 beam geometry. Results indicate that the mean slope residuals between MABEL and ATM remain small (< 0.05◦) through a wide range of localized slopes using ICESat-2 beam geometry. Furthermore, when MABEL data are subsampled by a factor of 4 to mimic the planned ICESat-2 transmit-energy configuration, the results are indistinguishable from the full-data-rate analysis. Results from MABEL suggest that ICESat-2 beam geometry and transmit-energy configuration are appropriate for the determination of slope on ∼90-m spatial scales, a measurement that will be fundamental to deconvolving the effects of surface slope from the ice-sheet surface change derived from ICESat-2.

Greenland↗

Intercomparison of surface meltwater routing models for the Greenland ice sheet and influence on subglacial effective pressures

Each summer, large volumes of surface meltwater drain off the Greenland ice sheet (GrIS) surface through moulins to the bed, impacting subglacial hydrology and ice flow dynamics. Supraglacial surface routing delays may propagate to englacial and subglacial hydrologic systems, requiring accurate assessment to correctly estimate subglacial effective pressures. We compare hourly supraglacial moulin discharge simulations from three surface meltwater routing models – the synthetic unit hydrograph (SUH), the bare-ice component of surface routing and lake filling (SRLF), and the rescaled width function (RWF) – for four internally drained catchments on the southwestern Greenland ice sheet surface. The routing models are forced identically using surface runoff from the Modèle Atmosphérique Régionale regional climate model (RCM). For each catchment, simulated moulin hydrographs are input to the SHAKTI subglacial hydrologic model to simulate diurnally varying subglacial effective-pressure variations in the vicinity of a single moulin. Overall, all three routing models produce more realistic moulin discharges than simply using RCM runoff outputs without surface routing but produce significant differences in peak moulin discharge and time to peak. In particular, the RWF yields later, smaller peak moulin discharges than the SUH or SRLF due to its representation of slow interfluve flow between supraglacial meltwater channels, and it can readily accommodate the seasonal evolution of supraglacial stream and river networks. Differences among the three routing models are reflected in a series of simple idealized subglacial hydrology simulations that yield different diurnal effective-pressure amplitudes; however, the supraglacial hydrologic system acts as short-term storage for surface meltwater, and the temporal mean effective pressure is relatively consistent across routing models.

Kang Yang↗

Determination of Local Slope on the Greenland Ice Sheet Using a Multibeam Photon-Counting Lidar in Preparation for the ICESat-2 Mission

The greatest changes in elevation in Greenland and Antarctica are happening along the margins of the ice sheets where the surface frequently has significant slopes. For this reason, the upcoming Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) mission utilizes pairs of laser altimeter beams that are perpendicular to the flight direction in order to extract slope information in addition to elevation. The Multiple Altimeter Beam Experimental Lidar (MABEL) is a high-altitude airborne laser altimeter designed as a simulator for ICESat-2. The MABEL design uses multiple beams at fixed angles and allows for local slope determination. Here, we present local slopes as determined by MABEL and compare them to those determined by the Airborne Topographic Mapper (ATM) over the same flight lines in Greenland. We make these comparisons with consideration for the planned ICESat-2 beam geometry. Results indicate that the mean slope residuals between MABEL and ATM remain small (< 0.05 degrees) through a wide range of localized slopes using ICESat-2 beam geometry. Furthermore, when MABEL data are subsampled by a factor of 4 to mimic the planned ICESat-2 transmit-energy configuration, the results are indistinguishable from the full-data-rate analysis. Results from MABEL suggest that ICESat-2 beam geometry and transmit-energy configuration are appropriate for the determination of slope on approx. 90-m spatial scales, a measurement that will be fundamental to deconvolving the effects of surface slope from the ice-sheet surface change derived from ICESat-2.

Photon-Counting↗

Getting around Antarctica: New High-Resolution Mappings of the Grounded and Freely-Floating Boundaries of the Antarctic Ice Sheet Created for the International Polar Year

Two ice-dynamic transitions of the Antarctic ice sheet - the boundary of grounded ice features and the freely-floating boundary - are mapped at 15-m resolution by participants of the International Polar Year project ASAID using customized software combining Landsat-7 imagery and ICESat/GLAS laser altimetry. The grounded ice boundary is 53 610 km long; 74% abuts to floating ice shelves or outlet glaciers, 19% is adjacent to open or sea-ice covered ocean, and 7% of the boundary ice terminates on land. The freely-floating boundary, called here the hydrostatic line, is the most landward position on ice shelves that expresses the full amplitude of oscillating ocean tides. It extends 27 521 km and is discontinuous. Positional (one-sigma) accuracies of the grounded ice boundary vary an order of magnitude ranging from +/- 52m for the land and open-ocean terminating segments to +/- 502m for the outlet glaciers. The hydrostatic line is less well positioned with errors over 2 km. Elevations along each line are selected from 6 candidate digital elevation models based on their agreement with ICESat elevation values and surface shape inferred from the Landsat imagery. Elevations along the hydrostatic line are converted to ice thicknesses by applying a firn-correction factor and a flotation criterion. BEDMAP-compiled data and other airborne data are compared to the ASAID elevations and ice thicknesses to arrive at quantitative (one-sigma) uncertainties of surface elevations of +/-3.6, +/-9.6, +/-11.4, +/-30 and +/-100m for five ASAID-assigned confidence levels. Over one-half of the surface elevations along the grounded ice boundary and over one-third of the hydrostatic line elevations are ranked in the highest two confidence categories. A comparison between ASAID-calculated ice shelf thicknesses and BEDMAP-compiled data indicate a thin-ice bias of 41.2+/-71.3m for the ASAID ice thicknesses. The relationship between the seaward offset of the hydrostatic line from the grounded ice boundary only weakly matches a prediction based on beam theory. The mapped products along with the customized software to generate them and a variety of intermediate products are available from the National Snow and Ice Data Center.

Bindschadler, R.↗

A Satellite-Derived Climate-Quality Data Record of the Clear-Sky Surface Temperature of the Greenland Ice Sheet

We have developed a climate-quality data record of the clear-sky surface temperature of the Greenland Ice Sheet using the Moderate-Resolution Imaging Spectroradiometer (MODIS) ice-surface temperature (1ST) algorithm. A climate-data record (CDR) is a time series of measurements of sufficient length, consistency, and continuity to determine climate variability and change. We present daily and monthly MODIS ISTs of the Greenland Ice Sheet beginning on 1 March 2000 and continuing through 31 December 2010 at 6.25-km spatial resolution on a polar stereographic grid. This record will be elevated in status to a CDR when at least nine more years of data become available either from MODIS Terra or Aqua, or from the Visible Infrared Imager Radiometer Suite (VIIRS) to be launched in October 2011. Our ultimate goal is to develop a CDR that starts in 1981 with the Advanced Very High Resolution (AVHRR) Polar Pathfinder (APP) dataset and continues with MODIS data from 2000 to the present, and into the VIIRS era. Differences in the APP and MODIS cloud masks have so far precluded the current 1ST records from spanning both the APP and MODIS time series in a seamless manner though this will be revisited when the APP dataset has been reprocessed. The complete MODIS 1ST daily and monthly data record is available online.

Hall, Dorothy K.↗

A Satellite-Derived Climate-Quality Data Record of the Clear-Sky Surface Temperature of the Greenland Ice Sheet

We have developed a climate-quality data record of the clear-sky surface temperature of the Greenland Ice Sheet using the Moderate-Resolution Imaging Spectroradiometer (MODIS) Terra ice-surface temperature (1ST) algorithm. A climate-data record (CDR) is a time series of measurements of sufficient length, consistency, and continuity to determine climate variability and change. We present daily and monthly Terra MODIS ISTs of the Greenland Ice Sheet beginning on 1 March 2000 and continuing through 31 December 2010 at 6.25-km spatial resolution on a polar stereographic grid within +/-3 hours of 17:00Z or 2:00 PM Local Solar Time. Preliminary validation of the ISTs at Summit Camp, Greenland, during the 2008-09 winter, shows that there is a cold bias using the MODIS IST which underestimates the measured surface temperature by approximately 3 C when temperatures range from approximately -50 C to approximately -35 C. The ultimate goal is to develop a CDR that starts in 1981 with the Advanced Very High Resolution (AVHRR) Polar Pathfinder (APP) dataset and continues with MODIS data from 2000 to the present. Differences in the APP and MODIS cloud masks have so far precluded the current IST records from spanning both the APP and MODIS IST time series in a seamless manner though this will be revisited when the APP dataset has been reprocessed. The Greenland IST climate-quality data record is suitable for continuation using future Visible Infrared Imager Radiometer Suite (VIIRS) data and will be elevated in status to a CDR when at least 9 more years of climate-quality data become available either from MODIS Terra or Aqua, or from the VIIRS. The complete MODIS IST data record will be available online in the summer of 2011.

Hall, Dorothy K.↗

Free oscillations in a climate model with ice-sheet dynamics

A study of stable periodic solutions to a simple nonlinear model of the ocean-atmosphere-ice system is presented. The model has two dependent variables: ocean-atmosphere temperature and latitudinal extent of the ice cover. No explicit dependence on latitude is considered in the model. Hence all variables depend only on time and the model consists of a coupled set of nonlinear ordinary differential equations. The globally averaged ocean-atmosphere temperature in the model is governed by the radiation balance. The reflectivity to incoming solar radiation, i.e., the planetary albedo, includes separate contributions from sea ice and from continental ice sheets. The major physical mechanisms active in the model are (1) albedo-temperature feedback, (2) continental ice-sheet dynamics and (3) precipitation-rate variations. The model has three-equilibrium solutions, two of which are linearly unstable, while one is linearly stable. For some choices of parameters, the stability picture changes and sustained, finite-amplitude oscillations obtain around the previously stable equilibrium solution. The physical interpretation of these oscillations points to the possibility of internal mechanisms playing a role in glaciation cycles.

Kallen, E.↗

Anomalous Circulation in July 2019 Resulting in Mass Loss on the Greenland Ice Sheet

Current mass loss on the Greenland Ice Sheet (GrIS) includes a significant contribution from surface runoff. The circumstances associated with melt events are important for understanding the global sea level contribution of the GrIS. In late July 2019, surface melt occurred over 62% of the GrIS, including Summit Station. The general circulation leading to the event is found to be dissimilar to 2012 and other events documented in the 21st century, with warm air associated with remote atmospheric blocking over western Europe eventually transiting west to the GrIS. Gravimetric data indicate that the 2019 summer mass loss was 137 Gt more than the 2004–2010 median, or about 92% of the 2012 record. Mass loss during the event was significant in GrIS northeastern regions in 2019. As compared to 2012, the southwest did not fully participate. Similar circulation patterns have not previously been associated with significant melt.

Richard I Cullather↗

Annual Sea Level Variability Induced by Changes in Sea Ice Extent and Accumulation on Ice Sheets: An Assessment Based on Remotely Sensed Data

Changes of mean annual net accumulation at the surface on the grounded ice sheets of East Antarctica, West Antarctica, and Greenland in response to variations in sea ice extent are estimated using grid-point values 100 km apart. The data bases are assembled principally by bilinear interpolation of remotely sensed brightness temperature (Nimbus-5 ESMR, Nimbus-7 SMMR), surface temperature (Nimbus-7 THIR), and surface elevation (ERS-1 radar altimeter). These data, complemented by field data where remotely sensed data are not available, are used in multivariate analyses in which mean annual accumulation (derived from firn emissivity) is the dependent variable; the independent variables are latitude, surface elevation, mean annual surface temperature, and mean annual distance to open ocean (as a source of energy and moisture). The last is the shortest distance measured between a grid point and the mean annual position of the 10% sea ice concentration boundary, and is used as an index of changes in sea ice extent as well as of mean concentration. Stepwise correlation analyses indicate that variations in sea ice extent of +/-50 km would lead to changes in accumulation inversely of +/-4% on East Antarctica, +/- 10% on West Antarctica, and +4% on Greenland. These results are compared with those obtained in a previous study using visually interpolated values from contoured compilations of field data; they substantiate the findings for the Antarctic ice sheets (+/-4% on East Antarctica, +/-9% in West Antarctica), and suggest a reduction by one half of the probable change of accumulation on Greenland (from +/-8%). The results also suggest a reduction of the combined contribution to sea level variability to +/- 0.19 mm/a (from +/- 0.22 mm/a).

Zwally, H. J.↗