一个海气耦合模式模拟的云辐射过程

利用NCC/IAP T63海气耦合模式进行了20 a积分,详细分析了模式对云量及其辐射影响的模拟能力。结果表明,模式能够模拟出云量分布的基本特征,但同ISCCP卫星观测资料及ERA再分析资料相比还存在较大的差距,总体表现为模拟的云量偏小,尤其是海洋上部分地区出现了异常的低值区。通过对云量方案的改进,明显改善了两大洋东岸、夏半球副热带到中纬度海洋上空低云的模拟。但模式对热带印度洋到西太平洋地区云量的模拟仍然存在明显的偏差,这主要是由于模式对该地区强对流云模拟能力差,造成该地区高云模拟存在较大的误差。对辐射及其云辐射强迫的分析表明,模式对长波及其云辐射强迫的模拟要明显好于短波。短波辐射模拟的偏差主要是由于短波云辐射强迫模拟过小、耦合模式对积雪和海冰模拟较差、以及未考虑气溶胶的影响等原因共同引起的;而长波辐射模拟的差距主要是云量以及下垫面温度模拟不足造成的。相应于云量方案的改进,两大洋东岸、夏半球副热带到中纬度海洋上辐射(尤其是短波辐射)的模拟有了明显的改善,这也明显改进了这些地区的净辐射模拟。     

中国东部和印度季风区云辐射特性的比较

基于 ＩＳＣＣＰ和 ＥＱＢＥ资料,本文比较了中国东部和印度季风区的云和云辐射强迫的气候特征。虽然它们同属于亚洲季风区,并且有相似的降水季节特征,但它们各自的云和云辐射强迫特征差异很大。在印度区域,所有的云量有着相同的季节变化,最大云量分布都出现在夏季,且总云量中以高云量为主。而中国东部云量的季节变化都比较复杂,在总云量中以中、低云量为主,最大总云量出现在春季。冬季的总云量和中、低云量要大于夏季。在全球云量分布中,中国东部最典型的特征是：该地区为全球最大的雨层云覆盖区。与云的分布和变化相关,印度季风区最大的负短波云辐射强迫,最大的正的长波辐射强迫和最大的负的净云辐射强迫发生在夏季,而在中国东部,大的负的短波云辐射强迫发生在春夏之交。年平均的负的短波云辐射强迫在中国东部地区明显要大于在印度季风区。     

副热带东太平洋低云对LASG耦合环流模式FGCM-0中"双ITCZ"的影响

与其他耦合环流模式一样,LASG耦合模式FGCM-0也存在虚假的“双ITCZ”。为了认识FGCM-0中“双ITCZ”,首先研究了FGCM-0的大气分量模式,剧INCAR(美国国家大气研究中心)的公用气候模式CCM3对秘鲁和加利福尼亚沿岸低云以及低层大气整体稳定度的模拟能力。发现：尽管CCM3模拟的低层大气整体稳定度与利用NCEP(美国国家环境预报中心)再分析资料分析的结果较一致,但模拟的低云量比ISCCP(国际卫星云气候计划)观测值显偏少。利用ISCCP低云量与由NCEP再分析温度场分析的低层整体稳定度之间的回归关系,修改了CCM3中低云参数化方案,并用于敏感性试验,以研究副热带东太平洋低云对FGCM-0中“双ITCZ”的影响。结果发现,修改的方案能显增强对低云量的模拟,秘鲁沿岸冷海域低云量增加能显减弱赤道以南热带东太平洋海表面温度(SST)的暖偏差,但同时也将使赤道冷舌增强、向西伸展更远;加利福尼亚沿岸低云量增加可以有效减弱赤道以北ITCZ区SST暖偏差。为了检验秘鲁沿岸SST与低云间的正反馈,又实施了一个控制秘鲁沿岸SST的敏感性试验,结果表明：控制秘鲁沿岸SST抑制其增暖,对自东南太平洋向西北至中、西赤道太平洋广大区域产生的影响,与增加秘鲁沿岸低云量产生的影响相似。     

