Improved simulation of the East Asian winter monsoon interannual variation by IAP/LASG AGCMs

JIN Chen-Xi， ZHOU Tin-Jun， GUO Zhun， WU Bond CHEN Xio-Long

aState Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics， Institute of Atmospheric Physics，Chinese Academy of Sciences， Beijing， China；bCollege of Earth Science， University of Chinese Academy of Sciences， Beijing， China；cClimate Change Research Center， Chinese Academy of Sciences， Beijing， China

Improved simulation of the East Asian winter monsoon interannual variation by IAP/LASG AGCMs

JIN Chen-Xia，b， ZHOU Tian-Juna，c， GUO Zhuna， WU Boaand CHEN Xiao-Longa

aState Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics， Institute of Atmospheric Physics，Chinese Academy of Sciences， Beijing， China；bCollege of Earth Science， University of Chinese Academy of Sciences， Beijing， China；cClimate Change Research Center， Chinese Academy of Sciences， Beijing， China

The simulation of the East Asian winter monsoon （EAWM） has been a challenge for climate models. In this study， the performances of two versions of the AGCM developed at the IAP， versions 1 and 2 of the Grid-point Atmospheric Model of the IAP/LASG （GAMIL1 and GAMIL2）， are evaluated in the context of mean state and interannual variation. Signifcant improvements are shown for GAMIL2 in comparison to GAMIL1. The simulated interannual variability of the EAWM， measured by the regional average of 1000 hPa meridional wind over East Asia， has evidently improved； the correlation coefcient with reanalysis data changes from 0.37 in GAMIL1 to 0.71 in GAMIL2. The associated interannual precipitation anomalies are also improved， in terms of both spatial pattern and magnitude. Analysis demonstrates that the improvements result from the better simulation of the El Niño-related Philippine Sea anticyclone （PSAC） in GAMIL2. The improved moist processes，including the stratiform condensation and evaporation in GAMIL2， lead to a reasonable atmospheric heating associated with El Niño in the tropical Pacifc， which further drives the PSAC as a Rossbywave response.

ARTICLE HISTORY

Revised 19 January 2016

Accepted 20 January 2016

East Asian winter monsoon；El Niño； Grid-point

Atmospheric Model of the IAP/LASG； Philippine Sea anticyclone； atmospheric heating

东亚冬季风的模拟性能一直是检验气候模式的重要指标。本文对比分析了LASG新旧两个版本大气环流模式对东亚冬季风年际变率的模拟性能。结果表明，GAMIL2模拟的东亚冬季风指数与观测的相关系数为0.71，相较于GAMIL1（0.37）有明显改进。主要原因是GAMIL2改进了云参数化方案，能够更合理的再现加热率对厄尔尼诺的响应，有助于激发菲律宾海反气旋，从而改进东亚冬季风的模拟。

1. Introduction

The East Asian winter monsoon （EAWM） is one of the most active regional climate systems in boreal winter. The intense variability of the EAWM can bring severe cold surges and heavy snowfall， which may cause serious damage over East Asia， such as the persistent low temperature and icy weather over southern China in January 2008 （Zhou et al. 2009）. In addition to regional impacts， as a coupled extratropical—tropical system， the EAWM also has substantial impacts on far away regions through changes in large-scale deep convection over the equatorial western Pacifc， or changes in the Hadley and Walker circulations （Chang and Lau 1980， 1982； Huang， Zhou，and Chen 2003）. Therefore， understanding the mechanisms of EAWM changes and improving its predictability have been important topics for the Asian climate research community（Ding 1994； Chen， Yang， and Huang 2005； Huang et al. 2012）.

