The global response to Younger Dryas boundary conditions in an AGCM simulation

 Geological evidence points to a global Younger Dryas (YD) climatic oscillation during the last glacial/ present interglacial transition phase. A convincing mechanism to explain this global YD climatic oscillation is not yet available. Nevertheless, a profound understanding of the mechanism behind the YD climate would lead to a better understanding of climate variability. Therefore, the Hamburg atmospheric circulation model was used to perform four numerical experiments on the YD climate. The objective of this study is to improve the understanding of different forcings influencing climate during the last glacial/interglacial transition and to investigate to what extent the model response agrees with global geological evidence of YD climate change. The following boundary conditions were altered: sea surface conditions, ice sheets, insolation and atmospheric CO₂ concentration. Sea surface temperatures based on foraminiferal assemblages proved to produce insufficient winter cooling in the N Atlantic Ocean in two experiments. It is proposed that this discrepancy is caused by uncertainties in the reconstruction method of sea surface temperatures. Therefore, a model-derived set of Atlantic surface ocean conditions was prescribed in a subsequent simulation. However, the latter set represented an Atlantic Ocean without a thermohaline circulation, which is not in agreement with evidence from ocean cores. The global response to the boundary conditions was analysed using three variables, namely surface temperature, zonal wind speed and precipitation. The statistical significance of the changes was tested with a two-tailed t-test. Moreover, the significant responses to cooled oceans were compared with geological evidence of a YD oscillation. This comparison revealed a good match in Europe, Greenland, Atlantic Canada and the N Pacific region, explaining the YD oscillation in these regions as a response to cooled N Atlantic and N Pacific Oceans. However, the results leave the YD climate in other regions completely unexplained. This reflects either an insufficient set of boundary conditions or the important role played by feedbacks within the coupled atmosphere-ocean-ice system. These feedbacks are poorly represented in the used atmospheric model, since ice sheets and the ocean surface conditions have to be prescribed. Received: 30 July 1996 / Accepted: 12 February 1997     

The impact of cold North Atlantic sea surface temperatures on climate: implications for the Younger Dryas cooling (11–10 k)

The sensitivity of global climate to colder North Atlantic sea surface temperatures is in vestigated with the use of the GISS general circulation model. North Atlantic ocean temperatures 18,000 B.P., resembling those prevalent during the Younger Dryas, were incorporated into the model of the present climate and also into an experiment using orbital parameters and land ice characteristic of 11,000 B.P. The results show that with both 11,000 B.P. and present conditions the colder ocean temperatures produce cooling over western and central Europe, in good agreement with Younger Dryas paleoclimatic evidence. Cooling also occurs over extreme eastern North America, although the precise magnitude and location depends upon the specification of ocean temperature change in the western Atlantic. Despite the presence of increased land ice and colder ocean temperatures, the Younger Dryas summer air temperatures at Northern Hemisphere midlatitudes in the model are warmer than those of today due to changes in the orbital parameters, chiefly precession, and atmospheric subsidence at the perimeter of the ice sheets.     

Coupled climate modelling of ocean circulation changes during ice age inception

Freshening of high latitude surface waters can change the large-scale oceanic transport of heat and salt. Consequently, atmospheric and sea ice perturbations over the deep water production sites excite a large-scale response establishing an oceanic "teleconnection" with time scales of years to centuries. To study these feedbacks, a coupled atmosphere-ocean-sea ice model consisting of a two dimensional atmospheric energy and moisture balance model (EMBM) coupled to a thermodynamic sea ice model and an ocean general circulation model is utilised. The coupled model reproduces many aspects of the present oceanic circulation. We also investigate the climate impact of changes in fresh water balance during an ice age initiation. In this experiment part of the precipitation over continents is stored within continental ice sheets. During the buildup of ice sheets the oceanic stratification in the North Atlantic is weakened by a reduced continental run-off leading to an enhanced thermohaline circulation. Under these conditions salinity is redistributed such that deep water is more saline than under present conditions. Once the ice sheets built up, we simulate an ice age climate without net fresh water storage on the continents. In this case the coupled model reproduces the shallow and weak overturning cell, an ice edge advance insulating the upper ocean, and many other aspects of the glacial circulation.     

Nonlinear variability of the climatic system from singular and power spectra of Late Quaternary records

Stable-isotope records from seven marine cores and one ice core provide invaluable information on the intricate behavior of the climatic system over time scales of 104 to 105 years. These records, in conjunction with a simple coupled climate model, help us understand major mechanisms of paleoclimatic variability. The time intervals covered by the records include the last glacial-interglacial cycle. In spite of the difference in the nature of the records, common features are revealed by advanced spectral-analysis tools. The dominant features are the presence of orbital frequencies, on the one hand, and a low number of internal degrees of freedom, on the other. The climatic system appears therefore to act on the Quaternary time scales considered as a forced nonlinear oscillator. The internal mechanisms giving rise to the aperiodic oscillations include ice-albedo feedback, precipitation-temperature feedback, and interactions between the ice sheets and the bedrock.     

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.     

