On the dependence of ENSO simulation on the coupled model mean state

Systematic model error remains a difficult problem for seasonal forecasting and climate predictions. An error in the mean state could affect the variability of the system. In this paper, we investigate the impact of the mean state on the properties of ENSO. A set of coupled decadal integrations have been conducted, where the mean state and its seasonal cycle have been modified by applying flux correction to the momentum-flux and a combination of heat and momentum fluxes. It is shown that correcting the mean state and the seasonal cycle improves the amplitude of SST inter-annual variability and also the penetration of the ENSO signal into the troposphere and the spatial distribution of the ENSO teleconnections are improved. An analysis of a multivariate PDF of ENSO shows clearly that the flux correction affects the mean, variance, skewness and tails of the distribution. The changes in the tails of the distribution are particularly noticeable in the case of precipitation, showing that without the flux correction the model is unable to reproduce the frequency of large events. For the inter-annual variability the momentum-flux correction alone has a large impact, while the additional heat-flux correction is important for the teleconnections. These results suggest that the current forecasts practices of removing the forecast bias a-posteriori or anomaly initialisation are by no means optimal, since they can not deal with the strong nonlinear interactions. A consequence of the results presented here is that the predictability on annual time-ranges could be higher than currently achieved. Whether or not the correction of the model mean state by some sort of flux correction leads to better forecasts needs to be addressed. In any case, flux correction may be a powerful tool for diagnosing coupled model errors and predictability studies.     

Sensitivity of ENSO characteristics to a new interactive flux correction scheme in a coupled GCM

A fast coupled global climate model (CGCM) is used to study the sensitivity of El Ni?o Southern Oscillation (ENSO) characteristics to a new interactive flux correction scheme. With no flux correction applied our CGCM reveals typical bias in the background state: for instance, the cold tongue in the tropical east Pacific becomes too cold, thus degrading atmospheric sensitivity to variations of sea surface temperature (SST). Sufficient atmospheric sensitivity is essential to ENSO. Our adjustment scheme aims to sustain atmospheric sensitivity by counteracting the SST drift in the model. With reduced bias in the forcing of the atmosphere, the CGCM displays ENSO-type variability that otherwise is absent. The adjustment approach employs a one-way anomaly coupling from the ocean to the atmosphere: heat fluxes seen by the ocean are based on full SST, while heat fluxes seen by the atmosphere are based on anomalies of SST. The latter requires knowledge of the model??s climatological SST field, which is accumulated interactively in the spin-up phase (??training??). Applying the flux correction already during the training period (by utilizing the evolving SST climatology) is necessary for efficiently reducing the bias. The combination of corrected fluxes seen by the atmosphere and uncorrected fluxes seen by the ocean implies a restoring mechanism that counteracts the bias and allows for long stable integrations in our CGCM. A suite of sensitivity runs with varying training periods is utilized to study the effect of different levels of bias in the background state on important ENSO properties. Increased duration of training amplifies the coupled sensitivity in our model and leads to stronger amplitudes and longer periods of the Nino3.4 index, increased emphasis of warm events that is reflected in enhanced skewness, and more pronounced teleconnections in the Pacific. Furthermore, with longer training durations we observe a mode switch of ENSO in our model that closely resembles the observed mode switch related to the mid-1970s ??climate shift??.     

The influence of systematic errors in the Southeast Pacific on ENSO variability and prediction in a coupled GCM

The impact of the warm SST bias in the Southeast Pacific (SEP) on the quality of seasonal and interannual variability and ENSO prediction in a coupled GCM is investigated. The reduction of this bias is achieved by means of empirical heat flux correction that is constant in time. It leads to a wide range of changes in the tropical Pacific climate including enhanced southeast trades, well-defined dry zone in the SEP, better simulation of the South Pacific Convergence Zone and stronger cross-equatorial asymmetry of the mean state in the eastern Pacific. As a result of the mean climate correction, significant improvements in the simulation of the seasonal cycle of the oceanic and atmospheric states are also observed both at the equator and basin-wide. Due to more realistic simulation of the seasonal evolution of the cold tongue, tropical convection and surface winds in the corrected version of the model, phase-lock of ENSO to the annual cycle looses its strong semi-annual component and becomes quite similar to the observed, although the amplitude of ENSO is reduced. Zonal wind stress response to the SST anomalies in the central-eastern Pacific also becomes more realistic. ENSO retrospective forecast experiments conducted with the directly coupled and the flux-corrected versions of the model demonstrate that deficiencies in the seasonal evolution of the cold tongue/Inter-Tropical Convergence Zone complex (that were largely due to the SEP bias in this model) and the related errors in the ENSO phase-lock to the annual cycle can seriously degrade ENSO prediction. By reducing these errors, ENSO predictive skill in the coupled model was substantially enhanced.     

