Progress in the Development and Application of Climate Ocean Models and Ocean-Atmosphere Coupled Models in China

A review is presented about the development and application of climate ocean models and oceanatmosphere coupled models developed in China as well as a review of climate variability and climate change studies performed with these models. While the history of model development is briefly reviewed, emphasis has been put on the achievements made in the last five years. Advances in model development are described along with a summary on scientific issues addressed by using these models. The focus of the review is the climate ocean models and the associated coupled models, including both global and regional models, developed at the Institute of Atmospheric Physics, Chinese Academy of Sciences. The progress of either coupled model development made by other institutions or climate modeling using internationally developed models also is reviewed.     

Comparison Between GAMIL, and CAM2 on Interannual Variability Simulation

Recently, a new atmospheric general circulation model （GAMIL： Grid-point Atmospheric Model of IAP LASG） has been developed at the Institute of Atmospheric Physics （IAP）, Chinese Academy of Sciences （CAS）, which is based on the Community Atmospheric Model Version 2 （CAM2） of the National Center for Atmospheric Research （NCAR）. Since the two models have the same physical processes but different dynamical cores, the interannual variability simulation performances of the two models are compared. The ensemble approach is used to reduce model internal variability. In general, the simulation performances of the two models are similar. Both models have good per- formance in simulating total space-time variability and the Southern Oscillation Index. GAMIL performs better in the Eastern Asian winter circulation simulation than CAM2, and the model internal variability of GAMIL has a better response to external forcing than that of CAM2. These indicate that the improvement of the dynamic core is very important. It is also verified that there is less predictability in the middle and high latitudes than in the low latitudes.     

The Performance of Atmospheric Component Model R42L9 of GOALS／LASG

This paper examines the performance of an atmospheric general circulation model (AGCM) developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics (LASG/IAP). It is a spectral model truncated at R42(2.8125°long×1.66°lat) resolution and with nine vertical levels, and referred to as R42L9/LASG hereafter. It is also the new version of atmospheric component model R15L9 of the global ocean-atmosphere-land system (GOALS/LASG). A 40-year simulation in which the model is forced with the climatological monthly mean sea surface temperature is compared with the 40-year (1958-97) U.S. National Center for Environmental Prediction (NGEP) global reanalysis and the 22-year (1979-2000) Xie-Arkin monthly precipitation climatology. The mean DJF and JJA geographical distributions of precipitation, sea level pressure, 500-hPa geopotential height, 850-hPa and 200-hPa zonal wind, and other fields averaged for the last 30-year integrat     

The Vertical Structures of Atmospheric Temperature Anomalies Associated with El Nio Simulated by the LASG/IAP AGCM： Sensitivity to Convection Schemes

The vertical structures of atmospheric temperature anomalies associated with El Nio are simulated with a spectrum atmospheric general circulation model developed by LASG/IAP (SAMIL). Sensitivity of the model’s response to convection scheme is discussed. Two convection schemes, i.e., the revised Zhang and Macfarlane (RZM) and Tiedtke (TDK) convection schemes, are employed in two sets of AMIP-type (Atmospheric Model Intercomparison Project) SAMIL simulations, respectively. Despite some deficiencies in the up...     

一个灵活的海洋——大气耦合环流模式

Based on the National Center for Atmospheric Research (NCAR) Climate System Model version 1(CSM-1), a Flexible coupled General Circulation Model version 0 (FGCM-0) is developed in this study through replacing CSM-1＇s oceanic component model with IAP L30T63 global oceanic general circulation model and some necessary modifications of the other component models. After the coupled model FGCM--0 is spun up for dozens of years, it has been run for 60 years without flux correction. The model does not only show the reasonable long-term mean climatology, but also reproduce a lot of features of the interannual variability of climate, e.g. the ENSO-like events in the tropical Pacific Ocean and the dipole mode pattern in the tropical Indian Ocean. Comparing FGCM-0 with the NCAR CSM-1, some common features are found, e.g. the overestimation of sea ice in the North Pacific and the simulated double ITCZ etc.The further analyses suggest that they may be attributed to errors in the atmospheric model.     

A fast version of LASG/IAP climate system model and its 1000-year control integration

