Role of the upper ocean structure in the response of ENSO-like SST variability to global warming

The response of El Niño and Southern Oscillation (ENSO)-like variability to global warming varies comparatively between the two different climate system models, i.e., the Meteorological Research Institute (MRI) and Geophysical Fluid Dynamics Laboratory (GFDL) Coupled General Circulation Models (CGCMs). Here, we examine the role of the simulated upper ocean temperature structure in the different sensitivities of the simulated ENSO variability in the models based on the different level of CO₂ concentrations. In the MRI model, the sea surface temperature (SST) undergoes a rather drastic modification, namely a tendency toward a permanent El Niño-like state. This is associated with an enhanced stratification which results in greater ENSO amplitude for the MRI model. On the other hand, the ENSO simulated by GFDL model is hardly modified although the mean temperature in the near surface layer increases. In order to understand the associated mechanisms we carry out a vertical mode decomposition of the mean equatorial stratification and a simplified heat balance analysis using an intermediate tropical Pacific model tuned from the CGCM outputs. It is found that in the MRI model the increased stratification is associated with an enhancement of the zonal advective feedback and the non-linear advection. In the GFDL model, on the other hand, the thermocline variability and associated anomalous vertical advection are reduced in the eastern equatorial Pacific under global warming, which erodes the thermocline feedback and explains why the ENSO amplitude is reduced in a warmer climate in this model. It is suggested that change in stratification associated with global warming impacts the equatorial wave dynamics in a way that enhances the second baroclinic mode over the gravest one, which leads to the change in feedback processes in the CGCMs. Our results illustrate that the upper ocean vertical structure simulated in the CGCMs is a key parameter of the sensitivity of ENSO-like SST variability to global warming.     

Quantitative assessment of the climate components driving the pacific decadal oscillation in climate models

The Pacific decadal oscillation (PDO) is defined as the first empirical orthogonal function (EOF) mode of the North Pacific sea surface temperature anomalies. In this study, we reconstructed the PDO using the first-order autoregressive model from various climate indices representing the El Niño-Southern oscillation (ENSO), Aleutian Low (AL), sea surface height (SSH), and thermocline depth over the Kuroshio–Oyashio extension (KOE) region. The climate indices were obtained from observation and twentieth-century simulations of the eight coupled general circulation models (CGCMs) participating in the Climate Model Intercomparison Project Phase III (CMIP3). In this manner, we quantitatively assessed the major climate components generating the PDO using observation and models. Based on observations, the PDO pattern in the central to eastern North Pacific was accurately reconstructed by the AL and ENSO indices, and that in the western North Pacific was best reconstructed by the SSH and thermocline indices. In the CMIP3 CGCMs, the relative contribution of each component to the generation of the PDO varied greatly from model to model, and observations, although the PDO patterns from most of the models were similar to the pattern observed. In the models, the PDO pattern in the eastern and western North Pacific were well reconstructed using the AL and SSH indices, respectively. However, the PDO pattern reconstructed by the ENSO index was quite different from the observed pattern, which was possibly due to the model's common deficiency in simulating the amplitude and location of the ENSO. Furthermore, the differences in the contribution of the KOE thermocline index between the observed pattern and most of the models indicated that the PDO pattern associated with ocean wave dynamics is not properly simulated by most models. Therefore, the virtually well simulated PDO pattern by models is a result of physically inconsistent processes.     

On the subarctic North Atlantic cooling due to global warming

