Effects of the surface wave-induced mixing on circulation in an isopycnal-coordinate oceanic circulation model

The influence of the nonbreaking surface wave-induced mixing under the mixed layer on the oceanic cir- culation was investigated using an isopycnal-coordinate oceanic circulation model. The effect of the wave- induced mixing within the mixed layer was eliminated via a bulk mixed layer model. The results show that the wave-induced mixing can penetrate through the mixed layer and into the oceanic interior. The wave- induced mixing under the mixed layer has an important effect on the distribution of temperature of the upper ocean at middle and high latitudes in summer, especially the structure of the seasonal thermocline. Moreover, the wave-induced mixing can affect the oceanic circulation, such as western boundary currents and the North Equatorial Currents through changes of sea surface height associated with the variation of the thermal structure of the upper ocean.

Advances in Studying Oceanic Circulation from Hydrographic Data with Applications in the South China Sea

Methods for studying oceanic circulation from hydrographic data are reviewed in the context of their applications in the South China Sea. These methods can be classified into three types according to their different dynamics as follows: (1) descriptive methods, (2) diagnostic methods without surface and bottom forcing, and (3) diagnostic methods with the above boundary forcing. The paper discusses the progress made in the above methods together with the advancement of study in the South China Sea circulation.

Numerical simulations of Atlantic meridional overturning circulation(AMOC)from OMIP experiments and its sensitivity to surface forcing

Atlantic meridional overturning circulation(AMOC)plays an important role in transporting heat meridionally in the Earth’s climate system and is also a key metrical tool to verify oceanic general circulation models.Two OMIP(Ocean Model Intercomparison Project phase 1 and 2)simulations with LICOM3(version 3 of the LASG/IAP Climate System Ocean Model)developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics(LASG),Institute of Atmospheric Physics(IAP),are compared in this study.Both simulations well reproduce the fundamental characteristics of the AMOC,but the OMIP1 simulation shows a significantly stronger AMOC than the OMIP2 simulation.Because the LICOM3 configurations are identical between these two experiments,any differences between them must be attributed to the surface forcing data.Further analysis suggests that sea surface salinity(SSS)differences should be mainly responsible for the enhanced AMOC in the OMIP1 simulation,but sea surface temperature(SST)also play an unignorable role in modulating AMOC.In the North Atlantic,where deep convection occurs,the SSS in OMIP1 is more saline than that in OMIP1.We find that in the major region of deep convection,the change of SSS has more significant effect on density than the change of SST.As a result,the SSS was more saline than that in OMIP2,leading to stronger deep convection and subsequently intensify the AMOC.We conduct a series of numerical experiments with LICOM3,and the results confirmed that the changes in SSS have more significant effect on the strength of AMOC than the changes in SST.

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.

Low-temperature alteration of oceanic island basalts and their contribution to transition metal circulation of the ocean

The major elements, rare earth elements （REE） and trace elements of four basalt samples from central and western Pacific ferro- manganese crust provinces have been analyzed using chemical methods and ICP - MS, respectively. The results indicate that the samples have been extensively altered and that the contents of their major elements have changed significantly. However, the similarity of REE partition patterns and trace element contents of basalt samples to those of fresh oceanic island basalts （OIB） indicate that the basalt samples originated as OIB. Because of low-temperature alteration, the contents of A1203 , Fe203 , MnO, K20 and P205 increased, while MgO and FeO decreased. Active components, such as magnesium and iron, were leached from OIB resulting in the relative enrichment of SiO2. The leaching of active components can cause the relative enrichment of REE, while the precipitation of LREE-rich ferromanganese oxides in vesicles and fissures not only causes an increase of REE contents, but also induces ＂fractionation＂ of LREE and HREE. Based on the enrichment mechanism of REE contents, the theoretical quantities of precipitated ferromanganese oxides and the depleted quantities of active components are calculated ： the depleted quantities of active components for the unit mass of fresh basalts vary in the range of 0.15 ～ 0. 657, and the precipitated quantities of ferromanganese oxides for the unit mass of fresh basahs vary in the range of 0. 006 ～ 0. 042. Of the major elements, the two most depleted are iron, and magnesium, with 18.28% ～ 70.95% of iron and 44.50% ～ 93.94% of magnesium in the fresh basalts was leached out. Theoretical calculation and geochemistry results both indicate that low-temperature alteration of basalts can supply abundant amount of metals to seawater, and may play an important role in ocean metal circulation.

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

Numerical Simulation of the Regional Ocean Circulation in the Coastal Areas of China

