全球纬向平均大气环流基本模态年代际变化的凝聚小波分析

应用1871—2008年NCEP／NCAR月平均再分析资料,研究了1948—2008年期间全球纬向平均大气环流基本模态的年代际变化.小波凝聚谱的结果表明全球纬向平均大气环流基本模态存在显著的20年左右周期的年代际变化.小波凝聚位相的结果清楚地显示了纬向平均大气环流基本模态的变化顺序.在20年左右的年代际变化时间尺度上,全球纬向平均温度超前纬向平均位势高度2个月,同时超前纬向平均流10个月出现变化;全球纬向平均位势高度又超前纬向平均流8个月出现变化.全球温度上升(下降), 将使高纬度的纬向平均位势高度降低(升高),中低纬度的纬向平均位势高度升高(降低);进而使得中高纬和热带的纬向平均西风加(减)速或东风减(加)速,同时使极地和副热带的西风减(加)速或东风加(减)速.20世纪70年代末期以来全球显著增暖的异常信号最早出现在南半球对流层顶附近,其次出现在南半球对流层低层、北半球对流层顶附近和北半球对流层低层.     

北半球夏季副热带西风急流变异及其对东亚夏季风气候异常的影响

探讨了东亚地区副热带西风急流(EAJS)位置的年际 变化特征、影响及其可能机制. EAJS南北变动的影响主要集中在亚澳季风区和气候平均的北 半球副热带西风急流轴的南北两侧,这与SOI或ENSO相联系的大气环流响应很不相同,后者的 影响主要表现在中、东太平洋上. 北半球副热带西风急流存在着两个显著不同的模态,其中 一个模态反映的是亚太尤其是东亚地区的西风急流的南北变异,另一个模态出现在150°E～1 20°W的中、东太平洋上. 它们分别联系着不同的太平洋海温异常分布,但都能对夏季200hPa 南亚高压的强度产生影响,尤其是南亚高压的东部,从而可以对我国东部夏季旱涝灾害的形成 产生作用.     

中国气候干湿变率与ENSO的关系及其稳定性

利用1951-01~2000-10中国160站气温和降水月平均资料, 计算了自修正PDSI指数. PDSI指数EOF分析第一模态空间场分布和1951~2000年PDSI指数的变化趋势分布十分相似, 第一模态时间系数反映了空间场随时间的演变情况. 研究发现, EOF分析所揭示的中国气候干湿变率和ENSO有着很好的关系. 这种关系表明, 在典型的ENSO暖状态, 中国大部分地区都偏干, 特别是华北地区更易偏干, 长江以南地区和西北容易偏湿, 而长江中下游地区处于变干和湿的过渡区, 变干或湿不明显. 在典型的ENSO冷状态则情况相反. 而中国气候干湿变率年际和年代际变化都对应着强El Niño事件; 反过来当发生强El Niño事件时, 中国气候干湿变率在年际和年代际尺度上有可能发生剧烈变化. 最近20~30 a中国气候干湿的年代际变化, 特别是华北自20世纪70年代末的变干和西北自80年代中期的变湿, 与ENSO朝更暖的状态变化及全球变暖有着紧密的联系. 1951~2000年中国气候干湿变率和ENSO关系的稳定性分析表明, 中国气候干湿变率和ENSO之间在3~8 a变化周期上存着很好的相关关系, 但这种相关关系不稳定, 存在着年代际变化: 1951~1962和1976~1991年两个时间段两者相关关系很高, 而在1963~1975和1992~2000年两时段内, 两者相关关系较差.     

北太平洋副热带高压年际变异与ENSO 循环之间的选择性相互作用

