A coupled climate model simulation of the Last Glacial Maximum, Part 2: approach to equilibrium

Abstract.The climate of the last glacial maximum (LGM) is simulated with a coupled climate model. The simulated climate undergoes a rapid adjustment during the first several decades after imposition of LGM boundary conditions, as described in Part 1, and then evolves toward equilibrium over 900 model years. The climate simulated by the coupled model at this period is compared with observationally-based LGM reconstructions and with LGM results obtained with an atmosphere-mixed layer (slab) ocean version of the model in order to investigate the role of ocean dynamics in the LGM climate. Global mean surface air temperature and sea surface temperature (SST) decrease by about 10 °C and 5.6 °C in the coupled model which includes ocean dynamics, compared to decreases of 6.3 and 3.8 °C in slab ocean case. The coupled model simulates a cooling of about 6.5 °C over the tropics, which is larger than that of the CLIMAP reconstruction (1.7 °C) and larger than that of the slab ocean simulation (3.3 °C), but which is in reasonable agreement with some recent proxy estimates. The ocean dynamics of the coupled model captures features found in the CLIMAP reconstructions such as a relative maximum of ocean cooling over the tropical Pacific associated with a mean La Niña-like response and lead to a more realistic SST pattern than in the slab model case. The reduction in global mean precipitation simulated in the coupled model is larger (15%) than that simulated with the slab ocean model (∼10%) in conjunction with the enhanced cooling. Some regions, such as the USA and the Mediterranean region, experience increased precipitation in accord with proxy paleoclimate evidence. The overall much drier climate over the ocean leads to higher sea surface salinity (SSS) in most ocean basins except for the North Atlantic where SSS is considerably lower due to an increase in the supply of fresh water from the Mississippi and Amazon rivers and presumably a decrease in salt transport by the weakened North Atlantic overturning circulation. The North Atlantic overturning stream function weakens to less than half of the control run value. The overturning is limited to a shallower depth (less than 1000 m) and its outflow is confined to the Northern Hemisphere. In the Southern Ocean, convection is much stronger than in the control run leading to a stronger overturning stream function associated with enhanced Antarctic Bottom Water formation. As a result, Southern Ocean water masses fill the entire deep ocean. The Antarctic Circumpolar Current (ACC) transport through the Drake Passage increases by about 25%. The ACC transport, despite weaker zonal winds, is enhanced due to changes in bottom pressure torque. The weakening of the overturning circulation in the North Atlantic and the accompanying 30% decrease in the poleward ocean heat transport contrasts with the strengthening of the overturning circulation in the Southern Ocean and a 40% increase in heat transport. As a result, sea ice coverage and thickness are affected in opposite senses in the two hemispheres. The LGM climate simulated by the coupled model is in reasonable agreement with paleoclimate proxy evidence. The dynamical response of the ocean in the coupled model plays an important role in determining the simulated, and undoubtedly, the actual, LGM climate.

[1]  S. Manabe,et al.  A global ocean-atmosphere climate model with seasonal variation for future studies of climate sensitivity , 1979 .

[2]  P. Gent,et al.  Isopycnal mixing in ocean circulation models , 1990 .

[3]  William R. Holland,et al.  Baroclinic and topographic influences on the transport in western boundary currents , 1973 .

[4]  E. Yu,et al.  Similar rates of modern and last-glacial ocean thermohaline circulation inferred from radiochemical data , 1996, Nature.

[5]  J. Beck,et al.  Abrupt changes in early Holocene tropical sea surface temperature derived from coral records , 1997, Nature.

[6]  G. Boer.  Climate change and the regulation of the surface moisture and energy budgets , 1993 .

[7]  E. Sikes,et al.  Equatorial Atlantic sea surface temperature for the last 30 kyr: a comparison of U37k', δ18O and foraminiferal assemblage temperature estimates , 1994 .

[8]  S. Philander,et al.  The climate of the Last Glacial Maximum: Results from a coupled atmosphere-ocean general circulation model , 1999 .

