Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ice sheet model”

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 145 records · Page 8

Backscatter model for the unusual radar properties of the Greenland Ice Sheet

A number of planetary objects exhibit large radar reflectivity and polarization ratios, and more recently, a similar behavior has been observed over a vast portion of the Earth's surface: the percolation facies of the Greenland Ice Sheet. Surface-based ranging radar data and snow stratigraphy studies demonstrated that the radar properties of that portion of Greenland are caused by enhanced scattering from massive, large, solid-ice bodies buried in the top few meters of the dry, cold, clean snowy surface of the ice sheet and created by seasonal melting and refreezing events. Here, we model the icy inclusions as randomly oriented, discrete, noninteracting, dielectric cylinders embedded in a transparent snow medium. An exact analytical solution is used to compute the scattered field from the cylinders. Using this model, we correctly predict the polarimetric radar observations gathered by an airborne imaging system at three wavelengths (5.6, 24, and 68 cm), between 19 deg and 65 deg incidence angle. The diameter and number density of the cylinders that are inferred from the radar data using the model are consistent with in situ observations of the icy inclusions. The large radar reflectivity and polarization ratios are interpreted as arising internal reflections of the radar signals in the icy inclusions that first-order external reflection models fail to predict. The results compare favorably with predictions from the coherent backscatter or weak localization theory and may provide a complementary framework for interpreting exotic radar echoes from other planetary objects.

Rignot, Eric↗

Spatial patterns in backscatter strength across the Greenland Ice Sheet

The relationship between the physical properties of the Greenland ice sheet and Synthetic Aperture Radar (SAR) data collected from aircraft and from ERS-1 is addressed. Limited aircraft data are combined with a description of the glacier surface to predict qualitatively the spatial and seasonal variation in backscatter strength across the ice sheet. In particular the model predicts relatively low backscatter near the ice edge where scattering is dominated by rough surface effects. Backscatter increases through the lake zone as volume scattering becomes important. Strongest backscatter is found in the percolation facies where volume scatter from snow grains and volume scatter from large, buried ice bodies becomes important. Backscatter weakens in the interior ice sheet where fine grained snow is the only mechanism producing backscatter.

Jezek, K. C.↗

How Much, How Fast?: A Science Review and Outlook for Research on the Instability of Antarctica's Thwaites Glacier in the 21st Century

Constraining how much and how fast the West Antarctic Ice Sheet (WAIS) will change in the coming decades has recently been identified as the highest priority in Antarctic research (National Academies, 2015). Here we review recent research on WAIS and outline further scientific objectives for the area now identified as the most likely to undergo near-term significant change: Thwaites Glacier and the adjacent Amundsen Sea. Multiple lines of evidence point to an ongoing rapid loss of ice in this region in response to changing atmospheric and oceanic conditions. Models of the ice sheets dynamic behavior indicate a potential for greatly accelerated ice loss as ocean-driven melting at the Thwaites Glacier grounding zone and nearby areas leads to thinning, faster flow, and retreat. A complete retreat of the Thwaites Glacier basin would raise global sea level by more than three meters by entraining ice from adjacent catchments. This scenario could occur over the next few centuries, and faster ice loss could occur through processes omitted from most ice flow models such as hydrofracture and ice cliff failure, which have been observed in recent rapid ice retreats elsewhere. Increased basal melt at the grounding zone and increased potential for hydrofracture due to enhanced surface melt could initiate a more rapid collapse of Thwaites Glacier within the next few decades.

Thwaites Glacier↗

The diatom record from beneath the West Antarctic Ice Sheet and the global proxy perspective

Recent glaciological evaluation and modeling of the marine-based West Antarctic Ice Sheet (WAIS) support the possibility that the WAIS disintegrated during one or more Pleistocene interglacial period(s). The magnitude of sea level and oxygen isotope variation during certain late-Pleistocene interglacial periods is also consistent with the possibility of major retreat of the WAIS. Although oxygen isotopes from deep-sea sediments provide the best available proxy record for global ice volume (despite the ambiguities in the record), the source of ice volume changes must be hypothesized. Based on the intensity of interglacial isotopic shifts recorded in Southern Ocean marine sedimentary records, stage 11 (400,000 years ago) is the strongest candidate for WAIS collapse, but the records for stages 9, 7, and 5.5 are all consistent with the possibility of multiple late-Pleistocene collapses. Seismic reflection studies through the WAIS have revealed thick successions of strata with seismic characteristics comparable to upper Tertiary marine sediments. Small samples of glacial diamictons from beneath the ice sheet have been collected via hot-water drilled access holes. These sediments include mixed diatom assemblages of varying ages. Late-Miocene diatoms dominate many samples, probably reflecting marine deposition in West Antarctic basins prior to development of a dominantly glacial phase in West Antarctica. In addition to late-Miocene diatoms, samples from Upstream B (1988/89) contain rare post-Miocene diatoms, many of which imply deposition in the West Antarctic interior during one or more Pleistocene deglaciation periods. Age-diagnostic fossils in glacial sediments beneath ice sheets provide relatively coarse chronostratigraphic control, but they do contain direct evidence of regional deglaciation. Thus, sub-glacial till samples provide the evidence regarding the source of ice sheet variability seen in well-dated proxy records. Combined, these independent data sets can provide a more comprehensive and less speculative interpretation of the history of past glacial minima in currently glaciated polar regions.

