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.

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.

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.

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.

SEASONAL HEAT BUDGET IN THE TROPICAL WESTERN PACIFIC OCEAN IN A GLOBAL GCM I. GENERAL DESCRIPTION*

This paper describes the large scale aspects of the seasonal surface heat budget and discusses its main forcing mechanisms in the tropical Western Pacific Ocean. The high - resolution general circulation model (Semtner & Chervin, 1992) used in this study reproduced well the observed upper-layer thermal structure and circulation. It is shown that at least on the average of the study region (20 °N -20°N, west boundary-160 °E) the semiannual variation is a dominant signal for all heat budget components and is presumably due to the sun's passing across the equator twice a year; but that the components have substantial differences in amplitude. The local Ekman divergence in the region does not change significantly through the year. As a result, the change in surface heat content is roughly half due to ocean ?atmosphere heat exchange and half due to heat advection by remotely forced verti-cal motion. Horizontal currents do not play a significant role directly by advection, because the wat-er which enters the region is not very much different in temperature from the water which leaves it.

Bottom water temperature measurements in the South China Sea,eastern Indian Ocean and western Pacific Ocean

文章报道了一批新的海底底水温度(BWT)数据,其中南海(SCS)158个站位、东印度洋(EIO)30个站位及西太平洋(WPO)37个站位。基于这批新的BWT数据,获得南海和西太平洋海域底水温度与水深经验关系,可为地球物理和物理海洋提供准确、可靠的海底温度边界。这将有助于海底油气资源调查与评估。同时,这批实测数据表明:1)水深超过3500m的海域,其底水温度在南海约为2.47℃,比东印度洋(～1.34℃)和西太平洋(～1.60℃)稍微偏高。这与大洋传送带模式所预测的情况比较吻合。该模式认为:低温高盐的海水,从北大西洋格陵兰岛和冰岛附近海域下沉到深层,然后向南流动,再与南极洲周围海域的低温高盐海水一同向北进入印度洋和太平洋。而南海是一个相对比较封闭的热带边缘海,其内部海水与印度洋和菲律宾海交换有限,导致海水温度整体高于印度洋和太平洋。2)台西南盆地水深在2700～3000m的部分站位,其底水温高达约3.00℃,明显高于其周边同水深海域底水温度(平均值约为2.33℃)。这可能是台西南盆地海底水热活动导致的结果。3)在东印度洋和西太平洋水深超过4800m海域,底水温度随着水压增大稍有升高,其升高率分别为10.6mK·MPa^(-1)和12.0mK·MPa^(-1)。这与理论估算的深层底水绝热压力温度梯度范围较为吻合。这也意味着东印度洋和西太平洋深层底水,主要由绝热自压作用导致其温度随着深度的增大而升高。

SEASONAL HEAT BUDGET IN THE TROPICAL WESTERN PACIFIC OCEAN IN A GLOBAL GCM Ⅱ. IN FIVE SUBREGIONS

The general features of the seasonal suuface heat budget in the tropical western Pacific Ocean,20°S-20°N, western boundary-160°E, were documented by Qu (1995) using a high-resolution generalcirculation model (GCM, Semtner & Chervin,1992) ard existing observations.Close inspection of thesmaller areas, with the whole region further partitioned into six parts, showed different mechanisms balancethe seasonal surface heat budget in different parts of the region The results of study on five subregionsare detailed in this article. In the equatorial (3°S - 3°N) aed North Equatorial Countercurrent(3°N-9°N) region, the surface the flux the does not change significantly throughout the year, so the surface heat content is determined largely by vertical motion near the equator and roughly helf due to horizontal and halfdue to vertical circulation in the region of the North Equatorial Countercurrent(NECC). In the othersubregions (9°N-20°N, 20°S -11°S aed 11°S -3°S ), however, in addition to ocean

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

