两个大气环流模式SAMIL和BCC_AGCM对北半球冬季极涡振荡的模拟对比

基于与NCEP再分析资料的比较,本文利用大气环流模式SAMIL和北京气候中心大气环流模式BCC_AGCM的1950～1999年的AMIP试验模拟数据,对北半球冬季环流及平流层极涡振荡的模拟性能进行了评估分析。结果表明两个模式都可以再现北半球环流基本型以及环流振荡的主导模态。对冬季气候平均态的模拟,两个模式模拟的热带—热带外温度梯度均偏大,极夜急流偏强,极涡偏冷偏强;100～20hPa平均位势高度场谐波分析表明两个模式模拟的行星波偏弱;气候平均的10hPa极夜急流均存在1个月的季节漂移,200hPa副热带西风急流较NCEP偏弱。两个模式模拟的环流变化的主导模态均代表极涡振荡,SAMIL极涡振荡的强度大于BCC_AGCM,BCC_AGCM极涡振荡的频率要高于SAMIL。连续功率谱分析表明,NCEP资料中极涡振荡存在4.6个月的显著周期,相应地,SAMIL中为5.5个月的显著周期,BCC_AGCM中为4.8个月。NCEP资料中的极涡振荡主要发生在12～3月,SAMIL模拟的极涡振荡主要发生在2～3月,BCC_AGCM模拟的极涡振荡主要发生在2～4月。

Winter Season Stratospheric Circulation in the SAMIL/LASG General Circulation Model

In this paper, we examine the performance of the 26-level version of the SAMIL/LASG GCM （R42/L26） in simulating the seasonal cycle and perpetual winter mean stratospheric circulation as well as its variability by comparing them with the NCEP/NCAR reanalysis. The results show that the model is capable of repro- ducing many key features of the climatology and seasonal variation of the stratospheric circulation despite that the model＇s mean polar vortex is stronger and more zonally symmetric compared to the observation. ~rther diagnosis of the results from a perpetual-January-run of the SAMIL/LASG GCM indicates that the dominant winter-season oscillation mode in the model＇s stratosphere exhibits a similar inter-seasonal timescale with similar spatial patterns as those inferred from the NCEP/NCAR reanalysis. In particular, the simulated polar vortex oscillation mode exhibits a dominant inter-seasonal timescale of about 120 days, and is accompanied with the simultaneous poleward and downward propagation of temperature anomalies in the stratosphere and the equatorward propagation of temperature anomalies in the troposphere. More encouragingly, the 26-layer version of the SAMIL/LASG GCM is able to produce three strong Stratospheric Sudden Warming events during the 1825 days of perpetual-January integration, with the polar westerly jet completely reversed for a few weeks without imposing any prescribed anomalous forcing at the lower boundary.

潜热通量异常对西北太平洋热带气旋活动影响的机理研究

使用1958—2001年的中国气象局台风年鉴资料和欧洲中心再分析资料,初步揭示了潜热通量异常对西北太平洋(含南海)热带气旋活动的可能影响机理。统计分析发现,西北太平洋(含南海)热带气旋年频数与北太平洋副热带(简称关键区)的潜热通量在过去40余年中均表现为显著的减弱趋势。诊断分析表明,关键区的潜热通量通过低层信风向西的水汽输送—整个西北太平洋热带气旋活动区水汽低层辐合上升而凝结—潜热释放这一链条改变西北太平洋(含南海)大气环境场条件(包括中层湿度场、高低层涡度场和高低层散度场),进而调制热带气旋活动。使用SAMIL模式进行关键区内潜热通量加倍和减半的敏感性试验,进一步证实了关键区潜热通量异常对西北太平洋热带气旋活动调制作用的可能机理,而对南海热带气旋活动影响较小。由此可推断,在过去40余年中,北太平洋副热带中部区潜热通量的减弱趋势可能是造成西北太平洋(含南海)热带气旋频数下降的主要原因之一。

