Advanced Search
Article Contents

A Hybrid Coupled Ocean-Atmosphere Model and Its Simulation of ENSO and Atmospheric Responses


doi: 10.1007/s00376-019-8197-8

  • A new hybrid coupled model (HCM) is presented in this study, which consists of an intermediate tropical Pacific Ocean model and a global atmospheric general circulation model. The ocean component is the intermediate ocean model (IOM) of the intermediate coupled model (ICM) used at the Institute of Oceanology, Chinese Academy of Sciences (IOCAS). The atmospheric component is ECHAM5, the fifth version of the Max Planck Institute for Meteorology atmospheric general circulation model. The HCM integrates its atmospheric and oceanic components by using an anomaly coupling strategy. A 100-year simulation has been made with the HCM and its simulation skills are evaluated, including the interannual variability of SST over the tropical Pacific and the ENSO-related responses of the global atmosphere. The model shows irregular occurrence of ENSO events with a spectral range between two and five years. The amplitude and lifetime of ENSO events and the annual phase-locking of SST anomalies are also reproduced realistically. Despite the slightly stronger variance of SST anomalies over the central Pacific than observed in the HCM, the patterns of atmospheric anomalies related to ENSO, such as sea level pressure, temperature and precipitation, are in broad agreement with observations. Therefore, this model can not only simulate the ENSO variability, but also reproduce the global atmospheric variability associated with ENSO, thereby providing a useful modeling tool for ENSO studies. Further model applications of ENSO modulations by ocean-atmosphere processes, and of ENSO-related climate prediction, are also discussed.
    摘要: ENSO是气候系统中最强的年际信号且具有广泛的气候影响, 所以, 对ENSO的模拟和预测具有重要的意义. 海气耦合模式是实现模拟和预测ENSO的主要工具, 其中, 简单海气耦合模式的动力过程简单, 一般为仅限于热带太平洋海域的距平模式, 计算量小, 但无法直接模拟和预测ENSO引起的全球范围的气候异常; 而全球海气耦合模式的物理过程更为复杂和完善, 可以直接模拟和预测热带外的气候异常, 但是一直受到“气候漂移”问题的困扰. 所以, 希望发展一个耦合模式, 既没有气候漂移, 又能够预测全球气候变化. 因此, 本文将热带太平洋区域海洋模式和全球大气环流模式耦合, 构建了一个新的混合型海气耦合模式(HCM), 其中, 海洋部分为中国科学院海洋研究所海气耦合模式(IOCAS ICM)的简单海洋模式, 大气部分为德国马普气象研究所的第五代大气环流模式(ECHAM5). 海洋和大气之间通过距平的动量通量和热量通量耦合. 利用HCM 100年模拟积分的结果, 对模式中ENSO及其全球大气响应的模拟进行了分析. 结果表明, 与IOCAS ICM相比, HCM中的ENSO事件具有显著的不规则性, 主要周期为2-5年振荡, 模拟的ENSO强度, 循环和季节锁相都与观测更为一致. 另外, HCM的大气分量表现出显著的年际变化, 冬, 夏季海平面气压, 温度, 降水和异常环流的分布与观测基本相符. 所以, 该模式不仅能够模拟ENSO及全球大气响应, 而且能够避免气候漂移, 计算量也远远小于复杂的全球海气耦合模式, 进而为ENSO模拟和预测研究提供了新的工具.
  • 加载中
  • An S.-I.,B. Wang, 2001: Mechanisms of locking of the El NiÑo and La NiÑa mature phases to boreal winter. J. Climate, 14, 2164-2176, https://doi.org/10.1175/1520-0442(2001)014<2164:MOLOTE>2.0.CO;2
    Ashok K.,S. K. Behera, S. A. Rao, H. Y. Weng, and T. Yamagata, 2007: El NiÑo Modoki and its possible teleconnection. J. Geophys. Res., 112, C11007, https://doi.org/10.1029/2006JC003798
    Barnston A. G.,M. K. Tippett, M. L. L'Heureux, S. H. Li, and D. G. DeWitt, 2012: Skill of real-time seasonal ENSO model predictions during 2002-11: Is our capability increasing? Bull. Amer. Meteor. Soc., 93, 631-651, https://doi.org/10.1175/BAMS-D-11-00111.1
    Barnston A. G.,M. K. Tippett, M. Ranganathan, and M. L. L'Heureux, 2017: Deterministic skill of ENSO predictions from the North American Multimodel Ensemble. Climate Dyn., https://doi.org/10.1007/s00382-017-3603-3
    Behringer D.,Y. Xue, 2004: Evaluation of the global ocean data assimilation system at NCEP: The Pacific Ocean. Eighth Symposium on Integrated Observing and Assimilation Systems for Atmosphere, Oceans, and Land Surface, AMS 84th Annual Meeting, Washington State Convention and Trade Center, Seattle, Washington, 11- 15.
    Bjerknes J.,1969: Atmospheric teleconnections from the equatorial Pacific. Mon. Wea. Rev., 97, 163, https://doi.org/10.1175/1520-0493(1969)097<0163:ATFTEP>2.3.CO;2
    Cane M. A.,S. E. Zebiak, and S. C. Dolan, 1986: Experimental forecasts of El NiÑo. Nature, 321, 827-832, https://doi.org/10.1038/321827a0
    Chen D. K.,S. E. Zebiak, A. J. Busalacchi, and M. A. Cane, 1995: An improved procedure for El NiÑo forecasting: Implications for predictability. Science, 269, 1699-1702, https://doi.org/10.1126/science.269.5231.1699
    Chen D. K.,M. A. Cane, A. Kaplan, S. E. Zebiak, and D. J. Huang, 2004: Predictability of El NiÑo over the past 148 years. Nature, 428(6984), 733-736, https://doi.org/10.1038/nature02439
    Chen, D. K.,Coauthors, 2015: Strong influence of westerly wind bursts on El NiÑo diversity. Nature Geoscience, 8, 339-345, https://doi.org/10.1038/ngeo2399
    Davey, M., Coauthors, 2002: STOIC: A study of coupled model climatology and variability in tropical ocean regions. Climate Dyn., 18, 403-420, https://doi.org/10.1007/s00382-001-0188-6
    Delworth, T. L.,Coauthors, 2006: GFDL's CM2 Global Coupled Climate Models. Part I: Formulation and Simulation Characteristics. J. Climate, 19(5), 643-674, https://doi.org/10.1175/JCLI3629.1
    Diaz H. F.,M. P. Hoerling, and J. K. Eischeid, 2001: ENSO variability, teleconnections and climate change. Int. J. Climatol., 21, 1845-1862, https://doi.org/10.1002/joc.631
    Duan W. S.,B. Tian, and H. Xu, 2014: Simulations of two types of El NiÑo events by an optimal forcing vector approach. Climate Dyn., 43, 1677-1692, https://doi.org/10.1007/s00382-013-1993-4
    Duan W. S.,X. C. Liu, K. Y. Zhu, and M. Mu, 2009: Exploring the initial errors that cause a significant "spring predictability barrier" for El NiÑo events. J. Geophys. Res., 114, C04022, https://doi.org/10.1029/2008JC004925
    Fedorov A. V.,2002: The response of the coupled tropical ocean-atmosphere to westerly wind bursts. Quart. J. Roy. Meteor. Soc., 128(579), 1-23, https://doi.org/10.1002/qj.200212857901
    Fedorov A. V.,S. N. Hu, M. Lengaigne, and E. Guilyardi, 2015: The impact of Westerly wind bursts and ocean initial state on the development, and diversity of El NiÑo events. Climate Dyn., 44(5-6), 1381-1401, https://doi.org/10.1007/s00382-014-2126-4
    Gao C.,R. H. Zhang, 2017: The roles of atmospheric wind and entrained water temperature (Te) in the second-year cooling of the 2010-12 La NiÑa event. Climate Dyn., 48, 597-617, https://doi.org/10.1007/s00382-016-3097-4
    Gualdi S.,A. Navarra, E. Guilyardi, and P. Delecluse, 2003: Assessment of the tropical Indo-Pacific climate in the SINTEX CGCM. Annals of Geophysics, 46(1), 1-26, https://doi.org/10.4401/ag-3385
    Guilyardi E.,A. Wittenberg, A. Fedorov, M. Collins, C. Z. Wang, A. Capotondi, G. J. van Oldenborgh, and T. Stockdale, 2009: Understanding El NiÑo in ocean-atmosphere general circulation models: Progress and challenges. Bull. Amer. Meteor. Soc., 90, 325-340, https://doi.org/10.1175/2008BAMS2387.1
    Gupta A. S.,N. C. Jourdain, J. N. Brown, and D. Monselesan, 2013: Climate drift in the CMIP5 models. J. Climate, 26(21), 8597-8615, https://doi.org/10.1175/JCLI-D-12-00521.1
    Hayashi M.,M. Watanab, 2017: ENSO complexity induced by state dependence of westerly wind events. J. Climate, 30(9), 3401-3420, https://doi.org/10.1175/JCLI-D-16-0406.1
    Hou Z. L.,J. P. Li, R. Q. Ding, J. Feng, and W. S. Duan, 2017: The application of nonlinear local Lyapunov vectors to the Zebiak-Cane model and their performance in ensemble prediction. Climate Dyn., 51, 283-304, https://doi.org/10.1007/s00382-017-3920-6
    Hu S. N.,A. V. Fedorov, M. Lengaigne, and E. Guilyardi, 2014: The impact of westerly wind bursts on the diversity and predictability of El NiÑo events: An ocean energetics perspective. Geophys. Res. Lett., 41, 4654-4663, https://doi.org/10.1002/2014GL059573
    Huang, B. Y.,Coauthors, 2014: Extended reconstructed sea surface temperature version 4 (ERSST.v4): Part I: Upgrades and intercomparisons. J. Climate, 28, 911-930, https://doi.org/10.1175/JCLI-D-14-00006.1
    Huffman G. J.,R. F. Adler, D. T. Bolvin, and G. J. Gu, 2009: Improving the global precipitation record: GPCP version 2.1. Geophys. Res. Lett., 36, L17808, https://doi.org/10.1029/2009GL040000
    Hurrell, J. W.,Coauthors, 2013: The community earth system model: A framework for collaborative research. Bull. Amer. Meteor. Soc., 94(9), 1339-1360, https://doi.org/10.1175/BAMS-D-12-00121.1
    Jin, E. K.,Coauthors, 2008: Current status of ENSO prediction skill in coupled ocean-atmosphere models. Climate Dyn., 31(6), 647-664, https://doi.org/10.1007/s00382-008-0397-3
    Jin F.-F.,S.-I. An, A. Timmermann, and J. X. Zhao, 2003: Strong El NiÑo events and nonlinear dynamical heating. Geophys. Res. Lett., 30, 1120, https://doi.org/10.1029/2002GL016356
    Kalnay, E., Coauthors, 1996: The NCEP/NCAR 40-year reanalysis project. Bull. Amer. Meteor. Soc., 77, 437, https://doi.org/10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2
    Kao H.-Y.,J.-Y. Yu, 2009: Contrasting eastern-pacific and central-pacific types of ENSO. J. Climate, 22(3), 615-632, https://doi.org/10.1175/2008JCLI2309.1
    Keenlyside N.,2001: Improved modeling of zonal currents and SST in the tropical Pacific. PhD dissertation, Monash University, 193 pp.
    Keenlyside N.,R. Kleeman, 2002: Annual cycle of equatorial zonal currents in the Pacific. J. Geophys. Res., 107, 3093, https://doi.org/10.1029/2000JC000711
    Kessler W. S.,2002: Is ENSO a cycle or a series of events? Geophys. Res. Lett., 29(23), 2125, https://doi.org/10.1029/2002gl015924
    Kug J. S.,F. F. Jin, and S.-I. An, 2009: Two types of El NiÑo events: Cold tongue El NiÑo and warm pool El NiÑo. J. Climate, 22(6), 1499-1515, https://doi.org/10.1175/2008JCLI2624.1
    Latif, M., Coauthors, 2001: ENSIP: The El NiÑo simulation intercomparison project. Climate Dyn., 18, 255-276, https://doi.org/10.1007/s003820100174
    Lengaigne M.,E. Guilyardi, J.-P. Boulanger, C. Menkes, P. Delecluse, P. Inness, J. Cole, and J. Slingo, 2004: Triggering of El NiÑo by westerly wind events in a coupled general circulation model. Climate Dyn., 23(6), 601-620, https://doi.org/10.1007/s00382-004-0457-2
    Levitus S.,1982: Climatological atlas of the world ocean. NOAA Prof. Paper 13, 173 pp.10.1029/EO064i049p00962-02d42441ccfe5870dbf336ba3ba886dd0ahttp%3A%2F%2Fonlinelibrary.wiley.com%2Fdoi%2F10.1029%2FEO064i049p00962-02%2Ffullhttp://onlinelibrary.wiley.com/doi/10.1029/EO064i049p00962-02/fullA project to objectively analyze historical ocean temperature, salinity, oxygen, and percent oxygen saturation data for the world ocean has recently been completed at the National Oceanic and Atmospheric Administration's (NOAA) Geophysical Fluid Dynamics Laboratory, Princeton, New Jersey. The results of the project are being made available through distribution of the Climatological Atlas of the World Ocean (NOAA Professional Paper No. 13), and through distribution of magnetic tapes containing the objective analyses.The sources of data used in the project were the Station Data, Mechanical Bathythermograph, and Expendable Bathythermograph files of the National Oceanographic Data Center (NODC) in Washington, D.C., updated through 1977 1978. The raw data were subjected to quality control procedures, averaged by one-degree squares, and then used as input to an objective analysis procedure that fills in one-degree squares containing no data and smooths the results. Due to the lack of synoptic observations for the world ocean, the historical data are composited by annual, seasonal, and (for temperature) monthly periods.
    Lian T.,D. K. Chen, Y. M. Tang, and Q. Y. Wu, 2014: Effects of westerly wind bursts on El NiÑo: A new perspective. Geophys. Res. Lett., 41, 3522-3527, https://doi.org/10.1002/2014GL059989
    Liu, W., Coauthors, 2014: Extended Reconstructed Sea Surface Temperature version 4 (ERSST.v4): Part II. Parametric and structural uncertainty estimations. J. Climate, 28, 931-951, https://doi.org/10.1175/JCLI-D-14-00007.1
    Lopez H.,B. P. Kirtman, 2013: Westerly wind bursts and the diversity of ENSO in CCSM3 and CCSM4. Geophys. Res. Lett., 40, 4722-4727, https://doi.org/10.1002/grl.50913
    Lopez H.,B. P. Kirtman, 2014: WWBs, ENSO predictability, the spring barrier and extreme events. J. Geophys. Res., 119, 10 114-10 138, https://doi.org/10.1002/2014JD021908
    Luo J.-J.,S. Masson, E. Roeckner, G. Madec, and T. Yamagata, 2005: Reducing climatology bias in an ocean-atmosphere CGCM with improved coupling physics. J. Climate, 18(13), 2344-2360, https://doi.org/10.1175/JCLI3404.1
    McCreary J. P.,1981: A linear stratified ocean model of the equatorial undercurrent. Philos. Trans. Roy. Soc. London, 298, 603-635, https://doi.org/10.1098/rsta.1981.0002
    McPhaden M. J.,1999: Genesis and evolution of the 1997-98 El NiÑo. Science, 283(5404), 950-954, https://doi.org/10.1126/science.283.5404.950
    McPhaden M. J.,X. R. Yu, 1999: Equatorial waves and the 1997-98 El NiÑo. Geophys. Res. Lett., 26(19), 2961-2964, https://doi.org/10.1029/1999GL004901
    McPhaden M. J.,S. E. Zebiak, and M. H. Glantz, 2006: ENSO as an integrating concept in Earth science. Science, 314, 1740-1745, https://doi.org/10.1126/science.1132588
    Mu M.,W. S. Duan, and B. Wang, 2007a: Season-dependent dynamics of nonlinear optimal error growth and El NiÑo-Southern Oscillation predictability in a theoretical model. J. Geophys. Res., 112, D10113, https://doi.org/10.1029/2005JD006981
    Mu M.,H. Xu, and W. S. Duan, 2007b: A kind of initial errors related to "spring predictability barrier" for El NiÑo events in Zebiak-Cane model. Geophys. Res. Lett., 34, L03709, https://doi.org/10.1029/2006GL027412
    Philand er, S. G. H., 1990: El NiÑo, La NiÑa, and the Southern Oscillation. Academic Press, 293 pp.
    Reynolds R. W.,T. M. Smith, 1995: A high-resolution global sea surface temperature climatology. J. Climate, 8, 1571-1583, https://doi.org/10.1175/1520-0442(1995)008<1571:AHRGSS>2.0.CO;2
    Richter I.,2015: Climate model biases in the eastern tropical oceans: Causes, impacts and ways forward. WIREs Climate Change, 6(3), 345-358, https://doi.org/10.1002/wcc.338
    Roeckner, E., Coauthors, 2003: The atmospheric general circulation model ECHAM 5. Part I: Model description. Max Planck Institute for Meteorology Rep., 349, 127 pp.
    Ropelewski C. F.,M. S. Halpert, 1987: Global and regional scale precipitation patterns associated with the El NiÑo/ Southern Oscillation. Mon. Wea. Rev., 115(8), 1606-1626, https://doi.org/10.1175/1520-0493(1987)115<1606:GARSPP>2.0.CO;2
    Stan C.,D. M. Straus, J. S. Frederiksen, H. Lin, E. D. Maloney, and C. Schumacher, 2017: Review of tropical-extratropical teleconnections on intraseasonal time scales. Rev. Geophys., 55, 902-937, https://doi.org/10.1002/2016RG000538
    Stricherz., J. N., D. M. Legler, and J. J. O'Brien, 1995: TOGA Pseudo-Stress Atlas 1985-1994. Volume II, Pacific Ocean, Florida State University, 173 pp.
    Tao L.-J.,R.-H. Zhang, and C. Gao, 2017: Initial error-induced optimal perturbations in ENSO predictions, as derived from an intermediate coupled model. Adv. Atmos. Sci., 34(6), 791-803, https://doi.org/10.1007/s00376-017-6266-4
    Tao L.-J.,C. Gao, and R.-H. Zhang, 2018: ENSO predictions in an intermediate coupled model influenced by removing initial condition errors in sensitive areas: A target observation perspective. Adv. Atmos. Sci., 35(7), 853-867, https://doi.org/10.1007/s00376-017-7138-7
    Timmermann, A., Coauthors, 2018: El NiÑo-southern oscillation complexity. Nature, 559, 535-545, https://doi.org/10.1038/s41586-018-0252-6
    Trenberth K. E.,G. W. Branstator, D. Karoly, A. Kumar, N. C. Lau, and C. Ropelewski, 1998: Progress during TOGA in understanding and modeling global teleconnections associated with tropical sea surface temperatures. J. Geophys. Res., 103, 14 291-14 324, https://doi.org/10.1029/97JC01444
    Tziperman E.,M. A. Cane, S. E. Zebiak, Y. Xue, and B. Blumenthal, 1998: Locking of El NiÑo's peak time to the end of the calendar year in the delayed oscillator picture of ENSO. J. Climate, 11, 2191-2199, https://doi.org/10.1175/1520-0442(1998)011<2191:LOENOS>2.0.CO;2
    Vega-Westhoff B.,R. L. Sriver, 2017: Analysis of ENSO's response to unforced variability and anthropogenic forcing using CESM. Scientific Reports, 7, 18047, https://doi.org/10.1038/s41598-017-18459-8
    Wang C. Z.,L. P. Zhang, S.-K. Lee, L. X. Wu, and C. R. Mechoso, 2014: A global perspective on CMIP5 climate model biases. Nat. Clim. Change, 4(3), 201-205, https://doi.org/10.1038/nclimate2118
    Webster P. J.,S. Yang, 1992: Monsoon and ENSO: Selectively interactive systems. Quart. J. Roy. Meteor. Soc., 118(507), 877-926, https://doi.org/10.1002/qj.49711850705
    Wittenberg A. T.,2009: Are historical records sufficient to constrain ENSO simulations? Geophys. Res. Lett., 36, L12702, https://doi.org/10.1029/2009GL038710
    Zebiak S. E.,M. A. Cane, 1987: A model El NiÑo-Southern Oscillation. Mon. Wea. Rev., 115(10), 2262, https://doi.org/10.1175/1520-0493(1987)115<2262:AMENO>2.0.CO;2
    Zhang R.-H.,2015: A hybrid coupled model for the Pacific ocean-atmosphere system. Part I: Description and basic performance. Adv. Atmos. Sci., 32(3), 301-318, https://doi.org/10.1007/s00376-014-3266-5
    Zhang R.-H.,C. Gao, 2016a: Role of subsurface entrainment temperature (Te) in the onset of El NiÑo events, as represented in an intermediate coupled model. Climate Dyn., 46, 1417-1435, https://doi.org/10.1007/s00382-015-2655-5
    Zhang R.-H.,C. Gao, 2016b: The IOCAS intermediate coupled model (IOCAS ICM) and its real-time predictions of the 2015-2016 El NiÑo event. Science Bulletin, 61(13), 1061-1070, https://doi.org/10.1007/s11434-016-1064-4
    Zhang R.-H.,A. J. Busalacchi, and D. G. DeWitt, 2008: The roles of atmospheric stochastic forcing (SF) and oceanic entrainment temperature (Te) in decadal modulation of ENSO. J. Climate, 21, 674-704, https://doi.org/10.1175/2007JCLI1665.1
    Zhang R.-H.,S. E. Zebiak, R. Kleeman, and N. Keenlyside, 2003: A new intermediate coupled model for El NiÑo simulation and prediction. Geophys. Res. Lett., 30, 2012, https://doi.org/10.1029/2003GL018010
    Zhang R.-H.,F. Zheng, J. Zhu, and Z. G. Wang, 2013: A successful real-time forecast of the 2010-11 La NiÑa event. Scientific Reports, 3, 1108, https://doi.org/10.1038/srep01108
    Zhang R.-H.,A. J. Busalacchi, R. G. Murtugudde, E. C. Hackert, and J. Ballabrera-Poy, 2004: A new approach to improved SST anomaly simulations using altimeter data: Parameterizing entrainment temperature from sea level. Geophys. Res. Lett., 31, L10304, https://doi.org/10.1029/2003GL019237
    Zhang R.-H.,S. E. Zebiak, R. Kleeman, and N. Keenlyside, 2005b: Retrospective El NiÑo forecasts using an improved intermediate coupled model. Mon. Wea. Rev., 133, 2777-2802, https://doi.org/10.1175/MWR3000.1
    Zhang R.-H.,R. Kleeman, S. E. Zebiak, N. Keenlyside, and S. Raynaud, 2005a: An empirical parameterization of subsurface entrainment temperature for improved SST anomaly simulations in an intermediate ocean model. J. Climate, 18, 350-371, https://doi.org/10.1175/JCLI-3271.1
    Zheng F.,J. Zhu, R.-H. Zhang, and G.-Q. Zhou, 2006: Ensemble hindcasts of SST anomalies in the tropical Pacific using an intermediate coupled model. Geophys. Res. Lett., 33, L19604, https://doi.org/10.1029/2006GL026994
    Zhu Y. C.,R.-H. Zhang, 2017: Scaling wind stirring effects in an oceanic bulk mixed layer model with application to an OGCM of the tropical Pacific. Climate Dyn., 51, 1927-1946, https://doi.org/10.1007/s00382-017-3990-5
    Zhu Y. C.,R.-H. Zhang, 2018: An Argo-derived background diffusivity parameterization for improved ocean simulations in the tropical pacific. Geophys. Res. Lett., 45, 1509-1517, https://doi.org/10.1002/2017GL076269
  • [1] Junya HU, Hongna WANG, Chuan GAO, Rong-Hua ZHANG, 2024: Different El Niño Flavors and Associated Atmospheric Teleconnections as Simulated in a Hybrid Coupled Model, ADVANCES IN ATMOSPHERIC SCIENCES, 41, 864-880.  doi: 10.1007/s00376-023-3082-x
    [2] ZHANG Rong-Hua, 2015: A Hybrid Coupled Model for the Pacific Ocean-Atmosphere System. Part I: Description and Basic Performance, ADVANCES IN ATMOSPHERIC SCIENCES, 32, 301-318.  doi: 10.1007/s00376-014-3266-5
    [3] Ben TIAN, Wansuo DUAN, 2016: Comparison of Constant and Time-variant Optimal Forcing Approaches in El Niño Simulations by Using the Zebiak-Cane Model, ADVANCES IN ATMOSPHERIC SCIENCES, 33, 685-694.  doi: 10.1007/s00376-015-5174-8
    [4] Guidi ZHOU, Xuhua CHENG, 2022: Impacts of Oceanic Fronts and Eddies in the Kuroshio-Oyashio Extension Region on the Atmospheric General Circulation and Storm Track, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 22-54.  doi: 10.1007/s00376-021-0408-4
    [5] ZHU Jieshun, SUN Zhaobo, ZHOU Guangqing, 2007: A Note on the Role of Meridional Wind Stress Anomalies and Heat Flux in ENSO Simulations, ADVANCES IN ATMOSPHERIC SCIENCES, 24, 729-738.  doi: 10.1007/s00376-007-0729-y
    [6] Wu Aiming, Ni Yunqi, 1999: A Hybrid Coupled Ocean-Atmosphere Model and ENSO Prediction Study, ADVANCES IN ATMOSPHERIC SCIENCES, 16, 405-418.  doi: 10.1007/s00376-999-0019-y
    [7] Zhou Tianjun, Yu Rucong, Li Zhaoxin, 2002: ENSO-Dependent and ENSO-Independent Variability over the Mid-Latitude North Pacific: Observation and Air-Sea Coupled Model Simulation, ADVANCES IN ATMOSPHERIC SCIENCES, 19, 1127-1147.  doi: 10.1007/s00376-002-0070-4
    [8] Ni Yunqi, Zou Li, Wu Aiming, 1996: Response Process of Ocean to Atmospheric Forcing and Optimal Response Frequency in the CZ Ocean Model, ADVANCES IN ATMOSPHERIC SCIENCES, 13, 43-58.  doi: 10.1007/BF02657027
    [9] HUANG Ping, HUANG Ronghui, 2009: Delayed Atmospheric Temperature Response to ENSO SST: Role of High SST and the Western Pacific, ADVANCES IN ATMOSPHERIC SCIENCES, 26, 343-351.  doi: 10.1007/s00376-009-0343-2
    [10] Chang LIU, Shaoqing ZHANG, Shan LI, Zhengyu LIU, 2017: Impact of the Time Scale of Model Sensitivity Response on Coupled Model Parameter Estimation, ADVANCES IN ATMOSPHERIC SCIENCES, 34, 1346-1357.  doi: 10.1007/s00376-017-6272-6
    [11] Lu ZHOU, Rong-Hua ZHANG, 2022: A Hybrid Neural Network Model for ENSO Prediction in Combination with Principal Oscillation Pattern Analyses, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 889-902.  doi: 10.1007/s00376-021-1368-4
    [12] Zhou Guangqing, Zeng Qingcun, 2001: Predictions of ENSO with a Coupled Atmosphere-Ocean General Circulation Model, ADVANCES IN ATMOSPHERIC SCIENCES, 18, 587-603.  doi: 10.1007/s00376-001-0047-8
    [13] Jiangbo JIN, Xiao DONG, Juanxiong HE, Yi YU, Hailong LIU, Minghua ZHANG, Qingcun ZENG, He ZHANG, Xin GAO, Guangqing ZHOU, Yaqi WANG, 2022: Ocean Response to a Climate Change Heat-Flux Perturbation in an Ocean Model and Its Corresponding Coupled Model, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 55-66.  doi: 10.1007/s00376-021-1167-y
    [14] Chuan GAO, Rong-Hua ZHANG, Xinrong WU, Jichang SUN, 2018: Idealized Experiments for Optimizing Model Parameters Using a 4D-Variational Method in an Intermediate Coupled Model of ENSO, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 410-422.  doi: 10.1007/s00376-017-7109-z
    [15] Zhang Renjian, Wang Mingxing, Zeng Qingcun, 2000: Global Two-Dimensional Chemistry Model and Simulation of Atmospheric Chemical Composition, ADVANCES IN ATMOSPHERIC SCIENCES, 17, 72-82.  doi: 10.1007/s00376-000-0044-3
    [16] Jiawen ZHU, Juanxiong HE, Duoying JI, Yangchun LI, He ZHANG, Minghua ZHANG, Xiaodong ZENG, Kece FEI, Jiangbo JIN, 2024: CAS-ESM2.0 Successfully Reproduces Historical Atmospheric CO2 in a Coupled Carbon−Climate Simulation, ADVANCES IN ATMOSPHERIC SCIENCES, 41, 572-580.  doi: 10.1007/s00376-023-3172-9
    [17] Bin MU, Juhui REN, Shijin YUAN, Rong-Hua ZHANG, Lei CHEN, Chuan GAO, 2019: The Optimal Precursors for ENSO Events Depicted Using the Gradient-definition-based Method in an Intermediate Coupled Model, ADVANCES IN ATMOSPHERIC SCIENCES, 36, 1381-1392.  doi: 10.1007/s00376-019-9040-y
    [18] Xianghui FANG, Fei ZHENG, 2018: Simulating Eastern- and Central-Pacific Type ENSO Using a Simple Coupled Model, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 671-681.  doi: 10.1007/s00376-017-7209-9
    [19] Chuan GAO, Xinrong WU, Rong-Hua ZHANG, 2016: Testing a Four-Dimensional Variational Data Assimilation Method Using an Improved Intermediate Coupled Model for ENSO Analysis and Prediction, ADVANCES IN ATMOSPHERIC SCIENCES, 33, 875-888.  doi: 10.1007/s00376-016-5249-1
    [20] Ling-Jiang TAO, Rong-Hua ZHANG, Chuan GAO, 2017: Initial Error-induced Optimal Perturbations in ENSO Predictions, as Derived from an Intermediate Coupled Model, ADVANCES IN ATMOSPHERIC SCIENCES, 34, 791-803.  doi: 10.1007/s00376-017-6266-4

