Interdecadal change in the Northern Hemisphere seasonal climate prediction skill: part I. The leading forced mode of atmospheric circulation

Using observations and 1-month lead hindcast data from six coupled atmosphere–ocean climate models, this study investigates the interdecadal change in the leading maximum covariance analysis mode (MCA1) of atmospheric circulation in response to the changes in the El Niño and Southern Oscillation (ENSO) occurred around late 1970s. We focus on boreal winter climate variability and predictability over the North Pacific–North American (NPNA) region using December–January–February prediction initiated from November 1st in the period of 1960–1980 (P1) and 1981–2001 (P2). Observed analysis reveals that ENSO variability, the related tropical convective activity, and thus the MCA1 are considerably enhanced from P1 to P2. As a result, surface climate anomalies over the NPNA are more significantly correlated with the MCA1 in P2 than P1, particularly over North America. The six coupled models and their multi-model ensemble not only are capable of capturing the interdecadal change of the MCA1 and its relationship with surface air temperature and precipitation over the NPNA regions but also have significantly higher forecast skills for the MCA1 and the surface climate anomalies in P2 than P1. However, models have systematic biases in the spatial distribution of the MCA1. It is demonstrated that the interdecadal change in the MCA1 should contribute to the improved forecast skill of the NPNA climate during recent epoch.     

Decadal predictability and forecast skill

The “potential predictability” of the climate system is the upper limit of available forecast skill and can be characterized by the ratio p of the predictable variance to the total variance. While the potential predictability of the actual climate system is unknown its analog q may be obtained for a model of the climate system. The usual correlation skill score r and the mean square skill score M are functions of p in the case of actual forecasts and potential correlation ρ and potential mean square skill score $\mathcal{M}$ are the same functions of q in the idealized model context. In the large ensemble limit the connection between model-based potential predictability and skill scores is particularly straightforward with $q=\rho^{2}=\mathcal{M}.$ Decadal predictions of annual mean temperature produced with the Canadian Centre for Climate Modelling and Analysis coupled climate model are analyzed for information on decadal climate predictability and actual forecast skill. Initialized forecast results are compared with the results of uninitialized climate simulations. Model-based values of potential predictability q and potential correlation skill ρ are obtained and ρ is compared with the actual forecast correlation skill r. The skill of externally forced and internally generated components of the variability are separately estimated. As expected, ρ > r and both decline with forecast range τ, at least for the first five years. The decline of skill is associated mainly with the decline of the skill of the internally generated component. The potential and actual skill of a forecast of time-averaged temperature depends on the averaging period. The skill of uninitialized simulations is low for short averaging times and increases as averaging time increases. By contrast, skill is high at short averaging times for forecasts initialized from observations and declines as averaging times increase to about three years, then increases somewhat at longer averaging times. The skills of the initialized forecasts and uninitialized simulations begin to converge for longer averaging times. The potential correlation skill ρ of the externally forced component of temperature is largest at tropical latitudes and the skill of the internally generated component is largest over the North Atlantic, parts of the Southern Ocean and to some extent the North Pacific. Potential skill over extratropical land is somewhat weaker than over oceans. The distribution of actual correlation skill r is broadly similar to that of potential skill for the externally forced component but less so for the internally generated component. Differences in potential and actual skill suggest where improvements in the forecast system might be found.     

Seasonal prediction skill of ECMWF System 4 and NCEP CFSv2 retrospective forecast for the Northern Hemisphere Winter