Various ocean reanalysis data reveal that the subarctic Atlantic sea surface temperature (SST) has been cooling during the twentieth century. A similar cooling pattern is found in the doubling CO₂ experiment obtained from the CMIP3 (coupled model intercomparison project third phase) compared to the pre-industrial experiment. Here, in order to investigate the main driver of this cooling, we perform the heat budget analysis on the subarctic Atlantic upper ocean temperature. The net surface heat flux associated with the increased concentration of greenhouse gases heats the subarctic ocean surface. In the most of models, the longwave radiation, latent heat flux, and sensible heat flux exert a warming effect, and the shortwave radiation exerts a cooling effect. On the other hand, the thermal advection by the meridional current reduces the subarctic upper ocean temperature in all models. This cold advection is attributed to the weakening of the meridional overturning circulation, which is related to the reduction in the ocean surface density. In particular, greater warming of the surface air than of the sea surface results in the reduction of surface evaporation and thereby enhanced freshening of the ocean surface water, while precipitation change was smaller than evaporation change. The thermal advections by both the wind-driven Ekman current and the density-driven geostrophic current contribute to cooling in most of the models, where the heat transport by the geostrophic current tends to be larger than that by the Ekman current.     

An OGCM Simulation of Seasonal and Interannual Variabilities in the Surface—Layer Pacific of the Equatorial Band

The heat budget is analyzed in the surface-layer (0-50 m) Pacific of the equatorial band (10°S-10°N),using the simulation of an ocean general circulation model from 1945 to 1993. The analysis indicates that downward net surface heat flux from the atmosphere and ocean advective heat fluxes play distinct roles in seasonal and interannual variabilities of surface-layer ocean temperature. The surface heat flux dominantly determines the ocean temperature in the seasonal time-scale. But, it has a negative feedback to the ocean temperature in the interannual time-scale. The interannual variability of ocean temperature is largely associated with the cold advection from off-equatorial divergent flow in the central Pacific and from upwelling in the cold tongue. Both the surface heat flux and ocean advective heat fluxes are important to the ocean temperature during an El Nino event. The ocean advective heat fluxes are further associated with local westward trade wind in the central Pacific. These results are largely consistent with some regional observational analyses.     

Interannual Salinity Variability in the Tropical Pacific in CMIP5 Simulations

Salinity variability and its causes in the tropical Pacific are analyzed using observations, reanalysis products and model simulations. The mixed-layer salinity(MLS) budget analyses from observations and reanalysis products indicate that its interannual evolution is closely related to ENSO and is predominantly governed by surface forcing and surface advection in the western-central equatorial Pacific. It is found that the observed MLS tendency leads Nin?o3.4 by about 12 months due to the effect of negative freshwater flux(evaporation minus precipitation). These observation-based analyses are used to evaluate the corresponding simulation using GFDL-ESM2 M. It is evident that the model can simulate the spatiotemporal variations of MLS with some discrepancies compared to observations. In the warm pool of the equatorial Pacific the MLS tendency in the model is sensitive to ocean dynamics, however model biases cause the tendency to be underestimated. In particular, the freshwater flux is overestimated while the ocean surface zonal current and vertical velocity at the base of the mixed layer are underestimated. Due to model biases in representing the related physics, the effects of surface forcing on the simulated MLS budget are overestimated and those of subsurface and surface advection are relatively weak. Due to weaker surface advection and subsurface forcing than observed, the simulated compensations for surface forcing are suppressed, and the simulated MLS tendency that leads Nin?o3.4 by 8–10 months, which is shorter than the observed lead time. These results are useful for the interpretation of observational analyses and other model simulations in the tropical Pacific.     

Interdecadal variability of the Pacific Ocean: model response to observed heat flux and wind stress anomalies

