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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The effect of sub-surface volume scattering on the accuracy of ice-sheet altimeter retracking algorithms

The NASA and ESA retracking algorithms are compared with an algorithm based upon a combined surface and volume (S/V) scattering model. First, the S/V, NASA, and ESA algorithms were used to retrack over 400,000 altimeter return waveforms from the Greenland and Antarctic ice sheets. The surface elevations from the S/V algorithm were compared with the elevations produced by the NASA and ESA algorithms to determine the relative accuracy of these algorithms when subsurface volume-scattering occurs. The results show that the NASA algorithm produced surface elevations within 35 to 50 cm of the S/V algorithm, while the performance of the ESA algorithm was slightly worse. Next, by analyzing several thousand satellite crossover points from the Antarctic data set, we determined the retracking algorithm that produced the most repeatable surface elevations. The elevations derived from the S/V algorithm had the smallest RMS error for the region of the East Antarctic plateau examined here. The ESA algorithm produced erroneous estimates of elevation change when seasonal variations were present; it measured 0.7 to 1.6-m change in elevation over a 6-month period on the East Antarctic plateau where accumulation rates are only 10 cm/year.

Davis, Curt H.↗
NASA NTRS2011

A New Global Mascon Solution Tuned for High-Latitude Ice Studies

A new global mascon solution has been developed with I-arc-degree spatial and IO-day temporal sampling. The global mas cons are estimated from the reduction of nearly 8 years of GRACE K-band range-rate data. Temporal and anisotropic spatial constraints have been applied for land, ocean and ice regions. The solution construction and tuning is focused towards the Greenland and Antarctic ice sheets (GIS and AIS) as well as the Gulf of Alaska mountain glaciers (GoA). Details of the solution development will be discussed, including the mascon parameter definitions, constraints, and the tuning of the constraint damping factor. Results will be presented, exploring the spatial and temporal variability of the ice sheets and GoA regions. A detailed error analysis will be discussed, including solution dependence on iteration, damping factor, forward modeling, and multitechnique comparisons. We also investigate the fundamental resolution of the solution and the spatial correlation of ice sheet inter-annual change. Finally, we discuss future improvements, including specific constraint application for the rest of the major land ice regions and improvements in solution regularization.

Luthcke, S. B.↗
NASA NTRS2014

Assessing the Impact of Laurentide Ice-sheet Topography on Glacial Climate

Simulations of past climates require altered boundary conditions to account for known shifts in the Earth system. For the Last Glacial Maximum (LGM) and subsequent deglaciation, the existence of large Northern Hemisphere ice sheets caused profound changes in surface topography and albedo. While ice-sheet extent is fairly well known, numerous conflicting reconstructions of ice-sheet topography suggest that precision in this boundary condition is lacking. Here we use a high-resolution and oxygen-isotopeenabled fully coupled global circulation model (GCM) (GISS ModelE2-R), along with two different reconstructions of the Laurentide Ice Sheet (LIS) that provide maximum and minimum estimates of LIS elevation, to assess the range of climate variability in response to uncertainty in this boundary condition.We present this comparison at two equilibrium time slices: the LGM, when differences in ice-sheet topography are maximized, and 14 ka, when differences in maximum ice-sheet height are smaller but still exist. Overall, we find significant differences in the climate response to LIS topography, with the larger LIS resulting in enhanced Atlantic Meridional Overturning Circulation and warmer surface air temperatures, particularly over northeastern Asia and the North Pacific. These up- and downstream effects are associated with differences in the development of planetary waves in the upper atmosphere, with the larger LIS resulting in a weaker trough over northeastern Asia that leads to the warmer temperatures and decreased albedo from snow and sea-ice cover. Differences between the 14 ka simulations are similar in spatial extent but smaller in magnitude, suggesting that climate is responding primarily to the larger difference in maximum LIS elevation in the LGM simulations. These results suggest that such uncertainty in ice-sheet boundary conditions alone may significantly impact the results of paleoclimate simulations and their ability to successfully simulate past climates, with implications for estimating climate sensitivity to greenhouse gas forcing utilizing past climate states.

glaciology↗
NASA NTRS1995

The Effect of Antarctic Ice Mass Changes on Crustal Motion and Global Geodetic Observables

Glaciological estimates of the present-day ice mass balance of Antarctica vary widely, indicating the need for additional data to constrain mass-balance models. For example, recent studies find both a positive and a negative mass balance of the Antarctic ice sheet. Analysis of studies suggest that observations of present-day crustal motion, as obtained from a GPS survey could assist in defining models.

glaciology ice mass changes crustal motion geodeti↗
DOE OSTI2024

A high-resolution pseudo-polygon discrete element model for regional sea ice

