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.     

Decadal climate variability simulated in a coupled general circulation model

A 1000 year integration of the CSIRO coupled ocean-atmosphere general circulation model is used to study low frequency (decadal to centennial) climate variability in precipitation and temperature. The model is shown to exhibit sizeable decadal variability for these fields, generally accounting for approximately 20 to 40% of the variability (greater than one year) in precipitation and up to 80% for temperature. An empirical orthogonal function (EOF) analysis is applied to the model output to show some of the major statistical modes of low frequency variability. The first EOF spatial pattern looks very much like that of the interannual ENSO pattern. It bears considerable resemblance to observational estimates and is centred in the Pacific extending into both hemispheres. It modulates both precipitation and temperature globally. The EOF has a time evolution that appears to be more than just red noise. Finally, the link between SST in the Pacific with Australian rainfall variability seen in observations is also evident in the model. Received: 29 August 1998 / Accepted: 31 July 1999     

Stratospheric climate and variability from a general circulation model and observations

The climate and natural variability of the large-scale stratospheric circulation simulated by a newly developed general circulation model are evaluated against available global observations. The simulation consisted of a 30-year annual cycle integration performed with a comprehensive model of the troposphere and stratosphere. The observations consisted of a 15-year dataset from global operational analyses of the troposphere and stratosphere. The model evaluation concentrates on the simulation of the evolution of the extratropical stratospheric circulation in both hemispheres. The December–February climatology of the observed zonal mean winter circulation is found to be reasonably well captured by the model, although in the Northern Hemisphere upper stratosphere the simulated westerly winds are systematically stronger and a cold bias is apparent in the polar stratosphere. This Northern Hemisphere stratospheric cold bias virtually disappears during spring (March–May), consistent with a realistic simulation of the spring weakening of the mean westerly winds in the model. A considerable amount of monthly interannual variability is also found in the simulation in the Northern Hemisphere in late winter and early spring. The simulated interannual variability is predominantly caused by polar warmings of the stratosphere, in agreement with observations. The breakdown of the Northern Hemisphere stratospheric polar vortex appears therefore to occur in a realistic way in the model. However, in early winter the model severely underestimates the interannual variability, especially in the upper troposphere. The Southern Hemisphere winter (June–August) zonal mean temperature is systematically colder in the model, and the simulated winds are somewhat too strong in the upper stratosphere. Contrary to the results for the Northern Hemisphere spring, this model cold bias worsens during the Southern Hemisphere spring (September–November). Significant discrepancies between the model results and the observations are therefore found during the breakdown of the Southern Hemisphere polar vortex. For instance, the simulated Southern Hemisphere stratosphere westerly jet continuously decreases in intensity more or less in situ from June to November, while the observed stratospheric jet moves downward and poleward.This paper was presented at the Third International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 4–8 Sept. 1995 under the auspice of the Max Planck Institute for Meteorology, Hamburg. Editor for these papers is L. Dümenil.     

Modes of climate variability as simulated by a coupled general circulation model. Part I: ENSO-like climate variability and its low-frequency modulation

The El Niño-Southern Oscillation (ENSO) is investigated in a multicentury integration conducted with the coupled general circulation model (CGCM) ECHAM3/LSG. The quasiperiodic interannual oscillations of the simulated equatorial Pacific climate system are due to subsurface temperature anomaly propagation and a positive atmosphere-ocean feedback. The gravest internal wave modes contribute to the generation of these anomalies. The simulated ENSO has a characteristic period of 5–8 years. Due to the coarse resolution of the ocean model the ENSO amplitude is underestimated by a factor of three as compared to observations. The model ENSO is associated with the typical atmospheric teleconnection patterns. Using wavelet statistics two characteristic interdecadal modulations of the ENSO variance are identified. The origins of a 22 and 35?y ENSO modulation as well as the characteristic ENSO response to greenhouse warming simulated by our model are discussed.     

An upper ocean general circulation model for climate studies: global simulation with seasonal cycle