利用NCEP/NCAR大气再分析资料以及Hadley中心海表温度资料,针对北太平洋副热带高压(简称副高)的完整系统,通过分析超前于ENSO事件的海平面副高年际异常特征及其对ENSO事件的触发作用以及ENSO事件对500 hPa副高和海平面副高的滞后影响,结果表明了北太平洋副热带高压年际变异和ENSO循环之间存在选择性相互作用.即在大多数情况下,一方面,前期海平面副高减弱会导致热带西太平洋表面西风异常,通过海洋平流过程触发El Nino事件在夏季发生发展,在秋冬季成熟; 而另一方面,El Nino事件在秋冬季发展成熟后,增强了赤道中太平洋的对流性热源,通过对异常热源的动力响应,同期和次年夏季500 hPa副高增强,又通过增强的Hadley环流作用,副热带地区下沉运动增强,从而使得次年夏季海平面副高增强,增强的海平面副高又有利于触发下一个La Nina事件.副高年际变异和ENSO循环之间相互作用的选择性主要取决于副高异常是否接近于赤道以及ENSO事件本身的持续性.这种相互作用有利于在热带太平洋海气系统产生准两年振荡.     

现代气候条件下降水变化的“西风模态”空间范围及其影响因子初探

全新世轨道-千年-百年-年代际尺度上中纬度亚洲内陆干旱区的湿度/降水变化与东亚季风区呈现出错/反位相的变化,简称为中纬度内陆干旱区气候变化的"西风模态",但至今对气候变化"西风模态"的适用空间范围及其影响因子缺乏深入研究.针对整个中纬度欧亚大陆(30°~60°N,0°~130°E),使用1960~2010年GPCC,CRU和CPC三个降水数据集的逐月降水资料,分别对冬、夏季降水的年际、年代际信号进行EOF分析,发现现代气候背景下的降水变化"西风模态"在夏季年代际尺度表现最为显著,亚洲中部内陆干旱区部分区域的年代际夏季降水与其东部的中纬度东亚季风区和西部的地中海周边地区均呈现出相反变化的空间格局,在中纬度欧亚大陆表现出"-+-"的空间变化模态.据此划分出气候变化"西风模态"核心区域,即西部以里海(约50°E)为界,东部到河西走廊西界(约90°E),南北界限与亚洲中部干旱区南北界线吻合(约36°~54°N),整体相当于中亚和中国新疆干旱区,本文将其命名为"西风模态核心区".依据逐月NCEP/NCAR再分析数据,探讨了亚洲中部干旱区气候变化"西风模态"的成因,发现其不但与中纬度大气环流的纬向波动传播及印度夏季风降水异常之间的共同作用有直接联系,而且与大西洋多年代际振荡(AMO)也密切相关.     

决定北半球冬季哈得莱环流年际变率的三维大气环流图像:观测分析和数值模拟

本文利用NCEP/NCAR再分析资料,分析了1979~2008年北半球冬季哈得莱(Hadley)环流年际变化的特征,在此基础上,讨论了在观测海温驱动下大气环流模式的模拟结果.观测分析表明,近30年北半球冬季哈得莱环流年际变率的主导模态呈现出空间上的非均匀变化,哈得莱环流圈位于热带部分与其位于副热带部分的强度变化符号相反,这在表征其年际变化特征的另一指标——经向风垂直切变中亦有显著体现.大气环流模式AMIP积分试验结果表明,北半球冬季哈得莱环流强度的上述年际变化源于海温强迫.分析发现,热带中东太平洋和南印度洋暖海温距平强迫导致了哈得莱环流强度年际变化的主导模态呈现出空间上的非均匀变化.ElNio的局地作用和大气桥作用激发的太平洋局地哈得莱环流(30°S~30°N,150°E~90°W)和大西洋局地哈得莱环流(30°S~30°N,90°W~10°W)并非呈现出整体一致的变化,尽管二者纬向平均后分别使气候平均的哈得莱环流圈强度加强和减弱.ElNio遥强迫作用激发的西北太平洋反气旋(0°~30°N,100°E~150°E)使北半球Hadley环流圈强度减弱,ElNio和南印度洋暖海温距平共同强迫出的南印度洋反气旋(30°S~0°,60°E~100°E)使南半球Hadley环流圈的强度亦减弱.上述局地哈得莱环流的变化叠加后,因纬向平均的太平洋局地哈得莱环流强度在(副)热带部分的增强大(小)于纬向平均的大西洋局地哈得莱环流和西北太平洋、南印度洋局地哈得莱环流在(副)热带地区的减弱,结果使得哈得莱环流圈的强度在(副)热带部分偏强(弱);较之南半球,北半球强度变化稍强.因此,北半球冬季哈得莱环流年际变率的主导模态在空间上呈现出非均匀变化.     

