Dynamics of nonlinear error growth and season-dependent predictability of El Nino events in the Zebiak-Cane model

With the intermediate-complexity Zebiak-Cane model, we investigate the 'spring predictability barrier' (SPB) problem for El Nino events by tracing the evolution of conditional nonlinear optimal perturbation (CNOP), where CNOP is superimposed on the El Nino events and acts as the initial error with the biggest negative effect on the El Nino prediction. We show that the evolution of CNOP-type errors has obvious seasonal dependence and yields a significant SPB, with the most severe occurring in predictions made before the boreal spring in the growth phase of El Nino. The CNOP-type errors can be classified into two types: one possessing a sea-surface-temperature anomaly pattern with negative anomalies in the equatorial central-western Pacific, positive anomalies in the equatorial eastern Pacific, and a thermocline depth anomaly pattern with positive anomalies along the Equator, and another with patterns almost opposite to those of the former type. In predictions through the spring in the growth phase of El Nino, the initial error with the worst effect on the prediction tends to be the latter type of CNOP error, whereas in predictions through the spring in the decaying phase, the initial error with the biggest negative effect on the prediction is inclined to be the former type of CNOP error. Although the linear singular vector (LSV)-type errors also have patterns similar to the CNOP-type errors, they cover a more localized area than the CNOP-type errors and cause a much smaller prediction error, yielding a less significant SPB. Random errors in the initial conditions are also superimposed on El Nino events to investigate the SPB. We find that, whenever the predictions start, the random errors neither exhibit an obvious season-dependent evolution nor yield a large prediction error, and thus may not be responsible for the SPB phenomenon for El Nino events. These results suggest that the occurrence of the SPB is closely related to particular initial error patterns. The two kinds of CNOP-type error are most likely to cause a significant SPB. They have opposite signs and, consequently, opposite growth behaviours, a result which may demonstrate two dynamical mechanisms of error growth related to SPB: in one case, the errors grow in a manner similar to El Nino; in the other, the errors develop with a tendency opposite to El Nino. The two types of CNOP error may be most likely to provide the information regarding the 'sensitive area' of El Nino-Southern Oscillation (ENSO) predictions. If these types of initial error exist in realistic ENSO predictions and if a target method or a data assimilation approach can filter them, the ENSO forecast skill may be improved. Copyright (C) 2009 Royal Meteorological Society

El Nino phenomenon in simple ocean data assimilation data

To study how the air and sea interact with each other during El Nino/La Nina onsets, extended associate pattern analysis (EAPA) is adopted with the simple ocean data assimilation (SODA) data. The results show that as El Nino/La Nina's parents their behaviors are quite different, there does not exist a relatively independent tropical atmosphere but does exist a relatively independent tropical Pacific Ocean because the air is heated from the bottom surface instead of the top surface and of much stronger baroclinic instability than the sea and has a very large inter-tropical convergence zone covering the most tropical Pacific Ocean. The idea that it is the wester burst and wind convergence, coming from middle latitudes directly that produce the seawater eastward movement and meridional convergence in the upper levels and result in the typical El Nino sea surface temperature warm signal is confirmed again.

The Obduction of Equatorial 13 degrees C Water in the Pacific Identified by a Simulated Passive Tracer

The obduction of equatorial 13 degrees C Water in the Pacific is investigated using a simulated passive tracer of the Consortium for Estimating the Circulation and Climate of the Ocean (ECCO). The result shows that the 13 degrees C Water initialized in the region 8 degrees N-8 degrees S, 130 degrees-90 degrees W enters the surface mixed layer in the eastern tropical Pacific, mainly through upwelling near the equator, in the Costa Rica Dome, and along the coast of Peru. Approximately two-thirds of this obduction occurs within 10 years after the 13 degrees C Water being initialized, with the upper portion of the water mass reaching the surface mixed layer in only about a month. The obduction of the 13 degrees C Water helps to maintain a cool sea surface temperature year-round, equivalent to a surface heat flux of about -6.0 W m(-2) averaged over the eastern tropical Pacific (15 degrees S-15 degrees N, 130 degrees W-eastern boundary) for the period of integration (1993-2006). During El Nino years, when the thermocline deepens as a consequence of the easterly wind weakening, the obduction of the 13 degrees C Water is suppressed, and the reduced vertical entrainment generates a warming anomaly of up to 10 W m(-2) in the eastern tropical Pacific and in particular along the coast of Peru, providing explanations for the warming of sea surface temperature that cannot be accounted for by local winds alone. The situation is reversed during La Nina years.

