Applying the energy- and water balance model for incorporation of the cryospheric component into a climate model. Part II. Modeled mass balance on the green land ice sheet surface

The Greenland ice sheet is a very important potential source of fresh water inflow to the World Ocean under warming climate conditions. Apparently, it was the same during the Last Interglacial 130-115 thousand years ago. In order to quantify input of the Greenland ice sheet to the rise of the global mean sea level in the past or in the future, we include a surface mass balance model block into the Earth System Model. The computational algorithm is based on the calculation of energy balance on the ice sheet surface. The key tuning parameter of the model is the daily amplitude of air surface temperature. It defines the area and the rate of snow or ice melting. The range of possible values of this parameter is determined during a series of numerical experiments. High sensitivity of meltwater runoff volume to surface air temperature amplitude is revealed.     

Snowpack model estimates of the mass balance of the Greenland ice sheet and its changes over the twentyfirst century

The performance of a snow cover model in capturing the ablation on the Greenland ice sheet is evaluated. This model allows an explicit calculation of the formation of melt water, of the fraction of melt water which re-freezes, and of runoff in the ablation region. The input climate variables to the snowpack model come from two climate models. While the higher resolution general circulation model (ECHAM 4), is closest to observations in its estimate of accumulation, it fails to give accurate results in its predictions of runoff, primarily in the southern half of the ice sheet. The two-dimensional low-resolution climate model (MIT 2D LO) produces estimates of runoff from the Greenland ice sheet within the range of uncertainty of the Inter governmental Panel on Climate Change (IPCC1) 1995 estimates. Both models reproduce some of the characteristics of the extent of the wet snow zone observed with satellite remote sensing; the MIT model is closer to observations in terms of areal extent and intensity of the melting in the southern half of the ice-sheet in July and August while the ECHAM model reproduces melting in the northern half of the ice sheet well. Changes in runoff from Greenland and Antarctica are often cited as one of the major concerns linked to anthropogenic changes in climate. Because it is based on physical principles and relies on the surface energy balance as input, the snow cover model can respond to the current climatic forcing as well as to future changes in climate on the century time scale without the limitations inherent in empirical parametrizations. For a reference climate scenario similar to the IPCC's IS92a, the model projects that the Greenland ice sheet does not contribute significantly to changes in the level of the ocean over the twenty-first century. Increases in accumulation over the central portion of the ice sheet offset most of the increase in melting and runoff, which takes place along the margins of the ice sheet. The range of uncertainty in the predictions of sea-level rise is estimated by repeating the calculation with the MIT model for seven climate change scenarios. The range is –0.5 to 1.7 cm.     

Effects of a melted greenland ice sheet on climate, vegetation, and the cryosphere

This paper investigates the possible implications for the earth-system of a melting of the Greenland ice-sheet. Such a melting is a possible result of increased high latitude temperatures due to increasing anthropogenic greenhouse gas emissions. Using an atmosphere-ocean general circulation model (AOGCM), we investigate the effects of the removal of the ice sheet on atmospheric temperatures, circulation, and precipitation. We find that locally over Greenland, there is a warming associated directly with the altitude change in winter, and the altitude and albedo change in summer. Outside of Greenland, the largest signal is a cooling over the Barents sea in winter. We attribute this cooling to a decrease in poleward heat transport in the region due to changes to the time mean circulation and eddies, and interaction with sea-ice. The simulated climate is used to force a vegetation model and an ice-sheet model. We find that the Greenland climate in the absence of an ice sheet supports the growth of trees in southern Greenland, and grass in central Greenland. We find that the ice sheet is likely to regrow following a melting of the Greenland ice sheet, the subsequent rebound of its bedrock, and a return to present day atmospheric CO₂ concentrations. This regrowth is due to the high altitude bedrock in eastern Greenland which allows the growth of glaciers which develop into an ice sheet.     

Impact of Greenland and Antarctic ice sheet interactions on climate sensitivity

