Book review

Summary and Conclusions The Last Great Ice Sheets provides a valuable compilation for current estimates of the extent and timing of the maximum ice extent of the last ice age period around the globe. The areas of uncertainty are indicated and emphasis is given to points of controversy where further research is needed to resolve the most important problems still remaining. The work thereby presents plausible upper and lower limit estimates for the maximum extent of the ice corresponding to about 18 ka BP. These limits can be used for the boundary conditions required for atmospheric climate modelling studies. For this type of application the differences in the upper and lower limits are not serious and consequently the possible errors in the reconstruction estimates are not so important. The greatest uncertainties in the ice cover reconstructions occur for the northerly limits of the grounded ice which tend to be in off-shore regions where reliable data is sparse. This signals the requirement for a much greater research effort to collect off-shore sediment sequences and other data to help clarify the ice cover fluctuation record. The possibility of an extensive arctic ice shelf system and marine ice sheet cover interconnecting with the various grounded ice sheets is put forward as a working hypothesis along with arguments in favour of the upper limit estimates for the maximum ice extent. The extent of the ice age ice shelves is an important question which should be examined by dynamic ice sheet modelling with explicit ice shelf formulation. The extensive growth of the ice sheets to seaward margins leads naturally to ice shelf formation. The major questions are: how extensive were the ice shelves, and to what extent did coalescence occur? The further hypothesis that the ice shelves may have had a major role in the growth and decay of the ice sheets is more difficult to support since the results of the dynamic modelling indicate that the ice shelves form readily as a result of extensive growth of the ice sheets but it is difficult to start extensive ice growth near sea level. The information presented for the timing of the advance towards the maximum ice cover and the subsequent pattern of retreat provides further valuable material to test the dynamic ice sheet and climate models including the reactions to the Earth's orbital radiation changes. These results for the rates of change of the ice cover taken together with the modelling results of Budd and Smith indicate that neither the ice sheets nor the bedrock depression had time to reach equilibrium states. The non-equilibrium nature of the ice sheets with only short periods between relatively rapid advance and retreat phases is an important factor that needs to be taken into account in assessing the reconstructions of the ice sheets. This emphasises the need for further work with dynamic ice sheet models and coupled global atmosphere ocean models to determine more clearly the sequence of changes during the growth and decay of the large ice sheets. ‘The Last Great Ice Sheets’ provides a timely data base and compilation to support these studies.     

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     

Radiative forcing and response of a GCM to ice age boundary conditions: cloud feedback and climate sensitivity

 A general circulation model is used to examine the effects of reduced atmospheric CO₂, insolation changes and an updated reconstruction of the continental ice sheets at the Last Glacial Maximum (LGM). A set of experiments is performed to estimate the radiative forcing from each of the boundary conditions. These calculations are used to estimate a total radiative forcing for the climate of the LGM. The response of the general circulation model to the forcing from each of the changed boundary conditions is then investigated. About two-thirds of the simulated glacial cooling is due to the presence of the continental ice sheets. The effect of the cloud feedback is substantially modified where there are large changes to surface albedo. Finally, the climate sensitivity is estimated based on the global mean LGM radiative forcing and temperature response, and is compared to the climate sensitivity calculated from equilibrium experiments with atmospheric CO₂ doubled from present day concentration. The calculations here using the model and palaeodata support a climate sensitivity of about 1 Wm^-2 K^-1 which is within the conventional range. Received: 8 February 1997 / Accepted: 4 June 1997     

Two climatic states and feedbacks on thermohaline circulation in an Earth system model of intermediate complexity