A fast version of the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geo- physical Fluid Dynamics(LASG)/Institute of Atmospheric Physics(IAP)climate system model is briefly documented.The fast coupled model employs a low resolution version of the atmospheric component Grid Atmospheric Model of IAP/LASG(GAMIL),with the other parts of the model,namely an oceanic component LASG/IAP Climate Ocean Model(LICOM),land component Common Land Model(CLM),and sea ice component from National Center for Atmospheric Research Community Climate System Model (NCAR CCSM2),as the same as in the standard version of LASG/IAP Flexible Global Ocean Atmosphere Land System model(FGOALS g).The parameterizations of physical and dynamical processes of the at- mospheric component in the fast version are identical to the standard version,although some parameter values are different.However,by virtue of reduced horizontal resolution and increased time-step of the most time-consuming atmospheric component,it runs faster by a factor of 3 and can serve as a useful tool for long- term and large-ensemble integrations.A 1000-year control simulation of the present-day climate has been completed without flux adjustments.The final 600 years of this simulation has virtually no trends in global mean sea surface temperatures and is recommended for internal variability studies.Several aspects of the control simulation’s mean climate and variability are evaluated against the observational or reanalysis data. The strengths and weaknesses of the control simulation are evaluated.The mean atmospheric circulation is well simulated,except in high latitudes.The Asian-Australian monsoonal meridional cell shows realistic features,however,an artificial rainfall center is located to the eastern periphery of the Tibetan Plateau persists throughout the year.The mean bias of SST resembles that of the standard version,appearing as a"double ITCZ"(Inter-Tropical Convergence Zone)associated with a westward extension of the equatorial eastern Pacific cold tongue.The sea ice extent is acceptable but has a higher concentration.The strength of Atlantic meridional overturning is 27.5 Sv.Evidence from the 600-year simulation suggests a modulation of internal variability on ENSO frequency,since both regular and irregular oscillations of ENSO are found during the different time periods of the long-term simulation.     

Modeling the tropical Pacific Ocean using a regional coupled climate model

A high-resolution tropical Pacific general circulation model (GCM) coupled to a global atmospheric GCM is described in this paper. The atmosphere component is the 5°×4°global general circulation model of the Institute of Atmospheric Physics (IAP) with 9 levels in the vertical direction. The ocean component with a horizontal resolution of 0.5°, is based on a low-resolution model (2°×1°in longitude-latitude).Simulations of the ocean component are first compared with its previous version. Results show that the enhanced ocean horizontal resolution allows an improved ocean state to be simulated; this involves (1) an apparent decrease in errors in the tropical Pacific cold tongue region, which exists in many ocean models,(2) more realistic large-scale flows, and (3) an improved ability to simulate the interannual variability and a reduced root mean square error (RMSE) in a long time integration. In coupling these component models, a monthly "linear-regression" method is employed to correct the model's exchanged flux between the sea and the atmosphere. A 100-year integration conducted with the coupled GCM (CGCM) shows the effectiveness of such a method in reducing climate drift. Results from years 70 to 100 are described.The model produces a reasonably realistic annual cycle of equatorial SST. The large SSTA is confined to the eastern equatorial Pacific with little propagation. Irregular warm and cold events alternate with a broad spectrum of periods between 24 and 50 months, which is very realistic. But the simulated variability is weaker than the observed and is also asymmetric in the sense of the amplitude of the warm and cold events.     

Performance of a Reconfigured Atmospheric General Circulation Model at Low Resolution

Paleoclimate simulations usually require model runs over a very long time.The fast integration version of a state-of-the-art general circulation model (GCM),which shares the same physical and dynamical processes but with reduced horizontal resolution and increased time step,is usually developed.In this study,we configure a fast version of an atmospheric GCM (AGCM),the Grid Atmospheric Model of IAP/LASG (Institute of Atmospheric Physics/State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics),at low resolution (GAMIL-L,hereafter),and compare the simulation results with the NCEP/NCAR reanalysis and other data to examine its performance.GAMIL-L,which is derived from the original GAMIL,is a finite difference AGCM with 72×40 grids in longitude and latitude and 26 vertical levels.To validate the simulated climatology and variability,two runs were achieved.One was a 60-year control run with fixed climatological monthly sea surface temperature (SST) forcing,and the other was a 50-yr (1950-2000) integration with observational time-varying monthly SST forcing.Comparisons between these two cases and the reanalysis,including intra-seasonal and inter-annual variability are also presented.In addition,the differences between GAMIL-L and the original version of GAMIL are also investigated. The results show that GAMIL-L can capture most of the large-scale dynamical features of the atmosphere, especially in the tropics and mid latitudes,although a few deficiencies exist,such as the underestimated Hadley cell and thereby the weak strength of the Asia summer monsoon.However,the simulated mean states over high latitudes,especially over the polar regions,are not acceptable.Apart from dynamics,the thermodynamic features mainly depend upon the physical parameterization schemes.Since the physical package of GAMIL-L is exactly the same as the original high-resolution version of GAMIL,in which the NCAR Community Atmosphere Model (CAM2) physical package was used,there are only small differences between them in the precipitation and temperature fields.Because our goal is to develop a fast-running AGCM and employ it in the coupled climate system model of IAP/LASG for paleoclimate studies such as ENSO and Australia-Asia monsoon,particular attention has been paid to the model performances in the tropics.More model validations,such as those ran for the Southern Oscillation and South Asia monsoon, indicate that GAMIL-L is reasonably competent and valuable in this regard.     