Scherer, Reed P.↗

Simulations of Firn Processes Over the Greenland and Antarctic Ice Sheets: 1980-2021

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 modeled. Here, we present new simulations of firn processes over the Greenland and Antarctic ice sheets (GrIS and AIS) using the Community Firn Model and atmospheric reanalysis variables for more than four decades. A data set 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 1980–2020 modeled firn column runoff, when combined with atmospheric variables from MERRA-2, generates realistic mean integrated surface mass balance values for the Greenland (+390 Gt yr(exp −1)) and Antarctic (+2612 Gt yr(exp −1)) ice sheets when compared to published model-ensemble means. We find that seasonal volume changes associated with firn air content are on average approximately 2.5 times larger than those associated with mass fluxes from surface processes for the AIS and 1.5 times larger for the GrIS; 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 nearly 5% of its firn air content, indicating a reduction in the total meltwater retention capability. Nearly all (94 %) of the meltwater produced over the Antarctic Ice Sheet is retained within the firn column through infiltration and refreezing.

Greenland↗

A computationally efficient statistically downscaled 100 m resolution Greenland product from the regional climate model MAR

The Greenland Ice Sheet (GrIS) has been contributing directly to sea level rise, and this contribution is projected to accelerate over the next decades. A crucial tool for studying the evolution of surface mass loss (e.g., surface mass balance, SMB) consists of regional climate models (RCMs), which can provide current estimates and future projections of sea level rise associated with such losses. However, one of the main limitations of RCMs is the relatively coarse horizontal spatial resolution at which outputs are currently generated. Here, we report results concerning the statistical downscaling of the SMB modeled by the Modèle Atmosphérique Régional (MAR) RCM from the original spatial resolution of 6 km to 100 m building on the relationship between elevation and mass losses in Greenland. To this goal, we developed a geospatial framework that allows the parallelization of the downscaling process, a crucial aspect to increase the computational efficiency of the algorithm. Using the results obtained in the case of the SMB, surface and air temperature are assessed through the comparison of the modeled outputs with in situ and satellite measurement. The downscaled products show a considerable improvement in the case of the downscaled product with respect to the original coarse output, with the coefficient of determination (R 2 ) increasing from 0.868 for the original MAR output to 0.935 for the SMB downscaled product. Moreover, the value of the slope and intercept of the linear regression fitting modeled and measured SMB values shifts from 0.865 for the original MAR to 1.015 for the downscaled product in the case of the slope and from the value −235 mm w.e. yr -1 (original) to −57 mm w.e. yr -1 (downscaled) in the case of the intercept, considerably improving upon results previously published in the literature.

Greenland Ice Sheet↗

Measuring Geophysical Parameters of the Greenland Ice Sheet using Airborne Radar Altimetry

This paper presents radar-altimeter scattering models for each of the diagenetic zones of the Greenland ice sheet. AAFE radar- altimeter waveforms obtained during the 1991 and 1993 NASA multi-sensor airborne altimetry experiments over Greenland reveal that the Ku-band return pulse changes significantly with the different diagenetic zones. These changes are due to varying amounts of surface and volume scattering in the return waveform. In the ablation and soaked zones, where surface scattering dominates the AAFE return, geophysical parameters such as rms surface height and rms surface slope are obtained by fitting the waveforms to a surface-scattering model. Waveforms from the percolation zone show that the sub-surface ice features have a much more significant effect on the return pulse than the surrounding snowpack. Model percolation waveforms, created using a combined surface- and volume-scattering model and an ice-feature distribution obtained during the 1993 field season, agree well with actual AAFE waveforms taken in the same time period. Using a combined surface- and volume-scattering model for the dry-snow-zone return waveforms, the rms surface height and slope and the attenuation coefficient of the snowpack are obtained. These scattering models not only allow geophysical parameters of the ice sheet to be measured but also help in the understanding of satellite radar-altimeter data.