积云对流和云物理过程调整对气候模拟的影响

本文利用中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室(LASG/IAP)的大气环流谱模式SAMIL,结合观测与政府间气候变化专门委员会第四次评估报告(IPCCAR4)大气模式集合平均结果,以大气辐射通量为例,诊断分析了物理过程调整前后模式对气候模拟的影响。旧版本SAMIL对大气辐射通量的模拟存在较大偏差,经过大气辐射过程、积云对流和诊断云等物理过程的调整后,新版本SAMIL模拟的全球辐射通量的年平均结果与观测的偏差大幅减小,其中大气顶能量收支的年变化及其平均值与观测更为接近。在积云对流方案调整基础上,通过对诊断云物理方案的进一步调整,新版本SAMIL对云物理量模拟更为合理,在赤道辐合带等区域,在很大程度上克服了单一积云对流物理过程调整引起的云宏观和微观属性不匹配问题,能模拟出夏季气候平均辐射通量的全球分布特征,尤其在东亚区域有较好的模拟能力。研究还表明,在热带和副热带对流活跃区域,当前SAMIL中积云对流过程偏差对辐射通量的模拟偏差有很大影响,而模式中较为简单的诊断云方案也会将云宏观物理量模拟偏差带入云微观量模拟中,也是主要偏差源之一。本文结果表明,要继续提高SAMIL的模拟性能,急需更新云物理参数化方案以改进云辐射过程的模拟,同时也需要有针对性的研究积云对流和云物理过程之间相互作用,并作进一步协同调整。

热带加热异常影响冬季平流层极涡强度的数值模拟

本文利用大气环流模式SAMIL/LASG,通过选择两种对流参数化方案,研究了热带加热异常对热带外平流层模拟的影响。结果表明,因不同对流参数化方案引起的热带对流加热状况的差异,可显著影响模式对北半球冬季平流层极涡强度的模拟偏差。与采用Manabe对流参数化方案相比,采用Tiedtke参数化方案可以显著改善对平流层极涡强度的模拟,使平流层极涡"过强"及极区"过冷"的模拟偏差得到明显改善。研究其中的影响过程发现,由于Manabe方案最大凝结潜热加热高度过低,在对流层中低层;而Tiedtke方案的最大凝结潜热加热位置在对流层中上层,因而Tiedtke(Manabe)方案时热带大气温度在对流层中上层较为偏暖(偏冷),在平流层低层较为偏冷(偏暖)。自秋季开始,与热带对流层高层温度的暖偏差相联系,热带外对流层高层以及热带平流层低层出现伴随的温度冷偏差;与之对应,平流层中纬度从秋季开始也出现持续的温度暖偏差。另外,随着秋冬季节平流层行星波活动的出现,Tiedtke方案时热带外地区行星波1波的强度也明显强于Manabe方案,使得秋冬季节涡动引起的向极热通量在Tiedtke方案时明显偏强,从而造成了冬季平流层极区温度偏暖、极涡强度偏弱。

基于LASG/IAP大气环流谱模式的气候系统模式

文章扼要介绍了基于LASG/IAP大气环流谱模式(SAMIL)的气候系统模式的新版本FGOALS-s的发展和结构。出于发展一个在东亚季风模拟方面有一定优势的气候系统模式之目的,FGOALS-s的大气模式分量SAMIL采用了较高的水平分辨率R42,这相当于2.8125°(经度)×1.66°(纬度),高于三角截断T42的分辨率。对FGOALS-s在模拟大气、陆面、海洋和海冰的气候平均态,以及主要的年际变率信号方面的能力进行了检验。分析表明,FGOALS-s成功地控制了气候漂移趋势,能够较为真实地模拟大气、海洋和陆面的气候平均态,特别是受益于大气模式的较高分辨率,由中国西南向东北延伸的夏季风雨带的分布,在模式中得到较为真实的再现,表明该模式在东亚夏季风的模拟上具有较强能力。耦合模式能够成功再现El Ni^no事件的非规则周期变化,但是其年际变化的振幅较之观测要弱。赤道中西太平洋年际变率的强度较之赤道中东太平洋要强。在中高纬度,模式模拟的北大西洋涛动模态,在空间分布上与观测接近。FGOALS-s模式存在的主要问题,是模拟的热带海温偏冷、而中纬度海温则偏暖,原因是模式模拟的云量分布存在偏差,它直接影响到海表的净热通量收支。模式模拟的北大西洋高纬度地区的海温明显偏冷,令该地区的年平均海冰分布的范围明显偏大;然而受南极周边海温偏高影响,南极洲周围的海冰范围则偏少。FGOALS-s的未来工作重点,宜放在大气模式的云过程、海洋模式的经向能量输送过程、以及海洋与大气的淡水通量耦合方案的改进方面。