Get Citation+

Export:  

Share Article

Manuscript History

Manuscript received: 27 September 2018
Manuscript revised: 23 January 2019
Manuscript accepted: 25 February 2019
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

A Hybrid Coupled Ocean-Atmosphere Model and Its Simulation of ENSO and Atmospheric Responses

    Corresponding author: Rong-Hua ZHANG, rzhang@qdio.ac.cn
  • 1. Chinese Academy of Sciences Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
  • 2. Pilot National Laboratory for Marine Science and Technology, Qingdao 266237, China
  • 3. Center for Ocean Mega-Science, Chinese Academy of Sciences, 7 Nanhai Road, Qingdao 266071, China
  • 4. University of Chinese Academy of Sciences, Beijing 10029, China

Abstract: A new hybrid coupled model (HCM) is presented in this study, which consists of an intermediate tropical Pacific Ocean model and a global atmospheric general circulation model. The ocean component is the intermediate ocean model (IOM) of the intermediate coupled model (ICM) used at the Institute of Oceanology, Chinese Academy of Sciences (IOCAS). The atmospheric component is ECHAM5, the fifth version of the Max Planck Institute for Meteorology atmospheric general circulation model. The HCM integrates its atmospheric and oceanic components by using an anomaly coupling strategy. A 100-year simulation has been made with the HCM and its simulation skills are evaluated, including the interannual variability of SST over the tropical Pacific and the ENSO-related responses of the global atmosphere. The model shows irregular occurrence of ENSO events with a spectral range between two and five years. The amplitude and lifetime of ENSO events and the annual phase-locking of SST anomalies are also reproduced realistically. Despite the slightly stronger variance of SST anomalies over the central Pacific than observed in the HCM, the patterns of atmospheric anomalies related to ENSO, such as sea level pressure, temperature and precipitation, are in broad agreement with observations. Therefore, this model can not only simulate the ENSO variability, but also reproduce the global atmospheric variability associated with ENSO, thereby providing a useful modeling tool for ENSO studies. Further model applications of ENSO modulations by ocean-atmosphere processes, and of ENSO-related climate prediction, are also discussed.

