ENSO representation in climate models: from CMIP3 to CMIP5

We analyse the ability of CMIP3 and CMIP5 coupled ocean–atmosphere general circulation models (CGCMs) to simulate the tropical Pacific mean state and El Niño-Southern Oscillation (ENSO). The CMIP5 multi-model ensemble displays an encouraging 30 % reduction of the pervasive cold bias in the western Pacific, but no quantum leap in ENSO performance compared to CMIP3. CMIP3 and CMIP5 can thus be considered as one large ensemble (CMIP3 + CMIP5) for multi-model ENSO analysis. The too large diversity in CMIP3 ENSO amplitude is however reduced by a factor of two in CMIP5 and the ENSO life cycle (location of surface temperature anomalies, seasonal phase locking) is modestly improved. Other fundamental ENSO characteristics such as central Pacific precipitation anomalies however remain poorly represented. The sea surface temperature (SST)-latent heat flux feedback is slightly improved in the CMIP5 ensemble but the wind-SST feedback is still underestimated by 20–50 % and the shortwave-SST feedbacks remain underestimated by a factor of two. The improvement in ENSO amplitudes might therefore result from error compensations. The ability of CMIP models to simulate the SST-shortwave feedback, a major source of erroneous ENSO in CGCMs, is further detailed. In observations, this feedback is strongly nonlinear because the real atmosphere switches from subsident (positive feedback) to convective (negative feedback) regimes under the effect of seasonal and interannual variations. Only one-third of CMIP3 + CMIP5 models reproduce this regime shift, with the other models remaining locked in one of the two regimes. The modelled shortwave feedback nonlinearity increases with ENSO amplitude and the amplitude of this feedback in the spring strongly relates with the models ability to simulate ENSO phase locking. In a final stage, a subset of metrics is proposed in order to synthesize the ability of each CMIP3 and CMIP5 models to simulate ENSO main characteristics and key atmospheric feedbacks.     

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.     

Role of the upper ocean structure in the response of ENSO-like SST variability to global warming

The response of El Niño and Southern Oscillation (ENSO)-like variability to global warming varies comparatively between the two different climate system models, i.e., the Meteorological Research Institute (MRI) and Geophysical Fluid Dynamics Laboratory (GFDL) Coupled General Circulation Models (CGCMs). Here, we examine the role of the simulated upper ocean temperature structure in the different sensitivities of the simulated ENSO variability in the models based on the different level of CO₂ concentrations. In the MRI model, the sea surface temperature (SST) undergoes a rather drastic modification, namely a tendency toward a permanent El Niño-like state. This is associated with an enhanced stratification which results in greater ENSO amplitude for the MRI model. On the other hand, the ENSO simulated by GFDL model is hardly modified although the mean temperature in the near surface layer increases. In order to understand the associated mechanisms we carry out a vertical mode decomposition of the mean equatorial stratification and a simplified heat balance analysis using an intermediate tropical Pacific model tuned from the CGCM outputs. It is found that in the MRI model the increased stratification is associated with an enhancement of the zonal advective feedback and the non-linear advection. In the GFDL model, on the other hand, the thermocline variability and associated anomalous vertical advection are reduced in the eastern equatorial Pacific under global warming, which erodes the thermocline feedback and explains why the ENSO amplitude is reduced in a warmer climate in this model. It is suggested that change in stratification associated with global warming impacts the equatorial wave dynamics in a way that enhances the second baroclinic mode over the gravest one, which leads to the change in feedback processes in the CGCMs. Our results illustrate that the upper ocean vertical structure simulated in the CGCMs is a key parameter of the sensitivity of ENSO-like SST variability to global warming.     

CMIP5耦合模式对太平洋年代际振荡的模拟与预估