Interannual horizontal heat advection in the surface mixed layer over the equatorial Pacific Ocean: assimilation versus TAO analysis

The interannual anomalies of horizontal heat advection in the surface mixed layer over the equatorial Pacific Ocean in an assimilation experiment are studied and compared with existing observational analyses. The assimilation builds upon a hindcast study that has produced a good simulation of the observed equatorial currents and optimizes the simulation of the Reynolds sea surface temperature (SST) data. The comparison suggests that the assimilation has improved the simulation of the interannual horizontal heat advection of the surface mixed layer significantly. During periods of interrupted current measurements, the assimilation is shown to produce more meaningful anomalies of the heat advection than the interpolation of the observational data does. The assimilation also shows that the eddy heat flux due to the correlation between high-frequency current and SST variations, which is largely overlooked by the existing observational analyses, is important for the interannual SST balance over the equatorial Pacific. The interannual horizontal heat advection anomalies are found to be sensitive to SST errors where oceanic currents are strong, which is a challenge for ENSO prediction. The study further suggests that the observational analyses of the tropical SST balance based on the TAO and the Reynolds SST data contain significant errors due to the large gradient errors in the Reynolds SST data, which are amplified into the advection anomalies by the large equatorial currents.

Progress of large-scale air-sea interaction studies in China

This paper summarizes the progress of large-scale air-sea interaction studies that has been achieved in China in the four-year period from July 1998 to July 2002, including seven aspects in the area of the air-sea interaction, namely air-sea interaction related to the tropical Pacific Ocean, monsoon-related air-sea interaction, air-sea interaction in the north Pacific Ocean, air-sea interaction in the Indian Ocean, air-sea interactions in the global oceans, field experiments, and oceanic cruise surveys. However more attention has been paid to the first and the second aspects because a large number of papers in the reference literature for preparing and organizing this paper are concentrated in the tropical Pacific Ocean, such as the ENSO process with its climatic effects and dynamics, and the monsoon-related air-sea interaction. The literature also involves various phenomena with their different time and spatial scales such as intraseasonal, annual, interannual, and interdecadal variabilities in the atmosphere/ocean interaction system, reflecting the contemporary themes in the four-year period at the beginning of an era from the post-TOGA to CLIVAR studies. Apparently, it is a difficult task to summarize the great progress in this area, as it is extracted from a large quantity of literature, although the authors tried very hard.

Equatorial pacific SSTA-related decadal variations of potential predictability of ENSO and interannual climate

Based on analysis of NCEP reanalysis data and SST indices of the recent 50 years, decadal changes of the potential predictability of ENSO and interannual climate anomalies were investigated. Autocorrelation of Nino3 SST anomalies (SSTA) and correlation between atmospheric anomalies fields and Nino3 SSTA exhibit obvious variation in different decades, which indicates that Nino3 SSTA-related potential predictability of ENSO and interannual climate anomalies has significant decadal changes. Time around 1977 is not only a shift point of climate on the interdecadal time scale but also a catastrophe point of potential predictability of ENSO and interannual climate. As a whole, ENSO and the PNA pattern in boreal winter are more predictable in 1980s than in 1960s and 1970s, while the Nino3 SSTA-related potential predictability of the Indian monsoon and the East Asian Monsoon is lower in 1980s than in 1960s and 1970s.

A hybrid coupled ocean-atmosphere model and ENSO prediction study