The impact of atmospheric initialisation on seasonal prediction of tropical Pacific SST

The impact of realistic atmospheric initialisation on the seasonal prediction of tropical Pacific sea surface temperatures is explored with the Predictive Ocean–Atmosphere Model for Australia (POAMA) dynamical seasonal forecast system. Previous versions of POAMA used data from an Atmospheric Model Intercomparison Project (AMIP)-style simulation to initialise the atmosphere for the hindcast simulations. The initial conditions for the hindcasts did not, therefore, capture the true intra-seasonal atmospheric state. The most recent version of POAMA has a new Atmosphere and Land Initialisation scheme (ALI), which captures the observed intra-seasonal atmospheric state. We present the ALI scheme and then compare the forecast skill of two hindcast datasets, one with AMIP-type initialisation and one with realistic initial conditions from ALI, focussing on the prediction of El Niño. For eastern Pacific (Niño3) sea surface temperature anomalies (SSTAs), both experiments beat persistence and have useful SSTA prediction skill (anomaly correlations above 0.6) at all lead times (forecasts are 9 months duration). However, the experiment with realistic atmospheric initial conditions from ALI is an improvement over the AMIP-type initialisation experiment out to about 6 months lead time. The improvements in skill are related to improved initial atmospheric anomalies rather than an improved initial mean state (the forecast drift is worse in the ALI hindcast dataset). Since we are dealing with a coupled system, initial atmospheric errors (or differences between experiments) are amplified though coupled processes which can then lead to long lasting errors (or differences).     

A comparison of full-field and anomaly initialization for seasonal to decadal climate prediction

There are two main approaches for dealing with model biases in forecasts made with initialized climate models. In full-field initialization, model biases are removed during the assimilation process by constraining the model to be close to observations. Forecasts drift back towards the model’s preferred state, thereby re-establishing biases which are then removed with an a posterior lead-time dependent correction diagnosed from a set of historical tests (hindcasts). In anomaly initialization, the model is constrained by observed anomalies and deviates from its preferred climatology only by the observed variability. In theory, the forecasts do not drift, and biases may be removed based on the difference between observations and independent model simulations of a given period. Both approaches are currently in use, but their relative merits are unclear. Here we compare the skill of each approach in comprehensive decadal hindcasts starting each year from 1960 to 2009, made using the Met Office decadal prediction system. Both approaches are more skilful than climatology in most regions for temperature and some regions for precipitation. On seasonal timescales, full-field initialized hindcasts of regional temperature and precipitation are significantly more skilful on average than anomaly initialized hindcasts. Teleconnections associated with the El Niño Southern Oscillation are stronger with the full-field approach, providing a physical basis for the improved precipitation skill. Differences in skill on multi-year timescales are generally not significant. However, anomaly initialization provides a better estimate of forecast skill from a limited hindcast set.     

Flux correction: tests with a simple ocean-atmosphere model

 Flux correction schemes are used in order to suppress the drift of coupled ocean atmosphere models. This technique is tested for a simple box model of the climate system. Two “perfect” models of the ocean and the atmosphere are available. These are coupled to form an ocean-atmosphere model representing the true climate system. This climate system is simulated by a climate model which is also constructed by coupling those two perfect models. This time, however, both models are run first separately as models of the atmosphere and the ocean. In that case, “observations” from the climate system are prescribed at the ocean surface in the uncoupled models. It is assumed that these observations are imperfect. A drift results, when these models are coupled to form an ocean-atmosphere stimulation model. A flux adjustment scheme is implemented to remove this drift. It is argued that the merits and shortcomings of the flux correction technique can be assessed more clearly this way than by coupling imperfect models as is done normally. Sensitivity tests are performed where either radiation parameters are changed or a salt anomaly is implanted. Model parameters are chosen such that the ocean has a thermally direct circulation in the unperturbed climate state. It is found that the flux correction technique is performing satisfactorily as long as the imposed perturbations are small enough so that the ocean circulation does not change its sense. If, however, the model climate is close to the transition to an indirect circulation, then the flux correction technique is unreliable. The predictions of the coupled model with flux correction may deviate substantially from the response of the climate system in that case. Received: 4 December 1995/Accepted: 15 October 1996     