利用COSMIC掩星资料研究对流层/下平流层大气比湿对ONI指数的响应

本文利用2006年6月至2014年6月COSMIC掩星观测的水汽廓线, 分析了对流层/下平流层(TLS)比湿信号对ENSO的响应.在数据处理中, 将COSMIC掩星水汽廓线计算得到的全球比湿数据内插为1000~30 hPa区间水平分辨率为5°×5°的三维格网, 在各等压面上求取各格网点去除年/月际信号后的比湿月异常值.然后在对比湿月异常时间序列低通滤波的基础上, 进行经验正交分解(EOF)得到比湿主成分, 并对该主成分信号进行二项式平滑; 接下来将平滑后的主成分信号与反映ENSO活动的ONI指数进行相关处理, 得到各等压面主成分信号相对于ONI指数的相关系数及对应的时间延迟.论文分析了包括Niño-3.4的5个代表性区域的TLS比湿异常主成分信号, 结果表明:在各区域, 采用本文先低通滤波再EOF分解的处理方法获得的TLS比湿异常主成分信号与ONI指数均有很强相关性, 对流层相关系数绝对值达0.8以上, 低平流层高于0.7, 在300~200 hPa的上对流层达到峰值; 各等压面上比湿异常主成分信号相对于ONI指数的时间延迟不尽相同, 在对流层中比湿异常主成分信号普遍滞后于ONI指数1~6个月; 在各区域, 比湿异常主成分与ONI指数相关系数绝对值达到最大的等压面都接近250 hPa, 最大相关系数绝对值均达到0.9以上.进一步对全球250 hPa等压面比湿异常主成分与ONI指数相关性的分析表明:两者强相关的区域主要集中在热带; 在这些强相关区域, 比湿异常主成分相对于ONI主要表现为滞后, 且相关系数越大, 相应的时间延迟越短.

     

南海夏季风爆发与南大洋海温变化之间的联系

利用1979-2009年NCEP第二套大气再分析资料和ERSST海温资料,分析南海夏季风爆发时间的年际和年代际变化特征,考察南海夏季风爆发早晚与南大洋海温之间的联系.主要结果为:(1)南海夏季风爆发时间年际和年代际变化明显,1979-1993年与1994-2009年前后两个阶段爆发时间存在阶段性突变;(2)南海夏季风爆发时间与前期冬季(12-1月)印度洋-南大洋(0-80°E,75°S-50°S)海温、春季(2-3月)太平洋-南大洋(170°E -80°W,75°S-50°S)海温都存在正相关关系,当前期冬、春季南大洋海温偏低(高)时,南海夏季风爆发偏早(晚).南大洋海温信号,无论是年际还是年代际变化,都对南海夏季风爆发具有一定的预测指示作用;(3)南大洋海温异常通过海气相互作用和大气遥相关影响南海夏季风爆发的迟早.当南大洋海温异常偏低(偏高)时,冬季南极涛动偏强(偏弱),同时通过遥相关作用使热带印度洋-西太平洋地区位势高度偏低(偏高)、纬向风加强(减弱),热带大气这种环流异常一直维持到春季4、5月份,位势高度和纬向风异常范围逐步向北扩展并伴随索马里越赤道气流的加强(减弱),从而为南海夏季风爆发偏早(偏晚)提供有利的环流条件.初步分析认为,热带大气环流对南大洋海气相互作用的遥响应与半球际大气质量重新分布引起的南北涛动有关.     

