Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice sheet modeling”

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 55 records · Page 3

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↗

Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt

Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.

54 ENVIRONMENTAL SCIENCES↗

Development of a Multilayer MODIS IST-Albedo Product of Greenland

A new multilayer IST-albedo Moderate Resolution Imaging Spectroradiometer (MODIS) product of Greenland was developed to meet the needs of the ice sheet modeling community. The multiple layers of the product enable the relationship between IST and albedo to be evaluated easily. Surface temperature is a fundamental input for dynamical ice sheet models because it is a component of the ice sheet radiation budget and mass balance. Albedo influences absorption of incoming solar radiation. The daily product will combine the existing standard MODIS Collection-6 ice-surface temperature, derived melt maps, snow albedo and water vapor products. The new product is available in a polar stereographic projection in NetCDF format. The product will ultimately extend from March 2000 through the end of 2017.

MODIS↗

Uncertainties in the Modelled CO2 Threshold for Antarctic Glaciation

frequently cited atmospheric CO2 threshold for the onset of Antarctic glaciation of approximately780 parts per million by volume is based on the study of DeConto and Pollard (2003) using an ice sheet model and the GENESIS climate model. Proxy records suggest that atmospheric CO2 concentrations passed through this threshold across the Eocene-Oligocene transition approximately 34 million years. However, atmospheric CO2 concentrations may have been close to this threshold earlier than this transition, which is used by some to suggest the possibility of Antarctic ice sheets during the Eocene. Here we investigate the climate model dependency of the threshold for Antarctic glaciation by performing offline ice sheet model simulations using the climate from 7 different climate models with Eocene boundary conditions (HadCM3L, CCSM3, CESM1.0, GENESIS, FAMOUS, ECHAM5 and GISS_ER). These climate simulations are sourced from a number of independent studies, and as such the boundary conditions, which are poorly constrained during the Eocene, are not identical between simulations. The results of this study suggest that the atmospheric CO2 threshold for Antarctic glaciation is highly dependent on the climate model used and the climate model configuration. A large discrepancy between the climate model and ice sheet model grids for some simulations leads to a strong sensitivity to the lapse rate parameter.

simulation↗

Ice Flow in the North East Greenland Ice Stream

Early observations with ERS-1 SAR image data revealed a large ice stream in North East Greenland (Fahnestock 1993). The ice stream has a number of the characteristics of the more closely studied ice streams in Antarctica, including its large size and gross geometry. The onset of rapid flow close to the ice divide and the evolution of its flow pattern, however, make this ice stream unique. These features can be seen in the balance velocities for the ice stream (Joughin 1997) and its outlets. The ice stream is identifiable for more than 700 km, making it much longer than any other flow feature in Greenland. Our research goals are to gain a greater understanding of the ice flow in the northeast Greenland ice stream and its outlet glaciers in order to assess their impact on the past, present, and future mass balance of the ice sheet. We will accomplish these goals using a combination of remotely sensed data and ice sheet models. We are using satellite radar interferometry data to produce a complete maps of velocity and topography over the entire ice stream. We are in the process of developing methods to use these data in conjunction with existing ice sheet models similar to those that have been used to improve understanding of the mechanics of flow in Antarctic ice streams.

Joughin, Ian↗

Sensitivity of an energy balance climate model with predicted snowfall rates

A snowfall parameterization and a polar-ice-sheet model are developed and applied to the two-level zonally averaged seasonal energy-balance climate model of Held and Suarez (1979), and sensitivity experiments involving changes in insolation are performed both with and without ice sheets. The results are presented in tables and graphs, and the hydrological-cycle response to insolation changes is found to be similar to that predicted by global-circulation models employing prescribed precipitation levels, with a somewhat higher sensitivity in the snow line. The area covered by ice sheets in the ice-sheet models is shown to be greater than that covered by permanent snow in the models without ice sheets, an effect attributed to lower surface temperatures over the ice. It is inferred that an increase in the solar constant can cause increased high-latitude precipitation but not an ice age.

Bowman, K. P.↗