CoPIVEP: a theory-based analysis of coupled processes and interannual variability in the Equatorial Pacific in four coupled GCMs

 The interannual variability over the tropical Pacific and a possible link with the mean state or the seasonal cycle is examined in four coupled ocean-atmosphere general circulation models (GCM). Each model is composed of a high-resolution ocean GCM of either the tropical Pacific or near-global oceans coupled to a moderate-resolution atmospheric GCM, without using flux correction. The oceanic subsurface is considered to describe the mean state or the seasonal cycle through the analytical formulations of some potential coupled processes. These coupled processes characterise the zonal gradient of sea surface temperature (hereafter SST), the oceanic vertical gradient of temperature and the equatorial upwelling. The simulated SST patterns of the mean state and the interannual signals are generally too narrow. The grid of the oceanic model could control the structure of the SST interannual signals while the behaviour of the atmospheric model could be important in the link between the oceanic surface and the subsurface. The first SST EOFs are different between the coupled models, however, the second SST EOFs are quite similar and could correspond to the return to the normal state while that of the observations (COADS) could favour the initial anomaly. All the models seem to simulate a similar equatorial wave-like dynamics to return to the normal state. The more the basic state is unstable from the coupled processes point of view, the more the interannual signal are high. It seems that the basic state could control the intensity of the interannual variability. Two models, which have a significant seasonal variation of the interannual variance, also have a significant seasonal variation of the instability with a few months lag. The potential seasonal phase locking of the interannual fluctuations need to be examined in more models to confirm its existence in current tropical GCMs. Received: 30 July 1999 / Accepted: 25 April 2000     

Reduction of the thermocline feedback associated with mean SST bias in ENSO simulation

Associated with the double Inter-tropical convergence zone problem, a dipole SST bias pattern (cold in the equatorial central Pacific and warm in the southeast tropical Pacific) remains a common problem inherent in many contemporary coupled models. Based on a newly-developed coupled model, we performed a control run and two sensitivity runs, one is a coupled run with annual mean SST correction and the other is an ocean forced run. By comparison of these three runs, we demonstrated that a serious consequence of this SST bias is to severely suppress the thermocline feedback in a realistic simulation of the El Ni?o/Southern Oscillation. Firstly, the excessive cold tongue extension pushes the anomalous convection far westward from the equatorial central Pacific, prominently diminishing the convection-low level wind feedback and thus the air-sea coupling strength. Secondly, the equatorial surface wind anomaly exhibits a relatively uniform meridional structure with weak gradient, contributing to a weakened wind-thermocline feedback. Thirdly, the equatorial cold SST bias induces a weakened upper-ocean stratification and thus yields the underestimation of the thermocline-subsurface temperature feedback. Finally, the dipole SST bias underestimates the mean upwelling through (a) undermining equatorial mean easterly wind stress, and (b) enhancing convective mixing and thus reducing the upper ocean stratification, which weakens vertical shear of meridional currents and near-surface Ekman-divergence.     