The McGill Paleoclimate Model-2 (MPM-2) is employed to study climate–thermohaline circulation (THC) interactions in a pre -industrial climate, with a special focus on the feedbacks on the THC from other climate system components. The MPM-2, a new version of the MPM, has an extended model domain from 90S to 90N, active winds and no oceanic heat and freshwater flux adjustments. In the MPM-2, there are mainly two stable modes for the Atlantic meridional overturning circulation (MOC) under the ‘present-day’ forcing (present-day solar forcing and the pre-industrial atmospheric CO₂ level of 280 ppm). The ‘on’ mode has an active North Atlantic deep water formation, while the ‘off’ mode has no such deep water formation. By comparing the ‘off’ mode climate state with its ‘on’ mode analogue, we find that there exist many large differences between the two climate states, which originate from large changes in the oceanic meridional heat transports. By suppressing or isolating each process associated with a continental ice sheet over North America, sea ice, the atmospheric hydrological cycle and vegetation, feedbacks from these components on the Atlantic MOC are investigated. Sensitivity studies investigating the role of varying continental ice growth and sea ice meridional transport in the resumption of the Atlantic MOC are also carried out. The results show that a fast ice sheet growth and an enhanced southward sea ice transport significantly favor the resumption of the Atlantic MOC in the MPM-2. In contrast to this, the feedback from the atmospheric hydrological cycle is a weak positive one. The vegetation-albedo feedback could enhance continental ice sheet growth and thus could also favor the resumption of the Atlantic MOC. However, before the shut-down of the Atlantic MOC, feedbacks from these components on the Atlantic MOC are very weak.     

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.     

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     

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.     

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.     

Relative impacts of insolation changes, meltwater fluxes and ice sheets on African and Asian monsoons during the Holocene

In order to better understand the evolution of the Afro-Asian monsoon in the early Holocene, we investigate the impact on boreal summer monsoon characteristics of (1) a freshwater flux in the North Atlantic from the surrounding melting ice sheets and (2) a remnant ice sheet over North America and Europe. Sensitivity experiments run with the IPSL_CM4 model show that both the meltwater flux and the remnant ice sheets induce a cooling of similar amplitude of the North Atlantic leading to a southward shift of the Inter-Tropical Convergence Zone over the tropical Atlantic and to a reduction of the African monsoon. The two perturbations have different impacts in the Asian sector. The meltwater flux results in a weakening of the Indian monsoon and no change in the East Asian monsoon, whereas the remnant ice sheets induce a strengthening of the Indian monsoon and a strong weakening of the East Asian monsoon. Despite the similar coolings in the Atlantic Ocean, the ocean heat transport is reduced only in the meltwater flux experiment, which induces slight differences between the two experiments in the role of the surface latent heat flux in the tropical energetics. In the meltwater experiment, the southward shift of the subtropical jet acts to cool the upper atmosphere over the Tibetan Plateau and hence to weaken the Indian monsoon. In the ice sheet experiment this effect is overwhelmed by the changes in extratropical stationary waves induced by the ice sheets, which are associated with a larger cooling over the Eurasian continent than in the meltwater experiment. However these sensitivity experiments suggest that insolation is the dominant factor explaining the relative changes of the African, Indian and East Asian monsoons from the early to the mid-Holocene.     

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.     

Impact of a realistic river routing in coupled ocean–atmosphere simulations of the Last Glacial Maximum climate

The presence of large ice sheets over North America and North Europe at the Last Glacial Maximum (LGM) strongly impacted Northern hemisphere river pathways. Despite the fact that such changes may significantly alter the freshwater input to the ocean, modified surface hydrology has never been accounted for in coupled ocean–atmosphere general circulation model simulations of the LGM climate. To reconstruct the LGM river routing, we use the ICE-5G LGM topography. Because of the uncertainties in the extent of the Fennoscandian ice sheet in the Eastern part of the Kara Sea, we consider two more realistic river routing scenarios. The first scenario is characterised by the presence of an ice dammed lake south of the Fennoscandian ice sheet, and corresponds to the ICE-5G topography. This lake is fed by the Ob and Yenisei rivers. In the second scenario, both these rivers flow directly into the Arctic Ocean, which is more consistent with the latest QUEEN ice sheet margin reconstructions. We study the impact of these changes on the LGM climate as simulated by the IPSL_CM4 model and focus on the overturning thermohaline circulation. A comparison with a classical LGM simulation performed using the same model and modern river basins as designed in the PMIP2 exercise leads to the following conclusions: (1) The discharge into the North Atlantic Ocean is increased by 2,000 m³/s between 38° and 54°N in both simulations that contain LGM river routing, compared to the classical LGM experiment. (2) The ice dammed lake is shown to have a weak impact, relative to the classical simulation, both in terms of climate and ocean circulation. (3) In contrast, the North Atlantic deep convection and meridional overturning are weaker than during the classical LGM run if the Ob and Yenisei rivers flow directly into the Arctic Ocean. The total discharge into the Arctic Ocean is increased by 31,000 m³/s, relative to the classical LGM simulation. Consequentially, northward ocean heat transport is weaker, and sea ice more extensive, in better agreement with existing proxy data.     

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.     

