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.     

一个海洋-海冰耦合模式的北极气候模拟

This paper evaluates the simulation of Arctic sea ice states using an ocean-ice coupled model that employs LASG/IAP(the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/the Institute of Atmospheric Physics) Climate Ocean Model(LICOM) and the sea-ice model from the Bergen Climate Model(BCM).It is shown that the coupled model can reasonably reproduce the major characteristics of the mean state,annual cycle,and interannual variability of the Arctic sea ice concentration.The coupled model also shows biases that were generally presented in other models,such as the underestimation of summer sea ice concentration and thickness as well as the unsatisfactory sea ice velocity.Sensitivity experiments indicate that the insufficient performance of the ocean model at high latitudes may be the main reason for the biases in the coupled model.The smoother and the fake "island",which had to be used due to the model’s grid in the North Pole region,likely caused the ocean model’s weak performance.Sea ice model thermodynamics are also responsible for the sea ice simulation biases.Therefore,both the thermodynamic module of the sea ice component and the model grid of the ocean component need to be further improved.     

Numerical Simulation of Atmosphere-Ocean-Sea Ice Interaction During Interannual Cycle in High Northern Latitudes

The interannual atmosphere-ocean-sea ice interaction （AOSI） in high northern latitudes is studied with a global atmosphere-ocean-sea ice coupled model system, in which the model components of atmosphere and land surface are from China National Climate Center and that of ocean and sea ice are from LASG, Institute of Atmospheric Physics, Chinese Academy of Sciences. A daily flux anomaly correction scheme is employed to couple the atmosphere model and the ocean model with the effect of inhomogenity of sea ice in high latitudes is considered. The coupled model system has been run for 50 yr and the results of the last 30 years are analyzed. After the sea level pressure （SLP）, surface air temperature （SAT）, sea surface temperature （SST）, sea ice concentration （SIC）, and sea surface sensible heat flux （SHF） are filtered with a digital filter firstly, their normalized anomalies are used to perform the decomposition of combined complex empirical orthogonal function （CCEOF） and then they are reconstructed with the leading mode. The atmosphere-ocean-sea ice interactions in high northern latitudes during a periodical cycle （approximately 4 yr） are analyzed. It is shown that： （1） When the North Atlantic Oscillation （NAO） is in its positive phase, the southerly anomaly appears in the Greenland Sea, SAT increases, the sea loses less SHF, SST increases and SIC decreases accordingly; when the NAO is in its negative phase, the northerly anomaly appears in the Greenland Sea, SAT decreases, the sea loses more SHF, SST decreases and SIC increases accordingly. There are similar features in the Barents Sea, but the phase of evolution in the Barents Sea is different from that in the Greenland Sea. （2） For an average of multi-years, there is a cold center in the inner part of the Arctic Ocean near the North Pole. When there is an anomaly of low pressure, which is closer to the Pacific Ocean, in the inner part of the Arctic Ocean, anomalies of warm advection appear in the region near the Pacif     

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.     

An Overview of BCC Climate System Model Development and Application for Climate Change Studies

This paper reviews recent progress in the development of the Beijing Climate Center Climate System Model(BCC-CSM) and its four component models(atmosphere,land surface,ocean,and sea ice).Two recent versions are described:BCC-CSM1.1 with coarse resolution(approximately 2.8125°×2.8125°) and BCC-CSM1.1(m) with moderate resolution(approximately 1.125°×1.125°).Both versions are fully coupled climate-carbon cycle models that simulate the global terrestrial and oceanic carbon cycles and include dynamic vegetation.Both models well simulate the concentration and temporal evolution of atmospheric CO_2 during the 20th century with anthropogenic CO2 emissions prescribed.Simulations using these two versions of the BCC-CSM model have been contributed to the Coupled Model Intercomparison Project phase five(CMIP5) in support of the Intergovernmental Panel on Climate Change(IPCC) Fifth Assessment Report(AR5).These simulations are available for use by both national and international communities for investigating global climate change and for future climate projections.Simulations of the 20th century climate using BCC-CSMl.l and BCC-CSMl.l(m) are presented and validated,with particular focus on the spatial pattern and seasonal evolution of precipitation and surface air temperature on global and continental scales.Simulations of climate during the last millennium and projections of climate change during the next century are also presented and discussed.Both BCC-CSMl.l and BCC-CSMl.l(m) perform well when compared with other CMIP5 models.Preliminary analyses indicate that the higher resolution in BCC-CSM1.1(m) improves the simulation of mean climate relative to BCC-CSMl.l,particularly on regional scales.     

