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     

Stabilization of thermohaline circulation by wind-driven and vertical diffusive salt transport

 The stability of the thermohaline circulation is investigated using an ocean general circulation model coupled to a simple atmospheric model. The atmospheric model is so developed that it represents the wind stress and the freshwater flux more realistically than existing energy balance models. The coupled model can reproduce the realistic deep ocean circulation without any flux adjustment. Effects of the wind stress and the vertical diffusion on the thermohaline circulation are studied by conducting various experiments with the coupled model. The Ekman upwelling between 60^∘N and 90^∘N brings up salt to the sea surface, while the compensation flow of the Ekman transport and the wind-driven gyre circulation between 30^∘N and 60^∘N carry salt horizontally to the high latitudes. By carrying out experiments where the wind stress is completely or partly removed, it is demonstrated that either of the vertical or the horizontal salt transport prevents the halocline formation at high latitudes and maintains the thermohaline circulation. For an experiment in which the vertical diffusivity is enhanced at high latitudes, it is shown that the vertical diffusion at high latitudes also prevents the halocline formation and stabilizes the thermohaline circulation. It is also shown that the value of the vertical diffusivity at high latitude affects the existence of the multiple equilibria of the thermohaline circulation. Received: 26 April 2000 / Accepted: 10 January 2001     

Multi-decadal thermohaline variability in an ocean–atmosphere general circulation model

A century scale integration of a near-global atmosphere–ocean model is used to study the multi-decadal variability of the thermohaline circulation (THC) in the Atlantic. The differences between the coupled and two supplementary ocean-only experiments suggest that a significant component of this variability is controlled by either a collective behavior of the ocean and the atmosphere, particularly in the form of air-sea heat exchange, or sub-monthly random noise present in the coupled system. Possible physical mechanisms giving rise to the mode of this THC variability are discussed. The SST anomaly associated with the THC variability resembles an interdecadal SST pattern extracted from observational data, as well as a pattern associated with the 50–60 year THC variability in the GFDL coupled model. In each case, a warming throughout the subpolar North Atlantic but concentrated along the Gulf Stream and its extension is indicated when the THC is strong. Concomitantly, surface air temperature has positive anomalies over the warmer ocean, with the strongest signal located downwind of the warmest SST anomalies and intruding into the western Eurasian Continent. In addition to the thermal response, there are also changes in the atmospheric flow pattern. More specifically, an anomalous northerly wind develops over the Labrador Sea when the THC is stronger than normal, suggesting a local primacy of the atmospheric forcing in the thermohaline perturbation structure.     

Predictability of SST forced climate signals in two atmospheric general circulation models

 The predictability of atmospheric responses to global sea surface temperature (SST) anomalies is evaluated using ensemble simulations of two general circulation models (GCMs): the GENESIS version 1.5 (GEN) and the ECMWF cycle 36 (ECM). The integrations incorporate observed SST variations but start from different initial land and atmospheric states. Five GEN 1980–1992 and six ECM 1980–1988 realizations are compared with observations to distinguish predictable SST forced climate signals from internal variability. To facilitate the study, correlation analysis and significance evaluation techniques are developed on the basis of time series permutations. It is found that the annual mean global area with realistic signals is variable dependent and ranges from 3 to 20% in GEN and 6 to 28% in ECM. More than 95% of these signal areas occur between 35 °S–35 °N. Due to the existence of model biases, robust responses, which are independent of initial condition, are identified over broader areas. Both GCMs demonstrate that the sensitivity to initial conditions decreases and the predictability of SST forced responses increases, in order, from 850 hPa zonal wind, outgoing longwave radiation, 200 hPa zonal wind, sea-level pressure to 500 hPa height. The predictable signals are concentrated in the tropical and subtropical Pacific Ocean and are identified with typical El Ni?o/ Southern Oscillation phenomena that occur in response to SST and diabatic heating anomalies over the equatorial central Pacific. ECM is less sensitive to initial conditions and better predicts SST forced climate changes. This results from (1) a more realistic basic climatology, especially of the upper-level wind circulation, that produces more realistic interactions between the mean flow, stationary waves and tropical forcing; (2) a more vigorous hydrologic cycle that amplifies the tropical forcing signals, which can exceed internal variability and be more efficiently transported from the forcing region. Differences between the models and observations are identified. For GEN during El Ni?o, the convection does not carry energy to a sufficiently high altitude, while the spread of the tropospheric warming along the equator is slower and the anomaly magnitude smaller than observed. This impacts model ability to simulate realistic responses over Eurasia and the Indian Ocean. Similar biases exist in the ECM responses. In addition, the relationships between upper and lower tropospheric wind responses to SST forcing are not well reproduced by either model. The identification of these model biases leads to the conclusion that improvements in convective heat and momentum transport parametrizations and basic climate simulations could substantially increase predictive skill. Received: 25 April 1996 / Accepted: 9 December 1996     