The seasonal prediction skill for the Northern Hemisphere winter is assessed using retrospective predictions (1982–2010) from the ECMWF System 4 (Sys4) and National Center for Environmental Prediction (NCEP) CFS version 2 (CFSv2) coupled atmosphere–ocean seasonal climate prediction systems. Sys4 shows a cold bias in the equatorial Pacific but a warm bias is found in the North Pacific and part of the North Atlantic. The CFSv2 has strong warm bias from the cold tongue region of the eastern Pacific to the equatorial central Pacific and cold bias in broad areas over the North Pacific and the North Atlantic. A cold bias in the Southern Hemisphere is common in both reforecasts. In addition, excessive precipitation is found in the equatorial Pacific, the equatorial Indian Ocean and the western Pacific in Sys4, and in the South Pacific, the southern Indian Ocean and the western Pacific in CFSv2. A dry bias is found for both modeling systems over South America and northern Australia. The mean prediction skill of 2 meter temperature (2mT) and precipitation anomalies are greater over the tropics than the extra-tropics and also greater over ocean than land. The prediction skill of tropical 2mT and precipitation is greater in strong El Nino Southern Oscillation (ENSO) winters than in weak ENSO winters. Both models predict the year-to-year ENSO variation quite accurately, although sea surface temperature trend bias in CFSv2 over the tropical Pacific results in lower prediction skill for the CFSv2 relative to the Sys4. Both models capture the main ENSO teleconnection pattern of strong anomalies over the tropics, the North Pacific and the North America. However, both models have difficulty in forecasting the year-to-year winter temperature variability over the US and northern Europe.     

How predictable is the northern hemisphere summer upper-tropospheric circulation?

The retrospective forecast skill of three coupled climate models (NCEP CFS, GFDL CM2.1, and CAWCR POAMA 1.5) and their multi-model ensemble (MME) is evaluated, focusing on the Northern Hemisphere (NH) summer upper-tropospheric circulation along with surface temperature and precipitation for the 25-year period of 1981–2005. The seasonal prediction skill for the NH 200-hPa geopotential height basically comes from the coupled models’ ability in predicting the first two empirical orthogonal function (EOF) modes of interannual variability, because the models cannot replicate the residual higher modes. The first two leading EOF modes of the summer 200-hPa circulation account for about 84% (35.4%) of the total variability over the NH tropics (extratropics) and offer a hint of realizable potential predictability. The MME is able to predict both spatial and temporal characteristics of the first EOF mode (EOF1) even at a 5-month lead (January initial condition) with a pattern correlation coefficient (PCC) skill of 0.96 and a temporal correlation coefficient (TCC) skill of 0.62. This long-lead predictability of the EOF1 comes mainly from the prolonged impacts of El Niño-Southern Oscillation (ENSO) as the EOF1 tends to occur during the summer after the mature phase of ENSO. The second EOF mode (EOF2), on the other hand, is related to the developing ENSO and also the interdecadal variability of the sea surface temperature over the North Pacific and North Atlantic Ocean. The MME also captures the EOF2 at a 5-month lead with a PCC skill of 0.87 and a TCC skill of 0.67, but these skills are mainly obtained from the zonally symmetric component of the EOF2, not the prominent wavelike structure, the so-called circumglobal teleconnection (CGT) pattern. In both observation and the 1-month lead MME prediction, the first two leading modes are accompanied by significant rainfall and surface air temperature anomalies in the continental regions of the NH extratropics. The MME’s success in predicting the EOF1 (EOF2) is likely to lead to a better prediction of JJA precipitation anomalies over East Asia and the North Pacific (central and southern Europe and western North America).     

The seasonal predictability of the Asian summer monsoon in a two-tiered forecast system

An extensive set of boreal summer seasonal hindcasts from a two tier system is compared with corresponding seasonal hindcasts from two other coupled ocean–atmosphere models for their seasonal prediction skill (for precipitation and surface temperature) of the Asian summer monsoon. The unique aspect of the two-tier system is that it is at relatively high resolution and the SST forcing is uniquely bias corrected from the multi-model averaged forecasted SST from the two coupled ocean–atmosphere models. Our analysis reveals: (a) The two-tier forecast system has seasonal prediction skill for precipitation that is comparable (over the Southeast Asian monsoon) or even higher (over the South Asian monsoon) than the coupled ocean–atmosphere. For seasonal anomalies of the surface temperature the results are more comparable across models, with all of them showing higher skill than that for precipitation. (b) Despite the improvement from the uncoupled AGCM all models in this study display a deterministic skill for seasonal precipitation anomalies over the Asian summer monsoon region to be weak. But there is useful probabilistic skill for tercile anomalies of precipitation and surface temperature that could be harvested from both the coupled and the uncoupled climate models. (c) Seasonal predictability of the South Asian summer monsoon (rainfall and temperature) does seem to stem from the remote ENSO forcing especially over the Indian monsoon region and the relatively weaker seasonal predictability in the Southeast Asian summer monsoon could be related to the comparatively weaker teleconnection with ENSO. The uncoupled AGCM with the bias corrected SST is able to leverage this teleconnection for improved seasonal prediction skill of the South Asian monsoon relative to the coupled models which display large systematic errors of the tropical SST’s.     

