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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Evaluating a Regional Climate Model Simulation of Greenland Ice Sheet Snow and Firn Density for Improved Surface Mass Balance Estimates

Modeling vertical profiles of snow and firn density near the surface of the Greenland ice sheet (GrIS) is key to estimating GrIS mass balance, and by extension, global sea level change. To understand sources of error in simulated GrIS density, we compare GrIS density profiles from a leading regional climate model with coincident in situ measurements. We identify key contributors to model density and mass balance biases, including underestimated simulated fresh snow density (which leads to underestimation of density in the top 1 m of snow by ~10%). In areas undergoing frequent melting, positive density biases (of 7% in the top 1 m, and 10% between 1 and 10 m) are likely associated with errors in representing meltwater production, retention and refreezing. The results highlight the importance of accurately capturing fresh snow density and meltwater processes in models used to estimate GrIS mass balance change.

P M Alexander↗

NASA NTRS2020

Enhancing New York City's resilience to sea level rise and increased coastal flooding

Accelerating Greenland and Antarctic Ice Sheet ice mass losses and potential West Antarctic Ice Sheet instability may lead to higher than previously anticipated future sea levels. The New York City Panel on Climate Change Antarctic Rapid Ice Melt (ARIM) upper-end, low probability sea level rise (SLR) scenario, which incorporates recent ice loss trends, improved ice sheet-ocean-atmosphere modeling, and potential ice sheet destabilization, projects SLR of up to 2.1 m by the 2080s and up to 2.9 m by 2100, at high greenhouse gas emissions (NPCC, 2019). These results exceed previous high-end SLR projections (90th percentile) of 1.5 m by the 2080s and 1.9 m by 2100, relative to 2000–2004 (NPCC, 2015). By 2100, the 1% annual chance (100-year) floodplain could cover 1/3 of the city's total area under ARIM; around 1/5 of the area could be flooded during monthly high tides. Some low-lying locations could become permanently inundated by late century. Will New York City coastal resiliency initiatives, guided, in part by NPCC findings, suffice for very high sea levels? Additional research is needed to determine technological, environmental, or economic limitations to coastal protection and to decide when and where strategic relocation may become necessary.

Sea level rise↗

NASA NTRS2012

The Annual Glaciohydrology Cycle in the Ablation Zone of the Greenland Ice Sheet: Part 1. Hydrology Model

We apply a novel one-dimensional glacier hydrology model that calculates hydraulic head to the tidewater-terminating Sermeq Avannarleq flowline of the Greenland ice sheet. Within a plausible parameter space, the model achieves a quasi-steady-state annual cycle in which hydraulic head oscillates close to flotation throughout the ablation zone. Flotation is briefly achieved during the summer melt season along a approx.17 km stretch of the approx.50 km of flowline within the ablation zone. Beneath the majority of the flowline, subglacial conduit storage closes (i.e. obtains minimum radius) during the winter and opens (i.e. obtains maximum radius) during the summer. Along certain stretches of the flowline, the model predicts that subglacial conduit storage remains open throughout the year. A calculated mean glacier water residence time of approx.2.2 years implies that significant amounts of water are stored in the glacier throughout the year. We interpret this residence time as being indicative of the timescale over which the glacier hydrologic system is capable of adjusting to external surface meltwater forcings. Based on in situ ice velocity observations, we suggest that the summer speed-up event generally corresponds to conditions of increasing hydraulic head during inefficient subglacial drainage. Conversely, the slowdown during fall generally corresponds to conditions of decreasing hydraulic head during efficient subglacial drainage.

Colgan, William↗

NASA NTRS2014

Evaluation of the Surface Representation of the Greenland Ice Sheet in a General Circulation Model

Simulated surface conditions of the Goddard Earth Observing System model, version 5 (GEOS 5) atmospheric general circulation model (AGCM) are examined for the contemporary Greenland Ice Sheet (GrIS). A surface parameterization that explicitly models surface processes including snow compaction, meltwater percolation and refreezing, and surface albedo is found to remedy an erroneous deficit in the annual net surface energy flux and provide an adequate representation of surface mass balance (SMB) in an evaluation using simulations at two spatial resolutions. The simulated 1980-2008 GrIS SMB average is 24.7+/-4.5 cm yr(- 1) water-equivalent (w.e.) at.5 degree model grid spacing, and 18.2+/-3.3 cm yr(- 1) w.e. for 2 degree grid spacing. The spatial variability and seasonal cycle of the simulation compare favorably to recent studies using regional climate models, while results from 2 degree integrations reproduce the primary features of the SMB field. In comparison to historical glaciological observations, the coarser resolution model overestimates accumulation in the southern areas of the GrIS, while the overall SMB is underestimated. These changes relate to the sensitivity of accumulation and melt to the resolution of topography. The GEOS-5 SMB fields contrast with available corresponding atmospheric models simulations from the Coupled Model Intercomparison Project (CMIP5). It is found that only a few of the CMIP5 AGCMs examined provide significant summertime runoff, a dominant feature of the GrIS seasonal cycle. This is a condition that will need to be remedied if potential contributions to future eustatic change from polar ice sheets are to be examined with GCMs.

