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.     

Simulated decadal oscillations of the Atlantic meridional overturning circulation in a cold climate state

The significance of the Atlantic meridional overturning circulation (MOC) for regional and hemispheric climate change requires a complete understanding using fully coupled climate models. Here we present a persistent, decadal oscillation in a coupled atmosphere–ocean general circulation model. While the present study is limited by the lack of comparisons with paleo-proxy records, the purpose is to reveal a new theoretically interesting solution found in the fully-coupled climate model. The model exhibits two multi-century-long stable states with one dominated by decadal MOC oscillations. The oscillations involve an interaction between anomalous advective transport of salt and surface density in the North Atlantic subpolar gyre. Their time scale is fundamentally determined by the advection. In addition, there is a link between the MOC oscillations and North Atlantic Oscillation (NAO)-like sea level pressure anomalies. The analysis suggests an interaction between the NAO and an anomalous subpolar gyre circulation in which sea ice near and south of the Labrador Sea plays an important role in generating a large local thermal anomaly and a meridional temperature gradient. The latter induces a positive feedback via synoptic eddy activity in the atmosphere. In addition, the oscillation only appears when the Nordic Sea is completely covered by sea ice in winter, and deep convection is active only near the Irminger Sea. Such conditions are provided by a substantially colder North Atlantic climate than today.     

Factors influencing anthropogenic carbon dioxide uptake in the North Atlantic in models of the ocean carbon cycle

The uptake and storage of anthropogenic carbon in the North Atlantic is investigated using different configurations of ocean general circulation/carbon cycle models. We investigate how different representations of the ocean physics in the models, which represent the range of models currently in use, affect the evolution of CO₂ uptake in the North Atlantic. The buffer effect of the ocean carbon system would be expected to reduce ocean CO₂ uptake as the ocean absorbs increasing amounts of CO₂. We find that the strength of the buffer effect is very dependent on the model ocean state, as it affects both the magnitude and timing of the changes in uptake. The timescale over which uptake of CO₂ in the North Atlantic drops to below preindustrial levels is particularly sensitive to the ocean state which sets the degree of buffering; it is less sensitive to the choice of atmospheric CO₂ forcing scenario. Neglecting physical climate change effects, North Atlantic CO₂ uptake drops below preindustrial levels between 50 and 300 years after stabilisation of atmospheric CO₂ in different model configurations. Storage of anthropogenic carbon in the North Atlantic varies much less among the different model configurations, as differences in ocean transport of dissolved inorganic carbon and uptake of CO₂ compensate each other. This supports the idea that measured inventories of anthropogenic carbon in the real ocean cannot be used to constrain the surface uptake. Including physical climate change effects reduces anthropogenic CO₂ uptake and storage in the North Atlantic further, due to the combined effects of surface warming, increased freshwater input, and a slowdown of the meridional overturning circulation. The timescale over which North Atlantic CO₂ uptake drops to below preindustrial levels is reduced by about one-third, leading to an estimate of this timescale for the real world of about 50 years after the stabilisation of atmospheric CO₂. In the climate change experiment, a shallowing of the mixed layer depths in the North Atlantic results in a significant reduction in primary production, reducing the potential role for biology in drawing down anthropogenic CO₂.     

Atlantic multidecadal oceanic variability and its influence on the atmosphere in a climate model

The mechanisms controlling the decadal to multidecadal variability of the Atlantic Meridional Overturning Circulation (MOC) and its influence on the atmosphere are investigated using a control simulation with the IPSL-CM4 climate model. The multidecadal fluctuations of the MOC are mostly driven by deep convection in the subpolar gyre, which occurs south of Iceland in the model. The latter is primarily influenced by the anomalous advection of salinity due to changes in the East Atlantic Pattern (EAP), which is the second mode of atmospheric variability in the North Atlantic region. The North Atlantic Oscillation is the dominant mode, but it plays a secondary role in the MOC fluctuations. During summer, the MOC variability is shown to have a significant impact on the atmosphere in the North Atlantic–European sector. The MOC influence is due to an interhemispheric sea surface temperature (SST) anomaly with opposite signs in the two hemispheres but largest amplitude in the northern one. The SST pattern driven by the MOC mostly resembles the model Atlantic Multidecadal Oscillation (AMO) and bears some similarity with the observed one. It is shown that the AMO reflects both the MOC influence and the local atmospheric forcing. Hence, the MOC influence on climate is best detected using lagged relations between climatic fields. The atmospheric response resembles the EAP, in a phase that might induce a weak positive feedback on the MOC.     