[9]  Richard G. Fairbanks,et al.  Tropical Temperature Variations Since 20,000 Years Ago: Modulating Interhemispheric Climate Change , 1994, Science.

[10]  R. Fairbanks.  A 17,000-year glacio-eustatic sea level record: influence of glacial melting rates on the Younger Dryas event and deep-ocean circulation , 1989, Nature.

[11]  A. Broccoli.  Tropical Cooling at the Last Glacial Maximum: An Atmosphere-Mixed Layer Ocean Model Simulation , 2000 .

[12]  William B. Curry,et al.  Weaker Gulf Stream in the Florida Straits during the Last Glacial Maximum , 1999, Nature.

[13]  G. Boer,et al.  The roles of radiation and dynamical processes in the El Niño-like response to global warming , 2002 .

[14]  Sandy P. Harrison,et al.  How well can we simulate past climates? Evaluating the models using global palaeoenvironmental datasets , 2000 .

[15]  N. Barkov,et al.  Aerosol concentrations over the last climatic cycle (160 kyr) from an Antarctic ice core , 1987, Nature.

[16]  Denis-Didier Rousseau,et al.  High-resolution record of the last Interglacial–glacial cycle in the Nussloch loess–palaeosol sequences, Upper Rhine Area, Germany , 2001 .

[17]  Laurent Labeyrie,et al.  Deepwater source variations during the last climatic cycle and their impact on the global deepwater circulation , 1988 .

[18]  Chris W. Hughes,et al.  Why Western Boundary Currents in Realistic Oceans are Inviscid: A Link between Form Stress and Bottom Pressure Torques , 2001 .

[19]  Michael Sarnthein,et al.  Global variations of surface ocean productivity in low and mid latitudes: Influence on CO2 reservoirs of the deep ocean and atmosphere during the last 21,000 years , 1988 .

[20]  R. Sausen,et al.  Periodically synchronously coupled integrations with the atmosphere-ocean general circulation model ECHAM3/LSG , 1998 .

[21]  Norman A. McFarlane,et al.  Greenhouse Gas–induced Climate Change Simulated with the CCC Second-Generation General Circulation Model , 1992 .

[22]  J. Kutzbach,et al.  The Influence of Changing Orbital Parameters and Surface Boundary Conditions on Climate Simulations for the Past 18 000 Years , 1986 .

[23]  G. Boer,et al.  Warming asymmetry in climate change simulations , 2001 .

[24]  William D. Hibler,et al.  Modeling Pack Ice as a Cavitating Fluid , 1992 .

[25]  T. Fichefet,et al.  A model study of the Atlantic thermohaline circulation during the last glacial maximum , 1994, Nature.

[26]  D. Robinson,et al.  Reappraisal of sea-ice distribution in Atlantic and Pacific sectors of the Southern Ocean at 18,000 yr BP , 1982, Nature.

[27]  S. Baum,et al.  Effect of vegetation on an ice-age climate model simulation , 1997 .

[28]  M. Horowitz,et al.  Glacial deep water geometry: South Atlantic benthic foraminiferal Cd/Ca and δ13C evidence , 2000 .

[29]  Anthony J. Broccoli,et al.  A coupled model study of the Last Glacial Maximum: Was part of the North Atlantic relatively warm? , 2001 .

[30]  D. Jolly,et al.  Tropical climates at the Last Glacial Maximum: a new synthesis of terrestrial palaeoclimate data. I. Vegetation, lake-levels and geochemistry , 1999 .

[31]  R. Sausen,et al.  Techniques for asynchronous and periodically synchronous coupling of atmosphere and ocean models , 1996 .

[32]  D. Schrag,et al.  Pore Fluid Constraints on the Temperature and Oxygen Isotopic Composition of the Glacial Ocean , 1996, Science.

[33]  W. Peltier,et al.  Ice Age Paleotopography , 1994, Science.

[34]  L. Harvey,et al.  Computational efficiency and accuracy of methods for asynchronously coupling atmosphere-ocean climate models. Part I. Testing with a mean annual model , 1986 .