摘要: ENSO是气候系统中最强的年际信号且具有广泛的气候影响, 所以, 对ENSO的模拟和预测具有重要的意义. 海气耦合模式是实现模拟和预测ENSO的主要工具, 其中, 简单海气耦合模式的动力过程简单, 一般为仅限于热带太平洋海域的距平模式, 计算量小, 但无法直接模拟和预测ENSO引起的全球范围的气候异常; 而全球海气耦合模式的物理过程更为复杂和完善, 可以直接模拟和预测热带外的气候异常, 但是一直受到“气候漂移”问题的困扰. 所以, 希望发展一个耦合模式, 既没有气候漂移, 又能够预测全球气候变化. 因此, 本文将热带太平洋区域海洋模式和全球大气环流模式耦合, 构建了一个新的混合型海气耦合模式(HCM), 其中, 海洋部分为中国科学院海洋研究所海气耦合模式(IOCAS ICM)的简单海洋模式, 大气部分为德国马普气象研究所的第五代大气环流模式(ECHAM5). 海洋和大气之间通过距平的动量通量和热量通量耦合. 利用HCM 100年模拟积分的结果, 对模式中ENSO及其全球大气响应的模拟进行了分析. 结果表明, 与IOCAS ICM相比, HCM中的ENSO事件具有显著的不规则性, 主要周期为2-5年振荡, 模拟的ENSO强度, 循环和季节锁相都与观测更为一致. 另外, HCM的大气分量表现出显著的年际变化, 冬, 夏季海平面气压, 温度, 降水和异常环流的分布与观测基本相符. 所以, 该模式不仅能够模拟ENSO及全球大气响应, 而且能够避免气候漂移, 计算量也远远小于复杂的全球海气耦合模式, 进而为ENSO模拟和预测研究提供了新的工具.