Interdecadal climate variability in the subpolar North Atlantic

The statistical relationships between various components of the subpolar North Atlantic air-sea-ice climate system are reexamined in order to investigate potential processes involved in interdecadal climate variability. It is found that sea surface temperature anomalies concentrated in the Labrador Sea region have a strong impact upon atmospheric sea level pressure anomalies over Greenland, which in turn influence the transport of freshwater and ice anomalies out of the Arctic Ocean, via Fram Strait. These freshwater and ice anomalies are advected around the subpolar gyre into the Labrador Sea affecting convection and the formation of Labrador Sea Water. This has an impact upon the transport of North Atlantic Current water into the subpolar gyre and thus, also upon sea surface temperatures in the region. An interdecadal negative feedback loop is therefore proposed as an internal source of climate variability within the subpolar North Atlantic. Through the lags associated with the correlations between different climatic components, observed horizontal advection time scales, and the use of Boolean delay equation models, the time scale for one cycle of this feedback loop is determined to have a period of about 21 years.     

The UVic earth system climate model: Model description,climatology, and applications to past,present and future climates

Abstract

A new earth system climate model of intermediate complexity has been developed and its climatology compared to observations. The UVic Earth System Climate Model consists of a three‐dimensional ocean general circulation model coupled to a thermodynamic/dynamic sea‐ice model, an energy‐moisture balance atmospheric model with dynamical feedbacks, and a thermomechanical land‐ice model. In order to keep the model computationally efficient a reduced complexity atmosphere model is used. Atmospheric heat and freshwater transports are parametrized through Fickian diffusion, and precipitation is assumed to occur when the relative humidity is greater than 85%. Moisture transport can also be accomplished through advection if desired. Precipitation over land is assumed to return instantaneously to the ocean via one of 33 observed river drainage basins. Ice and snow albedo feedbacks are included in the coupled model by locally increasing the prescribed latitudinal profile of the planetary albedo. The atmospheric model includes a parametrization of water vapour/planetary longwave feedbacks, although the radiative forcing associated with changes in atmospheric CO₂ is prescribed as a modification of the planetary longwave radiative flux. A specified lapse rate is used to reduce the surface temperature over land where there is topography. The model uses prescribed present‐day winds in its climatology, although a dynamical wind feedback is included which exploits a latitudinally‐varying empirical relationship between atmospheric surface temperature and density. The ocean component of the coupled model is based on the Geophysical Fluid Dynamics Laboratory (GFDL) Modular Ocean Model 2.2, with a global resolution of 3.6° (zonal) by 1.8° (meridional) and 19 vertical levels, and includes an option for brine‐rejection parametrization. The sea‐ice component incorporates an elastic‐viscous‐plastic rheology to represent sea‐ice dynamics and various options for the representation of sea‐ice thermodynamics and thickness distribution. The systematic comparison of the coupled model with observations reveals good agreement, especially when moisture transport is accomplished through advection.

Global warming simulations conducted using the model to explore the role of moisture advection reveal a climate sensitivity of 3.0°C for a doubling of CO₂, in line with other more comprehensive coupled models. Moisture advection, together with the wind feedback, leads to a transient simulation in which the meridional overturning in the North Atlantic initially weakens, but is eventually re‐established to its initial strength once the radiative forcing is held fixed, as found in many coupled atmosphere General Circulation Models (GCMs). This is in contrast to experiments in which moisture transport is accomplished through diffusion whereby the overturning is reestablished to a strength that is greater than its initial condition.

When applied to the climate of the Last Glacial Maximum (LGM), the model obtains tropical cooling (30°N‐30°S), relative to the present, of about 2.1°C over the ocean and 3.6°C over the land. These are generally cooler than CLIMAP estimates, but not as cool as some other reconstructions. This moderate cooling is consistent with alkenone reconstructions and a low to medium climate sensitivity to perturbations in radiative forcing. An amplification of the cooling occurs in the North Atlantic due to the weakening of North Atlantic Deep Water formation. Concurrent with this weakening is a shallowing of, and a more northward penetration of, Antarctic Bottom Water.

