Fundamental Framework and Experiments of the ThirdGeneration of IAP/LASG World OceanGeneral Circulation Model

A new generation of the IAP / LASG world ocean general circulation model is designed and presented based on the previous 20-layer model, with enhanced spatial resolutions and improved parameterizations. The model uses a triangular-truncated spectral horizontal grid system with its zonal wave number of 63 (T63) to match its atmospheric counterpart of a T63 spectral atmosphere general circulation model in a planned coupled ocean-atmosphere system. There are 30 layers in vertical direction, of which 20 layers are located above 1000 m for better depicting the permanent thermocline. As previous ocean models developed in IAP / LASG, a free surface (rather than “rigid-lid” approximation) is included in this model. Compared with the 20-layer model, some more detailed physical parameterizations are considered, including the along / cross isopycnal mixing scheme adapted from the Gent-MacWilliams scheme. The model is spun up from a motionless state. Initial conditions for temperature and salinity are taken from the three-dimensional distributions of Levitus’ annual mean observation. A preliminary analysis of the first 1000-year integration of a control experiment shows some encouraging improvements compared with the twenty-layer model, particularly in the simulations of permanent thermocline, thermohaline circulation, meridional heat transport, etc. resulted mainly from using the isopycnal mixing scheme. However, the use of isopycnal mixing scheme does not significantly improve the simulated equatorial thermocline. A series of numerical experiments show that the most important contribution to the improvement of equatorial thermocline and the associated equatorial under current comes from reducing horizontal viscosity in the equatorial regions. It is found that reducing the horizontal viscosity in the equatorial Atlantic Ocean may slightly weaken the overturning rate of North Atlantic Deep Water.

A Nine-layer Atmospheric General Circulation Model and Its Performance

Different versions of a new nine-layer general circulation model which is rhomboidally truncated at zonal wavenumber 15(L9R15)are introduced in this paper.On using the observed global monthly sea surfaCe temperature(SST)and sea ice(SI)data from 1979 to 1988 offered by the internahonal Atmospheric Model Iute-comparison Program(AMIP),these different model versions were integrated for the ten-year AMIP period. Results show that the model iscapable of simulating the basic states of the atmosphere and its interannual variability,and in performing reasonablesensitivity experiments.

An Eddy-Permitting Oceanic General Circulation Model and Its Preliminary Evaluation

An eddy-permitting, quasi-global oceanic general circulation model, LICOM (LASG/IAP (State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics) Climate System Ocean Model), with a uniform grid of 0.5? × 0.5? is established. Forced by wind stresses from Hellerman and Rosenstain (1983), a 40-yr integration is conducted with sea surface temperature and salinity being restored to the Levitus 94 datasets. The evaluation of the annual mean climatology of the LICOM control run shows that the large-scale circulation can be well reproduced. A comparison between the LICOM control run and a parallel integration of L30T63, which has the same framework but a coarse resolution, is also made to con?rm the impact of resolution on the model performance. On account of the reduction of horizontal viscosity with the enhancement of the horizontal resolution, LICOM improves the simulation with respect to not only the intensity of the large scale circulations, but also the magnitude and structure of the Equatorial Undercurrent and South Equatorial Current. Taking advantage of the ?ne grid size, the pathway of the Indonesian Through?ow (ITF) is better represented in LICOM than in L30T63. The transport of ITF in LICOM is more convergent in the upper layer. As a consequence, the Indian Ocean tends to get warmer in LICOM. The poleward heat transports for both the global and individual basins are also signi?cantly improved in LICOM. A decomposed analysis indicates that the transport due to the barotropic gyre, which primarily stands for the barotropic e?ect of the western boundary currents, plays a crucial role in making the di?erence.

Dynamical Framework of IAP Nine-Level Atmospheric General Circulation Model

