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 CO2-induced thickening/thinning of the Greenland and Antarctic ice sheets as simulated by a GCM (CCMI) and an ice-sheet model

Scaling analysis shows that the mean thickness of an ice sheet depends on the product of two poorly known quantities, the ice viscosity and the net snow accumulation rate. We adjust the viscosity of an ice sheet in order to get a consistent value of this product for the present-day ice sheet volume and area given the net snow accumulation rate calculated by an atmospheric general circulation model (GCM). We then hold this artificial rheology constant in further numerical experiments. We hope that in doing so we can partially compensate for systematic GCM errors in simulating the snow accumulation rate, and, therefore, thickening/thinning of ice sheets will depend mostly on the tendency in the net accumulation change rather than on its absolute value. Using this approach, the response of the Greenland and Antarctic ice sheets to doubling CO₂ concentration is simulated and the horizontal distribution of possible thickening/thinning of polar ice obtained. We find that, initially, the region of thickening ice is close to the area of increased snowfall rate, but later it significantly changes under the influence of internal ice flow dynamics. The sea-level changes predicted by our experiments agree with some empirical estimates. The sensitivity experiment with assigned basal sliding does not show significant changes in the large-scale ice topography, meaning, for example, that there is no indication of a possible disintegration of the West Antarctic ice sheet. At the same time, the regional thickening/thinning of ice (and consequently the sea-level change) depends strongly on processes at the ice sheet bottom.     

Greenland��s contribution to global sea-level rise by the end of the 21st century

The Greenland ice sheet holds enough water to raise the global sea level with ??7 m. Over the last few decades, observations manifest a substantial increase of the mass loss of this ice sheet. Both enhanced melting and increase of the dynamical discharge, associated with calving at the outlet-glacier fronts, are contributing to the mass imbalance. Using a dynamical and thermodynamical ice-sheet model, and taking into account speed up of outlet glaciers, we estimate Greenland??s contribution to the 21st-century global sea-level rise and the uncertainty of this estimate. Boundary fields of temperature and precipitation extracted from coupled climate-model projections used for the IPCC Fourth Assessment Report, are applied to the ice-sheet model. We implement a simple parameterization for increased flow of outlet glaciers, which decreases the bias of the modeled present-day surface height. It also allows for taking into account the observed recent increase in dynamical discharge, and it can be used for future projections associated with outlet-glacier speed up. Greenland contributes 0?C17?cm to global sea-level rise by the end of the 21st century. This range includes the uncertainties in climate-model projections, the uncertainty associated with scenarios of greenhouse-gas emissions, as well as the uncertainties in future outlet-glacier discharge. In addition, the range takes into account the uncertainty of the ice-sheet model and its boundary fields.     

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.     

Simulation of glacial-interglacial cycles by simple relaxation models: consistency with observational results

The model of Paillard and Parrenin (Earth Planet Sci Lett 227:263–271, 2004) was modified to obtain a closer fit to δ¹⁸O and CO₂ time series for the last 800 kyr. The model performance can be improved if its CO₂ sensitivity to I₆₅ insolation is eliminated and if different response times are assumed for ablation/accumulation of ice. Correlations between simulated and experimental time series for CO₂ and ice volume V increase from 0.59 and 0.63 to 0.79 and 0.88, respectively. According to these models, terminations are produced by I₆₅ amplification through CO₂-T and T-CO₂ feedbacks, in synergy with an extra CO₂ contribution from the deep ocean. This contribution is strongly dependent on ice-sheet extent and ice volume (or alternatively, CO₂ concentration, which is a good proxy of Antarctic temperature) but is insensitive to Southern Ocean (SO) insolation on 21 February (I₆₀). Change of deep SO state may be the “order parameter” for nonlinear deglacial changes. According to these models, 100 kyr periodicity of glacial cycles arises from the characteristic time of Antarctic ice sheet advance to the continental slope.     

A 3-D model for the Antarctic ice sheet: a sensitivity study on the glacial-interglacial contrast