Here, this work presents a pseudo-polygon discrete element model for high-resolution sea ice simulations. A scale-invariant bonded particle contact model is proposed to model joints between sea ice floes based on the smeared fracture model and a lattice spring beam model, and the Mohr–Coulomb failure criterion is implemented to represent the shearing failure mechanism of sea ice packings under complex loadings. All mechanical parameters of the bond model can be directly determined from laboratory tests. Validations of the proposed model are made by investigations of mechanical response and failure criteria of field sea ice sheets. Compared with the field observations of sea ice from satellite radar and in situ stress sensors, the proposed model is capable of reproducing the typical constitutive behavior and the Coulomb friction envelope of field sea ice. Finally, the proposed discrete element sea ice model is used to study the effect of loading rates on mechanical behavior including failure strength of regional sea ice.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2019

Improved Earth Oblateness Rate Reveals Increased Ice Sheet Losses and Mass-Driven Sea Level Rise

Satellite laser ranging (SLR) observations are routinely applied toward the estimation of dynamic oblateness, C(sub 20), which is the largest globally integrated component of Earth's time-variable gravity field. Since 2002, GRACE and GRACE Follow-On have revolutionized the recovery of higher spatial resolution features of global time-variable gravity, with SLR continuing to provide the most reliable estimates of C (sub 20).We quantify the effect of various SLR processing strategies on estimating C(sub 20) and demonstrate better signal recovery with the inclusion of GRACE-derived low-degree gravity information in the forward model. This improved SLR product modifies the Antarctic and Greenland Ice Sheet mass trends by -15.4 and -3.5 Gt/year, respectively, as compared to CSR TN11, and improves global mean sea level budget closure by modifying sea level rise by +0.08 mm/year. We recommend that this new C(sub 20) product be applied to RL06 GRACE data products for enhanced accuracy and scientific interpretation.

Loomis, B. D.↗
NASA NTRS1990

Comparison of measurements and theory for backscatter from bare and snow-covered saline ice

C-band radar backscatter measurements were made on artificially grown sea ice during the winters of 1987-1988 and 1988-1989. These measurements were made on smooth, rough, and snow-covered saline ice. The measured sigma-deg(theta) of smooth saline ice (rms height less than 0.05 cm) disagreed with small perturbation method (SPM) surface scattering predictions. Using physical parameters of the ice in a simple layer model, it us shown that this discrepancy can be explained by scattering from beneath the surface. A thin (7-cm) dry snow cover had a significant influence on backscatter from the smooth ice sheet. This influence was due to scattering from particles within the snow, and can be predicted by a commonly used empirical layer model for snow. The results of backscatter measurements of a moderately rough saline ice sheet were found to agree with SPM predictions.

Bredow, Jonathan W.↗
NASA NTRS2005

Verification, Validation, and Solution Quality in Computational Physics: CFD Methods Applied to Ice Sheet Physics

Procedures and methods for veri.cation of coding algebra and for validations of models and calculations used in the aerospace computational fluid dynamics (CFD) community would be ef.cacious if used by the glacier dynamics modeling community. This paper presents some of those methods, and how they might be applied to uncertainty management supporting code veri.cation and model validation for glacier dynamics. The similarities and differences between their use in CFD analysis and the proposed application of these methods to glacier modeling are discussed. After establishing sources of uncertainty and methods for code veri.cation, the paper looks at a representative sampling of veri.cation and validation efforts that are underway in the glacier modeling community, and establishes a context for these within an overall solution quality assessment. Finally, a vision of a new information architecture and interactive scienti.c interface is introduced and advocated.

Thompson, David E.↗
DOE OSTI2025

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗
NASA NTRS2019

A New Regional Climate Model for POLAR‐CORDEX: Evaluation of a 30‐Year Hindcast with COSMO‐CLM^2 Over Antarctica

Continent‐wide climate information over the Antarctic Ice Sheet (AIS) is important to obtain accurate information of present climate and reduce uncertainties of the ice sheet mass balance response and resulting global sea level rise to future climate change. In this study, the COSMO‐CLM^2 Regional Climate Model is applied over the AIS and adapted for the specific meteorological and climatological conditions of the region. A 30‐year hindcast was performed and evaluated against observational records consisting of long‐term ground‐based meteorological observations, automatic weather stations, radiosoundings, satellite records, stake measurements and ice cores. Reasonable agreement regarding the surface and upper‐air climate is achieved by the COSMO‐CLM^2 model, comparable to the performance of other state‐of‐the‐art climate models over the AIS. Meteorological variability of the surface climate is adequately simulated, and biases in the radiation and surface mass balance are small. The presented model therefore contributes as a new member to the COordinated Regional Downscaling EXperiment project over the AIS (POLAR‐CORDEX) and the CORDEX‐CORE initiative.