The global ocean circulation with a seasonal cycle has been simulated with a two-and-a-half layer upper-ocean model. This model was developed for the purpose of coupling to an atmospheric general circulation model for climate studies on decadal time scales. The horizontal resolution is 4° latitude by 5° longitude and is thus not eddy-resolving. Effects of bottom topography are neglected. In the vertical, the model resolves the oceanic mixed layer and the thermocline. A thermodynamic sea-ice model is coupled to the mixed layer. The model is forced at the surface with seasonally varying (a) observed wind stress, (b) heat fluxes, as defined by an atmospheric equilibrium temperature, and (c) Newtonian-type surface salt fluxes. The second layer is coupled to the underlying deep ocean through Newtonian-type diffusive heat and salt fluxes, convective overturning, and mass entrainment in the upwelling regions of the subpolar gyres. The overall global distributions of mixed layer temperature, salinity and thickness are favorably reproduced. Inherent limitations due to coarse horizontal resolution result in large mixed-layer temperature errors near continental boundaries and in weak current systems. Sea ice distributions agree well with observations except in the interiors of the Ross and Weddell Seas. A realistic time rate of change of heat storage is simulated. There is also realistic heat transport from low to high latitudes.     

Assessing the role of North Atlantic freshwater forcing in millennial scale climate variability: a tropical Atlantic perspective

This study analyzes a three-member ensemble of experiments, in which 0.1 Sv of freshwater was applied to the North Atlantic for 100 years in order to address the potential for large freshwater inputs in the North Atlantic to drive abrupt climate change. The model used is the GFDL R30 coupled ocean–atmosphere general circulation model. We focus in particular on the effects of this forcing on the tropical Atlantic region, which has been studied extensively by paleoclimatologists. In response to the freshwater forcing, North Atlantic meridional overturning circulation is reduced to roughly 40% by the end of the 100 year freshwater pulse. Consequently, the North Atlantic region cools by up to 8°C. The extreme cooling of the North Atlantic increases the pole-to-equator temperature gradient and requires more heat be provided to the high latitude Atlantic from the tropical Atlantic. To accommodate the increased heat requirement, the ITCZ shifts southward to allow for greater heat transport across the equator. Accompanying this southward ITCZ shift, the Northeast trade winds strengthen and precipitation patterns throughout the tropical Atlantic are altered. Specifically, precipitation in Northeast Brazil increases, and precipitation in Africa decreases slightly. In addition, we find that surface air temperatures warm over the tropical Atlantic and over Africa, but cool over northern South America. Sea-surface temperatures in the tropical Atlantic warm slightly with larger warm anomalies developing in the thermocline. These responses are robust for each member of the ensemble, and have now been identified by a number of freshwater forcing studies using coupled OAGCMs. The model responses to freshwater forcing are generally smaller in magnitude, but have the same direction, as paleoclimate data from the Younger Dryas suggest. In certain cases, however, the model responses and the paleoclimate data directly contradict one another. Discrepancies between the model simulations and the paleoclimate data could be due to a number of factors, including inaccuracies in the freshwater forcing, inappropriate boundary conditions, and uncertainties in the interpretation of the paleoclimate data. Despite these discrepancies, it is clear from our results that abrupt climate changes in the high latitude North Atlantic have the potential to significantly impact tropical climate. This warrants further model experimentation into the role of freshwater forcing in driving climate change.     

An approximate energy cycle for inter-member variability in ensemble simulations of a regional climate model

The presence of internal variability (IV) in ensembles of nested regional climate model (RCM) simulations is now widely acknowledged in the community working on dynamical downscaling. IV is defined as the inter-member spread between members in an ensemble of simulations performed by a given RCM driven by identical lateral boundary conditions (LBC), where different members are being initialised at different times. The physical mechanisms responsible for the time variations and structure of such IV have only recently begun to receive attention. Recent studies have shown empirical evidence of a close parallel between the energy conversions associated with the time fluctuations of IV in ensemble simulations of RCM and the energy conversions taking place in weather systems. Inspired by the classical work on global energetics of weather systems, we sought a formulation of an energy cycle for IV that would be applicable for limited-area domain. We develop here a novel formalism based on local energetics that can be applied to further our understanding IV. Prognostic equations for ensemble-mean kinetic energy and available enthalpy are decomposed into contributions due to ensemble-mean variables (EM) and those due to deviations from the ensemble mean (IV). Together these equations constitute an energy cycle for IV in ensemble simulations of RCM. Although the energy cycle for IV was developed in a context entirely different from that of energetics of weather systems, the exchange terms between the various reservoirs have a rather similar mathematical form, which facilitates some interpretations of their physical meaning.     