The dynamical framework of the nine-level version of the IAP AGCM is presented in this paper. The emphasis of the model's description is put on the following two aspects:(1) A model's standard atmosphere, which is a satisfactory approximation to the observed troposphere and lower stratosphere standard atmosphere, is introduced into the equations of the model to permit a more accurate calculation of the vertical transport terms, especially near the tropopause; (2) The vertical levels of the model are carefully selected to guarantee a smooth dependence of layer thickness upon pressure in order to reduce the truncation error involved in the unequal interval vertical finite-differencing. For testing the model, two kinds of linear baroclinic Rossby-Haurwitz waves, one of which has a dynamically stable vertical structure and the other has a relatively unstable one, are constructed to provide initial conditions for numerical experiments. The two waves have been integrated for more than 300 days and 100 days respectively by using the model and both of them are propagating westward with almost identical phase-speed during the time period of the integrations. No obvious change of the wave patterns is found at the levels in the model's troposphere. The amplitudes of both two waves at the uppermost level, however, exhibit rather significant oscillation with time, of which the periods are exactly 20 days and 25 days espectively.The explanation of this interesting phenomena is still under investigation.

A Numerical World Ocean General Circulation Model

This paper describes a numerical model of the world ocean based on the fully primitive equations. A 'Standard' ocean state is introduced into the equations of the model and the perturbed thermodynamic variables are used in the modlc's calculations. Both a free upper surface and a bottom topography are included in the model and a sigma coordinate is used to normalize the model's vertical component. The model has four unevenly-spaced layers and 4 × 5 horizontal resolution based on C-grid system. The finite-difference scheme of the model is designed to conserve the gross available energy in order to avoid fictitious energy generation or decay.The model has been tested in response to the annual mean surface wind stress, sea level air pressure and sea level air temperature as a preliminary step to its further improvement and its coupling with a global atmospheric general circulation model. Some of results, including currents, temperature and sea surface elevation simulated by the mode! arc presented.

The Dynamic Plateau Monsoon Index and Its Association with General Circulation Anomalies

Based on monthly ECMWF reanalysis-Interim （ERA-Interim） reanalysis data, along with monthly precipitation and temperature data, the Dynamic Plateau Monsoon Index （DPMI） is defined. The results of a contrast analysis of the DPMI versus the Traditional Plateau Monsoon Index （TPMI） are described. The response of general circulation to northern Qinghai-Xizang Plateau summer monsoon anomalies and the correlation of the DPMI with general circulation anomalies are investigated. The results show that, the DPMI reflected meteorological elements better and depicted climate variation more accurately than the TPMI. In years when the plateau summer monsoon is strong, the low over the plateau and the trough near the eastern coast of Asia are deeper and higher than normal over South China. This correlation corresponds to two anomalous cyclones over the plateau and the eastern coast of Asia and an anomalous anticyclone in South China. The plateau and its adjacent regions are affected by anomalous southwesterly winds that transport more moisture to South China and cause more precipitation. The lower reaches of the Yangtze River appear to receive more precipitation by means of the strong westerly water vapor flow transported from the ＂large triangle affecting the region＂. In years when the plateau summer monsoon is weak, these are opposite. The plateau monsoon is closely related to the intensity and position of the South Asian high, and the existence of a teleconnection pattern in the mid-upper levels suggests a possible linkage of the East Asian monsoon and the Indian monsoon to the plateau summer monsoon.

A Study on Sulfate Optical Properties and Direct Radiative Forcing Using LASG-IAP General Circulation Model

The direct radiative forcing （DRF） of sulfate aerosols depends highly on the atmospheric sulfate loading and the meteorology, both of which undergo strong regional and seasonal variations. Because the optical properties of sulfate aerosols are also sensitive to atmospheric relative humidity, in this study we first examine the scheme for optical properties that considers hydroscopic growth. Next, we investigate the seasonal and regional distributions of sulfate DRF using the sulfate loading simulated from NCAR CAM-Chem together with the meteorology modeled from a spectral atmospheric general circulation model （AGCM） developed by LASG-IAP. The global annual-mean sulfate loading of 3.44 mg m-2 is calculated to yield the DRF of -1.03 and -0.57 W m-2 for clear-sky and all-sky conditions, respectively. However, much larger values occur on regional bases. For example, the maximum all-sky sulfate DRF over Europe, East Asia, and North America can be up to -4.0 W m-2. The strongest all-sky sulfate DRF occurs in the Northern Hemispheric July, with a hemispheric average of -1.26 W m-2. The study results also indicate that the regional DRF are strongly affected by cloud and relative humidity, which vary considerably among the regions during different seasons. This certainly raises the issue that the biases in model-sinmlated regional meteorology can introduce biases into the sulfate DRF. Hence, the model processes associated with atmospheric humidity and cloud physics should be modified in great depth to improve the simulations of the LASG-IAP AGCM and to reduce the uncertainty of sulfate direct effects on global and regional climate in these simulations.