Variability of the Pacific Ocean is examined in numerical simulations with an ocean general circulation model forced by observed anomalies of surface heat flux, wind stress and turbulent kinetic energy (TKE) over the period 1970-88. The model captures the 1976-77 winter time climate shift in sea surface temperature, as well as its monthly, seasonal and longer term variability as evidenced in regional time series and empirical orthogonal function analyses. Examination of the surface mixed-layer heat budget reveals that the 1976-77 shift was caused by a unique concurrance of sustained heat flux input anomalies and very strong horizontal advection anomalies during a multi-month period preceding the shift in both the central Pacific region (where cooling occurred) and the California coastal region (where warming occurred). In the central Pacific, the warm conditions preceding and the cold conditions following the shift tend to be maintained by anomalous vertical mixing due to increases in the atmospheric momentum flux (TKE input) into the mixed layer (which deepens in the model after the shift) from the early 1970s to the late 1970s and 1980s. Since the ocean model does not contain feedback to the atmosphere and it succeeds in capturing the major features of the 1976-77 shift, it appears that the midlatitude part of the shift was driven by the atmosphere, although effects of midlatitude ocean-atmosphere feedback are still possible. The surface mixed-layer heat budget also reveals that, in the central Pacific, the effects of heat flux input and vertical mixing anomalies are comparable in amplitude while horizontal advection anomalies are roughly half that size. In the California coastal region, in contrast, where wind variability is much weaker than in the central Pacific, horizontal advection and vertical mixing effects on the mixed layer heat budget are only one-quarter the size of typical heat flux input anomalies.This paper was presented at the Second International Conference on Modelling of Global Climate Variability, held in Hamburg 7–11 September 1992 under the auspices of the Max Planck Institute for Meteorology. Guest Editor for these papers is L. Dümenil     

Ocean Response to a Climate Change Heat-Flux Perturbation in an Ocean Model and Its Corresponding Coupled Model

State-of-the-art coupled general circulation models(CGCMs)are used to predict ocean heat uptake(OHU)and sealevel change under global warming.However,the projections of different models vary,resulting in high uncertainty.Much of the inter-model spread is driven by responses to surface heat perturbations.This study mainly focuses on the response of the ocean to a surface heat flux perturbation F,as prescribed by the Flux-Anomaly-Forced Model Intercomparison Project(FAFMIP).The results of ocean model were compared with those of a CGCM with the same ocean component.On the global scale,the changes in global mean temperature,ocean heat content(OHC),and steric sea level(SSL)simulated in the OGCM are generally consistent with CGCM simulations.Differences in changes in ocean temperature,OHC,and SSL between the two models primarily occur in the Arctic and Atlantic Oceans(AA)and the Southern Ocean(SO)basins.In addition to the differences in surface heat flux anomalies between the two models,differences in heat exchange between basins also play an important role in the inconsistencies in ocean climate changes in the AA and SO basins.These discrepancies are largely due to both the larger initial value and the greater weakening change of the Atlantic meridional overturning circulation(AMOC)in CGCM.The greater weakening of the AMOC in the CGCM is associated with the atmosphere–ocean feedback and the lack of a restoring salinity boundary condition.Furthermore,differences in surface salinity boundary conditions between the two models contribute to discrepancies in SSL changes.     

The response of the North Pacific Decadal Variability to strong tropical volcanic eruptions

In this study, the effects of volcanic forcing on North Pacific climate variability, on interannual to decadal time scales, are examined using climate model simulations covering the last 600?years. The model used is the Bergen Climate Model, a fully coupled atmosphere–ocean general circulation model. It is found that natural external forcings, such as tropical strong volcanic eruptions (SVEs) and variations in total solar irradiance, play an important role in regulating North Pacific Decadal Variability (NPDV). In response to tropical SVEs the lower stratospheric pole–to–equator temperature gradient is enhanced. The North polar vortex is strengthened, which forces a significant positive Arctic Oscillation. At the same time, dipole zonal wind anomalies associated with strong polar vortex propagate downward from the lower stratosphere. Through positive feedbacks in the troposphere, the surface westerly winds across the central North Pacific are significantly weakened, and positive sea level pressure anomalies are formed in the North Pacific. This anomalous surface circulation results in changes in the net heat fluxes and the oceanic advection across the North Pacific. As a result of this, warm water converges in the subtropical western North Pacific, where the surface waters in addition are heated by significantly reduced latent and sensible heat fluxes from the ocean. In the eastern and high–latitude North Pacific the ocean loses more heat, and large–scale decreases in sea surface temperatures are found. The overall response of this chain of events is that the North Pacific enters a negative phase of the Pacific decadal oscillation (PDO), and this negative phase of the PDO is maintained for several years. It is thus concluded that the volcanic forcing plays a key role in the phasing of the PDO. The model results furthermore highlight the important role of troposphere–stratosphere coupling, tropical–extratropical teleconnections and extratropical ocean–atmosphere interactions for describing NPDV.     