We use the Earth system model of intermediate complexity LOVECLIM to show the effect of coupling interactive ice sheets on the climate sensitivity of the model on a millennial time scale. We compare the response to a 2×CO₂ warming scenario between fully coupled model versions including interactive Greenland and Antarctic ice sheet models and model versions with fixed ice sheets. For this purpose an ensemble of different parameter sets have been defined for LOVECLIM, covering a wide range of the model??s sensitivity to greenhouse warming, while still simulating the present-day climate and the climate evolution over the last millennium within observational uncertainties. Additional freshwater fluxes from the melting ice sheets have a mitigating effect on the model??s temperature response, leading to generally lower climate sensitivities of the fully coupled model versions. The mitigation is effectuated by changes in heat exchange within the ocean and at the sea?Cair interface, driven by freshening of the surface ocean and amplified by sea?Cice-related feedbacks. The strength of the effect depends on the response of the ice sheets to the warming and on the model??s climate sensitivity itself. The effect is relatively strong in model versions with higher climate sensitivity due to the relatively large polar amplification of LOVECLIM. With the ensemble approach in this study we cover a wide range of possible model responses.     

Long-term ice sheet–climate interactions under anthropogenic greenhouse forcing simulated with a complex Earth System Model

Several multi-century and multi-millennia simulations have been performed with a complex Earth System Model (ESM) for different anthropogenic climate change scenarios in order to study the long-term evolution of sea level and the impact of ice sheet changes on the climate system. The core of the ESM is a coupled coarse-resolution Atmosphere–Ocean General Circulation Model (AOGCM). Ocean biogeochemistry, land vegetation and ice sheets are included as components of the ESM. The Greenland Ice Sheet (GrIS) decays in all simulations, while the Antarctic ice sheet contributes negatively to sea level rise, due to enhanced storage of water caused by larger snowfall rates. Freshwater flux increases from Greenland are one order of magnitude smaller than total freshwater flux increases into the North Atlantic basin (the sum of the contribution from changes in precipitation, evaporation, run-off and Greenland meltwater) and do not play an important role in changes in the strength of the North Atlantic Meridional Overturning Circulation (NAMOC). The regional climate change associated with weakening/collapse of the NAMOC drastically reduces the decay rate of the GrIS. The dynamical changes due to GrIS topography modification driven by mass balance changes act first as a negative feedback for the decay of the ice sheet, but accelerate the decay at a later stage. The increase of surface temperature due to reduced topographic heights causes a strong acceleration of the decay of the ice sheet in the long term. Other feedbacks between ice sheet and atmosphere are not important for the mass balance of the GrIS until it is reduced to 3/4 of the original size. From then, the reduction in the albedo of Greenland strongly accelerates the decay of the ice sheet.     

Sensitivity of the last glacial inception to initial and surface conditions

We investigate the sensitivity of simulations of the last glacial inception (LGI) with respect to initial (size of the Greenland ice sheet) and surface (state of ocean/vegetation) conditions and two different CO₂ reconstructions. Utilizing the CLIMBER-2 Earth system model, we obtain the following results: (a) ice-sheet expansion in North America at the end of the Eemian can be reduced or even completely suppressed when pre-industrial or Eemian ocean/vegetation is prescribed. (b) A warmer surrounding ocean and, in particular, a large Laurentide ice sheet reduce the size of the Greenland ice sheet before and during the LGI. (c) A changing ocean contributes much stronger to the expansion of the Laurentide ice sheet when we apply the CO₂ reconstruction according to Barnola et al. (Nature 329:408–414, 1987) instead of Petit et al. (Nature 399:429–436, 1999). (d) In the fully coupled model, the CO₂ reconstruction used has only a small impact on the simulated ice sheets but it does impact the course of the climatic variables. (e) For the Greenland ice sheet, two equilibrium states exist under the insolation and CO₂ forcing at 128,000 years before present (128 kyear BP); the one with an ice sheet reduced by about one quarter as compared to its simulated pre-industrial size and the other with nearly no inland ice in Greenland. (f) Even the extreme assumption of no ice sheet in Greenland at the beginning of our transient simulations does not alter the simulated expansion of northern hemispheric ice sheets at the LGI.     

Thresholds for irreversible decline of the Greenland ice sheet