ENERGY BUDGET BIAS IN GLOBAL COUPLED OCEAN-ATMOSPHERE-LAND MODEL*

The energy budget of the two versions of the GOALS model (GOALS-1.1 and GOALS-2) is described and compared to observational estimates.The results illustrate that the simulated surface net shortwave radiation flux is underestimated in the high-latitude regions while the surface net longwave radiation flux is substantially overestimated in that region,which results in the lower surface air temperature (SAT) of the polar region and the stronger negative sensible heat flux in high latitudes.The overestimated sensible heat flux from surface to atmosphere in the continents causes the much warmer SAT centers,which may be the reason for the bias of the model SAT. The bias that the simulated precipitation is less than observation in most regions is closely related to the underestimated latent heat flux over most of the Eurasian Continent and the oceans, especially over the subtropical oceans.It can be seen that the bias in the OLR of the two models lies in low and middle latitudes,where the absorbed solar shortwave radiation flux at the top of the atmosphere is comparable to the NCEP reanalysis,but much less than ERBE data.This indicates that the improvement of cloud-radiation parameterization scheme in low and middle latitudes is of critical importance to the simulation of global energy budget.The simulated cloud cover from the GOALS-2 model with diagnosed cloud scheme is generally less except at equatorial areas, especially in the mid-latitude areas,which causes the large bias of energy budget there.It is suggested that the refinement of cloud parameterization is one of the most important tasks in the model's future development.     

A hybrid technique for computing the monthly mean net longwave surface radiation over oceanic areas

Summary Global maps of the monthly mean net upward longwave radiation flux at the ocean surface have been obtained for April, July, October 1985 and January 1986. These maps were produced by blending information obtained from a combination of general circulation model cloud radiative forcing fields, the top-of-the-atmosphere cloud radiative forcing from ERBE and TOVS profiles and sea surface temperature on ISCCP C1 tapes. The fields are compatible with known meteorological regimes of atmospheric water vapor content and cloudiness. There is a vast area of high net upward longwave radiation flux (> 80 W m^–2) in the eastern Pacific Ocean throughout most of the year. Areas of low net upward longwave radiation flux (< 40 Wm^–2) are the tropical convective regions and extra tropical regions that tend to have persistent low cloud cover. The technique used in this study relies on GCM simulations and so is subject to some of the uncertainties associated with the model. However, all input information regarding temperature, moisture and cloud cover is from satellite data having near global coverage. This feature of the procedure alone warrants its consideration for further use in compiling global maps of the net longwave radiation at the surface over the oceans.With 9 Figures     

Comments on ``Direct Radiative Forcing of Anthropogenic Aerosols over Oceans from Satellite Observation"

Previous observational studies have estimated anthropogenic aerosol direct radiative forcing over oceans without due consideration of cloudy-sky aerosols.However,when interaction between clouds and aerosols located below or above the cloud level is taken into account,the aerosol direct radiative forcing is larger by as much as 5 W m-2 in most mid-latitude regions in the Northern Hemisphere.     

Evaluation of cloud vertical structure simulated by recent BCC_AGCM versions through comparison with CALIPSO-GOCCP data