The global heat balance: heat transports in the atmosphere and ocean

The heat budget has been computed locally over the entire globe for each month of 1988 using compatible top-of-the-atmosphere radiation from the Earth Radiation Budget Experiment combined with European Centre for Medium Range Weather Forecasts atmospheric data. The effective heat sources and sinks (diabatic heating) and effective moisture sources and sinks for the atmosphere are computed and combined to produce overall estimates of the atmospheric energy divergence and the net flux through the Earth's surface. On an annual mean basis, this is directly related to the divergence of the ocean heat transport, and new computations of the ocean heat transport are made for the ocean basins. Results are presented for January and July, and the annual mean for 1988, along with a comprehensive discussion of errors. While the current results are believed to be the best available at present, there are substantial shortcomings remaining in the estimates of the atmospheric heat and moisture budgets. The issues, which are also present in all previous studies, arise from the diurnal cycle, problems with atmospheric divergence, vertical resolution, spurious mass imbalances, initialized versus uninitialized atmospheric analyses, and postprocessing to produce the atmospheric archive on pressure surfaces. Over land, additional problems arise from the complex surface topography, so that computed surface fluxes are more reliable over the oceans. The use of zonal means to compute ocean transports is shown to produce misleading results because a considerable part of the implied ocean transports is through the land. The need to compute the heat budget locally is demonstrated and results indicate lower ocean transports than in previous residual calculations which are therefore more compatible with direct ocean estimates. A Poisson equation is solved with appropriate boundary conditions of zero normal heat flux through the continental boundaries to obtain the ocean heat transport. Because of the poor observational data base, adjustments to the surface fluxes are necessary over the southern oceans. Error bars are estimated based on the large-scale spurious residuals over land of 30 W m^–2 over 1000 km scales (10¹² m²). In the Atlantic Ocean, a northward transport emerges at all latitudes with peak values of 1.1±0.2 PW (1 standard error) at 20 to 30°N. Comparable values are achieved in the Pacific at 20°N, so that the total is 2.1±0.3 PW. The peak southward transport is at 15 to 20°S of 1.9±0.3 PW made up of strong components from both the Pacific and Indian Oceans and with a heat flux from the Pacific into the Indian Ocean in the Indonesian throughflow. The pattern of poleward heat fluxes is suggestive of a strong role for Ekman transports in the tropical regions.     

北太平洋海水表层铅直热通量收支的季节变化及与海温的联系

分析表明，北太平洋中纬度地区海水表层铅直热通量收支的分布特征与海流密切有关。暖流区的受热相对较小或失热相对较大，冷流区则反之。海水表层铅直热通量收支的季节变化分别具有一年、半年和四个月三种周期，并分别与太阳辐射、海流以及大气环流等相系。海温变化的滞后时间，基本场为２个月，扰动场为１个月。     

耦合海-气环流模式中双热带辐合带现象及其热收支分析

文中研究了耦合海-气环流模式中的双热带辐合带(Double ITCZ)现象,并对这一现象的成因从海洋热量收支的角度进行了初步分析。Double ITCZ现象是在热带太平洋赤道南北两侧各出现一条ITCZ的现象,这是耦合海-气环流模式中的较为普遍的一种异常现象,与实际气候中出现的Double ITCZ现象并非指同一问题。文中对比观测和模式结果,指出了Double ITCZ现象的主要特征,针对它的出现过程进行细致分析,再利用模式输出的热量收支各项进行统计,得到了从海洋热收支分析得到的海表温度变化原因。与观测到的正常模态相比,Double ITCZ是一个异常的模态,它的特征突出地表现为西太平洋暖池区的降温和东南太平洋10°S附近的升温。海洋热量收支分析表明,暖池区的降温主要是受到扩散的作用,而表层平流的异常增暖在决定异常辐合带区升温过程中占第一位的作用。     