Decadal predictability and prediction skill of sea surface temperatures in the South Pacific region

The South Pacific Ocean is a key driver of climate variability within the Southern Hemisphere at different time scales. Previous studies have characterized the main mode of interannual sea surface temperature (SST) variability in that region as a dipolar pattern of SST anomalies that cover subtropical and extratropical latitudes (the South Pacific Ocean Dipole, or SPOD), which is related to precipitation and temperature anomalies over several regions throughout the Southern Hemisphere. Using that relationship and the reported low predictive skill of precipitation anomalies over the Southern Hemisphere, this work explores the predictability and prediction skill of the SPOD in near-term climate hindcasts using a set of state-of-the-art forecast systems. Results show that predictability greatly benefits from initializing the hindcasts beyond the prescribed radiative forcing, and is modulated by known modes of climate variability, namely El Niño-Southern Oscillation and the Interdecadal Pacific Oscillation. Furthermore, the models are capable of simulating the spatial pattern of the observed SPOD even without initialization, which suggests that the key dynamical processes are properly represented. However, the hindcast of the actual phase of the mode is only achieved when the forecast systems are initialized, pointing at SPOD variability to not be radiatively forced but probably internally generated. The comparison with the performance of an empirical prediction based on persistence suggests that initialization may provide skillful information for SST anomalies, outperforming damping processes, up to 2–3 years into the future.     

西北太平洋热带气旋频次的延伸期动力-统计预报方法和评估

介绍了西北太平洋热带气旋(TC)频次的延伸期预报方法,比较了新构建的动力-统计和统计预报模型的预测技巧,并探讨了预报误差来源及改进方向。动力-统计预报模型是基于动力模式预测的热带季节内振荡(ISO)信号及ISO-TC生成的同期统计关系来进行预报;统计预报模型则是基于TC生成的前兆ISO信号建模预报。预报评估结果显示,动力-统计混合预报模型的预报技巧高于统计预报模型,原因在于影响TC次季节变化的前兆信号并不稳定,且随着预报超前时间迅速消散,无法提供有效且稳定的可预报源;相反地,TC生成与同期的ISO背景场显著相关,动力模式对ISO(预报因子)有较好的预报能力,因此动力-统计相结合的预报方法为TC延伸期预报提供了有效途径。虽然目前动力-统计预报模型的预报技巧可达5～6周,但仍有进一步改进和提高的空间。通过对不同类型TC预报技巧检验和误差分析,研究认为年际和年代际背景场对ISO调控TC活动的影响不可忽略,且热带外ISO信号(如罗斯贝波破碎和西风急流强度等)对TC频次和轨迹也有显著影响,这些因子为TC延伸期预报提供了潜在可预报源。     

How are seasonal prediction skills related to models’ performance on mean state and annual cycle?