Potential climatic transitions with profound impact on Europe

We discuss potential transitions of six climatic subsystems with large-scale impact on Europe, sometimes denoted as tipping elements. These are the ice sheets on Greenland and West Antarctica, the Atlantic thermohaline circulation, Arctic sea ice, Alpine glaciers and northern hemisphere stratospheric ozone. Each system is represented by co-authors actively publishing in the corresponding field. For each subsystem we summarize the mechanism of a potential transition in a warmer climate along with its impact on Europe and assess the likelihood for such a transition based on published scientific literature. As a summary, the ‘tipping’ potential for each system is provided as a function of global mean temperature increase which required some subjective interpretation of scientific facts by the authors and should be considered as a snapshot of our current understanding.     

Simulation of climate phase lags in response to precession and obliquity forcing and the role of vegetation

Long (130,000 years) transient simulations with a coupled model of intermediate complexity (CLIMBER-2) have been performed. The main objective of this study is to examine leads and lags in the response to the climate system to separate obliquity and precession-induced insolation changes. Focus is on the role of internal feedbacks in the coupled atmosphere/ocean/sea-ice/vegetation system. No interactive ice sheets were used. The results show that leads and lags occur in response to the African/Asian monsoon, temperatures at high latitudes and the Atlantic thermohaline circulation. For the monsoon, leads and lags of the monthly precipitation with respect to the precession parameter were found, which are strongly modified by vegetation. In contrast, no lag was observed for the annual precipitation. At high latitudes during late winter/early spring a vegetation-induced lag with respect to the precession parameter was found in surface air temperatures. Again, no annual lag was detected. The lag in the monthly surface air temperatures induces a lag in the annual overturning in the Atlantic Ocean by changing the strength of the deep convection. The lag is several thousand years. The obliquity-related forcing does not give rise to lags in the climate system. We conclude that lags in monthly climatic variables, which are due to vegetation feedbacks, can result in an annual lag when a climatic process (like deep water formation) acts as a filter for certain months.     

Response in atmospheric circulation and sources of Greenland precipitation to glacial boundary conditions

The response in northern hemisphere atmospheric circulation and the resulting changes in moisture sources for Greenland precipitation to glacial boundary conditions are studied in NCAR’s CCM3 atmospheric general circulation model fitted with a moisture tracking functionality. We employ both the CLIMAP SST reconstruction and a modification thereto with reconstructions of glacial ice sheets and land masks. The individual components of the boundary conditions are added first one at a time and, finally, together. These steps show the atmospheric circulation to respond approximately linearly to the boundary condition changes, and the full glacial change may thus be decomposed into contributions from SST and topography changes, respectively. We find that using the CLIMAP SST reconstruction leads to a shift from Atlantic toward Pacific source regions not found with the modified reconstruction having cooler tropics and less sea ice. The occurrence of such a shift depends chiefly on the SST reconstruction and not on the existence of the large northern hemisphere glacial ice sheets. The influence of these circulation changes on important factors for ice core interpretation such as precipitation seasonality, condensation temperatures and source temperatures are assessed.     

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.     

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     

Multiple instabilities and modes of glacial rhythmicity in the plio-Pleistocene: a general theory of late Cenozoic climatic change

It has been noted that several distinct modes of glacial oscillation have existed during the past few million years, ranging from low-amplitude, high-frequency oscillations in the early Pliocene, through relatively high amplitude, predominantly near-40 ky period, oscillations in the late Pliocene and early Pleistocene, to the major near-100 ky period oscillations of the late Pleistocene. In addition to other plausible mechanisms suggested previously to explain aspects of this multirhythmic phenomenon, we now illustrate another possible contributor to this type of behavior based on the hypothesis that the slow-response climatic system is bistable and that two kinds of internal instability may be operative along with externally imposed forcing due to earth-orbital (Milankovitch) radiation changes and slow, tectonically-induced changes in atmospheric carbon dioxide. These two instabilities have been discussed previously: one is due to positive feedback in the global carbon cycle leading to near-100 ky free oscillations of the ice sheets, and the other is due to the potential for ice-calving catastrophes associated with bedrock variations that can lead to oscillations of a period near 40 ky, independent of obliquity forcing. Within the framework of a dynamical model containing the possibility for these two instabilities, as well as for stable modes, we show (1) how Milankovitch radiative changes or stochastic forcing influencing ice sheets can induce aperiodic (chaotic) transitions between the possible stable and unstable modes, and more significantly, (2) how progressive, long-term, tectonically-induced, changes in carbon dioxide, acting in concert with earth-orbital radiative variations in high Northern Hemisphere latitudes, can force systematic transitions between the modes. Such systematic changes can result in an ice mass chronology for the past 5 My that is qualitatively similar to the observed record of global ice mass. In essence, we have constructed a minimum dynamical model of the late Cenozoic climatic changes, containing what are believed to be the main physical factors determining these changes: ice mass, bedrock depression, atmospheric carbon dioxide concentration, deep ocean thermohaline state, Milankovitch radiation forcing, and slow tectonically-induced carbon dioxide forcing. This model forms the basis for a coherent theory for the complex climatic events of this long period.     

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.     