The EAWM is characterized by a cold Siberian—Mongolian high， strong northerlies at the surface， a broad East Asian trough at 500 hPa， and an enhanced East Asian jet stream in the upper troposphere. It is difcult to measure all of the EAWM's characteristics using a single index because of the diference in the EAWM variability between extratropical and tropical East Asia （Wang and Chen 2010；Wang et al. 2010）. The interannual variability of the EAWM can be separated into a northern mode and southern mode using EOF analysis of the 2-m air temperature （Wang et al. 2010）. The northern mode is characterized by enhanced surface pressure over central Siberia， related to the Arctic Oscillation （Gong， Wang， and Zhu 2001； Wang et al. 2010；Chen et al. 2013）. The southern mode is mainly attributed to the Philippine Sea anticyclone （PSAC）， associated with El Niño forcing （Zhang， Sumi， and Kimoto 1996； Jin and Zhou2014）. The main focus of the present study is the EAWM in southern East Asia.

Numerical climate models are useful tools to study the variability of the EAWM. Wei et al. （2014） examined the climatology of the EAWM in 41 fully coupled atmosphere—ocean models from CMIP5 and found that the biases of surface air temperature and precipitation are improved in CMIP5 models compared with CMIP3 models. The interannual variabilities of the EAWM-related circulations are successfully reproduced in CMIP5 multimodel ensembles； for example， about half of CMIP5 models successfully capture the observed ENSO—EAWM relationship （Gong et al. 2014）. Jiang et al. （2015） evaluated 77 coupled global climate models in CMIP3/CMIP5 and revealed an improvement in the performance of models from CMIP3 to CMIP5 with respect to temperature in China， but an underestimation of the strength of the EAWM as manifested by southerly wind anomalies.

Much efort has been devoted to the development of climate models at the LASG， IAP， Chinese Academy of Sciences （see Zhou， et al. 2014b for a summary）. FGOALS-g2 （Li， et al. 2013a） has shown reasonable performance in many aspects among both Chinese （Zhou，et al. 2014a） and international （Bellenger et al. 2013） CMIP5 models. However， the performance of FGOALS-g2 in simulating the EAWM remains unknown. In this study， we begin to address this knowledge gap by assessing the AGCM component of FGOALS-g2， the Grid-point Atmospheric Model of the IAP/LASG （GAMIL）. By comparing AMIP runs of GAMIL driven by observational SST， we aim to evaluate the performance of two versions of GAMIL in simulating the EAWM mean state and interannual variation. We also hope to understand the physical processes that have led to the improvements.

2. Model and data description

Two versions of the IAP/LASG AGCM （GAMIL1 and GAMIL2）are used in this study. They employ the same horizontal（128 × 60） and vertical （26 σ-layers） resolution （Li et al. 2012）. The main diferences between the two versions are their cloud-related processes， including the deep convective parameterization， convective cloud fraction， and microphysical schemes （Li， et al. 2013a， 2013b）. For the deep convective scheme， the CAPE adjustment closure in GAMIL1（Zhang and McFarlane 1995） is replaced by a free-tropospheric quasi-equilibrium closure in GAMIL2 （Zhang 2002）. For nonconvective cloud processes， a two-moment cloud microphysical scheme （Morrison and Gettelman 2008； Shi et al. 2010） is applied in GAMIL2 instead of the one-moment scheme （Rasch and Kristjánsson 1998） in GAMIL1. The AMIP experiments from GAMIL are performed with the same forcing data recommended by CMIP5， such as the solar constant， greenhouse gases， and HadISST data-set. Only the frst member， with the same period from 1979 to 1998， is used for analysis. All the datasets are interpolated into a 2.5° × 2.5° common grid by bilinear interpolation.

The observational and reanalysis datasets include: （1）atmospheric circulation and air temperature from NCEP-2（Kanamitsu et al. 2002）； （2） heating rate in deep convection， shallow convection， and large-scale condensation from the NCEP CFS Reanalysis （CFSR） （Saha et al. 2010）；（3） precipitation data from the GPCP （Adler et al. 2003）； and（4） SST data from ERSST.v3b （Smith et al. 2008）.

3. Results

3.1. Climatology of the EAWM simulated by the IAP/ LASG AGCMs