How Are El Nio and La Nia Events Improved in an Eddy-Resolving Ocean General Circulation Model?

The present study compares the performance of two versions of the LASG/IAP(State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics) Climate System Ocean Model(LICOM) in reproducing the interannual variability associated with El Nio and La Nia events in the tropical Pacific. Both versions are forced with the identical boundary conditions from observed or reanalysis data, in which one version has a finer spatial resolution of(1/10)° in the horizontal domain and 55 vertical layers, and the other version has a coarse resolution of 1° in the horizontal domain and 30 vertical layers. ENSO simulations form the two versions are compared with observations and, in particular, the improvements with regard to ENSO by the finer resolution ocean model are emphasized. As a result of the finer spatial resolution, both the vertical temperature gradient and vertical velocity are better represented in the equatorial Pacific than they are by the coarse resolution model; and thus, the corresponding vertical advections of temperature are more reasonable. Besides the mean climatology, simulated ENSO events and relevant feedbacks are much improved in the finer resolution model. A heat budget analysis suggests that both thermocline feedback and Ekman feedback are mainly responsible for the rapid increase in temperature anomalies during the developing and mature phases of ENSO events.     

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).     

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.     

Impact of Surface Exchange Coefficients and Sea Surface Cooling on Tropical Cyclone Simulation

The surface flux exchange associated with the exchange coefficients and upper ocean conditions is essential to the development of tropical cyclones (TCs). Using the Weather Research and Forecasting (WRF) model, the present study has investigated the impact of exchange coefficients and ocean coupling during Super Typhoon Saomai (2006). Firstly, two experiments with different formula of roughness are conducted. The experiment with the Donelan formula for drag coefficient (Cd) and ramped formula for enthalpy coefficient (Ck) can simulate stronger intensity compared to other experiments due to the increased surface wind and enthalpy fluxes. That is because the new formulas allows for a smaller Cd and larger Ck in the high wind regime than the former formulas did. Moreover, two coupled simulations between WRF and a one-dimensional ocean model are conducted to examine the feedback of sea surface cooling to the TC. In the experiments with a horizontal uniform mixed layer depth of 70 m, the sea surface cooling is too weak to change the evolution of TC. While in the experiment with an input mixed layer calculated using the Hybrid Coordinate Ocean Model (HYCOM) data, the significant sea surface cooling induces obvious impact on TC intensity and structure. Under the negative feedback of sea surface cooling, the sensible and latent heat fluxes decreases, especially in the right part of Saomai (2006). The negative feedback with coupled ocean model plays a vital role in simulating the intensity and structure of TC.     

Simulated Spatiotemporal Response of Ocean Heat Transport to Freshwater Enhancement in North Atlantic and Associated Mechanisms

The Atlantic Meridional Overturning Circulation(AMOC)transports a large amount of heat to northern high latitudes,playing an important role in the global climate change.Investigation of the freshwater perturbation in North Atlantic(NA)has become one of the hot topics in the recent years.In this study,the mechanism and pathway of meridional ocean heat transport(OHT)under the enhanced freshwater input to the northern high latitudes in the Atlantic are investigated by an ocean-sea ice-atmosphere coupled model.The results show that the anomalous OHT in the freshwater experiment(FW)is dominated by the meridional circulation kinetic and ocean thermal processes.In the FW,OHT drops down during the period of weakened AMOC while the upper tropical ocean turns warmer due to the retained NA warm currents.Conversely,OHT recovers as the AMOC recovers,and the mechanism can be generalized as:1)increased ocean heat content in the tropical Southern Ocean during the early integration provides the thermal condition for the recovery of OHT in NA;2)the OHT from the Southern Ocean enters the NA through the equator alongthe deep Ekman layer;3)in NA,the recovery of OHT appears mainly along the isopycnic layers of 24.70-25.77 below the mixing layer.It is then transported into the mixing layer from the "outcropping points"innorthern high latitudes,and finally released to the atmosphere by the ocean-atmosphere heat exchange.     