The global distribution of sources of atmospheric decadal variability and mechanisms over the tropical Pacific and southern North America

 Understanding natural atmospheric decadal variability is an important element of climate research, and here we investigate the geographic and seasonal diversity in the balance between its competing sources. Data are provided by an ensemble of multi-decadal atmospheric general circulation model experiments, forced by observed sea surface temperatures (SSTs), and verified against observations. First, the nature of internal atmospheric variability is studied. By assessing its spectral character, we refute the idea that internal modes may persist or oscillate on multi-annual time-scales, either through mechanisms purely internal to the atmosphere, or via coupling to the land surface; instead, they behave as a white noise process. Second, and more importantly, the role of oceanic forcing, relative to internal variability, is investigated by extending the ‘analysis of variance’ technique to the frequency domain. Significance testing and confidence intervals are also discussed. In the tropics, atmospheric decadal variability is usually dominated by oceanic forcing, although for some regions less so than at interannual time-scales. A moderate oceanic impact is also found for some extratropical regions in some seasons. Verification against observed mean sea-level pressure (MSLP) data suggests that many of these influences are realistic, although some model errors are also revealed. In other mid- and high-latitude regions, local simulated decadal variability is dominated by random processes, i.e. the integrated effects of chaotic weather systems. Third, we focus on the mechanisms of decadal variability in two specific regions (where the model is well behaved). Over the tropical Pacific, the relative impact of SSTs on decadal MSLP is strongly seasonal such that it peaks in September to November (SON). This is explained by noting that the model atmosphere is responsive to SSTs a little farther west in SON than it is in other seasons, and here it picks up relatively more decadal power from the ocean (the western Pacific being less dominated by ENSO time-scales), causing atmospheric ‘signal-to-noise ratios’ to be enhanced at decadal timescales in SON. Over southern North America, a strong SST impact is found in summer and autumn, resulting in an upward trend of MSLP over recent decades. We suggest this is caused by decadal SST variability in the Caribbean (and to some extent the tropical northeast Pacific in summer), which induces anomalous convective heating over these regions and hence the wider MSLP response. Received: 30 November 1998 / Accepted: 22 April 1999     

The seasonal cycle in a coupled experiment with an atmospheric GCM and a two-layer equatorial ocean model

The mean state and the seasonal cycle in the tropical Pacific are studied, using a new coupled tropical ocean-global atmosphere model. The atmospheric component is a general circulation model and the oceanic component is a two and a half layer model of the tropical Pacific. The coupling is based on delocalized physics: the spatial resolution of the physics of the atmospheric component is the same as the spatial resolution of the oceanic model. No flux corrections are applied. A 31 year experiment has been made with the climatological observed sea surface temperature outside the area of coupling. We observe a quick drift of the model which, after three years, reaches a warm mean state. The temperature bias varies geographically between 1?^°C and 2?^°C, but, in spite of this default, the eastern part of the basin remains colder than the west. This contrast is shown to be dependent on the shoaling of the thermocline east of 160^°W. There is a significant seasonal cycle with an amplitude and phase of the seasonal variations which are well reproduced with respect to many other models. It is shown that interactions between the ocean and the atmosphere in the central and eastern Pacific are sufficient to explain the gross features of its evolution. In July, easterlies intensify in the Southern Hemisphere and lead to a strong upwelling and an enhanced evaporation in the eastern part of the basin. This induces a cooling throughout the area. The cooling reaches a first maximum in October in the easternmost part of the basin, then propagates westward along the equator with a decreasing amplitude. In January it is reinforced in the central part of the basin because of a divergence of the current, which is too strong. The mechanisms found here emphasize the role of the upwelling in maintaining the equatorial Pacific climate, and are in agreement with those found in other simplified coupled models.     

On the mechanisms in a tropical ocean–global atmosphere coupled general circulation model. Part I: mean state and the seasonal cycle

 The mechanisms responsible for the mean state and the seasonal and interannual variations of the coupled tropical Pacific-global atmosphere system are investigated by analyzing a thirty year simulation, where the LMD global atmospheric model and the LODYC tropical Pacific model are coupled using the delocalized physics method. No flux correction is needed over the tropical region. The coupled model reaches its regime state roughly after one year of integration in spite of the fact that the ocean is initialized from rest. Departures from the mean state are characterized by oscillations with dominant periodicites at annual, biennial and quadriennial time scales. In our model, equatorial sea surface temperature and wind stress fluctuations evolved in phase. In the Central Pacific during boreal autumn, the sea surface temperature is cold, the wind stress is strong, and the Inter Tropical Convergence Zone (ITCZ) is shifted northwards. The northward shift of the ITCZ enhances atmospheric and oceanic subsidence between the equator and the latitude of organized convention. In turn, the stronger oceanic subsidence reinforces equatorward convergence of water masses at the thermocline depth which, being not balanced by equatorial upwelling, deepens the equatorial thermocline. An equivalent view is that the deepening of the thermocline proceeds from the weakening of the meridional draining of near-surface equatorial waters. The inverse picture prevails during spring, when the equatorial sea surface temperatures are warm. Thus temperature anomalies tend to appear at the thermocline level, in phase opposition to the surface conditions. These subsurface temperature fluctuations propagate from the Central Pacific eastwards along the thermocline; when reaching the surface in the Eastern Pacific, they trigger the reversal of sea surface temperature anomalies. The whole oscillation is synchronized by the apparent meridional motion of the sun, through the seasonal oscillation of the ITCZ. This possible mechanism is partly supported by the observed seasonal reversal of vorticity between the equator and the ITCZ, and by observational evidence of eastward propagating subsurface temperature anomalies at the thermocline level. Received: 7 April 1997 / Accepted: 15 July 1998     