A predictability study of simulated North Atlantic multidecadal variability

 The North Atlantic is one of the few places on the globe where the atmosphere is linked to the deep ocean through air–sea interaction. While the internal variability of the atmosphere by itself is only predictable over a period of one to two weeks, climate variations are potentially predictable for much longer periods of months or even years because of coupling with the ocean. This work presents details from the first study to quantify the predictability for simulated multidecadal climate variability over the North Atlantic. The model used for this purpose is the GFDL coupled ocean-atmosphere climate model used extensively for studies of global warming and natural climate variability. This model contains fluctuations of the North Atlantic and high-latitude oceanic circulation with variability concentrated in the 40–60 year range. Oceanic predictability is quantified through analysis of the time-dependent behavior of large-scale empirical orthogonal function (EOF) patterns for the meridional stream function, dynamic topography, 170 m temperature, surface temperature and surface salinity. The results indicate that predictability in the North Atlantic depends on three main physical mechanisms. The first involves the oceanic deep convection in the subpolar region which acts to integrate atmospheric fluctuations, thus providing for a red noise oceanic response as elaborated by Hasselmann. The second involves the large-scale dynamics of the thermohaline circulation, which can cause the oceanic variations to have an oscillatory character on the multidecadal time scale. The third involves nonlocal effects on the North Atlantic arising from periodic anomalous fresh water transport advecting southward from the polar regions in the East Greenland Current. When the multidecadal oscillatory variations of the thermohaline circulation are active, the first and second EOF patterns for the North Atlantic dynamic topography have predictability time scales on the order of 10–20 y, whereas EOF-1 of SST has predictability time scales of 5–7 y. When the thermohaline variability has weak multidecadal power, the Hasselmann mechanism is dominant and the predictability is reduced by at least a factor of two. When the third mechanism is in an extreme phase, the North Atlantic dynamic topography patterns realize a 10–20 year predictability time scale. Additional analysis of SST in the Greenland Sea, in a region associated with the southward propagating fresh water anomalies, indicates the potential for decadal scale predictability for this high latitude region as well. The model calculations also allow insight into regional variations of predictability, which might be useful information for the design of a monitoring system for the North Atlantic. Predictability appears to break down most rapidly in regions of active convection in the high-latitude regions of the North Atlantic. Received: 28 October 1996 / Accepted: 21 March 1997     

Interpretation of interbasin exchange in an isopycnal ocean model

 This work concerns an analysis of inter-basin and inter-layer exchanges in the component ocean part of the coupled ECHAM4/OPYC3 general circulation model, aimed at documenting the simulation of North Atlantic Deep Water (NADW) and related thermohaline circulations in the Indian and Pacific Oceans. The modeled NADW is formed mainly in the Greenland– Iceland–Norwegian Seas through a composite effect of deep convection and downward cross-isopycnal transport. The modeled deep-layer outflow of NADW can reach 16 Sv near 30 °S in the South Atlantic, with the corresponding upper-layer return flow mainly coming from the “cold water path” through Drake Passage. Less than 4 Sv of the Agulhas “leakage” water contributes to the replacement of NADW along the “warm water path”. In the South Atlantic Ocean, the model shows that some intermediate isopycnal layers with potential densities ranging between 27.0 and 27.5 are the major water source that compensate the NADW return flow and enhance the Circumpolar Deep Water (CDW) flowing from the Atlantic into Indian Ocean. The modeled thermohaline circulations in the Indian and Pacific Oceans indicate that the Indian Ocean may play the major role in converting deep water into intermediate water. About 16 Sv of the CDW-originating deep water enters the Indian Ocean northward of 31 °S, of which more than 13 Sv “upwell” mainly near the continental boundaries of Africa, South Asia and Australia through inter-layer exchanges and return to the Antarctic Circumpolar Current (ACC) as intermediate-layer water. As a contrast, only 4 Sv of Pacific intermediate water is connected to “upwelling” flow southward across 31 °S while the magnitude of northward deep flow across 31 °S in the Pacific Ocean is significantly greater than that in the Indian Ocean. The model suggests that a large portion of the deep waters entering the Pacific Ocean (about 14 Sv) “upwells” continually into some upper layers through the thermocline, and becomes the source of the Indonesian throughflow. Uncertainties in these results may be related to the incomplete adjustment of the model’s isopycnal layers and the sensitivity of the Indonesian throughflow to the model’s geography and topography. Received: 12 August 1997/Accepted: 12 March 1998     

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.     