利用第五次耦合模式比较计划（CMIP5）40个模式的模拟资料和分类集合的方法,评估了耦合模式对20世纪太平洋年代际振荡（PDO）特征的模拟能力。结果表明,CMIP5多数模式对PDO周期有着较好的刻画能力,能模拟出PDO的年代际变化周期。模式对PDO模态空间特征的模拟能力存在较大差异,小部分模式模拟效果较差。进一步的分析表明,对PDO模态模拟较好的第1类模式,能较好地再现热带太平洋与北太平洋海表温度异常（SSTA）年代际变化间的关系,而且热带太平洋SSTA通过大气遥相关影响北太平样海表温度的过程也模拟的较成功。对PDO模态模拟差的模式,不能合理模拟出热带太平洋SSTA对北太平洋海表温度影响的遥相关过程。以上研究也证实了热带太平洋地区海表温度的年代际变率对北太平洋海表温度年代际变率的重要影响,热带太平洋SSTA对北太平洋SSTA的影响是通过大气遥相关实现的。利用CMIP5中等排放情景模拟结果,分析了第1类模式预估的北太平洋年代际变率的特征,发现21世纪北太平洋年代际变率的主要模态为一致的正异常分布且呈现明显的上升趋势,第二模态则表现为类似于20世纪典型PDO的马蹄型SSTA分布。     

太平洋年代际振荡冷、暖背景下ENSO循环的特征

利用英国气象局哈德莱中心的月平均海温距平资料、美国Scripps海洋研究所联合环境分析中心(JEDAC)的海表和次表层海温观测资料以及NCEP/NCAR再分析资料,研究了太平洋年代际振荡(PDO)不同背景下ENSO循环的特征.结果表明,PDO为ENSO循环提供了一个年代际气候背景,在PDO的暖位相时期,El Nino事件发生的频率较高,强度较强;反之,在PDO的冷位相时期,La Nina事件发生的频率较高,强度较强.而且在不同的太平洋年代际振荡背景下,ENSO循环表现出不同的特征.在PDO冷位相时期,发生El Nino(La Nina)事件时,正(负)的SOTA从西太平洋沿温跃层向东传播,正(负)的SSTA从赤道东太平洋向西扩展到中太平洋,ENSO事件先在赤道东太平洋爆发.在PDO的暖位相时期,发生El Nino(La Nina)事件时,正(负)的SOTA首先出现在赤道中太平洋,然后沿温跃层向东传播,正(负)的SSTA从赤道中太平洋向东扩展到东太平洋,ENSO事件首先在中太平洋爆发.这为ENSO预测提供了新的线索.     

ENSO stability in coupled climate models and its association with mean state

In this study, using the Bjerknes stability (BJ) index analysis, we estimate the overall linear El Niño-Southern Oscillation (ENSO) stability and the relative contribution of positive feedbacks and damping processes to the stability in historical simulations of Coupled Model Intercomparison Project Phase 5 (CMIP5) models. When compared with CMIP3 models, the ENSO amplitudes and the ENSO stability as estimated by the BJ index in the CMIP5 models are more converged around the observed, estimated from the atmosphere and ocean reanalysis data sets. The reduced diversity among models in the simulated ENSO stability can be partly attributed to the reduced spread of the thermocline feedback and Ekman feedback terms among the models. However, a systematic bias persists from CMIP3 to CMIP5. In other words, the majority of the CMIP5 models analyzed in this study still underestimate the zonal advective feedback, thermocline feedback and thermodynamic damping terms, when compared with those estimated from reanalysis. This discrepancy turns out to be related with a cold tongue bias in coupled models that causes a weaker atmospheric thermodynamical response to sea surface temperature changes and a weaker oceanic response (zonal currents and zonal thermocline slope) to wind changes.     

ENSO teleconnections in projections of future climate in ECHAM5/MPI-OM

The teleconnections of the El Niño/Southern Oscillation (ENSO) in future climate projections are investigated using results of the coupled climate model ECHAM5/MPI-OM. For this, the IPCC SRES scenario A1B and a quadrupled CO₂ simulation are considered. It is found that changes of the mean state in the tropical Pacific are likely to condition ENSO teleconnections in the Pacific North America (PNA) region and in the North Atlantic European (NAE) region. With increasing greenhouse gas emissions the changes of the mean states in the tropical and sub-tropical Pacific are El Niño-like in this particular model. Sea surface temperatures in the tropical Pacific are increased predominantly in its eastern part and redistribute the precipitation further eastward. The dynamical response of the atmosphere is such that the equatorial east–west (Walker) circulation and the eastern Pacific inverse Hadley circulation are decreased. Over the subtropical East Pacific and North Atlantic the 200 hPa westerly wind is substantially increased. Composite maps of different climate parameters for positive and negative ENSO events are used to reveal changes of the ENSO teleconnections. Mean sea level pressure and upper tropospheric zonal winds indicate an eastward shift of the well-known teleconnection patterns in the PNA region and an increasing North Atlantic oscillation (NAO) like response over the NAE region. Surface temperature and precipitation underline this effect, particularly over the North Pacific and the central North Atlantic. Moreover, in the NAE region the 200 hPa westerly wind is increasingly related to the stationary wave activity. Here the stationary waves appear NAO-like.     