Climate models are usually evaluated by spinning them up under perpetual present‐day forcing and comparing the model results with present‐day observations. Implicit in this approach is the assumption that the present‐day observations are in equilibrium with the present‐day radiative forcing. The comparison of a long transient integration (starting at 6 KBP), forced by changing radiative forcing (solar, CO₂, orbital), with an equilibrium integration reveals substantial differences. Relative to the climatology from the present‐day equilibrium integration, the global mean surface air and sea surface temperatures (SSTs) are 0.74°C and 0.55°C colder, respectively. Deep ocean temperatures are substantially cooler and southern hemisphere sea‐ice cover is 22% greater, although the North Atlantic conveyor remains remarkably stable in all cases. The differences are due to the long timescale memory of the deep ocean to climatic conditions which prevailed throughout the late Holocene. It is also demonstrated that a global warming simulation that starts from an equilibrium present‐day climate (cold start) underestimates the global temperature increase at 2100 by 13% when compared to a transient simulation, under historical solar, CO₂ and orbital forcing, that is also extended out to 2100. This is larger (13% compared to 9.8%) than the difference from an analogous transient experiment which does not include historical changes in solar forcing. These results suggest that those groups that do not account for solar forcing changes over the twentieth century may slightly underestimate (～3% in our model) the projected warming by the year 2100.     

Quantifying the AMOC feedbacks during a 2×CO2 stabilization experiment with land-ice melting

The response of the Atlantic Meridional Overturning Circulation (AMOC) to an increase in atmospheric CO₂ concentration is analyzed using the IPSL-CM4 coupled ocean–atmosphere model. Two simulations are integrated for 70 years with 1%/year increase in CO₂ concentration until 2×CO₂, and are then stabilized for further 430 years. The first simulation takes land-ice melting into account, via a simple parameterization, which results in a strong freshwater input of about 0.13 Sv at high latitudes in a warmer climate. During this scenario, the AMOC shuts down. A second simulation does not include this land-ice melting and herein, the AMOC recovers after 200 years. This behavior shows that this model is close to an AMOC shutdown threshold under global warming conditions, due to continuous input of land-ice melting. The analysis of the origin of density changes in the Northern Hemisphere convection sites allows an identification as to the origin of the changes in the AMOC. The processes that decrease the AMOC are the reduction of surface cooling due to the reduction in the air–sea temperature gradient as the atmosphere warms and the local freshening of convection sites that results from the increase in local freshwater forcing. Two processes also control the recovery of the AMOC: the northward advection of positive salinity anomalies from the tropics and the decrease in sea-ice transport through the Fram Strait toward the convection sites. The quantification of the AMOC related feedbacks shows that the salinity related processes contribute to a strong positive feedback, while feedback related to temperature processes is negative but remains small as there is a compensation between heat transport and surface heat flux in ocean–atmosphere coupled model. We conclude that in our model, AMOC feedbacks amplify land-ice melting perturbation by 2.5.     

Role of the upper ocean structure in the response of ENSO-like SST variability to global warming

The response of El Niño and Southern Oscillation (ENSO)-like variability to global warming varies comparatively between the two different climate system models, i.e., the Meteorological Research Institute (MRI) and Geophysical Fluid Dynamics Laboratory (GFDL) Coupled General Circulation Models (CGCMs). Here, we examine the role of the simulated upper ocean temperature structure in the different sensitivities of the simulated ENSO variability in the models based on the different level of CO₂ concentrations. In the MRI model, the sea surface temperature (SST) undergoes a rather drastic modification, namely a tendency toward a permanent El Niño-like state. This is associated with an enhanced stratification which results in greater ENSO amplitude for the MRI model. On the other hand, the ENSO simulated by GFDL model is hardly modified although the mean temperature in the near surface layer increases. In order to understand the associated mechanisms we carry out a vertical mode decomposition of the mean equatorial stratification and a simplified heat balance analysis using an intermediate tropical Pacific model tuned from the CGCM outputs. It is found that in the MRI model the increased stratification is associated with an enhancement of the zonal advective feedback and the non-linear advection. In the GFDL model, on the other hand, the thermocline variability and associated anomalous vertical advection are reduced in the eastern equatorial Pacific under global warming, which erodes the thermocline feedback and explains why the ENSO amplitude is reduced in a warmer climate in this model. It is suggested that change in stratification associated with global warming impacts the equatorial wave dynamics in a way that enhances the second baroclinic mode over the gravest one, which leads to the change in feedback processes in the CGCMs. Our results illustrate that the upper ocean vertical structure simulated in the CGCMs is a key parameter of the sensitivity of ENSO-like SST variability to global warming.     

