Resuspension and advection processes affecting suspended particulate matter concentrations in the central English Channel

Suspended particulate matter (SPM) measurements obtained along a cross-section in the central English Channel (Wight-Cotentin transect) indicate that the area may be differentiated into: (1) an English coastal zone, associated with the highest concentrations; (2) a French coastal zone, with intermediate concentrations; and (3) the offshore waters of the Channel, characterised by a very low suspended-sediment load. The SPM particle-size distribution was modal close to the English coast (main mode 10-12 mu m); the remainder of the area was characterised by flat SPM distributions. Examination of the diatom communities in the SPM suggest:; that material resuspended in the intertidal zone and the estuarine environments was advected towards the offshore waters of the English Channel. Considerable variations in SPM concentrations occurred during a tidal cycle: maximum concentrations were sometimes up to 3 times higher than the minimum concentrations, Empirical orthogonal function (EOF) analysis of the SPM concentration time series indicates that, although the bottom waters were more turbid than the surficial waters, this was not likely to be the result of in situ sediment resuspension. Instead, the observed variations appear to be controlled mainly by advective mechanisms. The limited resuspension was probably caused by: (1) the limited availability of fine-grained material within the bottom sediments, and (2) 'bed-armouring' processes which protect the finer-grained fractions of the seabed material from erosion and entrainment within the overlying flow during the less energetic stages of the tide.

Interannual horizontal heat advection in the surface mixed layer over the equatorial Pacific Ocean: assimilation versus TAO analysis

The interannual anomalies of horizontal heat advection in the surface mixed layer over the equatorial Pacific Ocean in an assimilation experiment are studied and compared with existing observational analyses. The assimilation builds upon a hindcast study that has produced a good simulation of the observed equatorial currents and optimizes the simulation of the Reynolds sea surface temperature (SST) data. The comparison suggests that the assimilation has improved the simulation of the interannual horizontal heat advection of the surface mixed layer significantly. During periods of interrupted current measurements, the assimilation is shown to produce more meaningful anomalies of the heat advection than the interpolation of the observational data does. The assimilation also shows that the eddy heat flux due to the correlation between high-frequency current and SST variations, which is largely overlooked by the existing observational analyses, is important for the interannual SST balance over the equatorial Pacific. The interannual horizontal heat advection anomalies are found to be sensitive to SST errors where oceanic currents are strong, which is a challenge for ENSO prediction. The study further suggests that the observational analyses of the tropical SST balance based on the TAO and the Reynolds SST data contain significant errors due to the large gradient errors in the Reynolds SST data, which are amplified into the advection anomalies by the large equatorial currents.

A Lagrangian lattice Boltzmann method for Euler equations

A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.

Simulations of oil spread on a water surface

Numerical simulation of an oil slick spreading on still and wavy surfaces is described in this paper. The so-called sigma transformation is used to transform the time-varying physical domain into a fixed calculation domain for the water wave motions and, at the same time, the continuity equation is changed into an advection equation of wave elevation. This evolution equation is discretized by the forward time and central space scheme, and the momentum equations by the projection method. A damping zone is set up in front of the outlet boundary coupled with a Sommerfeld-Orlanski condition at that boundary to minimize the wave reflection. The equations for the oil slick are depth-averaged and coupled with the water motions when solving numerically. As examples, sinusoidal and solitary water waves, the oil spread on a smooth plane and on still and wavy water surfaces are calculated to examine the accuracy of simulating water waves by Navier-Stokes equations, the effect of damping zone on wave reflection and the precise structures of oil spread on waves.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.

Shallow Water Model On Cubed-Sphere By Multi-Moment Finite Volume Method

A global numerical model for shallow water flows on the cubed-sphere grid is proposed in this paper. The model is constructed by using the constrained interpolation profile/multi-moment finite volume method (CIP/MM FVM). Two kinds of moments, i.e. the point value (PV) and the volume-integrated average (VIA) are defined and independently updated in the present model by different numerical formulations. The Lax-Friedrichs upwind splitting is used to update the PV moment in terms of a derivative Riemann problem, and a finite volume formulation derived by integrating the governing equations over each mesh element is used to predict the VIA moment. The cubed-sphere grid is applied to get around the polar singularity and to obtain uniform grid spacing for a spherical geometry. Highly localized reconstruction in CIP/MM FVM is well suited for the cubed-sphere grid, especially in dealing with the discontinuity in the coordinates between different patches. The mass conservation is completely achieved over the whole globe. The numerical model has been verified by Williamson's standard test set for shallow water equation model on sphere. The results reveal that the present model is competitive to most existing ones. (C) 2008 Elsevier Inc. All rights reserved.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

