The Effects of Model Resolution on the Simulation of Regional Climate Extreme Events

The fifth-generation Pennsylvania State University/NCAR Mesoscale Model Version 3 (MM5V3) was used to simulate extreme heavy rainfall events over the Yangtze River Basin in June 1999. The effects of model's horizontal and vertical resolution on the extreme climate events were investigated in detail. In principle, the model was able to characterize the spatial distribution of monthly heavy precipitation. The results indicated that the increase in horizontal resolution could reduce the bias of the modeled heavy rain and reasonably simulate the change of daily precipitation during the study period. A finer vertical resolution led to obviously improve rainfall simulations with smaller biases, and hence, better resolve heavy rainfall events. The increase in both horizontal and vertical resolution could produce better predictions of heavy rainfall events. Not only the rainfall simulation altered in the cases of different horizontal and vertical grid spacing, but also other meteorological fields demonstrated diverse variations in terms of resolution change in the model. An evident improvement in the simulated sea level pressure resulted from the increase of horizontal resolution, but the simulation was insensitive to vertical grid spacing. The increase in vertical resolution could enhance the simulation of surface temperature as well as atmospheric circulation at low levels, while the simulation of circulation at middle and upper levels were found to be much less dependent on changing resolution. In addition, cumulus parameterization schemes showed high sensitivity to horizontal resolution. Different convective schemes exhibited large discrepancies in rainfall simulations with regards to changing resolution. The percentage of convective precipitation in the Grell scheme increased with increasing horizontal resolution. In contrast, the Kain-Fritsch scheme caused a reduced ratio of convective precipitation to total rainfall accumulations corresponding to increasing horizontal resolution.     

Modeling the tropical Pacific Ocean using a regional coupled climate model

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

What Determines the Amplitude of ENSO Events?

In this paper, the dynamic effect of oceanic upwelling on the intensity of El Nio-Southern Oscillation (ENSO) is studied using a simple coupled model (Zebiak-Cane Model). The term balance analysis in the temperature variability equation shows that the anomalous upwelling of the mean vertical temperature gradient and the mean advection of the anomalous meridional temperature gradient are the two of most important factors that determine the intensity of ENSO events, in which the "vertical oceanic heat flux" in the eastern equatorial Pacific (EEP) is the primary influencing factor. The lag correlation between "vertical heat flux (VHF)" and ENSO intensity shows that the highest correlation occurs when the former leads the latter by one to two weeks. The VHF is positively correlated with the background thermocline strength in the EEP, and an increase of both could result in strong ENSO variability. Comparison of the forced and coupled experiments suggests that the coupled process can affect both the intensity and frequency of ENSO.     

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

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

Nonlinearity Modulating Intensities and Spatial Structures of Central Pacific and Eastern Pacific El Nio Events

This paper compares data from linearized and nonlinear Zebiak–Cane model, as constrained by observed sea surface temperature anomaly(SSTA), in simulating central Pacific(CP) and eastern Pacific(EP) El Nio. The difference between the temperature advections(determined by subtracting those of the linearized model from those of the nonlinear model),referred to here as the nonlinearly induced temperature advection change(NTA), is analyzed. The results demonstrate that the NTA records warming in the central equatorial Pacific during CP El Nio and makes fewer contributions to the structural distinctions of the CP El Nio, whereas it records warming in the eastern equatorial Pacific during EP El Nio, and thus significantly promotes EP El Nio during El Nio–type selection. The NTA for CP and EP El Nio varies in its amplitude,and is smaller in CP El Nio than it is in EP El Nio. These results demonstrate that CP El Nio are weakly modulated by small intensities of NTA, and may be controlled by weak nonlinearity; whereas, EP El Nio are significantly enhanced by large amplitudes of NTA, and are therefore likely to be modulated by relatively strong nonlinearity. These data could explain why CP El Nio are weaker than EP El Nio. Because the NTA for CP and EP El Nio differs in spatial structures and intensities, as well as their roles within different El Nio modes, the diversity of El Nio may be closely related to changes in the nonlinear characteristics of the tropical Pacific.     

A New Approach for Classifying Two Types of El Nio Events