扰动位能与大气环流异常的耦合关系及机理研究

本文在前期局地扰动位能理论工作基础上,进一步着眼于扰动位能与大气环流异常的关系和机理问题,采用奇异值分解(SVD)方法分析了扰动位能与高、中、低层大气动能的耦合相关模态以及年际变率情况,同时考察它与作为大气质量分布的海平面气压之间的相关关系,并从物理上初步探讨了扰动位能与大气动能以及质量的相关机理问题.结果表明,扰动位能自身变率的第一模态同时也是与大气环流异常场耦合变化的主要模态,与之相对应的大气动能和质量场的主模态也是其自身变率的支配模态,说明它们之间存在相互制约的物理机制.而且,热带外地区大气环流主模态特征的形成,与扰动位能兼具全球尺度纬向对称结构和局地尺度纬向非对称性密不可分.另外,冬季南北半球环状模指数与扰动位能的前两阶矩存在非常好的相关关系,而在夏季这种关系明显削弱,仅存在于南半球.     

秋季北极海冰与欧亚冬季气温在年代际和年际尺度上的不同联系

北极海冰的急剧消融在近年来欧亚大陆频发的低温事件中扮演着关键角色.秋季北极海冰的偏少对应着冬季欧亚大陆的低温天气,然而二者的联系在年代际和年际两种时间尺度上存在显著区别.本文运用1979—2012年哈德莱中心第一套海冰覆盖率(HadISST1)、欧洲中心(ERA_Interim)的2 m温度、风场、海平面气压场、高度场等资料,分别研究了年代际和年际时间尺度上前期秋季北极海冰与欧亚冬季气温的联系.结果表明,欧亚和北极地区(0°—160°E,15°N—90°N)的冬季气温具有显著的年代际和年际变化.在年代际尺度上,温度异常分布在21世纪初由北极冷-大陆暖转为北极暖-大陆冷.这一年代际转折与前期秋季整个北极地区的海冰年代际减少联系密切.秋季北极全区海冰年代际偏少对应冬季欧亚大陆中高纬地区的高压异常,有利于北大西洋的暖湿气流北上和北极的冷空气南侵,造成北极暖-大陆冷的温度分布;在年际时间尺度上,温度异常分布主要由第一模态的年际变化部分和第二模态组成,且第一模态包含的年际变率信号也存在显著的年代际变化.年际尺度上全区北极海冰对欧亚冬季气温的影响远不及位于北冰洋西南边缘的巴伦支海、喀拉海和拉普捷夫海西部(30°E—120°E,75°N—85°N)的关键区海冰影响显著.关键区内海冰的偏少会引发冬季的北大西洋涛动负位相,导致北大西洋吹往欧亚大陆的暖湿气流减弱和欧亚大陆中高纬地区的气温偏低.     

Precipitation variability and its relationship to hydrologic variability in Alberta

By applying wavelet‐based empirical orthogonal function (WEOF) analysis to gridded precipitation (P) and empirical orthogonal function (EOF) analysis to gridded air temperature (T), potential evapotranspiration (PET), net precipitation (P‐PET) and runoff (Q), this paper examines the spatial, temporal and frequency patterns of Alberta's climate variability. It was found that only WEOF‐based precipitation patterns, possibly modulated by El Nino Southern Oscillation (ENSO) and Pacific Decadal Oscillation(PDO), delineated Alberta into four major regions which geographically represent northern Alberta Boreal forests, southern Alberta grasslands and Aspen Parklands and the Rocky Mountains and Foothills. The leading mode of wavelet‐based precipitation variability WPC1 showed that between 1900 and 2000, a wet climate dominated northern Alberta with significant 4–8, 11 and 25‐year periodic cycles, while the second mode WPC2 showed that between 1960 and 2000, southern Alberta grasslands were characterized by decreasing precipitation, dominated by 11‐year cycles, and the last two modes WPC3 and WPC4 were characterized by 4–7 and 25‐year cycles and both delineated regions where moisture from the Pacific Ocean penetrated the Rocky Mountains, accounted for much of the sub‐alpine climate. These results show that WEOF is superior to EOF in delineating Alberta precipitation variability to sub‐regions that more closely agree with its eco‐climate regions. Further, it was found that while WPC2 could not explain runoff variations in southern Alberta, WPC1, WPC3 and WPC4 accounted for runoff variability in their respective sub‐regions. Copyright © 2009 John Wiley & Sons, Ltd.     

