Observed and simulated multidecadal variability in the Northern Hemisphere

 Analyses of proxy based reconstructions of surface temperatures during the past 330 years show the existence of a distinct oscillatory mode of variability with an approximate time scale of 70 years. This variability is also seen in instrumental records, although the oscillatory nature of the variability is difficult to assess due to the short length of the instrumental record. The spatial pattern of this variability is hemispheric or perhaps even global in scale, but with particular emphasis on the Atlantic region. Independent analyses of multicentury integrations of two versions of the GFDL coupled atmosphere-ocean model also show the existence of distinct multidecadal variability in the North Atlantic region which resembles the observed pattern. The model variability involves fluctuations in the intensity of the thermohaline circulation in the North Atlantic. It is our intent here to provide a direct comparison of the observed variability to that simulated in a coupled ocean-atmosphere model, making use of both existing instrumental analyses and newly available proxy based multi-century surface temperature estimates. The analyses demonstrate a substantial agreement between the simulated and observed patterns of multidecadal variability in sea surface temperature (SST) over the North Atlantic. There is much less agreement between the model and observations for sea level pressure. Seasonal analyses of the variability demonstrate that for both the model and observations SST appears to be the primary carrier of the multidecadal signal. Received: 8 June 1999 / Accepted: 11 February 2000     

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     

Multidecadal variability of summer temperature over Romania and its relation with Atlantic Multidecadal Oscillation

We investigate the multidecadal variability of summer temperature over Romania as measured at 14 meteorological stations with long-term observational records. The dominant pattern of summer temperature variability has a monopolar structure and shows pronounced multidecadal variations. A correlation analysis reveals that these multidecadal variations are related with multidecadal variations in the frequency of four daily atmospheric circulation patterns from the North Atlantic region. It is found that on multidecadal time scales, negative summer mean temperature (TT) anomalies are associated with positive sea level pressure (SLP) anomalies centered over the northern part of the Atlantic Ocean and Scandinavia and negative SLP anomalies centered over the northern part of Africa. It is speculated that a possible cause of multidecadal fluctuations in the frequency of these four patterns are the sea surface temperature (SST) anomalies associated to the Atlantic Multidecadal Oscillation (AMO). These results have implications for predicting the evolution of summer temperature over Romania on multidecadal time scales.     

On the time-varying trend in global-mean surface temperature

The Earth has warmed at an unprecedented pace in the decades of the 1980s and 1990s (IPCC in Climate change 2007: the scientific basis, Cambridge University Press, Cambridge, 2007). In Wu et al. (Proc Natl Acad Sci USA 104:14889–14894, 2007) we showed that the rapidity of the warming in the late twentieth century was a result of concurrence of a secular warming trend and the warming phase of a multidecadal (~65-year period) oscillatory variation and we estimated the contribution of the former to be about 0.08°C per decade since ~1980. Here we demonstrate the robustness of those results and discuss their physical links, considering in particular the shape of the secular trend and the spatial patterns associated with the secular trend and the multidecadal variability. The shape of the secular trend and rather globally-uniform spatial pattern associated with it are both suggestive of a response to the buildup of well-mixed greenhouse gases. In contrast, the multidecadal variability tends to be concentrated over the extratropical Northern Hemisphere and particularly over the North Atlantic, suggestive of a possible link to low frequency variations in the strength of the thermohaline circulation. Depending upon the assumed importance of the contributions of ocean dynamics and the time-varying aerosol emissions to the observed trends in global-mean surface temperature, we estimate that up to one third of the late twentieth century warming could have been a consequence of natural variability.     

Seawater density variations in the North Atlantic and the Atlantic meridional overturning circulation

Seawater property changes in the North Atlantic Ocean affect the Atlantic meridional overturning circulation (AMOC), which transports warm water northward from the upper ocean and contributes to the temperate climate of Europe, as well as influences climate globally. Previous observational studies have focused on salinity and freshwater variability in the sinking region of the North Atlantic, since it is believed that a freshening North Atlantic basin can slow down or halt the flow of the AMOC. Here we use available data to show the importance of how density patterns over the upper ocean of the North Atlantic affect the strength of the AMOC. For the long-term trend, the upper ocean of the subpolar North Atlantic is becoming cooler and fresher, whereas the subtropical North Atlantic is becoming warmer and saltier. On a multidecadal timescale, the upper ocean of the North Atlantic has generally been warmer and saltier since 1995. The heat and salt content in the subpolar North Atlantic lags that in the subtropical North Atlantic by about 8–9 years, suggesting a lower latitude origin for the temperature and salinity anomalies. Because of the opposite effects of temperature and salinity on density for both long-term trend and multidecadal timescales, these variations do not result in a density reduction in the subpolar North Atlantic for slowing down the AMOC. Indeed, the variations in the meridional density gradient between the subpolar and subtropical North Atlantic Ocean suggest that the AMOC has become stronger over the past five decades. These observed results are supported by and consistent with some oceanic reanalysis products.     

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.     

