THE VARIATION AND MOVEMENT OF A KIND OF EDDIES IN SHELF SLOPE WATER

Considering the characteristics of the time and space scales of the eddies we established a quasi-static and quasi-geostrophic model to describe their variation and movement in shelf slope water. The analytical solution revealed the main properties of the variation: slow expansion and fast stagnation processes and the law of the eddy motion affected under the background field. All theoretical results are proved by satellite image measurements.

Anti-cyclonic eddies northwest of Luzon in summer-fall observed by satellite altimeters

Anti-cyclonic eddies northwest of Luzon of the Philippines in summer-fall are identified in the merged data products of satellite altimeters of Topex/Poseidon, Jason-1 and European Research Satellites. The generation and propagation of the anti-cyclonic eddies, which are confirmed by satellite ocean color data, are found to be a seasonal phenomenon that is phase-locked to the onset of the southwesterly monsoon and the relaxation of the cyclonic wind curl in the northeastern South China Sea. The eddies originate from northwest of Luzon in summer, move across the northeastern South China Sea to reach the China continental slope in fall, and propagate southwestward along the continental slope in fall-winter, inducing shelfbreak current variations in the western South China Sea in fall-winter. The anti-cyclonic eddy discovered by Li et al. (1998) in the northern South China Sea is found to originate from northwest of Luzon and carry primarily the South China Sea waters. It does not appear to be an eddy shed from the Kuroshio in the Luzon Strait area as alluded by Li et al. (1998) and others.

Application of altimetry data assimilation on mesoscale eddies simulation

Mesoscale eddy plays an important role in the ocean circulation. In order to improve the simulation accuracy of the mesoscale eddies, a three-dimensional variation (3DVAR) data assimilation system called Ocean Variational Analysis System (OVALS) is coupled with a POM model to simulate the mesoscale eddies in the Northwest Pacific Ocean. In this system, the sea surface height anomaly (SSHA) data by satellite altimeters are assimilated and translated into pseudo temperature and salinity (T-S) profile data. Then, these profile data are taken as observation data to be assimilated again and produce the three-dimensional analysis T-S field. According to the characteristics of mesoscale eddy, the most appropriate assimilation parameters are set up and testified in this system. A ten years mesoscale eddies simulation and comparison experiment is made, which includes two schemes: assimilation and non-assimilation. The results of comparison between two schemes and the observation show that the simulation accuracy of the assimilation scheme is much better than that of non-assimilation, which verified that the altimetry data assimilation method can improve the simulation accuracy of the mesoscale dramatically and indicates that it is possible to use this system on the forecast of mesoscale eddies in the future.

Numerical Study Of Heat Transfer Mechanism In Turbulent Supercritical Co2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Multiscale Method for Turbulent-Flow Computations

A main method of predicting turbulent flows is to solve LES equations, which was called traditional LES method. The traditional LES method solves the motions of large eddies of size larger than filtering scale An while modeling unresolved scales less than Delta_n. Hughes et al argued that many shortcomings of the traditional LES approaches were associated with their inabilities to successfully differentiate between large and small scales. One may guess that a priori scale-separation would be better, because it can predict scale-interaction well compared with posteriori scale-separation. To this end, a multi-scale method was suggested to perform scale-separation computation. The primary contents of the multiscale method are l) A space average is used to differentiate scale. 2) The basic equations include the large scale equations and fluctuation equations. 3) The large-scale equations and fluctuation equations are coupled through turbulent stress terms. We use the multiscale equations of n=2, i.e., the large and small scale (LSS) equations, to simulate 3-D evolutions of a channel flow and a planar mixing layer flow Some interesting results are given.

Numerical Study of Heat Transfer Mechanism in Turbulent Supercritical CO2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Computation of turbulence-generated noise by large-eddy simulation

The hybrid method of large eddy simulation (LES) and the Lighthill analogy is being developed to compute the sound radiated from turbulent flows. The results obtained from the hybrid method are often contaminated by the absence of small scales in LES, since the energy level of sound is much smaller than that of turbulent flows. Previous researches investigate the effects of subgrid sacle (SGS) eddies on the frequency spectra of sound radiated by isotropic turbulence and suggest a SGS noise model to represent the SGS contributions to the frequency spectra. Their investigations are conducted in physical space and are unavoidably influenced by boundary conditions. In this paper, we propose to perform such calculations in Fourier space so that the effects of boundary conditions can be correctly treated. Posteriori tests are carried out to investigate the SGS contribution to the sound. The results obtained recover the -7/2 law within certain wave-number ranges, but under-estimate the amplitudes of the frequency spectra. The reason for the underestimation is also discussed.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

The Yangtze River dolphin or baiji (Lipotes vexillifer): population status and conservation issues in the Yangtze River, China