On the subarctic North Atlantic cooling due to global warming

Various ocean reanalysis data reveal that the subarctic Atlantic sea surface temperature (SST) has been cooling during the twentieth century. A similar cooling pattern is found in the doubling CO₂ experiment obtained from the CMIP3 (coupled model intercomparison project third phase) compared to the pre-industrial experiment. Here, in order to investigate the main driver of this cooling, we perform the heat budget analysis on the subarctic Atlantic upper ocean temperature. The net surface heat flux associated with the increased concentration of greenhouse gases heats the subarctic ocean surface. In the most of models, the longwave radiation, latent heat flux, and sensible heat flux exert a warming effect, and the shortwave radiation exerts a cooling effect. On the other hand, the thermal advection by the meridional current reduces the subarctic upper ocean temperature in all models. This cold advection is attributed to the weakening of the meridional overturning circulation, which is related to the reduction in the ocean surface density. In particular, greater warming of the surface air than of the sea surface results in the reduction of surface evaporation and thereby enhanced freshening of the ocean surface water, while precipitation change was smaller than evaporation change. The thermal advections by both the wind-driven Ekman current and the density-driven geostrophic current contribute to cooling in most of the models, where the heat transport by the geostrophic current tends to be larger than that by the Ekman current.     

The role of salinity in the decadal variability of the North Atlantic meridional overturning circulation

An OGCM hindcast is used to investigate the linkages between North Atlantic Ocean salinity and circulation changes during 1963–2003. The focus is on the eastern subpolar region consisting of the Irminger Sea and the eastern North Atlantic where a careful assessment shows that the simulated interannual to decadal salinity changes in the upper 1,500 m reproduce well those derived from the available record of hydrographic measurements. In the model, the variability of the Atlantic meridional overturning circulation (MOC) is primarily driven by changes in deep water formation taking place in the Irminger Sea and, to a lesser extent, the Labrador Sea. Both are strongly influenced by the North Atlantic Oscillation (NAO). The modeled interannual to decadal salinity changes in the subpolar basins are mostly controlled by circulation-driven anomalies of freshwater flux convergence, although surface salinity restoring to climatology and other boundary fluxes each account for approximately 25% of the variance. The NAO plays an important role: a positive NAO phase is associated with increased precipitation, reduced northward salt transport by the wind-driven intergyre gyre, and increased southward flows of freshwater across the Greenland–Scotland ridge. Since the NAO largely controlled deep convection in the subpolar gyre, fresher waters are found near the sinking region during convective events. This markedly differs from the active influence on the MOC that salinity exerts at decadal and longer timescales in most coupled models. The intensification of the MOC that follows a positive NAO phase by about 2 years does not lead to an increase in the northward salt transport into the subpolar domain at low frequencies because it is cancelled by the concomitant intensification of the subpolar gyre which shifts the subpolar front eastward and reduces the northward salt transport by the North Atlantic Current waters. This differs again from most coupled models, where the gyre intensification precedes that of the MOC by several years.     

Abrupt millennial variability and interdecadal-interstadial oscillations in a global coupled model: sensitivity to the background climate state