On the basis of Zeng's theorehcal design, a coupled general circulation model(CGCM) is develO￣ with itscharacteristics 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 30W and 30'S with horizontal grid spacing of ic in latitude and 2°in longitude,and with 14 vertical layers.The atmospheric component is a global GCM with low-resolution of 4°in lahtude and 5°in longitude,and tWo layers of equal mass in the verhcal between the surfaCe and 200 hFa.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 correchon,and to so-called predictor-corrector monthly coupling(PCMC),are conducted tO show the esistence 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 anneal mean and its seasollal cycle of the atmospheric and￣ante circulahon.The model also ProduCeS the coherent intermnual variations of the climate system, manifesting the observed EI Nifio/Southern OSCillation(ENSO).

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.

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

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.

An extended variable-grid global ocean circulation model and its preliminary results of the equatorial Pacific circulation

To investigate the interaction between the tropical Pacific and China seas a variable-grid global ocean circulation model with fine grid covering the area from 20°S to 50°N and from 99° to 150°E is developed. Numerical computation of the annually cyclic circulation fields is performed. The results of the annual mean zonal currents and deep to abyssal western boundary currents in the equatorial Pacific Ocean are reported. The North Equatorial Current,the North Equatorial Countercurrent, the South Equatorial Current and the Equatorial Undercurrent are fairly well simulated. The model well reproduces the northward flowing abyssal western boundary current.From the model results a lower deep western boundary current east of the Bismarck-Solomon-New Hebrides Island chain at depths around 2 000 m has been found. The model results also show that the currents in the equatorial Pacific Ocean have multi-layer structures both in zonal currents and western boundary currents, indicating that the global ocean overturning thermohaline circulation appears of multi-layer pattern.

Review on observational studies of western tropical Pacific Ocean circulation and climate

The Western Tropical Pacific(WTP) Ocean holds the largest area of warm water(>28℃) in the world ocean referred to as the Western Pacific Warm Pool(WPWP),which modulates the regional and global climate through strong atmospheric convection and its variability.The WTP is unique in terms of its complex 3-D ocean circulation system and intensive multiscale variability,making it crucial in the water and energy cycle of the global ocean.Great advances have been made in understanding the complexity of the WTP ocean circulation and associated climate impact by the international scientific community since the 1960 s through field experiments.In this study,we review the evolving insight to the 3-D structure and multi-scale variability of the ocean circulation in the WTP and their climatic impacts based on in-situ ocean observations in the past decades,with emphasis on the achievements since 2000.The challenges and open que stions remaining are reviewed as well as future plan for international study of the WTP ocean circulation and climate.

A THREE-LAYER OCEAN CIRCULATION MODEL APPLICATION TO NUMERICAL STUDIES ON THE NORTH PACIFIC OCEAN CIRCULATION

The formulation and justification of a three-layer baroclinic ocean model developed to simulate thegeneral circulation of the ocean are described in this paper.Test of the model in simulating the annualmean circulation patterns in the North Pacific under the prescribed atmospheric forcing,which consists ofthe climatological surface wind stress and sea surface heat flux,and comparison of the results withobservations showed that the model basically simulated the large scale features of the annual meancirculation patterns in the North Pacific Ocean such as those of the intensified western boundary currentsand the North Equatorial Currents and Undercurrents.But due to the coarse resolution of the model,some details of these currents were poorly reproduced.The seasonal variations of the North Pacific Oceancirculation driven by the seasonal mean sea surface wind stress was calculated,the different aspects of theseresults were analyzed and the main current(the intensified western boundary currents)transports

A quasi-global 1/10° eddy-resolving ocean general circulation model and its preliminary results

In this study,a quasi-global(excluding the Arctic Ocean) eddy-resolving ocean general circulation model(OGCM) is established based on the latest version of the LASG/IAP Climate system Ocean Model(LICOM2.0).The horizontal resolution and vertical resolution are increased to 1/10° and 55 layers,respectively.Forced by the surface fluxes from the reanalysis and observed data,the model has been integrated for approximately 20 model years(20 a).During the last 8 a,the model is driven by daily mean wind stresses from QuikSCAT and heat fluxes from reanalysis data from 2000 to 2007.The basic performance of the OGCM is analyzed using the last 8 a simulation output.Compared with the simulation of the coarse-resolution OGCM,the eddy-resolving OGCM not only better simulates the spatial-temporal features of mesoscale eddies and the paths and positions of western boundary currents but also reproduces the large meander of the Kuroshio Current and its interannual variability.Another aspect,namely,the complex structures of equatorial Pacific currents and currents in the coastal ocean of China,are better captured due to the increased horizontal and vertical resolution.

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.

