The role of atmosphere feedbacks during ENSO in the CMIP3 models

Several studies using ocean-atmosphere GCMs suggest that the atmospheric component plays a dominant role in the modelled ENSO. To help elucidate these findings, the two main atmosphere feedbacks relevant to ENSO, the Bjerknes positive feedback (µ) and the heat flux negative feedback (), are analysed here in 12 coupled GCMs. We find that the models generally underestimate both feedbacks, leading to an error compensation. The strength of is inversely related to the ENSO amplitude in the models and the latent heat and shortwave flux components of this feedback dominate. Furthermore, the shortwave component could help explain the model diversity in both overall and ENSO amplitude.

North Atlantic subtropical mode waters and ocean memory in HadCM3

A study of the formation and propagation of volume anomalies in North Atlantic Mode Waters is presented, based on 100 yr of monthly mean fields taken from the control run of the Third Hadley Centre Coupled Ocean-Atmosphere GCM (HadCM3). Analysis of the temporal and. spatial variability in the thickness between pairs of isothermal surfaces bounding the central temperature of the three main North Atlantic subtropical mode waters shows that large-scale variability in formation occurs over time scales ranging from 5 to 20 yr. The largest formation anomalies are associated with a southward shift in the mixed layer isothermal distribution, possibly due to changes in the gyre dynamics and/or changes in the overlying wind field and air-sea heat fluxes. The persistence of these anomalies is shown to result from their subduction beneath the winter mixed layer base where they recirculate around the subtropical gyre in the background geostrophic flow. Anomalies in the warmest mode (18 degrees C) formed on the western side of the basin persist for up to 5 yr. They are removed by mixing transformation to warmer classes and are returned to the seasonal mixed layer near the Gulf Stream where the stored heat may be released to the atmosphere. Anomalies in the cooler modes (16 degrees and 14 degrees C) formed on the eastern side of the basin persist for up to 10 yr. There is no clear evidence of significant transformation of these cooler mode anomalies to adjacent classes. It has been proposed that the eastern anomalies are removed through a tropical-subtropical water mass exchange mechanism beneath the trade wind belt (south of 20 degrees N). The analysis shows that anomalous mode water formation plays a key role in the long-term storage of heat in the model, and that the release of heat associated with these anomalies suggests a predictable climate feedback mechanism.

Elimination of the Greenland Ice Sheet in a High CO2 Climate

Projections of future global sea level depend on reliable estimates of changes in the size of polar ice sheets. Calculating this directly from global general circulation models (GCMs) is unreliable because the coarse resolution of 100 km or more is unable to capture narrow ablation zones, and ice dynamics is not usually taken into account in GCMs. To overcome these problems a high-resolution (20 km) dynamic ice sheet model has been coupled to the third Hadley Centre Coupled Ocean-Atmosphere GCM (HadCM3). A novel feature is the use of two-way coupling, so that climate changes in the GCM drive ice mass changes in the ice sheet model that, in turn, can alter the future climate through changes in orography, surface albedo, and freshwater input to the model ocean. At the start of the main experiment the atmospheric carbon dioxide concentration was increased to 4 times the preindustrial level and held constant for 3000 yr. By the end of this period the Greenland ice sheet is almost completely ablated and has made a direct contribution of approximately 7 m to global average sea level, causing a peak rate of sea level rise of 5 mm yr-1 early in the simulation. The effect of ice sheet depletion on global and regional climate has been examined and it was found that apart from the sea level rise, the long-term effect on global climate is small. However, there are some significant regional climate changes that appear to have reduced the rate at which the ice sheet ablates.

The physical properties of the atmosphere in the new Hadley Centre Global Environmental Model (HadGEM1). Part II: Aspects of variability and regional climate

The performance of the atmospheric component of the new Hadley Centre Global Environmental Model (HadGEM1) is assessed in terms of its ability to represent a selection of key aspects of variability in the Tropics and extratropics. These include midlatitude storm tracks and blocking activity, synoptic variability over Europe, and the North Atlantic Oscillation together with tropical convection, the Madden-Julian oscillation, and the Asian summer monsoon. Comparisons with the previous model, the Third Hadley Centre Coupled Ocean-Atmosphere GCM (HadCM3), demonstrate that there has been a considerable increase in the transient eddy kinetic energy (EKE), bringing HadGEM1 into closer agreement with current reanalyses. This increase in EKE results from the increased horizontal resolution and, in combination with the improved physical parameterizations, leads to improvements in the representation of Northern Hemisphere storm tracks and blocking. The simulation of synoptic weather regimes over Europe is also greatly improved compared to HadCM3, again due to both increased resolution and other model developments. The variability of convection in the equatorial region is generally stronger and closer to observations than in HadCM3. There is, however, still limited convective variance coincident with several of the observed equatorial wave modes. Simulation of the Madden-Julian oscillation is improved in HadGEM1: both the activity and interannual variability are increased and the eastward propagation, although slower than observed, is much better simulated. While some aspects of the climatology of the Asian summer monsoon are improved in HadGEM1, the upper-level winds are too weak and the simulation of precipitation deteriorates. The dominant modes of monsoon interannual variability are similar in the two models, although in HadCM3 this is linked to SST forcing, while in HadGEM1 internal variability dominates. Overall, analysis of the phenomena considered here indicates that HadGEM1 performs well and, in many important respects, improves upon HadCM3. Together with the improved representation of the mean climate, this improvement in the simulation of atmospheric variability suggests that HadGEM1 provides a sound basis for future studies of climate and climate change.