A hybrid coupled ocean-atmosphere model is designed, which consists of a global AGCM and a simple anomaly ocean model in the tropical Pacific. Retroactive experimental predictions initiated in each season from 1979 to 1994 are performed. Analyses indicate that: (1) The overall predictive capability of this model for SSTA over the central-eastern tropical Pacific can reach one year, and the error is not larger than 0.8 degrees C. (2) The prediction skill depends greatly on the season when forecasts start. However, the phenomenon of SPB (spring prediction barrier) is not found in the model. (3) The ensemble forecast method can effectively improve prediction results. A new initialization scheme is discussed.

Intraseasonal variability of Indian Ocean sea surface temperature during boreal winter: Madden-Julian Oscillation versus submonthly forcing and processes

[ 1] Intraseasonal variability of Indian Ocean sea surface temperature (SST) during boreal winter is investigated by analyzing available data and a suite of solutions to an ocean general circulation model for 1998 - 2004. This period covers the QuikSCAT and Tropical Rainfall Measuring Mission (TRMM) observations. Impacts of the 30 - 90 day and 10 - 30 day atmospheric intraseasonal oscillations (ISOs) are examined separately, with the former dominated by the Madden-Julian Oscillation (MJO) and the latter dominated by convectively coupled Rossby and Kelvin waves. The maximum variation of intraseasonal SST occurs at 10 degrees S - 2 degrees S in the wintertime Intertropical Convergence Zone (ITCZ), where the mixed layer is thin and intraseasonal wind speed reaches its maximum. The observed maximum warming ( cooling) averaged over ( 60 degrees E - 85 degrees E, 10 degrees S - 3 degrees S) is 1.13 degrees C ( - 0.97 degrees C) for the period of interest, with a standard deviation of 0.39 degrees C in winter. This SST change is forced predominantly by the MJO. While the MJO causes a basin-wide cooling ( warming) in the ITCZ region, submonthly ISOs cause a more complex SST structure that propagates southwestward in the western-central basin and southeastward in the eastern ocean. On both the MJO and submonthly timescales, winds are the deterministic factor for the SST variability. Short-wave radiation generally plays a secondary role, and effects of precipitation are negligible. The dominant role of winds results roughly equally from wind speed and stress forcing. Wind speed affects SST by altering turbulent heat fluxes and entrainment cooling. Wind stress affects SST via several local and remote oceanic processes.

IMPACTS OF THE TROPICAL PACIFIC COUPLED PROCESS ON THE INTERANNUAL VARIABILITY IN THE INDIAN OCEAN

The basic features of climatology and interannual variations of tropical Pacific and Indian Oceans were analyzed using a coupled general circulation model (CGCM), which was constituted with an intermediate 2.5-layer ocean model and atmosphere model ECHAM4. The CGCM well captures the spatial and temporal structure of the Pacific El Nino-Southern Oscillation (ENSO) and the variability features in the tropical Indian Ocean. The influence of Pacific air-sea coupled process on the Indian Ocean variability was investigated carefully by conducting numerical experiments. Results show that the occurrence frequency of positive/negative Indian Ocean Dipole (IOD) event will decrease/increase with the presence/absence of the coupled process in the Pacific Ocean. Further analysis demonstrated that the air-sea coupled process in the Pacific Ocean affects the IOD variability mainly by influencing the zonal gradient of thermocline via modulating the background sea surface wind.

Variability of surface circulation in the South China Sea from satellite altimeter data

11-year satellite altimeter sea surface height (SSH) anomaly data from January 1993 to December 2003 are used to present the dominant spatial patterns and temporal variations of the South China Sea (SCS) surface circulation through Empirical Orthogonal Function (EOF) analysis. The first three EOF modes show the obvious seasonal variations of SSH in the SCS. EOF mode one is generally characterized by a basin-wide circulation. Mode two describes the double-cell basin scale circulation structure. The two cells were located off west of the Luzon Island and southeast of Vietnam, respectively. EOF mode three presents the mesoscale eddy structure in the western SCS, which develops into a strong cyclonic eddy rapidly from July to September. EOF mode one and mode three are also embedded with interannual signals, indicating that the SCS surface circulation variation is influenced by El Nino events prominently. The strong El Nino of 1997/98 obviously changed the SCS circulation structure. This study also shows that there existed a series of mesoscale eddies in the western SCS, and their temporal variation indicates intra-seasonal and interannual signals.