Easterly wave regimes and associated convection over West Africa and tropical Atlantic: results from the NCEP/NCAR and ECMWF reanalyses

 NCEP/NCAR and ECMWF daily reanalyses are used to investigate the synoptic variability of easterly waves over West Africa and tropical Atlantic at 700 hPa in northern summer between 1979–1995 (1979–1993 for ECMWF). Spectral analysis of the meridional wind component at 700 hPa highlighted two main periodicity bands, between 3 and 5 days, and 6 and 9 days. The 3–5-day easterly wave regime has already been widely investigated, but only on shorter datasets. These waves grow both north and south of the African Easterly Jet (AEJ). The two main tracks, noted over West Africa at 5 °N and 15 °N, converge over the Atlantic on latitude 17.5 °N. These waves are more active in August–September than in June–July. Their average wavelength/phase speed varies from about 3000 km/8 m s^-1 north of the jet to 5000 km/12 m s^-1 south of the jet. Rainfall, convection and monsoon flux are significantly modulated by these waves, convection in the Inter-Tropical Convergence Zone (ITCZ) being enhanced in the trough and ahead of it, with a wide meridional extension. Compared to the 3–5-day waves, the 6–9-day regime is intermittent and the corresponding wind field pattern has both similar and contrasting characteristics. The only main track is located north of the AEJ along 17.5 °N both over West Africa and the Atlantic. The mean wavelength is higher, about 5000 km long, and the average phase speed is about 7 m s^-1. Then the wind field perturbation is mostly evident at the AEJ latitude and north of it. The perturbation structure is similar to that of 3–5-days in the north except that the more developed circulation centers, moving more to the north, lead to a large modulation of the jet zonal wind component. South of the AEJ, the wind field perturbation is weaker and quite different. The zonal wind core of the jet appears to be an almost symmetric axis in the 6–9-day wind field pattern, a clockwise circulation north of the AEJ being associated with a counter-clockwise circulation south of the jet, and vice versa. These 6–9-day easterly waves also affect significantly rainfall, convection and monsoon flux but in a different way, inducing large zonal convective bands in the ITCZ, mostly in the trough and behind it. As opposed to the 3–5-day wave regime, these rainfall anomalies are associated with anomalies of opposite sign over the Guinea coast and the Sahelian regions. Over the continent, these waves are more active in June–July, and in August–September over the ocean. GATE phase I gave an example of such an active 6–9-day wave pattern. Considered as a sequence of weak easterly wave activity, this phase was also a sequence of high 6–9-day easterly wave activity. We suggest that the 6–9-day regime results from an interaction between the 3–5-day easterly wave regime (maintained by the barotropic/baroclinic instability of the AEJ), and the development of strong anticyclonic circulations, north of the jet over West Africa, and both north and south of the jet over the Atlantic, significantly affecting the jet zonal wind component. The permanent subtropical anticyclones (Azores, Libya, St Helena) could help initiation and maintenance of such regime over West Africa and tropical Atlantic. Based on an a priori period-band criterion, our synoptic classification has enabled us to point out two statistical and meteorological easterly wave regimes over West Africa and tropical Atlantic. NCEP/NCAR and ECMWF reanalyses are in good agreement, the main difference being a more developed easterly wave activity in the NCEP/NCAR reanalyses, especially for the 3–5-day regime over the Atlantic. Received: 28 May 1998 / Accepted: 2 May 1999     

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

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