Granger causality of coupled climate processes: Ocean feedback on the North Atlantic oscillation

This study uses a Granger causality time series modeling approach to quantitatively diagnose the feedback of daily sea surface temperatures (SSTs) on daily values of the North Atlantic Oscillation (NAO) as simulated by a realistic coupled general circulation model (GCM). Bivariate vector autoregressive time series models are carefully fitted to daily wintertime SST and NAO time series produced by a 50-yr simulation of the Third Hadley Centre Coupled Ocean-Atmosphere GCM (HadCM3). The approach demonstrates that there is a small yet statistically significant feedback of SSTs oil the NAO. The SST tripole index is found to provide additional predictive information for the NAO than that available by using only past values of NAO-the SST tripole is Granger causal for the NAO. Careful examination of local SSTs reveals that much of this effect is due to the effect of SSTs in the region of the Gulf Steam, especially south of Cape Hatteras. The effect of SSTs on NAO is responsible for the slower-than-exponential decay in lag-autocorrelations of NAO notable at lags longer than 10 days. The persistence induced in daily NAO by SSTs causes long-term means of NAO to have more variance than expected from averaging NAO noise if there is no feedback of the ocean on the atmosphere. There are greater long-term trends in NAO than can be expected from aggregating just short-term atmospheric noise, and NAO is potentially predictable provided that future SSTs are known. For example, there is about 10%-30% more variance in seasonal wintertime means of NAO and almost 70% more variance in annual means of NAO due to SST effects than one would expect if NAO were a purely atmospheric process.

El Nino-mean state-seasonal cycle interactions in a multi-model ensemble

The modelled El Nino-mean state-seasonal cycle interactions in 23 coupled ocean-atmosphere GCMs, including the recent IPCC AR4 models, are assessed and compared to observations and theory. The models show a clear improvement over previous generations in simulating the tropical Pacific climatology. Systematic biases still include too strong mean and seasonal cycle of trade winds. El Nino amplitude is shown to be an inverse function of the mean trade winds in agreement with the observed shift of 1976 and with theoretical studies. El Nino amplitude is further shown to be an inverse function of the relative strength of the seasonal cycle. When most of the energy is within the seasonal cycle, little is left for inter-annual signals and vice versa. An interannual coupling strength (ICS) is defined and its relation with the modelled El Nino frequency is compared to that predicted by theoretical models. An assessment of the modelled El Nino in term of SST mode (S-mode) or thermocline mode (T-mode) shows that most models are locked into a S-mode and that only a few models exhibit a hybrid mode, like in observations. It is concluded that several basic El Nino-mean state-seasonal cycle relationships proposed by either theory or analysis of observations seem to be reproduced by CGCMs. This is especially true for the amplitude of El Nino and is less clear for its frequency. Most of these relationships, first established for the pre-industrial control simulations, hold for the double and quadruple CO2 stabilized scenarios. The models that exhibit the largest El Nino amplitude change in these greenhouse gas (GHG) increase scenarios are those that exhibit a mode change towards a T-mode (either from S-mode to hybrid or hybrid to T-mode). This follows the observed 1976 climate shift in the tropical Pacific, and supports the-still debated-finding of studies that associated this shift to increased GHGs. In many respects, these models are also among those that best simulate the tropical Pacific climatology (ECHAM5/MPI-OM, GFDL-CM2.0, GFDL-CM2.1, MRI-CGM2.3.2, UKMO-HadCM3). Results from this large subset of models suggest the likelihood of increased El Nino amplitude in a warmer climate, though there is considerable spread of El Nino behaviour among the models and the changes in the subsurface thermocline properties that may be important for El Nino change could not be assessed. There are no clear indications of an El Nino frequency change with increased GHG.

Multidecadal modulation of El Nino-Southern Oscillation (ENSO) variance by Atlantic Ocean sea surface temperatures