辐射和积云对流过程对大气辐射通量的影响

基于中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室发展的全球大气环流谱模式(SAMIL-R42L26),研究了澳大利亚气象局研究中心(BMRC)新辐射方案和新Zhang-McFarlane积云对流方案对大气辐射通量模拟的影响。新辐射方案相比原辐射方案在辐射计算光谱分辨率、气体吸收和计算效率等方面作了很多改进,其对大气辐射通量的模拟能力相应提高。在晴空条件下,大气顶出射长波、大气吸收短波和地表入射短波等与观测的偏差较原辐射方案明显减小,尤其是在对流活跃区域。在云天条件下大气辐射通量与观测的偏差也较原辐射方案减小,但其偏差依然较大,这与模式中积云对流参数化方案模拟能力不足引起的辐射通量偏差有关。为此,换用了新Zhang-McFarlane积云对流方案,其结果表明,对流活跃区水汽含量显著增加,原对流方案中偏强的"双赤道辐合带"现象明显减弱,赤道辐合带地区的大气辐射通量偏差有明显减小,在海洋地区晴空大气顶出射长波和地表入射短波的量值及空间分布均接近观测结果,同时大气顶全球平均能量收支的年变化和观测结果趋于一致,其中模拟的全球年平均大气顶能量收支和观测的偏差不到0.6W/m2。试验结果同时表明,在未来研究中引入气溶胶分布、调整相关的云物理和陆面过程等物理参数化方案是进一步提高SAMIL-R42L26辐射通量模拟性能的关键。

Water Vapor and Cloud Radiative Forcings over the Pacific Ocean Simulated by the LASG/IAP AGCM:Sensitivity to Convection Schemes

Characteristics of the total clear-sky greenhouse effect （GA） and cloud radiative forcings （CRFs）, along with the radiative-related water vapor and cloud properties simulated by the Spectral Atmospheric Model developed by LASGIAP （SAMIL） are evaluated. Impacts of the convection scheme on the simulation of CRFs are discussed by using two AMIP （Atmospheric Model Inter-comparison Project） type simulations employing different convection schemes： the new Zhang-McFarlane （NZH） and Tiedtke （TDK） convection schemes. It shows that both the climatological GA and its response to El Nio warming are simulated well, both in terms of spatial pattern and magnitude. The impact of the convection scheme on GA is not significant. The climatological longwave CRF （LWCRF） and its response to El Nio warming are simulated well, but with a prominently weaker magnitude. The simulation of the climatology （response） of LWCRF in the NZH （TDK） run is slightly more realistic than in the TDK （NZH） simulation, indicating significant impacts of the convection scheme. The shortwave CRF （SWCRF） shows large biases in both spatial pattern and magnitude, and the results from the TDK run are better than those from the NZH run. A spuriously excessive negative climatological SWCRF over the southeastern Pacific and an insufficient response of SWCRF to El Nio warming over the tropical Pacific are seen in the NZH run. These two biases are alleviated in the TDK run, since it produces vigorous convection, which is related to the low threshold for convection to take place. Also, impacts of the convection scheme on the cloud profile are discussed.