The impact of combined ENSO and PDO on the PNA climate: a 1,000-year climate modeling study

This study analyzes the atmospheric response to the combined Pacific interannual ENSO and decadal–interdecadal PDO variability, with a focus on the Pacific-North American (PNA) sector, using a 1,000-year long integration of the Canadian Center for Climate Modelling and Analysis (CCCma) coupled climate model. Both the tropospheric circulation and the North American temperature suggest an enhanced PNA-like climate response and impacts on North America when ENSO and PDO variability are in phase. The anomalies of the centers of action for the PNA-like pattern are significantly different from zero and the anomaly pattern is field significant. In association with the stationary wave anomalies, large stationary wave activity fluxes appear in the mid-high latitudes originating from the North Pacific and flowing downstream toward North America. There are significant Rossby wave source anomalies in the extratropical North Pacific and in the subtropical North Pacific. In addition, the axis of the Pacific storm track shifts southward with the positive PNA. Atmospheric heating anomalies associated with ENSO variability are confined primarily to the tropics. There is an anomalous heating center over the northeast Pacific, together with anomalies with the same polarity in the tropical Pacific, for the PDO variability. The in-phase combination of ENSO and PDO would in turn provide anomalous atmospheric energy transports towards North America from both the Tropical Pacific and the North Pacific, which tends to favor the occurrence of stationary wave anomalies and would lead to a PNA-like wave anomaly structure. The modeling results also confirm our analysis based on the observational record in the twentieth century.     

Quantifying internal variability in a regional climate model: a case study for Southern Africa

This study aims at presenting various methodologies to separate the reproducible and irreproducible components of seasonal and intraseasonal climate variability simulated by a regional climate model over Southern Africa (south of 15°S), during an austral summer rainy season representative of the climatology. To that end, a 30-member ensemble simulation is performed using WRF laterally forced by the ERA40 reanalyses. Retained metrics include the analysis of weather regimes, signal-to-noise ratio, inter-member standard deviation and coefficient of variation. At the seasonal timescale, simulated rainfall amounts generally show a strong reproducibility, except in the subtropics and over the southern part of the Mozambique Channel. There, the number of rainy days is roughly similar in all members, while their average intensity varies extensively. At the intraseasonal timescale, the chronology of weather regimes, derived from the 500?hPa geopotential height, is highly reproducible. Rainfall variability is much less reproducible, especially in the central parts of the domain and near its outflow boundaries. Analysis of a South African regional index nonetheless indicates that both wet and dry spells tend to be accurately simulated and occur in phase in most members, demonstrating that they are embedded in large-scale variability patterns. Internal variability is lastly related to the lateral forcings along the domain boundaries. An objective classification of inflow/outflow mass fluxes allows identification of the recurrent synoptic configurations that favor strong or weak regional reproducibility. The main uncertainties concern the basic morphological features of rain-bearing systems (i.e., their spatial extension, location and propagation speed). Consequences for tropical-temperate interactions are then discussed.     

The impact of resolution on the adjustment and decadal variability of the Atlantic meridional overturning circulation in a coupled climate model

Variations in the Atlantic meridional overturning circulation (MOC) exert an important influence on climate, particularly on decadal time scales. Simulation of the MOC in coupled climate models is compromised, to a degree that is unknown, by their lack of fidelity in resolving some of the key processes involved. There is an overarching need to increase the resolution and fidelity of climate models, but also to assess how increases in resolution influence the simulation of key phenomena such as the MOC. In this study we investigate the impact of significantly increasing the (ocean and atmosphere) resolution of a coupled climate model on the simulation of MOC variability by comparing high and low resolution versions of the same model. In both versions, decadal variability of the MOC is closely linked to density anomalies that propagate from the Labrador Sea southward along the deep western boundary. We demonstrate that the MOC adjustment proceeds more rapidly in the higher resolution model due the increased speed of western boundary waves. However, the response of the Atlantic sea surface temperatures to MOC variations is relatively robust—in pattern if not in magnitude—across the two resolutions. The MOC also excites a coupled ocean-atmosphere response in the tropical Atlantic in both model versions. In the higher resolution model, but not the lower resolution model, there is evidence of a significant response in the extratropical atmosphere over the North Atlantic 6?years after a maximum in the MOC. In both models there is evidence of a weak negative feedback on deep density anomalies in the Labrador Sea, and hence on the MOC (with a time scale of approximately ten years). Our results highlight the need for further work to understand the decadal variability of the MOC and its simulation in climate models.     