Observations suggest a possible link between the Atlantic Multidecadal Oscillation (AMO) and El Nino Southern Oscillation (ENSO) variability, with the warm AMO phase being related to weaker ENSO variability. A coupled ocean-atmosphere model is used to investigate this relationship and to elucidate mechanisms responsible for it. Anomalous sea surface temperatures (SSTs) associated with the positive AMO lead to change in the basic state in the tropical Pacific Ocean. This basic state change is associated with a deepened thermocline and reduced vertical stratification of the equatorial Pacific ocean, which in turn leads to weakened ENSO variability. We suggest a role for an atmospheric bridge that rapidly conveys the influence of the Atlantic Ocean to the tropical Pacific. The results suggest a non-local mechanism for changes in ENSO statistics and imply that anomalous Atlantic ocean SSTs can modulate both mean climate and climate variability over the Pacific.

Ocean circulation and climate in an idealised Pangean OAGCM

An idealised Pangean configuration is integrated in a coupled ocean atmosphere general circulation model to investigate the form of the ocean circulation and its impacts on the large scale climate system. A vigorous, hemispherically symmetric overturning is found, driven by deep water formation at high latitudes. Whilst the peak mass transport is around 100Sv, a low vertical temperature gradient in the ocean means that the maximum heat transport is only 1.2PW. The geographical change in the coupled model is found to produce a global average warming of 2°C, despite an increase in global surface albedo. This occurs through changes in the atmospheric water vapour and cloud distributions. There is also reduction in the equator-pole temperature gradient, largely attributable to the same causes, avoiding the paradox of low meridional temperature gradients without increased polar heat transport.

The impact of global freshwater forcing on the thermohaline circulation: adjustment of North Atlantic convection sites in a CGCM

On the time scale of a century, the Atlantic thermohaline circulation (THC) is sensitive to the global surface salinity distribution. The advection of salinity toward the deep convection sites of the North Atlantic is one of the driving mechanisms for the THC. There is both a northward and a southward contributions. The northward salinity advection (Nsa) is related to the evaporation in the subtropics, and contributes to increased salinity in the convection sites. The southward salinity advection (Ssa) is related to the Arctic freshwater forcing and tends on the contrary to diminish salinity in the convection sites. The THC changes results from a delicate balance between these opposing mechanisms. In this study we evaluate these two effects using the IPSL-CM4 ocean-atmosphere-sea-ice coupled model (used for IPCC AR4). Perturbation experiments have been integrated for 100 years under modern insolation and trace gases. River runoff and evaporation minus precipitation are successively set to zero for the ocean during the coupling procedure. This allows the effect of processes Nsa and Ssa to be estimated with their specific time scales. It is shown that the convection sites in the North Atlantic exhibit various sensitivities to these processes. The Labrador Sea exhibits a dominant sensitivity to local forcing and Ssa with a typical time scale of 10 years, whereas the Irminger Sea is mostly sensitive to Nsa with a 15 year time scale. The GIN Seas respond to both effects with a time scale of 10 years for Ssa and 20 years for Nsa. It is concluded that, in the IPSL-CM4, the global freshwater forcing damps the THC on centennial time scales.

The vertical structure of ocean heat transport

One of the most important contributions the ocean makes to Earth's climate is through its poleward heat transport: about 1.5 PW or more than 30% of that accomplished by the ocean-atmosphere system (Trenberth and Caron, 2001). Recently, concern has arisen over whether global warming could affect this heat transport (Watson et al., 2001), for example, reducing high latitude convection and triggering a collapse of the deep overturning circulation (Rahmstorf, 1995). While the consequences of abrupt changes in oceanic circulation should be of concern, we argue that the attention devoted to deep circulations is disproportionate to their role in heat transport. For this purpose, we introduce a heat function which identifies the contribution to the heat transport by different components of the oceanic circulation. A new view of the ocean emerges in which a shallow surface intensified circulation dominates the poleward heat transport.

Simulated changes in the relationship between tropical ocean temperatures and the western African monsoon during the mid-Holocene

Results from nine coupled ocean-atmosphere simulations have been used to investigate changes in the relationship between the variability of monsoon precipitation over western Africa and tropical sea surface temperatures (SSTs) between the mid-Holocene and the present day. Although the influence of tropical SSTs on the African monsoon is generally overestimated in the control simulations, the models reproduce aspects of the observed modes of variability. Thus, most models reproduce the observed negative correlation between western Sahelian precipitation and SST anomalies in the eastern tropical Pacific, and many of them capture the positive correlation between SST anomalies in the eastern tropical Atlantic and precipitation over the Guinea coastal region. Although the response of individual model to the change in orbital forcing between 6 ka and present differs somewhat, eight of the models show that the strength of the teleconnection between SSTs in the eastern tropical Pacific and Sahelian precipitation is weaker in the mid-Holocene. Some of the models imply that this weakening was associated with a shift towards longer time periods (from 3–5 years in the control simulations toward 4–10 years in the mid-Holocene simulations). The simulated reduction in the teleconnection between eastern tropical Pacific SSTs and Sahelian precipitation appears to be primarily related to a reduction in the atmospheric circulation bridge between the Pacific and West Africa but, depending on the model, other mechanisms such as increased importance of other modes of tropical ocean variability or increased local recycling of monsoonal precipitation can also play a role.