ABSTRACT The abilities of BCC-AGCM2.1 and BCC_AGCM2.2 to simulate the annual-mean cloud vertical structure （CVS） were evaluated through comparison with GCM-Oriented CALIPSO Cloud Product （CALIPSO-GOCCP） data. BCC-AGCM2.2 has a dynamical core and physical processes that are consistent with BCC-AGCM2.1, but has a higher horizontal resolution. Results showed that both BCC-AGCM versions underestimated the global-mean total cloud cover （TCC）, middle cloud cover （MCC） and low cloud cover （LCC）, and that BCC_AGCM2.2 underestimated the global-mean high cloud cover （HCC）. The global-mean cloud cover shows a systematic decrease from BCCA-GCM2.1 to BCC_AGCM2.2, especially for HCC. Geographically, HCC is significantly overestimated in the tropics, particularly by BCC_AGCM2,1, while LCC is generally overestimated over extra-tropical lands, but significantly underestimated over most of the oceans, especially for subtropical marine stratocumulus clouds. The leading EOF modes of CVS were extracted. The BCC_AGCMs perform well in reproducing EOF1, but with a larger variance explained. The two models also capture the basic features of EOF3, except an obvious deficiency in eigen- vector peaks. EOF2 has the largest simulation biases in both position and strength of eigenvector peaks. Furthermore, we investigated the effects of CVS on relative shortwave and longwave cloud radiative forcing （RSCRF and RLCRF）. Both BCC_AGCM versions successfully reproduce the sign of regression coefficients, except for RLCRF in PC1. However, the RSCRF relative contributions from PC1 and PC2 are overestimated, while the relative contribution from PC3 is underes timated in both BCC_AGCM versions. The RLCRF relative contribution is underestimated for PC2 and overestimated for PC3.     

NCC/IAP T63海气耦合模式对云滴有效半径的敏感性分析

利用NCC/IAP T63海气耦合模式研究了模拟对云滴有效半径的敏感性。结果表明,可变水云粒子有效半径的引入（REL和REL_CAM3方案）导致对流层中下部云消光光学厚度的明显变化,但REL和REL_CAM3之间的变化位置不同、符号相反;而引入可变的冰云粒子有效半径后（REI方案）,对流层中上部的云消光光学厚度明显减小。云辐射强迫的变化主要受REI方案的影响,长波和短波云辐射强迫均以减小为主;从全球平均来看,大气顶云辐射强迫的变化对冰云有效半径的变化更加敏感,尽管量级还有明显差距,各种方案均正确模拟出了大气顶净云辐射强迫的符号。大气顶净辐射的变化与净云辐射强迫的变化基本一致,总体以减少为主。对于地表气温的变化,REL和REI方案总体上有着较为一致的分布,夏季以降温为主,冬季的欧亚大陆中部和北美北部增温明显。此外,大陆地表气温的变化主要受地表热量输送的影响,而热带海洋上空却表现出明显的不一致性。同时引入可变的水云和冰云粒子（RIW和RIW_CAM3方案）,云的消光光学厚度和辐射强迫相对单独引入主要呈线性响应,但地表气温表现出明显的非线性特征,非线性最明显的地区主要位于大陆地区。       

Evaluating the Impacts of Cloud Microphysical and Overlap Parameters on Simulated Clouds in Global Climate Models

The improvement of the accuracy of simulated cloud-related variables, such as the cloud fraction, in global climate models (GCMs) is still a challenging problem in climate modeling. In this study, the influence of cloud microphysics schemes (one-moment versus two-moment schemes) and cloud overlap methods (observation-based versus a fixed vertical decorrelation length) on the simulated cloud fraction was assessed in the BCC_AGCM2.0_CUACE/Aero. Compared with the fixed decorrelation length method, the observation-based approach produced a significantly improved cloud fraction both globally and for four representative regions. The utilization of a two-moment cloud microphysics scheme, on the other hand, notably improved the simulated cloud fraction compared with the one-moment scheme; specifically, the relative bias in the global mean total cloud fraction decreased by 42.9%–84.8%. Furthermore, the total cloud fraction bias decreased by 6.6% in the boreal winter (DJF) and 1.64% in the boreal summer (JJA). Cloud radiative forcing globally and in the four regions improved by 0.3%?1.2% and 0.2%?2.0%, respectively. Thus, our results showed that the interaction between clouds and climate through microphysical and radiation processes is a key contributor to simulation uncertainty.     

Impacts of an improved low-level cloud scheme on the eastern Pacific ITCZ-cold tongue complex