Closure of the Panama Seaway during the Pliocene: implications for climate and Northern Hemisphere glaciation

The “Panama Hypothesis” states that the gradual closure of the Panama Seaway, between 13 million years ago (13 Ma) and 2.6 Ma, led to decreased mixing of Atlantic and Pacific water Masses, the formation of North Atlantic Deep water and strengthening of the Atlantic thermohaline circulation, increased temperatures and evaporation in the North Atlantic, increased precipitation in Northern Hemisphere (NH) high latitudes, culminating in the intensification of Northern Hemisphere Glaciation (NHG) during the Pliocene, 3.2–2.7 Ma. Here we test this hypothesis using a fully coupled, fully dynamic ocean-atmosphere general circulation model (GCM) with boundary conditions specific to the Pliocene, and a high resolution dynamic ice sheet model. We carry out two GCM simulations with “closed” and “open” Panama Seaways, and use the simulated climatologies to force the ice sheet model. We find that the models support the “Panama Hypothesis” in as much as the closure of the seaway results in a more intense Atlantic thermohaline circulation, enhanced precipitation over Greenland and North America, and ultimately larger ice sheets. However, the volume difference between the ice sheets in the “closed” and “open” configurations is small, equivalent to about 5 cm of sea level. We conclude that although the closure of the Panama Seaway may have slightly enhanced or advanced the onset of NHG, it was not a major forcing mechanism. Future work must fully couple the ice sheet model and GCM, and investigate the role of orbital and CO₂ effects in controlling NHG.     

The influence of continental ice,atmospheric CO2, and land albedo on the climate of the last glacial maximum

The contributions of expanded continental ice, reduced atmospheric CO₂, and changes in land albedo to the maintenance of the climate of the last glacial maximum (LGM) are examined. A series of experiments is performed using an atmosphere-mixed layer ocean model in which these changes in boundary conditions are incorporated either singly or in combination. The model used has been shown to produce a reasonably realistic simulation of the reduced temperature of the LGM (Manabe and Broccoli 1985b). By comparing the results from pairs of experiments, the effects of each of these environmental changes can be determined.Expanded continental ice and reduced atmospheric CO₂ are found to have a substantial impact on global mean temperature. The ice sheet effect is confined almost exclusively to the Northern Hemisphere, while lowered CO₂ cools both hemispheres. Changes in land albedo over ice-free areas have only a minor thermal effect on a global basis. The reduction of CO₂ content in the atmosphere is the primary contributor to the cooling of the Southern Hemisphere. The model sensitivity to both the ice sheet and CO₂ effects is characterized by a high latitude amplification and a late autumn and early winter maximum.Substantial changes in Northern Hemisphere tropospheric circulation are found in response to LGM boundary conditions during winter. An amplified flow pattern and enhanced westerlies occur in the vicinity of the North American and Eurasian ice sheets. These alterations of the tropospheric circulation are primarily the result of the ice sheet effect, with reduced CO₂ contributing only a slight amplification of the ice sheet-induced pattern.     

Long-term effects of anthropogenic CO<Subscript>2</Subscript> emissions simulated with a complex earth system model