A simple method for integrating a complex model into an ensemble data assimilation system using MPI

This paper details a strategy for modifying the source code of a complex model so that the model may be used in a data assimilation context, {and gives the standards for implementing a data assimilation code to use such a model}. The strategy relies on keeping the model separate from any data assimilation code, and coupling the two through the use of Message Passing Interface (MPI) {functionality}. This strategy limits the changes necessary to the model and as such is rapid to program, at the expense of ultimate performance. The implementation technique is applied in different models with state dimension up to $2.7 \times 10^8$. The overheads added by using this implementation strategy in a coupled ocean-atmosphere climate model are shown to be an order of magnitude smaller than the addition of correlated stochastic random errors necessary for some nonlinear data assimilation techniques.

Glacial lake drainage in Patagonia (13-8 kyr) and response of the adjacent Pacific Ocean

Large freshwater lakes formed in North America and Europe during deglaciation following the Last Glacial Maximum. Rapid drainage of these lakes into the Oceans resulted in abrupt perturbations in climate, including the Younger Dryas and 8.2 kyr cooling events. In the mid-latitudes of the Southern Hemisphere major glacial lakes also formed and drained during deglaciation but little is known about the magnitude, organization and timing of these drainage events and their e ect on regional climate. We use 16 new single-grain optically stimulated luminescence (OSL) dates to de ne three stages of rapid glacial lake drainage in the Lago General Carrera/Lago Buenos Aires and Lago Cohrane/ Pueyrredón basins of Patagonia and provide the rst assessment of the e ects of lake drainage on the Paci c Ocean. Lake drainage occurred between 13 and 8 kyr ago and was initially gradual eastward into the Atlantic, then subsequently reorganized westward into the Paci c as new drainage routes opened up during Patagonian Ice Sheet deglaciation. Coupled ocean-atmosphere model experiments using HadCM3 with an imposed freshwater surface “hosing” to simulate glacial lake drainage suggest that a negative salinity anomaly was advected south around Cape Horn, resulting in brief but signi cant impacts on coastal ocean vertical mixing and regional climate.

Height variability from the MIROC - IPCC model for the 20th century compared to that of the TOPEX/POSEIDON altimeter

Planetary waves are key to large-scale dynamical adjustment in the global ocean as they transfer energy from the east to the west side of oceanic basins; they connect the forcing in the ocean interior with the variability at its boundaries: and they change the local heat content, thus coupling oceanic, atmospheric, and biological processes. Planetary waves, mostly of the first baroclinic mode, are observed as distinctive patterns in global time series of sea surface height anomaly (SSHA) and heat storage. The goal of this study is to compare and validate large-scale SSHA signals from coupled ocean-atmosphere general circulation Model for Interdisciplinary Research on Climate (MIROC) with TOPEX/POSEIDON satellite altimeter observations. The last decade of the models` time series is selected for comparison with the altimeter data. The wave patterns are separated from the meso- and large-scale SSHA signals by digital filters calibrated to select the same spectral bands in both model and altimeter data. The band-wise comparison allows for an assessment of the model skill to simulate the dynamical components of the observed wave field. Comparisons regarding both the seasonal cycle and the Rossby wave Held differ significantly among basins. When carried within the same basin, differences can occur between equal latitudes in opposite hemispheres. Furthermore, at some latitudes the MIROC reproduces biannual, annual and semiannual planetary waves with phase speeds and average amplitudes similar to those observed by the altimeter, but with significant differences in phase. (C) 2008 Elsevier Ltd. All rights reserved.

The Southern Ocean and Its Climate in CCSM4

The new Community Climate System Model, version 4 (CCSM4), provides a powerful tool to understand and predict the earth's climate system. Several aspects of the Southern Ocean in the CCSM4 are explored, including the surface climatology and interannual variability, simulation of key climate water masses (Antarctic Bottom Water, Subantarctic Mode Water, and Antarctic Intermediate Water), the transport and structure of the Antarctic Circumpolar Current, and interbasin exchange via the Agulhas and Tasman leakages and at the Brazil-Malvinas Confluence. It is found that the CCSM4 has varying degrees of accuracy in the simulation of the climate of the Southern Ocean when compared with observations. This study has identified aspects of the model that warrant further analysis that will result in a more comprehensive understanding of ocean-atmosphere-ice dynamics and interactions that control the earth's climate and its variability.