Lunar and Planetary Science XXXV: Special Session: Mars Climate Change

The session "Mars Climate Change" contained the following reports:Geological Evidence for Climate Change on Mars; A New Astronomical Solution for the Long Term Evolution of the Insolation Quantities of Mars; Interpreting Martian Paleoclimate with a Mars General Circulation Model; History and Progress of GCM Simulations on Recent Mars Climate Change; Northern and Southern Permafrost Regions on Mars with High Content of Water Ice: Similarities and Differences; Periods of Active Permafrost Layer Formation in the Recent Geological History of Mars; Microclimate Zones in the Dry Valleys of Antarctica: Implications for Landscape; Evolution and Climate Change on Mars; Geomorphic Evidence for Martian Ground Ice and Climate Change; Explaining the Mid-Latitude Ice Deposits with a General Circulation Model; Tharsis Montes Cold-based Glaciers: Observations and Constraints for Modeling and Preliminary Results; Ice Sheet Modeling: Terrestrial Background and Application to Arsia Mons Lobate Deposit, Mars; Enhanced Water-Equivalent Hydrogen on the Western Flanks of the Tharsis Montes and Olympus Mons: Remnant Subsurface Ice or Hydrate Minerals?; and New Age Mars.

Source record↗

Lunar and Planetary Science XXXV: Special Session: Mars Climate Change

The session "Mars Climate Change" included the following topics:Geological Evidence for Climate Change on Mars; A New Astronomical Solution for the Long Term Evolution of the Insolation Quantities of Mars; Interpreting Martian Paleoclimate with a Mars General Circulation Model; History and Progress of GCM Simulations on Recent Mars Climate Change; Northern and Southern Permafrost Regions on Mars with High Content of Water Ice: Similarities and Differences; Periods of Active Permafrost Layer Formation in the Recent Geological History of Mars; Microclimate Zones in the Dry Valleys of Antarctica: Implications for Landscape Evolution and Climate Change on Mars; Geomorphic Evidence for Martian Ground Ice and Climate Change; Explaining the Mid-Latitude Ice Deposits with a General Circulation Model; Tharsis Montes Cold-based Glaciers: Observations and Constraints for Modeling and Preliminary Results; Ice Sheet Modeling: Terrestrial Background and Application to Arsia Mons Lobate Deposit, Mars; Enhanced Water-Equivalent Hydrogen on the Western Flanks of the Tharsis Montes and Olympus Mons: Remnant Subsurface Ice or Hydrate Minerals?; and New Age Mars.

Source record↗

Ice Flow in the Northeast Greenland Ice Stream

As the use of satellite data for icesheet research has accelerated, new discoveries have changed the way that we think about ice discharge, and new tools have changed the way that we measure it. We have begun to apply a combination of remote sensing and ice sheet modeling techniques to study ice flow in the recently identified northeast Greenland ice stream. This ice stream shows evidence of organized flow far into the interior and has low-slope areas of rapid flow and regions of enhanced shear that resemble the ice streams of West Antarctica. The onset of rapid flow close to the ice divide and the evolution of its flow pattern, however, make this ice stream unique. Our investigation is aimed at increasing our understanding of this ice stream and its origins and role in the discharge of ice from the ice sheet and the effects that discharge may have on the history and evolution of the ice in northern Greenland. The ice flow data comes from satellite radar interferometry tied to reference points. We have applied this technique to measure the velocity over the entire ice stream, which is more than 700 km in length. The velocity contours reveal distinct, well-developed shear margins, which are coincident with linear features in the Synthetic Aperture Radar (SAR) amplitude imagery. It is surprising to see such defined, enhanced flow so far inland and so close to the ice divide. The velocity data and topography data are being used to feed a modeling effort that allows us to study the dynamics and character rapid flow. Additional information is contained in the original.

Joughin, Ian↗

The Annual Glaciohydrology Cycle in the Ablation Zone of the Greenland Ice Sheet: Part 2. Observed and Modeled Ice Flow

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 (46km 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.

ice flow↗