Antarctica↗
NASA NTRS2023

A New Model of Dry Firn Densification Constrained By Continuous Strain Measurements Near South Pole

Converting measurements of ice-sheet surface elevation change to mass change requires measurements of accumulation and knowledge of the evolution of the density profile in the firn. Most firn-densification models are tuned using measured depth-density profiles, a method which is based on an assumption that the density profile in the firn is invariant through time. Here we present continuous measurements of firn-compaction rates in 12 boreholes near the South Pole over a two-year period. To our knowledge, these are the first continuous measurements of firn compaction on the Antarctic Plateau. We use the data to derive a new firn-densification algorithm framed as a constitutive relationship. We also compare our measurements to compaction rates predicted by several existing firn densification models. Results indicate that an activation energy of 60 kJ mol (−1), a value within the range used by current models, best predicts the seasonal cycle in compaction rates on the Antarctic Plateau. Our results suggest models can predict firn-compaction rates with at best 7% uncertainty and cumulative firn compaction on a two year timescale with at best 8% uncertainty.

C. Max Stevens↗
NASA NTRS1995

Polarimetric scattering from layered media with multiple species of scatterers

Geophysical media are usually heterogeneous and contain multiple species of scatterers. In this paper a model is presented to calculate effective permittivities and polarimetric backscattering coefficients of multispecies-layered media. The same physical description is consistently used in the derivation of both permittivities and scattering coefficients. The strong permittivity fluctuation theory is extended to account for the multiple species of scatterers with a general ellipsoidal shape whose orientations are randomly distributed. Under the distorted Born approximation, polarimetric scattering coefficients are obtained. These calculations are applicable to the special cases of spheroidal and spherical scatterers. The model is used to study effects of scatterer shapes and multispecies mixtures on polarimetric signatures of heterogeneous media. The multispecies model accounts for moisture content in scattering media such as snowpack in an ice sheet. The results indicate a high sensitivity of backscatter to moisture with a stronger dependence for drier snow and ice grain size is important to the backscatter. For frost-covered saline ice, model results for bare ice are compared with measured data at C band and then the frost flower formation is simulated with a layer of fanlike ice crystals including brine infiltration over a rough interface. The results with the frost cover suggest a significant increase in scattering coefficients and a polarimetric signature closer to isotropic characteristics compared to the thin saline ice case.

Nghiem, S. V.↗
NASA NTRS1992

Modeling orbital changes on tectonic time scales

Geologic time series indicate significant 100 ka and 400 ka pre-Pleistocene climate fluctuations, prior to the time of such fluctuations in Pleistocene ice sheets. The origin of these fluctuations must therefore depend on phenomena other than the ice sheets. In a previous set of experiments, we tested the sensitivity of an energy balance model to orbital insolation forcing, specifically focusing on the filtering effect of the Earth's geography. We found that in equatorial areas, the twice-yearly passage of the sun across the equator interacts with the precession index to generate 100 ka and 400 ka power in our modeled time series. The effect is proportional to the magnitude of land in equatorial regions. We suggest that such changes may reflect monsoonal variations in the real climate system, and the subsequent wind and weathering changes may transfer some of this signal to the marine record. A comparison with observed fluctuations of Triassic lake levels is quite favorable. A number of problems remain to be studied or clarified: (1) the EBM experiments need to be followed up by a limited number of GCM experiments; (2) the sensitivity to secular changes in orbital forcing needs to be examined; (3) the possible modifying role of sedimentary processes on geologic time series warrants considerably more study; (4) the effect of tectonic changes on Earth's rotation rate needs to be studied; and (5) astronomers need to make explicit which of their predictions are robust and geologists and astronomers have to agree on which of the predictions are most testable in the geologic record.

Crowley, Thomas J.↗
NASA NTRS2024

A Physics-Based Antarctic Melt Detection Technique: Combining Advanced Microwave Scanning Radiometer 2, Radiative-Transfer Modeling, and Firn Modeling