A transient climate change simulation with greenhouse gas and aerosol forcing: projected climate to the twenty-first century

 The potential climatic consequences of increasing atmospheric greenhouse gas (GHG) concentration and sulfate aerosol loading are investigated for the years 1900 to 2100 based on five simulations with the CCCma coupled climate model. The five simulations comprise a control experiment without change in GHG or aerosol amount, three independent simulations with increasing GHG and aerosol forcing, and a simulation with increasing GHG forcing only. Climate warming accelerates from the present with global mean temperatures simulated to increase by 1.7 °C to the year 2050 and by a further 2.7 °C by the year 2100. The warming is non-uniform as to hemisphere, season, and underlying surface. Changes in interannual variability of temperature show considerable structure and seasonal dependence. The effect of the comparatively localized negative radiative forcing associated with the aerosol is to retard and reduce the warming by about 0.9 °C at 2050 and 1.2 °C at 2100. Its primary effect on temperature is to counteract the global pattern of GHG-induced warming and only secondarily to affect local temperatures suggesting that the first order transient climate response of the system is determined by feedback processes and only secondarily by the local pattern of radiative forcing. The warming is accompanied by a more active hydrological cycle with increases in precipitation and evaporation rates that are delayed by comparison with temperature increases. There is an “El Nino-like” shift in precipitation and an overall increase in the interannual variability of precipitation. The effect of the aerosol forcing is again primarily to delay and counteract the GHG-induced increase. Decreases in soil moisture are common but regionally dependent and interannual variability changes show considerable structure. Snow cover and sea-ice retreat. A PNA-like anomaly in mean sea-level pressure with an enhanced Aleutian low in northern winter is associated with the tropical shift in precipitation regime. The interannual variability of mean sea-level pressure generally decreases with largest decreases in the tropical Indian ocean region. Changes to the ocean thermal structure are associated with a spin-down of the Atlantic thermohaline circulation together with a decrease in its variability. The effect of aerosol forcing, although modest, differs from that for most other quantities in that it does not act primarily to counteract the GHG forcing effect. The barotropic stream function in the ocean exhibits modest change in the north Pacific but accelerating changes in much of the Southern Ocean and particularly in the north Atlantic where the gyre spins down in conjunction with the decrease in the thermohaline circulation. The results differ in non-trivial ways from earlier equilibrium 2 × CO₂ results with the CCCma model as a consequence of the coupling to a fully three-dimensional ocean model and the evolving nature of the forcing. Received: 24 September 1998 / Accepted: 8 October 1999     

Coupled ocean-atmosphere surface variability and its climate impacts in the tropical Atlantic region

 This study examines time evolution and statistical relationships involving the two leading ocean-atmosphere coupled modes of variability in the tropical Atlantic and some climate anomalies over the tropical 120 °W–60 °W region using selected historical files (75-y near global SSTs and precipitation over land), more recent observed data (30-y SST and pseudo wind stress in the tropical Atlantic) and reanalyses from the US National Centers for Environmental Prediction (NCEP/NCAR) reanalysis System on the period 1968–1997: surface air temperature, sea level pressure, moist static energy content at 850 hPa, precipitable water and precipitation. The first coupled mode detected through singular value decomposition of the SST and pseudo wind-stress data over the tropical Atlantic (30 °N–20 °S) expresses a modulation in the thermal transequatorial gradient of SST anomalies conducted by one month leading wind-stress anomalies mainly in the tropical north Atlantic during northern winter and fall. It features a slight dipole structure in the meridional plane. Its time variability is dominated by a quasi-decadal signal well observed in the last 20–30 ys and, when projected over longer-term SST data, in the 1920s and 1930s but with shorter periods. The second coupled mode is more confined to the south-equatorial tropical Atlantic in the northern summer and explains considerably less wind-stress/SST cross-covariance. Its time series features an interannual variability dominated by shorter frequencies with increased variance in the 1960s and 1970s before 1977. Correlations between these modes and the ENSO-like Nino3 index lead to decreasing amplitude of thermal anomalies in the tropical Atlantic during warm episodes in the Pacific. This could explain the nonstationarity of meridional anomaly gradients on seasonal and interannual time scales. Overall the relationships between the oceanic component of the coupled modes and the climate anomaly patterns denote thermodynamical processes at the ocean/atmosphere interface that create anomaly gradients in the meridional plane in a way which tends to alter the north–south movement of the seasonal cycle. This appears to be consistent with the intrinsic non-dipole character of the tropical Atlantic surface variability at the interannual time step and over the recent period, but produces abnormal amplitude and/or delayed excursions of the intertropical convergence zone (ITCZ). Connections with continental rainfall are approached through three (NCEP/NCAR and observed) rainfall indexes over the Nordeste region in Brazil, and the Guinea and Sahel zones in West Africa. These indices appear to be significantly linked to the SST component of the coupled modes only when the two Atlantic modes+the ENSO-like Nino3 index are taken into account in the regressions. This suggests that thermal forcing of continental rainfall is particularly sensitive to the linear combinations of some basic SST patterns, in particular to those that create meridional thermal gradients. The first mode in the Atlantic is associated with transequatorial pressure, moist static energy and precipitable water anomaly patterns which can explain abnormal location of the ITCZ particularly in northern winter, and hence rainfall variations in Nordeste. The second mode is more associated with in-phase variations of the same variables near the southern edge of the ITCZ, particularly in the Gulf of Guinea during the northern spring and winter. It is primarily linked to the amplitude and annual phase of the ITCZ excursions and thus to rainfall variations in Guinea. Connections with Sahel rainfall are less clear due to the difficulty for the model to correctly capture interannual variability over that region but the second Atlantic mode and the ENSO-like Pacific variability are clearly involved in the Sahel climate interannual fluctuations: anomalous dry (wet) situations tend to occur when warmer (cooler) waters are present in the eastern Pacific and the gulf of Guinea in northern summer which contribute to create a northward (southward) transequatorial anomaly gradient in sea level pressure over West Africa. Received: 14 April 1998 / Accepted: 24 December 1998     

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).     