Effects of variability of sea ice transport through the Fram Strait on the intensity of the Atlantic deep circulation

Results from an ice-ocean coupled model are used to investigate the impact of long-term variability in sea ice transport at the Fram Strait on the intensity of the Atlantic deep circulation. An increase (or decrease) in sea ice transport through the Fram Strait leads to a stronger (or weaker) deep circulation in the Atlantic. Change in the sea ice transport is accompanied by a salinity anomaly in the surface layer of the Arctic Ocean. Such an anomaly could inversely affect the Atlantic circulation once it reaches deep water formation regions. If the Canadian Archipelago is closed, the anomaly is subsequently transported through the Fram Strait, and counters the initial changes in the Atlantic deep circulation. On the other hand, if the Canadian Archipelago is open, some of the anomaly is transported to the Canadian Archipelago, and the initial change in the Atlantic deep circulation persists. In the Arctic Ocean basin, the time scale and path of the salinity anomalys propagation depends on the large-scale flow at the surface of the Arctic Ocean. Our results suggest that the salinity anomaly transport and its propagation pathway out of the Arctic Ocean are important determinants of the role of sea ice transport variability through the Fram Strait in controlling the intensity of the Atlantic deep circulation.     

Simulation of the Atlantic meridional overturning circulation in an atmosphere-ocean global coupled model. Part II: weakening in a climate change experiment: a feedback mechanism

Most state-of-the art global coupled models simulate a weakening of the Atlantic meridional overturning circulation (MOC) in climate change scenarios but the mechanisms leading to this weakening are still being debated. The third version of the CNRM (Centre National de Recherches Météorologiques) global atmosphere-ocean-sea ice coupled model (CNRM-CM3) was used to conduct climate change experiments for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4). The analysis of the A1B scenario experiment shows that global warming leads to a slowdown of North Atlantic deep ocean convection and thermohaline circulation south of Iceland. This slowdown is triggered by a freshening of the Arctic Ocean and an increase in freshwater outflow through Fram Strait. Sea ice melting in the Barents Sea induces a local amplification of the surface warming, which enhances the cyclonic atmospheric circulation around Spitzberg. This anti-clockwise circulation forces an increase in Fram Strait outflow and a simultaneous increase in ocean transport of warm waters toward the Barents Sea, favouring further sea ice melting and surface warming in the Barents Sea. Additionally, the retreat of sea ice allows more deep water formation north of Iceland and the thermohaline circulation strengthens there. The transport of warm and saline waters toward the Barents Sea is further enhanced, which constitutes a second positive feedback.     

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.     

On the freshwater forcing and transport of the Atlantic thermohaline circulation

The 'conveyor belt' circulation of the Atlantic Ocean transports large amounts of heat northward, acting as a heating system for the northern North Atlantic region. It is widely thought that this circulation is driven by atmospheric freshwater export from the Atlantic catchment region, and that it transports freshwater northward to balance the loss to the atmosphere. Using results from a simple conceptual model and a global circulation model, it is argued here that the freshwater loss to the atmosphere arises mainly in the subtropical South Atlantic and is balanced by northward freshwater transport in the wind-driven subtropical gyre, while the thermohaline circulation transports freshwater southward. It is further argued that the direction of freshwater transport is closely linked to the dynamical regime and stability of the 'conveyor belt': if its freshwater transport is indeed southward, then its flow is purely thermally driven and inhibited by the freshwater forcing. In this case the circulation is not far from Stommel's saddle-node bifurcation, and a circulation state without NADW formation would also be stable. Received: 10 February 1996 / Accepted: 30 May 1996     

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.     

Climatic fluctuations on the century time scale: A review of high-resolution proxy data and possible mechanisms