Simulations of the East Asian Subtropical Westerly Jet by LASG/IAP AGCMs

Performances of two LASG/IAP （State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics） Atmospheric General Circulation Models （AGCMs）, namely GAMIL and SAMIL, in simulating the major characteristics of the East Asian subtropical westerly jet （EASWJ） in the upper troposphere are examined in this paper. The mean vertical and horizontal structures and the correspondence of the EASWJ location to the meridional temperature gradient in the upper troposphere are well simulated by two models. However, both models underestimate the EASWJ intensity in winter and summer, and are unable to simulate the bimodal distribution of the major EASWJ centers in mid-summer, relative to the observation, especially for the SAMIL model. The biases in the simulated EASWJ intensity are found to be associated with the biases of the meridional temperature gradients in the troposphere, and furthermore with the surface sensible heat flux and condensation latent heating. The models capture the major characteristics of the seasonal evolution of the diabatic heating rate averaged between 30°-45°N, and its association with the westerly jet. However, the simulated maximum diabatic heating rate in summer is located westward in comparison with the observed position, with a relatively strong diabatic heating intensity, especially in GAMIL. The biases in simulating the diabatic heating fields lead to the biases in simulating the temperature distribution in the upper troposphere, which may further affect the EASWJ simulations. Therefore, it is necessary to improve the simulation of the meridional temperature gradient as well as the diabatic heating field in the troposphere for the improvement of the EASWJ simulation by GAMIL and SAMIL models.

The Application of Flux-Form Semi-Lagrangian Transport Scheme in a Spectral Atmosphere Model

A flux-form semi-Lagrangian transport scheme （FFSL） was implemented in a spectral atmospheric GCM developed and used at IAP/LASG. Idealized numerical experiments show that the scheme is good at shape preserving with less dissipation and dispersion, in comparison with other conventional schemes, hnportantly, FFSL can automatically maintain the positive definition of the transported tracers, which was an underlying problem in the previous spectral composite method （SCM）. To comprehensively investigate the impact of FFSL on GCM results, we conducted sensitive experiments. Three main improvements resulted： first, rainfall simulation in both distribution and intensity was notably improved, which led to an improvement in precipitation frequency. Second, the dry bias in the lower troposphere was significantly reduced compared with SCM simulations. Third, according to the Taylor diagram, the FFSL scheme yields simulations that are superior to those using the SCM： a higher correlation between model output and observation data was achieved with the FFSL scheme, especially for humidity in lower troposphere. However, the moist bias in the middle and upper troposphere was more pronounced with the FFSL scheme. This bias led to an over-simulation of precipitable water in comparison with reanalysis data. Possible explanations, as well as solutions, are discussed herein.

The Wavenumber-Frequency Characteristics of the Tropical Waves in an Aqua-Planet GCM

Based on the aqua-planet experiments, the wavenumber-frequency characteristics of tropical waves and their influencing factors in SST distribution and the convective parameterization scheme are investigated using the spectral atmospheric general circulation model （SAMIL）. Space-time spectral analysis is used to obtain the variance of convectively coupled tropical waves. In the Control experiment with maximum SST located at the equator the simulated tropical-wave behaviors are in agreement with those in observations and theoretical solutions. When the maximum SST is located at 5°N, the symmetric and antisymmetric waves are much weaker than those in the control experiment, suggesting that tropical wave activities are very sensitive to the SST distributions. Importantly, the variance maximum of Madden-Julian oscillation （MJO） is found to occur around 5°N, which suggests that the development of the MJO depends largely on the latitude of maximum SST. Furthermore, the seasonal variations of MJO may be mainly caused by the seasonal variations of the maximum SST. The experiment results with two different cumulus schemes the Manabe moist convective adjustment and Zhang-McFarlane （ZM） convective scheme, were also compared to examine the impacts of convective parameterization. Weakened variances of each individual tropical wave when the ZM scheme is used suggest that the ZM scheme is not favorable for the tropical wave activities. However, the wave characteristics are different when the ZM scheme is used in different models, which may imply that the simulated basic state is important to the meridional distributions of the waves. The MJO signals suggest that the parameterization scheme may have great influence on the strength, but have less direct impact on the MJO distribution. The frequency of the tropical waves may be associated with the moisture control of convection and the large-scale condensation scheme used in the model.