Performance of Convective Parameterization Schemes in Asia Using RegCM:Simulations in Three Typical Regions for the Period 1998–2002

This study discusses the sensitivity of convective parameterization schemes(CPSs) in the Regional Climate Model(version 4.3)(Reg CM4.3) over East/South Asia. The simulations using different CPSs in Reg CM are compared to discover a suitable scheme for this region, as the performance of different schemes is greatly influenced by region and seasonality. Over Southeast China and the Bay of Bengal, the Grell scheme exhibits the lowest RMSEs of summer precipitation compared to observed data. Moreover, the Emanuel over land and Grell over ocean(ELGO) scheme enhances the simulation, in comparison with any single CPS(Grell/Emanuel) over Western Ghats, Sri Lanka, and Southeast India. Over the Huang–Huai–Hai Plain(3H) and Tibetan Plateau(TP) regions of China, the Tiedtke scheme simulates the more reasonable summer precipitation with high correlation coefficient and comparable amplitude. Especially, it reproduces a minimum convective precipitation bias of 8 mm d-1and the lowest RMSEs throughout the year over East/South Asia. Furthermore, for seasonal variation of precipitation, the Tiedtke scheme results are closer to the observed data over the 3H and TP regions. However, none of the CPSs is able to simulate the seasonal variation over North Pakistan(NP). In comparison with previous research, the results of this study support the Grell scheme over South Asia. However, the Tiedtke scheme shows superiority for the 3H, TP and NP regions. The thicker PBL, less surface latent heat flux, the unique ability of deep convection and the entrainment process in the Tiedtke scheme are responsible for reducing the wet bias.     

Evaluation of Haney-Type Surface ThermalBoundary Conditions Using a CoupledAtmosphere and Ocean Model

A coupled atmosphere-ocean model developed at the Institute for Space Studies at NASA Goddard Space Flight Center (Russell et al., 1995) was used to verify the validity of Haney-type surface thermal boundary condition, which linearly connects net downward surface heat flux Q to air / sea temperature difference △T by a relaxation coefficient k. The model was initiated from the National Centers for Environmental Prediction (NCEP) atmospheric observations for 1 December 1977, and from the National Ocean Data Center (NODC) global climatological mean December temperature and salinity fields at 1° ×1° resolution. The time step is 7.5 minutes. We integrated the model for 450 days and obtained a complete model-generated global data set of daily mean downward net surface flux Q, surface air temperature TA,and sea surface temperature To. Then, we calculated the cross-correlation coefficients (CCC) between Q and △T. The ensemble mean CCC fields show (a) no correlation between Q and △T in the equatorial regions, and (b) evident correlation (CCC≥ 0.7) between Q and △T in the middle and high latitudes.Additionally, we did the variance analysis and found that when k= 120 W m-2K-1, the two standard deviations, σQ and σk△T, are quite close in the middle and high latitudes. These results agree quite well with a previous research (Chu et al., 1998) on analyzing the NCEP re-analyzed surface data, except that a smaller value of k (80 W m-2K-1) was found in the previous study.     

SEASONAL CHANGES OF SOIL MOISTURE INDUCED BY INCREASED CO2 AS SIMULATED BY THE ATMOSPHERIC GENERAL CIRCULATION AND MIXED LAYER OCEAN MODEL