Given observed initial conditions, how well do coupled atmosphere–ocean models predict precipitation climatology with 1-month lead forecast? And how do the models’ biases in climatology in turn affect prediction of seasonal anomalies? We address these questions based on analysis of 1-month lead retrospective predictions for 21 years of 1981–2001 made by 13 state-of-the-art coupled climate models and their multi-model ensemble (MME). The evaluation of the precipitation climatology is based on a newly designed metrics that consists of the annual mean, the solstitial mode and equinoctial asymmetric mode of the annual cycle, and the rainy season characteristics. We find that the 1-month lead seasonal prediction made by the 13-model ensemble has skills that are much higher than those in individual model ensemble predictions and approached to those in the ERA-40 and NCEP-2 reanalysis in terms of both the precipitation climatology and seasonal anomalies. We also demonstrate that the skill for individual coupled models in predicting seasonal precipitation anomalies is positively correlated with its performances on prediction of the annual mean and annual cycle of precipitation. In addition, the seasonal prediction skill for the tropical SST anomalies, which are the major predictability source of monsoon precipitation in the current coupled models, is closely link to the models’ ability in simulating the SST mean state. Correction of the inherent bias in the mean state is critical for improving the long-lead seasonal prediction. Most individual coupled models reproduce realistically the long-term annual mean precipitation and the first annual cycle (solstitial mode), but they have difficulty in capturing the second annual (equinoctial asymmetric) mode faithfully, especially over the Indian Ocean (IO) and Western North Pacific (WNP) where the seasonal cycle in SST has significant biases. The coupled models replicate the monsoon rain domains very well except in the East Asian subtropical monsoon and the tropical WNP summer monsoon regions. The models also capture the gross features of the seasonal march of the rainy season including onset and withdraw of the Asian–Australian monsoon system over four major sub-domains, but striking deficiencies in the coupled model predictions are observed over the South China Sea and WNP region, where considerable biases exist in both the amplitude and phase of the annual cycle and the summer precipitation amount and its interannual variability are underestimated.     

QBO influence on extratropical predictive skill

The quasi-biennial oscillation (QBO) in the zonal wind in the tropical stratosphere is one of the most predictable aspects of the circulation anywhere in the atmosphere and can be accurately forecast for many months in advance. If the stratospheric QBO systematically (and significantly) affects the tropospheric circulation, it potentially provides a predictable signal useful for seasonal forecasting. The stratospheric QBO itself is generally not well represented in current numerical models, however, including those used for seasonal prediction and this potential may not be exploited by current numerical-model based forecast systems. The purpose of the present study is to ascertain if a knowledge of the state of the QBO can contribute to extratropical boreal winter seasonal forecast skill and, if so, to motivate further research in this area. The investigation is in the context of the second Historical Forecasting Project (HFP2), a state-of-the-art multimodel two-tier ensemble seasonal forecasting system. The first tier, consisting of a prediction of sea surface temperature anomalies (SSTAs), is followed by the second tier which is a prediction of the state of the atmosphere and surface using an AGCM initialized from atmospheric analyses and using the predicted SSTs as boundary conditions. The HFP2 forecasts are successful in capturing the extratropical effects of sea surface temperature anomalies in the equatorial Pacific to the extent that a linear statistical correction based on the NINO3.4 index does not provide additional extratropical skill. By contrast, knowledge of the state of the stratospheric QBO can be used statistically to add extratropical skill centred in the region of the North Atlantic Oscillation. Although the additional skill is modest, the result supports the contention that taking account of the QBO could improve extratropical seasonal forecasting skill. This might be done statistically after the fact, by forcing the QBO state into the forecast model as it runs or, preferably, by using models which correctly represent the physical processes and behaviour of the QBO.     

Decadal prediction skill in a multi-model ensemble

Decadal climate predictions may have skill due to predictable components in boundary conditions (mainly greenhouse gas concentrations but also tropospheric and stratospheric aerosol distributions) and initial conditions (mainly the ocean state). We investigate the skill of temperature and precipitation hindcasts from a multi-model ensemble of four climate forecast systems based on coupled ocean-atmosphere models. Regional variations in skill with and without trend are compared with similarly analysed uninitialised experiments to separate the trend due to monotonically increasing forcings from fluctuations around the trend due to the ocean initial state and aerosol forcings. In temperature most of the skill in both multi-model ensembles comes from the externally forced trends. The rise of the global mean temperature is represented well in the initialised hindcasts, but variations around the trend show little skill beyond the first year due to the absence of volcanic aerosols in the hindcasts and the unpredictability of ENSO. The models have non-trivial skill in hindcasts of North Atlantic sea surface temperature beyond the trend. This skill is highest in the northern North Atlantic in initialised experiments and in the subtropical North Atlantic in uninitialised simulations. A similar result is found in the Pacific Ocean, although the signal is less clear. The uninitialised simulations have good skill beyond the trend in the western North Pacific. The initialised experiments show some skill in the decadal ENSO region in the eastern Pacific, in agreement with previous studies. However, the results in this study are not statistically significant (p?≈?0.1) by themselves. The initialised models also show some skill in forecasting 4-year mean Sahel rainfall at lead times of 1 and 5?years, in agreement with the observed teleconnection from the Atlantic Ocean. Again, the skill is not statistically significant (p?≈?0.2). Furthermore, uninitialised simulations that include volcanic aerosols have similar skill. It is therefore still an open question whether initialisation improves predictions of Sahel rainfall. We conclude that the main source of skill in forecasting temperature is the trend forced by rising greenhouse gas concentrations. The ocean initial state contributes to skill in some regions, but variations in boundary forcings such as aerosols are as important in decadal forecasting.     

