▪ Abstract  The 100 kyr quasiperiodic variation of continental ice cover, which has been a persistent feature of climate system evolution throughout the most recent 900 kyr of Earth history, has occurred as a consequence of changes in the seasonal insolation regime forced by the influence of gravitational n-body effects in the Solar System on the geometry of Earth's orbit around the Sun. The impacts of the changing surface ice load upon both Earth's shape and gravitational field, as well as upon sea-level history, have come to be measurable using a variety of geological and geophysical techniques. These observations are invertible to obtain useful information on both the internal viscoelastic structure of the solid Earth and on the detailed spatiotemporal characteristics of glaciation history. This review focuses upon the most recent advances that have been achieved in each of these areas, advances that have proven to be central to the construction of the refined model of the global process of glacial isostatic adjustment, denoted ICE-5G (VM2). A significant test of this new global model will be provided by the global measurement of the time dependence of the gravity field of the planet that will be delivered by the GRACE satellite system that is now in space.

A new model of the last deglaciation event of the Late Quaternary ice age is here described and denoted as ICE-6G_C (VM5a). It differs from previously published models in this sequence in that it has been explicitly refined by applying all of the available Global Positioning System (GPS) measurements of vertical motion of the crust that may be brought to bear to constrain the thickness of local ice cover as well as the timing of its removal. Additional space geodetic constraints have also been applied to specify the reference frame within which the GPS data are described. The focus of the paper is upon the three main regions of Last Glacial Maximum ice cover, namely, North America, Northwestern Europe/Eurasia, and Antarctica, although Greenland and the British Isles will also be included, if peripherally, in the discussion. In each of the three major regions, the model predictions of the time rate of change of the gravitational field are also compared to that being measured by the Gravity Recovery and Climate Experiment satellites as an independent means of verifying the improvement of the model achieved by applying the GPS constraints. Several aspects of the global characteristics of this new model are also discussed, including the nature of relative sea level history predictions at far-field locations, in particular the Caribbean island of Barbados, from which especially high-quality records of postglacial sea level change are available but which records were not employed in the development of the model. Although ICE-6G_C (VM5a) is a significant improvement insofar as the most recently available GPS observations are concerned, comparison of model predictions with such far-field relative sea level histories enables us to identify a series of additional improvements that should follow from a further stage of model iteration.

Ocean Drilling Program (ODP) Site 677 provided excellent material for high resolution stable isotope analysis of both benthonic and planktonic foraminifera through the entire Pleistocene and upper Pliocene. The oxygen isotope record is readily correlated with the SPECMAP stack (Imbrie et al. 1984) and with the record from DSDP 607 (Ruddiman et al. 1986) but a significantly better match with orbital models is obtained by departing from the timescale proposed by these authors below Stage 16 (620 000 years). It is the stronger contribution from the precession signal in the record from ODP Site 677 that provides the basis for the revised timescale. Our proposed modification to the timescale would imply that the currently adopted radiometric dates for the Matuyama–Brunhes boundary, the Jaramillo and Olduvai Subchrons and the Gauss–Matuyama boundary underestimate their true astronomical ages by between 5 and 7%.

The Atlantic thermohaline circulation (THC) is an important part of the earth's climate system. Previous research has shown large uncertainties in simulating future changes in this critical system. The simulated THC response to idealized freshwater perturbations and the associated climate changes have been intercompared as an activity of World Climate Research Program (WCRP) Coupled Model Intercomparison Project/Paleo-Modeling Intercomparison Project (CMIP/PMIP) committees. This intercomparison among models ranging from the earth system models of intermediate complexity (EMICs) to the fully coupled atmosphere-ocean general circulation models (AOGCMs) seeks to document and improve understanding of the causes of the wide variations in the modeled THC response. The robustness of particular simulation features has been evaluated across the model results. In response to 0.1-Sv (1 Sv ϵ 10 6 m 3 s Ϫ1) freshwater input in the northern North Atlantic, the multimodel ensemble mean THC weakens by 30% after 100 yr. All models simulate some weakening of the THC, but no model simulates a complete shutdown of the THC. The multimodel ensemble indicates that the surface air temperature could present a complex anomaly pattern with cooling south of Greenland and warming over the Barents and Nordic Seas. The Atlantic ITCZ tends to shift southward. In response to 1.0-Sv freshwater input, the THC switches off rapidly in all model simulations. A large cooling occurs over the North Atlantic. The annual mean Atlantic ITCZ moves into the Southern Hemisphere. Models disagree in terms of the reversibility of the THC after its shutdown. In general, the EMICs and AOGCMs obtain similar THC responses and climate changes with more pronounced and sharper patterns in the AOGCMs.

A gravitationally self-consistent theory of postglacial relative sea level change is used to infer the variation of surface ice and water cover since the Last Glacial Maximum (LGM). The results show that LGM ice volume was approximately 35 percent lower than suggested by the CLIMAP reconstruction and the maximum heights of the main Laurentian and Fennoscandian ice complexes are inferred to have been commensurately lower with respect to sea level. Use of these Ice Age boundary conditions in atmospheric general circulation models will yield climates that differ significantly from those previously inferred on the basis of the CLIMAP data set.