An Updated Coupled Model for Land-Atmosphere Interaction.Part II:Simulations of Biological Processes

In Part I, the authors succeeded in coupling the spectral atmospheric model （SAMIL_R42L9） developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences （LASG/IAP/CAS） with the land surface model, Atmosphere-Vegetation-Interaction-Model （AVIM） and analyzed the climate basic state and land surface physical fluxes simulated by R42_AVIM. In this Part Ⅱ, we further evaluate the simulated results of the biological processes, including leaf area index （LAI）, biomass and net primary productivity （NPP） etc. Results indicate that R42_AVIM can simulate the global distribution of LAI and has good consistency with the monthly mean LAI provided by Max Planck Institute for Meteorology. The simulated biomass corresponds reasonably to the vegetation classifications. In addition, the simulated annual mean NPP has a consistent distribution with the data provided by IGBP and MODIS, and compares well with the work in literature. This land-atmosphere coupled model will offer a new experiment tool for the research on the two-way interaction between climate and biosphere, and the global terrestrial ecosystem carbon cycle.

An Updated Coupled Model for Land-Atmosphere Interaction.Part I:Simulations of Physical Processes

A new two-way land-atmosphere interaction model （R42_AVIM） is fulfilled by coupling the spectral atmospheric model （SAMIL_R42L9） developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences （LASG/IAP/CAS） with the land surface model, Atmosphere-Vegetation-Interaction-Model （AVIM）. In this coupled model, physical and biological components of AVIM are both included. Climate base state and land surface physical fluxes simulated by R42_AVIM are analyzed and compared with the results of R42_SSIB [which is coupled by SAMIL_R42L9 and Simplified Simple Biosphere （SSIB） models]. The results show the performance of the new model is closer to the observations. It can basically guarantee that the land surface energy budget is balanced, and can simulate June-July-August （JJA） and December-January- February （DJF） land surface air temperature, sensible heat flux, latent heat flux, precipitation, sea level pressure and other variables reasonably well. Compared with R42_SSIB, there are obvious improvements in the JJA simulations of surface air temperature and surface fluxes. Thus, this land-atmosphere coupled model will offer a good experiment platform for land-atmosphere interaction research.

Impact of Convective Momentum Transport by Deep Convection on Simulation of Tropical Intraseasonal Oscillation

This paper focuses on the impacts of convective momentum transport(CMT) on simulations of the tropical intraseasonal oscillation(TIO) in SAMIL. Two sets of experiments are performed, which give different reality of the Madden-Julian Oscillation(MJO). The Tiedtke cumulus parameterization scheme is used for all experiments. It is found that simulations of the TIO can be influenced by CMT, and the impacts on the simulated TIO depend on the model capability in simulating the MJO. CMT tends to have large influences to the model that can simulate the eastward propagation of the MJO. CMT can further influence the long-term mean of zonal wind and its vertical shear. Zonal wind suffers from easterlies biases at low level and westerlies biases at upper level when CMT is introduced. Such easterlies biases at low level reduce the reality of the simulated tropical intraseasonal oscillation. When CMT is introduced in the model, MJO signals disappear but the model's mean state improves. Therefore, a more appropriate way is needed to introduce CMT to the model to balance the simulated mean state and TIO signals.

两种对流参数化方案下降水和潜热加热空间结构的模拟及其影响