Heat budget of the south-central equatorial Pacific in CMIP3 models

ABSTRACT Using data from 17 coupled models and nine sets of corresponding Atmospheric Model Intercomparison Project （AMIP） results, we investigated annual and seasonal variation biases in the upper 50 m of the south-central equatorial Pacific, with a focus on the double-ITCZ bias, and examined the causes for the amplitude biases by using heat budget analysis. The results showed that, in the research region, most of the models simulate SSTs that are higher than or similar to observed. The simulated seasonal phase is close to that observed, but the amplitudes of more than half of the model results are larger than or equal to observations. Heat budget analysis demonstrated that strong shortwave radiation in individual atmospheric models is the main factor that leads to high SST values and that weak southward cold advection is an important mechanism for maintaining a high SST. For seasonal circulation, large surface shortwave radiation amplitudes cause large SST amplitudes.     

The Impact of Global Warming on the Pacific Decadal Oscillation and the Possible Mechanism简

The response of the Pacific Decadal Oscillation （PDO） to global warming according to the Fast Ocean Atmosphere Model （FOAM） and global warming comparison experiments of 11 IPCC AR4 models is investigated. The results show that North Pacific ocean decadal variability, its dominant mode （i.e., PDO）, and atmospheric decadal variability, have become weaker under global warming, but with PDO shifting to a higher frequency. The SST decadal variability reduction maximum is shown to be in the subpolar North Pacific Ocean and western North Pacific （PDO center）. The atmospheric decadal variability reduction maximum is over the PDO center. It was also found that oceanic baroclinic Rossby waves play a key role in PDO dynamics, especially those in the subpolar ocean. As the frequency of ocean buoyancy increases under a warmer climate, oceanic baroclinic Rossby waves become faster, and the increase in their speed ratio in the high latitudes is much larger than in the low latitudes. The faster baroclinic Rossby waves can cause the PDO to shift to a higher frequency, and North Pacific decadal variability and PDO to become weaker.     

Role of horizontal density advection in seasonal deepening of the mixed layer in the subtropical Southeast Pacific

The mechanisms behind the seasonal deepening of the mixed layer(ML) in the subtropical Southeast Pacific were investigated using the monthly Argo data from 2004 to 2012. The region with a deep ML(more than 175 m) was found in the region of(22?–30?S, 105?–90?W), reaching its maximum depth(～200 m) near(27?–28?S, 100?W) in September. The relative importance of horizontal density advection in determining the maximum ML location is discussed qualitatively. Downward Ekman pumping is key to determining the eastern boundary of the deep ML region. In addition, zonal density advection by the subtropical countercurrent(STCC) in the subtropical Southwest Pacific determines its western boundary, by carrying lighter water to strengthen the stratification and form a "shallow tongue" of ML depth to block the westward extension of the deep ML in the STCC region. The temperature advection by the STCC is the main source for large heat loss from the subtropical Southwest Pacific. Finally, the combined effect of net surface heat flux and meridional density advection by the subtropical gyre determines the northern and southern boundaries of the deep ML region: the ocean heat loss at the surface gradually increases from 22?S to 35?S, while the meridional density advection by the subtropical gyre strengthens the stratification south of the maximum ML depth and weakens the stratification to the north. The freshwater flux contribution to deepening the ML during austral winter is limited. The results are useful for understanding the role of ocean dynamics in the ML formation in the subtropical Southeast Pacific.     

水平非均匀对流边界层热量平衡和平流输送作用的大涡模拟