Multidecadal variability in local growing season during 1901–2009

Global warming exerts a lengthening effect on the growing season, with observational evidences emerging from different regions over the world. However, the difficulty for a global overview of this effect for the last century arises from limited availability of the long-term daily observations. In this study, we find a good linear relationship between the start (end) date of local growing season (LGS) and the monthly mean temperature in April (October) using the global gridded daily temperature dataset for 1960–1999. Using homogenized daily temperature records from nine stations where the time series go back to the beginning of the twentieth century, we find that the rate of change in the start (end) date of the LGS for per degree warming in April (October) mean temperature keeps nearly constant throughout the time. This enables us to study LGS changes during the last century using global gridded monthly mean temperature data. The results show that during the period 1901–2009, averaged over the observation areas, the LGS length has increased by a rate of 0.89 days decade^?1, mainly due to an earlier start (?0.58 days decade^?1). This is smaller than those estimates for the late half of the twentieth century, because of multidecadal climate variability (MDV). A MDV component of the LGS index series is extracted by using Ensemble Empirical Mode Decomposition method. The MDV exhibits significant positive correlation with the Atlantic Multi–decadal Oscillation (AMO) over most of the Northern Hemisphere lands, but negative in parts of North America and Western Asia for start date of LGS. These are explained by analyzing differences in atmospheric circulation expressed by sea level pressure departures between the warm and cool phases of AMO. It is suggested that MDV in association with AMO accelerates the lengthening of LGS in Northern Hemisphere by 53 % for the period 1980–2009.     

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

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

Stochastically-forced multidecadal variability in the North Atlantic: a model study

Observations show a multidecadal signal in the North Atlantic ocean, but the underlying mechanism and cause of its timescale remain unknown. Previous studies have suggested that it may be driven by the North Atlantic Oscillation (NAO), which is the dominant pattern of winter atmospheric variability. To further address this issue, the global ocean general circulation model, Nucleus for European Modelling of the Ocean (NEMO), is driven using a 2,000 years long white noise forcing associated with the NAO. Focusing on key ocean circulation patterns, we show that the Atlantic Meridional Overturning Circulation (AMOC) and Sub-polar gyre (SPG) strength both have enhanced power at low frequencies but no dominant timescale, and thus provide no evidence for a oscillatory ocean-only mode of variability. Instead, both indices respond linearly to the NAO forcing, but with different response times. The variability of the AMOC at 30°N is strongly enhanced on timescales longer than 90 years, while that of the SPG strength starts increasing at 15 years. The different response characteristics are confirmed by constructing simple statistical models that show AMOC and SPG variability can be related to the NAO variability of the previous 53 and 10 winters, respectively. Alternatively, the AMOC and the SPG strength can be reconstructed with Auto-regressive (AR) models of order seven and five, respectively. Both statistical models reconstruct interannual and multidecadal AMOC variability well, while on the other hand, the AR(5) reconstruction of the SPG strength only captures multidecadal variability. Using these methods to reconstruct ocean variables can be useful for prediction and model intercomparision.     

Naturally forced multidecadal variability of the Atlantic meridional overturning circulation

The mechanisms by which natural forcing factors alone could drive simulated multidecadal variability in the Atlantic meridional overturning circulation (AMOC) are assessed in an ensemble of climate model simulations. It is shown for a new state-of-the-art general circulation model, HadGEM2-ES, that the most important of these natural forcings, in terms of the multidecadal response of the AMOC, is solar rather than volcanic forcing. AMOC strengthening occurs through a densification of the North Atlantic, driven by anomalous surface freshwater fluxes due to increased evaporation. These are related to persistent North Atlantic atmospheric circulation anomalies, driven by forced changes in the stratosphere, associated with anomalously weak solar irradiance during the late nineteenth and early twentieth centuries. Within a period of approximately 100 years the 11-year smoothed ensemble mean AMOC strengthens by 1.5 Sv and subsequently weakens by 1.9 Sv, representing respectively approximately 3 and 4 standard deviations of the 11-year smoothed control simulation. The solar-induced variability of the AMOC has various relevant climate impacts, such as a northward shift of the intertropical convergence zone, anomalous Amazonian rainfall, and a sustained increase in European temperatures. While this model has only a partial representation of the atmospheric response to solar variability, these results demonstrate the potential for solar variability to have a multidecadal impact on North Atlantic climate.     