The origin and bifurcation structure of abrupt millennial-scale climate transitions under steady external solar forcing and in the absence of atmospheric synoptic variability is studied by means of a global coupled model of intermediate complexity. We show that the origin of Dansgaard-Oeschger type oscillations in the model is caused by the weaker northward oceanic heat transport in the Atlantic basin. This is in agreement with previous studies realized with much simpler models, based on highly idealized geometries and simplified physics. The existence of abrupt millennial-scale climate transitions during glacial times can therefore be interpreted as a consequence of the weakening of the negative temperature-advection feedback. This is confirmed through a series of numerical experiments designed to explore the sensitivity of the bifurcation structure of the Atlantic meridional overturning circulation to increased atmospheric CO₂ levels under glacial boundary conditions. Contrasting with the cold, stadial, phases of millennial oscillations, we also show the emergence of strong interdecadal variability in the North Atlantic sector during warm interstadials. The instability driving these interdecadal-interstadial oscillations is shown to be identical to that found in ocean-only models forced by fixed surface buoyancy fluxes, that is, a large-scale baroclinic instability developing in the vicinity of the western boundary current in the North Atlantic. Comparisons with modern observations further suggest a physical mechanism similar to that driving the 30–40?years time scale associated with the Atlantic multidecadal oscillation.     

The sensitivity of the meridional overturning circulation to modelling uncertainty in a perturbed physics ensemble without flux adjustment

Ocean dynamics play a key role in the climate system, by redistributing heat and freshwater. The uncertainty of how these processes are represented in climate models, and how this uncertainty affects future climate projections can be investigated using perturbed physics ensembles of global circulation models (GCMs). Techniques such as flux adjustments should be avoided since they can impact the sensitivity of the ensemble to the imposed forcing. In this study a method for developing an coupled ensemble with a GCM that does not use flux adjustment is presented. The ensemble is constrained by using information from a prior ensemble with a mixed layer ocean coupled to an atmosphere GCM, to reduce drifts in the coupled ensemble. Constraints on parameter perturbations are derived by using observational constraints on surface temperature, and top of the atmosphere radiative fluxes. As an example of such an ensemble developed with this methodology, uncertainty in response of the meridional overturning circulation (MOC) to increased CO₂ concentrations is investigated. The ensemble mean MOC strength is 17.1?Sv and decreases by 2.1?Sv when greenhouse gas concentrations are doubled. No rapid changes or shutdown of the MOC are seen in any of the ensemble members. There is a strong negative relationship between global mean temperature and MOC strength across the ensemble which is not seen in a multimodel ensemble. A positive relationship between climate sensitivity and the decrease of MOC strength is also seen.     

Sensitivity of North Atlantic subpolar gyre and overturning to stratification-dependent mixing: response to global warming

We use a reduced complexity climate model with a three-dimensional ocean component and realistic topography to investigate the effect of stratification-dependent mixing on the sensitivity of the North Atlantic subpolar gyre (SPG), and the Atlantic meridional overturning circulation (AMOC), to idealized CO₂ increase and peaking scenarios. The vertical diffusivity of the ocean interior is parameterized as κ ∼ N ^−α, where N is the local buoyancy frequency. For all parameter values 0 ≤ α ≤ 3, we find the SPG, and subsequently the AMOC, to weaken in response to increasing CO₂ concentrations. The weakening is significantly stronger for α ≥ α_cr ≈ 1.5. Depending on the value of α, two separate model states develop. These states remain different after the CO₂ concentration is stabilized, and in some cases even after the CO₂ concentration has been decreased again to the pre-industrial level. This behaviour is explained by a positive feedback between stratification and mixing anomalies in the Nordic Seas, causing a persistent weakening of the SPG.     

CLIMBER-2: a climate system model of intermediate complexity. Part II: model sensitivity

 A set of sensitivity experiments with the climate system model of intermediate complexity CLIMBER-2 was performed to compare its sensitivity to changes in different types of forcings and boundary conditions with the results of comprehensive models (GCMs). We investigated the climate system response to changes in freshwater flux into the Northern Atlantic, CO₂ concentration, solar insolation, and vegetation cover in the boreal zone and in the tropics. All these experiments were compared with the results of corresponding experiments performed with different GCMs. Qualitative, and in many respects, quantitative agreement between the results of CLIMBER-2 and GCMs demonstrate the ability of our climate system model of intermediate complexity to address diverse aspects of the climate change problem. In addition, we used our model for a series of experiments to assess the impact of some climate feedbacks and uncertainties in model parameters on the model sensitivity to different forcings. We studied the role of freshwater feedback and vertical ocean diffusivity for the stability properties of the thermohaline ocean circulation. We show that freshwater feedback plays a minor role, while changes of vertical diffusivity in the ocean considerably affect the circulation stability. In global warming experiments we analysed the impact of hydrological sensitivity and vertical diffusivity on the long-term evolution of the thermohaline circulation. In the boreal and tropical deforestation experiments we assessed the role of an interactive ocean and showed that for both types of deforestation scenarios, an interactive ocean leads to an additional cooling due to albedo and water vapour feedbacks. Received: 28 May 2000 / Accepted: 9 November 2000     

