Effect of seasonal forcing on global circulation in a world ocean general circulation model

 Effects of the seasonal variation in thermohaline and wind forcing on the abyssal circulation are investigated by using an ocean general circulation model. To isolate effects of the seasonality in the thermohaline forcing from those in the wind forcing, we carry out three experiments with (1) annual-mean wind forcing and perpetual-winter thermohaline forcing, (2) annual-mean wind forcing and seasonal thermohaline forcing, and (3) seasonal wind forcing and seasonal thermohaline forcing. The deep water under the seasonal thermohaline forcing becomes warmer than under the perpetual-winter thermohaline forcing. Although the perpetual-winter thermohaline forcing is widely used and believed to reproduce the deep water better than the annual-mean forcing, the difference between the results of the perpetual-winter and the seasonal thermohaline forcing is significant. The seasonal variation of the Ekman convergence and divergence produces meridional overturning cells extending to the bottom because the period of seasonal cycle is shorter than the adjustment timescale by baroclinic Rossby waves. The heat transport owing to those Ekman flows and temperature anomalies makes the upper water (0–200 m) colder at low to mid-latitudes (40S–40N) and warmer at high latitudes. Also the deep water becomes warmer owing to the warming of the northern North Atlantic, the main source region of North Atlantic Deep Water. The model is also synchronously (i.e., without acceleration) integrated with seasonal forcing for 5400 y. A past study suggested that under seasonal forcing, a sufficient equilibrium state can be achieved after only decades of synchronous integration following more than 10 000 y of accelerated integration. Here, the result so obtained is compared with that of the 5400-y synchronous integration. The difference in the global average temperature is as small as 0.12 °C, and most of the difference is confined to the Southern Ocean. Received: 1 May 1998 / Accepted: 5 January 1999     

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.     

Equilibrium response of thermohaline circulation to large changes in atmospheric CO2 concentration

This study evaluates the equilibrium response of a coupled ocean–atmosphere model to the doubling, quadrupling, and halving of CO₂ concentration in the atmosphere. Special emphasis in the study is placed upon the response of the thermohaline circulation in the Atlantic Ocean to the changes in CO₂ concentration of the atmosphere. The simulated intensity of the thermohaline circulation (THC) is similar among three quasi-equilibrium states with the standard, double the standard, and quadruple the standard amounts of CO₂ concentration in the atmosphere. When the model atmosphere has half the standard concentration of CO₂, however, the THC is very weak and shallow in the Atlantic Ocean. Below a depth of 3 km, the model oceans maintain very thick layer of cold bottom water with temperature close to –2 °C, preventing the deeper penetration of the THC in the Atlantic Ocean. In the Circumpolar Ocean of the Southern Hemisphere, sea ice extends beyond the Antarctic Polar front, almost entirely covering the regions of deepwater ventilation. In addition to the active mode of the THC, there exists another stable mode of the THC for the standard, possibly double the standard (not yet confirmed), and quadruple the standard concentration of atmospheric carbon dioxide. This second mode is characterized by the weak, reverse overturning circulation over the entire Atlantic basin, and has no ventilation of the entire subsurface water in the North Atlantic Ocean. At one half the standard CO₂ concentration, however, the intensity of the first mode is so weak that it is not certain whether there are two distinct stable modes or not. The paleoceanographic implications of the results obtained here are discussed as they relate to the signatures of the Cenozoic changes in the oceans.An erratum to this article can be found at     

On the horizontal density ratio in the upper ocean

The horizontal density ratio in the upper ocean is examined using SeaSoar data collected over the last 15 years in the Pacific, Atlantic, and Indian Oceans. The horizontal density ratio R is defined to be the ratio of the relative effect of temperature and salinity on density. A front with a horizontal density ratio of 1 is said to be compensated since temperature and salinity gradients compensate in their effect on density. The statistics of density ratio are examined through calculation of conditional probability density functions. Case studies from each of the oceans elucidate processes affecting the density ratio. Global distributions of density ratio are calculated as functions of mixed-layer depth, distance below the mixed layer, and magnitude of thermohaline variability. Compensation is found in all oceans, on 3–4 km horizontal scales, when the mixed layer is deep and significant thermohaline variability exists. The tendency for compensation is stronger as mixed-layer depth increases. Conversely, compensation is not typical in shallow mixed layers, or when thermohaline variability is weak. The thermocline density ratio is found to be 2, in agreement with previous observational studies, and consistent with the process of salt fingering. The transition from R=1 in the mixed layer to R=2 in the thermocline is sharp when the mixed layer is deep. The ubiquity of compensation in the mixed layer is consistent with recent theory that suggests horizontal eddy diffusivity is a growing function of density gradient.     