正In recent decades, the typical El Nio events with the warmest SSTs in the tropical eastern Pacific have become less common, and a different of El Nio with the warmest SSTs in the central Pacific, which is flanked on the east and west by cooler SSTs, has become more fre-quent. The more recent type of El Nio was referred to as central Pacific El Nio, warm pool El Nio, or dateline El Nio, or the El Nio Modoki. Central Pacific El Nio links to a different tropical-to-extratropical teleconnection and exerts different impacts on climate, and several classification approaches have been proposed. In this study, a new classification approach is proposed, which is based on the linear combination (sum or difference) of the two leading Empirical Orthogonal Functions (EOFs) of tropical Pacific Ocean sea surface temperature anomaly (SSTA), and the typical El Ni o index (TENI) and the central El Nio index (CENI) are able to be derived by projecting the observed SSTA onto these combinations. This classification not only reflects the characteristics of non-orthogonality between the two types of events but also yields one period peaking at approximate two to seven years. In particular, this classification can distin-guish the different impacts of the two types of events on rainfall in the following summer in East China. The typical El Nio events tend to induce intensified rainfall in the Yangtze River valley, whereas the central Pacific El Nio tends to induce intensified rainfall in the Huaihe River valley. Thus, the present approach may be appropriate for studying the impact of different types of El Nio on the East Asian climate.     

Influence of positive and negative Indian Ocean Dipoles on ENSO via the Indonesian Throughflow: Results from sensitivity experiments

The role of the Indonesian Throughflow(ITF) in the influence of the Indian Ocean Dipole(IOD) on ENSO is investigated using version 2 of the Parallel Ocean Program(POP2) ocean general circulation model. We demonstrate the results through sensitivity experiments on both positive and negative IOD events from observations and coupled general circulation model simulations. By shutting down the atmospheric bridge while maintaining the tropical oceanic channel, the IOD forcing is shown to influence the ENSO event in the following year, and the role of the ITF is emphasized. During positive IOD events,negative sea surface height anomalies(SSHAs) occur in the eastern Indian Ocean, indicating the existence of upwelling.These upwelling anomalies pass through the Indonesian seas and enter the western tropical Pacific, resulting in cold anomalies there. These cold temperature anomalies further propagate to the eastern equatorial Pacific, and ultimately induce a La Nia-like mode in the following year. In contrast, during negative IOD events, positive SSHAs are established in the eastern Indian Ocean, leading to downwelling anomalies that can also propagate into the subsurface of the western Pacific Ocean and travel further eastward. These downwelling anomalies induce negative ITF transport anomalies, and an El Nio-like mode in the tropical eastern Pacific Ocean that persists into the following year. The effects of negative and positive IOD events on ENSO via the ITF are symmetric. Finally, we also estimate the contribution of IOD forcing in explaining the Pacific variability associated with ENSO via ITF.     

Sensitivity Difference in the Extratropical Atmosphere to Two Types of El Nio Events

A comparison of sensitivity in extratropical circulation in the Northern Hemisphere(NH)and Southern Hemisphere(SH)is conducted through observational analyses and diagnostic linear model experiments for two types of El Nio events,the traditional El Nio with the strongest warmth in the eastern tropical Pacific(EP El Nio)and the El Nio Modoki with the strongest warmth in the central tropical Pacific(CP El Nio).It is shown that CP El Nio favors the occurrence of a negative-phase Northern Annular Mode(NAM),while EP El Nio favors that of the Pacific-North American(PNA)pattern.In SH,both EP and CP El Nio induce a negative phase Southern Annular Mode(SAM).However,the former has a greater amplitude,which is consistent with the stronger sea surface temperature(SST)warmth.The difference in the two types of El Nio events in NH may originate from the dependence of heating-induced extratropical response on the location of initial heating,which may be associated with activity of the stationary wave.In SH,the lack of sensitivity to the location of heating can be associated with weaker activity of the stationary wave therein.     

THE DESIGN OF A LIMITED AREA MESOSCALE MODEL AND PREDICTION EXPERIMENTS

Based on the dynamic frame of a 6-layer primitive equations model, a mesoscale primitive equations model is designed that includes larger range of physical processes. It is run in (x, y, σ) system with the model atmosphere topping at 10 hPa and availability of multiple conditions of horizontal boundaries, both horizontal and vertical resolution are adjustable; the precipitation scheme includes large scale and deep cumulus convective precipitation; the ground temperature is computed using surface heat budget equations; exchanges between land atmosphere and between ocean and atmosphere are considered, the Liouis format is used in the computation of vertical exchange budget; a scheme that combines the second and fourth order is employed in horizontal diffusion in which the coefficient is the function of the location of grid points and wind fields; the integration scheme is in the form of economic central difference.With the resolution that horizontal grids are spaced at intervals of 80 km and vertical length is unequally spaced into 16 layers, the model is experimented with 26 cases of forecast. The result has shown stable model computation, good prediction of major synoptic patterns and close reproduction of real precipitation. Statistics for a number of assessment indexes are given in this paper and comparisons are made to the original 6-layer model in respect of the forecasting ability and model properties.     