The atmospheric general circulation model coupled to the mixed layer ocean model has been used to simulate thechanges of the global soil moisture.Comparing the simulated results with observations,it is shown that the model is ca-pable of doing sensitive experiments about the carbon dioxide change.The 2×CO_2/1×CO_2 comparison shows that there are the obvious changes of the soil moisture in the global forfour seasons.There are the wet soil moisture in the lower latitudes of both hemispheres and dry soil moisture in the mid-dle latitudes of both hemispheres for four seasons.The dry soil moisture in summer and wet in other seasons are foundin the northern higher latitudes.The analyses of the physical feedbacks responsible for the CO_2-induced changes of soil moisture show that the bud-gets of the surface water and heat are the important factors.     

Dependence of the agcm climatology on the method of prescribing surface boundary conditions and its climatological implication

By using IAP 9L AGCM, two sets of long-term climatological integration have been per-formed with the two different interpolation procedures for generating the daily surface boundary conditions. One interpolation procedure is the so-called “traditional” scheme, for which the daily surface boundary conditions are obtained by linearly interpolating between the observed monthly mean values, however the observed monthly means cannot be preserved after interpolation. The other one is the “new” scheme, for which the daily surface boundary conditions are obtained by linearly interpolating between the "artificial" monthly mean values which are based on, but are dif-ferent from the observed ones, after interpolating with this new scheme, not only the observed monthly mean values are preserved, the time series of the new generated daily values is also more consistent with the observation. Comparison of the model results shows that the differences of the globally or zonally averaged fields between these two integrations are quite small, and this is due to the compensating effect between the different regions. However, the differences of the two patterns (the global or regional geographical distributions), are quite significant, for example, the magni-tude of the difference in the JJA mean rainfall between these two integrations can exceed 2 mm/day over Asian monsoon regions, and the difference in DJF mean surface air temperature can also exceed 2oC over this region. The fact that the model climatology depends quite strongly on the method of prescribing the daily surface boundary conditions suggests that in order to validate the climate model or to predict the short-term climate anomalies, either the " new* interpolation scheme or the high frequency surface boundary conditions (e.g., daily or weekly data instead of the monthly data) should be introduced. Meanwhile, as for the coupled model, the daily coupling scheme between the different component cli?mate models (e.g., atmospheric and oceanic general circulation models) is preferred in order to partly eliminate the “climate drift” problem which may appear during the course of direct coupling.     

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.     

Improved ENSO forecasts by assimilating sea surface temperature observations into an intermediate coupled model

A simple method for initializing intermediate coupled models (ICMs) using only sea surface temperature (SST) anomaly data is comprehensively tested in two sets of hindcasts with a new ICM. In the initialization scheme, both the magnitude of the nudging parameter and the duration of the assimilation are considered, and initial conditions for both atmosphere and ocean are generated by running the coupled model with SST anomalies nudged to the observations. A comparison with the observations indicates that the scheme can generate realistic thermal fields and surface dynamic fields in the equatorial Pacific through hindcast experiments. An ideal experiment is performed to get the optimal nudging parameters which include the nudging intensity and nudging time length. Twelve-month-long hindcast experiments are performed with the model over the period 1984–2003 and the period 1997–2003. Compared with the original prediction results, the model prediction skills are significantly improved by the nudging method especially beyond a 6-month lead time during the two different periods. Potential problems and further improvements are discussed regarding the new coupled assimilation system.     

Development of a Single-Column Model in RegCM4 and Its Preliminary Application for Evaluating PBL Schemes in Simulating the Dry Convection Boundary Layer

A single-column model (SCM) is developed in the regional climate model RegCM4. The evolution of a dry convection boundary layer (DCBL) is used to evaluate this SCM. Moreover, four planetary boundary layer (PBL) schemes, namely the Holtslag-Boville scheme (HB), Yonsei University scheme (YSU), and two University of Washington schemes (UW01, Grenier-Bretherton-McCaa scheme and UW09, Bretherton-Park scheme), are compared by using the SCM approach. A large-eddy simulation (LES) of the DCBL is performed as a benchmark to examine how well a PBL parameterization scheme reproduces the LES results, and several diagnostic outputs are compared to evaluate the schemes. The results show that the SCM is proper constructed. In general, with the DCBL case, the YSU scheme performs best for reproducing the LES results, which include well-mixed features and vertical sensible heat fluxes; the simulated wind speed, turbulent kinetic energy, entrainment flux, and height of the entrainment zone are all underestimated in the UW09; the UW01 has all those biases of the UW09 but larger, and the simulated potential temperature is not well mixed; the HB is the least skillful scheme, by which the PBL height, entrainment flux, height of the entrainment zone, and the vertical gradients within the mixed layer are all overestimated, and a inversion layer near the top of the surface layer is wrongly simulated.Although more cases and further testing are required, these simulations show encouraging results towards the use of this SCM framework for evaluating the simulated physical processes by the RegCM4.     