Dependence of NAO variability on coupling with sea ice

The variance of the North Atlantic Oscillation index (denoted N) is shown to depend on its coupling with area-averaged sea ice concentration anomalies in and around the Barents Sea (index denoted B). The observed form of this coupling is a negative feedback whereby positive N tends to produce negative B, which in turn forces negative N. The effects of this feedback in the system are examined by modifying the feedback in two modeling frameworks: a statistical vector autoregressive model (F _VAR) and an atmospheric global climate model (F _CAM) customized so that sea ice anomalies on the lower boundary are stochastic with adjustable sensitivity to the model??s evolving N. Experiments show that the variance of N decreases nearly linearly with the sensitivity of B to N, where the sensitivity is a measure of the negative feedback strength. Given that the sea ice concentration field has anomalies, the variance of N goes down as these anomalies become more sensitive to N. If the sea ice concentration anomalies are entirely absent, the variance of N is even smaller than the experiment with the most sensitive anomalies. Quantifying how the variance of N depends on the presence and sensitivity of sea ice anomalies to N has implications for the simulation of N in global climate models. In the physical system, projected changes in sea ice thickness or extent could alter the sensitivity of B to N, impacting the within-season variability and hence predictability of N.     

The role of the individual air-sea flux components in CO2-induced changes of the ocean's circulation and climate

 In this study we investigate the role of heat, freshwater and momentum fluxes in changing the oceanic climate and thermohaline circulation as a consequence of increasing atmospheric CO₂ concentration. Two baseline integrations with a fully coupled ocean atmosphere general circulation model with either fixed or increasing atmospheric CO₂ concentrations have been performed. In a set of sensitivity experiments either freshwater (precipitation, evaporation and runoff from the continents) and/or momentum fluxes were no longer simulated, but prescribed according to one of the fully coupled baseline experiments. This approach gives a direct estimate of the contribution from the individual flux components. The direct effect of surface warming and the associated feedbacks in ocean circulation are the dominant processes in weakening the Atlantic thermohaline circulation in our model. The relative contribution of momentum and freshwater fluxes to the total response turned out to be less than 25%, each. Changes in atmospheric water vapour transport lead to enhanced freshwater input into middle and high latitudes, which weakens the overturning. A stronger export of freshwater from the Atlantic drainage basin to the Indian and Pacific ocean, on the other hand, intensifies the Atlantic overturning circulation. In total the modified freshwater fluxes slightly weaken the Atlantic thermohaline circulation. The contribution of the modified momentum fluxes has a similar magnitude, but enhances the formation of North Atlantic deep water. Salinity anomalies in the Atlantic as a consequence of greenhouse warming stem in almost equal parts from changes in net freshwater fluxes and from changes in ocean circulation caused by the surface warming due to atmospheric heat fluxes. Important effects of the momentum fluxes are a poleward shift of the front between Northern Hemisphere subtropical and subpolar gyres and a southward shift in the position of the Antarctic circumpolar current, with a clear signal in sea level. Received: 3 May 1999 / Accepted: 11 December 1999     

Impacts of interruption of the Agulhas leakage on the tropical Atlantic in coupled ocean�Catmosphere simulations