On the longer climatic time scales, changes in the elevation and extent of the Antarctic ice sheet have an important role in modulating global atmospheric and oceanographic processes, and contribute significantly to world-wide sea levels. In this paper, a 3-D time-dependent thermomechanical model for the entire ice sheet is presented, that is subsequently used to examine the effects of glacial-interglacial shifts in environmental boundary conditions on its geometry. The model takes into account a coupled ice shelf, grounding-line dynamics, basal sliding and isostatic bed adjustment and considers the fully coupled velocity and temperature fields. Ice flow is calculated on a fine mesh (40 km horizontal grid size and 10 layers in the vertical) for grounded and floating ice and a stress transition zone in between at the grounding line, where all stress components contribute in the effective stress in the flow law. There is free interaction between ice sheet and ice shelf, so that the entire geometry is internally generated. A simulation of the present ice sheet reveals that the model is able to yield realistic results. A series of sensitivity experiments are then performed, in which lower temperatures, reduced accumulation rates and lower global sea level stands are imposed, either singly or in combination. By comparing results of pairs of experiments, the effects of each of these environmental changes can be determined. In agreement with glacial-geological evidence, we found that the most pronounced changes show up in the West Antarctic ice sheet configuration. They appear to be essentially controlled by variations in eustatic sea level, whereas typical glacial-interglacial changes in temperature and ice deposition rates tend to balance one another. These findings support the hypothesis that the Antarctic ice sheet basically follows glacial episodes in the northern hemisphere by means of sea-level teleconnections. Grounding occurs more readily in the Weddell sea than in the Ross sea and long time scales appear to be involved: it may take up to 30–40000 years for these continental shelf areas to become completely grounded after an initial stepwise perturbation in boundary conditions. According to these reconstructions, a steady state Antarctic ice sheet may contribute some 16 m to global sea level lowering at maximum glaciation.     

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.     

一个新的干沉积速度型及其在高分辨区域酸沉积模式中的应用

为研究区域尺度干沉积过程，发展了一新的干沉积速度型（ＮＤＤＶＰ）。ＮＤＤＶＰ与１０２２次于沉积速度（Ｖｄ）实测资料之间平均比值为１．０６±０．８２，表明ＮＤＤＶＰ与实测资料有好的一致性。用分别具有ＮＤＤＶＰ和原Ｖｄ型的高分进区域酸沉积模式（ＥＭ３）作实例预测，发现仅由Ｖｄ型的不同所引起的ＳＯ＝４中心浓度之间最大比值可达２．４，用ＮＤＤＶＰ预测的物质浓度和干沉积量比用原Ｖｄ型有明显的改善。     

On the Response of the Greenland Ice Sheet to Greenhouse Climate Change

Numerical computations are performed with the three-dimensional polythermal ice-sheet model SICOPOLIS in order to investigate the possible impact of a greenhouse-gas-induced climate change on the Greenland ice sheet. The assumed increase of the mean annual air temperature above the ice covers a range from T = 1°C to 12°C, and several parameterizations for the snowfall and the surface melting are considered. The simulated shrinking of the ice sheet is a smooth function of the temperature rise, indications for the existence of critical thresholds of the climate input are not found. Within 1000 model years, the ice-volume decrease is limited to 10% of the present volume for T 3°C, whereas the most extreme scenario, T = 12°C, leads to an almost entire disintegration, which corresponds to a sea-level equivalent of 7 m. The different snowfall and melting parameterizations yield an uncertainty range of up to 20% of the present ice volume after 1000 model years.     

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.     

An investigation of the small ice cap instability in the Southern Hemisphere with a coupled atmosphere-sea ice-ocean-terrestrial ice model

 A simple climate model has been developed to investigate the existence of the small ice cap instability in the Southern Hemisphere. The model consists of four coupled components: an atmospheric energy balance model, a thermodynamic snow-sea ice model, an oceanic mixed layer model and a terrestrial ice model. Results from a series of experiments involving different degrees of coupling in the model show that the instability appears only in those cases when an explicit representation of the Antarctic ice sheet is not included in the model. In order to determine which physical processes in the ice sheet model lead to a stabilization of the system we have conducted several sensitivity experiments in each of which a given ice sheet process has been removed from the control formulation of the model. Results from these experiments suggest that the feedback between the elevation of the ice sheet and the snow accumulation-ice ablation balance is responsible for the disappearance of the small ice cap instability in our simulation. In the model, the mass balance of the ice sheet depends on the air temperature at sea level corrected for altitude and it is, therefore, a function of surface elevation. This altitude-mass balance feedback effectively decouples the location of the ice edge from any specific sea level isotherm, thus decreasing the model sensitivity to the albedo-temperature feedback, which is responsible for the appearance of the instability. It is also shown that the elevation-radiative cooling feedback tends to stabilize the ice sheet, although its effect does not seem to be strong enough to remove the instability. Another interesting result is that for those simulations which include the terrestrial ice model with elevation-dependent surface mass balance, hysteresis is exhibited, where for a given level of external forcing, two stable solutions with different, non-zero ice-sheet volume and area and different air and ocean temperature fields occur. However, no unstable transition between the two solutions is ever observed. Our results suggest that the small ice cap instability mechanism could be unsuitable for explaining the inception of glaciation in Antarctica. Received: 14 April 1997 / Accepted: 22 October 1997     