A statistically-based low-level cloud parameterization scheme is introduced, modified, and applied in the Flexible coupled General Circulation Model (FGCM-0). It is found that the low-level cloud scheme makes improved simulations of low-level cloud fractions and net surface shortwave radiation fluxes in the subtropical eastern oceans off western coasts in the model. Accompanying the improvement in the net surface shortwave radiation fluxes, the simulated distribution of SSTs is more reasonably asymmetrical about the equator in the tropical eastern Pacific, which suppresses, to some extent, the development of the double ITCZ in the model. Warm SST biases in the ITCZ north of the equator are more realistically reduced, too. But the equatorial cold tongue is strengthened and extends further westward, which reduces the precipitation rate in the western equatorial Pacific but increases it in the ITCZ north of the equator in the far eastern Pacific. It is demonstrated that the low-level cloud-radiation feedback would enhance the cooperative feedback between the equatorial cold tongue and the ITCZ. Based on surface layer heat budget analyses, it is demonstrated that the reduction of SSTs is attributed to both the thermodynamic cooling process modified by the increase of cloud fractions and the oceanic dynamical cooling processes associated with the strengthened surface wind in the eastern equatorial Pacific, but it is mainly attributed to oceanic dynamical cooling processes associated with the strengthening of surface wind in the central and western equatorial Pacific.     

A source of hail embryos

Abstract

The sensitivity of the annual cycle of ice cover in Baffin Bay to short‐wave radiation is investigated. The Princeton Ocean Model (POM) is used and is coupled with a multi‐category, dynamic‐thermodynamic sea‐ice model in which the surface energy balance governs the growth rates of ice of varying thickness. During spring and summer the short‐wave radiation flux dominates other surface heat fluxes and thus has the greatest effect on the ice melt. The sensitivity of model results to short‐wave radiation is tested using several, commonly used, shortwave parameterizations under climatological, as well as short‐term, atmospheric forcing. The focus of this paper is short‐term and annual variability. It is shown that simulated ice cover is sensitive to the short‐wave radiation formulation during the melting phase. For the Baffin Bay simulation, the differences in the resulting ice area and volume, integrated from May to November, can be as large as 45% and 70%, respectively. The parameterization of the effect of cloud cover on the short‐wave radiation can result in the sea‐ice area and volume changes reaching 20% and 30%, respectively. The variation of the cloud amount represents cloud data error, and has a relatively small effect (less then ±4%) on the simulated ice conditions. This is due to the fact that the effect of cloud cover on the short‐wave radiation flux is largely compensated for by its effect on the net near‐surface long‐wave radiation flux.     

辐射参数化方案对一个海气耦合模式云和辐射模拟的影响

比较Morcrette辐射方案和Fu_Liou辐射方案对NCC/IAP T63海气耦合模式云和辐射模拟的影响, 结果表明:两种方案模拟的大气顶入射辐射存在明显的差异; 晴空大气Fu_Liou方案的短波吸收能力在全球普遍较Morcrette方案低; 在60°S~60°N之间, Fu_Liou方案模拟的行星反照率更接近于ERBE卫星观测; 在对大气顶净辐射的模拟上, 除了冬季的太平洋和大西洋东岸云量明显减少的部分地区外, Fu_Liou方案对大气顶净辐射的模拟总体上较Morcrette方案有了较为明显的改善; Fu_Liou方案模拟的海洋低层云显著减少, 而热带地区高云的模拟明显增加; 由于采用了“二元云量”算法, 尽管云量有所减少, Fu_Liou方案模拟的云短波吸收作用仍有所增强, 一定程度上改进了Morcrette方案云的短波吸收作用偏弱的现象。     

BCC_AGCM2.1对中国东部地区云辐射特征模拟的偏差分析