AGCM↗

NASA NTRS2023

An Evaluation of A Physics-Based Firn Model and A Semi-Empirical Firn Model Across the Greenland Ice Sheet (1980–2020)

The Greenland Ice Sheet's (GrIS) firn layer buffers the ice sheet's contribution to sea level rise by storing meltwater in its pore space. However, available pore space and meltwater retention capability is lost due to ablation of the firn layer and refreezing of meltwater as near-surface ice slabs in the firn. Understanding how firn properties respond to climate is important for constraining the GrIS's future contribution to sea level rise in a warming climate. Observations of firn density provide detailed information about firn properties, but they are spatially and temporally limited. Here we use two firn models, the physics-based SNOWPACK model and the Community Firn Model configured with a semi-empirical densification equation (CFM-GSFC), to quantify firn properties across the GrIS from 1980 through 2020. We use an identical forcing (Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) atmospheric reanalysis) for SNOWPACK and the CFM-GSFC in order to isolate firn model differences. To evaluate the models, we compare simulated firn properties, including firn air content (FAC), to measurements from the Surface Mass Balance and Snow on Sea Ice Working Group (SUMup) dataset of snow and firn density. Both models perform well (mean absolute percentage errors of 14 % in SNOWPACK and 16 % in the CFM-GSFC), though their performance is hindered by the spatial resolution of the atmospheric forcing. In the ice-sheet-wide simulations, the 1980–1995 average spatially integrated FAC (i.e., air volume in the firn) for the upper 100 m is 34 645 km 3 from SNOWPACK and 28 581 km 3 from the CFM-GSFC. The discrepancy in the magnitude of the modeled FAC stems from differences in densification with depth and variations in the sensitivity of the models to atmospheric forcing. In more recent years (2005–2020), both models simulate substantial depletion of pore space. During this period, the spatially integrated FAC across the entire GrIS decreases by 3.2 % (−66.6 km 3 yr −1 ) in SNOWPACK and 1.5 % (−17.4 km 3 yr −1 ) in the CFM-GSFC. These differing magnitudes demonstrate how model differences propagate throughout the FAC record. Over the full modeled record (1980–2020), SNOWPACK simulates a loss of pore space equivalent to 3 mm of sea level rise buffering, while the CFM-GSFC simulates a loss of 1 mm. The greatest depletion in FAC is along the margins and especially along the western margin where observations and models show the formation of near-surface, low-permeability ice slabs that may inhibit meltwater storage.

Megan Thompson-Munson↗

NASA NTRS2018

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↗

DOE OSTI2024

Direct Geologic Constraints on the Timing of Late Holocene Ice Thickening in the Amundsen Sea Embayment, Antarctica

Abstract Constraining past West Antarctic Ice Sheet (WAIS) change helps validate numerical models simulating future ice sheet dynamics. Following rapid deglaciation during the mid‐Holocene, ice near Thwaites Glacier was ∼35 m thinner than present; however, the timing of ice regrowth to its present configuration remains unknown. To fill this knowledge gap, we present cosmogenic nuclide exposure ages of cobbles from the surface of a moraine situated between Thwaites and Pope glaciers. We infer that the moraine formed and stabilized in the Late Holocene (∼1.4 ka) when a small glacier thickened. We also present a novel reconstruction of WAIS volume constrained by sea‐level data, which demonstrates that moraine formation coincided with a large‐scale WAIS readvance. Our new geologic constraints will help inform models of the solid Earth response to surface mass loading, improving our understanding of ice sheet dynamics in a vulnerable part of WAIS.