Modelling the short-term response of the Greenland ice-sheet to global warming

 A two-dimensional vertically integrated ice flow model has been developed to test the importance of various processes and concepts used for the prediction of the contribution of the Greenland ice-sheet to sea-level rise over the next 350 y (short-term response). The mass balance is modelled by the degree-day method and the energy-balance method. The lithosphere is considered to respond isostatically to a point load and the time evolution of the bedrock follows from a viscous asthenosphere. According to the IPCC-IS92a scenario (with constant aerosols after 1990) the Greenland ice-sheet is likely to cause a global sea level rise of 10.4 cm by 2100 AD. It is shown, however, that the result is sensitive to precise model formulations and that simplifications as used in the sea-level projection in the IPCC-96 report yield less accurate results. Our model results indicate that, on a time scale of a hundred years, including the dynamic response of the ice-sheet yields more mass loss than the fixed response in which changes in geometry are not incorporated. It appears to be important to consider sliding, as well as the fact that climate sensitivity increases for larger perturbations. Variations in predicted sea-level change on a time scale of hundred years depend mostly on the initial state of the ice-sheet. On a time scale of a few hundred years, however, the variability in the predicted melt is dominated by the variability in the climate scenarios. Received: 21 August 1996/Accepted: 12 May 1997     

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.     

Recent precipitation trends over the polar ice sheets

Summary Meteorological and glaciological analyses are integrated to examine the precipitation trends during the last three decades over the ice sheets covering Antarctica and Greenland. For Antarctica, the best data source is provided by glaciologically-measured trends of snow accumulation, and for limited sectors of East Antarctica consistency with precipitation amounts calculated from the atmospheric water balance equation is obtained. For Greenland, precipitation rates parameterized from atmospheric analyses yield the only comprehensive depiction. The precipitation rate over Antarctica appears to have increased by about 5% over a time period spanning the accumulation means for the 1955–65 to 1965–75 periods, while over Greenland it has decreased by about 15% since 1983 with a secondary increase over the southern part of the ice sheet starting in 1977. At the end of the 10-year overlapping period, the global sea-level impact of the precipitation changes over Antarctica dominates that for Greenland and yields a net ice-sheet precipitation contribution of roughly 0.02 mm yr^–1. These changes are likely due to marked variations in the cyclonic forcing affecting the ice sheets, but are only weakly reflected in the temperature regime, consistent with the episodic nature of cyclonic precipitation. These conclusions are not founded on high quality data bases. The importance of such changes for understanding global sea-level variations argues for a modest research effort to collect simultaneous meteorological and glaciological observations in order to describe and understand the current precipitation variations over both ice sheets. Some suggestions are offered for steps that could be taken.With 8 Figures     

A one-dimensional model for the interaction between continental-scale ice sheets and atmospheric stationary waves

 The great continental ice sheets of the Pleistocene represented a significant topographic obstacle to the westerly winds in northern midlatitudes. This work explores how consequent changes in the atmospheric stationary wave pattern might have affected the shape and growth of the ice sheets themselves. A one dimensional (1-D) model is developed which permits an examination of the types and magnitudes of the feedbacks that might be expected. When plausible temperature perturbations are introduced at the ice-sheet margin which are proportional to the stationary wave amplitude, the equilibrium shape of the ice sheet is significantly altered, and depends on the sign of the perturbation. The proposed feedback also affects the response of the ice sheet to time-varying climate forcing. The results suggest that the evolution of a continental-scale ice sheet with a land-based margin may be significantly determined by the changes it induces in the atmospheric circulation. Received: 1 October 1999 / Accepted: 17 July 2000     

Projecting twenty-first century regional sea-level changes

We present regional sea-level projections and associated uncertainty estimates for the end of the 21 ^st century. We show regional projections of sea-level change resulting from changing ocean circulation, increased heat uptake and atmospheric pressure in CMIP5 climate models. These are combined with model- and observation-based regional contributions of land ice, groundwater depletion and glacial isostatic adjustment, including gravitational effects due to mass redistribution. A moderate and a warmer climate change scenario are considered, yielding a global mean sea-level rise of 0.54 ±0.19 m and 0.71 ±0.28 m respectively (mean ±1σ). Regionally however, changes reach up to 30 % higher in coastal regions along the North Atlantic Ocean and along the Antarctic Circumpolar Current, and up to 20 % higher in the subtropical and equatorial regions, confirming patterns found in previous studies. Only 50 % of the global mean value is projected for the subpolar North Atlantic Ocean, the Arctic Ocean and off the western Antarctic coast. Uncertainty estimates for each component demonstrate that the land ice contribution dominates the total uncertainty.     

Response of the Antarctic ice sheet to future greenhouse warming