A new complex earth system model consisting of an atmospheric general circulation model, an ocean general circulation model, a three-dimensional ice sheet model, a marine biogeochemistry model, and a dynamic vegetation model was used to study the long-term response to anthropogenic carbon emissions. The prescribed emissions follow estimates of past emissions for the period 1751–2000 and standard IPCC emission scenarios up to the year 2100. After 2100, an exponential decrease of the emissions was assumed. For each of the scenarios, a small ensemble of simulations was carried out. The North Atlantic overturning collapsed in the high emission scenario (A2) simulations. In the low emission scenario (B1), only a temporary weakening of the deep water formation in the North Atlantic is predicted. The moderate emission scenario (A1B) brings the system close to its bifurcation point, with three out of five runs leading to a collapsed North Atlantic overturning circulation. The atmospheric moisture transport predominantly contributes to the collapse of the deep water formation. In the simulations with collapsed deep water formation in the North Atlantic a substantial cooling over parts of the North Atlantic is simulated. Anthropogenic climate change substantially reduces the ability of land and ocean to sequester anthropogenic carbon. The simulated effect of a collapse of the deep water formation in the North Atlantic on the atmospheric CO₂ concentration turned out to be relatively small. The volume of the Greenland ice sheet is reduced, but its contribution to global mean sea level is almost counterbalanced by the growth of the Antarctic ice sheet due to enhanced snowfall. The modifications of the high latitude freshwater input due to the simulated changes in mass balance of the ice sheet are one order of magnitude smaller than the changes due to atmospheric moisture transport. After the year 3000, the global mean surface temperature is predicted to be almost constant due to the compensating effects of decreasing atmospheric CO₂ concentrations due to oceanic uptake and delayed response to increasing atmospheric CO₂ concentrations before.     

The Pangean ice age: studies with a coupled climate-ice sheet model

 Application of an ice sheet model developed for the Pleistocene to the extensive Carboniferous glaciation on Gondwana yields an ice sheet which has several features consistent with observations. While complete deglaciation is not achieved without CO₂ changes, the Milankovich-induced fluctuations in ice sheet volume are comparable to Pleistocene glacial/ interglacial signals. This result is shown to hold for a large fraction of physically reasonable parameter space. The model also exhibits multiple equilibria and sharp bifurcations, as infinitesimal changes in the solar constant or precipitation can lead to a qualitatively different climate. The success of the model in predicting ice location in an environment quite different from the Pleistocene provides additional support for the robustness of the basic model physics and suggests that the model can be applied with some confidence to other pre-Pleistocene glaciations. Received: 30 June 1998 / Accepted: 5 January 1999     

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.     

Impact of ice sheet induced North Atlantic oscillation on East Asian summer monsoon during an interglacial 500,000?years ago

Marine Isotope Stage (MIS) 13, an interglacial about 500,000?years ago, is unique due to an exceptionally strong East Asia summer monsoon (EASM) occurring in a relatively cool climate with low greenhouse gas concentrations (GHG). This paper attempts to find one of the possible mechanisms for this seeming paradox. Simulations with an Earth System model LOVECLIM show that the presence of ice sheets over North America and Eurasia during MIS-13 induces a positive phase of the winter North Atlantic Oscillation (NAO) like feature. The ocean having a longer memory than the atmosphere, the oceanic anomalies associated with NAO persists until summer. The signals of summer NAO are transmitted to East Asia to reinforce the monsoon there through the stationary waves excited at the Asian Jet entrance. The geopotential height shows clearly a mid-latitude wave train with positive anomalies over the eastern Mediterranean/Caspian Sea and the Okhotsk Sea and a negative anomaly over Lake Baikal. This reinforces the effect of the high-latitude wave train induced independently by the Eurasian ice sheet topography as shown in previous study. These features reinforce the Meiyu front and enhance the precipitation over East Asia. The results obtained from LOVECLIM are further confirmed by an atmospheric general circulation model, ARPEGE.     

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.