Primary Reasoning behind the Double ITCZ Phenomenon in a Coupled Ocean-Atmosphere General Circulation Model

This paper investigates the processes behind the double ITCZ phenomenon, a common problem in Coupled ocean-atmosphere General Circulation Models (CGCMs), using a CGCM—FGCM-0 (Flexible General Circulation Model, version 0). The double ITCZ mode develops rapidly during the ?rst two years of the integration and becomes a perennial phenomenon afterwards in the model. By way of Singular Value Decomposition (SVD) for SST, sea surface pressure, and sea surface wind, some air-sea interactions are analyzed. These interactions prompt the anomalous signals that appear at the beginning of the coupling to develop rapidly. There are two possible reasons, proved by sensitivity experiments: (1) the overestimated east-west gradient of SST in the equatorial Paci?c in the ocean spin-up process, and (2) the underestimated amount of low-level stratus over the Peruvian coast in CCM3 (the Community Climate Model, Version Three). The overestimated east-west gradient of SST brings the anomalous equatorial easterly. The anomalous easterly, a?ected by the Coriolis force in the Southern Hemisphere, turns into an anomalous westerly in a broad area south of the equator and is enhanced by atmospheric anomalous circulation due to the underestimated amount of low-level stratus over the Peruvian coast simulated by CCM3. The anomalous westerly leads to anomalous warm advection that makes the SST warm in the southeast Paci?c. The double ITCZ phenomenon in the CGCM is a result of a series of nonlocal and nonlinear adjustment processes in the coupled system, which can be traced to the uncoupled models, oceanic component, and atmospheric component. The zonal gradient of the equatorial SST is too large in the ocean component and the amount of low-level stratus over the Peruvian coast is too low in the atmosphere component.

A parameterization scheme of vertical mixing due to inertial internal wave breaking in the ocean general circulation model

Based on the theoretical spectral model of inertial internal wave breaking （fine structure） proposed previ- ously, in which the effects of the horizontal Coriolis frequency component f-tilde on a potential isopycnal are taken into account, a parameterization scheme of vertical mixing in the stably stratified interior be- low the surface mixed layer in the ocean general circulation model （OGCM） is put forward preliminarily in this paper. Besides turbulence, the impact of sub-mesoscale oceanic processes （including inertial internal wave breaking product） on oceanic interior mixing is emphasized. We suggest that adding the inertial inter- hal wave breaking mixing scheme （F-scheme for short） put forward in this paper to the turbulence mixing scheme of Canuto et al. （T-scheme for short） in the OGCM, except the region from 15°S to 15°N. The numeri- cal results ofF-scheme by usingWOA09 data and an OGCM （LICOM, LASG/IAP climate system ocean model） over the global ocean are given. A notable improvement in the simulation of salinity and temperature over the global ocean is attained by using T-scheme adding F-scheme, especially in the mid- and high-latitude regions in the simulation of the intermediate water and deep water. We conjecture that the inertial internal wave breaking mixing and inertial forcing of wind might be one of important mechanisms maintaining the ventilation process. The modeling strength of the Atlantic meridional overturning circulation （AMOC） by using T-scheme adding F-scheme may be more reasonable than that by using T-scheme alone, though the physical processes need to be further studied, and the overflow parameterization needs to be incorporated. A shortcoming in F-scheme is that in this paper the error of simulated salinity and temperature by using T-scheme adding F-scheme is larger than that by using T-scheme alone in the subsurface layer.

Influences of Freshwater from Major Rivers on Global Ocean Circulation and Temperatures in the MIT Ocean General Circulation Model