Possible future changes in land ice volume are mentioned frequently as an important aspect of the greenhouse problem. This paper deals with the response of the Antarctic ice sheet and presents a tentative projection of changes in global sea level for the next few hundred years, due to changes in its surface mass balance. We imposed a temperature scenario, in which surface air temperature rises to 4.2° C in the year 2100 AD and is kept constant afterwards. As GCM studies seem to indicate a higher temperature increase in polar latitudes, the response to a more extreme scenario (warming doubled) has also been investigated. The mass balance model, driven by these temperature perturbations, consists of two parts: the accumulation rate is derived from present observed values and is consequently perturbed in proportion to the saturated vapour pressure at the temperature above the inversion layer. The ablation model is based on the degree-day method. It accounts for the daily temperature cycle, uses a different degree-day factor for snow and ice melting and treats refreezing of melt water in a simple way. According to this mass balance model, the amount of accumulation over the entire ice sheet is presently 24.06 × 10¹¹ m³ of ice, and no runoff takes place. A 1°C uniform warming is then calculated to increase the overall mass balance by an amount of 1.43 × 10¹¹ m³ of ice, corresponding to a lowering of global sea level with 0.36 mm/yr. A temperature increase of 5.3°C is needed for the increase in ablation to become more important than the increase in accumulation and the temperature would have to rise by as much as 11.4°C to produce a zero surface mass balance. Imposing the Bellagio-scenario and accumulating changes in mass balance forward in time (static response) would then lower global sea level by 9 cm by 2100 AD. In a subsequent run with a high-resolution 3-D thermomechanic model of the ice sheet, it turns out that the dynamic response of the ice sheet (as compared to the direct effect of the changes in surface mass balance) becomes significant after 100 years or so. Ice-discharge across the grounding-line increases, and eventually leads to grounding-line retreat. This is particularly evident in the extreme case scenario and is important along the Antarctic Peninsula and the overdeepened outlet glaciers along the East Antarctic coast. Grounding-line retreat in the Ross and Ronne-Filchner ice shelves, on the other hand, is small or absent.     

Influence of Antarctic ice sheet lowering on the Southern Hemisphere climate: modeling experiments mimicking the mid-Miocene

A coupled global atmosphere-ocean model is used to study the influence of the Antarctica ice sheet in a configuration that mimics that of the early Miocene on the atmospheric and oceanic circulations. Based on different climate simulations of the present day (CTR) and conducted with distinct Antarctic ice sheet topography (AIS-EXP), it is found that the reduction of the Antarctic ice sheet topography (AIS) induces warming of the Southern Hemisphere and reduces the meridional thermal gradient. Consequently, the atmospheric transient low level eddy heat flux $[(\overline{v^{\prime}T^{\prime}})]$ and the eddy momentum flux $[(\overline{u^{\prime}v^{\prime}})]$ are reduced causing the reduced transport of heat from the mid-latitudes to the pole. The stationary flow and transient wave anomalies generate changes in the SSTs which modify the rate of deep water formation, strengthening the formation of the Antarctic Bottom Water. Substantial changes are predicted to occur in the atmospheric and oceanic heat transport and a comparison between the total heat transport of the atmosphere-ocean system, as simulated by the AIS-EXP and the CTR runs, shows that the reduction of the AIS height leads to reduced Southern Hemisphere poleward and increased equatorward heat transport. These results are in agreement with reduced storm track activities and baroclinicity.     

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.     

Simulations of the Last Glacial Maximum climates using a general circulation model: prescribed versus computed sea surface temperatures

 The climate during the Last Glacial Maximum (LGM) has been simulated using the UK Universities Global Atmospheric Modelling Programme (UGAMP) general circulation model (GCM) with both prescribed sea surface temperatures (SSTs) based on the CLIMAP reconstruction and computed SSTs with a simple thermodynamic slab ocean. Consistent with the Paleoclimate Modelling Intercomparison Project (PMIP), the other boundary conditions include the large changes in ice-sheet topography and geography, a lower sea level, a lower concentration of CO₂ in the atmosphere, and a slightly different insolation pattern at the top of the atmosphere. The results are analysed in terms of changes in atmospheric circulation. Emphasis is given to the changes in surface temperatures, planetary waves, storm tracks and the associated changes in distribution of precipitation. The model responds in a similar manner to the changes in boundary conditions to previous studies in global mean statistics, but differs in its treatment of regional climates. Results also suggest that both the land ice sheets and sea ice introduce significant changes in planetary waves and transient eddy activity, which in turn affect regional climates. The computed SST simulations predict less sea ice and cooler tropical temperatures than those based on CLIMAP SSTs. It is unclear as to whether this is a model and/or a data problem, but the resulting changes in land temperatures and precipitation can be large. Snow mass budget analysis suggests that there is net ice loss along the southern edges of the Laurentide and Fennoscandian ice sheets and net ice gain over other parts of the two ice sheets. The net accumulation is mainly due to the decrease in ablation in the cold climate rather than to the changes in snowfall. The characteristics of the Greenland ice-sheet mass balance in the LGM simulations is also quite different from those in the present-day (PD) simulations. The ablation in the LGM simulations is negligible while it is a very important process in the ice mass budget in the PD simulations. Received: 10 January 1997 / Accepted: 11 December 1997