The numerical simulation of Pacific Ocean circulation with five-level model

A five-level oceanic primitive equation model has been developed. This model is integrated numerically with annual mean wind stress and heat flux at sea surface for 30 a. The ocean circulations tend to quasi-stability. The simulated results show that the computed annual mean currents and sea surface temperature agree well with the observations.

The influence of explicit tidal forcing in a climate ocean circulation model

The eight main tidal constituents have been implemented in the global ocean general circulation model with approximate 1° horizontal resolution.Compared with the observation data,the patterns of the tidal amplitudes and phases had been simulated fairly well.The responses of mean circulation,temperature and salinity are further investigated in the global sense.When implementing the tidal forcing,wind-driven circulations are reduced,especially those in coastal regions.It is also found that the upper cell transport of the Atlantic meridional overturning circulation（AMOC） reduces significantly,while its deep cell transport is slightly enhanced from 9×106m3/s to 10×106 m3/s.The changes of circulations are all related to the increase of a bottom friction and a vertical viscosity due to the tidal forcing.The temperature and salinity of the model are also significantly affected by the tidal forcing through the enhanced bottom friction,mixing and the changes in mean circulation.The largest changes occur in the coastal regions,where the water is cooled and freshened.In the open ocean,the changes are divided into three layers：cooled and freshened on the surface and below 3 000 m,and warmed and salted in the middle in the open ocean.In the upper two layers,the changes are mainly caused by the enhanced mixing,as warm and salty water sinks and cold and fresh water rises;whereas in the deep layer,the enhancement of the deep overturning circulation accounts for the cold and fresh changes in the deep ocean.

Tide simulation in a global eddy-resolving ocean model

The tide plays a pivotal role in the ocean,affecting the global ocean circulation and supplying the bulk of the energy for the global meridional overturning circulation.To further investigate internal tides and their impacts on circulation,it is imperative to incorporate tidal forcing into the eddy-resolving global ocean circulation model.In this study,we successfully incorporated explicit tides(eight major constituents)into a global eddy-resolving general ocean circulation model and evaluated its tidal simulation ability.We obtained harmonic constants by analyzing sea surface height through tidal harmonic analysis and compared them with the analysis data Topex Poseidon Cross-Overs v9(TPXO9),the open ocean tide dataset from 102 open-ocean tide observations,and tide gauge stations from World Ocean Circulation Experiment.The results demonstrated that the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics(LASG/IAP)Climate System Ocean Model 3.0(LICOM3.0)effectively simulated tides,with errors predominantly occurring in nearshore regions.The tidal amplitude simulated in LICOM3.0 was greater than that of TPXO9,and these high-amplitude areas exhibited greater errors.The amplitude error of the M_(2) constituent was larger,while the phase error of the K_(1) constituent was more significant.Furthermore,we further compared our results with those from other models.

Seasonal variability of the isopycnal surface circulation in the South China Sea derived from a variable-grid global ocean circulation model

In this study, we develop a variable-grid global ocean general circulation model （OGCM） with a fine grid （1/6）° covering the area from 20°S-50°N and from 99°-150°E, and use the model to investigate the isopycnal surface circulation in the South China Sea （SCS）. The simulated results show four layer structures in vertical： the surface and subsurface circulation of the SCS are characterized by the monsoon driven circulation, with basin-scaled cyclonic gyre in winter and anti-cyclonic gyre in summer. The intermediate layer circulation is opposite to the upper layer, showing anti-cyclonic gyre in winter but cyclonic gyre in summer. The circulation in the deep layer is much weaker in spring and summer, with the maximum velocity speed below 0.6 cm/s. In fall and winter, the SCS deep layer circulation shows strong east boundary current along the west coast of Philippine with the velocity speed at 1.5 m/s, which flows southward in fall and northward in winter. The results have also revealed a fourlayer vertical structure of water exchange through the Luzon Strait. The dynamics of the intermediate and deep circulation are attributed to the monsoon driving and the Luzon Strait transport forcing.