Responses of global ocean circulation and temperature to freshwater runoff from major rivers were studied by blocking regional runoff in the global ocean general circulation model （OGCM） developed at the Massachusetts Institute of Technology. Runoff into the tropical Atlantic, the western North Pacific, and the Bay of Bengal and northern Arabian Sea were selectively blocked. The blocking of river runoff first resulted in a salinity increase near the river mouths （2 practical salinity units）. The saltier and, therefore, denser water was then transported to higher latitudes in the North Atlantic, North Pacific, and southern Indian Ocean by the mean currents. The subsequent density contrasts between northern and southern hemispheric oceans resulted in changes in major ocean currents. These anomalous ocean currents lead to significant temperature changes （I^C -2~C） by the resulting anomalous heat transports. The current and temperature anomalies created by the blocked river runoff propagated from one ocean basin to others via coastal and equatorial Kelvin waves. This study suggests that river runoff may be playing an important role in oceanic salinity, temperature, and circulations; and that partially or fully blocking major rivers to divert freshwater for societal purposes might significantly change ocean salinity, circulations, temperature, and atmospheric climate. Further studies are necessary to assess the role of river runoff in the coupled atmosphere-ocean system.

THE EFFECTS OF GENERAL CIRCULATION ANOMALIES ON PRECIPITATION IN THE RAINING SEASON OF SHANDONG

Using the monthly mean reanalysis data of height, wind and vertical velocity of NCEP/NCAR from 1958 to 2003, the characteristics of 500hPa height anomalies in the Northern Hemisphere, convergence and divergence on the levels of 200hPa and 925hPa and vertical stream over the Shandong province were studied and compared with that over North China. The study shows that the teleconnection pattern called East Asia-Pacific pattern (EAP) plays important roles in the summer rainfall of Shandong. Weaker Okhotsk sea high, stronger upward stream over tropical areas and Shandong and northward location of subtropical high and stronger Indian low, are likely to lead to more rain in Shandong province in summer, but there are some flood (drought) years in which these characteristics such as the location and intensity of subtropical high change with the least extent. The difference between Shandong and North China in 500hPa height is that there is a teleconnection pattern called Europe / Asia pattern in North China, while there is a negative East Asia Pattern in Shandong during drought years. There is a teleconnection relation between Shandong summer rainfall and convective activities happening in tropical areas in spring and summer accomplished by the Rossby waves. During the flood years in Shandong, the convective activities happening in tropical India Ocean and tropical west Pacific and Walker circulation are much stronger than that during the drought years, but the convective activities happening in tropical east Pacific are weaker than those during the drought years.

The African Climate as Predicted by the IAP Grid-Point Nine-Layer Atmospheric General Circulation Model (IAP-9L-AGCM)

ＴｈｅＡｆｒｉｃａｎＣｌｉｍａｔｅａｓＰｒｅｄｉｃｔｅｄｂｙｔｈｅＩＡＰＧｒｉｄ－ＰｏｉｎｔＮｉｎｅ－ＬａｙｅｒＡｔｍｏｓｐｈｅｒｉｃＧｅｎｅｒａｌＣｉｒｃｕｌａｔｉｏｎＭｏｄｅｌ（ＩＡＰ－９Ｌ－ＡＧＣＭ）ＣｈｉｎｅｋｅＴｈｅｏＣｈｉｄｉｅｚｉｅ①，...

EXPERIMENTS OF A REDUCED GRID IN LASG/IAP WORLD OCEAN GENERAL CIRCULATION MODELS (OGCMs)

Due to the decrease in grid size associated with the convergence of meridians toward the poles inspherical coordinates, the time steps in many global climate models with finite-difference method are restrictedto be unpleasantly small. To overcome the problem, a reduced grid is introduced to LASG/IAP world oceangeneral circulation models. The reduced grid is implemented successfully in the coarser resolutions versionmodel L30T63 at first. Then, it is carried out in the improved version model LICOM with finer resolutions. Inthe experiment with model L30T63, under time step unchanged though, execution time per single model run isshortened significantly owing to the decrease of grid number and filtering execution in high latitudes. Resultsfrom additional experiments with L30T63 show that the time step of integration can be quadrupled at most inreduced grid with refinement ratio 3. In the experiment with model LICOM and with the model’s original timestep unchanged, the model covered area is extended to the whole globe from its original case with the grid pointof North Pole considered as an isolated island and the results of experiment are shown to be acceptable.