CMIP5多模式对阿留申低压气候特征的模拟检验与预估

利用观测的海温资料和海平面气压资料,检验了CMIP5(Coupled Model Intercomparison Project 5,CMIP5)多模式对阿留申低压(Aleutian Low,AL)特征指数的时空分布和变化的模拟能力;从AL周期及变化趋势等方面,分析了CMIP5模式预估的未来AL的变化特征。结果表明,CMIP5模式及其集合平均能够很好地模拟AL的环流结构,对AL的气候态有着较强的模拟能力,尤其是模式对于东太平洋海表温度的模拟能力直接影响其对于AL的模拟效果。模式的集合平均对变率强度的模拟偏强,且对于变率的模拟效果逊于对气候态的模拟。22个模式中的16个模式能模拟出AL强度指数的年代际变化周期,对年代际周期有着较好的刻画能力。Historical试验下对于AL的变化趋势存在着较大的不确定性,而相对于两种不同排放情景,随着排放的增加,AL更加偏北,强度增强,年际、年代际周期变得更加显著。在两种排放情景下模式的集合平均以及多数模式模拟出AL有着向北和增强的趋势。     

Ocean mixed layer processes in the Pacific Decadal Oscillation in coupled general circulation models

It is investigated how the Pacific Decadal Oscillation (PDO) is simulated differently among various coupled general circulation models (CGCMs), and how it is related to the heat budget of the simulated ocean mixed layer, which includes the surface heat flux and ocean heat transport. For this purpose the dataset of the climate of the 20th Century experiment (20C3M) from nine CGCMs reported to IPCC’s AR4 are used, while the MRI and MIROC models are examined in detail. Detailed analyses of these two CGCMs reveal that the PDO is mainly affected by ocean heat transport rather than surface heat flux, in particular in the MRI model which has a larger contribution of ocean heat transport to the heat budget. It is found that the ocean heat transport due to Ekman advection versus geostrophic advection contributes differently to the PDO in the western and central North Pacific. Specifically, the strength of PDO tends to be larger for CGCMs with a larger ocean heat transport in the region.     

Weak ENSO Asymmetry Due to Weak Nonlinear Air–Sea Interaction in CMIP5 Climate Models

State-of-the-art climate models have long-standing intrinsic biases that limit their simulation and projection capabilities.Significantly weak ENSO asymmetry and weakly nonlinear air–sea interaction over the tropical Pacific was found in CMIP5(Coupled Model Intercomparison Project, Phase 5) climate models compared with observation. The results suggest that a weak nonlinear air–sea interaction may play a role in the weak ENSO asymmetry. Moreover, a weak nonlinearity in air–sea interaction in the models may be associated with the biases in the mean climate—the cold biases in the equatorial central Pacific. The excessive cold tongue bias pushes the deep convection far west to the western Pacific warm pool region and suppresses its development in the central equatorial Pacific. The deep convection has difficulties in further moving to the eastern equatorial Pacific, especially during extreme El Ni o events, which confines the westerly wind anomaly to the western Pacific. This weakens the eastern Pacific El Ni o events, especially the extreme El Ni o events, and thus leads to the weakened ENSO asymmetry in climate models. An accurate mean state structure(especially a realistic cold tongue and deep convection) is critical to reproducing ENSO events in climate models. Our evaluation also revealed that ENSO statistics in CMIP5 climate models are slightly improved compared with those of CMIP3. The weak ENSO asymmetry in CMIP5 is closer to the observation. It is more evident in CMIP5 that strong ENSO activities are usually accompanied by strong ENSO asymmetry, and the diversity of ENSO amplitude is reduced.     

Interdecadal Variability of the ENSO–North Pacific Atmospheric Circulation in Winter