Application of the characteristic CIP method to a shallow water model on the sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the efficiency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.

A model study on carbon cycle and phytoplankton dynamical processes in the Bohai Sea

The carbon cycle of lower trophic level in the Bohai Sea is studied with a three-dimension-al biological and physical coupled model. The influences of the processes (including horizontal advection,river nutrient load, active transport etc. ) on the phytoplankton biomass and its evolution are estimated.The Bohai Sea is a weak sink of the CO2 in the atmosphere. During the cycle, 13.7% of the gross pro-duction of the phytoplankton enter the higher trophic level and 76.8 % of it are consumed by the respira-tion itself. The nutrient reproduction comes mainly from the internal biogeochemical loop and the rem-ineralization is an important mechanism of the nutrient transfer from organic form to inorganic. Horizon-tal advection decreases the total biomass and the eutrophication in some sea areas. Change in the nutrientload of a river can only adjust the local system near its estuary. Controlling the input of the nutrient,which limits the alga growth, can be very useful in lessening the phytoplankton biomass.

Geophysical signatures associated with fluid flow and gas hydrate occurrence in a tectonically quiescent sequence, Qiongdongnan Basin, South China Sea

Interpretation of high-resolution two-dimensional (2D) and three-dimensional (3D) seismic data collected in the Qiongdongnan Basin, South China Sea reveals the presence of polygonal faults, pockmarks, gas chimneys and slope failure in strata of Pliocene and younger age. The gas chimneys are characterized by low-amplitude reflections, acoustic turbidity and low P-wave velocity indicating fluid expulsion pathways. Coherence time slices show that the polygonal faults are restricted to sediments with moderate-amplitude, continuous reflections. Gas hydrates are identified in seismic data by the presence of bottom simulating reflectors (BSRs), which have high amplitude, reverse polarity and are subparallel to seafloor. Mud diapirism and mounded structures have variable geometry and a great diversity regarding the origin of the fluid and the parent beds. The gas chimneys, mud diapirism, polygonal faults and a seismic facies-change facilitate the upward migration of thermogenic fluids from underlying sediments. Fluids can be temporarily trapped below the gas hydrate stability zone, but fluid advection may cause gas hydrate dissociation and affect the thickness of gas hydrate zone. The fluid accumulation leads to the generation of excess pore fluids that release along faults, forming pockmarks and mud volcanoes on the seafloor. These features are indicators of fluid flow in a tectonically-quiescent sequence, Qiongdongnan Basin. Geofluids (2010) 10, 351-368.

Characteristics and flux of settling particulate matter in neritic waters: The southern Yellow Sea and the East China Sea

The typology and flux of settling particulate matter (SPM) were investigated based on sediment trap sampling at six typical stations in the Yellow Sea and the East China Sea. The settling particulate matter in the neritic seas was sorted into three categories, lithogenic particles, living organisms, and particle aggregates. The mass of individual organisms is an important portion of particulate matter in the neritic waters. The aggregates contain six types, mucus aggregates, fecal pellets, diatom aggregates, silicoflagellate aggregates, tintinnids, and miscellaneous aggregates, of which the silicoflagellate aggregates and tintinnids are the most abundant in the Yellow Sea and the East China Sea. High particle fluxes, such as 215 to 874 g m(-2). day(-1) SPM in the bottom layer, were found at three stations where the water was well mixed, and the maximum flux was detected in the boundary area between the Yellow Sea and the East China Sea, where a wide nepheloid layer was present. Hence, particle flux in neritic waters can be easily shifted by water turbulence. The net vertical flux (123 to 961 mg C day(-1)), the contribution of lateral advection to resuspension flux (5 to 76%), and the particulate organic carbon export ratio (18 to 60%) were estimated for the other three stations where the water was stratified. The highest values were all found in the upwelling area off the Zhejiang coast, suggesting that the area of high productivity provides a high net vertical flux of SPM. (C) 2010 Elsevier Ltd. All rights reserved.

Effect of low-frequency Rossby wave on thermal structure of the upper southwestern tropical Indian Ocean