采用大涡模拟所获的数据结果，分析地面热通量沿平均风方向存在 突变的条件下对流边界层的热量平衡和平流输送作用。分析表明边界层内模拟所得结果 可以很好地满足热量平衡关系。除地面热通量项以外，平流项（包括水平平流和垂直平 流）对边界层加热率的作用可达地面热通量不均匀性差值的大小，是影响边界层内热量 平衡的最重要因子，平均速度散度项对热量平衡的作用也不可忽略，但湍流通量散度项 的作用则很小。     

Effects of tropical instability wave (TIW)-induced surface wind feedback in the tropical Pacific Ocean

Tropical instability waves (TIWs) arise from oceanic instability in the eastern tropical Pacific and Atlantic Oceans, having a clear atmospheric signature that results in coupled atmosphere–ocean interactions at TIW scales. In this study, the extent to which TIW-induced surface wind feedback influences the ocean is examined using an ocean general circulation model (OGCM). The TIW-induced wind stress (τ_TIW) part is diagnostically determined using an empirical τ_TIW model from sea surface temperature (SST) fields simulated in the OGCM. The interactively represented TIW wind tends to reduce TIW activity in the ocean and influence the mean state, with largest impacts during TIW active periods in fall and winter. In December, the interactive τ_TIW forcing induces a surface cooling (an order of ?0.1 to ?0.3 °C), an increased heat flux into the ocean, a shallower mixed layer and a weakening of the South Equatorial Current in the eastern equatorial Pacific. Additionally, the TIW wind effect yields a pronounced latitudinal asymmetry of sea level field across the equator, and a change to upper thermal structure, characterized by a surface cooling and a warming below in the thermocline, leading to a decreased temperature gradient between the mixed layer and the thermocline. Processes responsible for the τ_TIW–induced cooling effects are analyzed. Vertical mixing and meridional advection are the two terms in the SST budget that are dominantly affected by the TIW wind feedback: the cooling effect from the vertical mixing on SST is enhanced, with the maximum induced cooling in winter; the warming effect from the meridional advection is reduced in July–October, but enhanced in November–December. Additional experiments are performed to separate the relative roles the affected surface momentum and heat fluxes play in the cooling effect on SST. This ocean-only modeling work indicates that the effect of TIW-induced wind feedback is small but not negligible, and may need to be adequately taken into account in large-scale climate modeling.     

A model study on oceanic processes during the Indian Ocean Dipole termination

In this study, the processes affecting the temperature variability over the Southeastern Tropical Indian Ocean (STIO) during 1958–2000, accomplishing the positive and negative Indian Ocean Dipole (IOD) events are analyzed. The upper ocean heat budget analysis of the STIO has been carried out to understand the oceanic process during the termination of the recent strongest IOD events. The three recent strongest positive IOD events (1961, 1994 and 1997) and a strong negative IOD event (1996) are considered for detailed analysis. The heat budget analysis revealed that the positive net-surface heat flux and vertical advection played dominant roles in the termination of 1997 IOD event, whereas horizontal and vertical advections are responsible for the termination of IOD events during 1961 and 1994. The anomalous negative surface heat flux and horizontal advection caused the dipole termination during the negative dipole year 1996. The findings are well supported by the analysis of anomaly correlation between model upper ocean heat content tendency and heat budget components. Significant intra-seasonal oscillations (ISOs) in sea surface temperature (SST) anomaly are seen during the initial phase of termination in the eastern equatorial Indian Ocean during 1961 and 1994 IOD events. The influence of ISOs in SST is not so evident during the IOD termination in 1997. It is found that the termination processes have started more than a month prior to the actual IOD termination.     

On the seasonal sea surface temperature variations in the tropical Pacific

