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.

A modification to the Munk wind-driven ocean circulation theory

In order to fulfill the no-slip condition at the western and eastern boundaries of the ocean basin, introduced ＂effective wind stress＂, which has much larger spatial variations towards the boundaries than in the ocean interior. The effective wind stress can thus be decomposed into spatially slow-varying and fast varying components. Careful scale analysis on the classical Munk winddriven ocean circulation theory, which consists of the interior Sverdrup flow and the western boundary current but of no eastern boundary current, shows that the wind stress curl appearing in the Sverdrup equation must have negligible spatial variations. In the present model the spatially slow-varying component of the wind stress appears in the Sverdrup equation, and the spatially fastvarying component becomes the forcing term of the boundary equations. As a result, in addition to the classical Munk solution the present model has an extra term at the western boundary which （Northern Hemisphere） increases the northward transport as well as the southward return transport, and has a term at the eastern boundary corresponding to the eastern boundary current.

Wind-Driven,Double-Gyre,Ocean Circulation in a Reduced-Gravity,2.5-Layer,Lattice Boltzmann Model

A coupled lattice Boltzmann （LB） model with second-order accuracy is applied to the reduced-gravity, shallow water, 2.5-layer model for wind-driven double-gyre ocean circulation. By introducing the secondorder integral approximation for the collision operator, the model becomes fully explicit. The Coriolis force and other external forces are included in the model with second-order accuracy, which is consistent with the discretization accuracy of the LB equation. The feature of the multiple equilibria solutions is found in the numerical experiments under different Reynolds numbers based on this LB scheme. With the Reynolds number increasing from 3000 to 4000, the solution of this model is destabilized from the anti-symmetric double-gyre solution to the subtropic gyre solution and then to the subpolar gyre solution. The transitions between these equilibria states are also found in some parameter ranges. The time-dependent variability of the circulation based on this LB simulation is also discussed for varying viscosity regimes. The flow of this model exhibits oscillations with different timescales varying from subannual to interannual. The corresponding statistical oscillation modes are obtained by spectral analysis. By analyzing the spatiotemporal structures of these modes, it is found that the subannual oscillation with a 9-month period originates from the barotropic Rossby basin mode, and the interarmual oscillations with periods ranging from 1.5 years to 4.6 years originate from the recirculation gyre modes, which include the barotropic and the baroclinic recirculation gyre modes.

AN ACCURATE SOLUTION OF THE LINEAR THEORY OF THE WIND-DRIVEN OCEAN CIRCULATION-I. THE GENERALIZED SOLUTION

To model the wind-driven ocean circulation of the isobath rectangular basin,the linear vorticity equation with the meridional friction terra was used compared to the Munk'stheory on the ocean circulation. The generalized solution of the vorlicity equation was thus workedout in the sense of Fourier averaging by using the corrected Fourier expansion. The method to obtainthe undetermined coefficients was presented using the viscous boundary conditions.

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.

Tides and Wind-Driven Circulation in the Tropical and Southern Atlantic Ocean:The BRAZCOAST System

The Brazilian coast is characterized by different tidal regimes and distinct meteorological influences. The northern part has larger tidal amplitudes and is permanently affected by trade winds and tropical disturbances; the southern portion has smaller tidal amplitudes and is frequently influenced by extratropical cyclone activity. Besides these aspects, many features regarding current structure and behavior are also present, such as the equatorial system of currents, the subtropical gyre and the corresponding western boundary currents, and the Brazil-Malvinas confluence region. Within this context, efforts were made to develop the BRAZCOAST system, capable of describing the processes that determine the oceanic circulation from large to coastal scales. A customized version of the Princeton Ocean Model(POM) was implemented in a basin-scale domain covering the whole of the tropical and southern Atlantic Ocean, with 0.5° spatial resolution, as well as three nested grids with(1/12)° resolution covering the different parts of the Brazilian shelf, in a one-way procedure. POM was modified to include tidal potential generator terms and a partially-clamped boundary condition for tidal elevations. The coarse grid captured large-scale features, while the nested grids detailed local circulations affected by bathymetry and coastal restrictions. An interesting aspect at the coarse grid level was the relevance of the Weddell Sea to the location of the tidal amphidromic systems.

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.

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.

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.

A Note on the South China Sea Shallow Interocean Circulation

The existing estimates of the volume transport from the Pacific Ocean to the South China Sea are summarized, showing an annual mean westward transport, with the Taiwan Strait outflow subtracted, of 3.5±2.0 Sv （1 Sv=-0^6 ma s^-1）. Results of a global ocean circulation model show an annual mean transport of 3.9 Sv from the Pacific to the Indian Ocean through the South China Sea. The boreal winter transport is larger and exhibits a South China Sea branch of the Pacific-to-Indian Ocean throughflow, which originates from the western Philippine Sea toward the Indonesian Seas through the South China Sea, as well as through the Karimata and Mindoro Straits. The southwestward current near the continental slope of the northern South China Sea is shown to be a combination of this branch and the interior circulation gyre. This winter branch can be confirmed by trajectories of satellite-tracked drifters, which clearly show a flow from the Luzon Strait to the Karimata Strait in winter. In summer, the flow in the Karimata Strait is reversed. Numerical model results indicate that the Pacific water can enter the South China Sea and exit toward the Sulu Sea, but no observational evidence is available. The roles of the throughiiow branch in the circulation, water properties and air-sea exchange of the South China Sea, and in enhancing and regulating the volume transport and reducing the heat transport of the Indonesian Throughflow, are discussed.

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.