The tropical intraseasonal oscillation in SAMIL coupled and uncoupled general circulation models

Simulations of tropical intraseasonal oscillation(TISO) in SAMIL,the Spectral Atmospheric Model from the Institute of Atmospheric Physics(IAP) State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics(LASG) coupled and uncoupled general circulation models were comprehensively evaluated in this study.Compared to the uncoupled model,the atmosphere-ocean coupled model improved the TISO simulation in the following aspects:(1) the spectral intensity for the 30-80-day peak eastward periods was more realistic;(2) the eastward propagation signals over western Pacific were stronger;and(3) the variance distribution and stronger signals of Kelvin waves and mixed Rossby gravity waves were more realistic.Better performance in the coupled run was assumed to be associated with a better mean state and a more realistic relationship between precipitation and SST.In both the coupled and uncoupled runs,the unrealistic simulation of the eastward propagation over the equatorial Indian Ocean might have been associated with the biases of the precipitation mean state over the Indian Ocean,and the unrealistic split of maximum TISO precipitation variance over the Pacific might have corresponded to the exaggeration of the double Intertropical Convergence Zone(ITCZ) structure in precipitation mean state.However,whether a better mean state leads to better TISO activity remains questionable.Notably,the northward propagation over the Indian Ocean during summer was not improved in the mean lead-lag correlation analysis,but case studies have shown some strong cases to yield remarkably realistic northward propagation in coupled runs.     

Impacts of a GCM＇s Resolution on MJO Simulation

Long-term integrations are conducted using the Spectral Atmospheric Model （referred to as SAMIL）, which was developed in the Laboratory for Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics （LASG） in the Institute of Atmospheric Physics （IAP）, with different resolutions to inves-tigate sensitivity of the Madden-Julian Oscillation （MJO） simulations to the model＇s resolution （horizontal and vertical）. Three resolutions of the model, R15L9, R42L9 and R42L26, with identical physical processes, all produced the basic observed features of the MJO, including the spatiotemporal space-time spectra and eastward propagation. No fundamental differences among these simulations were found. This indicates that the model resolution is not a determining factor for simulating the MJO. Detailed differences among these modeling results suggest, however, that model resolution can substantially affect the simulated MJO in certain aspects. For instance, at a lower horizontal resolution, high frequency disturbances were weaker and the structures of the simulated MJO were better defined to a certain extent. A higher vertical resolution led to a more realistic spatiotemporal spectrum and spatial distribution of MJO precipitation. Meanwhile, increasing the model＇s resolution improved simulation of the climatology. However, increasing the resolution should be based on improving the cumulus parameterization scheme.     

SAMIL模式中Tiedtke 积云对流方案对热带降水模拟的影响

目前, 大多数全球耦合模式及大气环流模式在降水模拟中普遍存在不同程度的“热带偏差”问题, 中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室所发展的全球大气环流谱模式SAMIL-R42L26也存在这一现象, 主要表现在SPCZ (南太平洋辐合区) 降水过强且过分东伸、 赤道附近降水偏少等方面。本文通过修改SAMIL中的积云对流方案有效地削弱了这一偏差, 并进一步探讨其原因, 发现对流方案修改后, 改变了对流层低层至地面的温度分布状态, 进而影响了风速及散度场的模拟, 最终通过垂直速度的调整反作用于对流过程。比较修改前后对流过程云底质量通量, 发现修改后的方案主要通过削弱浅对流来提高热带降水的模拟性能。     

区域海气耦合模式FROALS的发展及其应用