[35]  D. Rind,et al.  Terrestrial Conditions at the Last Glacial Maximum and CLIMAP Sea-Surface Temperature Estimates: Are They Consistent? , 1985, Quaternary Research.

[36]  Edward A. Boyle,et al.  North Atlantic thermohaline circulation during the past 20,000 years linked to high-latitude surface temperature , 1987, Nature.

[37]  G. Boer,et al.  Climate sensitivity and response , 2003 .

[38]  R. Fairbanks,et al.  Variability in the deep and intermediate water circulation of the Atlantic Ocean during the past 25,000 years: Northern Hemisphere modulation of the Southern Ocean , 1987 .

[39]  P. Webster,et al.  Late Quaternary Ice Age Climates of Tropical Australasia: Interpretations and Reconstructions , 1978, Quaternary Research.

[40]  A. Mcintyre,et al.  Seasonal reconstructions of the earth's surface at the last glacial maximum , 1981 .

[41]  Alan C. Mix,et al.  Foraminiferal faunal estimates of paleotemperature: Circumventing the No‐analog problem yields cool Ice Age tropics , 1999 .

[42]  Y. Rosenthal,et al.  Cd/Ca changes in a Deep Cape Basin Core over the past 730,000 years: Response of circumpolar deepwater variability to northern hemisphere ice sheet melting? , 1994 .

[43]  A. Weaver,et al.  Simulated influence of carbon dioxide, orbital forcing and ice sheets on the climate of the Last Glacial Maximum , 1998, Nature.

[44]  C. Hillaire‐Marcel,et al.  Reconstruction of sea-surface temperature, salinity, and sea-ice cover in the northern North Atlantic during the last glacial maximum based on dinocyst assemblages , 2000 .

[45]  B. McAvaney,et al.  Influence of vegetation changes during the Last Glacial Maximum using the BMRC atmospheric general circulation model , 2001 .

[46]  T. Crowley,et al.  Late Pleistocene variations in Antarctic sea ice I: effect of orbital insolation changes , 1988 .

[47]  W. Peltier,et al.  Global water balance and atmospheric water vapour transport at last glacial maximum: climate simulations with the Canadian Climate Centre for Modelling and Analysis atmospheric general circulation model , 2000 .

[48]  S. Goldstein,et al.  Reduced North Atlantic Deep Water flux to the glacial Southern Ocean inferred from neodymium isotope ratios , 2000, Nature.

[49]  S. Manabe,et al.  Response of a Coupled Ocean–Atmosphere Model to Increasing Atmospheric Carbon Dioxide: Sensitivity to the Rate of Increase , 1999 .

[50]  S. Levitus,et al.  A diagnosis of interpentadal circulation changes in the North Atlantic , 1991 .

[51]  D. Rea,et al.  Quaternary Fluctuations in the Northern Hemisphere Trade Winds and Westerlies , 1985, Quaternary Research.

[52]  J. Hansen,et al.  Radiative forcing and climate response , 1997 .

[53]  M. Ledbetter.  Bottom-current speed in the Vema Channel recorded by particle size of sediment fine-fraction , 1984 .

[54]  Syukuro Manabe,et al.  Two Stable Equilibria of a Coupled Ocean-Atmosphere Model , 1988 .

[55]  Stefan Rahmstorf,et al.  Simulation of modern and glacial climates with a coupled global model of intermediate complexity , 1998, Nature.

[56]  G. Boer,et al.  A transient climate change simulation with greenhouse gas and aerosol forcing: experimental design and comparison with the instrumental record for the twentieth century , 2000 .

[57]  R. Sausen,et al.  Techniques for asynchronous and periodically synchronous coupling of atmosphere and ocean models , 1996 .

[58]  A. Molina-Cruz.  The Relation of the Southern Trade Winds to Upwelling Processes During the Last 75,000 years , 1977, Quaternary Research.

[59]  E. Boyle.  Cadmium and delta13C Paleochemical Ocean Distributions During the Stage 2 Glacial Maximum , 1992 .