The regional ocean circulation in the coastal areas of China (including a part of the western PacificOcean, the South China Sea and the Bay of Bengal et al.) is simulated by using the improved Princeton University ocean circulation model (POM). Compared with the modeling results obtained by the large-scaleocean general circulation model (OGCM), the basic ocean circulation features simulated by the regionalocean circulation model al-e in good agreement with that simulated by OGCM and some detailed characteristics such as the regional ocean circulation, sea temperature, salinity and flee sea surface height have alsobeen obtained which are in good accord with the observations. These results indicate that the regional oceancirculation model has good capability to produce the legional ocean circulation characteristics and it can beused to develop coupled legional ocean-atmospheric model systems.

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.

Abyssal Circulation in the Philippine Sea

The abyssal circulation in the Philippine Sea(PS)is investigated,with outputs from the Simple Ocean Data Assimilation version 2.2.4(SODA224).The deep-water currents in SODA224 are carefully evaluated,with sparse in situ observations in the North Pacific Ocean.In the upper deep layer(20003000 m)of the PS,a strong westward current,which originates from the Northeast Pacific Basin and enters the PS through the Yap-Mariana Junction,exists along 1114 N.This strong westward current bifurcates into two western boundary currents off the Philippines.The northward-flowing current flows out of the PS around 2021 N,whereas the southward-flowing current transports deep water from the northern hemisphere to the southern hemisphere.In the lower deep layer(30004500 m),the inflow water first flows northward to the east of the Western Mariana Basin and then turns westward at approximately 18 N.The inflow water mainly enters the Philippine Basin(PB),with a small part turning southward to constitute a weak cyclonic circulation.The water entering the PB mainly merges into a strong southward western boundary current in the south-ern PB.In the bottom layer(below 4500 m),both the northeast and northwest PB show single cyclonic gyres,whereas the south PB shows a single anticyclonic gyre.Moreover,comparisons with the observations indicate the possible existence of a cyclonic sense of circulation over the Philippine Trench.The current study provides the implications for future observations,which are needed to fur-ther investigate the temporospatial variations of the abyssal circulation in the PS on multiple scales.

Oceanic Climatology in the Coupled Model FGOALS-g2: Improvements and Biases

The present study examines simulated oceanic climatology in the Flexible Global Ocean-Atmosphere- Land System model, Grid-point Version 2 （FGOALS-g2） forced by historical external forcing data. The oceanic temperatures and circulations in FGOALS-g2 were found to be comparable to those observed, and substantially improved compared to those simulated by the previous version, FGOALS-gl.0. Compared with simulations by FGOALS-gl.0, the shallow mixed layer depths were better captured in the eastern Atlantic and Pacific Ocean in FGOALS-g2. In the high latitudes of the Northern Hemisphere, the cold biases of SST were about 1℃-5℃ smaller in FGOALS-g2. The associated sea ice distributions and their seasonal cycles were more realistic in FGOALS-g2. The pattern of Atlantic Meridional Overturning Circulation （AMOC） was better simulated in FGOALS-g2, although its magnitude was larger than that found in observed data. The simulated Antarctic Circumpolar Current （ACC） transport was about 140 Sv through the Drake Passage, which is close to that observed. Moreover, Antarctic Intermediate Water （AAIW） was better captured in FGOALS-g2. However, large SST cold biases （〉3℃） were still found to exist around major western boundary currents and in the Barents Sea, which can be explained by excessively strong oceanic cold advection and unresolved processes owing to the coarse resolution. In the Indo-Pacific warm pool, the cold biases were partly related to the excessive loss of heat from the ocean. Along the eastern coast in the Atlantic and Pacific Oceans, the warm biases were due to overestimation of shortwave radiation. In the Indian Ocean and Southern Ocean, the surface fresh biases were mainly due to the biases of precipitation. In the tropical Pacific Ocean, the surface fresh biases （〉2 psu） were mainly caused by excessive precipitation and oceanic advection. In the Indo-Pacific Ocean, fresh biases were also found to dominate in the upper 1000 m, except in the northeastern Indian Ocean. There were warm and salty biases （3℃-4℃ and 1-2 psu） from the surface to the bottom in the Labrador Sea, which might be due to large amounts of heat transport and excessive evaporation, respectively. For vertical structures, the maximal biases of temperature and salinity were found to be located at depths of 〉600 m in the Arctic Ocean, and their values exceeded 4℃ and 2 psu, respectively.

Intercomparison of three South China Sea circulation models