The Greenland ice sheet will decline in volume in a warmer climate. If a sufficiently warm climate is maintained for a few thousand years, the ice sheet will be completely melted. This raises the question of whether the decline would be reversible: would the ice sheet regrow if the climate cooled down? To address this question, we conduct a number of experiments using a climate model and a high-resolution ice-sheet model. The experiments are initialised with ice sheet states obtained from various points during its decline as simulated in a high-CO₂ scenario, and they are then forced with a climate simulated for pre-industrial greenhouse gas concentrations, to determine the possible trajectories of subsequent ice sheet evolution. These trajectories are not the reverse of the trajectory during decline. They converge on three different steady states. The original ice-sheet volume can be regained only if the volume has not fallen below a threshold of irreversibility, which lies between 80 and 90% of the original value. Depending on the degree of warming and the sensitivity of the climate and the ice-sheet, this point of no return could be reached within a few hundred years, sooner than CO₂ and global climate could revert to a pre-industrial state, and in that case global sea level rise of at least 1.3 m would be irreversible. An even larger irreversible change to sea level rise of 5 m may occur if ice sheet volume drops below half of its current size. The set of steady states depends on the CO₂ concentration. Since we expect the results to be quantitatively affected by resolution and other aspects of model formulation, we would encourage similar investigations with other models.     

Multistability of the Greenland ice sheet and the effects of an adaptive mass balance formulation

The effect of a warmer climate on the Greenland ice sheet as well as its ability to regrow from a reduced geometry is important knowledge when studying future climate. Here we use output from a general circulation model to construct adaptive temperature and precipitation patterns to force an ice flow model off-line taking into consideration that the patterns change in a non-uniform way (both spatially and temporally) as the geometry of the ice sheet evolves and as climate changes. In a series of experiments we investigate the retreat from the present day configuration, build-up from ice free conditions of the ice sheet during a warmer-than-present climate and how the ice sheet moves between states. The adaptive temperature and accumulation patterns as well as two different constant-pattern formulations are applied and all experiments are run to steady state. All results fall into four different groups of geometry regardless of the applied accumulation pattern and initial state. We find that the ice sheet is able to survive and build up at higher temperatures using the more realistic adaptive patterns compared to the classic constant patterns. In contrast, decay occurs at considerably higher temperatures than build-up when the other formulations are used. When studying the motion between states it is clear that the initial state is crucial for the result. The ice sheet is thus multistable at least for certain temperature forcings, and this implies that the ice sheet not does not necessarily return to its initial configuration after a temperature excursion.     

Greenland Ice Sheet Contribution to Future Global Sea Level Rise based on CMIP5 Models简

Sea level rise （SLR） is one of the major socioeconomic risks associated with global warming. Mass losses from the Greenland ice sheet （GrIS） will be partially responsible for future SLR, although there are large uncertainties in modeled climate and ice sheet behavior. We used the ice sheet model SICOPOLIS （Simulation COde for POLythermal Ice Sheets） driven by climate projections from 20 models in the fifth phase of the Coupled Model Intercomparison Project （CMIP5） to estimate the GrlS contribution to global SLR. Based on the outputs of the 20 models, it is estimated that the GrIS will contribute 0-16 （0-27） cm to global SLR by 2100 under the Representative Concentration Pathways （RCP） 4.5 （RCP 8.5） scenarios. The projected SLR increases further to 7-22 （7-33） cm with 2~basal sliding included. In response to the results of the multimodel ensemble mean, the ice sheet model projects a global SLR of 3 cm and 7 cm （10 cm and 13 cm with 2~basal sliding） under the RCP 4.5 and RCP 8.5 scenarios, respectively. In addition, our results suggest that the uncertainty in future sea level projection caused by the large spread in climate projections could be reduced with model-evaluation and the selective use of model outputs.     

Quantifying Arctic contributions to climate predictability in a regional coupled ocean-ice-atmosphere model

The relative importance of regional processes inside the Arctic climate system and the large scale atmospheric circulation for Arctic interannual climate variability has been estimated with the help of a regional Arctic coupled ocean-ice-atmosphere model. The study focuses on sea ice and surface climate during the 1980s and 1990s. Simulations agree reasonably well with observations. Correlations between the winter North Atlantic Oscillation index and the summer Arctic sea ice thickness and summer sea ice extent are found. Spread of sea ice extent within an ensemble of model runs can be associated with a surface pressure gradient between the Nordic Seas and the Kara Sea. Trends in the sea ice thickness field are widely significant and can formally be attributed to large scale forcing outside the Arctic model domain. Concerning predictability, results indicate that the variability generated by the external forcing is more important in most regions than the internally generated variability. However, both are in the same order of magnitude. Local areas such as the Northern Greenland coast together with Fram Straits and parts of the Greenland Sea show a strong importance of internally generated variability, which is associated with wind direction variability due to interaction with atmospheric dynamics on the Greenland ice sheet. High predictability of sea ice extent is supported by north-easterly winds from the Arctic Ocean to Scandinavia.     