[60]  A. Semtner.  A MODEL FOR THE THERMODYNAMIC GROWTH OF SEA ICE IN NUMERICAL INVESTIGATIONS OF CLIMATE , 1975 .

[61]  S. Manabe,et al.  Climate Calculations with a Combined Ocean-Atmosphere Model , 1969 .

[62]  Timothy T. Barrows,et al.  Sea-surface temperatures of the southwest Pacific Ocean during the Last Glacial Maximum , 2000 .

[63]  G. P. Lohmann,et al.  Carbon Isotopic Changes in Benthic Foraminifera from the Western South Atlantic: Reconstruction of Glacial Abyssal Circulation Patterns , 1982, Quaternary Research.

[64]  J. Duplessy,et al.  Deep water formation in the North Atlantic Ocean during the last ice age , 1980, Nature.

[65]  Seong-Joong Kim,et al.  A coupled climate model simulation of the Last Glacial Maximum, Part 1: transient multi-decadal response , 2002 .

[66]  J. Duplessy,et al.  Changes in the distribution of δ13C of deep water ΣCO2 between the Last Glaciation and the Holocene , 1988 .

[67]  M. Sarnthein,et al.  Glacial and interglacial wind regimes over the eastern subtropical Atlantic and North-West Africa , 1981, Nature.

[68]  Syukuro Manabe,et al.  Transient responses of a coupled ocean-atmosphere model to gradual changes of atmospheric CO2 , 1991 .

[69]  R. Pierrehumbert,et al.  Huascaran δ18O as an indicator of tropical climate during the Last Glacial Maximum , 1999 .

[70]  A. Weaver,et al.  The Canadian Centre for Climate Modelling and Analysis global coupled model and its climate , 2000 .

[71]  S. Manabe,et al.  The Influence of Continental Ice Sheets on the Climate of an Ice Age , 1985 .

[72]  W. Broecker,et al.  Cooling of Tropical Brazil (5�C) During the Last Glacial Maximum , 1995, Science.

[73]  E. Boyle,et al.  Deep Circulation of the North Atlantic over the Last 200,000 Years: Geochemical Evidence , 1982, Science.

[74]  N. Shackleton,et al.  Intensified deep Pacific inflow and ventilation in Pleistocene glacial times , 2001, Nature.

[75]  C. Lorius,et al.  Vostok ice core provides 160,000-year record of atmospheric CO2 , 1987, Nature.

[76]  A. Kitoh,et al.  A simulation of the Last Glacial Maximum with a coupled atmosphere‐ocean GCM , 2001 .

[77]  E. Bard,et al.  Interhemispheric synchrony of the last deglaciation inferred from alkenone palaeothermometry , 1997, Nature.

[78]  W. Broecker.  Oxygen Isotope Constraints on Surface Ocean Temperatures , 1986, Quaternary Research.

[79]  C. Weyhenmeyer,et al.  Cool glacial temperatures and changes in moisture source recorded in oman groundwaters , 2000, Science.

[80]  P. Colinvaux,et al.  A Long Pollen Record from Lowland Amazonia: Forest and Cooling in Glacial Times , 1996, Science.

[81]  J. Blanchet,et al.  The Canadian Climate Centre Second-Generation General Circulation Model and Its Equilibrium Climate , 1992 .

[82]  A. Royer,et al.  Ice age aerosol content from East Antarctic ice core samples and past wind strength , 1981, Nature.

[83]  H. Elderfield,et al.  Past temperature and δ18O of surface ocean waters inferred from foraminiferal Mg/Ca ratios , 2000, Nature.

[84]  Laurent Labeyrie,et al.  Changes in east Atlantic deepwater circulation over the last 30 , 1994 .

[85]  S. J. Lambert,et al.  Second-order space-time climate difference statistics , 2001 .

[86]  E. Mosley‐Thompson,et al.  Late Glacial Stage and Holocene Tropical Ice Core Records from Huascar�n, Peru , 1995, Science.

[87]  U. Beyerle,et al.  Palaeotemperature reconstruction from noble gases in ground water taking into account equilibration with entrapped air , 2000, Nature.