Glacial Thermohaline Circulation and Climate：Forcing from the North or South？

Based on the evidence available from both observations and model simulations, the author proposes a view that may provide a unified interpretation of the North Atlantic thermohaline variability. Because of the slow response time of the Southern Ocean （millennia） and the relatively faster response time of the North Atlantic （centuries）, the North Atlantic thermohaline circulation is controlled predominantly by the climate forcing over the Southern Ocean at the long glacial cycle timescales, but by the North Atlantic climate forcing at the short millennial timescaies.     

Moisture transport from the Atlantic to the Pacific basin and its response to North Atlantic cooling and global warming

Atmospheric moisture transport from the Atlantic to the Pacific basin plays an important role in regulating North Atlantic salinity and thus the strength of the thermohaline circulation. Potential changes in the strength of this moisture transport are investigated for two different climate-change scenarios: North Atlantic cooling representative of Heinrich events, and increased greenhouse gas (GHG) forcing. The effect of North Atlantic cooling is studied using a coupled regional model with comparatively high resolution that successfully simulates Central American gap winds and other important aspects of the region. Cooler North Atlantic sea surface temperature (SST) in this model leads to a regional decrease of atmospheric moisture but also to an increase in wind speed across Central America via an anomalous pressure gradient. The latter effect dominates, resulting in a 0.13 Sv (1 Sv = 10⁶ m³ s^?1) increase in overall moisture transport to the Pacific basin. In fresh water forcing simulations with four different general circulation models, the wind speed effect is also present but not strong enough to completely offset the effect of moisture decrease except in one model. The influence of GHG forcing is studied using simulations from the Intergovernmental Panel on Climate Change archive. In these simulations atmospheric moisture increases globally, resulting in an increase of moisture transport by 0.25 Sv from the Atlantic to Pacific. Thus, in both scenarios, moisture transport changes act to stabilize the thermohaline circulation. The notion that the Andes effectively block moisture transport from the Atlantic to the Pacific basin is not supported by the simulations and atmospheric reanalyses examined here. This indicates that such a blocking effect does not exist or else that higher resolution is needed to adequately represent the steep orography of the Andes.     

A zonally averaged,three-basin ocean circulation model for climate studies

A two-dimensional, three-basin ocean model suitable for long-term climate studies is developed. The model is based on the zonally averaged form of the primitive equations written in spherical coordinates. The east-west density difference which arises upon averaging the momentum equations is taken to be proportional to the meridional density gradient. Lateral exchanges of heat and salt between the basins are explicitly resolved. Moreover, the model includes bottom topography and has representations of the Arctic Ocean and of the Weddell and Ross seas. Under realistic restoring boundary conditions, the model reproduces the global conveyor belt: deep water is formed in the Atlantic between 60 and 70°N at a rate of about 17 Sv (1 Sv=10⁶ m3 s^–1) and in the vicinity of the Antarctic continent, while the Indian and Pacific basins show broad upwelling. Superimposed on this thermohaline circulation are vigorous wind-driven cells in the upper thermocline. The simulated temperature and salinity fields and the computed meridional heat transport compare reasonably well with the observational estimates. When mixed boundary conditions (i.e., a restoring condition on sea-surface temperature and flux condition on sea-surface salinity) are applied, the model exhibits an irregular behavior before reaching a steady state characterized by self-sustained oscillations of 8.5-y period. The conveyor-belt circulation always results at this stage. A series of perturbation experiments illustrates the ability of the model to reproduce different steady-state circulations under mixed boundary conditions. Finally, the model sensitivity to various factors is examined. This sensitivity study reveals that the bottom topography and the presence of a submarine meridional ridge in the zone of the Drake Passage play a crucial role in determining the properties of the model bottom-water masses. The importance of the seasonality of the surface forcing is also stressed.     