1. Baiji were sighted 17 times during three recent simultaneous multi-vessel surveys in the Yangtze River, China (November 4-10, 1997; December 4-9, 1998; October 31-November 5, 1999). There were 11 sightings in 1997 (consisting of 17 animals), five in 1998 (seven animals), and two in 1999 (four animals). It was concluded that 13 individuals Could be considered as a minimum number of the baiji currently in the Yangtze River. 2. An annual rate of population decrease was roughly estimated as 10%. From the body sizes observed, the proportions of old, adult and immature individuals were approximately estimated at 57, 26, and 17% respectively. 3. Baiji showed a significant attraction to confluences and sand bars with large eddies. The present distribution range of the baiji is less than 1400 km in length in the Yangtze main river. Distances between the two nearest groups of baiji appear to be increasing. 4. Two typical sightings are described, in which surfacing and movements of baiji were recorded. Baiji were often found swimming together with finless porpoise. In the surveys they occurred in the same group in 63% of occurrences. Interactions between baiji and finless porpoise are described and discussed. 5. Human activities are the main threats to the baiji. Illegal electrical fishing accounted for 40% of known mortalities during the 1990s. Engineering explosions for maintaining navigation channels have become another main cause of baiji deaths. The last hope of saving the species may be to translocate the remaining baiji into a semi-captive reserve. known as the 'Baiji Semi-natural Reserve'. Copyright (C) 2003 John Wiley Sons, Ltd.

A quasi-synoptic interpretation of water mass distribution and circulation in the western North Pacific II: Circulation

Using the data of conductivity-temperature-depth (CTD) intensive observations conducted during Oct.-Nov. 2005, this study provides the first three-dimension quasi-synoptic description of the circulation in the western North Pacific. Several novel phenomena are revealed, especially in the deep ocean where earlier observations were very sparse. During the observations, the North Equatorial Current (NEC) splits at about 12A degrees N near the sea surface. This bifurcation shifts northward with depth, reaching about 20A degrees N at 1 000 m, and then remains nearly unchanged to as deep as 2 000 m. The Luzon Undercurrent (LUC), emerging below the Kuroshio from about 21A degrees N, intensifies southward, with its upper boundary surfacing around 12A degrees N. From there, part of the LUC separates from the coast, while the rest continues southward to join the Mindanao Current (MC). The MC extends to 2 000 m near the coast, and appears to be closely related to the subsurface cyclonic eddies which overlap low-salinity water from the North Pacific. The Mindanao Undercurrent (MUC), carrying waters from the South Pacific, shifts eastward upon approaching the Mindanao coast and eventually becomes part of the eastward undercurrent between 10A degrees N and 12A degrees N at 130A degrees E. In the upper 2 000 dbar, the total westward transport across 130A degrees E between 7.5A degrees N and 18A degrees N reaches 65.4 Sv (1 Sv = 10(-6) m(3)s(-1)), the northward transport across 18A degrees N from Luzon coast to 130A degrees E is up to 35.0 Sv, and the southward transport across 7.5A degrees N from Mindanao coast to 130A degrees E is 27.9 Sv.

Eddy generation and evolution in the North Pacific Subtropical Countercurrent (NPSC) zone

The seasonal generation and evolution of eddies in the region of the North Pacific Subtropical Countercurrent remain poorly understood due to the scarcity of available data. We used TOPEX/POSEIDON altimetry data from 1992 to 2007 to study the eddy field in this zone. We found that velocity shear between this region and the neighboring North Equatorial Current contributes greatly to the eddy generation. Furthermore, the eddy kinetic energy level (EKE) shows an annual cycle, maximum in April/May and minimum in December/January. Analyses of the temporal and spatial distributions of the eddy field revealed clearly that the velocity shear closely related to baroclinic instability processes. The eddy field seems to be more zonal than meridional, and the energy containing length scale shows a surprising lag of 2-3 months in comparison with the 1-D and 2-D EKE level. A similar phenomenon is observed in individual eddies in this zone. The results show that in this eddy field band, the velocity shear may drive the EKE level change so that the eddy field takes another 2-3 months to grow and interact to reach a relatively stable state. This explains the seasonal evolution of identifiable eddies.

Variability of surface circulation in the South China Sea from satellite altimeter data