Predictability of the western North Pacific summer climate demonstrated by the coupled models of ENSEMBLES

The Asian monsoon system, including the western North Pacific (WNP), East Asian, and Indian monsoons, dominates the climate of the Asia-Indian Ocean-Pacific region, and plays a significant role in the global hydrological and energy cycles. The prediction of monsoons and associated climate features is a major challenge in seasonal time scale climate forecast. In this study, a comprehensive assessment of the interannual predictability of the WNP summer climate has been performed using the 1-month lead retrospective forecasts (hindcasts) of five state-of-the-art coupled models from ENSEMBLES for the period of 1960–2005. Spatial distribution of the temporal correlation coefficients shows that the interannual variation of precipitation is well predicted around the Maritime Continent and east of the Philippines. The high skills for the lower-tropospheric circulation and sea surface temperature (SST) spread over almost the whole WNP. These results indicate that the models in general successfully predict the interannual variation of the WNP summer climate. Two typical indices, the WNP summer precipitation index and the WNP lower-tropospheric circulation index (WNPMI), have been used to quantify the forecast skill. The correlation coefficient between five models’ multi-model ensemble (MME) mean prediction and observations for the WNP summer precipitation index reaches 0.66 during 1979–2005 while it is 0.68 for the WNPMI during 1960–2005. The WNPMI-regressed anomalies of lower-tropospheric winds, SSTs and precipitation are similar between observations and MME. Further analysis suggests that prediction reliability of the WNP summer climate mainly arises from the atmosphere–ocean interaction over the tropical Indian and the tropical Pacific Ocean, implying that continuing improvement in the representation of the air–sea interaction over these regions in CGCMs is a key for long-lead seasonal forecast over the WNP and East Asia. On the other hand, the prediction of the WNP summer climate anomalies exhibits a remarkable spread resulted from uncertainty in initial conditions. The summer anomalies related to the prediction spread, including the lower-tropospheric circulation, SST and precipitation anomalies, show a Pacific-Japan or East Asia-Pacific pattern in the meridional direction over the WNP. Our further investigations suggest that the WNPMI prediction spread arises mainly from the internal dynamics in air–sea interaction over the WNP and Indian Ocean, since the local relationships among the anomalous SST, circulation, and precipitation associated with the spread are similar to those associated with the interannual variation of the WNPMI in both observations and MME. However, the magnitudes of these anomalies related to the spread are weaker, ranging from one third to a half of those anomalies associated with the interannual variation of the WNPMI in MME over the tropical Indian Ocean and subtropical WNP. These results further support that the improvement in the representation of the air–sea interaction over the tropical Indian Ocean and subtropical WNP in CGCMs is a key for reducing the prediction spread and for improving the long-lead seasonal forecast over the WNP and East Asia.     

Global seasonal climate predictability in a two tiered forecast system: part I: boreal summer and fall seasons

This paper shows demonstrable improvement in the global seasonal climate predictability of boreal summer (at zero lead) and fall (at one season lead) seasonal mean precipitation and surface temperature from a two-tiered seasonal hindcast forced with forecasted SST relative to two other contemporary operational coupled ocean–atmosphere climate models. The results from an extensive set of seasonal hindcasts are analyzed to come to this conclusion. This improvement is attributed to: (1) The multi-model bias corrected SST used to force the atmospheric model. (2) The global atmospheric model which is run at a relatively high resolution of 50 km grid resolution compared to the two other coupled ocean–atmosphere models. (3) The physics of the atmospheric model, especially that related to the convective parameterization scheme. The results of the seasonal hindcast are analyzed for both deterministic and probabilistic skill. The probabilistic skill analysis shows that significant forecast skill can be harvested from these seasonal hindcasts relative to the deterministic skill analysis. The paper concludes that the coupled ocean–atmosphere seasonal hindcasts have reached a reasonable fidelity to exploit their SST anomaly forecasts to force such relatively higher resolution two tier prediction experiments to glean further boreal summer and fall seasonal prediction skill.     