We assess the performance of the models in simulating the EAWM with a focus on the surface climate. The climatology of 1000-hPa wind， SLP and 2-m temperature derived from the simulations and NCEP-2 reanalysis （hereafter referred to as ‘observation' to facilitate the discussion） are shown in Figure 1. In the observation （Figures 1（a））， a salient structure of the Siberian—Mongolian high is evident at the surface， which induces strong northerlies from 60°N down to the equator. Northwesterlies prevail toward the subtropical Pacifc， north of 30°N， and turn westward along the coast of East Asia into the South China Sea. The above features are reproduced well by the two versions of GAMIL（Figure 1（b） and （c））. We further assess the performance of the models in simulating the 3D features （0°—60°N， 100—140°E） of the EAWM by using a Taylor diagram （Figure 1（d））. Generally， the performance of the models in simulating the 2-m temperature， 500-hPa geopotential height， and 300-hPa zonal wind is better than that of precipitation，SLP and 1000-hPa meridional wind （V1000）. Quantitative comparison suggests that GAMIL2 performs better than GAMIL1 in simulating the climatology of precipitation， SLP，and V1000.

3.2. Interannual variability of the EAWM simulated by the IAP/LASG AGCMs

Following Ji et al. （1997）， an EAWM index is defned as the 1000-hPa meridional wind averaged within the box（10—30°N， 115—130°E） over East Asia. The interannual variation of the EAWM is reproduced well in both GAMIL1 and GAMIL2 （Figure 2（a））. The correlation coefcients of the EAWM index between the observation and the two models indicate signifcant improvement from GAMIL1 （0.37） to GAMIL2 （0.71）. The interannual variability of the EAWM is dominated by El Niño， as demonstrated by the signifcant correlation coefcient of -0.75 between the EAWM indexand Niño3.4 index in the observation. Note that in Table 1 the correlation coefcient between the Niño3.4 index and EAWM index is -0.59 for GAMIL1 and -0.64 for GAMIL2，demonstrating El Niño's forcing of the interannual variability of the EAWM.

Figure 1.Climatology of 1000-hPa winds （vectors； m s-1）， SLP （contours； hPa）， and 2-m temperature （shading； °C） in winter （December—January—February）， derived from （a） the observation （NCEP-2）， （b） GAMIL1， and （c） GAMIL2. （d） Taylor diagram: The observation is considered as the reference （REF）； the angular （vertical） coordinate is the pattern correlation coefcient （ratio of standard deviation）between the model and observation. The nearer the distance between a number and REF， the better the performance of the corresponding model.

To show more robust evidence， we further examine the interannual EAWM anomalous patterns by regressing the observed EAWM index onto the precipitation and 1000-hPa wind anomalies （Figure 2）. In the observation（Figure 2（b））， associated with a high EAWM index， there is a cyclone over the Philippine Sea， with enhanced northeasterly winds over South China. Correspondingly， the precipitation exhibits a dipole pattern， with negative anomalies extending from continental South China to the south of Japan and positive anomalies over the western Pacifc.

Both models reproduce the cyclone and dipole rainfall patterns well. However， the Philippine Sea cyclone in GAMIL1 is weaker than that in the observation and is shifted more eastward （Figure 2（c））. Accordingly， the anomalous precipitation dipole pattern is weak and shifts eastward. Compared with GAMIL1， the performance of GAMIL2 is improved in this regard， along with a better simulation of the precipitation anomaly pattern and its magnitude （Figures 2（d） and （e））. Since the EAWM index is signifcantly negatively correlated with Niño3.4 index（Table 1）， the cyclonic circulations associated with the positive EAWM index shown in Figure 2 are actually indicative of the impact of anticyclonic circulation over the Philippine Sea in El Niño-year winters （Zhang， Sumi， and Kimoto 1996）. We further examine the PSAC simulated by the two versions of GAMIL.