In this paper we use a coupled ocean?Catmosphere model to investigate the impact of the interruption of Agulhas leakage of Indian ocean water on the tropical Atlantic, a region where strong coupled ocean?Catmosphere interactions occur. The effect of a shut down of leakage of Indian ocean water is isolated from the effect of a collapse of the MOC. In our experiments, the ocean model is forced with boundary conditions in the southeastern corner of the domain that correspond to no interocean exchange of Indian ocean water into the Atlantic. The southern boundary condition is taken from the Levitus data and ensures an MOC in the Atlantic. Within this configuration, instead of warm and salty Indian ocean water temperature (cold) and salinity (fresh) anomalies of southern ocean origin propagate into the South Atlantic and eventually reach the equatorial region, mainly in the thermocline. This set up mimics the closure of the ??warm water path?? in favor of the ??cold water path??. As part of the atmospheric response, there is a northward shift of the intertropical convergence zone (ITCZ). The changes in trade winds lead to reduced Ekman pumping in the equatorial region. This leads to a freshening and warming of the surface waters along the equator. Especially in the Cold Tongue region, the cold and fresh subsurface anomalies do not reach the surface due to the reduced upwelling. The anomaly signals are transported by the equatorial undercurrent and spread away from the equator within the thermocline. Part of the anomaly eventually reaches the Tropical North Atlantic, where it affects the Guinea Dome. Surprisingly, the main effect at the surface is small on the equator and relatively large at the Guinea Dome. In the atmosphere, the northward shift of the ITCZ is associated with a band of negative precipitation anomalies and higher salinities over the Tropical South Atlantic. An important implication of these results is that the modified water characteristics due to a shut down of the Agulhas leakage remain largely unaffected when crossing the equatorial Atlantic and therefore can affect the deepwater formation in the North Atlantic. This supports the hypothesis that the Agulhas leakage is an important source region for climate change and decadal variability of the Atlantic.     

An ENSO stability analysis. Part II: results from the twentieth and twenty-first century simulations of the CMIP3 models

In this study, a Bjerknes stability (BJ) index, proposed by Jin et al. (2006), is adopted to assess the overall stability of El Niño and Southern Oscillation (ENSO) in state-of-the-art coupled models. The twentieth and twenty-first century simulations of 12 coupled models among the coupled model intercomparison project phase 3 models used in the intergovernmental panel on climate change forth assessment report demonstrate a significant positive correlation between ENSO amplitude and ENSO stability as measured by the BJ index. The simulations also show a diversity of behavior regarding the ENSO stability among the coupled models, which can be attributed to different mean state and sensitivity of an oceanic and atmospheric response to wind and SST forcing from model to model. When respective components of the BJ index obtained from the coupled models are compared with those from observations, it is revealed that most coupled models underestimate the thermodynamic damping effect and the positive effect of the zonal advective and thermocline feedback. Under increased CO₂ induced warm climate, changes, relative to the twentieth century simulations, in the damping and feedback terms responsible for the ENSO stability measured by the BJ index can be linked to mean state changes and associated atmospheric and oceanic response sensitivity changes. There is a clear multi-model trend in the damping terms and positive zonal advective feedback, thermocline feedback, and Ekman feedback terms under enhanced greenhouse gas conditions. However, the various behavior among the coupled models in competition between the positive feedback and negative damping terms in the BJ index formula prevent the formation of a definitive conclusion regarding future projections of ENSO stability using the current coupled models.     

The effect of sea-ice extent in the North Atlantic on the stability of the thermohaline circulation in global warming experiments

Different climate models simulate different behavior of the Atlantic meridional overturning circulation (MOC) under the same global warming scenario. We propose a plausible explanation for this and argue that a proper simulation of the present-day climate in the subpolar North Atlantic is important. This is illustrated using results from idealized global warming experiments, in which both the radiative forcing scenario and the model employed are the same, with the only major difference being the initial subpolar North Atlantic climate. The initial conditions are made progressively colder, with more extensive sea-ice cover in the northern North Atlantic.The key result is that starting from conditions which are too cold in the North Atlantic and with sea-ice that is too extensive leads to an MOC that is more stable to the radiative forcing. Furthermore, under considerably underestimated sea surface temperatures in subpolar regions, the MOC can even intensify. A reduction of freshwater flux associated with the reduction of sea-ice melt is shown to be important for such unusual behavior of the MOC. Other mechanisms are also considered, but not deemed as important in explaining published inter-model differences.