11-year satellite altimeter sea surface height (SSH) anomaly data from January 1993 to December 2003 are used to present the dominant spatial patterns and temporal variations of the South China Sea (SCS) surface circulation through Empirical Orthogonal Function (EOF) analysis. The first three EOF modes show the obvious seasonal variations of SSH in the SCS. EOF mode one is generally characterized by a basin-wide circulation. Mode two describes the double-cell basin scale circulation structure. The two cells were located off west of the Luzon Island and southeast of Vietnam, respectively. EOF mode three presents the mesoscale eddy structure in the western SCS, which develops into a strong cyclonic eddy rapidly from July to September. EOF mode one and mode three are also embedded with interannual signals, indicating that the SCS surface circulation variation is influenced by El Nino events prominently. The strong El Nino of 1997/98 obviously changed the SCS circulation structure. This study also shows that there existed a series of mesoscale eddies in the western SCS, and their temporal variation indicates intra-seasonal and interannual signals.

Vertical structure and evolution of the Luzon Warm Eddy

Eddies are frequently observed in the northeastern South China Sea (SCS). However, there have been few studies on vertical structure and temporal-spatial evolution of these eddies. We analyzed the seasonal Luzon Warm Eddy (LWE) based on Argo float data and the merged data products of satellite altimeters of Topex/Poseidon, Jason-1 and European Research Satellites. The analysis shows that the LWE extends vertically to more than 500 m water depth, with a higher temperature anomaly of 5A degrees C and lower salinity anomaly of 0.5 near the thermocline. The current speeds of the LWE are stronger in its uppermost 200 m, with a maximum speed of 0.6 m/s. Sometimes the LWE incorporates mixed waters from the Kuroshio Current and the SCS, and thus has higher thermohaline characteristics than local marine waters. Time series of eddy kinematic parameters show that the radii and shape of the LWE vary during propagation, and its eddy kinetic energy follows a normal distribution. In addition, we used the empirical orthogonal function (EOF) here to analyze seasonal characteristics of the LWE. The results suggest that the LWE generally forms in July, intensifies in August and September, separates from the coast of Luzon in October and propagates westward, and weakens in December and disappears in February. The LWE's westward migration is approximately along 19A degrees N latitude from northwest of Luzon to southeast of Hainan, with a mean speed of 6.6 cm/s.

冬季东中国海环流及中尺度涡旋数值模拟

该文在POM模式基础上,建立一个б坐标系下的三维斜压预报模式,考虑了海底地形、外来流、长江径流、海面风应力、海面热通量等多方面因素的影响,模拟结果较好地体现了冬季东中国海环流及涡旋的特征.在此数值模拟的基础上,分别在改变风场、以及黑潮流量的条件下进行数值实验,以考察它们对冬季东中国海环流和中尺度涡旋的影响.数值模拟和数值实验的结果表明,黑潮在流经东海时沿东海陆坡流动,在其两侧出现一些涡旋.冬季台湾暖流的表层水主要来自台湾东北的黑潮水,少部分来自台湾海峡;次表层,台湾东北的黑潮水和台湾海峡的入流对台湾暖流的贡献相当;深层,台湾暖流主要来源于台湾东北的黑潮水.冬季对马暖流的主要来源是黑潮分支和台湾暖流的分支,朝鲜沿岸流也是对马暖流的来源之一.

大洋中尺度涡旋与源区黑潮的相互作用研究

黑潮是西北太平洋一支高温、高盐、高流速的西边界流，它起源于菲律宾以东海域，沿吕宋岛、台湾东岸、东海大陆架外缘流动，并穿越吐噶喇海峡返回太平洋。其中吕宋岛至台湾岛这段西边界流命名为源区黑潮。源区黑潮在流经吕宋海峡以东时，由于失去岸界支撑会发生形变，并通过吕宋海峡与南海环流体系进行质量、动量和能量的交换，并影响着南海的环流特征。因此研究黑潮通过吕宋海峡入侵南海的关键动力过程和机制，对我们国家的军事、生产、环境等具有十分重大的意义，一直是国内外海洋学家关注的焦点。 本文利用高度计资料和涡旋精度的高分辨率数值模式（OFES）输出结果分析了南海东北部和吕宋海峡东侧海域的上层环流特征，并独立建立1.5层浅水约化重力模式下的理想模型，分析了黑潮从吕宋海峡入侵南海的非线性特征和西传的大洋中尺度涡旋与黑潮在吕宋海峡处的相互作用。 对高度计资料和模式结果资料分析的结果表明季风是影响南海北部上层环流的主要因素，黑潮以各种方式入侵南海也是影响南海北部环流特征的重要组成部分。通过数值模拟，我们得到在不同的参数下，黑潮入侵南海会存在稳定态、分叉、周期解等非线性特征，对应着黑潮南海流套、分离流环、无入侵等多种流态。通过讨论大洋中尺度涡旋与黑潮在吕宋海峡处的相互作用，发现涡旋在西传遇到黑潮后基本局地耗散或随黑潮继续向北移动，但当黑潮较弱时，涡旋可能将大洋中的能量通过吕宋海峡传入南海。