Figure 2.（a） Normalized EAWM index （1000-hPa meridional wind averaged within the box （10—30°N， 115—130°E）） derived from the observation （NCEP-2； black line）， GAMIL1 （red line）， and GAMIL2 （blue line）. Green line: normalized （-1） × Niño3.4 index. （b—e） Regression coefcients of the December—January—February 1000-hPa winds （vectors； m s-1） and precipitation （shaded； mm d-1） with respect to the EAWM index derived from （b） the observation （NCEP-2）， （c） GAMIL1， （d） GAMIL2， and （e） the diference between GAMIL2 and GAMIL1.

Table 1.Time correlation coefcients among the EAWM index （EAWMI） derived from GAMIL and the observation （NCEP-2）， and the Niño3.4 index.

The establishment of the PSAC in boreal winter associated with El Niño events is remotely driven by the warm SST anomalies in the central-eastern Pacifc and enhanced by the local cold SST via wind—evaporation—SST feedback（Wang， Wu， and Fu 2000）. To understand the response of the models to El Niño-related SST forcing， the anomalies of precipitation， 1000-hPa wind and 500-hPa vertical velocity are regressed onto the Niño3.4 index （Figure 3）. In the observation （Figures 3（a） and （b））， in response to the warming SST， positive precipitation anomalies and anomalous upward motion are seen over the equatorial central Pacifc， whereas the western North Pacifc witnesses negative precipitation anomalies along with anomalous sinking motion. The PSAC induced by the Rossby-wave response in the lower troposphere （Wang， Wu， and Fu 2000） and local anomalous sinking motion signifcantly afects the EAWM. In comparison to the observation， the positive （negative）responses of precipitation and upward （sinking） motion over the equatorial central Pacifc （western North Pacifc） to El Niño are weakly simulated by GAMIL1 （Figures 3（c） and（d））， whereas signifcant improvement is seen in GAMIL2（Figures 3（e） and （f））. The improvement of the simulation of the PSAC from GAMIL1 to GAMIL2 ultimately results in a better simulation of the EAWM interannual variation.

Figure 3.Regression coefcients of the December—January—February 1000-hPa winds （vectors； m s-1）， precipitation （left column； shaded；mm d-1）， 500-hPa omega （right column； shaded； pa s-1）， and SST （contours； °C） with respect to the Niño3.4 index derived from （a， b） the observation （NCEP-2）， （c， d） GAMIL1， and （e， f） GAMIL2.

Figure 4.Regression coefcients of the December—January—February heating rate （K d-1） vertically integrated from 1000 to 200 hPa in total moist processes with respect to the Niño3.4 index derived from （a） the observation （CFSR）， （b） GAMIL1， （c） GAMIL2， and （d） the diference between GAMIL2 and GAMIL1.

But why does GAMIL2 better reproduce the PSAC？ Since in the AMIP-type simulation both the central equatorial Pacifc warming and the local colder SST in the western Pacifc are prescribed as the observation， the improved response should result from the improved sensitivity of the model to the SST forcing. In comparison to GAMIL1， cloud-related processes， including deep convective parameterization， the convective cloud fraction and microphysical schemes， are much improved in GAMIL2 （Li， et al. 2013a）. Associated with the enhanced stratiform condensation and evaporation in GAMIL2， the climatology of annual mean water vapor （or relative humidity）， cloud fraction， and in-cloud liquid water path， are also improved （Li， Wang， and Zhang 2014）. In Figure 4， we examine the regression coefcients of the vertically integrated heating rate with respect to the Niño3.4 index. Clearly， the heating rate simulated by GAMIL1 is too weak in comparison to the observation， while the result of GAMIL2 isgreatly improved and closer to the observation. The improved moist physical processes in GAMIL2 lead to a more reasonable anomalous atmospheric heating associated with El Niño forcing in the tropical Pacifc， which further drives a better PSAC as the Rossby-wave response. The improved simulation of the PSAC further results in more reasonable wind anomalies over the East Asian continent.