Interdecadal variability of the ENSO teleconnection to the wintertime North Pacific

The El Niño/Southern Oscillation (ENSO) strongly influences the large-scale atmospheric circulation over the extratropical North Pacific during boreal winter, which has an important impact on North American winter climate. This study analyses the interdecadal variability of the ENSO teleconnection to the wintertime extratropical North Pacific, over the period 1900–2010, using a range of observationally derived datasets and an ensemble of atmospheric model simulations. The observed teleconnection strength is found to vary substantially over the 20th century. Specifically, 31-year periods in the early-century (1912–1942), mid-century (1946–1976) and the late-century (1980–2010) are identified in the observations when the ENSO teleconnection to the North Pacific circulation are found to be particularly strong, weak and strong respectively. The ENSO teleconnection to the North Pacific in the atmospheric model ensemble is weak in the mid-century period and substantially stronger in the late-century, closely following the variability in the observed ENSO-North Pacific teleconnection. In the early-century, however, the atmospheric model also exhibits a weak teleconnection to the North Pacific, unlike in observations. In a subset of the model realisations that exhibit similar ENSO-North Pacific teleconnection as in observations during the early-century period there are large differences in extratropical circulation but not in equatorial Pacific precipitation anomalies, in contrast to the late-century period. This suggests that the high correlation in the early century period is largely due to internal extratropical variability. The important implications of these results for seasonal predictability and the assessment of seasonal forecasting systems are discussed.

     

Potential predictability of Eurasian snow cover

Potential predictability and skill of simulated Eurasian snow cover are explored using a suite of seasonal ensemble hindcasts (i.e. retrospective forecasts), an ensemble climate simulation (spanning the years 1982–1998) and observations. Using remotely sensed observations of snow cover, we find significant point-wise correlation over the North Atlantic and North Pacific between winter and spring averaged sea-surface temperatures and Eurasian snow cover area. The observed correlation shows no discernible pattern related to the El Niño-Southern Oscillation (ENSO). The hindcasts show correlation patterns similar to the observations. However, the climate simulation shows an exaggerated ENSO pattern. The results underscore the importance of initialization in seasonal climate forecasts, and that the observed potential predictability of Eurasian snowcover cannot be solely attributed to ENSO.     

Hindcast skill and predictability for precipitation and two-meter air temperature anomalies in global circulation models over the Southeast United States

This paper presents an assessment of the seasonal prediction skill of current global circulation models, with a focus on the two-meter air temperature and precipitation over the Southeast United States. The model seasonal hindcasts are analyzed using measures of potential predictability, anomaly correlation, Brier skill score, and Gerrity skill score. The systematic differences in prediction skill of coupled ocean–atmosphere models versus models using prescribed (either observed or predicted) sea surface temperatures (SSTs) are documented. It is found that the predictability and the hindcast skill of the models vary seasonally and spatially. The largest potential predictability (signal-to-noise ratio) of precipitation anywhere in the United States is found in the Southeast in the spring and winter seasons. The maxima in the potential predictability of two-meter air temperature, however, reside outside the Southeast in all seasons. The largest deterministic hindcast skill over the Southeast is found in wintertime precipitation. At the same time, the boreal winter two-meter air temperature hindcasts have the smallest skill. The large wintertime precipitation skill, the lack of corresponding two-meter air temperature hindcast skill, and a lack of precipitation skill in any other season are features common to all three types of models (atmospheric models forced with observed SSTs, atmospheric models forced with predicted SSTs, and coupled ocean–atmosphere models). Atmospheric models with observed SST forcing demonstrate a moderate skill in hindcasting spring-and summertime two-meter air temperature anomalies, whereas coupled models and atmospheric models forced with predicted SSTs lack similar skill. Probabilistic and categorical hindcasts mirror the deterministic findings, i.e., there is very high skill for winter precipitation and none for summer precipitation. When skillful, the models are conservative, such that low-probability hindcasts tend to be overestimates, whereas high-probability hindcasts tend to be underestimates.     