Nichols, Keir A. [Imperial College London London U↗

DOE OSTI2025

Improving Projections of Antarctic Ice Sheet Contribution to Sea‐Level Change Through 2300 by Capturing Gravitational, Rotational, and Deformational Effects

Ice-mass change induces regionally varying patterns of sea-level change due to gravitational, rotational, and deformational (GRD) effects, which in turn influence marine-based ice stability in Antarctica. For improved projection of the Antarctic Ice Sheet (AIS), there is a need for including GRD effects in modeling and improving understanding of basin-by-basin sensitivity of ice evolution to GRD effects under a range of climate scenarios. We couple a high-resolution, higher-order ice-sheet model with a 1D global sea-level model that fully captures GRD effects, and simulate ice evolution in Antarctica under the Ice Sheet Model Intercomparison Project for CMIP6 experiments. We perform two sets of coupled simulations incorporating 1D Maxwell solid Earth structure suitable for West and East Antarctica and show that the Amundsen Sea Embayment (ASE) in West Antarctica has the highest sensitivity to GRD effects—in high-emission scenarios, grounding-line retreat accelerates by hundreds of kilometers by 2300 without GRD effects, but GRD effects delay this retreat on a timescale of decades. However, we find that delay times do not show a clear relationship to the strength of climate forcing alone. Furthermore, GRD effects can influence ice-sheet dynamics more than the choice of climate model for a given emissions scenario. In contrast, East Antarctica exhibits minimal sensitivity to GRD effects throughout the study period. These findings underscore the critical role of GRD effects in shaping future West AIS evolution, highlighting the importance of constraining the regional 3D Earth structure and bed topography in West Antarctica, particularly the ASE.

58 GEOSCIENCES↗

NASA NTRS2014

Future Antarctic Bed Topography and Its Implications for Ice Sheet Dynamics

The Antarctic bedrock is evolving as the solid Earth responds to the past and ongoing evolution of the ice sheet. A recently improved ice loading history suggests that the Antarctic Ice Sheet (AIS) has generally been losing its mass since the Last Glacial Maximum. In a sustained warming climate, the AIS is predicted to retreat at a greater pace, primarily via melting beneath the ice shelves.We employ the glacial isostatic adjustment (GIA) capability of the Ice Sheet System Model (ISSM) to combine these past and future ice loadings and provide the new solid Earth computations for the AIS.We find that past loading is relatively less important than future loading for the evolution of the future bed topography. Our computations predict that the West Antarctic Ice Sheet (WAIS) may uplift by a few meters and a few tens of meters at years AD 2100 and 2500, respectively, and that the East Antarctic Ice Sheet is likely to remain unchanged or subside minimally except around the Amery Ice Shelf. The Amundsen Sea Sector in particular is predicted to rise at the greatest rate; one hundred years of ice evolution in this region, for example, predicts that the coastline of Pine Island Bay will approach roughly 45mmyr−1 in viscoelastic vertical motion. Of particular importance, we systematically demonstrate that the effect of a pervasive and large GIA uplift in the WAIS is generally associated with the flattening of reverse bed slope, reduction of local sea depth, and thus the extension of grounding line (GL) towards the continental shelf. Using the 3-D higher-order ice flow capability of ISSM, such a migration of GL is shown to inhibit the ice flow. This negative feedback between the ice sheet and the solid Earth may promote stability in marine portions of the ice sheet in the future.

ice dynamics↗

DOE OSTI2025

Emulation With Uncertainty Quantification of Regional Sea‐Level Change Caused by the Antarctic Ice Sheet

Abstract Projecting regional sea‐level change under various climate‐change scenarios typically involves running forward simulations of the Earth's gravitational, rotational and deformational (GRD) response to ice‐mass change, which requires substantial computational cost if applied to probabilistic frameworks requiring thousands to millions of samples. Here we build emulators of regional sea‐level change at 27 coastal locations, due to the GRD effects associated with future Antarctic Ice Sheet mass change over the 21st century. The emulators are evaluated against a numerical sea‐level model applied to an ensemble of ice‐sheet model simulations of the Antarctic Ice Sheet through 2100. We build a physics‐based emulator using a recent sensitivity kernel approach and compare it to machine learning based emulators (neural network and conditional variational autoencoder methods). In order to quantify uncertainty, we derive well‐calibrated prediction intervals for regional sea‐level change via split‐conformal inference and linear regression, and show that Monte Carlo dropout does not yield well‐calibrated uncertainties in this instance. We also demonstrate substantial gains in computational efficiency using both the physics‐based emulator and neural networks in comparison to the numerical model for the complete regional sea‐level solution. Overall, we find the physics‐based emulator modestly outperforms the machine learning emulators for this problem.

58 GEOSCIENCES↗

NASA NTRS2017

Determining Ice-Sheet Uplift Surrounding Subglacial Lakes with a Viscous Plate Model

We develop a viscous model of plate bending suitable for studying ice-sheet flexure due to subglacial lake filling and draining, and apply this model to determine the area of ice-sheet uplift surrounding a subglacial lake. The choice of a viscous model reflects our interest in Antarctic subglacial lakes, which can fill and drain on time scales of months to decades. Experiments with idealized lake shapes show that the size of the uplift area relative to lake area depends on subglacial water pressure and ice-sheet thickness, with the viscous material parameters scaling the magnitude of uplift rate within this area. The water pressure therefore has a strong control on the evolution of the lake shape and related subglacial hydrological development, but is not yet well constrained by observations. Due to the likelihood that ice flexure will affect subglacial lake filling and draining, we suggest that the insights of this study should be applied to development of a realistic ice sheet-hydrological coupled model.

ice sheets↗

NASA NTRS1997

A Comparison of Geoid Undulations for West Central Greenland

The accuracy of a new local gravity field model, GEOID94A, is examined at a site on the western Greenland ice sheet. The model, developed by the Danish National Survey and Cadastre, incorporates several new gravity data sets including an extensive amount of airborne gravity data. Model-derived geoid undulations were compared to independently determined undulations found by differencing the elevations from Global Positioning System controlled airborne laser altimetry and optical leveling surveys. Differences between the two sets of undulations were less than +/- 6 cm RMS. The comparison improved (+/- 5 cm RMS) when GEOID94A undulations were adjusted by local gravity observations also acquired at the site. Our comparisons demonstrate that GEOID94A adequately models the long to intermediate wavelengths of the gravity field. We conclude that GEOID94A constitutes a reliable reference model for studies of Greenland's gravity field.

Roman, Daniel R.↗

NASA NTRS1999

Interdisciplinary Determination of Ice Sheet Accumulation Patterns: Combined Atmospheric Modeling and Field and Remote Sensing Studies

This research is focusing on two related areas that are fundamental to the NASA PARCA (Program for Arctic Regional Climate Assessment) program. The primary research area is the determination of the amount, rate, and timing of accumulation at distributed sites in the dry snow zone of Greenland and evaluation of these results in light of accumulation modeling results. The secondary research area is the calibration of the isotope "thermometer" at these ice sheet sites as well as the determination of long-term temperature trends in Greenland.

Shuman, C. A.↗

NASA NTRS1984

Origin of the Greenland Ice Sheet

A climate model resolving the seasonal cycle and the two horizontal dimensions was developed and applied to several problems of current interest. Models of this type are useful when for various reasons a general circulation model experiment is not warranted or not feasible. For example, in cases where the signal to natural variability is small it may be advantageous to first consider such a statistical dynamical model because extremely long runs may be necessary in the application. In this case the simpler statistical dynamical model serves as a pilot study device. The model developed is a thermodynamic model whose solution yields the equilibrium seasonal cycle for the surface temperature field over the globe. The model is essentially a statement of the conservation of heat energy for individual columns of the earth atmosphere system. Various terms such as the infrared radiation flux to space are parameterized with earth radiation budget data from satellites such as Nimbus 6. The primary agent modulating the seasonal cycle amplitude is the heat capacity per unit area which is a strong function of surface type -- ocean surface can store 60 times more heat per unit time than land. By adjusting its few empirical parameters the model can be brought into remarkable agreement with the observed seasonal cycle. The model is then very useful for looking at the dependence of the seasonal cycle of the temperature on such externally defined variables as the Earth's orbital elements (eccentricity, tilt, precession of equinoxes) or the configuration of land-sea geography which can be changed by continental drift.

North, G. R.↗

NASA NTRS1998

Nonlinear Insolation Forcing: A Physical Mechanism for Climate Change

This paper focuses on recent advances in the understanding of nonlinear insolation forcing for climate change. The amplitude-frequency resonances in the insolation variations induced by the Earth's changing obliquity are emergent and may provide a physical mechanism to drive the glaciation cycles. To establish the criterion that nonlinear insolation forcing is responsible for major climate changes, the cooperative phenomena between the frequency and amplitude of the insolation are defined as insolation pulsation. Coupling of the insolation frequency and amplitude variations has established an especially new and interesting series of insolation pulses. These pulses would modulate the insolation in such a way that the mode of insolation variations could be locked to generate the 100-kyr ice age cycle which is a long-time geophysical puzzle. The nonlinear behavior of insolation forcing is tested by energy balance and ice sheet climate models and the physical mechanism behind this forcing is explained in terms of pulse duration in the incoming solar radiation. Calculations of the solar energy flux at the top of the atmosphere show that the duration of the negative and positive insolation pulses is about 2 thousand years which is long enough to prolong glaciation into deep ice ages and cause rapid melting of large ice sheets in the high latitudes of the northern hemisphere. We have performed numerical simulations of climate response to nonlinear insolation forcing for the past 2 million years. Our calculated results of temperature fluctuations are in good agreement with the climate cycles as seen in the terrestrial biogenic silica (BDP-96-2) data as well as in the marine oxygen isotope (delta(sup 18)O) records.

Liu, H. S.↗

NASA NTRS2001

Seasonal and Interannual Variations of Ice Sheet Surface Elevation at the Summit of Greenland: Observed and Modeled

Observed seasonal and interannual variations in the surface elevation over the summit of the Greenland ice sheet are modeled using a new temperature-dependent formulation of firn-densification and observed accumulation variations. The observed elevation variations are derived from ERS (European Remote Sensing)-1 and ERS-2 radar altimeter data for the period between April 1992 and April 1999. A multivariate linear/sine function is fitted to an elevation time series constructed from elevation differences measured by radar altimetry at orbital crossovers. The amplitude of the seasonal elevation cycle is 0.25 m peak-to-peak, with a maximum in winter and a minimum in summer. Inter-annually, the elevation decreases to a minimum in 1995, followed by an increase to 1999, with an overall average increase of 4.2 cm a(exp -1) for 1992 to 1999. Our densification formulation uses an initial field-density profile, the AWS (automatic weather station) surface temperature record, and a temperature-dependent constitutive relation for the densification that is based on laboratory measurements of crystal growth rates. The rate constant and the activation energy commonly used in the Arrhenius-type constitutive relation for firn densification are also temperature dependent, giving a stronger temperature and seasonal amplitudes about 10 times greater than previous densification formulations. Summer temperatures are most important, because of the strong non-linear dependence on temperature. Much of firn densification and consequent surface lowering occurs within about three months of the summer season, followed by a surface build-up from snow accumulation until spring. Modeled interannual changes of the surface elevation, using the AWS measurements of surface temperature and accumulation and results of atmospheric modeling of precipitation variations, are in good agreement with the altimeter observations. In the model, the surface elevation decreases about 20 cm over the seven years due to more compaction driven by increasing summer temperatures. The minimum elevation in 1995 is driven mainly by a temporary accumulation decrease and secondarily by warmer temperatures. However, the overall elevation increase over the seven years is dominated by the accumulation increase in the later years.

Zwally, H. Jay↗

NASA NTRSDate not supplied

The Impact of Ice Sheet Surface Representation on Surface Mass Balance in the Goddard Earth Observing System

Surface conditions and their impacts on surface mass balance (SMB) play important roles in ice sheet dynamics and ice sheet interactions with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) – an ecosystem of models and reanalyses – represent ice sheet and glacier SMB components, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. We find the successful representation of ice sheet surface mass balance in MERRA-2, and similar systems, is partly due to compensating biases in surface energy budget. In the atmospheric system, subsequent changes to cloud microphysics and in the longwave radiative transfer model are found to have reduced biases in surface net longwave radiation fluxes while exacerbating surface net shortwave radiation biases, producing erroneously high near-surface temperatures and surface melt in summer months. These issues can be exacerbated by poor experiment initialization and the use of two-moment cloud microphysics within the GEOS ocean-atmosphere coupled system. Here, we document the spatial and temporal extent and causes of these biases in the ice sheet surface energy budget across GEOS systems and implement a range of model improvements to mitigate these issues. We examine the combined impact on these changes on regional energy budget and SMB using both free running and replay experiments (to simulate the impact in reanalyses). Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet SMB in GMAO’s future reanalyses and forecasting systems.

Lauren C Andrews↗

NASA NTRSDate not supplied

Ice Sheet and Ice Shelf Freshwater Forcing of the NASA GISS ModelE GCM

Recent mass loss from ice sheets and ice shelves is now persistent and prolonged enough that it impacts downstream oceanographic conditions. To demonstrate this, we use an ensemble of coupled GISS-E2.1-G (a CMIP6 model) simulations forced with historical estimates of anomalous freshwater, in addition to other climate forcings, from 1990 through 2019. In this ensemble there are detectable differences in zonal-mean sea surface temperatures (SST) and sea ice in the Southern Ocean, and in regional sea level around Antarctica and in the western North Atlantic. These impacts mostly improve the model's representation of historical changes, including reversing the forced trends in Southern Ocean surface temperature and Antarctic sea ice. The changes in SST may have implications for estimates of the SST pattern effect on climate sensitivity and for cloud feedbacks. We conclude that the changes are sufficiently large that model groups should strive to include more accurate estimates of these drivers in all-forcing historical simulations in future coupled model intercomparisons.

Ice sheet↗