Asymmetry in the interhemispheric planetary wave-tide link between the two hemispheres

This study assesses the relation between the year-to-year variability of the semidiurnal tides (SDT) observed at high latitudes of both hemispheres and the global stratospheric stationary planetary wave (SPW) with zonal wavenumber S=1 (SPW1) derived from the UKMO temperature data. No significant positive correlation can be identified between the interannual variability of the Northern Hemisphere (NH) SDT and the Southern Hemisphere (SH) SPW1 for austral late-winter months. In contrast, a good consistency is evident for the interannual variations between the SDT observed at Rothera (68°S, 68°W) and the Arctic SPW1 for NH mid-winter months. Since it has been observed that during austral summer the non-migrating SDT often plays a significant role at the latitude of Rothera, a physical link between the SH SDT and the NH SPW is suggested. This asymmetry in the interhemispheric link is also noted in a recent study.     

Relationship and its instability of ENSO — Chinese variations in droughts and wet spells

Monthly data of Self-Calibrated Palmer Drought Severity Index (PDSI) from 1951 to 2000 are calculated using historical precipitation and temperature data for Chinese 160 stations. Temporal and spatial patterns of the first empirical orthogonal function (EOF) of the PDSI reveals a fairly linear trend resulting from trends in precipitation and surface temperature, which is similar to the linear PDSI trend during 1951–2000 calculated using all monthly data. The EOF analysis also reveals that the leading mode correlates significantly with ENSO events in time and space. The ENSO EOF shows that during the typical warm phase of ENSO, surface conditions are drier in most regions of China, especially North China, but wetter than normal in the southern regions of Changjiang River, and Northwest China. During the typical cold phase of ENSO, these anomalies reverse sign. From 1951 to 2000, there are large multi-year to decadal variations in droughts and wet spells over China, which are all closely related to strong El Niño events. In other words, when one strong El Niño event happens, there is a possible big variability in droughts and wet spells over China on the multi-year or decadal scale. Studies also suggest that during the last 2–3 decades climate changes over China, especially North China’s drying and northwest China’s wetting, are closely related to the shift in ENSO towards warmer events and global warming since the late 1970s. The instability of the relationship is also studied. It is revealed that there is a good correlation between ENSO and Chinese variations in droughts and wet spells in the 3–8-year band, but the correlation between ENSO and Chinese variations in droughts and wet spells is instable. Studies suggest that there are decadal changes in the correlation: the wavelet coherency between ENSO and Chinese variations in droughts and wet spells is high during 1951–1962 and 1976–1991, but low during 1963–1975 and 1992–2000.     

Evaluation of ocean-atmospheric indices as predictors for summer streamflow of the Yangtze River based on ROC analysis