We investigate the influence of low-frequency Rossby waves on the thermal structure of the upper southwestern tropical Indian Ocean (SWTIO) using Argo profiles, satellite altimetric data, sea surface temperature, wind field data and the theory of linear vertical normal mode decomposition. Our results show that the SWTIO is generally dominated by the first baroclinic mode motion. As strong downwelling Rossby waves reach the SWTIO, the contribution of the second baroclinic mode motion in this region can be increased mainly because of the reduction in the vertical stratification of the upper layer above thermocline, and the enhancement in the vertical stratification of the lower layer under thermocline also contributes to it. The vertical displacement of each isothermal is enlarged and the thermal structure of the upper level is modulated, which is indicative of strong vertical mixing. However, the cold Rossby waves increase the vertical stratification of the upper level, restricting the variability related to the second baroclinic mode. On the other hand, during decaying phase of warm Rossby waves, Ekman upwelling and advection processes associated with the surface cyclonic wind circulation can restrain the downwelling processes, carrying the relatively colder water to the near-surface, which results in an out-of-phase phenomenon between sea surface temperature anomaly (SSTA) and sea surface height anomaly (SSHA) in the SWTIO.

Modeling vertical heat transport in the wave affected surface layer of the ocean

In considering the vertical heat transport problems in the upper ocean, the flat upper boundary approximation for the free surface and the horizontal homogenous hypothesis are usually applied. However, due to the existence of the wave motion, the application of this approximation may result in some errors to the solar irradiation since it decays quickly in respect to the actual thickness of the water layer below the surface; on the other hand, due to the fluctuation of the water layer depth, it is improper to neglect the effects of the horizontal advection and turbulent diffusion since they also contribute to the vertical heat transport. A new model is constructed in this study to reflect these effects. The corresponding numerical simulations show that the wave motion may remarkably accelerate the vertical heat transferring process and the variation of the temperature in the wave affected layer appears in an oscillating manner.

An N-shape thermal front in the western South Yellow Sea in winter

An N-shape thermal front in the western South Yellow Sea (YS) in winter was detected using Advanced Very High Resolution Radiation (AVHRR) Sea Surface Temperature data and in-situ observations with a merged front-detecting method. The front, which exists from late October through early March, consists of western and eastern wings extending roughly along the northeast-southwest isobaths with a southeastward middle segment across the 20-50 m isobaths. There are north and south inflexions connecting the middle segment with the western and eastern wings, respectively. The middle segment gradually moves southwestward from November through February with its length increasing from 62 km to 107 km and the southern inflexion moving from 36.2A degrees N to 35.3A degrees N. A cold tongue is found to coexist with the N-shape front, and is carried by the coastal jet penetrating southward from the tip of the Shandong Peninsula into the western South YS as revealed by a numerical simulation. After departing from the coast, the jet flows as an anti-cyclonic recirculation below 10 m depth, trapping warmer water originally carried by the compensating Yellow Sea Warm Current (YSWC). A northwestward flowing branch of the YSWC is also found on the lowest level south of the front. The N-shape front initially forms between the cold tongue and warm water involved in the subsurface anti-cyclonical recirculation and extends upwards to the surface through vertical advection and mixing. Correlation analyses reveal that northerly and easterly winds tend to be favorable to the formation and extension of the N-shape front probably through strengthening of the coastal jet and shifting the YSWC pathway eastward, respectively.

STUDY ON THE FORMATIONAL MECHANISM OF THE NORTHERN YELLOW (HUANGHAI) SEA COLD WATER MASS .2. DISCUSSION ON THE SOLUTION OF THE MODEL

The theoretical solution of the model of the Northern Yellow (Huanghai) Sea Cold Water Mass (NYSCWM) reveals that the NYSCWM is mainly formed through the continuous temperature increase of the overwintered water body above the Northern Yellow Sea Depression (NYSD) after spring when heat is continuously conducted from the sea surface to the deeper layer. In the NYSCWM's growing period, (June-July), nonlinear vertical convection and advection effects continuously increase, and are gradually balanced by the heat diffusion effect as the temperature increases from the surface to the bottom, which leads to the formation of an intensive thermocline and lateral front. Meanwhile, the three-dimensional circulation correspondingly occurs. In the NYSCWM's entire growing period, the horizontal circulation is always in the cyclonic motion, while the vertical circulation passes through a transition from a period with the cold centre as downwelling to a period with the cold centre as upwelling.