The atmospheric sulfur cycle in ECHAM-4 and its impact on the shortwave radiation

 The atmospheric general circulation model ECHAM-4 is coupled to a chemistry model to calculate sulfate mass distribution and the radiative forcing due to sulfate aerosol particles. The model simulates the main components of the hydrological cycle and, hence, it allows an explicit treatment of cloud transformation processes and precipitation scavenging. Two experiments are performed, one with pre-industrial and one with present-day sulfur emissions. In the pre-industrial emission scenario SO₂ is oxidized faster to sulfate and the in-cloud oxidation via the reaction with ozone is more important than in the present-day scenario. The atmospheric sulfate mass due to anthropogenic emissions is estimated as 0.38 Tg sulfur. The radiative forcing due to anthropogenic sulfate aerosols is calculated diagnostically. The backscattering of shortwave radiation (direct effect) as well as the impact of sulfate aerosols on the cloud albedo (indirect effect) is estimated. The model predicts a direct forcing of −0.35 W m^-2 and an indirect forcing of −0.76 W m^-2. Over the continents of the Northern Hemisphere the direct forcing amounts to −0.64 W m^-2. The geographical distribution of the direct and indirect effect is very different. Whereas the direct forcing is strongest over highly polluted continental regions, the indirect forcing over sea exceeds that over land. It is shown that forcing estimates based on monthly averages rather than on instantaneous sulfate pattern overestimate the indirect effect but have little effect on the direct forcing. Received: 16 October 1996/Accepted: 24 October 1996     

On the role of flux adjustments in an idealized coupled climate model

 The effect of employing flux adjustments on the climatic response of an idealized coupled model to an imposed radiative forcing is investigated with two coupled models, one of which employs flux adjustments. A linear reduction (to the planetary longwave flux) of 4 W/m² is applied over a 70 y period and held constant thereafter. Similar model responses are found (during the initial 70 y period) for global-scale diagnostics of hemispheric air temperature due to the nearly linear surface-air temperature response to the radiative forcing. Significant regional scale differences do exist, however. As the perturbation away from the present climate grows, basin-scale diagnostics (such as meridional overturning rates) begin to diverge between flux adjusted and non-flux adjusted models. Once the imposed radiative forcing is held constant, differences in global mean air temperature of up to 0.5 °C are found, with large regional-scale differences in air temperature and overturning rates within the North Atlantic and Southern Ocean. Two additional experiments with the flux adjusted model (beginning from points further along the control integration) suggest that the elimination of much of the coupling shock before the radiative forcing is applied leads to results slightly closer to the non-flux adjusted case, although large differences still persist. In particular a dipole structure indicating an enhanced warming within the Pacific sector of the Southern Ocean, and cooling within the Atlantic sector is not reproduced by the flux adjusted models. This disparity is intimately linked to the Southern Ocean overturning cell along with the flux adjustments employed as well as the drift arising from coupling shock. If a similar form of sensitivity exists in more realistic coupled models, our results suggest: (1) perturbation experiments should not be undertaken until after the coupled model control experiment is carried out for several hundred years (thereby minimizing the coupling shock); (2) care should be exercised in the interpretation of regional-scale results (over the ocean) in coupled models which employ flux adjustments; (3) care should also be taken in interpreting even global-scale diagnostics in flux adjusted models for large perturbations about the present climate. Received: 15 November 1996 / Accepted: 4 June 1997     

