Diagnostic Calculation of the Oceanic Circulation

A 2-dimensional global free surface diagnostic model, combined with dynamic calculation, is used to investigate the world ocean circulation; the model has a horizontal resolution of 1/4°×1/4°. The simulated results agree well with the results of other modesl and observations. The distribution of Stream Function suggests that the main circulation systems in the wodd ocean have been represented, including oceanic currents strengthened in the oceanic western. Be close to the observed results, the net mass transport of the Kuroshio axes is estimated about 54Sv; The distribution of the horizontal circulation in each layer shows that the main circulation systems in the world ocean are well simulated, for example, the Kuroshio and the Antarctic Circumpolar Current can go down to the bottom layer, but the Gulf Stream cannot, and its direction reverses at the depths of 1000 to 2 000 m.

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.

Quantifying the Role of the Eddy Transfer Coefficient in Simulating the Response of the Southern Ocean Meridional Overturning Circulation to Enhanced Westerlies in a Coarse-resolution Model

This study assesses the capability of a coarse-resolution ocean model to replicate the response of the Southern Ocean Meridional Overturning Circulation(MOC) to intensified westerlies,focusing on the role of the eddy transfer coefficient(κ).κ is a parameter commonly used to represent the velocities induced by unresolved eddies.Our findings reveal that a stratification-dependent κ,incorporating spatiotemporal variability,leads to the most robust eddy-induced MOC response,capturing 82% of the reference eddy-resolving simulation.Decomposing the eddy-induced velocity into its vertical variation(VV) and spatial structure(SS) components unveils that the enhanced eddy compensation response primarily stems from an augmented SS term,while the introduced VV term weakens the response.Furthermore,the temporal variability of the stratification-dependent κ emerges as a key factor in enhancing the eddy compensation response to intensified westerlies.The experiment with stratification-dependent κ exhibits a more potent eddy compensation response compared to the constant κ,attributed to the structure of κ and the vertical variation of the density slope.These results underscore the critical role of accurately representing κ in capturing the response of the Southern Ocean MOC and emphasize the significance of the isopycnal slope in modulating the eddy compensation mechanism.

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.

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.

Simulation of Pacific Ocean Circulations Based on Global Warming from 1960 to 1999

The Pacific Ocean circulations were simulated based on the global warming from 1960 to 1999 by using the Non-Boussinesq POP model and the data of wind stress and temperature at 1 000 hPa from the NCEP. The results show that the circulation in the tropical Pacific Ocean was weakening during the past 40 years. The heat transported to the tropical western Pacific Ocean coast by the north equatorial current and the heat transported to middle and high latitudes in the southem hemisphere by the south equatorial current decreased with time due to the global warming, while the heat transported to middle and high latitudes in the northern hemisphere by the north equatorial current increased with time due to the global warming.

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.

The effect of Coriolis-Stokes forcing on upper ocean circulation in a two-way coupled wave-current model

We investigated the Stokes drift-driven ocean currents and Stokes drift-induced wind energy input into the upper ocean using a two-way coupled wave-current modeling system that consists of the Princeton Ocean Model generalized coordinate system (POMgcs), Simulating WAves Nearshore (SWAN) wave model, and the Model Coupling Toolkit (MCT). The Coriolis-Stokes forcing (CSF) computed using the wave parameters from SWAN was incorporated with the momentum equation of POMgcs as the core coupling process. Experimental results in an idealized setting show that under the steady state, the scale of the speed of CSF-driven current was 0.001 m/s and the maximum reached 0.02 m/s. The Stokes drift-induced energy rate input into the model ocean was estimated to be 28.5 GW, taking 14% of the direct wind energy rate input. Considering the Stokes drift effects, the total mechanical energy rate input was increased by approximately 14%, which highlights the importance of CSF in modulating the upper ocean circulation. The actual run conducted in Taiwan Adjacent Sea (TAS) shows that: 1) CSF-based wave-current coupling has an impact on ocean surface currents, which is related to the activities of monsoon winds; 2) wave-current coupling plays a significant role in a place where strong eddies present and tends to intensify the eddy's vorticity; 3) wave-current coupling affects the volume transport of the Taiwan Strait (TS) throughflow in a nontrivial degree, 3.75% on average.

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.

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.

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.