A previous study of GAMIL suggested that improved simulation of the tropical mean climate is associated with an improved representation of the tropical climate anomalies related to El Niño events （Li， Wang， and Zhang 2014）. Our analysis fnds that the EAWM simulated by GAMIL2 is better than that by GAMIL1， suggesting a connection between the EAWM mean states and the interannual variation. Further work based on multi-model intercomparisons is needed to clarify the issue.

4. Summary

The climatology and interannual variability of the EAWM simulated by two versions of the IAP/LASG AGCM， GAMIL1 and GAMIL2， are evaluated in this study. The physical processes that dominate the model performance are analysed. The main fndings can be summarized as follows:

（1） The two versions of GAMIL reasonably simulate the climatological pattern of the winter monsoon over the East Asian continent， including the 2-m temperature， SLP， and V1000. The performance for the upper troposphere （500-hPa geopotential height and 300-hPa zonal wind） are better than that for precipitation， SLP， and V1000. GAMIL2 shows better performance than GAMIL1 in simulating the climatology of precipitation， SLP， and V1000.

（2） The interannual variability of the EAWM simulated by GAMIL2 is signifcantly improved compared to GAMIL1， as evidenced by change in correlation coefcient from 0.37 to 0.71. The improvement is also evident in the surface wind anomalies and the dipole rainfall pattern and its magnitude. The improved simulation of monsoon anomalies derives from the better simulation of the PSAC associated with El Niño forcing. The improved moist processes in GAMIL2， including stratiform condensation and evaporation， have led to a reasonable atmospheric heating associated with El Niño forcing in the tropical Pacifc， which further drives the PSAC as a Rossby-wave response.

Funding

This work was supported by the National Natural Science Foundation of China ［grant numbers 41330423 and 41420104006］.

Notes on contributors

Chen-xi Jin is PhD candidate at Institute of Atmospheric Physics， Chinese Academy of Sciences. His main research interests focus on East Asian winter monsoon. Recent publications include papers in Atmospheric and Oceanic Sciences Letters， and Chinese Journal of Atmospheric Sciences.

Tian-jun Zhou is a scientist in the LASG， Institute of Atmospheric Physics， Chinese Academy of Sciences. His research interests include coupled atmosphere—ocean modeling and climate dynamics， East Asian climate change and climate variability， air—sea interaction and Asian-Australian monsoon predictability. Recent publications include papers in Journal of Climate， Scientifc Reports，Journal of Geophysical Research — Atmospheres， Climate Dynamics and other journals.

References

Adler， R. F.， G. J. Huffman， A. Chang， R. Ferraro， P. P. Xie， J. Janowiak， B. Rudolf， et al. 2003. “The Version-2 Global Precipitation Climatology Project （GPCP） Monthly Precipitation Analysis （1979—Present）.” Journal of Hydrometeorology 4: 1147—1167.

Bellenger， H.， E. Guilyardi， J. Leloup， M. Lengaigne， and J. Vialard. 2013. “ENSO Representation in Climate Models: From CMIP3 to CMIP5.” Climate Dynamics 42: 1999—2018.

Chang， C. P.， and K. M. Lau. 1980. “Northeasterly Cold Surges and near-Equatorial Disturbances over the Winter MONEX Area during December 1974. Part II: Planetary-Scale Aspects.”Monthly Weather Review 108: 298—312.

Chang， C. P.， and K. M. Lau. 1982. “Short-Term Planetary-Scale Interaction over the Tropics and Midlatitudes during Northern Winter. Part I: Contrast between Active and Inactive Periods.” Monthly Weather Review 110: 933—946.

Chen， W.， K. Wei， L. Wang， and Q. Zhou. 2013. “Climate Variability and Mechanisms of the East Asian Winter Monsoon and the Impact from the Stratosphere.” ［In Chinese.］ Chinese Journal of Atmospheric Sciences 37 （2）: 425—438.