Mechanisms of tropical Pacific interannual-to-decadal variability in the ARPEGE/ORCA global coupled model

Spatial and temporal structures of interannual-to-decadal variability in the tropical Pacific Ocean are investigated using results from a global atmosphere–ocean coupled general circulation model. The model produces quite realistic mean state characteristics, despite a sea surface temperature cold bias and a thermocline that is shallower than observations in the western Pacific. The periodicity and spatial patterns of the modelled El Niño Southern Oscillations (ENSO) compare well with those observed over the last 100 years, although the quasi-biennial timescale is dominant. Lag-regression analysis between the mean zonal wind stress and the 20°C isotherm depth suggests that the recently proposed recharge-oscillator paradigm is operating in the model. Decadal thermocline variability is characterized by enhanced variance over the western tropical South Pacific (~7°S). The associated subsurface temperature variability is primarily due to adiabatic displacements of the thermocline as a whole, arising from Ekman pumping anomalies located in the central Pacific, south of the equator. Related wind anomalies appear to be caused by SST anomalies in the eastern equatorial Pacific. This quasi-decadal variability has a timescale between 8 years and 20 years. The relationship between this decadal tropical mode and the low-frequency modulation of ENSO variance is also discussed. Results question the commonly accepted hypothesis that the low-frequency modulation of ENSO is due to decadal changes of the mean state characteristics.     

Multi-year predictability in a coupled general circulation model

Multi-year to decadal variability in a 100-year integration of a BMRC coupled atmosphere-ocean general circulation model (CGCM) is examined. The fractional contribution made by the decadal component generally increases with depth and latitude away from surface waters in the equatorial Indo-Pacific Ocean. The relative importance of decadal variability is enhanced in off-equatorial “wings” in the subtropical eastern Pacific. The model and observations exhibit “ENSO-like” decadal patterns. Analytic results are derived, which show that the patterns can, in theory, occur in the absence of any predictability beyond ENSO time-scales. In practice, however, modification to this stochastic view is needed to account for robust differences between ENSO-like decadal patterns and their interannual counterparts. An analysis of variability in the CGCM, a wind-forced shallow water model, and a simple mixed layer model together with existing and new theoretical results are used to improve upon this stochastic paradigm and to provide a new theory for the origin of decadal ENSO-like patterns like the Interdecadal Pacific Oscillation and Pacific Decadal Oscillation. In this theory, ENSO-driven wind-stress variability forces internal equatorially-trapped Kelvin waves that propagate towards the eastern boundary. Kelvin waves can excite reflected internal westward propagating equatorially-trapped Rossby waves (RWs) and coastally-trapped waves (CTWs). CTWs have no impact on the off-equatorial sub-surface ocean outside the coastal wave guide, whereas the RWs do. If the frequency of the incident wave is too high, then only CTWs are excited. At lower frequencies, both CTWs and RWs can be excited. The lower the frequency, the greater the fraction of energy transmitted to RWs. This lowers the characteristic frequency (reddens the spectrum) of variability off the equator relative to its equatorial counterpart. At low frequencies, dissipation acts as an additional low pass filter that becomes more effective, as latitude increases. At the same time, ENSO-driven off-equatorial surface heating anomalies drive mixed layer temperature responses in both hemispheres. Both the eastern boundary interactions and the accumulation of surface heat fluxes by the surface mixed layer act to low pass filter the ENSO-forcing. The resulting off-equatorial variability is therefore more coherent with low pass filtered (decadal) ENSO indices [e.g. NINO3 sea-surface temperature (SST)] than with unfiltered ENSO indices. Consequently large correlations between variability and NINO3 extend further poleward on decadal time-scales than they do on interannual time-scales. This explains why decadal ENSO-like patterns have a broader meridional structure than their interannual counterparts. This difference in appearance can occur even if ENSO indices do not have any predictability beyond interannual time-scales. The wings around 15–20°S, and sub-surface variability at many other locations are predictable on interannual and multi-year time-scales. This includes westward propagating internal RWs within about 25° of the equator. The slowest of these take up to 4 years to reach the western boundary. This sub-surface predictability has significant oceanographic interest. However, it is linked to only low levels of SST variability. Consequently, extrapolation of delayed action oscillator theory to decadal time-scales might not be justified.     