Variation in joint mode of the tropical Indian-Pacific Ocean thermodynamic anomaly

Previous research has defined the index of the Indian-Pacific thermodynamic anomaly joint mode (IPTAJM) and suggested that the winter IPTAJM has an important impact on summer rainfall over China. However, the possible causes for the interannual and decadal variability of the IPTAJM are still unclear. Therefore, this work investigates zonal displacements of both the western Pacific warm pool (WPWP) and the eastern Indian Ocean warm pool (EIOWP). The relationships between the WPWP and the EIOWP and the IPTAJM are each examined, and then the impacts of the zonal wind anomalies over the equatorial Pacific and Indian Oceans on the IPTAJM are studied. The WPWP eastern edge anomaly displays significant interannual and decadal variability and experienced a regime shift in about 1976 and 1998, whereas the EIOWP western edge exhibits only distinct interannual variability. The decadal variability of the IPTAJM may be mainly caused by both the zonal migration of the WPWP and the 850 hPa zonal wind anomaly over the central equatorial Pacific. On the other hand, the zonal migrations of both the WPWP and the EIOWP and the zonal wind anomalies over the central equatorial Pacific and the eastern equatorial Indian Ocean may be all responsible for the interannual variability of the IPTAJM.

Mixed-layer water oscillations in tropical Pacific for ENSO cycle

The main modes of interannal variabilities of thermocline and sea surface wind stress in the tropical Pacific and their interactions are investigated, which show the following results. (1) The thermocline anomalies in the tropical Pacific have a zonal dipole pattern with 160 W as its axis and a meridional seesaw pattern with 6-8 degrees N as its transverse axis. The meridional oscillation has a phase lag of about 90 to the zonal oscillation, both oscillations get together to form the El Nino/La Nina cycle, which behaves as a mixed layer water oscillates anticlockwise within the tropical Pacific basin between equator and 12 degrees N. (2) There are two main patterns of wind stress anomalies in the tropical Pacific, of which the first component caused by trade wind anomaly is characterized by the zonal wind stress anomalies and its corresponding divergences field in the equatorial Pacific, and the abnormal cross- equatorial flow wind stress and its corresponding divergence field, which has a sign opposite to that of the equatorial region, in the off-equator of the tropical North Pacific, and the second component represents the wind stress anomalies and corresponding divergences caused by the ITCZ anomaly. (3) The trade winds anomaly plays a decisive role in the strength and phase transition of the ENSO cycle, which results in the sea level tilting, provides an initial potential energy to the mixed layer water oscillation, and causes the opposite thermocline displacement between the west side and east side of the equator and also between the equator and 12 degrees N of the North Pacific basin, therefore determines the amplitude and route for ENSO cycle. The ITCZ anomaly has some effects on the phase transition. (4) The thermal anomaly of the tropical western Pacific causes the wind stress anomaly and extends eastward along the equator accompanied with the mixed layer water oscillation in the equatorial Pacific, which causes the trade winds anomaly and produces the anomalous wind stress and the corresponding divergence in favor to conduce the oscillation, which in turn intensifies the oscillation. The coupled system of ocean-atmosphere interactions and the inertia gravity of the mixed layer water oscillation provide together a phase-switching mechanism and interannual memory for the ENSO cycle. In conclusion, the ENSO cycle essentially is an inertial oscillation of the mixed layer water induced by both the trade winds anomaly and the coupled ocean-atmosphere interaction in the tropical Pacific basin between the equator and 12 degrees N. When the force produced by the coupled ocean-atmosphere interaction is larger than or equal to the resistance caused by the mixed layer water oscillation, the oscillation will be stronger or maintain as it is, while when the force is less than the resistance, the oscillation will be weaker, even break.