Sensitivity study of the seasonal mean circulation in the northern South China Sea

A study of the circulation in the northern South China Sea （SCS） is carried out with the aid of a three-dimensional, high-resolution regional ocean model. One control and two sensitivity experiments are performed to qualitatively investigate the effects of surface wind forcing, Kuroshio intrusion, and bottom topographic influence on the circulation in the northern SCS. The model results show that a branch of the Kuroshio in the upper layer can intrude into the SCS and have direct influence on the circulation over the continental shelf break in the northern SCS. There are strong southward pressure gradients along a zonal belt largely seaward of the continental slope. The pressure gradients are opposite in the southern and northern parts of the Luzon Strait, indicating inflow and outflow through the strait, respectively. The sensitivity experiments suggest that the Kuroshio intrusion is responsible for generating the imposed pressure head along the shelf break and has no obvious seasonal variations. The lateral forcing through the Luzon Strait and Taiwan Strait can induce the southwestward slope current and the northeastward SCS Warm Current in the northern SCS. Without the lateral forcing, there is the continental slope. The wind forcing mainly causes the The wind-induced water pile-up results in the southward no high-pressure-gradient zonal belt seaward of seasonal variation of the circulation in the SCS. high pressure gradient along the northwestern boundary of the basin. Without the blocking of the plateau around Dongsha Islands, the intruded Kuroshio tends to extend northwest and the SCS branch of the Kuroshio becomes wider and stronger. The analyses presented here are qualitative in nature but should lead to a better understanding of the oceanic responses in the northern SCS to these external influence factors.     

Formation and pathways of North Atlantic Deep Water in a coupled ice–ocean model of the Arctic–North Atlantic Oceans

We investigate the formation process and pathways of deep water masses in a coupled ice–ocean model of the Arctic and North Atlantic Oceans. The intent is to determine the relative roles of these water masses from the different source regions (Arctic Ocean, Nordic Seas, and Subpolar Atlantic) in the meridional overturning circulation. The model exhibits significant decadal variability in the deep western boundary current and the overturning circulation. We use detailed diagnostics to understand the process of water mass formation in the model and the resulting effects on the North Atlantic overturning circulation. Particular emphasis is given to the multiple sources of North Atlantic Deep Water, the dominant deep water masses of the world ocean. The correct balance of Labrador Sea, Greenland Sea and Norwegian Sea sources is difficult to achieve in climate models, owing to small-scale sinking and convection processes. The global overturning circulation is described as a function of potential temperature and salinity, which more clearly signifies dynamical processes and clarifies resolution problems inherent to the high latitude oceans. We find that fluxes of deep water masses through various passages in the model are higher than observed estimates. Despite the excessive volume flux, the Nordic Seas overflow waters are diluted by strong mixing and enter the Labrador Sea at a lighter density. Through strong subpolar convection, these waters along with other North Atlantic water masses are converted into the densest waters [similar density to Antarctic Bottom Water (AABW)] in the North Atlantic. We describe the diminished role of salinity in the Labrador Sea, where a shortage of buoyant surface water (or excess of high salinity water) leads to overly strong convection. The result is that the Atlantic overturning circulation in the model is very sensitive to the surface heat flux in the Labrador Sea and hence is correlated with the North Atlantic Oscillation. As strong subpolar convection is found in other models, we discuss broader implications.     