Climate modification by future ice sheet changes and consequences for ice sheet mass balance

The future evolution of global ice sheets under anthropogenic greenhouse forcing and its impact on the climate system, including the regional climate of the ice sheets, are investigated with a comprehensive earth system model consisting of a coupled Atmosphere–Ocean General Circulation Model, a dynamic vegetation model and an ice sheet model. The simulated control climate is realistic enough to permit a direct coupling of the atmosphere and ice sheet components, avoiding the use of anomaly coupling, which represents a strong improvement with respect to previous modelling studies. Glacier ablation is calculated with an energy-balance scheme, a more physical approach than the commonly used degree-day method. Modifications of glacier mask, topographic height and freshwater fluxes by the ice sheets influence the atmosphere and ocean via dynamical and thermodynamical processes. Several simulations under idealized scenarios of greenhouse forcing have been performed, where the atmospheric carbon dioxide stabilizes at two and four times pre-industrial levels. The evolution of the climate system and the ice sheets in the simulations with interactive ice sheets is compared with the simulations with passively coupled ice sheets. For a four-times CO₂ scenario forcing, a faster decay rate of the Greenland ice sheet is found in the non-interactive case, where melting rates are higher. This is caused by overestimation of the increase in near-surface temperature that follows the reduction in topographic height. In areas close to retreating margins, melting rates are stronger in the interactive case, due to changes in local albedo. Our results call for careful consideration of the feedbacks operating between ice sheets and climate after substantial decay of the ice sheets.     

Global Warming and the Greenland Ice Sheet

The Greenland coastal temperatures have followed the early 20th century global warming trend. Since 1940, however, the Greenland coastal stations data have undergone predominantly a cooling trend. At the summit of the Greenland ice sheet the summer average temperature has decreased at the rate of 2.2 °C per decade since the beginning of the measurements in 1987. This suggests that the Greenland ice sheet and coastal regions are not following the current global warming trend. A considerable and rapid warming over all of coastal Greenland occurred in the 1920s when the average annual surface air temperature rose between 2 and 4 °C in less than ten years (at some stations the increase in winter temperature was as high as 6 °C). This rapid warming, at a time when the change in anthropogenic production of greenhouse gases was well below the current level, suggests a high natural variability in the regional climate. High anticorrelations (r = ?0.84 to?0.93) between the NAO (North Atlantic Oscillation) index and Greenland temperature time series suggest a physical connection between these processes. Therefore, the future changes in the NAO and Northern Annular Mode may be of critical consequence to the future temperature forcing of the Greenland ice sheet melt rates.     

A continuous simulation of global ice volume over the past 1 million years with 3-D ice-sheet models

Sea-level records show large glacial-interglacial changes over the past million years, which on these time scales are related to changes of ice volume on land. During the Pleistocene, sea-level changes induced by ice volume are largely caused by the waxing and waning of the large ice sheets in the Northern Hemisphere. However, the individual contributions of ice in the Northern and Southern Hemisphere are poorly constrained. In this study, for the first time a fully coupled system of four 3-D ice-sheet models is used, simulating glaciations on Eurasia, North America, Greenland and Antarctica. The ice-sheet models use a combination of the shallow ice and shelf approximations to determine sheet, shelf and sliding velocities. The framework consists of an inverse forward modelling approach to derive a self-consistent record of temperature and ice volume from deep-sea benthic δ¹⁸O data over the past 1 million years, a proxy for ice volume and temperature. It is shown that for both eustatic sea level and sea water δ¹⁸O changes, the Eurasian and North American ice sheets are responsible for the largest part of the variability. The combined contribution of the Antarctic and Greenland ice sheets is about 10 % for sea level and about 20 % for sea water δ¹⁸O during glacial maxima. However, changes in interglacials are mainly caused by melt of the Greenland and Antarctic ice sheets, with an average time lag of 4 kyr between melt and temperature. Furthermore, we have tested the separate response to changes in temperature and sea level for each ice sheet, indicating that ice volume can be significantly influenced by changes in eustatic sea level alone. Hence, showing the importance of a simultaneous simulation of all four ice sheets. This paper describes the first complete simulation of global ice-volume variations over the late Pleistocene with the possibility to model changes above and below present-day ice volume, constrained by observations of benthic δ¹⁸O proxy data.     