Evolving Perspectives on Abrupt Seasonal Changes of the General Circulation

Professor Duzheng YE（Tu-cheng YEH） was decades ahead of his time in proposing a model experiment to investigate whether abrupt seasonal changes of the general circulation can arise through circulation feedbacks alone, unrelated to underlying inhomogeneities at the lower boundary. Here, we introduce Professor YEH＇s ideas during the 1950 s and 1960 s on the general circulation and summarize the results and suggestions of Yeh et al.（1959） on abrupt seasonal changes. We then review recent advances in understanding abrupt seasonal changes arising from model experiments like those proposed by Yeh et al.（1959）. The model experiments show that circulation feedbacks can indeed give rise to abrupt seasonal transitions.In these transitions, large-scale eddies that originate in midlatitudes and interact with the zonal mean flow and meridional overturning circulations in the tropics play central roles.

Simulating Tropical Instability Waves in the Equatorial Eastern Pacific with a Coupled General Circulation Model

Satellite observations of SSTs have revealed the existence of unstable waves in the equatorial eastern Pacific and Atlantic oceans. These waves have a 20-40-day periodicity with westward phase speeds of 0.4-0.6 m s^-1 and wavelengths of 1000-2000 km during boreal summer and fall. They are generally called tropical instability waves （TIWs）. This study investigates TIWs simulated by a high-resolution coupled atmosphere-ocean general circulation model （AOGCM）. The horizontal resolution of the model is 120 km in the atmosphere, and 30 km longitude by 20 km latitude in the ocean. Model simulations show good agreement with the observed main features associated with TIWs. The results of energetics analysis reveal that barotropic energy conversion is responsible for providing the main energy source for TIWs by extracting energy from the meridional shear of the climatological-mean equatorial currents in the mixed layer. This deeper and northward-extended wave activity appears to gain its energy through baroclinic conversion via buoyancy work, which further contributes to the asymmetric distribution of TIWs. It is estimated that the strong cooling effect induced by equatorial upwelling is partially （-30%-40%） offset by the equatorward heat flux due to TIWs in the eastern tropical Pacific during the seasons when TIWs are active. The atmospheric mixed layer just above the sea surface responds to the waves with enhanced or reduced vertical mixing. Furthermore, the changes in turbulent mixing feed back to sea surface evaporation, favoring the westward propagation of TIWs. The atmosphere to the south of the Equator also responds to TIWs in a similar way, although TIWs are much weaker south of the Equator.

PREVIOUS CONVECTION ANOMALY IN AUTUMN AND WINTER ASSOCIATED WITH THE GENERAL CIRCULATION OVER EAST ASIA IN WINTER AND SPRING AND APRIL PRECIPITATION IN SHANDONG

Using OLR and 850 hPa and 200 hPa wind fields data (1979 - 2006), this paper diagnoses the characteristics of convection over the tropical area in preceding autumns and winters in association with April precipitation anomalies in Shandong province. It is found that preceding convection anomalies over the Western Pacific Warm Pool in December have close relationships with the April precipitation in Shandong. Further analysis of the relationship with the general circulation over the East Asia shows that the convection anomaly over the Western Pacific Warm Pool has close relationships with the Main East Asian Trough, the Hadley cell over East Asia and the Walker cell. The characteristics of East Asian atmospheric circulation anomalies accompanied with stronger (weaker) convection are consistent with those of less (more) April precipitation anomalies in Shandong. Therefore, the convection anomaly over the tropics in December may be an important indicator for April precipitation in Shandong.

A three-dimensional ocean general circulation model for mesoscale eddies——Ⅱ Diagnostic analysis

A three-dimensional density field associated with mesoscaie unstable waves generated by the 3-D, primitive-equation model (Wang and Ikeda, 1996) is provided to the quasi-geostrophic pressure tendency and ω-equations, and to the (ageostrophic) Q-vector equation. Diagnostic analyses, analogous to the approaches in meteorology: ω-equation and Q-vector method, are for the first time developed to examine the mesoscaie dynamical processes and mechanisms of the unstable waves propagating in the mid-latitude ocean. The weaknesses and strengths of these two diagnostic approaches are evaluated and compared to the model results. The Q-vector method is then recommended to diagnose the vertical motion associated with the mesoscaie dynamics from a hydrographic CTD (conductivity-temperature-depth) array, while the quasi-geostrophic equations produce some small-scale features (errors) in the diagnosed fields.