The stability of the North Atlantic thermohaline circulation in a coupled ocean-atmosphere general circulation model

 The stability of the Atlantic thermohaline circulation against meltwater input is investigated in a coupled ocean-atmosphere general circulation model. The meltwater input to the Labrador Sea is increased linearly for 250 years to a maximum input of 0.625 Sv and then reduced again to 0 (both instantaneously and linearly decreasing over 250 years). The resulting freshening forces a shutdown of the formation of North Atlantic deepwater and a subsequent reversal of the thermohaline circulation of the Atlantic, filling the deep Atlantic with Antarctic bottom water. The change in the overturning pattern causes a drastic reduction of the Atlantic northward heat transport, resulting in a strong cooling with maximum amplitude over the northern North Atlantic and a southward shift of the sea-ice margin in the Atlantic. Due to the increased meridional temperature gradient, the intertropical convergence zone over the Atlantic is displaced southward and the westerlies in the Northern Hemisphere gain strength. We identify four main feedbacks affecting the stability of the thermohaline circulation: the change in the overturning circulation of the Atlantic leads to longer residence times of the surface water in high-northern latitudes, which allows them to accumulate more precipitation and runoff from the continents. As a consequence the stratification in the North Atlantic becomes more stable. This effect is further amplified by an enhanced northward atmospheric water vapour transport, which increases the freshwater input into the North Atlantic. The reduced northward oceanic heat transport leads to colder sea-surface temperatures and an intensification of the atmospheric cyclonic circulation over the Norwegian Sea. The associated Ekman transports cause increased upwelling and increased freshwater export with the East Greenland Current. Both the cooling and the wind-driven circulation changes largely compensate for the effects of the first two feedbacks. The wind-stress feedback destabilizes modes without deep water formation in the North Atlantic, but has been neglected in almost all studies so far. After the meltwater input stops, the North Atlantic deepwater formation resumed in all experiments and the meridional overturning returned within 200 years to a conveyor belt pattern. This happened although the formation of North Atlantic deep water was suppressed in one experiment for more than 300 years and the Atlantic overturning had settled into a circulation pattern with Antarctic bottom water as the only source of deep water. It is a clear indication that cooling and wind-stress feedback are more effective, at least in our model, than advection feedback and increased atmospheric water vapour transport. We conclude that the conveyor belt-type thermohaline circulation seems to be much more stable than hitherto assumed from experiments with simpler models. Received 31 January 1996/Accepted 22 August 1996     

Centennial oscillations in an ocean-ice coupled model

Irregular centennial oscillations, with a spectral peak at 106 years, were obtained from an ocean-ice coupled model for the North Atlantic with realistic coastline and bottom topography. The model’s thermohaline circulation is forced by mixed boundary conditions, i.e., a Haney-type relaxation condition for temperature, but an equivalent virtual salt flux condition for salinity. All forcing fields are taken from the observed monthly mean climatological wind stress and buoyancy fluxes. The oscillations appeared in the form of a surface-intensified tripole in both the sea surface temperature and salinity fields located in the vicinity of the Labrador Sea. The oscillations involve a delicate interplay between heat and fresh water advection by meridional overturning circulation, horizontal gyres, vertical convection, and the seasonal cycle. The oscillations are primarily control?led by the salinity component of the circulation; however, sea ice plays a minor role in driving the oscillations observed in the model. On the other hand, a regular seasonal cycle in the forcing fields is an important ingredient for the centennial oscillations.     

Reconciling theories of a mechanically driven meridional overturning circulation with thermohaline forcing and multiple equilibria

It has recently been suggested that the structure and strength of the meridional overturning circulation in the global ocean is governed by the input of mechanical energy to the system by winds and tides. However, it is not clear how this suggestion relates to the existence of multiple equilibria of the meridional overturning circulation, which depends on thermohaline feedbacks and is more consistent with a buoyancy-driven view of the circulation. Both theories have been illustrated by box models in the past (Stommel in Tellus 13:224–230, 1961; Gnanadesikan in Science 283:2077–2079, 1999). Here we incorporate these two theories into a single box model in an attempt to reconcile the roles of mechanical and buoyancy forcing in driving the meridional overturning circulation. The box model has two equilibrium solutions, one with sinking at high northern latitudes as in the present-day Atlantic, and one without. The circulation is mechanically driven, but the northern sinking can be thought of as a release valve which acts as a sink of potential energy when the surface water at high northern latitudes is dense enough to convect. While the source of energy comes from mechanical forcing, the presence or otherwise of multiple equilibria is therefore determined by thermohaline feedbacks. In some areas of parameter space an oscillation between the model’s two circulation regimes occurs, reminiscent of a bipolar seesaw.     