A numerical world ocean general circulation model: Part II. A baroclinic experiment

A new six-layer world ocean general circulation model based on the primitive system of equations is described in detail and its performance in the case of a homogeneous ocean is described. These test integrations show that the model is capable of reproducing the observed mean barotropic or vertically-integrated transport, as well as the seasonal variability of the major ocean gyres. The surface currents, however, are dominated by the Ekman transport, and such non-linear features as the western boundary currents and the equatorial countercurrents are poorly represented. The abyssal boundary countercurrents are also absent due to the lack of thermohaline forcing. The most conspicuous effect of the bottom topography on a homogeneous ocean is seen in the Southern ocean where the calculated Antarctic circumpolar transport through the Drake passage ( ≈ 10 Sv, with bathymetry included) greatly underestimates the observed transport (≈ 100 Sv).     

Holocene climatic change in Morocco: a quantitative reconstruction from pollen data

 Annual precipitation, July and January temperatures were reconstructed from a continuous Holocene pollen sequence from the Middle Atlas, Morocco, using the best modern analogues method. The reconstructions show a clear difference between the early and late Holocene: from ∼10 ka to ∼6.5 ka the climate was drier and warmer than during the period since 6.5 ka. The average value of annual precipitation was ∼870 mm until 6.5 ka, then rose to ∼940 mm. Between 10 ka and 6.5 ka January and July temperatures were about 4 °C higher than the present. Both temperatures show a marked decrease between 7 ka and 6 ka. After 6.5 ka July and January temperatures fluctuated between 21 and 23 °C, and 2.5 and 5 °C respectively. January temperatures show a period of intermediate values (∼3.5 °C) between 4 ka and 5.5 ka. The reconstructed climate values generally match palaeolimnological data from the same core, which show five intervals of low lake level during the Holocene. They are also consistent with regional-scale COHMAP simulated palaeoclimate that shows contrasting patterns of rainfall variation between the northwesternmost part of Africa and the intertropical band. Received: 7 July 1997 / Accepted: 28 May 1998     

Classification of the stratified atmospheric boundary layers at Molve (Croatia) based on the similarity theory

Summary The stability parameter μ is suggested as the one which is determinable with satisfying accuracy for routine application by means of commonly accessible meteorological data at the Molve location (Croatia). The similarity functions applied for vertical wind speed simulation in the planetary boundary layer (PBL) at Molve were useful for the determination of local stability classes. Universal similarity functions were applied for unstable and neutral stability, whereas local similarity functions were established for stable stratification. Wind speed simulations were performed using two types of wind models. The Monin-Obukhov similarity theory was included in both types. However, it turned out that for the operative determination of the stability of the 35 m deep lowest layer, the stability parameter μ was locally a better stability parameter than the Monin-Obukhov parameter z/L. That was possibly because 35 m deep lowest layer sometimes (depending upon stability) includes a large proportion of the Ekman layer and parameter μ is originally designed for the deeper part of PBL than z/L that is originally designed for the surface layer. At Molve, the input data for local wind models as well as for the stability parameter μ were wind speed at 35 m and temperature at 2 and 35 m above the ground.     

Statistical-dynamical analysis of wind observations usinga variational approach

Summary  In this paper, we present a variational analysis of wind and pressure data which takes into account statistical characteristics of the data and a linear model describing the dynamical relations between the analysed variables. The variational approach is used to adjust the given data to the dynamical model and it will be shown that this adjustment process can be controlled in an uncomplicated, comprehensible, and reproducible way by very few parameters only. The dependence of the analysis results on these parameters is investigated in theory and this theoretical conclusions are tested by an application. It is studied how effectively this method can be used to correct erroneous data. We apply the analysis to the COADS (Comprehensive Ocean-Atmosphere Data Set) wind and pressure data of the Januaries 1951–1993 over the North Atlantic and introduce the Ekman balance as a weak dynamical constraint. This specific data set is chosen because several previous investigations suggested that there is a spurious trend in the COADS wind speed of ∼ 1 ms⁻¹ since the mid-century. The results show that the control parameters can be effectively used to shift the wind field continuously between the identity with the input data and the exact consistency with the dynamical model. But it has to be admitted that the Ekman balance is inadequate in the tropics and that it overestimates the magnitude of the horizontal vector wind although this dynamical model is more suitable than the frequently used geostrophic balance. Received January 20, 2000 Revised April 10, 2000     

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.