Ferraro, Ellen J.↗

Obliquity, Ice Sheets, and Layered Sediments on Mars: What Spacecraft Observations and Climate Models are Telling Us

The Mars Odyssey Gamma-Ray Spectrometer (GRS) data present a quandary: On the one hand, large deposits of (inferred) water ice are located where thermal models suggest they would form and best be protected, e.g., if deposited during periods of higher obliquity. On the other hand, the volume mixing ratios (approx. 70%) are so high that diffusive deposition of water in regolith pore space (which is the process assumed by these models) cannot be the primary formation mechanism. Furthermore, given that the water is inferred to be so close to the surface (less than a few 10's of cm's), it must be in communication with the atmosphere on time scales that are geologically relatively short (10(exp 3)-10(exp 6) years); therefore the water cannot be archaic. Considering the GRS data, images of mantled, fretted, and disaggregated terrain, and new climate modeling of Mars orbital cycles, we are led to an alternate conclusion about the ice deposits: that they form as subaerial ice sheets. This scenario not only provides a simple explanation for these observations, but may also help explain the formation of globally distributed, sedimentary layered deposits.

Richardson, M. I.↗

Multi-model Comparison of the Volcanic Sulfate Deposition from the 1815 Eruption of Mt. Tambora

The eruption of Mt. Tambora in 1815 was the largest volcanic eruption of the past 500 years. The eruption had significant climatic impacts, leading to the 1816 "year without a summer", and remains a valuable event from which to understand the climatic effects of large stratospheric volcanic sulfur dioxide injections. The eruption also resulted in one of the strongest and most easily identifiable volcanic sulfate signals in polar ice cores, which are widely used to reconstruct the timing and atmospheric sulfate loading of past eruptions. As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), five state-of-the-art global aerosol models simulated this eruption. We analyze both simulated background (no Tambora) and volcanic (with Tambora) sulfate deposition to polar regions and compare to ice core records. The models simulate overall similar patterns of background sulfate deposition, although there are differences in regional details and magnitude. However, the volcanic sulfate deposition varies considerably between the models with differences in timing, spatial pattern and magnitude. Mean simulated deposited sulfate on Antarctica ranges from 19 to 264 kgkm-2 and on Greenland from 31 to 194 kgkm-2, as compared to the mean ice-core derived estimates of roughly 50 kgkm-2 for both Greenland and Antarctica. The ratio of the hemispheric atmospheric sulfate aerosol burden after the eruption to the average ice sheet deposited sulfate varies between models by up to a factor of 15. Sources of this inter-model variability include differences in both the formation and the transport of sulfate aerosol. Our results suggest that deriving relationships between sulfate deposited on ice sheets and atmospheric sulfate burdens from model simulations may be associated with greater uncertainties than previously thought.

Marshall, Lauren↗

New York City Panel on Climate Change 2019 Report Chapter 3: Sea Level Rise

The New York City Panel on Climate Change (NPCC, 2015) sea level rise projections provide the current scientific basis for New York City scientific decision making and planning, as reflected in, for example, the City's Climate Resiliency Design Guidelines. However, since the IPCC (2013) and NPCC (2015) reports, recent observations show mounting glacier and ice sheet losses leading to rising sea levels. Furthermore, new developments in modeling interactions between oceans, atmosphere, and ice sheets suggest the possibility of a significantly higher global mean sea level rise (GMSLR) by 2100 than previously anticipated, particularly under elevated greenhouse gas emission scenarios.

Gornitz, Vivien M.↗

A combined surface/volume scattering retracking algorithm for ice sheet satellite altimetry

An algorithm that is based upon a combined surface-volume scattering model is developed. It can be used to retrack individual altimeter waveforms over ice sheets. An iterative least-squares procedure is used to fit the combined model to the return waveforms. The retracking algorithm comprises two distinct sections. The first generates initial model parameter estimates from a filtered altimeter waveform. The second uses the initial estimates, the theoretical model, and the waveform data to generate corrected parameter estimates. This retracking algorithm can be used to assess the accuracy of elevations produced from current retracking algorithms when subsurface volume scattering is present. This is extremely important so that repeated altimeter elevation measurements can be used to accurately detect changes in the mass balance of the ice sheets. By analyzing the distribution of the model parameters over large portions of the ice sheet, regional and seasonal variations in the near-surface properties of the snowpack can be quantified.