Seasonal differences of model predictability and the impact of SST in the Pacific

Both seasonal potential predictability and the impact of SST in the Pacific on the forecast skill over China are investigated by using a 9-level global atmospheric general circulation model developed at the Institute of Atmospheric Physics under the Chinese Academy of Sciences (IAP9L-AGCM). For each year during 1970 to 1999, the ensemble consists of seven integrations started from consecutive observational daily atmospheric fields and forced by observational monthly SST. For boreal winter, spring and summer,the variance ratios of the SST-forced variability to the total variability and the differences in the spatial correlation coefficients of seasonal mean fields in special years versus normal years are computed respectively. It follows that there are slightly inter-seasonal differences in the model potential predictability in the Tropics. At northern middle and high latitudes, prediction skill is generally low in spring and relatively high either in summer for surface air temperature and middle and upper tropospheric geopotential height or in winter for wind and precipitation. In general, prediction skill rises notably in western China, especially in northwestern China, when SST anomalies (SSTA) in the Nino-3 region are significant. Moreover,particular attention should be paid to the SSTA in the North Pacific (NP) if one aims to predict summer climate over the eastern part of China, i.e., northeastern China, North China and southeastern China.     

Characteristics of tropical Pacific SST predictability in coupled GCM forecasts using the NCEP CFS

The limits of predictability of El Niño and the Southern Oscillation (ENSO) in coupled models are investigated based on retrospective forecasts of sea surface temperature (SST) made with the National Centers for Environmental Prediction (NCEP) coupled forecast system (CFS). The influence of initial uncertainties and model errors associated with coupled ENSO dynamics on forecast error growth are discussed. The total forecast error has maximum values in the equatorial Pacific and its growth is a strong function of season irrespective of lead time. The largest growth of systematic error of SST occurs mainly over the equatorial central and eastern Pacific and near the southeastern coast of the Americas associated with ENSO events. After subtracting the systematic error, the root-mean-square error of the retrospective forecast SST anomaly also shows a clear seasonal dependency associated with what is called spring barrier. The predictability with respect to ENSO phase shows that the phase locking of ENSO to the mean annual cycle has an influence on the seasonal dependence of skill, since the growth phase of ENSO events is more predictable than the decay phase. The overall characteristics of predictability in the coupled system are assessed by comparing the forecast error growth and the error growth between two model forecasts whose initial conditions are 1 month apart. For the ensemble mean, there is fast growth of error associated with initial uncertainties, becoming saturated within 2 months. The subsequent error growth follows the slow coupled mode related the model’s incorrect ENSO dynamics. As a result, the Lorenz curve of the ensemble mean NINO3 index does not grow, because the systematic error is identical to the same target month. In contrast, the errors of individual members grow as fast as forecast error due to the large instability of the coupled system. Because the model errors are so systematic, their influence on the forecast skill is investigated by analyzing the erroneous features in a long simulation. For the ENSO forecasts in CFS, a constant phase shift with respect to lead month is clear, using monthly forecast composite data. This feature is related to the typical ENSO behavior produced by the model that, unlike the observations, has a long life cycle with a JJA peak. Therefore, the systematic errors in the long run are reflected in the forecast skill as a major factor limiting predictability after the impact of initial uncertainties fades out.     

Interdecadal change of interannual variability and predictability of two types of ENSO