Chen， W.， S. Yang， and R. H. Huang. 2005. “Relationship between Stationary Planetary Wave Activity and the East Asian Winter Monsoon.” Journal of Geophysical Research 110: D14110. doi:10.1029/2004JD005669.

Ding， Y. H. 1994. Monsoon over China. Dordrecht: Kluwer Academic Publishers， 419pp.

Gong， H. N.， L. Wang， W. Chen， R. G. Wu， K. Wei， and X. F. Cui. 2014. “The Climatology and Interannual Variability of the East Asian Winter Monsoon in CMIP5 Models.” Journal of Climate 27: 1659—1678.

Gong， D. Y.， S. W. Wang， and J. H. Zhu. 2001. “East Asian Winter Monsoon and Arctic Oscillation.” Geophysical Research Letters 28: 2073—2076.

Huang， R. H.， J. L. Chen， L. Wang， and Z. D. Lin. 2012.“Characteristics， Processes， and Causes of the Spatio-Temporal Variabilities of the East Asian Monsoon System.”Advances in Atmospheric Sciences 29 （5）: 910—942.

Huang， R. H.， L. T. Zhou， and W. Chen. 2003. “The Progresses of Recent Studies on the Variabilities of the East Asian Monsoonand Their Causes.” Advances in Atmospheric Sciences 20 （1）: 55—69.

Ji， L. R.， S. Q. Sun， K. Arpe， and L. Bengtsson. 1997. “Model Study on the Interannual Variability of Asian Winter Monsoon and Its Infuence.” Advances in Atmospheric Sciences 14: 1—22.

Jiang， D. B.， Z. P. Tian， and X. M. Lang. 2015. “Reliability of Climate Models for China through the IPCC Third to Fifth Assessment Reports.” International Journal of Climatology. doi:10.1002/ joc.4406.

Jin， C. X.， and T. J. Zhou. 2014. “Analysis of the Interannual Variations of the East Asian Winter Monsoon Simulation by Four CMIP5 GCMs.” ［In Chinese.］ Chinese Journal of Atmospheric Sciences 38 （3）: 453—468.

Kanamitsu， M.， W. Ebisuzaki， J. Woollen， S.-K. Yang， J. J. Hnilo， M. Fiorino， and G. L. Potter. 2002. “NCEP—DOE AMIP-II Reanalysis（R-2）.” Bulletin of the American Meteorological Society 83: 1631—1643.

Li， L. J.， P. F. Lin， Y. Q. Yu， B. Wang， T. J. Zhou， L. Liu， and J. P. Liu. 2013a. “The Flexible Global Ocean-Atmosphere-Land System Model， Grid-Point Version 2: FGOALS-g2.” Advances in Atmospheric Sciences 30: 543—560.

Li， L. J.， B. Wang， L. Dong， L. Liu， S. Shen， N. Hu， and W. Q. Sun. 2013b. “Evaluation of Grid-Point Atmospheric Model of IAP LASG Version 2 （GAMIL2）.” Advances in Atmospheric Sciences 30: 855—867.

Li， L. J.， B. Wang， and G. J. Zhang. 2014.“The Role of Nonconvective Condensation Processes in Response of Surface Shortwave Cloud Radiative Forcing to El Niño Warming.” Journal of Climate 27: 6721—6736.

Li， L. J.， X. Xie， B. Wang， and L. Dong. 2012. “Evaluating the Performances of GAMIL1.0 and GAMIL2.0 during TWP-ICE with CAPT.” Atmospheric and Oceanic Science Letters 5: 38—42. Morrison， H.， and A. Gettelman. 2008. “A New Two-Moment Bulk Stratiform Cloud Microphysics Scheme in the Community Atmosphere Model， Version 3 （CAM3）. Part I: Description and Numerical Tests.” Journal of Climate 21: 3642—3659.