Three numerical oceanic circulation models: POM(Princeton ocean model), MICOM(Miami isopycnal coordinates ocean model) and GFDL model, which adopt sigma coordinate, isopycnal coordinate and depth coordinate respectively, are used in the South China Sea(SCS) circulation modeling. Model domain has the same topography, grid resolution, initial conditions and surface boundary conditions. The maximum ocean depth is set as 1 000 m. Grid resolution is 0.5o×0.5o.Initial conditions are supplied by climatological temperature and salinity data in January. Climatological wind stress, surface temperature and salinity are used as surface forcing. Lateral boundaries take enclosed boundary conditions artificially. Focusing on the common point of different ocean circulation models, the circulation pattern in winter and summer, sea surface height in the northern SCS, seasonal cycle of the mixed layer thickness in the southern SCS, barotropic stream function in winter are selected to carry out intercomparison. In winter, a strong cyclonic gyre occupies the whole SCS. In summer, a strong anticyclonic gyre occupies the southern SCS and a weak cyclonic gyre occupies the northern SCS. The thickness of the mixed layer shows bimodal features in the southern SCS. Sea surface height anomaly(SSHA) in the northern SCS has an eastward propagating feature, in agreement with the remote sensing observation. Barotropic stream functions indicate that the circulation of the upper ocean is mainly forced by inputting of wind stress curl under closed boundary conditions. In addition, three models also show distinct differences. The basin-scale circulation from MICOM is distinct. Output of POM has more mesoscale eddies than others. GFDL model seems good at simulating westward intensification.

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.

Design and Numerical Simulation of an Arctic Ocean Circulation and Thermodynamic Sea－Ice Model

In this paper, the first version of a new Arctic Ocean circulation and thermodynamic sea-ice model is presentedby the authors based on the framework of a twenty-layer World Oceanic general circulation model developed byZhang et al. in 1994. The model's domain covers the Arctic Ocean and Greenland-Norwegian Seas with the horizontal resolution of 200 km×200 km on a stereographic projection plane. In vertical, the model uses the Eta-coordinate(Sigma modified to have quasi-horizontal coordinate surfaces) and has ten unevenly-spaced layers to cover the deepest water column of 3000 m. Two 150-year integrations of coupling the ocean circulation model with the sea-icemodel have been performed with seasonally cyclic surface boundary conditions. The only difference between the tWoexperiments is in the model's geography. Some preliminary analyses of the experimental results have been done focused on the following aspects: (1) surface layer temperature, salinity and current; (2) the' Atlantic Layer'; (3)sea-ice cover and its seasonal variation. In comparison with the available observational data, these results are acceptable with reasonable accuracy.

A two-time-level split-explicit ocean circulation model (MASNUM) and its validation

A two-time-level, three-dimensional numerical ocean circulation model（named MASNUM） was established with a two-level, single-step Eulerian forward-backward time-differencing scheme. A mathematical model of large-scale oceanic motions was based on the terrain-following coordinated, Boussinesq, Reynolds-averaged primitive equations of ocean dynamics. A simple but very practical Eulerian forward-backward method was adopted to replace the most preferred leapfrog scheme as the time-differencing method for both barotropic and baroclinic modes. The forward-backward method is of second-order of accuracy, computationally efficient by requiring only one function evaluation per time step, and free of the computational mode inherent in the three-level schemes. This method is superior to the leapfrog scheme in that the maximum time step of stability is twice as large as that of the leapfrog scheme in staggered meshes thus the computational efficiency could be doubled. A spatial smoothing method was introduced to control the nonlinear instability in the numerical integration. An ideal numerical experiment simulating the propagation of the equatorial Rossby soliton was performed to test the amplitude and phase error of this new model. The performance of this circulation model was further verified with a regional（northwest Pacific） and a quasi-global（global ocean simulation with the Arctic Ocean excluded） simulation experiments. These two numerical experiments show fairly good agreement with the observations. The maximum time step of stability in these two experiments were also investigated and compared between this model and that model which adopts the leapfrog scheme.

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.

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.

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.

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.

Wind-Driven Ocean Circulation in Shallow Water Lattice Boltzmann Model

A lattice Boltzmann (LB) model with overall second-order accuracy is applied to the 1.5-layer shallow water equation for a wind-driven double-gyre ocean circulation. By introducing the second-order integral approximation for the collision operator, the model becomes fully explicit. In this case, any iterative technique is not needed. The Coriolis force and other external forces are included in the model with second-order accuracy, which is consistent with the discretized accuracy of the LB equation. The numerical results show correct physics of the ocean circulation driven by the double-gyre wind stress with different Reynolds numbers and different spatial resolutions. An intrinsic low-frequency variability of the shallow water model is also found. The wind-driven ocean circulation exhibits subannual and interannual oscillations, which are comparable to those of models in which the conventional numerical methods are used.