通过与观测及再分析资料的对比,评估了中国国家气候中心大气环流模式BCC_AGCM 2.1对中国东部地区云辐射特征的模拟性能,并着重分析了模拟偏差的原因.在云辐射特征的基本气候态模拟方面,模式能大致再现中国东部中纬度层状云大值带,以及层状云冷季多、暖季少的季节特征,模拟的短波云辐射强迫也具有与观测相对应的季节变化特征.在云辐射强迫和地面温度相互影响过程的模拟方面,模式也能模拟出与观测相近的相互作用过程,即地面温度降低伴随着层状云云量增多以及负的净云辐射强迫加强,升温时层状云云量减少和净云辐射强迫减弱.但模式模拟的大陆层状云云量系统性偏少(尤其在冷季),使得模式在该处的短波云辐射强迫明显偏弱.初步分析表明,造成层状云模拟差异的主要原因是在中国西南地区对流层低层模式模拟的偏南气流明显偏弱以及陆-气潜热通量偏小.偏南气流偏弱导致低层散度和垂直运动条件不利于中层云的形成.同时偏南气流偏弱也不利于向西南地区的水汽输送,再加上模式模拟地表向上潜热通量偏小,这二者都使得模式模拟中国西南区域对流层低层的水汽含量严重偏少,相对湿度偏低,同样不利于层状云生成和发展.水汽偏少进一步导致在冷异常情况下青藏高原下游云辐射-地表温度反馈模拟偏弱,即呈现冷异常时,水汽条件偏弱限制了云量增加,弱化了进一步降低温度的反馈过程.     

层状云微物理参数反演及其辐射效应的个例研究

利用毫米波云雷达、微波辐射计联合反演方法,对2015年11月11日安徽寿县的一次层状云过程的云参数进行了反演,将所得云参数加入到SBDART辐射传输模式中,进行辐射通量计算,并将计算的地面辐射通量与观测的地面辐射通量进行了对比分析。研究表明:1)利用毫米波雷达和微波辐射计数据联合反演的云参数比较可靠;2)利用SBDART模式并结合反演的云参数,可以准确实时地计算地面及其他高度层的长短波辐射通量;3)在反演的云参数中,光学厚度对地面各种辐射通量的影响是最大的,云层的光学厚度越大,到达地面的太阳短波辐射越小,地面反射短波辐射也越小。另外云底温度越高,云体向下发射的红外长波辐射越大。地面向上的长波辐射是地面温度的普朗克函数,随地面温度而变;4)云对地面的短波辐射强迫为负值,对地面有降温的作用。云对地面的长波辐射强迫是一个正值,对地面有一个增温的作用;5)云对地面的净辐射强迫随时间变化很大,它的正负与太阳高度角和云参数有关。     

东亚地区云和地表反照率对硫酸盐直接辐射强迫的影响

文中利用区域气候模式深入探讨了东亚地区云及地表反照率对硫酸盐直接辐射强迫的影响 ,同时定量估算了东亚地区云区气溶胶的直接辐射强迫 ,讨论了硫酸盐对地表和大气短波辐射平衡产生的不同影响。研究表明 :云对气溶胶的直接辐射强迫具有很强的减弱作用 ,这种减弱作用不仅取决于云覆盖份数 ,而且取决于云的光学厚度。就区域平均而言 ,文中模拟的东亚地区气溶胶直接辐射强迫为 - 0 .0 97W /m2 ,占总直接辐射强迫的 10 .4%左右。表明云对硫酸盐直接辐射具有很强的减弱作用 ,在估算其总的直接辐射强迫时 ,云区的贡献不可忽视。较高的地表反照率会减弱硫酸盐的直接辐射强迫 ,而较低地表反照率则会增加硫酸盐的辐射强迫。硫酸盐气溶胶对大气辐射平衡影响非常小 ,但对地表辐射平衡产生重要影响 ,影响程度与大气几乎一致。     

Seasonal evolution of subtropical anticyclones in the climate system model FGOALS-s2