Impacts of Seasonal Fossil and Ocean Emissions on the Seasonal Cycle of Atmospheric CO2

The seasonal cycle of atmospheric CO2 at surface observation stations in the northern hemisphere is driven primarily by net ecosystem production (NEP) fluxes from terrestrial ecosystems. In addition to NEP from terrestrial ecosystems, surface fluxes from fossil fuel combustion and ocean exchange also contribute to the seasonal cycle of atmospheric CO2. Here the authors use the Goddard Earth Observing System-Chemistry (GEOS-Chem) model (version 8-02-01), with modifications, to assess the impact of these fluxes on the seasonal cycle of atmospheric CO2 in 2005. Modifications include monthly fossil and ocean emission inventories. CO2 simulations with monthly varying and annual emission inventories were carried out separately. The sources and sinks of monthly averaged net surface flux are different from those of annual emission inventories for every month. Results indicate that changes in monthly averaged net surface flux have a greater impact on the average concentration of atmospheric CO2 in the northern hemisphere than on the average concentration for latitudes 30-90°S in July. The concentration values differ little between both emission inventories over the latitudinal range from the equator to 30°S in January and July. The accumulated impacts of the monthly averaged fossil and ocean emissions contribute to an increase of the total global monthly average of CO2 from May to December.An apparent discrepancy for global average CO2 concentration between model results and observation was because the observation stations were not sufficiently representative. More accurate values for monthly varying net surface flux will be necessary in future to run the CO2 simulation.     

CAS-FGOALS Datasets for the Two Interglacial Epochs of the Holocene and the Last Interglacial in PMIP4

Two versions of the Chinese Academy of Sciences Flexible Global Ocean–Atmosphere–Land System model(CASFGOALS), version f3-L and g3, are used to simulate the two interglacial epochs of the mid-Holocene and the Last Interglacial in phase 4 of the Paleoclimate Modelling Intercomparison Project(PMIP4), which aims to study the impact of changes in orbital parameters on the Earth’s climate. Following the PMIP4 experimental protocols, four simulations for the mid-Holocene and two simulations for the La...     

A Hybrid Coupled Model for the Pacific Ocean–Atmosphere System.Part I: Description and Basic Performance

A hybrid coupled model(HCM) is constructed for El Nino–Southern Oscillation(ENSO)-related modeling studies over almost the entire Pacific basin. An ocean general circulation model is coupled to a statistical atmospheric model for interannual wind stress anomalies to represent their dominant coupling with sea surface temperatures. In addition, various relevant forcing and feedback processes exist in the region and can affect ENSO in a significant way; their effects are simply represented using historical data and are incorporated into the HCM, including stochastic forcing of atmospheric winds, and feedbacks associated with freshwater flux, ocean biology-induced heating(OBH), and tropical instability waves(TIWs). In addition to its computational efficiency, the advantages of making use of such an HCM enable these related forcing and feedback processes to be represented individually or collectively, allowing their modulating effects on ENSO to be examined in a clean and clear way. In this paper, examples are given to illustrate the ability of the HCM to depict the mean ocean state, the circulation pathways connecting the subtropics and tropics in the western Pacific, and interannual variability associated with ENSO. As satellite data are taken to parameterize processes that are not explicitly represented in the HCM, this work also demonstrates an innovative method of using remotely sensed data for climate modeling. Further model applications related with ENSO modulations by extratropical influences and by various forcings and feedbacks will be presented in Part II of this study.