Improved Tropical Cyclone Forecasts with Increased Vertical Resolution in the TRAMS Model

The numerical simulation of typhoons has been found to be very sensitive to the vertical resolution of the model. During the updating of the TRAMS model from version 1.0 to 3.0, the horizontal resolution has been increased from 36 km to 9 km, while the vertical layer number only increased from 55 to 65 layers. The lack of high vertical resolution limits the performance of the TRAMS model in typhoon forecasting to a certain extent. In order to study the potential improvement of typhoon forecasting by increasing the vertical resolution, this paper increases the vertical resolution of the TRAMS model from 65 to 125 layers for the first time for a comparative simulation test. The results of the case study with Typhoon Hato (2017) show that the model with high vertical resolution can significantly enhance the warm structure caused by water vapor flux convergence and vertical transport, thus accurately simulating the rapid strengthening process of the typhoon. Meanwhile, the model with 125-layer vertical resolution can simulate the asymmetric structural characteristics of the wind field, which are closer to the observations and can help to reduce the bias in typhoon track forecasting. The improvement of vertical resolution is also trialed by using the batch test results of several landfalling typhoons in 2016-2017. The experimental results show that the typhoon forecast of the model becomes consistent with the observations only when the number of vertical layers of the model increases to about 125 layers, which in turn causes a large computational burden. In the next step, we will try to solve the computational burden problem caused by ultra-high vertical resolution with the top boundary nesting technique, and realize the application of high vertical resolution in the actual operation of the TRAMS model.     

Impacts of a GCM＇s Resolution on MJO Simulation

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

A Hybrid Coupled Ocean–Atmosphere Model and Its Simulation of ENSO and Atmospheric Responses

A new hybrid coupled model(HCM) is presented in this study, which consists of an intermediate tropical Pacific Ocean model and a global atmospheric general circulation model. The ocean component is the intermediate ocean model(IOM)of the intermediate coupled model(ICM) used at the Institute of Oceanology, Chinese Academy of Sciences(IOCAS). The atmospheric component is ECHAM5, the fifth version of the Max Planck Institute for Meteorology atmospheric general circulation model. The HCM integrates its atmospheric and oceanic components by using an anomaly coupling strategy. A100-year simulation has been made with the HCM and its simulation skills are evaluated, including the interannual variability of SST over the tropical Pacific and the ENSO-related responses of the global atmosphere. The model shows irregular occurrence of ENSO events with a spectral range between two and five years. The amplitude and lifetime of ENSO events and the annual phase-locking of SST anomalies are also reproduced realistically. Despite the slightly stronger variance of SST anomalies over the central Pacific than observed in the HCM, the patterns of atmospheric anomalies related to ENSO,such as sea level pressure, temperature and precipitation, are in broad agreement with observations. Therefore, this model can not only simulate the ENSO variability, but also reproduce the global atmospheric variability associated with ENSO, thereby providing a useful modeling tool for ENSO studies. Further model applications of ENSO modulations by ocean–atmosphere processes, and of ENSO-related climate prediction, are also discussed.     

Contrasting the Skills and Biases of Deterministic Predictions for the Two Types of El Nio

The tropical Pacific has begun to experience a new type of El Nio, which has occurred particularly frequently during the last decade, referred to as the central Pacific(CP) El Nio. Various coupled models with different degrees of complexity have been used to make real-time El Nio predictions, but high uncertainty still exists in their forecasts. It remains unknown as to how much of this uncertainty is specifically related to the new CP-type El Nio and how much is common to both this type and the conventional Eastern Pacific(EP)-type El Nio. In this study, the deterministic performance of an El Nio–Southern Oscillation(ENSO) ensemble prediction system is examined for the two types of El Nio. Ensemble hindcasts are run for the nine EP El Nio events and twelve CP El Nio events that have occurred since 1950. The results show that(1) the skill scores for the EP events are significantly better than those for the CP events, at all lead times;(2) the systematic forecast biases come mostly from the prediction of the CP events; and(3) the systematic error is characterized by an overly warm eastern Pacific during the spring season, indicating a stronger spring prediction barrier for the CP El Nio. Further improvements to coupled atmosphere–ocean models in terms of CP El Nio prediction should be recognized as a key and high-priority task for the climate prediction community.     

High resolution global modeling of the atmospheric circulation