Multidecadal-to-centennial SST variability in the MPI-ESM simulation ensemble for the last millennium

We assess the responses of North Atlantic, North Pacific, and tropical Indian Ocean Sea Surface Temperatures (SSTs) to natural forcing and their linkage to simulated global surface temperature (GST) variability in the MPI-Earth System Model simulation ensemble for the last millennium. In the simulations, North Atlantic and tropical Indian Ocean SSTs show a strong sensitivity to external forcing and a strong connection to GST. The leading mode of extra-tropical North Pacific SSTs is, on the other hand, rather resilient to natural external perturbations. Strong tropical volcanic eruptions and, to a lesser extent, variability in solar activity emerge as potentially relevant sources for multidecadal SST modes’ phase modulations, possibly through induced changes in the atmospheric teleconnection between North Atlantic and North Pacific that can persist over decadal and multidecadal timescales. Linkages among low-frequency regional modes of SST variability, and among them and GST, can remarkably vary over the integration time. No coherent or constant phasing is found between North Pacific and North Atlantic SST modes over time and among the ensemble members. Based on our assessments of how multidecadal transitions in simulated North Atlantic SSTs compare to reconstructions and of how they contribute characterizing simulated multidecadal regional climate anomalies, past regional climate multidecadal fluctuations seem to be reproducible as simulated ensemble-mean responses only for temporal intervals dominated by major external forcings.     

Non-uniform interhemispheric temperature trends over the past 550 years

The warming trend over the last century in the northern hemisphere (NH) was interrupted by cooling from ad 1940 to 1975, a period during which the southern hemisphere experienced pronounced warming. The cause of these departures from steady warming at multidecadal timescales are unclear; the prevailing explanation is that they are driven by non-uniformity in external forcings but recent models suggest internal climate drivers may play a key role. Paleoclimate datasets can help provide a long-term perspective. Here we use tree-rings to reconstruct New Zealand mean annual temperature over the last 550 years and demonstrate that this has frequently cycled out-of-phase with NH mean annual temperature at a periodicity of around 30–60 years. Hence, observed multidecadal fluctuations around the recent warming trend have precedents in the past, strongly implicating natural climate variation as their cause. We consider the implications of these changes in understanding and modelling future climate change.     

Bi-decadal variability excited in the coupled ocean–atmosphere system by strong tropical volcanic eruptions

Decadal and bi-decadal climate responses to tropical strong volcanic eruptions (SVEs) are inspected in an ensemble simulation covering the last millennium based on the Max Planck Institute—Earth system model. An unprecedentedly large collection of pre-industrial SVEs (up to 45) producing a peak annual-average top-of-atmosphere radiative perturbation larger than ?1.5 Wm^?2 is investigated by composite analysis. Post-eruption oceanic and atmospheric anomalies coherently describe a fluctuation in the coupled ocean–atmosphere system with an average length of 20–25 years. The study provides a new physically consistent theoretical framework to interpret decadal Northern Hemisphere (NH) regional winter climates variability during the last millennium. The fluctuation particularly involves interactions between the Atlantic meridional overturning circulation and the North Atlantic gyre circulation closely linked to the state of the winter North Atlantic Oscillation. It is characterized by major distinctive details. Among them, the most prominent are: (a) a strong signal amplification in the Arctic region which allows for a sustained strengthened teleconnection between the North Pacific and the North Atlantic during the first post-eruption decade and which entails important implications from oceanic heat transport and from post-eruption sea ice dynamics, and (b) an anomalous surface winter warming emerging over the Scandinavian/Western Russian region around 10–12 years after a major eruption. The simulated long-term climate response to SVEs depends, to some extent, on background conditions. Consequently, ensemble simulations spanning different phases of background multidecadal and longer climate variability are necessary to constrain the range of possible post-eruption decadal evolution of NH regional winter climates.     

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

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

Atmospheric response to the North Atlantic Ocean variability on seasonal to decadal time scales