Rasch， P. J.， and J. E. Kristjánsson. 1998. “A Comparison of the CCM3 Model Climate Using Diagnosed and Predicted Condensate Parameterizations.” Journal of Climate 11: 1587—1614.

Saha， S.， S. Moorthi， H. L. Pan， X. R. Wu， J. D. Wang， S. Nadiga， and P. Tripp. 2010. “The NCEP Climate Forecast System Reanalysis.”Bulletin of the American Meteorological Society 91: 1015—1057. Shi， X. J.， B. Wang， X. H. Liu， M. H. Wang， L. L. Juan， and L. Dong. 2010. “Aerosol Indirect Efects on Warm Clouds in the Grid-Point Atmospheric Model of IAP LASG （GAMIL）.” Atmospheric and Oceanic Science Letters 3 （4）: 237—241.

Smith， T.， R. Reynolds， T. Peterson， and J. Lawrimore. 2008.“Improvements to NOAA's Historical Merged Land-Ocean Surface Temperature Analysis （1880—2006）.” Journal of Climate 21: 2283—2296.

Wang， L.， and W. Chen. 2010. “How Well Do Existing Indices Measure the Strength of the East Asian Winter Monsoon？”Advances in Atmospheric Sciences 27 （4）: 855—870.

Wang， B.， Z. W. Wu， C. P. Chang， J. Liu， J. P. Li， and T. J. Zhou. 2010.“Another Look at Interannual-to-Interdecadal Variations of the East Asian Winter Monsoon: The Northern and Southern Temperature Modes.” Journal of Climate 23: 1495—1512.

Wang， B.， R. G. Wu， and X. H. Fu. 2000. “Pacifc-East Asian Teleconnection: How Does ENSO Afect East Asian Climate？”Journal of Climate 13: 1517—1536.

Wei， K.， T. Xu， Z. C. Du， H. N. Gong， and B. H. Xie. 2014. “How Well Do the Current State-of-the-Art CMIP5 Models Characterise the Climatology of the East Asian Winter Monsoon？” Climate Dynamics 43: 1241—1255.

Zhang， G. J. 2002. “Convective Quasi-Equilibrium in Midlatitude Continental Environment and Its Efect on Convective Parameterization.” Journal of Geophysical Research 107: 4220. doi:10.1029/2001JD001005.

Zhang， G. J.， and N. A. McFarlane. 1995. “Sensitivity of Climate Simulations to the Parameterization of Cumulus Convection in the Canadian Climate Centre General Circulation Model.”Atmosphere Ocean 33: 407—446.

Zhang， R. H.， A. Sumi， and M. Kimoto. 1996. “Impact of El Nino on the East Asian Monsoon: A Diagnostic Study of the ‘86/87 and‘91/92 Events.” Journal of the Meteorological Society of Japan 74 （1）: 49—62.

Zhou， W.， J. C. L. Chan， W. Chen， J. Ling， J. Pinto， and Y. P. Shao. 2009. “Synoptic-Scale Controls of Persistent Low Temperature and Icy Weather over Southern China in January 2008.”Monthly Weather Review 137 （11）: 3978—3991.

Zhou， T. J.， X. L. Chen， L. Dong， B. Wu， W. M. Man， L. X. Zhang， R. P. Lin，J. C. Yao， F. F. Song， and C. B. Zhao. 2014a. “Chinese Contribution to CMIP5: An Overview of Five Chinese Models' Performances.”Journal of Meteorological Research 28 （4）: 481—509.

Zhou， T. J.， Y. Q. Yu， Y. M. Liu， and B. Wang. 2014b. Flexible Global Ocean-Atmosphere-Land System Model. Heidelberg: Springer，483 pp.

28 December 2015

CONTACT ZHOU Tian-Jun zhoutj@lasg.iap.ac.cn

© 2016 The Author（s）. Published by Taylor & Francis.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License （http://creativecommons.org/licenses/by/4.0/）， which permits unrestricted use， distribution， and reproduction in any medium， provided the original work is properly cited.