区域海气耦合模式是研究局地海气相互作用过程影响气候变率的重要平台,也是对全球气候模式进行"动力降尺度"的重要工具.本文介绍了LASG(State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics)/IAP(Institute of Atmospheric Physics)发展的区域海气耦合模式FROALS(Flexible Regional Ocean-Atmosphere-Land System model),并总结了过去五年围绕该区域海气耦合模式开展的研究工 作.FROALS的特点之一是有两个完全不同的大气模式分量和海洋模式分量选项,可以适应不同的模拟研究需 求.针对区域海气耦合模式在西北太平洋地区的模拟偏差,通过分步骤考察不同大气模式分量和不同海洋模式分量对模式模拟性能的影响,指出大气模式是导致区域海气耦合偏差的主要分量.通过改进对流触发的相对湿度阈值标准,有效地改善了此前区域海气耦合模式在亚洲季风区普遍出现的"模拟海温冷偏差".改进的FROALS对西北太平洋地区的大气和海洋环境有较好的模拟能力,合理地再现了西北太平洋地区表层洋流气候态和年际变率.较之非耦合模式,考虑区域海气耦合过程后,改进了东亚和南亚地区的降水和热带气旋潜势年际变率的模拟.最后,针对东亚—西北太平洋地区,利用FROALS对IAP/LASG全球气候模式模拟和预估的结果进行了动力降尺 度,得到了东亚区域50 km高分辨率区域气候变化信息.分析显示,FROALS模拟得到的东亚区域气候较之全球气候模式和非耦合区域气候模式结果具有明显的"增值",显示出区域海气耦合模式在该区域良好的应用前景.     

SAMIL大气环流模式海面湍流通量参数化方案研究

将中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室发展的大气环流谱模式SAMIL_R42L26 2.08中的海气通量参数化方案——Louis方案, 与新发展的一种新型近海层湍流通量参数化方案——LGLC方案进行比较和分析。离线测试结果表明, LGLC方案计算的通量结果与观测数据具有更好的一致性, 且由于其区分热力粗糙度和动力粗糙度, 使得对热量通量的计算更加准确。在线测试则证明, 引入LGLC方案的SAMIL模式对洋面风应力、 感热通量、 潜热通量和降水率的模拟能力有了进一步的提高, 尤其对北半球夏季印度季风和南海季风区的降水改善明显。     

Impacts of cumulus momentum transport on MJO simulation

Vertical cumulus momentum transport is an important physical process in the tropical atmosphere and plays a key role in the evolution of the tropical atmospheric system. This paper focuses on the impact of the vertical cumulus momentum transport on Madden-Julian Oscillation (MJO) simulation in two global climate models (GCMs). The Tiedtke cumulus parameterization scheme is applied to both GCMs [CAM2 and Spectral Atmospheric general circulation Model of LASG/IAP (SAMIL)]. It is found that the MJO simulation ability might be influenced by the vertical cumulus momentum transport through the cumulus parameterization scheme. However, the use of vertical momentum transport in different models provides different results. In order to improve model's MJO simulation ability, we must introduce vertical cumulus momentum transport in a more reasonable way into models. Furthermore, the coherence of the parameterization and the underlying model also need to be considered.     

Development of global coupled ocean-atmosphere general circulation models

It has long been believed that a climate model capable of realistically simulating many features of global climate, variability, and climate change must interactively represent the major components of the dynamically coupled climate system, particularly the atmosphere, ocean, and cryosphere. This effort traditionally has been constrained by computing power, our understanding of the observed system, and climate modeling capability. With the advent of supercomputers, improved understanding of global climate processes, and computationally efficient general circulation climate models, we have witnessed a rapid increase in the simulation of global climate by coupling together various representations of atmosphere, ocean, and sea ice. Beginning in the late 1960s and continuing through the early 1980s, general circulation models (GCMs) of the atmosphere, ocean, and sea ice were coupled and run asynchronously to produce credible simulations of the global climate. Systematic errors in these component models later led some modeling groups to use flux correction or flux adjustment, whereby either one or several of the variables at the air-sea interface are adjusted to bring the simulations in closer agreement with observations. Further advances in computing power and climate modeling techniques in the past few years have allowed global coupled ocean-atmosphere GCMs to be run synchronously (i.e., atmosphere and ocean communicate at least once each model day). Computing constraints, combined with the need for multidecadal climate integrations, still only allow relatively coarse-grid ocean GCMs to be coupled to correspondingly coarse-grid atmospheric models (on the order of 500 km × 500 km). However, results from this current generation of global, coupled GCMs have revealed interesting characteristics associated with ocean dynamics and global climate in experiments with gradual increases of carbon dioxide. Another somewhat surprising aspect of the global-coupled GCM simulations is the appearance of some features associated with the El Niño-Southern Oscillation. Along with concurrent efforts with other types of limited-domain, dynamical coupled models, this has led to the realization that inherent unstable coupled modes exist in the climate system that are the unique product of the interactive coupling of the atmosphere and the ocean. All of these efforts are leading to the next generation of coupled ocean-atmosphere GCMs. These models will run on even faster and larger-memory computers and will have higher-resolution atmosphere and ocean components, more accurate sea-ice formulations, improved cloud-radiation schemes, and increasingly realistic land-surface processes.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. DümenilThe National Center for Atmospheric Research is sponsored by the National Science Foundation     

