Using seasonal hindcasts to understand the origin of the equatorial cold tongue bias in CGCMs and its impact on ENSO

The cold equatorial SST bias in the tropical Pacific that is persistent in many coupled OAGCMs severely impacts the fidelity of the simulated climate and variability in this key region, such as the ENSO phenomenon. The classical bias analysis in these models usually concentrates on multi-decadal to centennial time series needed to obtain statistically robust features. Yet, this strategy cannot fully explain how the models errors were generated in the first place. Here, we use seasonal re-forecasts (hindcasts) to track back the origin of this cold bias. As such hindcasts are initialized close to observations, the transient drift leading to the cold bias can be analyzed to distinguish pre-existing errors from errors responding to initial ones. A time sequence of processes involved in the advent of the final mean state errors can then be proposed. We apply this strategy to the ENSEMBLES-FP6 project multi-model hindcasts of the last decades. Four of the five AOGCMs develop a persistent equatorial cold tongue bias within a few months. The associated systematic errors are first assessed separately for the warm and cold ENSO phases. We find that the models are able to reproduce either El Niño or La Niña close to observations, but not both. ENSO composites then show that the spurious equatorial cooling is maximum for El Niño years for the February and August start dates. For these events and at this time of the year, zonal wind errors in the equatorial Pacific are present from the beginning of the simulation and are hypothesized to be at the origin of the equatorial cold bias, generating too strong upwelling conditions. The systematic underestimation of the mixed layer depth in several models can also amplify the growth of the SST bias. The seminal role of these zonal wind errors is further demonstrated by carrying out ocean-only experiments forced by the AOCGCMs daily 10-meter wind. In a case study, we show that for several models, this forcing is sufficient to reproduce the main SST error patterns seen after 1 month in the AOCGCM hindcasts.     

Mechanisms controlling warm water volume interannual variations in the equatorial Pacific: diabatic versus adiabatic processes

Variations of the volume of warm water above the thermocline in the equatorial Pacific are a good predictor of ENSO (El Ni?o/Southern Oscillation) and are thought to be critical for its preconditioning and development. In this study, the Warm Water Volume (WWV) interannual variability is analysed using forced general circulation model experiments and an original method for diagnosing processes responsible for WWV variations. The meridional recharge/discharge to higher latitudes drives 60% of the ENSO-related equatorial WWV variations, while diabatic processes in the eastern equatorial Pacific account for the remaining 40%. Interior meridional transport is partially compensated by western boundary transports, especially in the southern hemisphere. Diabatic equatorial WWV formation (depletions) during La Ni?a (El Ni?o) are explained by enhanced (reduced) diathermal transport through enhanced (reduced) vertical mixing and penetrating solar forcing at the 20°C isotherm depth. The respective contribution of diabatic and adiabatic processes during build-ups/depletions strongly varies from event-to-event. The WWV build-up during neutral ENSO phases (e.g. 1980–1982) is almost entirely controlled by meridional recharge, providing a text-book example for the recharge/discharge oscillator’s theory. On the other hand, diabatic processes are particularly active during the strongest La Ni?a events (1984, 1988, 1999), contributing to more than 70% of the WWV build-up, with heating by penetrative solar fluxes explaining as much as 30% of the total build-up due to a very shallow thermocline in the eastern Pacific. This study does not invalidate the recharge/discharge oscillator theory but rather emphasizes the importance of equatorial diabatic processes and western boundary transports in controlling WWV changes.     

Low and high frequency Madden–Julian oscillations in austral summer: interannual variations

The Madden–Julian oscillation (MJO) is the main component of intraseasonal variability of the tropical convection, with clear climatic impacts at an almost-global scale. Based on satellite observations, it is shown that there are two types of austral-summer MJO events (broadly defined as 30–120 days convective variability with eastward propagation of about 5 m/s). Equatorial MJO events have a period of 30–50 days and tend to be symmetric about the equator, whereas MJO events centered near 8°S tend to have a longer period of 55–100 days. The lower-frequency variability is associated with a strong upper-ocean response, having a clear signature in both sea surface temperature and its diurnal cycle. These two MJO types have different interannual variations, and are modulated by the Indian Ocean Dipole (IOD). Following a negative IOD event, the lower-frequency southern MJO variability increases, while the higher-frequency equatorial MJO strongly diminishes. We propose two possible explanations for this change in properties of the MJO. One possibility is that changes in the background atmospheric circulation after an IOD favour the development of the low-frequency MJO. The other possibility is that the shallower thermocline ridge and mixed layer depth, by enhancing SST intraseasonal variability and thus ocean–atmosphere coupling in the southwest Indian Ocean (the breeding ground of southern MJO onset), favour the lower-frequency southern MJO variability.     

Tidal mixing in the Indonesian Seas and its effect on the tropical climate system