The simulation characteristics of the seasonal evolution of subtropical anticyclones in the Northern Hemisphere are documented for the Flexible Global Ocean-Atmosphere-Land Systemmodel, Spectral Version 2 (FGOALS-s2), developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, the Institute of Atmospheric Physics. An understanding of the seasonal evolution of the subtropical anticyclones is also addressed. Compared with the global analysis established by the European Centre for Medium-Range Forecasts, the ERA-40 global reanalysis data, the general features of subtropical anticyclones and their evolution are simulated well in both winter and summer, while in spring a pronounced bias in the generation of the South Asia Anticyclone(SAA) exists. Its main deviation in geopotential height from the reanalysis is consistent with the bias of temperature in the troposphere. It is found that condensation heating (CO) plays a dominant role in the seasonal development of the SAA and the subtropical anticyclone over the western Pacific (SAWP) in the middle troposphere. The CO biases in the model account for the biases in the establishment of the SAA in spring and the weaker strength of the SAA and the SAWP from spring to summer. CO is persistently overestimated in the central-east tropical Pacific from winter to summer, while it is underestimated over the area from the South China Sea to the western Pacific from spring to summer. Such biases generate an illusive anticyclonic gyre in the upper troposphere above the middle Pacific and delay the generation of the SAA over South Asia in April. In midsummer, the simulated SAA is located farther north than in the ERA-40 data owing to excessively strong surface sensible heating (SE) to the north of the Tibetan Plateau. Whereas, the two surface subtropical anticyclones in the eastern oceans during spring to summer are controlled mainly by the surface SE over the two continents in the Northern Hemisphere, which are simulated reasonably well, albeit with their centers shifted westwards owing to the weaker longwave radiation cooling in the simulation associated with much weaker local stratiform cloud. Further improvements in the related parameterization of physical processes are therefore identified.     

Cloud Variations under Subtropical High Conditions

The cloud variations under subtropical high(STH) conditions during summers over a ten-year period are studied using combined data from the International Satellite Cloud Climatology Project and the National Centers for Environmental Prediction.The results reveal that clouds mainly experience an isolated evolution in the STHs,which is designated in this study by the 1540 gpm geopotential lines at 850 hPa.In the STH domain throughout the Northern Hemisphere,the average amount of total clouds exceeds 30%.Low clouds dominate in the STH domain,contributing over 60%of total cloud amount within the Pacific subtropical high and over 40%within the Atlantic subtropical high.The prevalence of low clouds in above regions is determined by the circulation pattern around 150°-180°E and 850 hPa,which suppresses both the upward development of the cloud tops and the water vapor divergences near the surface.Furthermore,clouds present great geographical incoherence within the STH domain.In the eastern STHs,the amount of middle and low clouds increases to peak in the early morning and decreases to a trough in the afternoon,while the amount of high clouds remains stable throughout the day.Conversely,in the western STHs,the diurnal amplitude of low and middle clouds is less than three,while high clouds dramatically reach the maximum in the afternoon and drop to the minimum in the evening.Among the nine cloud categories,stratocumulus clouds with greater optical thickness account for the most under STH conditions,no matter their occurrence or amount,causing more shortwave cloud radiative forcing to cool the local atmosphere and surface as a consequence.     

基于RegCM4.6的中国地区云对地表太阳辐射的影响研究

使用RegCM4.6区域气候模式，选取Emanuel和Mix（Grell+Emanuel）两种积云对流参数化方案，以2016年为例，分别对中国地区云短波辐射强迫及其涉及的物理量进行数值模拟，揭示其时空分布特征，并探究两种积云对流参数方案模拟效果的差异及其原因。结果表明：从季节平均来看，全国地表云短波辐射强迫均为负值，云对地表为冷却效应，冬季最小，春夏季较大。塔里木盆地四季均为辐射强迫低值区，夏季冷却效应最弱，辐射强迫绝对值低于40 W·m-2；全天空地表净短波辐射分布也呈显著季节差异，除夏季外均呈由南向北逐渐递减的分布趋势；晴空地表净短波辐射在横断山脉处和塔里木盆地处均比较低，其中春季最为明显；两个方案所得的季节空间分布特征大致相同，但在数值上存在差异。春季时全国大部分地区全天空地表净短波辐射通量差异最大，在55 W·m-2左右，晴空地表净短波辐射通量在青藏高原处差异在60W·m-2左右。     