Future climate of the Caribbean from a super-high-resolution atmospheric general circulation model

Present-day (1979–2003) and future (2075–2099) simulations of mean and extreme rainfall and temperature are examined using data from the Meteorological Research Institute super-high-resolution atmospheric general circulation model. Analyses are performed over the 20-km model grid for (1) a main Caribbean basin, (2) sub-regional zones, and (3) specific Caribbean islands. Though the model’s topography underestimates heights over the eastern Caribbean, it captures well the present-day spatial and temporal variations of seasonal and annual climates. Temperature underestimations range from 0.1 °C to 2 °C with respect to the Japanese Reanalysis and the Climatic Research Unit datasets. The model also captures fairly well sub-regional scale variations in the rainfall climatology. End-of-century projections under the Intergovernmental Panel on Climate Change SRES A1B scenario indicate declines in rainfall amounts by 10–20 % for most of the Caribbean during the early (May–July) and late (August–October) rainy seasons relative to the 1979–2003 baselines. The early dry season (November–January) is also projected to get wetter in the far north and south Caribbean by approximately 10 %. The model also projects a warming of 2–3 °C over the Caribbean region. Analysis of future climate extremes indicate a 5–10 % decrease in the simple daily precipitation intensity but no significant change in the number of consecutive dry days for Cuba, Jamaica, southern Bahamas, and Haiti. There is also indication that the number of hot days and nights will significantly increase over the main Caribbean basin.     

Storm tracks in a warmer climate: sensitivity studies with a simplified global circulation model

 A simplified global circulation model is used to analyse a greenhouse warming experiment simulated by a comprehensive general circulation model. The given GCM scenario and control climates are assimilated by the simplified model using a dynamical relaxation technique. Two sets of sensitivity experiments investigate the influence of upper and lower tropospheric changes in baroclinicity on the Northern Hemisphere winter storm tracks. The results show that the three-dimensional structure of both the background flow and the changes in baroclinicity are important for the behaviour of mid-latitude eddy activity in relation to modifications of the baroclinicity. In general, the mid-latitude eddy activity is more sensitive to lower than to upper level changes in baroclinicity. The results further suggest that the simulated storm track changes in the GCM scenario are dominated by local modes of baroclinic instability. Received: 17 December 1996 / Accepted: 14 May 1998     

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.     

Decadal-timescale changes of the Atlantic overturning circulation and climate in a coupled climate model with a hybrid-coordinate ocean component

A wide range of statistical tools is used to investigate the decadal variability of the Atlantic Meridional Overturning Circulation (AMOC) and associated key variables in a climate model (CHIME, Coupled Hadley-Isopycnic Model Experiment), which features a novel ocean component. CHIME is as similar as possible to the 3rd Hadley Centre Coupled Model (HadCM3) with the important exception that its ocean component is based on a hybrid vertical coordinate. Power spectral analysis reveals enhanced AMOC variability for periods in the range 15–30 years. Strong AMOC conditions are associated with: (1) a Sea Surface Temperature (SST) anomaly pattern reminiscent of the Atlantic Multi-decadal Oscillation (AMO) response, but associated with variations in a northern tropical-subtropical gradient; (2) a Surface Air Temperature anomaly pattern closely linked to SST; (3) a positive North Atlantic Oscillation (NAO)-like pattern; (4) a northward shift of the Intertropical Convergence Zone. The primary mode of AMOC variability is associated with decadal changes in the Labrador Sea and the Greenland Iceland Norwegian (GIN) Seas, in both cases linked to the tropical activity about 15 years earlier. These decadal changes are controlled by the low-frequency NAO that may be associated with a rapid atmospheric teleconnection from the tropics to the extratropics. Poleward advection of salinity anomalies in the mixed layer also leads to AMOC changes that are linked to processes in the Labrador Sea. A secondary mode of AMOC variability is associated with interannual changes in the Labrador and GIN Seas, through the impact of the NAO on local surface density.     