The sensitivity of the tropical climate to tidal mixing in the Indonesian Archipelago (IA) is investigated using a coupled general circulation model. It is shown that the introduction of tidal mixing considerably improves water masses properties in the IA, generating fresh and cold anomalies in the thermocline and salty and cold anomalies at the surface. The subsurface fresh anomalies are advected in the Indian Ocean thermocline and ultimately surface to freshen the western part of the basin whereas surface salty anomalies are advected in the Leuwin current to salt waters along the Australian coast. The ~0.5°C surface cooling in the IA reduces by 20% the overlying deep convection. This improves both the amount and structure of the rainfall and weakens the wind convergence over the IA, relaxes the equatorial Pacific trade winds and strengthens the winds along Java coast. These wind changes causes the thermocline to be deeper in the eastern equatorial Pacific and shallower in the eastern Indian Ocean. The El Nino Southern Oscillation (ENSO) amplitude is therefore slightly reduced while the Indian Ocean Dipole/Zonal Mode (IODZM) variability increases. IODZM precursors, related to ENSO events the preceding winter in this model, are also shown to be more efficient in promoting an IODZM thanks to an enhanced wind/thermocline coupling. Changes in the coupled system in response tidal mixing are as large as those found when closing the Indonesian Throughflow, emphasizing the key role of IA on the Indo-Pacific climate.     

Impact of partial steps and momentum advection schemes in a global ocean circulation model at eddy-permitting resolution

Series of sensitivity tests were performed with a z-coordinate, global eddy-permitting (1/4°) ocean/sea-ice model (the ORCA-R025 model configuration developed for the DRAKKAR project) to carefully evaluate the impact of recent state-of-the-art numerical schemes on model solutions. The combination of an energy–enstrophy conserving (EEN) scheme for momentum advection with a partial step (PS) representation of the bottom topography yields significant improvements in the mean circulation. Well known biases in the representation of western boundary currents, such as in the Atlantic the detachment of the Gulf Stream, the path of the North Atlantic Current, the location of the Confluence, and the strength of the Zapiola Eddy in the south Atlantic, are partly corrected. Similar improvements are found in the Pacific, Indian, and Southern Oceans, and characteristics of the mean flow are generally much closer to observations. Comparisons with other state-of-the-art models show that the ORCA-R025 configuration generally performs better at similar resolution. In addition, the model solution is often comparable to solutions obtained at 1/6 or 1/10° resolution in some aspects concerning mean flow patterns and distribution of eddy kinetic energy. Although the reasons for these improvements are not analyzed in detail in this paper, evidence is shown that the combination of EEN with PS reduces numerical noise near the bottom, which is likely to affect current–topography interactions in a systematic way. We conclude that significant corrections of the mean biases presently seen in general circulation model solutions at eddy-permitting resolution can still be expected from the development of numerical methods, which represent an alternative to increasing resolution.     

Skill assessment of three earth system models with common marine biogeochemistry

We have assessed the ability of a common ocean biogeochemical model, PISCES, to match relevant modern data fields across a range of ocean circulation fields from three distinct Earth system models: IPSL-CM4-LOOP, IPSL-CM5A-LR and CNRM-CM5.1. The first of these Earth system models has contributed to the IPCC 4th assessment report, while the latter two are contributing to the ongoing IPCC 5th assessment report. These models differ with respect to their atmospheric component, ocean subgrid-scale physics and resolution. The simulated vertical distribution of biogeochemical tracers suffer from biases in ocean circulation and a poor representation of the sinking fluxes of matter. Nevertheless, differences between upper and deep ocean model skills significantly point to changes in the underlying model representations of ocean circulation. IPSL-CM5A-LR and CNRM-CM5.1 poorly represent deep-ocean circulation compared to IPSL-CM4-LOOP degrading the vertical distribution of biogeochemical tracers. However, their representations of surface wind, wind stress, mixed-layer depth and geostrophic circulations (e.g., Antarctic Circumpolar Current) have been improved compared to IPSL-CM4-LOOP. These improvements result in a better representation of large-scale structure of biogeochemical fields in the upper ocean. In particular, a deepening of 20–40 m of the summer mixed-layer depth allows to capture the 0–0.5 μgChl L^?1 concentrations class of surface chlorophyll in the Southern Ocean. Further improvements in the representation of the ocean mixed-layer and deep-ocean ventilation are needed for the next generations of models development to better simulate marine biogeochemistry. In order to better constrain ocean dynamics, we suggest that biogeochemical or passive tracer modules should be used routinely for both model development and model intercomparisons.     

Internal and forced variability along a section between Greenland and Portugal in the CLIPPER Atlantic model

Numerical models are used to estimate the meridional overturning and transports along the paths of two hydrographic cruises, carried out in 1997 and 2002 from Greenland to Portugal. We have examined the influence of the different paths of the two cruises and found that it could explain 0.4 to 2 Sv of difference in overturning (the precise value is model-dependent). Models show a decrease in the overturning circulation between 1997 and 2002, with different amplitudes. The CLIPPER ATL6 model reproduces well the observed weakening of the overturning in density coordinates between the cruises; in the model, the change is due to the combination of interannual and high-frequency forcing and internal variability associated with eddies and meanders. Examination of the -coordinate overturning reveals model–data discrepancies: the vertical structure in the models does not change as much as the observed one. The East Greenland current variability is mainly wind-forced in the ATL6 model, while fluctuations due to eddies and instabilities explain a large part of the North Atlantic Current variability. The time-residual transport of dense water and heat due to eddy correlations between currents and properties is small across this section, which is normal to the direction of the main current.