Surface melt on ice shelves has been linked to hydrofracture and subsequent ice shelf breakup. Since the 1990s, scientists have been using microwave radiometers to detect melt on ice shelves and ice sheets by applying various statistical thresholding techniques to identify significant increases in brightness temperature that are associated with melt. In this study, instead of using a fixed threshold, we force the Snow Radiative Transfer model (SMRT) with outputs from the Community Firn Model (CFM) to create a dynamic, physics-based threshold for melt. In the process, we also combine our method with statistical thresholding techniques and produce microwave grain size information in the process. We run this “Hybrid Method” across the Larsen C ice shelf as well as thirteen sites on the Antarctic Ice Sheet. Melt and non-melt days from the Hybrid Method and three statistical thresholding techniques match with surface energy balance within 94±1%; the effect of melt on the passive microwave is mostly binary and thus largely detectable by statistical thresholding techniques as well as physics-based techniques. Rather than always replacing statistical thresholding techniques with the Hybrid Method, we recommend using the Hybrid Method in studies where melt volume or grain size is of interest. In this study, we show that the Hybrid Method can be used to (a) model dry snow brightness temperatures of Antarctic snow and (b) derive a measure of grain size; therefore, it is an important step forwards towards using firn and radiative transfer modeling to quantify melt rather than to simply detect melt days.

Marissa E Dattler↗
DOE OSTI2024

A multi-model assessment of the early last deglaciation (PMIP4 LDv1): a meltwater perspective

Abstract. The last deglaciation (∼20–11 ka BP) is a period of a major, long-term climate transition from a glacial to interglacial state that features multiple centennial- to decadal-scale abrupt climate variations whose root cause is still not fully understood. To better understand this time period, the Paleoclimate Modelling Intercomparison Project (PMIP) has provided a framework for an internationally coordinated endeavour in simulating the last deglaciation whilst encompassing a broad range of models. Here, we present a multi-model intercomparison of 17 transient simulations of the early part of the last deglaciation (∼20–15 ka BP) from nine different climate models spanning a range of model complexities and uncertain boundary conditions and forcings. The numerous simulations available provide the opportunity to better understand the chain of events and mechanisms of climate changes between 20 and 15 ka BP and our collective ability to simulate them. We conclude that the amount of freshwater forcing and whether it follows the ice sheet reconstruction or induces an inferred Atlantic meridional overturning circulation (AMOC) history, heavily impacts the deglacial climate evolution for each simulation rather than differences in the model physics. The course of the deglaciation is consistent between simulations except when the freshwater forcing is above 0.1 Sv – at least 70 % of the simulations agree that there is warming by 15 ka BP in most places excluding the location of meltwater input. For simulations with freshwater forcings that exceed 0.1 Sv from 18 ka BP, warming is delayed in the North Atlantic and surface air temperature correlations with AMOC strength are much higher. However, we find that the state of the AMOC coming out of the Last Glacial Maximum (LGM) also plays a key role in the AMOC sensitivity to model forcings. In addition, we show that the response of each model to the chosen meltwater scenario depends largely on the sensitivity of the model to the freshwater forcing and other aspects of the experimental design (e.g. CO2 forcing or ice sheet reconstruction). The results provide insight into the ability of our models to simulate the first part of the deglaciation and how choices between uncertain boundary conditions and forcings, with a focus on freshwater fluxes, can impact model outputs. We can use these findings as helpful insight in the design of future simulations of this time period.

Snoll, Brooke (ORCID:0000000230914394)↗
NASA NTRS2002

Comparative Results of Using Different Methods for Discovery of Microorganisms in very Ancient Layers of the Central Antarctic Glacier above the Lake Vostok

The ice sheet of the Central Antarctic is considered by the scientific community worldwide, as a model to elaborate on different methods to search for life outside Earth. This became especially significant in connection with the discovery of the underglacial lake in the vicinity of the Russian Antarctic Station Vostok. Lake Vostok is considered by many scientists as an analog of the ice covered seas of Jupiter's satellite Europa. According to the opinion of many researchers there is the possibility that relict forms of microorganisms, well preserved since the Ice Age, may be present in this lake. Investigations throughout the thickness of the ice sheet above Lake Vostok show the presence of microorganisms belonging to different well-known taxonomic groups, even in the very ancient horizons near close to floor of the glacier. Different methods were used to search for microorganisms that are rarely found in the deep ancient layers of an ice sheet. The method of aseptic sampling from the ice cores and the results of controlled sterile conditions in all stages when conducting these investigations, are described in detail in previous reports. Primary investigations tried the usual methods of sowing samples onto different nutrient media, and the result was that only a few microorganisms grew on the media used. The possibility of isolating the organisms obtained for further investigations, by using modern methods including DNA-analysis, appears to be the preferred method. Further investigations of the very ancient layers of the ice sheet by radioisotopic, luminescence, and scanning electron microscopy methods at different modifications, revealed the quantity and morphological diversity of the cells of microorganisms that were distributed on the different horizons. Investigations over many years have shown that the microflora in the very ancient strata of the Antarctic ice cover, nearest to the bedrock, support the effectiveness of using a combination of different methods to search for signs of life in ancient icy formations, which might play a role in the long-term preservation and transportation of microbial life throughout the Universe.

Abyzov, S. S.↗
NASA NTRS2020

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↗