A significant interdecadal climate shift of interannual variability and predictability of two types of the El Niño-Southern Oscillation (ENSO), namely the canonical or eastern Pacific (EP)-type and Modoki or central Pacific (CP) type, are investigated. Using the retrospective forecasts of six-state-of-the-art coupled models and their multi-model ensemble (MME) for December–January–February during the period of 1972–2005 along with corresponding observed and reanalyzed data, we examine the climate regime shift that occurred in the winter of 1988/1989 and how the shift affected interannual variability and predictability of two types of ENSO for the two periods of 1972–1988 (hereafter PRE) and 1989–2005 (hereafter POST). The result first shows substantial interdecadal changes of observed sea surface temperature (SST) in mean state and variability over the western and central Pacific attributable to the significant warming trend in the POST period. In the POST period, the SST variability increased (decreased) significantly over the western (eastern) Pacific. The MME realistically reproduces the observed interdecadal changes with 1- and 4-month forecast lead time. It is found that the CP-type ENSO was more prominent and predictable during the POST than the PRE period while there was no apparent difference in the variability and predictability of the EP-type ENSO between two periods. Note that the second empirical orthogonal function mode of the Pacific SST during the POST period represents the CP-type ENSO but that during the PRE period captures the ENSO transition phase. The MME better predicts the former than the latter. We also investigate distinctive regional impacts associated with the two types of ENSO during the two periods.     

A theoretical relationship between probabilistic relative operating characteristic skill and deterministic correlation skill in dynamical seasonal climate prediction

Understanding the relationship between probabilistic and deterministic prediction skills is of important significance for the study of seasonal forecasting and verification. Based on the Brier skill score methodology, we have previously found a theoretical relationship between the probabilistic resolution skill and the deterministic correlation (i.e., anomaly correlation; AC) skill and a lack of necessary or consistent relationship between the probabilistic reliability skill and the deterministic skill in dynamical seasonal prediction. Here, we further theoretically investigate the relationship between the probabilistic relative operating characteristic (ROC) skill and the deterministic skill. The ROC measures the discrimination attribute of probabilistic forecast quality, another important attribute besides the resolution and reliability. With some simplified assumptions, we first derive theoretical expressions for the hit and false-alarm rates that are basic ingredients for the ROC curve, then demonstrate a sole dependence of the ROC curve on the AC, and finally analytically derive a relationship between the related ROC score and the AC. Such a theoretically derived ROC-AC relationship is further examined using dynamical models’ ensemble seasonal hindcasts, which is well verified. The finding here along with our previous findings implies that the discrimination and resolution attributes of probabilistic seasonal forecast skill are intrinsically equivalent to the corresponding deterministic skill, while the reliability appears to be the fundamental attribute of the probabilistic skill that differs from the deterministic skill, which constitutes an understanding of the fundamental similarities and difference between the two types of seasonal forecasting skills and predictability and can offer important implications for the study of seasonal forecasting and verification.

     

Impact of in-consistency between the climate model and its initial conditions on climate prediction

We investigated the influence of dynamical in-consistency of initial conditions on the predictive skill of decadal climate predictions. The investigation builds on the fully coupled global model “Coupled GCM for Earth Simulator” (CFES). In two separate experiments, the ocean component of the coupled model is full-field initialized with two different initial fields from either the same coupled model CFES or the GECCO2 Ocean Synthesis while the atmosphere is initialized from CFES in both cases. Differences between both experiments show that higher SST forecast skill is obtained when initializing with coupled data assimilation initial conditions (CIH) instead of those from GECCO2 (GIH), with the most significant difference in skill obtained over the tropical Pacific at lead year one. High predictive skill of SST over the tropical Pacific seen in CIH reflects the good reproduction of El Niño events at lead year one. In contrast, GIH produces additional erroneous El Niño events. The tropical Pacific skill differences between both runs can be rationalized in terms of the zonal momentum balance between the wind stress and pressure gradient force, which characterizes the upper equatorial Pacific. In GIH, the differences between the oceanic and atmospheric state at initial time leads to imbalance between the zonal wind stress and pressure gradient force over the equatorial Pacific, which leads to the additional pseudo El Niño events and explains reduced predictive skill. The balance can be reestablished if anomaly initialization strategy is applied with GECCO2 initial conditions and improved predictive skill in the tropical Pacific is observed at lead year one. However, initializing the coupled model with self-consistent initial conditions leads to the highest skill of climate prediction in the tropical Pacific by preserving the momentum balance between zonal wind stress and pressure gradient force along the equatorial Pacific.