Transports and budgets of volume,heat, and salt from a global eddy-resolving ocean model

The results from an integration of a global ocean circulation model have been condensed into an analysis of the volume, heat, and salt transports among the major ocean basins. Transports are also broken down between the model's Ekman, thermocline, and deep layers. Overall, the model does well. Horizontal exchanges of mass, heat, and salt between ocean basins have reasonable values; and the volume of North Atlantic Deep Water (NADW) transport is in general agreement with what limited observations exist. On a global basis the zonally integrated meridional heat transport is poleward at all latitudes except for the latitude band 30°S to 45°S. This anomalous transport is most likely a signature of the model's inability to form Antarctic Intermediate (AAIW) and Antarctic bottom water (AABW) properly. Eddy heat transport is strong at the equator where its convergence heats the equatorial Pacific about twice as much as it heats the equatorial Atlantic. The greater heating in the Pacific suggests that mesoscale eddies may be a vital mechanism for warming and maintaining an upwelling portion of the global conveyor-belt circulation. The model's fresh water transport compares well with observations. However, in the Atlantic there is an excessive southward transport of fresh water due to the absence of the Mediterranean outflow and weak northward flow of AAIW. Eddies in the mid-latitudes act to redistribute heat and salt down the mean gradients. Residual fluxes calculated from a sum of the computed advective (including eddies), forced, and stored fluxes of heat and salt represent transport mostly due to vertical sub-grid scale mixing processes. Perhaps the model's greatest weakness is the lack of strong AAIW and AABW circulation cells. Accurate thermohaline forcing in the North Atlantic (based on numerous hydrographic observations) helps the model adequately produce NADW. In contrast, the southern ocean is an area of sparse observation. Better thermohaline observations in this area may be needed if models such as this are to produce the deep convection that will achieve more accurate simulations of the global 3-dimensional circulation.     

Century-scale variability in a randomly forced,two-dimensional thermohaline ocean circulation model

The response of a two-dimensional thermohaline ocean circulation model to a random freshwater flux superimposed on the usual mixed boundary conditions for temperature and salinity is considered. It is shown that for a wide range of vertical and horizontal diffusivities and a box geometry that approximates the Atlantic Ocean, 200–300 yr period oscillations exist in the basic-state, interhemispheric meridional overturning circulation with deep convection in the north. These fluctuations can also be described in terms of propagating salinity anomalies which travel in the direction of the thermohaline flow. For large horizontal (K _h = 15 × 103 m²/s) and small vertical (K _v = 0.5 × 10^–4 m²/s) diffusivities, the random forcing also excites deca-millennial oscillations in the basic structure of the thermohaline circulation. In this case, the meridional circulation pattern slowly oscillates between three different stages: a large positive cell, with deep convection in the North Atlantic and upwelling in the south; a symmetric two-cell circulation, with deep convection in both polar regions and upwelling near the equator; and a large negative cell, with deep convection in the South Atlantic and upwelling in the north. Each state can persist for 0 (10 kyr).     

温盐环流与全球增暖的数值模拟 (一)纬向平均温盐环流的模拟

本文是用简单海一气耦合模型模拟温盐环流在全球增暖事件中作用的研究工作的第一部分。为了建立一个简单海一气耦合模型，我们首先根据Ｗｒｉｇｈｔ和Ｓｔｏｋｅｒ等人的设计复制出一个包括大西洋、太平洋和南大洋在内的二维温盐环流模式，从等温、等盐和无运动的初始状态出发，在给定的年平均海表强迫下将模式积分了４０００年，模拟出了和原作相似的温盐环流。对模拟结果的分析表明，相对于北太平洋而言，北大西洋北部的高盐、低温特点（后者是由两大洋在地理上的差别决定的）是形成当代温盐环流的主要原因；从与温盐环流相联系的海表热通量来看，北大西洋北部是向大气提供热量的主要源地；模式温盐环流对于海表盐度通量的敏感性试验的结果表明，对于纬圈平均的二维模式而言，要想模拟出合理的温盐环流就必须人为地提高北大西洋北部的海表盐度，文章分析了这种作法的物理根据；模式中的对流过程对于温盐环流的维持是至关重要的，对比有无季节循环的试验结果可以看出，虽然温度场的明显的季节变化只出现在模式的最上面两层，但由于引进季节循环后冬季高纬海洋的对流活动加强，后者直接影响到温盐环流，使更多的深海热量上传并向大气释放。这是使海洋温跃层得以保持合理．厚度的一个重要原因。     