Interannual Variability of Late-spring Circulation and Diabatic Heating over the Tibetan Plateau Associated with Indian Ocean Forcing

The thermal forcing of the Tibetan Plateau （TP） during boreal spring, which involves surface sensible heating, latent heating released by convection and radiation flux heat, is critical for the seasonal and subseasonal variation of the East Asian summer monsoon. Distinct from the situation in March and April when the TP thermal forcing is modulated by the sea surface temperature anomaly （SSTA） in the North Atlantic, the present study shows that it is altered mainly by the SSTA in the Indian Ocean Basin Mode （IOBM） in May, according to in-situ observations over the TP and MERRA reanalysis data. In the positive phase of the IOBM, a local Hadley circulation is enhanced, with its ascending branch over the southwestern Indian Ocean and a descending one over the southeastern TP, leading to suppressed precipitation and weaker latent heat over the eastern TP. Meanwhile, stronger westerly flow and surface sensible heating emerges over much of the TP, along with slight variations in local net radiation flux due to cancellation between its components. The opposite trends occur in the negative phase of the IOBM. Moreover, the main associated physical processes can be validated by a series of sensitivity experiments based on an atmospheric general circulation model, FAMIL. Therefore, rather than influenced by the remote SSTAs of the northern Atlantic in the early spring, the thermal forcing of the TP is altered by the Indian Ocean SSTA in the late spring on an interannual timescale.

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.

The Shallow Meridional Overturning Circulation in the Northern Indian Ocean and Its Interannual Variability

The shallow meridional overturning circulation (upper 1000 m) in the northern Indian Ocean and its interannual variability are studied, based on a global ocean circulation model (MOM2) with an integration of 10 years (1987-1996). It is shown that the shallow meridional overturning circulation has a prominent seasonal reversal characteristic. In winter, the flow is northward in the upper layer and returns southward at great depth. In summer, the deep northward inflow upwells north of the equator and returns southward in the Ekman layer. In the annual mean, the northward inflow returns through two branches: one is a southward flow in the Ekman layer, the other is a flow that sinks near 10°N and returns southward between 500 m and 1000 m. There is significant interannual variability in the shallow meridional overturning circulation, with a stronger (weaker) one in 1989 (1991) and with a period of about four years. The interannual variability of the shallow meridional overturning circulation is intimately r

Effect of Langmuir circulation on upper ocean mixing in the South China Sea

Effect of Langmuir circulation （LC） on upper ocean mixing is investigated by a two-way wave-current coupled model. The model is coupled of the ocean circulation model ROMS （regional ocean modeling system） to the surface wave model SWAN （simulating waves nearshore） via the model-coupling toolkit. The LC already certified its importance by many one-dimensional （1D） research and mechanism analysis work. This work focuses on inducing LC＇s effect in a three-dimensional （3-D） model and applying it to real field modeling. In ROMS, the Mellor-Yamada turbulence closure mixing scheme is modified by including LC＇s effect. The SWAN imports bathymetry, free surface and current information from the ROMS while exports signifi- cant wave parameters to the ROMS for Stokes wave computing every 6 s. This coupled model is applied to the South China Sea （SCS） during September 2008 cruise. The results show that LC increasing turbulence and deepening mixed layer depth （MLD） at order of O （10 m） in most of the areas, especially in the north part of SCS where most of our measurements operated. The coupled model further includes wave break- ing which will brings more energy into water. When LC works together with wave breaking, more energy is transferred into deep layer and accelerates the MLD deepening. In the north part of the SCS, their effects are more obvious. This is consistent with big wind event in the area of the Zhujiang River Delta. The shallow water depth as another reason makes them easy to influence the ocean mixing as well.

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.

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.

Numerical simulation of wind-driven circulation and pollutant transport in Taihu Lake based on a quadtree grid

In this study,a two-dimensional flow-pollutant coupled model was developed based on a quadtree grid.This model was established to allow the accurate simulation of wind-driven flow in a large-scale shallow lake with irregular natural boundaries when focusing on important smallscale localized flow features.The quadtree grid was created by domain decomposition.The governing equations were solved using the finite volume method,and the normal fluxes of mass,momentum,and pollutants across the interface between cells were computed by means of a Godunov-type Osher scheme.The model was employed to simulate wind-driven flow in a circular basin with non-uniform depth.The computed values were in agreement with analytical data.The results indicate that the quadtree grid has fine local resolution and high efficiency,and is convenient for local refinement.It is clear that the quadtree grid model is effective when applied to complex flow domains.Finally,the model was used to calculate the flow field and concentration field of Taihu Lake,demonstrating its ability to predict the flow and concentration fields in an actual water area with complex geometry.

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.