An informal review is presented of recent developments in numerical simulation of the global atmospheric circulation with very fine numerical resolution models. The focus is on results obtained recently with versions of the GFDL SKYHI model and the Atmospheric Model for the Earth Simulator (AFES) global atmospheric models. These models have been run with effective horizontal grid resolution of 10–40 km and fine vertical resolution. The results presented demonstrate the utility of such models for the study of a diverse range of phenomena. Specifically the models are shown to simulate the development of tropical cyclones with peak winds and minimum central pressures comparable to those of the most intense hurricanes actually observed. More fundamentally, the spectrum of energy content in the mesoscale in the flow can be reproduced by these models down to near the smallest explicitly-resolved horizontal scales. In the middle atmosphere it is shown that increasing horizontal resolution can lead to significantly improved overall simulation of the global-scale circulation. The application of the models to two specific problems requiring very fine resolution global will be discussed. The spatial and temporal variability of the vertical eddy flux of zonal momentum associated with gravity waves near the tropopause is evaluated in the very fine resolution AFES model. This is a subject of great importance for understanding and modelling the flow in the middle atmosphere. Then the simulation of the small scale variations of the semidiurnal surface pressure oscillation is analyzed, and the signature of significant topographic modulation of the semidiurnal atmospheric tide is identified.     

Modeling Arctic Boundary Layer Cloud Streets at Grey-zone Resolutions

To better understand how model resolution affects the formation of Arctic boundary layer clouds,we investigated the influence of grid spacing on simulating cloud streets that occurred near Utqiaġvik(formerly Barrow),Alaska,on 2 May 2013 and were observed by MODIS(the Moderate Resolution Imaging Spectroradiometer).The Weather Research and Forecasting model was used to simulate the clouds using nested domains with increasingly fine resolution ranging from a horizontal grid spacing of 27 km in the boundary-layer-parameterized mesoscale domain to a grid spacing of 0.111 km in the large-eddy-permitting domain.We investigated the model-simulated mesoscale environment,horizontal and vertical cloud structures,boundary layer stability,and cloud properties,all of which were subsequently used to interpret the observed roll-cloud case.Increasing model resolution led to a transition from a more buoyant boundary layer to a more shear-driven turbulent boundary layer.The clouds were stratiform-like in the mesoscale domain,but as the model resolution increased,roll-like structures,aligned along the wind field,appeared with ever smaller wavelengths.A stronger vertical water vapor gradient occurred above the cloud layers with decreasing grid spacing.With fixed model grid spacing at 0.333 km,changing the model configuration from a boundary layer parameterization to a large-eddy-permitting scheme produced a more shear-driven and less unstable environment,a stronger vertical water vapor gradient below the cloud layers,and the wavelengths of the rolls decreased slightly.In this study,only the large-eddy-permitting simulation with gird spacing of 0.111 km was sufficient to model the observed roll clouds.     

Major modes of short-term climate variability in the newly developed NUIST Earth System Model (NESM)

A coupled earth system model(ESM) has been developed at the Nanjing University of Information Science and Technology(NUIST) by using version 5.3 of the European Centre Hamburg Model(ECHAM), version 3.4 of the Nucleus for European Modelling of the Ocean(NEMO), and version 4.1 of the Los Alamos sea ice model(CICE). The model is referred to as NUIST ESM1(NESM1). Comprehensive and quantitative metrics are used to assess the model's major modes of climate variability most relevant to subseasonal-to-interannual climate prediction. The model's assessment is placed in a multi-model framework. The model yields a realistic annual mean and annual cycle of equatorial SST, and a reasonably realistic precipitation climatology, but has difficulty in capturing the spring–fall asymmetry and monsoon precipitation domains. The ENSO mode is reproduced well with respect to its spatial structure, power spectrum, phase locking to the annual cycle, and spatial structures of the central Pacific(CP)-ENSO and eastern Pacific(EP)-ENSO; however, the equatorial SST variability,biennial component of ENSO, and the amplitude of CP-ENSO are overestimated. The model captures realistic intraseasonal variability patterns, the vertical-zonal structures of the first two leading predictable modes of Madden–Julian Oscillation(MJO), and its eastward propagation; but the simulated MJO speed is significantly slower than observed. Compared with the T42 version, the high resolution version(T159) demonstrates improved simulation with respect to the climatology, interannual variance, monsoon–ENSO lead–lag correlation, spatial structures of the leading mode of the Asian–Australian monsoon rainfall variability, and the eastward propagation of the MJO.     

Numerical Study for Characteristic Change of Asian Summer Monsoon Circulation and Its Influence Mechanism during the EI Nino Period