Summary In this paper, we investigated physical processes that control the seasonal variations of sea surface temperature in the tropical Pacific, using an intermediate ocean model. It is found that the westward propagation of sea surface temperature along the equator is attributed to dynamic response of the ocean to the wind (that consists of 3-dimensional temperature advection), whereas the northward propagation of sea surface temperature in the eastern Pacific results from the thermodynamic response of the ocean to the surface heat flux, primarily shortwave radiation that includes the effect of low-level stratus clouds. The remote response of the eastern Pacific sea surface temperature to seasonally varying wind in the western Pacific is of secondary importance, compared to the local wind forcing. The results suggest that the mechanism that controls the seasonal cycle of sea surface temperature is different from that associated with El Nino-Southern Oscillation.With 9 Figures     

ENSO representation in climate models: from CMIP3 to CMIP5

We analyse the ability of CMIP3 and CMIP5 coupled ocean–atmosphere general circulation models (CGCMs) to simulate the tropical Pacific mean state and El Niño-Southern Oscillation (ENSO). The CMIP5 multi-model ensemble displays an encouraging 30 % reduction of the pervasive cold bias in the western Pacific, but no quantum leap in ENSO performance compared to CMIP3. CMIP3 and CMIP5 can thus be considered as one large ensemble (CMIP3 + CMIP5) for multi-model ENSO analysis. The too large diversity in CMIP3 ENSO amplitude is however reduced by a factor of two in CMIP5 and the ENSO life cycle (location of surface temperature anomalies, seasonal phase locking) is modestly improved. Other fundamental ENSO characteristics such as central Pacific precipitation anomalies however remain poorly represented. The sea surface temperature (SST)-latent heat flux feedback is slightly improved in the CMIP5 ensemble but the wind-SST feedback is still underestimated by 20–50 % and the shortwave-SST feedbacks remain underestimated by a factor of two. The improvement in ENSO amplitudes might therefore result from error compensations. The ability of CMIP models to simulate the SST-shortwave feedback, a major source of erroneous ENSO in CGCMs, is further detailed. In observations, this feedback is strongly nonlinear because the real atmosphere switches from subsident (positive feedback) to convective (negative feedback) regimes under the effect of seasonal and interannual variations. Only one-third of CMIP3 + CMIP5 models reproduce this regime shift, with the other models remaining locked in one of the two regimes. The modelled shortwave feedback nonlinearity increases with ENSO amplitude and the amplitude of this feedback in the spring strongly relates with the models ability to simulate ENSO phase locking. In a final stage, a subset of metrics is proposed in order to synthesize the ability of each CMIP3 and CMIP5 models to simulate ENSO main characteristics and key atmospheric feedbacks.     

Versions g1.0 and g1.1 of the LASG/IAP Flexible Global Ocean–Atmosphere–Land System Model

The latest two versions of the IAP Flexible Global Ocean-Atmosphere-Land System (FGOALS) model- versions g1.0 and g1.1, are described in this study. Both two versions are fully coupled GCMs without any flux correction, major changes for g1.1 mainly lie in four aspects: (1) advection schemes for tracer in the ocean component model; (2) zonal filter scheme in high latitudes in the ocean component model; (3) coupling scheme for fresh water flux in high latitudes; and (4) an improved algorithm of airsea turbulent flux depending on the surface current of the ocean. As a result, the substantial cold biases in the tropical Pacific and high latitudes are improved by g1.1, especially g1.1 simulates more reasonable equatorial thermocline, poleward heat transport, zonal overturning stream function in the ocean and sea ice distribution than g1.0. Significant ENSO variability are simulated by both versions, however the ENSO behavior by g1.0 differs from the observed one in many aspects: about twice ENSO amplitude as observed, false ENSO asymmetry, only one peak period around 3 years, etc. Due to improved mean climate state by g1.1, many basic characteristics of ENSO are reproduced by g1.1, e.g., more reasonable ENSO amplitude, two peaks of power spectra for ENSO events, and positive SST skewness in the eastern Pacific as observed.