Are atmospheric biases responsible for the tropical Atlantic SST biases in the CNRM-CM5 coupled model?

In this study, the CNRM-CM5 model is shown to simulate too warm SSTs in the tropical Atlantic as most state-of-the-art CMIP5 models. The warm bias develops within 1 or 2 months in decadal experiments initialised in January using an observationally derived state. To better quantify the role of the atmospheric biases in initiating this warm SST bias, several sensitivity experiments have been performed. In a first set of experiments, the surface solar net heat flux sent to the ocean model is academically corrected over the southeastern tropical Atlantic Ocean. This correction locally reduces the warm SST bias by more than 50 % with some remote impacts over equatorial regions. In contrast, the solar heat flux correction has locally little impact on the spring cooling. A second set of experiments quantifies the role of surface winds, using a nudging technique. When applied in a narrow equatorial region, the wind correction mainly improves the SST annual cycle amplitude along the Equator. It promotes not only the spring cooling along the Equator in preconditioning the mixed-layer depth but also in the southeastern Atlantic along the African coast. These local and remote effects are attributed to the more realistic representation of the oceanic equatorial circulation, driven by corrected winds. These results are consistent with those reported by Wahl et al. (Clim Dyn 36:891–906, 2011) in a very similar study with the Kiel Climate Model. The solar and wind biases have comparable effects in their study, although the importance of off-equatorial winds is less clear in our study. Diagnosing the wind energy flux provides a physical understanding of the equatorial region. When combining the corrections of both the equatorial wind and the southeastern solar heat flux, no obvious feedback between them is evidenced. The present study also emphasizes the need to consider two time-scales, the annual mean and the seasonal cycle, as well as two regions, the equatorial and the southeastern Atlantic regions, to comprehensively address the Atlantic SST bias. As pointed out in Richter (Clim Dyn, doi:10.1007/s00382-012-1624-5, 2013), the need to improve the atmospheric component of the CNRM-CM model is emphasized, even though strong positive coupling feedbacks are highlighted.     

A systematic approach to identify the sources of tropical SST errors in coupled models using the adjustment of initialised experiments

Understanding the sources of systematic errors in climate models is challenging because of coupled feedbacks and errors compensation. The developing seamless approach proposes that the identification and the correction of short term climate model errors have the potential to improve the modeled climate on longer time scales. In previous studies, initialised atmospheric simulations of a few days have been used to compare fast physics processes (convection, cloud processes) among models. The present study explores how initialised seasonal to decadal hindcasts (re-forecasts) relate transient week-to-month errors of the ocean and atmospheric components to the coupled model long-term pervasive SST errors. A protocol is designed to attribute the SST biases to the source processes. It includes five steps: (1) identify and describe biases in a coupled stabilized simulation, (2) determine the time scale of the advent of the bias and its propagation, (3) find the geographical origin of the bias, (4) evaluate the degree of coupling in the development of the bias, (5) find the field responsible for the bias. This strategy has been implemented with a set of experiments based on the initial adjustment of initialised simulations and exploring various degrees of coupling. In particular, hindcasts give the time scale of biases advent, regionally restored experiments show the geographical origin and ocean-only simulations isolate the field responsible for the bias and evaluate the degree of coupling in the bias development. This strategy is applied to four prominent SST biases of the IPSLCM5A-LR coupled model in the tropical Pacific, that are largely shared by other coupled models, including the Southeast Pacific warm bias and the equatorial cold tongue bias. Using the proposed protocol, we demonstrate that the East Pacific warm bias appears in a few months and is caused by a lack of upwelling due to too weak meridional coastal winds off Peru. The cold equatorial bias, which surprisingly takes 30 years to develop, is the result of an equatorward advection of midlatitude cold SST errors. Despite large development efforts, the current generation of coupled models shows only little improvement. The strategy proposed in this study is a further step to move from the current random ad hoc approach, to a bias-targeted, priority setting, systematic model development approach.     