In this paper, the relation between Asian summer monsoon circulation and sea surface temperature anomalies over equatorial central-eastern Pacific is investigated by using a global spectral model. This model has nine layers in the vertical and the model variables are represented in the horizontal as truncated expansions of the surface spherical harmonics with rhomboidal truncation at wave number 15. The model involves comparatively complete physical processes and parameterizations with mountains.Using the above model, two experimental schemes are designed, namely control case and anomalous sea surface temperature case. The above two schemes are respectively integrated for forty days and the simulated results are obtained from the last 30-day averaged simulations.The simulations show that positive SST anomalies over equatorial central-eastern Pacific weakens Indian monsoon circulation,decreases precipitation in Indian sub-continent whereas it intensifies East Asian monsoon circulation and increases preci     

Distinguished Effects of Interannual Salinity Variability on the Development of the Central-Pacific El Ni o Events

El Nio events in the central equatorial Pacific (CP) are gaining increased attention,due to their increasing intensity within the global warming context.Various physical processes have been identified in the climate system that can be responsible for the modulation of El Nio,especially the effects of interannual salinity variability.In this work,a comprehensive data analysis is performed to illustrate the effects of interannual salinity variability using surface and subsurface salinity fields from the Met Office ENSEMBLES (EN3) quality controlled ocean dataset.It is demonstrated that during the developing phase of an El Nio event,a negative sea surface salinity (SSS) anomaly in the western-central basin acts to freshen the mixed layer (ML),decrease oceanic density in the upper ocean,and stabilize the upper layers.These related oceanic processes tend to reduce the vertical mixing and entrainment of subsurface water at the base of the ML,which further enhances the warm sea surface temperature (SST) anomalies associated with the El Nio event.However,the effects of interannually variable salinity are much more significant during the CP-El Nio than during the eastern Pacific (EP) El Nio,indicating that the salinity effect might be an important contributor to the development of CP-El Nio events.     

Connections Between Different Types of El Nio and Southern/Northern Oscillation

It has long been acknowledged that there are two types of El Nio events, i.e., the eastern Pacific El Nio (EE) and the central Pacific El Nio (CE), according to the initial position of the anomalous warm water and its propagation direction. In this paper, the oceanic and atmospheric evolutions and the possible mechanisms of the two types of El Nio events were examined. It is found that all the El Nio events, CE or EE, could be attributed to the joint impacts of the eastward advection of warm water from the western Pacific warm pool (WPWP) and the local warming in the equatorial eastern Pacific. Before the occurrence of CE events, WPWP had long been in a state of being anomalous warm, so the strength of eastward advection of warm water was much stronger than that of EE, which played a major role in the formation of CE. While for the EE events, most contribution came from the local warming of the equatorial eastern Pacific. It is further identified that the immediate cause leading to the difference of the two types of El Nio events was the asynchronous variations of the Southern Oscillation (SO) and the Northern Oscillation (NO) as defined by Chen in 1984. When the transition from the positive phase of the NO (NO+) to NO- was prior to that from SO+ to SO-, there would be eastward propagation of westerly anomalies from the tropical western Pacific induced by NO and hence the growth of warm sea surface temperature anomalies in WPWP and its eastward propagation. This was followed by lagged SO-induced weakening of southeast trade winds and local warming in the equatorial eastern Pacific. These were conducive to the occurrence of the CE. On the contrary, the transition from SO+ to SO- leading the transition of NO would favor the occurrence of EE type events.     

Comparison of Constant and Time-variant Optimal Forcing Approaches in El Nio Simulations by Using the Zebiak–Cane Model

Model errors offset by constant and time-variant optimal forcing vector approaches(termed COF and OFV, respectively)are analyzed within the framework of El Nio simulations. Applying the COF and OFV approaches to the well-known Zebiak–Cane model, we re-simulate the 1997 and 2004 El Nio events, both of which were poorly degraded by a certain amount of model error when the initial anomalies were generated by coupling the observed wind forcing to an ocean component. It is found that the Zebiak–Cane model with the COF approach roughly reproduced the 1997 El Nio, but the 2004 El Nio simulated by this approach defied an ENSO classification, i.e., it was hardly distinguishable as CP-El Nio or EP-El Nio. In both El Nio simulations, substituting the COF with the OFV improved the fit between the simulations and observations because the OFV better manages the time-variant errors in the model. Furthermore, the OFV approach effectively corrected the modeled El Nio events even when the observational data(and hence the computational time) were reduced.Such a cost-effective offset of model errors suggests a role for the OFV approach in complicated CGCMs.