Natural forcing of climate during the last millennium: fingerprint of solar variability

The variability of the climate during the last millennium is partly forced by changes in total solar irradiance (TSI). Nevertheless, the amplitude of these TSI changes is very small so that recent reconstruction data suggest that low frequency variations in the North Atlantic Oscillation (NAO) and in the thermohaline circulation may have amplified, in the North Atlantic sector and mostly in winter, the radiative changes due to TSI variations. In this study we use a state-of-the-art climate model to simulate the last millennium. We find that modelled variations of surface temperature in the Northern Hemisphere are coherent with existing reconstructions. Moreover, in the model, the low frequency variability of this mean hemispheric temperature is found to be correlated at 0.74 with the solar forcing for the period 1001?C1860. Then, we focus on the regional climatic fingerprint of solar forcing in winter and find a significant relationship between the low frequency TSI forcing and the NAO with a time lag of more than 40?years for the response of the NAO. Such a lag is larger than the around 20-year lag suggested in other studies. We argue that this lag is due, in the model, to a northward shift of the tropical atmospheric convection in the Pacific Ocean, which is maximum more than four decades after the solar forcing increase. This shift then forces a positive NAO through an atmospheric wave connection related to the jet-stream wave guide. The shift of the tropical convection is due to the persistence of anomalous warm SST forcing the anomalous precipitation, associated with the advection of warm SST by the North Pacific subtropical gyre in a few decades. Finally, we analyse the response of the Atlantic meridional overturning circulation to solar forcing and find that the former is weakened when the latter increases. Changes in wind stress, notably due to the NAO, modify the barotropic streamfunction in the Atlantic 50?years after solar variations. This implies a wind-driven modification of the oceanic circulation in the Atlantic sector in response to changes in solar forcing, in addition to the variations of the thermohaline circulation.     

全球海气耦合模式中热盐环流对大气强迫的响应

大气环流与热盐环流 (THC)变化之间的因果关系 ,是海气相互作用研究领域的一个悬而未决的问题。作者利用一个全球海气耦合模式 -挪威卑尔根气候模式 (BCM)的 3 0 0a积分结果 ,讨论了冬季北大西洋涛动 (NAO)对海洋的强迫与热盐环流的年际调整之间的关系。结果发现 ,在NAO活动的正位相 ,伴随着中纬度西风带的加强 ,北大西洋拉布拉多海热通量损失剧增 ,同时海表盐度出现正距平 ,二者的共同作用 ,令表层海水变沉、密度增大 ,海洋层结出现不稳定 ,导致深对流发生。在NAO活动达到最强劲状态之后 3个月 ,拉布拉多海对流也达到最深。北大西洋热盐环流强度变化对拉布拉多海对流活动的响应 ,要滞后 3a左右。而在年际尺度上 ,大西洋的极向热输送变化和热盐环流的变化则基本是同步的。对流活动对大气存在明显的反馈作用。在对流活动深度达到最大之后 1～4个月 ,对流热释放令拉布拉多海表层气温明显升高     

Weak oceanic heat transport as a cause of the instability of glacial climates