We review the century time scale climatic variability that is observed in high-resolution proxy data records covering the past 10 000 yr. Cyclic variations with time scales ranging from 50 to 400 yr occur in oxygen isotope ratios derived from ice cores, tree-ring index series, pollen records and sea-ice extents. Century time scale cycles can also be identified in some biological and historical records and in long-term instrumental observations. In order to appreciate the century scale cycles in the context of climatic variability in general, a brief survey of all climatic time scales is presented.The traditional interpretation that decadal-to-century scale fluctuations in the climate system are externally forced, e.g. by variations in solar properties, is questioned. A different mechanism for these fluctuations is proposed on the basis of recent findings of numerical models of the ocean's thermohaline circulation. The results indicate that this oceanic circulation exhibits natural variability on the century time scale which produces oscillations in the ocean-to-atmosphere heat flux. Although global in extent, these fluctuations are largest in the Atlantic Ocean.     

Evaluation of ocean and climate models using present‐day observations and forcing

Abstract

The most common method used to evaluate climate models involves spinning them up under perpetual present‐day forcing and comparing the model results with present‐day observations. This approach clearly ignores any potential long‐term memory of the model ocean to past climatic conditions. Here we examine the validity of this approach through the 6000‐year integration of a coupled atmosphere–ocean–sea‐ice model. The coupled model is initially spun‐up with atmospheric CO₂ concentrations and orbital parameters applicable for 6KBP. The model is then integrated forward in time to 2100. Results from this transient coupled model simulation are compared with the results from two additional simulations, in which the model is spun up with perpetual 1850 (preindustrial) and 1998 (present‐day) atmospheric CO₂ concentrations and orbital parameters. This comparison leads to substantial differences between the equilibrium climatologies and the transient simulation, even at 1850 (in weakly ventilated regions), prior to any significant changes in atmospheric CO₂. When compared to the present‐day equilibrium climatology, differences are very large: the global mean surface air and sea surface temperatures are ,0.5°C and ,0.4°C colder, respectively, deep ocean temperatures are substantially cooler, Southern Hemisphere sea‐ice cover is 38% larger, and the North Atlantic conveyor 16% weaker in the transient case. These differences are due to the long timescale memory of the deep ocean to climatic conditions which prevailed throughout the late Holocene, as well as to its large thermal inertia. It is also demonstrated that a ‘cold start’ global warming simulation (one that starts from a 1998 equilibrium climatology) underestimates the global temperature increase at 2100 by ,10%. Our results question the accuracy of current techniques for climate model evaluation and underline the importance of using paleoclimatic simulations in parallel with present‐day simulations in this evaluation process.     

Atmospheric water vapor flux,bifurcation of the thermohaline circulation,and climate change

Latitudinal heat transport in the ocean and atmosphere represents a fundamental process of the Earth's climate system. The ocean component of heat transport is effected by the thermohaline circulation. Changes in this circulation, and hence latitudinal heat transport, would have a significant effect on global climate. Paleoclimate evidence from the Greenland ice cores and deep sea sediment cores suggests that during much of glacial time the climate system oscillated between two different states. Bimodal equilibrium states of the thermohaline circulation have been demonstrated in climate models. We address the question of the role of the atmospheric hydrological cycle on the global thermohaline circulation and the feedback to the climate system through changes in the ocean's latitudinal heat transport, with a simple coupled ocean-atmosphere energy-salt balance model. Two components of the atmospheric hydrological cycle, i.e., latitudinal water vapor transport and the net flux of water vapor from the Atlantic to the Pacific Ocean appear to play separate roles. If the inter-basin transport is sufficiently large, small changes in water vapor transport over the North Atlantic can effect bifurcation or a rapid transition between two different equilibria in the global thermohaline circulation; maximum difference between the modes occurs in the North Atlantic. If the inter-basin transport is from the Pacific to the Atlantic and sufficiently large, latitudinal vapor transport in the North Pacific controls the bifurcations, with maximum changes occurring in the North Pacific. For intermediate values of inter-basin transport, no rapid transitions occur in either basin. In the regime with vapor flux from the Atlantic to the Pacific, the on mode has strong production of deep water in the North Atlantic and a large flux of heat to the atmosphere from the high latitude North Atlantic. The off mode has strong deep water production in the Southern Ocean and weak production in the North Pacific. Heat transport into the high latitude North Atlantic by the ocean is reduced to about 20% of the on mode value. For estimated values of water vapor transport for the present climate the model asserts that while water vapor transport from the Atlantic to the Pacific Ocean is sufficiently large to make the North Atlantic the dominant region for deep water production, latitudinal water vapor transport is sufficiently low that the thermohaline circulation appears stable, i.e., far from a bifurcation point. This conclusion is supported to some extent by the fact that the high latitude temperature of the atmosphere as recorded in the Greenland ice cores has changed little over the last 9000 years.