A Coupled General Circulation Model for the Tropical Pacific Ocean and Global Atmosphere

On the basis of Zeng’s theoretical design, a coupled general circulation model (CGCM) is developed with its characteristics different from other CGCMs such as the unified vertical coordinates and subtraction of the standard stratification for both atmosphere and ocean, available energy consideration, and so on. The oceanic component is a free surface tropical Pacific Ocean GCM between 30oN and 30oS with horizontal grid spacing of 1o in latitude and 2o in longitude, and with 14 vertical layers. The atmospheric component it a global GCM with low-resolution of 4o in latitude and 5o in longitude, and two layers or equal man in the vertical between the surface and 200 hPa. The atmospheric GCM includes comprehensive physical processes. The coupled model is subjected to seasonally-varying cycle. Several coupling experiments, ranging from straight forward coupling without flux correction to one with flux correction, and to so-called predictor-corrector monthly coupling (PCMC), are conducted to show the existence and final controlling of the climate drift in the coupled system. After removing the climate drift with the PCMC scheme, the coupled model is integrated for more than twenty years. The results show reasonable simulations of the annual mean and its seasonal cycle of the atmospheric and oceanic circulation. The model also produces the coherent interannual variations of the climate system, manifesting the observed El Ni?o / Southern Oscillation (ENSO).     

北半球高纬地区年际尺度循环过程中的气-海-冰相互作用关系

基于一个全球气-海-冰耦合模式数值模拟结果,对北半球高纬度地区年际尺度的气-海-冰相互作用进行了分析。在所使用的全球气-海-冰耦合模式中,大气环流模式和陆面过程模式来自国家气候中心,海洋环流模式和海冰模式来自中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室。采用一种逐日通量距平耦合方案实现次网格尺度海冰非均匀条件下大气环流模式和海洋环流模式在高纬地区的耦合。只对50 a模拟结果中的后30 a结果进行了分析。在分析中,首先对滤波后的北半球高纬度地区海平面气压、表面大气温度、海表面温度、海冰密集度及海表面感热通量的标准化距平做联合复经验正交函数分解,取第一模进行重建,然后讨论了在一个循环周期(约4 a)中北半球高纬度地区气-海-冰的作用关系。结果表明:(1)当北大西洋涛动处于正位相时,格陵兰海出现南风异常,使表面大气温度升高,海洋失去感热通量减少,海洋表面温度升高,海冰密集度减小;当北大西洋涛动处于负位相时,格陵兰海出现北风异常,使表面大气温度降低,海洋失去感热通量增多,海洋表面温度降低,海冰密集度增加。巴伦支海变化特点与格陵兰海相似,但在时间上并不完全一致。(2)多年平均而言,北冰洋内部靠近极点区域为冷中心。当北冰洋内部为低压异常时,因异常中心偏向太平洋一侧,使北冰洋内部靠近太平洋部分为暖平流异常,靠近大西洋一侧为冷平流异常。伴随着暖、冷平流异常,这两侧分别出现暖异常和冷异常,海表面给大气的感热通量分别偏少和偏多,上述海区海表面温度分别偏高和偏低,海冰密集度分别偏小和偏大。当北冰洋内部为高压异常时特点正好与上述相反。由上述分析结果可知,在海洋、大气年际循环中,大尺度大气环流变率起主导作用,海洋表面温度和海冰密集度变化主要是对大气环流变化的响应。     

Impacts of Indian Ocean SST biases on the Indian Monsoon: as simulated in a global coupled model