The stability of the thermohaline circulation of modern and glacial climates is compared with the help of a two dimensional ocean—atmosphere—sea ice coupled model. It turns out to be more unstable as less freshwater forcing is required to induce a polar halocline catastrophy in glacial climates. The large insulation of the ocean by the extensive sea ice cover changes the temperature boundary condition and the deepwater formation regions moves much further South. The nature of the instability is of oceanic origin, identical to that found in ocean models under mixed boundary conditions. With similar strengths of the oceanic circulation and rates of deep water formation for warm and cold climates, the loss of stability of the cold climate is due to the weak thermal stratification caused by the cooling of surface waters, the deep water temperatures being regulated by the temperature of freezing. Weaker stratification with similar overturning leads to a weakening of the meridional oceanic heat transport which is the major negative feedback stabilizing the oceanic circulation. Within the unstable regime periodic millennial oscillations occur spontaneously. The climate oscillates between a strong convective thermally driven oceanic state and a weak one driven by large salinity gradients. Both states are unstable. The atmosphere of low thermal inertia is carried along by the oceanic overturning while the variation of sea ice is out of phase with the oceanic heat content. During the abrupt warming events that punctuate the course of a millennial oscillation, sea ice variations are shown respectively to damp (amplify) the amplitude of the oceanic (atmospheric) response. This sensitivity of the oceanic circulation to a reduced concentration of greenhouse gases and to freshwater forcing adds support to the hypothesis that the millennial oscillations of the last glacial period, the so called Dansgaard—Oeschger events, may be internal instabilities of the climate system.     

The relation of meridional pressure gradients to North Atlantic deep water volume transport in an ocean general circulation model

We use a coarse resolution ocean general circulation model to study the relation between meridional pressure and density gradients in the Southern Ocean and North Atlantic and the Atlantic meridional overturning circulation. In several experiments, we artificially modify the meridional density gradients by applying different magnitudes of the Gent–McWilliams isopycnal eddy diffusion coefficients in the Southern Ocean and in the North Atlantic and investigate the response of the simulated Atlantic meridional overturning to such changes. The simulations are carried out close to the limit of no diapycnal mixing, with a very small explicit vertical diffusivity and a tracer advection scheme with very low implicit diffusivities. Our results reveal that changes in eddy diffusivities in the North Atlantic affect the maximum of the Atlantic meridional overturning, but not the outflow of North Atlantic Deep Water into the Southern Ocean. In contrast, changes in eddy diffusivities in the Southern Ocean affect both the South Atlantic outflow of North Atlantic Deep Water and the maximum of the Atlantic meridional overturning. Results from these experiments are used to investigate the relation between meridional pressure gradients and the components of the Atlantic meridional overturning. Pressure gradients and overturning are found to be linearly related. We show that, in our simulations, zonally averaged deep pressure gradients are very weak between 20°S and about 30°N and that between 30°N and 60°N the zonally averaged pressure grows approximately linearly with latitude. This pressure difference balances a westward geostrophic flow at 30–40°N that feeds the southbound deep Atlantic western boundary current. We extend our analysis to a large variety of experiments in which surface freshwater forcing, vertical mixing and winds are modified. In all experiments, the pycnocline depth, assumed to be the relevant vertical scale for the northward volume transport in the Atlantic, is found to be approximately constant, at least within the coarse vertical resolution of the model. The model behaviour hence cannot directly be related to conceptual models in which changes in the pycnocline depth determine the strength of Atlantic meridional flow, and seems conceptually closer to Stommel’s box model. In all our simulations, the Atlantic overturning seems to be mainly driven by Southern Ocean westerlies. However, the actual strength of the Atlantic meridional overturning is not determined solely by the Southern Ocean wind stress but as well by the density/pressure gradients created between the deep water formation regions in the North Atlantic and the inflow/outflow region in the South Atlantic.     

Revisiting meridional overturing bistability using a minimal set of state variables: stochastic theory

A stochastic analytical model of the Atlantic meridional overturning circulation (AMOC) is presented and tested against climate model data. AMOC stability is characterised by an underlying deterministic differential equation describing the evolution of the central state variable of the system, the average Atlantic salinity. Stability of an equilibrium implies that infinitesimal salinity perturbations are damped, and violation of this requirement yields a range of unoccupied salinity states. The range of states is accurately predicted by the analytical model for a coupled climate model of intermediate complexity. The introduction of climatic noise yields an equation describing the evolution of the probability density function of the state variable, and therefore the AMOC. Given the hysteresis behaviour of the steady AMOC versus surface freshwater forcing, the statistical model is able to describe the variability of the AMOC based on knowledge of the variability in the forcing. The method accurately describes the wandering between AMOC-On and AMOC-Off states in the climate model. The framework presented is a first step in relating the stability of the AMOC to more observable aspects of its behaviour, such as its transient response to variable forcing.