Antecedent anomalies of sea surface temperature and atmospheric circulation are important signals for making long-term streamflow forecasts. In this study, four groups of ocean-atmospheric indices, i.e, El Niño Southern Oscillation (ENSO), the Northern Hemisphere atmospheric circulation, the Southern Hemisphere atmospheric circulation (SAC), and the Western Pacific and Indian Ocean SST (WPI), are evaluated for forecasting summer streamflow of the Yangtze River. The gradient boosting regression tree (GBRT) is used to forecast streamflow based on each group of indices. The score based on receiver operating characteristics (ROC) curves, i.e., area under the ROC curve (AUC), is used to evaluate skills of models for identifying the high category and the low category of summer streamflow. It is found that the ENSO group and the SAC group show higher AUC values. Furthermore, both AUC values of GBRT models and individual indices show that the low flow years are easier to be identified than the high flow years. The result of this study highlights the skill from the Southern Hemisphere circulation systems for forecasting summer streamflow of the Yangtze River. Results of relative influences of predictors in GBRT models and AUC of individual indices indicate some key ocean-atmospheric indices, such as the Multivariate ENSO Index and the 500-hPa height of the east of Australia.     

南半球平流层极涡崩溃早晚年环流异常特征

极涡崩溃是平流层大气环流一个重要的变化过程,本文利用31年的再分析资料研究了南半球平流层极涡崩溃早晚年的异常特征.研究结果表明,南半球平流层极涡崩溃偏早年极涡崩溃前后平流层环流场异常表现为整层一致的变化,即都为正温度异常、正位势高度异常和负纬向风异常;而南半球平流层极涡崩溃偏晚年极涡崩溃前后平流层环流场异常的整层一致性的变化不典型,而在符号上与极涡崩溃偏早年的异常相反.与北半球平流层极涡崩溃前后环流异常相反明显不同,南半球平流层极涡崩溃偏早或偏晚年在极涡崩溃前后的环流异常保持相同的性质.进一步分析表明行星波活动在南极极涡的崩溃过程中起到了重要作用,极涡崩溃早年上传行星波比极涡崩溃晚年强,并且持续时间长.通过波流相互作用,行星波的异常使得极涡崩溃早年和晚年10月的平流层高纬地区分别为位势高度正异常和负异常,环流异常持续保持可能最终影响了南半球平流层极涡的崩溃时间.分析显示南半球极涡崩溃偏晚与La Niña事件之间可能存在一定的联系,但在极涡崩溃偏早年与赤道太平洋海表温度异常(SSTA)并无明显关系.     

Characterization of Global Mean Sea Level Variations Observed by TOPEX/POSEIDON Using Empirical Orthogonal Functions

The TOPEX/POSEIDON (T/P) satellite altimeter mission has provided estimates of global mean sea level since late 1992 with a precision of approximately 4 mm. Over the first 3.5 years of the mission, T/P has observed a mean sea level rise of +0.5 mm/year when on-board estimates of the instrument drift are employed (and after correcting for a recently discovered software error), and +2.8 mm/year when an additional external tide gauge-based calibration estimate is used. A preliminary estimate of the error in the latter estimate is 1.3 mm/year, however this issue requires more research. Characterization of the observed sea level variations using Empirical Orthogonal Functions (EOFs) indicates that most of the mean sea level rise can be described by a single mode of the EOF expansion. The spatial characteristics of this mode suggests it is related to the El Nino Southern Oscillation (ENSO) phenomena. EOF analysis of sea level variations from the Semtner/Chervin ocean circulation model reveal a nearly identical mode, although its effect on mean sea level is unknown due to a constant volume constraint used in the model. EOF analysis of measured sea surface temperature (SST) variations also show a mode with similar temporal and spatial structure. However, the concentration of the observed sea level rise in this mode does not preclude the possibility that multiple phenomena have contributed to this mode, thus a link between the observed sea level rise and the ENSO phenomena is only weakly suggested. The absolute value of the observed mean sea level rise will depend on refinements currently being made in the instrument calibration techniques. In addition, the possibility of interannual and decadal variations of global mean sea level requires that a much longer time series of satellite altimetry be collected before variations caused by climate change can be unambiguously detected.     

Interannual variability of Antarctic Oscillation and its influence on East Asian climate during boreal winter and spring