The 1988–2003 Greenland ice sheet melt extent using passive microwave satellite data and a regional climate model

Measurements from ETH-Camp and JAR1 AWS (West Greenland) as well as coupled atmosphere-snow regional climate simulations have highlighted flaws in the cross-polarized gradient ratio (XPGR) technique used to identify melt from passive microwave satellite data. It was found that dense clouds (causing notably rainfall) on the ice sheet severely perturb the XPGR melt signal. Therefore, the original XPGR melt detection algorithm has been adapted to better incorporate atmospheric variability over the ice sheet and an updated melt trend for the 1988–2003 period has been calculated. Compared to the original algorithm, the melt zone area increase is eight times higher (from 0.2 to 1.7% year⁻¹). The increase is higher with the improved XPGR technique because rainfall also increased during this period. It is correlated to higher atmospheric temperatures. Finally, the model shows that the total ice sheet runoff is directly proportional to the melt extent surface detected by satellites. These results are important for the understanding of the effect of Greenland melting on the stability of the thermohaline circulation.     

Climate sensitivity studies of the Greenland ice sheet using satellite AVHRR,SMMR, SSM/I and in situ data

Summary The feasibility of using satellite data for climate research over the Greenland ice sheet is discussed. In particular, we demonstrate the usefulness of Advanced Very High Resolution Radiometer (AVHRR) Local Area Coverage (LAC) and Global Area Coverage (GAC) data for narrow-band albedo retrieval. Our study supports the use of lower resolution AVHRR (GAC) data for process studies over most of the Greenland ice sheet. Based on LAC data time series analysis, we can resolve relative albedo changes on the order of 2–5%. In addition, we examine Scanning Multichannel Microwave Radiometer (SMMR) and Special Sensor Microwave Imager (SSM/I) passive microwave data for snow typing and other signals of climatological significance. Based on relationships between in situ measurements and horizontally polarized 19 and 37 GHz observations, wet snow regions are identified. The wet snow regions increase in aerial percentage from 9% of the total ice surface in June to a maximum of 26% in August 1990. Furthermore, the relationship between brightness temperatures and accumulation rates in the northeastern part of Greenland is described. We found a consistent increase in accumulation rate for the northeastern part of the ice sheet from 1981 to 1986.With 16 Figures     

GCM simulations of the Last Glacial Maximum surface climate of Greenland and Antarctica

 The LMDz variable grid GCM was used to simulate the Last Glacial Maximum (LGM, 21 ky Bp.) climate of Greenland and Antarctica at a spatial resolution of about 100 km.The high spatial resolution allows to investigate the spatial variability of surface climate change signals, and thus to address the question whether the sparse ice core data can be viewed as representative for the regional scale climate change. This study addresses primarily surface climate parameters because these can be checked against the, limited, ice core record. The changes are generally stronger for Greenland than for Antarctica, as the imposed changes of the forcing boundary conditions (e.g., sea surface temperatures) are more important in the vicinity of Greenland. Over Greenland, and to a limited extent also in Antarctica, the climate shows stronger changes in winter than in summer. The model suggests that the linear relationship between the surface temperature and inversion strength is modified during the LGM. The temperature dependency of the moisture holding capacity of the atmosphere alone cannot explain the strong reduction in snowfall over central Greenland; atmospheric circulation changes also play a crucial role. Changes in the high frequency variability of snowfall, atmospheric pressure and temperature are investigated and possible consequences for the interpretation of ice core records are discussed. Using an objective cyclone tracking scheme, the importance of changes of the atmospheric dynamics off the coasts of the ice sheets, especially for the high frequency variability of surface climate parameters, is illustrated. The importance of the choice of the LGM ice sheet topography is illustrated for Greenland, where two different topographies have been used, yielding results that differ quite strongly in certain nontrivial respects. This means that the paleo-topography is a significant source of uncertainty for the modelled paleoclimate. The sensitivity of the Greenland LGM climate to the prescribed sea surface conditions is examined by using two different LGM North Atlantic data sets. Received: 23 October 1997 / Accepted: 17 March 1998     

Linear response functions to project contributions to future sea level