The NCEP twentieth century reanalyis and a 500-year control simulation with the IPSL-CM5 climate model are used to assess the influence of ocean-atmosphere coupling in the North Atlantic region at seasonal to decadal time scales. At the seasonal scale, the air-sea interaction patterns are similar in the model and observations. In both, a statistically significant summer sea surface temperature (SST) anomaly with a horseshoe shape leads an atmospheric signal that resembles the North Atlantic Oscillation (NAO) during the winter. The air-sea interactions in the model thus seem realistic, although the amplitude of the atmospheric signal is half that observed, and it is detected throughout the cold season, while it is significant only in late fall and early winter in the observations. In both model and observations, the North Atlantic horseshoe SST anomaly pattern is in part generated by the spring and summer internal atmospheric variability. In the model, the influence of the ocean dynamics can be assessed and is found to contribute to the SST anomaly, in particular at the decadal scale. Indeed, the North Atlantic SST anomalies that follow an intensification of the Atlantic meridional overturning circulation (AMOC) by about 9 years, or an intensification of a clockwise intergyre gyre in the Atlantic Ocean by 6 years, resemble the horseshoe pattern, and are also similar to the model Atlantic Multidecadal Oscillation (AMO). As the AMOC is shown to have a significant impact on the winter NAO, most strongly when it leads by 9 years, the decadal interactions in the model are consistent with the seasonal analysis. In the observations, there is also a strong correlation between the AMO and the SST horseshoe pattern that influences the NAO. The analogy with the coupled model suggests that the natural variability of the AMOC and the gyre circulation might influence the climate of the North Atlantic region at the decadal scale.     

Leading patterns of the tropical Atlantic variability in a coupled general circulation model

This paper examines the mean annual cycle, interannual variability, and leading patterns of the tropical Atlantic Ocean simulated in a long-term integration of the climate forecast system (CFS), a state-of-the-art coupled general circulation model presently used for operational climate prediction at the National Centers for Environmental Prediction. By comparing the CFS simulation with corresponding observation-based analyses or reanalyses, it is shown that the CFS captures the seasonal mean climate, including the zonal gradients of sea surface temperature (SST) in the equatorial Atlantic Ocean, even though the CFS produces warm mean biases and underestimates the variability over the southeastern ocean. The seasonal transition from warm to cold phase along the equator is delayed 1 month in the CFS compared with the observations. This delay might be related to the failure of the model to simulate the cross-equatorial meridional wind associated with the African monsoon. The CFS also realistically simulates both the spatial structure and spectral distributions of the three major leading patterns of the SST anomalies in the tropical Atlantic Ocean: the south tropical Atlantic pattern (STA), the North tropical Atlantic pattern (NTA), and the southern subtropical Atlantic pattern (SSA). The CFS simulates the seasonal dependence of these patterns and partially reproduces their association with the El Niño-Southern Oscillation. The dynamical and thermodynamical processes associated with these patterns in the simulation and the observations are similar. The air-sea interaction processes associated with the STA pattern are well simulated in the CFS. The primary feature of the anomalous circulation in the Northern Hemisphere (NH) associated with the NTA pattern resembles that in the Southern Hemisphere (SH) linked with the SSA pattern, implying a similarity of the mechanisms in the evolution of these patterns and their connection with the tropical and extratropical anomalies in their respective hemispheres. The anomalies associated with both the SSA and NTA patterns are dominated by atmospheric fluctuations of equivalent-barotropic structure in the extratropics including zonally symmetric and asymmetric components. The zonally symmetric variability is associated with the annular modes, the Arctic Oscillation in the NH and the Antarctic Oscillation in the SH. The zonally asymmetric part of the anomalies in the Atlantic is teleconnected with the anomalies over the tropical Pacific. The misplaced teleconnection center over the southern subtropical ocean may be one of the reasons for the deformation of the SSA pattern in the CFS.     

Interdecadal variability of the meridional overturning circulation as an ocean internal mode

The meridional overturning circulation (MOC) in the coupled ECHAM5/MPIOM exhibits variability at periods of near 30 years and near 60 years. The 30-year variability, referred to as interdecadal variability (IDV), exist in an ocean model driven by climatological atmospheric forcing, suggesting that it is maintained by ocean dynamics; the 60-year variability, the multidecadal variability (MDV), is only observed in the fully coupled model and therefore is interpreted as an atmosphere–ocean coupled mode. The coexistence of the 30-year IDV and the 60-year MDV provides a possible explanation for the widespread time scales observed in climate variables. Further analyses of the climatologically forced ocean model shows that, the IDV is related to the interplay between the horizontal temperature-dominated density gradients and the ocean circulation: temperature anomalies move along the cyclonic subpolar gyre leading to fluctuations in horizontal density gradients and the subsequent weakening and strengthening of the MOC. This result is consistent with that from less complex models, indicating the robustness of the IDV. We further show that, along the North Atlantic Current path, the sea surface temperature anomalies are determined by the slow LSW advection at the intermediate depth.     

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.     

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

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