The interdecadal change in the relationship between the El Niño–Southern Oscillation (ENSO) and atmospheric circulation over the North Pacific is investigated using both observational data and an atmospheric general circulation model. There are two prominent modes of winter mid-latitude atmospheric variability in the North Pacific: the West Pacific (WP) teleconnection and the Aleutian Low (AL). The relationship between ENSO and the WP-AL patterns changed notably around the late 1970s. From 1957 to 1975, during the mature phase of ENSO, significant sea surface temperature anomalies (SSTAs) occurred, mainly in the equatorial eastern Pacific Ocean; the associated atmospheric circulation anomaly pattern resembles the negative phase of a WP teleconnection pattern. In contrast, for the 1978–2011 period, significant negative SSTAs were observed in the western and extratropical Pacific in both hemispheres, with some significant positive SSTAs appearing over the eastern Pacific. This is in agreement with the defined regions of a mega-ENSO, the associated atmospheric circulation anomaly pattern resembles the AL mode. Further analysis suggests that a negative–positive anomaly pattern in the 500?hPa geopotential height throughout the entire North Pacific, possibly enhanced by the SSTAs in the extratropical North Pacific associated with the mature phase of ENSO, is responsible for modulating the relationship between ENSO and the North Pacific atmospheric circulation.     

El-Nino Southern Oscillation simulated and predicted in SNU coupled GCMs

The characteristics of the El-Nino Southern Oscillation (ENSO) simulated in free integrations using two versions of the Seoul National University (SNU) ocean–atmosphere coupled global climate model (CGCM) are examined. A revised version of the SNU CGCM is developed by incorporating a reduced air–sea coupling interval (from 1?day to 2?h), a parameterization for cumulus momentum transport, a minimum entrainment rate threshold for convective plumes, and a shortened auto-conversion time scale of cloud water to raindrops. With the revised physical processes, lower tropospheric zonal wind anomalies associated with the ENSO-related sea surface temperature anomalies (SSTA) are represented with more realism than those in the original version. From too weak, the standard deviation of SST over the eastern Pacific becomes too strong in the revised version due to the enhanced air–sea coupling strength and intraseasonal variability associated with ENSO. From the oceanic side, the stronger stratification and the shallower-than-observed thermocline over the eastern Pacific also contribute to the excessive ENSO. The impacts of the revised physical processes on the seasonal predictability are investigated in two sets of the hindcast experiment performed using the two versions of CGCMs. The prediction skill measured by anomaly correlation coefficients of monthly-mean SSTA shows that the new version has a higher skill over the tropical Pacific regions compared to the old version. The better atmospheric responses to the ENSO-related SSTA in the revised version lead to the basin-wide SSTA maintained and developed in a manner that is closer to observations. The symptom of an excessively strong ENSO of the new version in the free integration is not prominent in the hindcast experiment because the thermocline depth over the eastern Pacific is maintained as initialized over the arc of time of the hindcast (7?months).     

Temporal Characteristics of Pacific Decadal Oscillation (PDO) and ENSO and Their Relationship Analyzed with Method of Empirical Mode Decomposition (EMD)