We propose linear response functions to separately estimate the sea-level contributions of thermal expansion and solid ice discharge from Greenland and Antarctica. The response function formalism introduces a time-dependence which allows for future rates of sea-level rise to be influenced by past climate variations. We find that this time-dependence is of the same functional type, R(t) ～ t ^α, for each of the three subsystems considered here. The validity of the approach is assessed by comparing the sea-level estimates obtained via the response functions to projections from comprehensive models. The pure vertical diffusion case in one dimension, corresponding to α =  ?0.5, is a valid approximation for thermal expansion within the ocean up to the middle of the twenty first century for all Representative Concentration Pathways. The approximation is significantly improved for α =  ? 0.7. For the solid ice discharge from Greenland we find an optimal value of α =  ?0.7. Different from earlier studies we conclude that solid ice discharge from Greenland due to dynamic thinning is bounded by 0.42 m sea-level equivalent. Ice discharge induced by surface warming on Antarctica is best captured by a positive value of α = 0.1 which reflects the fact that ice loss increases with the cumulative amount of heat available for softening the ice in our model.     

Heat,momentum and moisture budgets of the katabatic layer over the melting zone of the west Greenland ice sheet in summer

The budgets of momentum, heat and moisture of the atmospheric boundary layer overlying the melting zone of the west Greenland ice sheet during an 8-day period in summer are calculated. To do so, the governing budget equations are derived and presented in terms of vertically averaged quantities. Moreover, stationarity is assumed in the present study. Measurements collected during the GIMEX-91 experiment are used to calculate the contribution of the different terms in the equations to the budget.During summer, a well developed katabatic wind system is present over the melting zone of the Greenland ice sheet. The budgets show that advection in the katabatic layer is small for momentum, heat and humidity, when the horizontal length scale of the integration area is sufficiently large (>50 km). This indicates that in principle one-dimensional atmospheric models can be used to study the boundary layer over the melting zone of the Greenland ice sheet. The background stratification plays a crucial role in the heat and moisture budget. Vertical divergence of longwave radiation provides one-third and the turbulent flux of sensible heat the rest of the cooling of the boundary layer. Moisture is added to the boundary layer by evaporation which is a significant term in the moisture budget. Negative buoyancy (katabatic forcing) dominates the momentum budget in the downslope direction. Coriolis forcing is important, stressing the large spatial scale of the katabatic winds on the Greenland ice sheet.     

Reconstruction of Climate of the Eemian Interglacial Using an Earth System Model. Part 2. The Response of the Greenland Ice Sheet to Climate Change

In the framework of the study of the Eemian interglacial we consider the role of the Greenland ice sheet in the rise of the mean level of the World Ocean. Its contribution estimated as 2 m confirms the newest estimates based on the model results and on the proxy data analysis. In the beginning of the Eemian interglacial (earlier than 126 thousand years ago) mass lost occurs through the marine margin of the sheet. During the next five millennia, the negative surface mass balance plays the leading role. Taking into account the contribution of Greenland ice sheet, ocean thermal expansion, and the melting of mountain glaciers and ice caps, it is very probable that the West Antarctic ice sheet was the main source of the global sea level growth equal to 6–9 m the compared to the present.     

Late Pliocene to Pleistocene sensitivity of the Greenland Ice Sheet in response to external forcing and internal feedbacks

The timing and nature of ice sheet variations on Greenland over the last ～5 million years remain largely uncertain. Here, we use a coupled climate-vegetation-ice sheet model to determine the climatic sensitivity of Greenland to combined sets of external forcings and internal feedbacks operating on glacial-interglacial timescales. In particular, we assess the role of atmospheric pCO₂, orbital forcing, and vegetation dynamics in modifying thresholds for the onset of glaciation in late Pliocene and Pleistocene. The response of circum-Arctic vegetation to declining levels of pCO₂ (from 400 to 200 ppmv) and decreasing summer insolation includes a shift from boreal forest to tundra biomes, with implications for the surface energy balance. The expansion of tundra amplifies summer surface cooling and heat loss from the ground, leading to an expanded summer snow cover over Greenland. Atmospheric and land surface fields respond to forcing most prominently in late spring-summer and are more sensitive at lower Pleistocene-like levels of pCO₂. We find cold boreal summer orbits produce favorable conditions for ice sheet growth, however simulated ice sheet extents are highly dependent on both background pCO₂ levels and land-surface characteristics. As a result, late Pliocene ice sheet configurations on Greenland differ considerably from late Pleistocene, with smaller ice caps on high elevations of southern and eastern Greenland, even when orbital forcing is favorable for ice sheet growth.