Surface currents and equatorial thermocline in a coupled upper ocean-atmosphere GCM

The Oregon State University coupled upper ocean-atmosphere GCM is evaluated in terms of the simulated winds, ocean currents and thermocline depth variations. Although the zonal wind velocities in the model are underestimated by a factor of about three and the zonal current velocities are underestimated by a factor of about five, the model is seen to qualitatively simulate the major features of the gyral scale currents, and the phases of the seasonal variation of the principal equatorial currents are in reasonable agreement with observations. The simulated tropical currents are dominated by Ekman transport and the eastern boundary currents do not penetrate far enough equatorward, while the western boundary currents do not penetrate far enough poleward. The subtropical trade wind belt and the mid-latitude westerlies are displaced equatorward of observations; hence, the mid-latitude eastward currents, principally the Kuroshio-North Pacific Drift and the Gulf Stream-North Atlantic Current are displaced equatorward. In spite of these shortcomings the surface current simulation of this two-layer upper ocean model is comparable with that of other ocean GCMs of coarse resolution. The coupled model successfully simulates the deepening of the thermocline westward across Pacific as a consequence of the prevailing Walker circulation. The region of most intense simulated surface forcing is located in the western Pacific due to a southwestward displacement of the northeast trade winds relative to observations; hence the equatorial Pacific is dominated by eastward propagation of thermocline depth variations. The excessively strong Ekman divergence and upwelling in the western Pacific cools the local warm pool, while incorrectly simulated westerlies in the eastern Pacific suppress upwelling and inhibit cooling from below. These features reduce the simulated trans-Pacific sea-surface temperature gradient, weakening the Walker circulation and the anomalies associated with the simulated Southern Oscillation. Offprint requests to: KR Sperber     

Annual cycle and interannual variability in the tropical pacific as simulated by three versions of FGOALS

The seasonal cycle and interannual variability in the tropical oceans simulated by three versions of the Flexible Ocean-Atmosphere-Land System (FGOALS) model (FGOALS-g1.0, FGOALS-g2 and FGOALSs2), which have participated in phases 3 and 5 of the Coupled Model Intercomparison Project (CMIP3 and CMIP5), are presented in this paper. The seasonal cycle of SST in the tropical Pacific is realistically reproduced by FGOALS-g2 and FGOALSs2, while it is poorly simulated in FGOALS-g1.0. Three feedback mechanisms responsible for the SST annual cycle in the eastern Pacific are evaluated. The ocean-atmosphere dynamic feedback, which is successfully reproduced by both FGOALS-g2 and FGOALS-s2, plays a key role in determining the SST annual cycle, while the overestimated stratus cloud-SST feedback amplifies the annual cycle in FGOALS-s2. Because of the serious warm bias existing in FGOALS-g1.0, the ocean-atmosphere dynamic feedback is greatly underestimated in FGOALS-g1.0, in which the SST annual cycle is mainly driven by surface solar radiation. FGOALS-g1.0 simulates much stronger ENSO events than observed, whereas FGOALS-g2 and FGOALSs2 successfully simulate the observed ENSO amplitude and period and positive asymmetry, but with less strength. Further ENSO feedback analyses suggest that surface solar radiation feedback is principally responsible for the overestimated ENSO amplitude in FGOALS-g1.0. Both FGOALS-g1.0 and FGOALS-s2 can simulate two different types of El Ni-no events — with maximum SST anomalies in the eastern Pacific (EP) or in the central Pacific (CP) — but FGOALS-g2 is only able to simulate EP El Ni-no, because the negative cloud shortwave forcing feedback by FGOALS-g2 is much stronger than observed in the central Pacific.