1. IntroductionPacific Decadal Oscillation (PDO) is a long-termENSO-like variability of the North Pacific. It can becharacterized by the first principal component of EOFof the North Pacific SST (Zhu and Yang, 2003; Tren-berth, 1990; Yang and Zhang, 2003). ENSO is thestrongest signal of annular change of global climatesystem (Trenberth, 1997). The spatial pattern of PDOis a wedge similar to El Nino. In the cool (warm)phases of PDO, the central and northwest Pacific is ofwarm (co…     

ENSO phase-locking to the boreal winter in CMIP3 and CMIP5 models

In this study, the El Nino-Southern Oscillation (ENSO) phase-locking to the boreal winter in CMIP3 and CMIP5 models is examined. It is found that the models that are poor at simulating the winter ENSO peak tend to simulate colder seasonal-mean sea-surface temperature (SST) during the boreal summer and associated shallower thermocline depth over the eastern Pacific. These models tend to amplify zonal advection and thermocline depth feedback during boreal summer. In addition, the colder eastern Pacific SST in the model can reduce the summertime mean local convective activity, which tends to weaken the atmospheric response to the ENSO SST forcing. It is also revealed that these models have more serious climatological biases over the tropical Pacific, implying that a realistic simulation of the climatological fields may help to simulate winter ENSO peak better. The models that are poor at simulating ENSO peak in winter also show excessive anomalous SST warming over the western Pacific during boreal winter of the El Nino events, which leads to strong local convective anomalies. This prevents the southward shift of El Nino-related westerly during boreal winter season. Therefore, equatorial westerly is prevailed over the western Pacific to further development of ENSO-related SST during boreal winter. This bias in the SST anomaly is partly due to the climatological dry biases over the central Pacific, which confines ENSO-related precipitation and westerly responses over the western Pacific.     

北半球冬季阿留申低压－冰岛低压相关关系的年代际变化及其模拟

使用多种长期观测和再分析资料,分析了北半球冬季阿留申低压和冰岛低压相关关系的年代际变化。结果表明,两低压存在显著的负相关关系,使北太平洋和北大西洋海平面气压形成跷跷板式的变化(Aleutian Low-Icelandic Low Seesaw,AIS)。此外,AIS还存在显著的年代际变化,在1935~1949年和1980年后较为显著,其余时期并不显著。对1980年代的年代际转变分析表明,太平洋年代际振荡(Pacific Decadal Oscillation,PDO)在1970年代末的位相转变是AIS这次年代际转变的主要原因。PDO由负位相转变为正位相,使全球大部分大洋海表温度升高,而北太平洋海表温度降低,两低压显著变深,低压南部西风增强,从而通过Rossby波的频散效应使两低压强度形成显著负相关。1930年代中期的年代际转变与此类似,但强度较弱。同时,年代际背景的变化也影响到两低压的年际变化。在给定海表温度和海冰分布的驱动下,大气环流模式IAP AGCM4能基本模拟出AIS年代际转变的过程和机理,但仍存在一些偏差。     

Investigating the possibility of a human component in various pacific decadal oscillation indices

The pacific decadal oscillation (PDO) is a mode of natural decadal climate variability, typically defined as the principal component of North Pacific sea surface temperature (SST) anomalies. To remove any global warming signal present in the data, the traditional definition specifies that monthly-mean, global-average SST anomalies are subtracted from the local anomalies. Differences in the warming rates over the globe and the PDO region may therefore be aliased into the PDO index. Here, we examine the possibility of a human component in the PDO, considering three different definitions. The implications of these definitions are explored using SSTs from both observations and simulations of historical and future climate, all projected onto (definition-dependent) observed PDO patterns. In the twenty first century scenarios, a systematic anthropogenic component is found in all three PDO indices. Under the first definition??in which no warming signal is removed??this component is so large that it is also statistically detectable in the observed PDO. Using the second/traditional definition, this component is also large, and arises primarily from the differential warming rates predicted in the North Pacific and over global oceans. Removing the spatial average SST signal in the PDO region (in the third definition) partially solves this problem, but a human signal persists because the predicted pattern of SST response to human forcing projects strongly onto the PDO pattern. This illustrates the importance of separating internally-generated and externally-forced components in the PDO, and suggests that caution should be exercised in using PDO indices for statistical removal of ??natural variability?? effects from observational datasets.     

ENSO Frequency Asymmetry and the Pacific Decadal Oscillation in Observations and 19 CMIP5 Models

Using observational data and the pre-industrial simulations of 19 models from the Coupled Model Intercomparison Project Phase 5(CMIP5), the El Ni o(EN) and La Ni a(LN) events in positive and negative Pacific Decadal Oscillation(PDO) phases are examined. In the observational data, with EN(LN) events the positive(negative) SST anomaly in the equatorial eastern Pacific is much stronger in positive(negative) PDO phases than in negative(positive) phases. Meanwhile,the models cannot reasonably reproduce this difference. Besides, the modulation of ENSO frequency asymmetry by the PDO is explored. Results show that, in the observational data, EN is 300% more(58% less) frequent than LN in positive(negative)PDO phases, which is significant at the 99% confidence level using the Monte Carlo test. Most of the CMIP5 models exhibit results that are consistent with the observational data.     

CMIP5气候模式对中国东部夏季降水年代际变化的模拟性能评估

使用分类集合的方法评估了第五次耦合模式比较计划(CMIP5)多个耦合模式对中国东部夏季降水年代际变化的模拟性能.结果表明,在评估的38个模式中,仅有6个模式(第1类模式)可以成功再现1970年代末中国东部夏季降水年代际变化的主要特征,即长江流域降水偏多、而华北和华南偏少.这些模式模拟的成功归因于它们能较好再现1970年代末东亚夏季风的年代际减弱及相关的环流场的变化,包括东亚沿海的偏北风异常以及西太平洋副热带高压的偏向西南、强度增强等.而对降水年代际变化模拟很差的第2类模式,则模拟出不出东亚夏季风的这种减弱特征.进一步的分析表明,两类CMIP5模式对太平洋年代际振荡(PDO)空间分布特征都有较好的再现能力,但对PDO年代际转变特征的模拟能力则差异较大.第1类模式能很好地模拟出1970年代末热带海洋的增暖和相关的PDO位相由负到正的转换,而第2类模式所模拟的PDO位相转变与观测完全相反,且也不能模拟出热带中东太平洋海洋的年代际增暖及江淮流域夏季的变冷,因此导致该类模式对1970年代末东亚夏季风的减弱和中国东部夏季雨型的年代际转变没有模拟能力.由此也表明,对耦合模式来说,中国夏季降水年代际变化的模拟能力在很大程度上取决于模式对海洋年代际变化信号的模拟.     

SST and ENSO variability and change simulated in historical experiments of CMIP5 models

This work documents the diversity in Coupled Model Inter-comparison Project Phase 5 (CMIP5) models in simulating different aspects of sea surface temperature (SST) variability, particularly those associated with the El Niño–Southern Oscillation (ENSO), as well as the impact of low-frequency variations on the ENSO variability and its global teleconnection. The historical simulations (1870–2005) include 10 models with ensemble member ranging from 3 to 10 that are forced with observed atmospheric composition changes reflecting both natural and anthropogenic forcings. It is shown that the majority of the CMIP5 models capture the relative large SST anomaly variance in the tropical central and eastern Pacific, as well as in North Pacific and North Atlantic. The frequency of ENSO is not well captured by almost all models, particularly for the period of 5–6 years. The low-frequency variations in SST caused by external forcings affect the SST variability and also modify the global teleconnection of ENSO. The models reproduce the global averaged SST low-frequency variations, particularly since 1970s. However, majority of the models are unable to correctly simulate the spatial pattern of the observed SST trends. These results suggest that it is still a challenge to reproduce the features of global historical SST variations with the state-of-the-art coupled general circulation model.     

A brief introduction to the issue of climate and marine fisheries

Climatic variability has profound effects on the distribution, abundance and catch of oceanic fish species around the world. The major modes of this climate variability include the El Niño-Southern Oscillation (ENSO) events, the Pacific Decadal Oscillation (PDO) also referred to as the Interdecadal Pacific Oscillation (IPO), the Indian Ocean Dipole (IOD), the Southern Annular Mode (SAM) and the North Atlantic Oscillation (NAO). Other modes of climate variability include the North Pacific Gyre Oscillation (NPGO), the Atlantic Multidecadal Oscillation (AMO) and the Arctic Oscillation (AO). ENSO events are the principle source of interannual global climate variability, centred in the ocean–atmosphere circulations of the tropical Pacific Ocean and operating on seasonal to interannual time scales. ENSO and the strength of its climate teleconnections are modulated on decadal timescales by the IPO. The time scale of the IOD is seasonal to interannual. The SAM in the mid to high latitudes of the Southern Hemisphere operates in the range of 50–60 days. A prominent teleconnection pattern throughout the year in the Northern Hemisphere is the North Atlantic Oscillation (NAO) which modulates the strength of the westerlies across the North Atlantic in winter, has an impact on the catches of marine fisheries. ENSO events affect the distribution of tuna species in the equatorial Pacific, especially skipjack tuna as well as the abundance and distribution of fish along the western coasts of the Americas. The IOD modulates the distribution of tuna populations and catches in the Indian Ocean, whilst the NAO affects cod stocks heavily exploited in the Atlantic Ocean. The SAM, and its effects on sea surface temperatures influence krill biomass and fisheries catches in the Southern Ocean. The response of oceanic fish stocks to these sources of climatic variability can be used as a guide to the likely effects of climate change on these valuable resources.