1 Introduction Antarctic Oscillation is a major mode of Southern Hemispheric (SH) extratropical atmospheric circula- tion. The SH Annual Mode represents a zonally sym- metric exchange of mass between polarward of 60°S and 40°S[1], therefore AAO indicates both the intensity of circumpolar low and zonal mean west wind at mid-high latitudes in SH. Positive phase of AAO tends to deepen circumpolar low and enhance west wind at mid-high latitudes in SH. Many studies show that AAO has bar…     

Quasi-biennial and quasi-triennial oscillations in geomagnetic activity indices

Data for geomagnetic activity index aa for 1868–1994 were subjected to spectral analysis for 12 intervals each of 11 consecutive years. In each interval, QBO and QTO (quasi-biennial and quasi-triennial oscillations) were observed at ∼ 2.00, 2.15, 2.40, 2.70 y and ∼ 3.20, 3.40 y, but not all in all intervals. These fluctuations are absent near (2–3 y before and after) the sunspot minima and are present only as 2 or 3 peaks in aa indices, one near or before the sunspot maximum and the other (one or two, generally the larger ones) in the declining phase of the sunspot cycle. Comparison with the solar wind (1965 onwards) showed a fairly good match, indicating that the aa variations were mostly due to similar variations in the solar wind, which must have their origin in solar physical processes. A few aa variations did not match with solar wind. When compared with terrestrial phenomena, no match was found with stratospheric low-latitude zonal wind QBO; but some QTO in aa matched QTO in ENSO (El Nino/ Southern Oscillation). This may or may not be a chance coincidence and needs further exploration.     

Global upper ocean heat content and climate variability

Observational data from the Global Temperature and Salinity Profile Program were used to calculate the upper ocean heat content (OHC) anomaly. The thickness of the upper layer is taken as 300 m for the Pacific/Atlantic Ocean and 150 m for the Indian Ocean since the Indian Ocean has shallower thermoclines. First, the optimal spectral decomposition scheme was used to build up monthly synoptic temperature and salinity dataset for January 1990 to December 2009 on 1° × 1° grids and the same 33 vertical levels as the World Ocean Atlas. Then, the monthly varying upper layer OHC field (H) was obtained. Second, a composite analysis was conducted to obtain the total-time mean OHC field ([`([`(H)])] \bar{\bar{H}} ) and the monthly mean OHC variability ( [(\textH)\tilde] \widetilde{\text{H}} ), which is found an order of magnitude smaller than [^(\textH)] \widehat{\text{H}} . Third, an empirical orthogonal function (EOF) method is conducted on the residue data ( [^(\textH)] \widehat{\text{H}} ), deviating from [(\textH)\tilde] \widetilde{\text{H}}  +  [(\textH)\tilde] \widetilde{\text{H}} , in order to obtain interannual variations of the OHC fields for the three oceans. In the Pacific Ocean, the first two EOF modes account for 51.46% and 13.71% of the variance, representing canonical El Nino/La Nina (EOF-1) and pseudo-El Nino/La Nina (i.e., El Nino Modoki; EOF-2) events. In the Indian Ocean, the first two EOF modes account for 24.27% and 20.94% of the variance, representing basin-scale cooling/warming (EOF-1) and Indian Ocean Dipole (EOF-2) events. In the Atlantic Ocean, the first EOF mode accounts for 49.26% of the variance, representing a basin-scale cooling/warming (EOF-1) event. The second EOF mode accounts for 8.83% of the variance. Different from the Pacific and Indian Oceans, there is no zonal dipole mode in the tropical Atlantic Ocean. Fourth, evident lag correlation coefficients are found between the first principal component of the Pacific Ocean and the Southern Oscillation Index with a maximum correlation coefficient (0.68) at 1-month lead of the EOF-1 and between the second principal component of the Indian Ocean and the Dipole Mode Index with maximum values (around 0.53) at 1–2-month advance of the EOF-2. It implies that OHC anomaly contains climate variability signals.