末次冰消期以来东亚季风演化：来自东海和北黄海沉积孢粉记录

本论文通过对冲绳海槽南部的MD05-2908孔、北黄海的B-U35、B-U41、B-L44孔和北黄海的299份表层样中的孢粉组成、堆积速率、有机碳氮同位素的综合研究，分析了北黄海的孢粉来源、传播动力和途径，并与其它古气候记录对比，重建了末次冰消期以来东亚季风的演化。 MD05-2908孔孢粉组合特征显示6800-3800a BP间气候温暖干燥，其间存在着两次明显的温度降低、湿度增加阶段，即6800-6600a BP和5800-4850a BP；3800a BP以后，气候变得比较寒冷湿润，中世纪暖期和小冰期分别位于990-480a BP和480-230a BP。 北黄海表层孢粉百分含量等值线图显示北黄海不同海域孢粉来源不同，北黄 海东部海域(123.5°E以东的陆架区域)孢粉主要来源于鸭绿江；北黄海中部海域(123.5°E-122.6°E的陆架区域)孢粉来源比较复杂，既有来自辽东半岛、沿海岛屿、山东半岛，还有由黄海暖流携带而来的；北黄海西部海域(122.6°E以西的陆架区域)孢粉以黄海沿岸流携带而来的黄河来源的为主。 北黄海表层孢粉百分含量等值线图还显示，北黄海东部海域孢粉分布主要受 鸭绿江冲淡水影响，波浪和潮汐的作用较弱；北黄海中部海域孢粉分布受黄海暖流和辽东沿岸流的分支共同影响；黄海暖流与黄海沿岸流在山东半岛东北部相遇，形成一逆时针方向的涡旋，该涡旋控制着北黄海西部海域孢粉的分布。 B-U35孔、B-U41和B-L44孔孢粉组合特征显示，12830-12350a BP时气候 寒冷干燥，植被以旱生草本植物为主，时间上可能对应YD事件；12350-10100a BP时植被以针叶林为主，林下生长着大量蕨类植物，气候寒冷湿润，可能对应着YD事件；10100-6600a BP时植被以落叶阔叶林为主，伴生大量旱生草本植物，气候温暖干燥；6600-5000a BP时植被以针叶-落叶阔叶混交林为主，旱生草本植物含量下降，相对上一阶段温度有所降低，但气候更湿润，可能对应全新世最佳期；5000-4000a BP时针叶林所占比例扩大，温度再次降低，湿度变化不明显；4000-500a BP旱生草本植物含量增加，气候寒冷干燥；500a BP-至今，受人类活动影响，植被遭到严重破坏。 B-U35孔和B-U41孔淡水水生植物花粉百分含量曲线显示11750a BP以来黄河河口地区沼泽湿地面积几经变化。10100-6600a BP时沼泽湿地范围在上一阶段的基础上进一步扩大；6600-2850a BP间沼泽湿地大面积缩小，水生植物属种式微；2850a BP之后黄河河口地区沼泽湿地面积再次扩大，水生植物重新繁盛。 MD05-2908孔蕨类孢子百分含量被用来作为中国东南部季风演化的代用指标。结果显示6800-3800a BP时夏季风较强，其间有两次夏季风减弱阶段，分别位于6800-6600a BP和5800-4850a BP，3800a BP以来，夏季风明显减弱，此期的夏季风减弱除受太阳辐射量的减少影响外，可能还与El Nino的频繁爆发有关。 B-U35孔和B-U41孔蒿属、藜科、麻黄科百分含量被用来作为中国中北部夏季风演化的代用指标。研究表明，11750a BP以来东亚夏季风强度在10100-6600a BP时，显著加强，6600a BP左右，夏季风开始减弱，4000a BP以后夏季风显著减弱。 将MD05-2908孔、B-U35孔和B-U41孔夏季风代用指标进行对比发现，中晚全新世以来中国南-北湿度变化特征正好相反，东亚夏季风强度变化以及由此引起的雨锋的南进北撤可能是导致中国不同区域中晚全新世以来降雨量不同的主要原因。