Davis, Curt H.↗

Temperature and Snowfall in Western Queen Maud Land Increasing Faster than Climate Model Projections

East Antarctic Ice Sheet (EAIS) mass balance is largely driven by snowfall. Recently, increased snowfall in Queen Maud Land led to years of EAIS mass gain. It is difficult to determine whether these years of enhanced snowfall are anomalous or part of a longer-term trend, reducing our ability to assess the mitigating impact of snowfall on sea level rise. We determine that the recent snowfall increases in western Queen Maud Land (QML) are part of a long-term trend (+5.2 +/- 3.7% decade(exp -1)) and are unprecedented over the past two millennia. Warming between 1998 and 2016 is significant and rapid (+1.1 +/- 0.7 C decade(exp -1)). Using these observations, we determine that the current accumulation and temperature increases in QML from an ensemble of global climate simulations are too low, which suggests that projections of the QML contribution to sea level rise are potentially overestimated with a reduced mitigating impact of enhanced snowfall in a warming world.

WAIS↗

Evaluating Greenland Surface-Mass-Balance and Firn-Densification Data Using ICESat-2 Altimetry

Surface-mass-balance (SMB) and firn-densification (FD) models are widely used in altimetry studies as a tool to separate atmospheric-driven from ice-dynamics-driven ice-sheet mass changes and to partition observed volume changes into ice-mass changes and firn-air-content changes. Until now, SMB models have been principally validated based on comparison with ice core and weather station data or comparison with widely separated flight radar-survey flight lines. Firn-densification models have been primarily validated based on their ability to match net densification over decades, as recorded in firn cores, and the short-term time-dependent component of densification has rarely been evaluated at all. The advent of systematic ice-sheet-wide repeated ice-surface-height measurements from ICESat-2 (the Ice Cloud, and land Elevation Satellite, 2) allows us to measure the net surface-height change of the Greenland ice sheet at quarterly resolution and compare the measured surface-height differences directly with those predicted by three FD–SMB models: MARv3.5.11 (Modèle Atmosphérique Régional version 3.5.11) and GSFCv1.1 and GSFCv1.2 (the Goddard Space Flight Center FD–SMB models version 1.1 and 1.2). By segregating the data by season and elevation, and based on the timing and magnitude of modelled processes in areas where we expect minimal ice-dynamics-driven height changes, we investigate the models' accuracy in predicting atmospherically driven height changes. We find that while all three models do well in predicting the large seasonal changes in the low-elevation parts of the ice sheet where melt rates are highest, two of the models (MARv3.5.11 and GSFCv1.1) systematically overpredict, by around a factor of 2, the magnitude of height changes in the high-elevation parts of the ice sheet, particularly those associated with melt events. This overprediction seems to be associated with the melt sensitivity of the models in the high-elevation part of the ice sheet. The third model, GSFCv1.2, which has an updated high-elevation melt parameterization, avoids this overprediction.

Benjamin E. Smith↗

Topographic Correction of Geothermal Heat Flux in Greenland and Antarctica

Earth’s thick, polar ice sheets insulate the bedrock beneath them from the colder temperatures at the ice surface. Consequently, the small amount of geothermal heat that Earth’s bedrock releases can have a disproportionate role in controlling ice flow. Geothermal heat flux measurements beneath ice sheets are sparse, and models do not presently account for variable bed topography. Here, we present a simple statistical method for correcting geothermal heat flux models to make them consistent with known subglacial bed topography. This method is based on previous geothermal heat flux measurements across valleys and ridges, and it helps us explore the broader significance of geothermal heat flux variability beneath the Greenland and Antarctic ice sheets. Our method can also be used to estimate the influence of topography on geothermal heat flux in non-ice-covered areas.

William Colgan↗

Slope-induced errors in radar altimetry over continental ice sheets

Altimeter measurements of ranges to modelled irregular surfaces are simulated and two correction schemes are used to reconstruct the modelled surfaces from the simulated data. When the cross-track slope is negligible, the problem can be treated as two-dimensional and the relocation method is preferred to the slope correction method. In a two-dimensional test case, 85 percent of the slope-induced rms error was removed by the relocation method. The mean error along profiles of about 75 km or longer is usually reduced more than the rms error. An alternative slope correction scheme, which uses the local slope to calculate the expected error, is less effective in the two-dimensional case. Over a simulated three-dimensional surface, where groundtracks are widely spaced and cross-track slopes are significant, the slope-correction method must be used in at least the cross-track direction.

Brenner, A. C.↗