海气耦合模式FGOALS_gl模拟的水汽和云辐射反馈过程

本文分析了中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室(LASG)发展的快速耦合模式FGOALS_gl对低纬太平洋区域水汽温室效应(Ga)和云辐射强迫特征的模拟能力,讨论了模拟偏差的成因.结果表明,FGOALS_gl能合理再现Ga、云辐射强迫的气候态空间分布特征,但也存在明显的偏差.模式低估了冷...     

A regional response in mean westerly circulation and rainfall to projected climate warming over Tasmania, Australia

Coupled ocean–atmosphere general circulation models (GCMs) lack sufficient resolution to model the regional detail of changes to mean circulation and rainfall with projected climate warming. In this paper, changes in mean circulation and rainfall in GCMs are compared to those in a variable resolution regional climate model, the Conformal Cubic Atmospheric Model (CCAM), under a high greenhouse gas emissions scenario. The study site is Tasmania, Australia, which is positioned within the mid-latitude westerlies of the southern hemisphere. CCAM projects a different response in mean sea level pressure and mid-latitude westerly circulation to climate warming to the GCMs used as input, and shows greater regional detail of the boundaries between regions of increasing and decreasing rainfall. Changes in mean circulation dominate the mean rainfall response in western Tasmania, whereas changes to rainfall in the East Coast are less related to mean circulation changes. CCAM projects an amplification of the dominant westerly circulation over Tasmania and this amplifies the seasonal cycle of wet winters and dry summers in the west. There is a larger change in the strength than in the incidence of westerly circulation and rainfall events. We propose the regional climate model displays a more sensitive atmospheric response to the different rates of warming of land and sea than the GCMs as input. The regional variation in these results highlight the need for dynamical downscaling of coupled general circulation models to finely resolve the influence of mean circulation and boundaries between regions of projected increases and decreases in rainfall.     

耦合模式FGOALS_s模拟的东亚夏季风

本文评估了中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室的海气耦合模式FGOALS_s对东亚夏季风的模拟能力, 并通过与观测海温强迫下单独大气模式SAMIL试验结果的比较, 分析了海气耦合过程对模式性能的影响。结果表明, FGOALS_s基本能够模拟出东亚夏季风系统的气候态分布及其演变过程, 但也存在明显偏差, 主要表现为模拟的温度场在对流层中上层一致性偏冷, 导致模式中环流系统强度偏弱; 而温度经向梯度模拟的不足, 直接影响到东亚副热带西风急流的模拟。通过与观测海温强迫下SAMIL模拟结果的对比发现, SAMIL模拟的温度场、 环流场以及风场较之耦合模式结果更接近观测, 但也存在与FGOALS_s类似的模式偏差。因此, 大气模式固有的偏差对耦合模式的模拟偏差有重要影响。分析发现, 对于西太平洋降水的模拟而言, 耦合模式结果更加合理, 表明海气相互作用过程对模式性能有重要影响。本文的结果表明, 大气模式自身的误差是导致耦合模式误差的主要原因。通过更新云－辐射模块改进大气模式模拟的温度场, 应是FGOALS_s后续发展的首要工作。     

任意正交曲线坐标系下的海洋模式动力框架的发展与评估

本文发展了一个可以适用于任意水平正交曲线坐标系的海洋模式动力框架,并将其应用于中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室发展的气候系统海洋模式LICOM2.0（LASG/IAP Climate system Ocean Model,version2.0）。在经纬网格坐标系下,新的动力框架与LICOM2.0原有的动力框架模拟结果完全一致。基于新的动力框架,海洋模式可采用能够准确描述北冰洋地形的三极网格,克服了LICOM2.0经纬网格版本必须将北极点处理为孤岛的缺陷,从而显著改进了模式对于北冰洋环流和北大西洋经圈翻转流函数（AMOC）的模拟能力。此外,引进三极网格还可以避免模式网格距随纬度增加而急剧减小带来的计算不稳定,在LICOM2.0的三极网格版本中,模式不需要采用任何空间滤波方案仍然能够保证计算的稳定性,从而与LICOM2.0的经纬网格版本相比,极大地提高了模式的并行效率,这一点在当水平分辨率提高到0.1度时表现得尤为明显,海洋模式的并行加速比可以从经纬网格版本的5.8左右提高到三极网格版本的15.0左右。