In this study, the impact of the ocean–atmosphere coupling on the atmospheric mean state over the Indian Ocean and the Indian Summer Monsoon (ISM) is examined in the framework of the SINTEX-F2 coupled model through forced and coupled control simulations and several sensitivity coupled experiments. During boreal winter and spring, most of the Indian Ocean biases are common in forced and coupled simulations, suggesting that the errors originate from the atmospheric model, especially a dry islands bias in the Maritime Continent. During boreal summer, the air-sea coupling decreases the ISM rainfall over South India and the monsoon strength to realistic amplitude, but at the expense of important degradations of the rainfall and Sea Surface Temperature (SST) mean states in the Indian Ocean. Strong SST biases of opposite sign are observed over the western (WIO) and eastern (EIO) tropical Indian Ocean. Rainfall amounts over the ocean (land) are systematically higher (lower) in the northern hemisphere and the south equatorial Indian Ocean rainfall band is missing in the control coupled simulation. During boreal fall, positive dipole-like errors emerge in the mean state of the coupled model, with warm and wet (cold and dry) biases in the WIO (EIO), suggesting again a significant impact of the SST errors. The exact contributions and the distinct roles of these SST errors in the seasonal mean atmospheric state of the coupled model have been further assessed with two sensitivity coupled experiments, in which the SST biases are replaced by observed climatology either in the WIO (warm bias) or EIO (cold bias). The correction of the WIO warm bias leads to a global decrease of rainfall in the monsoon region, which confirms that the WIO is an important source of moisture for the ISM. On the other hand, the correction of the EIO cold bias leads to a global improvement of precipitation and circulation mean state during summer and fall. Nevertheless, all these improvements due to SST corrections seem drastically limited by the atmosphere intrinsic biases, including prominently the unimodal oceanic position of the ITCZ (Inter Tropical Convergence Zone) during summer and the enhanced westward wind stress along the equator during fall.     

热带太平洋年代际平均气候态变化与ENSO循环

文中用观测的热带太平洋海表温度资料、风应力资料和OLR资料，通过多时间尺度分析，将与ENSO有关的变化分为3个主要的分量，一是2～7a的ENSO循环尺度，二是8～20a的年代际尺度，三是20a以上的平均气候态变化。讨论了热带太平洋这种平均气候态变化的主要特征以及与ENSO循环的关系，并用耦合模式的数值试验来研究平均气候态的变化对ENSO循环的影响。结果表明热带太平洋的平均气候态在20世纪70年代后期发生了一次由冷态向暖态的变化，主要增暖区是沿赤道以及热带东太平洋的，海表温度变化最大中心可以达到0.6℃。伴随着海表温度的变化，赤道西太平洋的西风距平加强，赤道东太平洋的东风距平也加强，在赤道中太平洋形成了一个加强的辐合中心。年代际平均气候冷暖态的变化对ENSO最直接的线性影响是使ElNio位相增加，而形成ENSO冷位相和暖位相的不对称。另一方面较暖的平均气候态可能引起海洋和大气之间的耦合加强，导致ENSO循环振荡有所加强。     

印度洋对ENSO事件的响应:观测与模拟

观测事实显示,在El Ni(n～)o期间,伴随着赤道中东太平洋表层海温(SST)的升高,热带印度洋SST出现正距平.作者利用海气耦合模式模拟了印度洋对ENSO事件的上述响应,并进而讨论了其物理机制.所用模式为法国国家科研中心Pierre-Simon-Laplace 全球环境科学联合实验室(IPSL)发展的全球海气耦合模式.该模式成功地控制了气候漂移,能够合理再现印度洋的基本气候态.观测中与ENSO相关的热带印度洋SST变化,表现为全海盆一致的正距平,并且这种变化要滞后赤道中东太平洋SST变化大约一个季度,意味着它主要是对东太平洋SST强迫的一种遥响应,模式结果也支持这一机制,尽管模式中的南方涛动现象被夸大了,使得模拟的与ENSO相关联的SST正距平的位置南移,阿拉伯海和孟加拉湾被负距平(而不是正距平)所控制.研究表明,东太平洋主要通过大气桥影响潜热释放来影响印度洋SST变化.赤道东太平洋El Ni(n～)o事件的发展,导致印度洋上空风场异常自东而西传播;伴随着风场的变化,潜热发生相应变化,并最终导致SST异常的发生.非洲东海岸受索马里急流控制的海域,其SST的变化不能简单地利用热通量的变化来解释.证据显示,印度洋的增暖是ENSO事件发生的结果而不是其前期信号.