利用中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室(LASG/IAP)新版大气环流模式SAMIL(Spectral Atmosphere Model of IAP LASG),分析比较了Manabe和Tiedtke两种对流方案对降水和潜热加热空间结构的模拟性能.结果表明,Tiedtke方案模拟的对流降水偏强,层状降水偏弱,致使层状降水占总降水比例低于观测.与此相比,Manabe方案较为合理地再现了沿赤道的层状降水中心,但其模拟的孟加拉湾降水中心向东向北伸展,在高原下游有虚假强降水出现.通过比较两种对流方案模拟的潜热加热空间结构,结果表明Tiedtke方案的模拟结果远较Manabe方案合理,表现为与观测一致的位于对流层中层的对流加热峰值,以及0℃层之上的层状加热与其下的冷却.而Manabe方案模拟的对流加热峰值所在高度较观测偏低,大致位于700 hPa.另外对于层状加热,Manabe方案模拟不出0℃层以下的冷却,而表现为整层加热.温度、比湿等要素的Taylor图表明:Tiedtke方案的模拟结果与ERA40,NCEP再分析资料的相关系数更高,尤其是对流层高层200 hPa温度的模拟.文中同时提出Tiedtke方案下层状降水模拟失真的可能原因是层状降水方案忽略了积云对流过程中卷出水凝物的影响,该猜想通过敏感性试验得到证实.

Features of rainfall and latent heating structure simulated by two convective parameterization schemes

Using the latest version of SAMIL (Spectral Atmosphere Model of IAP LASG) developed by LASG/IAP,we evaluate the model performance by analyzing rainfall,latent heating structure and other basic fields with two different convective parameterization schemes:Manabe Scheme and Tiedtke Scheme.Results show that convective precipitation is excessively overestimated while stratiform precipitation is underestimated by Tiedtke scheme,thus causing less stratiform rainfall proportion compared with TRMM observation.In contrast,for Manabe scheme stratiform rainfall belt is well simulated,although precipitation center near Bay of Bengal (BOB) spreads eastward and northward associated with unrealistic strong rainfall downstream of the Tibet Plateau.The simulated latent heating structure indicates that Tiedtke scheme has an advantage over Manabe scheme,as the maximum convective latent heating near middle of troposphere is well reproduced.Moreover,the stratiform latent heating structure is also well simulated by Tiedtke scheme with warming above freezing level and cooling beneath freezing level.As for Manabe scheme,the simulated maximum convective latent heating lies near 700 hPa,lower than the observation.Additionally,the warming due to stratiform latent heating extends to the whole vertical levels,which is unreasonable compared with observation.Taylor diagram further indicates that Tiedtke scheme is superior to Manabe scheme as higher correlation between model output and observation data is achieved when Tiedtke scheme is employed,especially for the temperature near 200 hPa.Finally,a possible explanation is addressed for the unrealistic stratiform rainfall by Tiedtke scheme,which is due to the neglect of detrained cloud water and cloud ice during convective process.The speculation is verified through an established sensitivity experiment.

Development of a toy column model and its application in testing cumulus convection parameterizations

A single-column model is constructed based on parameterizations inherited from the Finite-volume/Spectral Atmospheric Model F/SAMIL and tested in simulations of tropical convective systems. Two representative convection schemes are compared in terms of their performances on precipitation types, individual physical tendencies, and temperature and moisture fields. The main difference between the two selected schemes is in their representation of entraining/detraining process. The Tiedtke scheme assumes bulk entrainment, while the Zhang–Mc Farlane scheme parameterizes entrainment/detrainment rates under the spectrum concept. Large-scale forcing and verification data are taken from the GATE phase III field campaign, during which abundant convective events were observed. Given the same triggering function and closure assumption, results show that entrainment/detrainment representation remains the dominant factor on the simulation of cumulus mass flux and of temperature and moisture fields. By analyzing sources and sinks of heat and moisture, this study reveals how parameterization components compensate for each other and make model results insensitive to parameterization changes in certain fields, thus suggesting the need to treat parameterizations as systems rather than individual components.