Primary study on pattern of general circulation of atmosphere before uplift of the Tibetan Plateau in eastern Asia

The Tibetan Plateau is a key factor for the pattern of the general circulation of the atmosphere (GCA) in eastern Asia. The pattern of the GCA after the uplift of the plateau is well known, while the pattern of the GCA before the uplift of the plateau is lack of direct evidences. Based on the knowability of desert, a section recording wind directions across the Cretaceous northern hemisphere mid-low latitude desert belt is measured and the pattern of the GCA in the Cretaceous is revealed. The result shows that the eastern Asia was really controlled by the planetary circulation before the uplift of the plateau, i.e. westerlies in the north and northeast trades in the south. The convert belt between westerlies and trades had drifted northwards and southwards. The possibility of existence of paleo-monsoon is also dealt with and a possibly imposed paleo-monsoon is suggested.

A Statistically-Based Low-Level Cloud Scheme and Its Tentative Application in a General Circulation Model

In this study, a statistical cloud scheme is first introduced and coupledwith a first-order turbulence scheme with second-order turbulence moments parameterized by thetimescale of the turbulence dissipation and the vertical turbulent diffusion coefficient. Then theability of the scheme to simulate cloud fraction at different relative humidity, verticaltemperature profile, and the timescale of the turbulent dissipation is examined by numericalsimulation. It is found that the simulated cloud fraction is sensitive to the parameter used in thestatistical cloud scheme and the timescale of the turbulent dissipation. Based on the analyses, theintroduced statistical cloud scheme is modified. By combining the modified statistical cloud schemewith a boundary layer cumulus scheme, a new statistically-based low-level cloud scheme is proposedand tentatively applied in NCAR (National Center for Atmospheric Research) CCM3 (Community ClimateModel version 3). It is found that the simulation of low-level cloud fraction is markedly improvedand the centers with maximum low-level cloud fractions are well simulated in the cold oceans off thewestern coasts with the statistically-based low-level cloud scheme applied in CCM3. It suggeststhat the new statistically-based low-level cloud scheme has a great potential in the generalcirculation model for improving the low-level cloud parameterization.

Seasonality of the Interaction between Convection over the Western Pacific and General Circulation in the Northern Hemisphere

The seasonality of the interaction between convection over the western Pacific and general circulation in the Northern Hemisphere (NH) is analyzed in the present paper with singular value decomposition (SVD) and empirical orthogonal function (EOF) analysis approaches, based on 500 hPa monthly mean geopotential height data and high-cloud amount data. The analyses demonstrate that coupled dominant patterns in the interaction between the convection over the western Pacific and the general circulation in NH are different in various seasons. In spring, the convection over the western Pacific is closely related with the western Atlantic (WA) and North Pacific (NP) like patterns of the general circulation in NH, and some associations between the WA and NP like patterns and the El Ni o /Southern Oscillation (ENSO) cycle are also existed. The Pacific Japan (PJ) pattern is the dominant pattern in the interaction between the interannual variabilities of the convection over the western Pacific and the general circulation in NH summer. The WA like pattern and 3-4 year period oscillation are also relatively obvious for the summer case. In autumn, the convection over the western Pacific is closely linked with the Eurasian (EU) like pattern and the Atlantic oscillation in the general circulation in NH, it is suggested that in autumn the variation of convective activity over the western Pacific is largely affected by the general circulation anomaly (cold air from high latitudes ) through EU like teleconnection pattern. Abrupt change happened by the end of 1980′s in the autumn interaction. The strong interaction between the western Pacific (WP) and EU like patterns in the general circulation in NH and the convection over the western Pacific and a linear trend of increasing of this interaction are also suggested in winter. It is also demonstrated that the interaction in summer and winter is stronger than in the transition seasons (spring and autumn).