Simulated climate change during the last 1,000 years: comparing the ECHO-G general circulation model with the MAGICC simple climate model

An intercomparison of eight climate simulations, each driven with estimated natural and anthropogenic forcings for the last millennium, indicates that the so-called “Erik” simulation of the ECHO-G coupled ocean-atmosphere climate model exhibits atypical behaviour. The ECHO-G simulation has a much stronger cooling trend from 1000 to 1700 and a higher rate of warming since 1800 than the other simulations, with the result that the overall amplitude of millennial-scale temperature variations in the ECHO-G simulation is much greater than in the other models. The MAGICC (Model for the Assessment of Greenhouse-gas-Induced Climate Change) simple climate model is used to investigate possible causes of this atypical behaviour. It is shown that disequilibrium in the initial conditions probably contributes spuriously to the cooling trend in the early centuries of the simulation, and that the omission of tropospheric sulphate aerosol forcing is the likely explanation for the anomalously large recent warming. The simple climate model results are used to adjust the ECHO-G Erik simulation to mitigate these effects, which brings the simulation into better agreement with the other seven models considered here and greatly reduces the overall range of temperature variations during the last millennium simulated by ECHO-G. Smaller inter-model differences remain which can probably be explained by a combination of the particular forcing histories and model sensitivities of each experiment. These have not been investigated here, though we have diagnosed the effective climate sensitivity of ECHO-G to be 2.39±0.11 K for a doubling of CO₂.     

Validation of synoptic circulation patterns simulated by the Canadian climate centre general circulation model for western north America: Research note

Abstract

In order to assess the ability of a GCM to simulate regional to synoptic scale atmospheric structures, a correlation‐based computer‐assisted gridded map typing procedure is used to compare daily pressure (MSL) and geopotential height fields (500 hPa) from a GCM simulation of the present climate to a decade of NMC analyses. The model is able to reproduce the entire range of synoptic circulation types. However, statistically significant differences in the seasonal frequencies and variances of the main circulation types are evident. These differences, which are most pronounced in the winter (at 500 hPa) and in spring and autumn in the MSL fields, are consistent with subtle errors in the predicted fields at the hemispheric scale. The lack of agreement between the NMC climatology and the “control” simulation precludes extension of this approach to investigation of climate change impacts in western north America, and to more meteorologically dynamic extra‐tropical regions. The map‐typing procedure is shown to be an appropriate GCM synoptic‐scale validation tool that permits direct comparison of GCM output and observed fields. As such, it has the potential to elucidate the regional‐scale impacts of global climatic change through established synoptic circulation environment relationships.     

Entropy production and coarse graining of the climate fields in a general circulation model

We extend the analysis of the thermodynamics of the climate system by investigating the role played by processes taking place at various spatial and temporal scales through a procedure of coarse graining. We show that the coarser is the graining of the climatic fields, the lower is the resulting estimate of the material entropy production. In other terms, all the spatial and temporal scales of variability of the thermodynamic fields provide a positive contribution to the material entropy production. This may be interpreted also as that, at all scales, the temperature fields and the heating fields resulting from the convergence of turbulent fluxes have a negative correlation, while the opposite holds between the temperature fields and the radiative heating fields. Moreover, we obtain that the latter correlations are stronger, which confirms that radiation acts as primary driver for the climatic processes, while the material fluxes dampen the resulting fluctuations through dissipative processes. We also show, using specific coarse-graining procedures, how one can separate the various contributions to the material entropy production coming from the dissipation of kinetic energy, the vertical sensible and latent heat fluxes, and the large scale horizontal fluxes, without resorting to the full three-dimensional time dependent fields. We find that most of the entropy production is associated to irreversible exchanges occurring along the vertical direction, and that neglecting the horizontal and time variability of the fields has a relatively small impact on the estimate of the material entropy production. The approach presented here seems promising for testing climate models, for assessing the impact of changing their parametrizations and their resolution, as well as for investigating the atmosphere of exoplanets, because it allows for evaluating the error in the estimate of their thermodynamical properties due to the lack of high-resolution data. The findings on the impact of coarse graining on the thermodynamic fields on the estimate of the material entropy production deserve to be explored in a more general context, because they provide a way for understanding the relationship between forced fluctuations and dissipative processes in continuum systems.