1. Introduction
  • The ENSO phenomenon is the dominant source of interannual variability on Earth and can cause severe climatic and socioeconomic effects around the globe (Ropelewski and Halpert, 1987; Philander, 1990; Diaz et al., 2001; McPhaden et al., 2006; Stan et al., 2017). Hence, ENSO prediction has received widespread attention and a great deal of modeling efforts have been made to simulate and predict ENSO-related anomalies in the tropical Pacific.

    Coupled ocean-atmosphere models are major tools to simulate and predict ENSO, since ENSO arises from the air-sea interactions in the tropical Pacific Ocean (Bjerknes, 1969). In the past several decades, various coupled models with different degrees of complexity have been developed for ENSO studies, including intermediate coupled models (ICMs), hybrid coupled models (HCMs), and comprehensive coupled general circulation models (CGCMs). In recent years, CGCMs have progressed remarkably (Gualdi et al., 2003; Delworth et al., 2006; Guilyardi et al., 2009; Hurrell et al., 2013), and most current CGCMs are capable of simulating a good resemblance of the observed ENSO and associated variability. However, many CGCMs still suffer from the problem of climate drift, with considerable errors in reproducing the climatology of the equatorial Pacific, such as a too-cold cold tongue and warm SST biases in the eastern boundaries of the Pacific and Atlantic basins (Latif et al., 2001; Davey et al., 2002; Luo et al., 2005; Gupta et al., 2013; Wang et al., 2014; Richter, 2015; Zhu and Zhang, 2017, Zhu and Zhang, 2018). Among the aforementioned models, ICMs are efficient in computing, and show certain advantages for ENSO simulation and prediction. For example, compared to a complex CGCM, they can lead to ways in which coupling mechanisms are more easily understood. Furthermore, ICMs are free from the climate drift problem because they are normally designed to be anomaly models. In addition, owing to their low computational cost compared to CGCMs, ICMs allow a large number of experiments to be performed, feasibly and affordably. Famously, the Zebiak-Cane model was the first dynamical coupled model to successfully predict ENSO events (Cane et al., 1986). It is an intermediate anomaly model that describes the evolution of anomalies with respect to a prescribed seasonally varying background state. Since 1986, the Zebiak-Cane model has been used to conduct realistic ENSO simulations, and is widely used in dynamics and predictability studies of ENSO (Zebiak and Cane, 1987; Chen et al., 1995, Chen et al., 2004; Mu et al., 2007a, b; Duan et al., 2014; Hou et al., 2017).

    (Zhang et al., 2003) developed an improved ICM on the basis of an intermediate ocean model (IOM) designed by (Keenlyside and Kleeman, 2002). This ICM has been used at the Institute of Oceanology, Chinese Academy of Sciences (IOCAS), and is referred to as IOCAS ICM. It is an anomaly model, including an IOM with an embedded nonlinear SST anomaly model describing the thermodynamics for the surface mixed layer, and a statistical atmosphere model. Being constructed by a singular value decomposition (SVD) analysis technique, the atmosphere model determines the wind stress interannual anomaly from the ocean model SST anomaly. The SVD analysis is conducted to derive the relationship between the historical wind stress anomalies and SST anomalies, in which the wind stress data are the ensemble mean of 24-member ECHAM4.5 simulations and hence filter out unrelated atmospheric noise. Another striking feature of the ICM is an explicit parameterization of the subsurface entrainment temperature (Te) in terms of the thermocline variability based on SVD analysis (Zhang and Gao, 2016a, b). IOCAS ICM not only realistically reproduces the interannual variability associated with ENSO, but also shows high prediction skill over the central equatorial Pacific. Therefore, it has been widely used for ENSO modeling and predictability studies (Zhang et al., 2005a, b, 2008, 2013; Zheng et al., 2006; Gao et al., 2017; Tao et al., 2017, Tao et al., 2018). In particular, it has been routinely used for ENSO real-time predictions since 2003, and the model results are available at the website of the International Research Institute for Climate and Society, Columbia University (http://iri.columbia.edu/climate/ENSO/currentinfo/update.html).

    However, ENSO events in IOCAS ICM are much more regular in terms of their spatiotemporal evolution compared to observations. This is mostly due to the lack of modulation from the stochastic wind forcing, which is excluded by the statistical atmosphere model of the ICM since it only represents the interannual wind response to SST anomalies (Zhang2003,Zhang2008). Actually, an increasing number of studies are highlighting the influence of high-frequency synoptic-scale atmospheric variations on the complexity and predictability of ENSO. For example, (Zhang et al., 2008) indicated that atmospheric stochastic forcing can cause significant irregularity of ENSO and modulate the ENSO period and amplitude randomly. Westerly wind bursts (WWBs), a form of tropical weather noise usually occurring in the western and central Pacific, have been argued to play an important role in modulating the strength, asymmetry, irregularity or timing of ENSO events (McPhaden, 1999; McPhaden and Yu, 1999; Fedorov, 2002; Kessler, 2002; Lengaigne et al., 2004). In particular, (Fedorov et al., 2015) demonstrated that WWBs, when acting on different ocean initial states, can result in different flavors of El NiÑo——namely, central Pacific (CP) and eastern Pacific (EP) El NiÑo——thereby contributing to ENSO diversity (Lopez and Kirtman, 2013; Lian et al., 2014; Hu et al., 2014; Chen et al., 2015; Hayashi and Watanab, 2017). (Lopez and Kirtman, 2014) indicated WWB activities can lead to a rapid spring decline in the signal-to-noise ratio of a coupled system that induces a spring predictability barrier of ENSO prediction (Webster and Yang, 1992; Mu et al., 2007a, b; Duan et al., 2009). For the coupled internal variability, (Wittenberg, 2009) demonstrated that there was a long-term modulation of ENSO within a 2000-year simulation of GFDL CM2.1. (Vega-Westhoff and Sriver, 2017) further showed that the internal variations within CESM dominate the SST changes in the eastern Pacific compared to the forcing of climate change. It can be seen that a good representation of atmospheric internal variability and the coupled internal variability in the model is important for simulating ENSO characteristics and improving ENSO predictions. In addition, although ICMs do not suffer from climate drift, they cannot directly simulate and predict climate anomalies in the extratropics associated with ENSO-related SST forcing, whereas this can be achieved by adopting a complete AGCM. Therefore, we set out to improve the atmospheric component of IOCAS ICM for ENSO simulation and prediction, despite the model already being competitive in terms of ENSO prediction compared to many CGCMs.

    In this study, we construct a new coupled model for ENSO studies by coupling an AGCM (namely, ECHAM5) with the ocean component of IOCAS ICM (i.e., the IOM). Typically and previously, an HCM connects an OGCM to a simplified atmosphere model (e.g., Zhang, 2015). Here, we couple a simple ocean model with a complex AGCM. Despite this difference, the general hybrid concept (i.e., including one complex component and one simple component) is similar, and thus the coupled model presented here is still referred to as an HCM. Such an HCM is computationally efficient to run compared to CGCMs, and yet physically realistic enough to depict the atmosphere, allowing the investigation of ENSO modulations from the inherent variability of the atmosphere and global atmospheric response to ENSO, which are missing in any type of ICM. The main purpose of this paper is to document and validate the HCM simulations in terms of ENSO variability and global atmospheric responses to ENSO.

    To facilitate comparison, the framework of IOCAS ICM is first introduced in section 2, followed by descriptions of the atmospheric component model and the coupling procedure of the HCM, as well as the datasets used to validate the model. Section 3 presents the simulated interannual variability of SST in the tropical Pacific associated with ENSO within the HCM, and a comparison between the HCM and IOCAS ICM is also discussed in this section. In section 4, the global atmospheric response to ENSO is examined. Finally, a summary and discussion are given in section 5.

2. Model descriptions and datasets
  • IOCAS ICM includes an IOM, an SST anomaly model with an empirical parameterization of Te, and a statistical atmosphere model in the tropical Pacific Ocean.

    The dynamical IOM was originally designed by (Keenlyside, 2001) and (Keenlyside and Kleeman, 2002), and consists of a linear and a nonlinear component. The linear component is extended from the (McCreary, 1981) baroclinic modal model with a horizontally varying background stratification. The first 10 baroclinic modes are resolved in the vertical layers, plus two surface layers governed by Ekman dynamics to simulate the combined effect of high-order baroclinic modes. The nonlinear component, a simplified model derived from the residual nonlinear momentum equations, is incorporated within the two surface layers to provide corrections to the linear component where the nonlinearity cannot be ignored.

    An SST anomaly model is embedded into the ocean dynamical framework to describe the evolution of interannual temperature anomalies in the surface mixed layer. The governing equation includes ocean horizontal advection and entrainment by both specified mean currents and model simulated anomalous currents. (Zhang et al., 2003) demonstrated that the performance of SST simulations in the equatorial Pacific is significantly affected by the parameterization of Te. It has been shown that the interannual variability of the sea level and Te are closely correlated in the tropical Pacific (Zhang et al., 2004). Thus, an empirical Te submodel is constructed from the historical data of the sea level and Te anomalies by using the SVD method (Zhang et al., 2005a). The surface heat flux anomaly is negatively proportional to the local SST anomaly, with a thermal damping coefficient of (100 d)-1.

    A statistical atmospheric model is also constructed based on the SVD analysis. The SVD is used to determine the relationship between the historical wind stress and SST anomaly, and thus the wind stress and SST anomaly fields are specifically related. As such, according to the constructed wind stress model, the interannual wind response can be calculated given an SST anomaly.

    IOCAS ICM spans the tropical Pacific and Atlantic oceans (only the Pacific basin is considered in this work). Its domain covers (33.5°S-33.5°N, 124°-30°E), with a realistic representation of the continents. The model has a zonal grid with 2° spacing and a meridional grid stretching from 0.5° within 10° of the equator to 3° at the meridional northern and southern boundaries. Vertically, the ocean is assumed to be flat-bottomed with a depth of 5500 m. The linear component has 33 levels and 8 levels are in the upper 125 m. The two surface layers, in which nonlinear effects and high-order baroclinic modes are simulated, span the upper 125 m and are divided by a surface mixed layer whose depth is determined by a stability criterion from the annual mean temperature and salinity data in (Levitus, 1982). The dynamical ocean model and the SSTA model have the same grids. The model time step is 4800 s. The model's climatological fields include the SST of (Reynolds and Smith, 1995), model currents generated using the Florida State University wind stress (Stricherz et al., 1995), and thermocline depth constructed from (Levitus, 1982). The climatological fields are updated once monthly. More detailed descriptions of IOCAS ICM are given by (Zhang and Gao, 2016b). A 50-year control run of the ICM is used for comparison with the HCM constructed in this study.

  • The atmosphere model used in this study is ECHAM5, which is a global spectral model based on the primitive equations. Prognostic variables consist of vorticity, divergence, temperature, surface pressure, cloud water and water vapor. The horizontal spectral resolution of ECHAM5 used in this work is T63 (1.875°× 1.875°), with 19 vertical hybrid levels up to a pressure level of 10 hPa. The model employs a semi-implicit leapfrog time-stepping scheme with a weak time filter to inhibit the spurious computational modes. The model time step is 1200 s for dynamics and physics, and the radiation is calculated at 2-h intervals. A more detailed model description of ECHAM5 is given by (Roeckner et al., 2003).

    Before being coupled to the ocean component, two experiments are carried out with ECHAM5 alone, in which the SST boundary data are from AMIP II. First, ECHAM5 is run for 50 years, forced by the climatological monthly mean of SST, from which the climatology of ECHAM5 can be obtained and is used to calculate the anomalous coupling flux in the HCM. Second, ECHAM5 is forced by monthly historical SST from 1956 to 2000 in the tropical Pacific Ocean and the SST beyond the tropical Pacific is set to be the climatological monthly mean. Our analysis focuses on the period between 1961 and 2000. This experiment, together with observations, is used to evaluate the global atmospheric simulations of the HCM.

  • In this work, an HCM is constructed by coupling the ocean component of IOCAS ICM (i.e., the IOM) to ECHAM5. The atmosphere and ocean models are only coupled in the tropical Pacific Ocean. Beyond the active coupling regions, the underlying SST of the atmosphere is specified as the climatological monthly mean. To maintain the continuity of the boundary forcing, sponge layers are introduced at the northern and southern boundaries of the tropical Pacific, acting to relaxing SSTs to the climatological monthly mean. The coupling frequency between the atmosphere and ocean model is once per day. Given an SST, the atmosphere produces a total wind stress field and a net surface heat flux field, which is the sum of the solar radiation, longwave radiation, and sensible and latent heat fluxes. By subtracting the atmospheric climatology taken from the uncoupled simulation mentioned in section 2.2, the wind stress and heat flux anomalies are calculated and passed to force the ocean model. The ocean sends daily mean SST anomalies, superimposed on the climatological mean, back to the atmosphere. In order to maintain the radiative-convective equilibrium, the total SST is limited to be no greater than 30°C (Jin et al., 2003). Before being transferred to force the ocean model, the magnitude of wind stress anomalies is adjusted by multiplying a scalar parameter (ατ). This parameter is called the relative coupling coefficient and represents the strength of the interannual wind forcing on the ocean. Several tuning experiments are performed with different values of ατ to examine the coupled behavior in the HCM. It is found that taking ατ=0.8 can produce a reasonable interannual variability in the tropical Pacific. Because the model grids are different in the atmosphere and ocean, the exchanged variables are interpolated during the transfer process. The initial atmospheric field is a restart file from the previous 50-year integration. The initial condition of the oceanic component is the steady state after a long-term integration of the ocean model. The coupled model is integrated for 200 years and the last 100 years are used in the following analyses.

  • Observational and reanalysis datasets used for evaluating the model simulation include the following: SST from ERSST.v4 from 1951 to 2000 (Huang et al., 2014; Liu et al., 2014); SSH and zonal wind stress from GODAS from 1981 to 2017 (Behringer and Xue, 2004); SLP, 2-m temperature, and 500-hPa geopotential height from the NCEP-NCAR reanalysis from 1951 to 2000 (Kalnay et al., 1996); and precipitation from the GPCP, version 2.2, combined precipitation dataset from 1979 to 2012 (Huffman et al., 2009). These data are referred to as "observations" in the following. Interannual anomalies of all variables are defined as the deviations from their corresponding mean seasonal cycle.

3. ENSO variability within the tropical Pacific Ocean
  • In this section, the tropical Pacific variability is described from a 100-year integration of the HCM and compared with available observations. The comparison between the HCM and IOCAS ICM is also documented in this section.

    A time series of monthly SST anomalies averaged in the NiÑo3.4 region (5°S-5°N, 120°-170°W) from the HCM, IOCAS ICM and observations are shown in Fig. 1. It shows that the HCM is able to realistically simulate the interannual variability of SST anomalies in the tropical Pacific. As in the observations, the simulated time series of NiÑo3.4 in the HCM exhibits pronounced irregularity. For example, there are extreme El NiÑo events that peak in 1983 and 1998 in the observations. Similar warm events appear during the model years of 20-40 and 65-75, with anomalies of up to 3°C. The model also shows a suppressed period with anomalies of less than 1°C for model years 1-20 and 40-50, similar to the observations during the years 1975-80. However, ENSO events in the ICM are quite regular (Fig. 1c), probably due to the absence of modulation from the stochastic forcing in the atmosphere (Zhang et al., 2008). In addition, as shown in Fig. 1, the simulated ENSO events in the ICM are significantly stronger than those in the HCM and observations. The standard deviation of NiÑo3.4 index is often used to measure the amplitude of ENSO, and the standard deviation derived from the observations, ICM and HCM are 0.85°C, 1.33°C and 0.95°C, respectively. Compared with observations, the amplitude of ENSO is substantially overestimated in the ICM, whereas the ENSO amplitude in the HCM is comparable to observation.

    Figure 1.  Time series of SST anomalies (units: °C) averaged over the NiÑo3.4 region from (a) observation, (b, d) the HCM and (c) IOCAS ICM. The results of the HCM and ICM are based on the model periods of 1-100 and 1-50 respectively, while the observations are during 1951-2000.

    It is well known that ENSO events are inclined to lock the peak times to boreal winter (Tziperman et al., 1998; An and Wang, 2001). In the ICM and HCM, the standard deviation of SST anomalies in the NiÑo3.4 region is minimum in late spring and maximum in winter (Fig. 2a). That is, the feature of phase-locking is effectively captured by both the ICM and HCM. However, the magnitude of the standard deviation is excessively large in the ICM, while it is more reasonable in the HCM.

    Figure 2.  (a) Standard deviation of SST anomalies (units: °C) over the NiÑo3.4 region as a function of calendar month. (b) Power spectra of NiÑo3.4 index in the models and observation.

    The power spectra of the NiÑo3.4 SST anomaly time series are illustrated in Fig. 2b. The dominant observed period of the NiÑo3.4 index is four years, with a broad spectrum between two and five years. In the ICM, the SST variation is characterized by a sharp peak at a period of three years, but the dominant frequency presents a much too strong power and a too narrow band width compared to observations. The ICM run generates another peak at a period of 18 months, while it is non-significant in observations. In contrast, the HCM shows a broader peak between two and five years, although the largest spectral peak in the simulation is at three years. In addition, the power of the dominant frequency in the HCM is also closer to observation. In other words, the HCM simulation compares better with the observed spectrum than the ICM in terms of dominant period and spectral range. The significant difference in spectra between the HCM and ICM is mainly due to the modulation of atmospheric stochastic forcing, which is implicitly present in the atmospheric component of HCM but missing in the ICM.

    In Fig. 3, the standard deviation of SST in the tropical Pacific is shown to further quantify the interannual SST variability. In the observations, the strongest interannual SST variability is located in the eastern equatorial Pacific and near the Peruvian coast. The HCM reproduces the observed maximum variance center near the Peruvian coast but overestimates the interannual variability of SST in the central Pacific. In contrast, the significant SST variability in the ICM tends to occur in the central equatorial Pacific, extending a little too far west than observed. Also, the interannual variance of SST in the ICM is much stronger than observed. Therefore, compared to the ICM, the HCM can capture the amplitude and structure of interannual variability in the tropical Pacific more effectively. By the way, both the HCM and ICM simulations cannot capture the so-called EP and CP El NiÑo events (Ashok et al., 2007; Kao and Yu, 2009; Kug et al., 2009; Timmermann et al., 2018).

    Figure 3.  Standard deviation of interannual variability of tropical Pacific SST (units: °C) in (a) observation, (b) the HCM, and (c) IOCAS ICM.

    The evolution of SST, zonal wind stress and sea level anomalies on the equator from the models and observations is shown in Fig. 4. As in observations, both the ICM and HCM exhibit a coherent interannual oscillation between corresponding atmospheric and oceanic anomaly fields. The SSTAs are characterized by a standing pattern with the largest anomalies centered in the central-eastern tropical Pacific (Figs. 4a-c), while the maximum wind stress anomalies are located near the date line (Figs. 4d-f). During the development of ENSO events, the zonal wind stress anomalies feature an eastward propagation from the western equatorial Pacific into the central Pacific, consistent with that found in nature. The sea level anomalies, closely related to the variation of the thermocline, show an obvious eastward phase propagation along the equator, initially emerging in the western equatorial Pacific and then propagating eastward as they amplify (Figs. 4g-i). The simulated anomaly fields in the ICM are much smoother than those in the HCM, while the HCM variability is quite irregular, as illustrated in the NiÑo3.4 time series (Fig. 1), and thus reproduces a more realistic interannual variability as in observations.

    Figure 4.  Evolution of SST (top panels; units: °C), zonal wind stress (middle panels; units: dyn cm-2), and sea level (bottom panels; units: cm) anomalies averaged between 2°S and 2°N in the (a, d, g) observation, (b, e, h) HCM and (c, f, i) IOCAS ICM. The results of the models are based on the model-year period of 11-25 and the observations during 1986-2000, respectively.

4. Atmospheric variability in response to ENSO-related SST forcing
  • It is well known that ENSO not only affects regions neighboring the tropical Pacific, but also induces changes in the atmospheric circulation far away from the forcing region and excites global-scale impacts (Ropelewski and Halpert, 1987; Trenberth et al., 1998; Diaz et al., 2001). In fact, ENSO-related short-term climate prediction can be considered to have two parts: one for the SST anomalies in the tropical Pacific and the other for global climate variability induced by these SST anomalies. Although ENSO is now predictable several months ahead (Jin et al., 2008; Barnston et al., 2012, Barnston et al., 2017), it has become even more challenging to understand the processes and mechanisms by which tropical SST-induced signals affect conditions in other remote regions and to further predict their impacts effectively. By coupling the tropical Pacific Ocean to the AGCM (i.e., ECHAM5), the HCM constructed in this study can directly simulate and predict the ENSO-related global atmospheric interannual variability. In this section, the simulation of the global atmospheric response to ENSO in the HCM is examined.

    To evaluate the model's performance, the correlations are examined between the SST anomalies averaged over the NiÑo3.4 region and anomalies of SLP, 200-hPa zonal wind, 2-m temperature, precipitation, and 500-hPa geopotential height, for boreal winter (December-January, February; DJF) and summer (June-July-August, JJA), when the global atmosphere exhibits the clearest response to the ENSO-related SST forcing.

    Results for the SLP and 200-hPa zonal wind are shown in Figs. 5 and 6 for DJF and JJA, respectively, for the model and observations. For DJF, the overall correlation patterns from the model and observations show a large degree of similarity, although the HCM yields somewhat weaker correlations outside the tropical Pacific than those found in the observations. For SLP, the Southern Oscillation signature dominates, with anomalies of opposite sign over the western and eastern tropical Pacific (Figs. 5a and c). Consistent with the SLP, a weakening of the upper-tropospheric westerlies in the eastern tropical Pacific is simulated, with a meridional wave-like structure straddling the equator (Figs. 5b and d). The model correlations for JJA also represent good approximations to those observed in the Pacific Ocean. However, there are obvious differences outside the Pacific. For example, the HCM is unable to reproduce the positive correlations of SLP in the tropical Atlantic (Fig. 6c). The correlations of 200-hPa zonal wind anomalies in the tropical Atlantic and Indian Ocean also seem to be much weaker than those inferred from observations (Fig. 6d). This can be ascribed to the fact that SST fields in these ocean basins are prescribed to the climatology without interannual variations.

    Figure 5.  Correlations of (a, c, e) SLP (units: hPa) and (b, d, f) 200-hPa zonal wind (U200; untis: m s-1) with NiÑo3.4 index for DJF from observation (top panels), the HCM (middle panels) and ECHAM5 forced by observed SST (bottom panels). The black dotted areas indicate that the correlations are statistically significant at the 5% level, as determined by a t-test.

    Figure 6.  As in Fig. 5 but for JJA.

    Figures 7a and c compare the observed and modeled correlations between the 2-m temperature and NiÑo3.4 index for DJF. In observations, significant positive correlations are in the central and eastern tropical Pacific. One can also find similar positive correlations in the tropical Indian Ocean and along the Pacific coast in North and South America. In contrast, significant negative correlations are seen over the western parts of the Pacific spreading eastward from the western equatorial Pacific in both hemispheres and showing a characteristic horseshoe pattern. The model correlations represent good approximations to those observed in the Pacific Ocean basin but the meridional band width of high correlations is smaller than in nature. Outside the Pacific, correlations are too weak compared to those observed over certain regions, such as the Indian Ocean and South Pacific. Figure 7 also shows the correlation patterns of precipitation. Strong positive correlations are observed in the central equatorial Pacific, with a narrower meridional width than those of the surface temperature (Figs. 7b and d). Significant negative correlations are generally found in the subtropical Pacific, Maritime Continent region, and northern Brazil. The model's responses closely match those observed, especially in the central equatorial Pacific and subtropical Pacific, where the correlations between the precipitation and ENSO are positive and negative, respectively. The response over the Maritime Continent is somewhat underestimated by the model, while stronger negative correlations are found in the north band of the subtropical Pacific. Positive precipitation correlation across the southern United States and negative correlation over Brazil are reproduced well by the model, whereas the response over Africa and the Indian Ocean is not well matched to what is observed.

    Figure 7.  Correlations of (a, c, e) 2-m temperature (TMP; units: °C) and (b, d, f) precipitation (PRECIP; units: mm d-1) with NiÑo3.4 index for DJF from observation (top panels), the HCM (middle panels) and ECHAM5 forced by observed SST (bottom panels). The black dotted areas indicate that the correlations are statistically significant at the 5% level, as determined by a t-test.

    Comparing the correlation patterns between JJA and DJF (Figs. 7 and 8), it clearly shows that, for the 2-m temperature, there is a significant warming signal throughout the year in the central-eastern tropical Pacific. Different from the temperature, the integrated global responses of precipitation in JJA are much weaker than those in DJF of the HCM, as well as in the observations. In JJA, the model simulations do well in reproducing the positive correlations of precipitation in the central tropical Pacific, but the response over the Maritime Continent is underestimated and overly strong negative correlations compared to those observed are found in the south and north bands of the subtropical Pacific (Fig. 8d).

    Figure 8.  As in Fig. 7 but for JJA.

    To evaluate the performance of the HCM in reproducing the observed teleconnections between ENSO and mid-tropospheric atmospheric circulation, the correlation fields between NiÑo3.4 and 500-hPa geopotential height are presented in Fig. 9 for DJF and JJA. Positive (negative) correlations imply high (low) heights corresponding to El NiÑo (La NiÑa) events. In the tropical and subtropical latitudes, positive correlations are high at all longitudes, indicating a tropospheric warming during El NiÑo events. A wave train pattern spans the Pacific-American region, which initially emerges in the western equatorial Pacific and then propagates into the extratropics. It consists of three centers of pressure anomalies, located over the northeastern Pacific, Canada, and southeastern USA, respectively. The strongest circulation impact is felt in DJF, when the correlation amplitude is nearly double compared to that in JJA. The model captures the high positive correlations throughout the tropics but with a smaller meridional extent than in nature.

    As mentioned in section 2.2, we conducted a stand-alone AGCM experiment, in which ECHAM5 was forced with observed interannual SST in the tropical Pacific Ocean. The results are shown in Figs. 5-9 in order to compare with the HCM properties. They show that ENSO-related global atmospheric anomalies in the surface climate and circulation simulated by the HCM are in a better agreement with those in ECHAM5 than observations, both with respect to the amplitude and spatial structure of correlations. For example, during JJA in Fig. 6c, the observed positive correlations of SLP in the tropical Atlantic are absent in the HCM, which is also the case for ECHAM5 (Fig. 6e). In the Indian and Atlantic oceans, both models show weaker correlation signals of 200-hPa zonal wind anomalies than observed. For the temperature and precipitation, the weak responses outside the tropical Pacific in the HCM are consistent with those in ECHAM5, while this is not the case for the observation (Figs. 7 and 8). The lack of agreement outside of the tropical Pacific between the HCM and observations can be due to the prescribed SST climatology without interannual variations beyond the tropical Pacific basin. This is evidenced by the good agreement between the HCM and ECHAM5, which is forced with observed interannual SST in the tropical Pacific.

    Figure 9.  Correlations between 500-hPa geopotential height (H500) and NiÑo3.4 index for (a, c, e) DJF and (b, d, f) JJA from observation (top panels), the HCM (middle panels) and ECHAM5 forced by observed SST (bottom panels). The black dotted areas indicate that the correlations are statistically significant at the 5% level, as determined by a t-test.

5. Summary and discussion
  • In this study, an HCM is presented that couples a tropical Pacific Ocean model (i.e., the ocean component of IOCAS ICM) to a global spectral AGCM. A 100-year HCM-based simulation is examined for ENSO and the related global atmospheric response, the results of which are compared to available observations and to a simulation using IOCAS ICM.

    The HCM is able to reproduce realistic and irregular occurrence of ENSO events with substantial active and inactive periods. Both the observed and simulated NiÑo3.4 SST anomaly shows a broad spectrum of oscillation periods between two and five years, although the dominant period is four years in the observations and three years in the model. The HCM successfully locates the maximum center of ENSO variance near the Peruvian coast, but somewhat overestimates the SST variability in the central tropical Pacific. The model ENSO indicates reasonable annual phase-locking with minimum variance in late spring and maximum variance in boreal winter. Time-longitude distributions of observed and the HCM-simulated SST, zonal wind stress, and sea level anomalies illustrate that these fields exhibit a coherent interannual variation. In the simulation, a standing pattern of SST anomalies exists, with the largest anomalies in the central-eastern tropical Pacific, while the largest variation of zonal wind stress anomalies is near the date line. Eastward propagation is apparent for zonal wind stress and sea level anomalies, similar to that in nature. In contrast, the ENSO cycle in IOCAS ICM is far too regular, dominated by a three-year oscillation; and the maximum interannual variance of SST is too far west and with a much stronger amplitude.

    It is also found that the HCM is capable of characterizing realistic atmospheric variations in response to ENSO-related SSTA forcing for DJF and JJA, respectively. Particularly, similar to the observations, a high correlation is seen between the tropical Pacific SST anomalies and atmospheric response in the winter. This is evidenced by examining the correlation fields between NiÑo3.4 index and anomalies of SLP, 200-hPa zonal wind, 2-m temperature, precipitation, and 500-hPa geopotential height. Within the tropics, the signature of Southern Oscillation is well reproduced, characterized by opposite anomalies of SLP over the western and eastern tropical Pacific, with a weakening of the upper-tropospheric westerlies. The model also reasonably captures the response of 2-m temperature and atmospheric circulation in the middle troposphere, including the tropospheric warming throughout the tropics and the remote impacts in the Pacific-American sector. The model-simulated correlation patterns of precipitation in DJF closely match the observations; however, the patterns in JJA deviate from the observations over some regions, such as the weaker signals in the Maritime Continent and stronger signals in the north and south bands of the subtropical Pacific. Compared with observations, the HCM is in better agreement with the stand-alone ECHAM5 forced with observed interannual SST in the tropical Pacific, both with respect to the amplitude and spatial structure of the atmospheric response. Therefore, the lack of agreement outside the tropical Pacific Ocean between the HCM and observations is mainly due to the fact that SST fields in these regions are prescribed to the climatology without interannual variations.

    Compared to IOCAS ICM, the HCM shows clear advantages for ENSO-related studies. For example, the simulated ENSO characteristics show substantial improvement in the HCM compared to those in IOCAS ICM, including a more reasonable irregularity, amplitude, periodicity and life cycle of ENSO. By adopting an AGCM, the HCM can realistically describe the atmospheric variability, and thus can represent the nonlinear processes of the coupled system better than the ICM. Also, the HCM can characterize the global atmospheric responses to ENSO, which are not possible in the tropics-only ICM. In addition, by employing an anomaly coupling strategy, this HCM not only is free from the problem of climate drift, but also can effectively reduce the computational cost compared to fully coupled general circulation models. Nevertheless, like other HCMs, this HCM also has its disadvantages. For example, the atmosphere and ocean are only actively coupled in the tropical Pacific Ocean in the HCM, which may significantly underestimate the atmospheric response outside the tropical Pacific, particularly at regional or local scales. The tropical oceans have long been recognized as a key player in global climate, in which large-scale convection over the warm water provides an important portion of the driving energy for the general circulation of the atmosphere. Therefore, as a step forward, the HCM constructed in this work can be further expanded to include the tropical Indian and Atlantic oceans, with the expectation of improving the simulation of global climate. In fact, the relevant work is currently underway and will hopefully be reported in the future.

    This paper presents our efforts in developing a new HCM and a preliminary analysis of the HCM-based simulation for ENSO variability and the associated global atmospheric response. There are still questions that remain to be answered, and further detailed analyses are clearly needed. For example, compared to IOCAS ICM, the HCM produces a more realistic irregularity of the ENSO cycle. But what model physics are responsible for the ENSO regularity? As mentioned in the introduction, an increasing number of studies are highlighting the potential influence of stochastic atmospheric noise (e.g., WWBs) on the complexity and predictability of ENSO events (Zhang et al., 2008; Lopez and Kirtman, 2013, Lopez and Kirtman, 2014; Hu et al., 2014; Lian et al., 2014; Chen et al., 2015). With the HCM constructed in this study, the impact of stochastic atmospheric forcing on ENSO needs to be further discussed. In addition, the diagnoses of the simulation suggest that the HCM may be suitable for the prediction of interannual climate anomalies associated with ENSO. But is the HCM more skillful in forecasting ENSO and the related climate variability compared with those in IOCAS ICM? These questions should be addressed by conducting hindcast experiments and more detailed diagnoses in the future.

Reference

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return