Vertical mixing in the bottom mixed layer near the continental shelf break in the East China Sea

The strength of the vertical mixing in the bottom mixed layer near the continental shelf break in the East China Sea was directly measured with the Micro-Scale Profiler (MSP). It has been shown that there is no significant statistical relation between the turbulent energy dissipation and the degree of the stratificationN ². It seems that the vigorous turbulence occurs not constantly but intermittently in the bottom mixed layer so that a large variation of is found depending on the time. In contrast to , the coefficient of the vertical eddy diffusivityK _z is mostly determined byN such thatK _z is large in the bottom mixed layer and small in the thermocline. Large value ofK _z in the bottom mixed layer is also found in the time series ofK _z estimated in terms of Richardson number calculated from the data obtained with electromagnetic current meters. The value ofK _z more than 10 cm²s^–1 frequently occur in the layer of 20–25 m thick just above the bottom.     

Estimation of phytoplankton pigments,total suspended matter and the concentration of its components by Secchi depth in Guinea's coastal zone

Relationships between the surface concentrations of phytoplankton pigments (C _ph), total suspended matter (C _sm), particulate organic carbon (C _poc), and total suspended phosphorus (C _sp), on the one hand, and the relative water transparency determined through the Secchi disc depth (z _d), on the other, are analysed using the data compiled in the Guinean coastal waters (Tabunsu and Tonkima river estuaries) during November–December 1990. The functions ofC _ph,C _sm=f(z _d) are matched up with the experimentally derived data, as well as with the model bio-optical state of seawater. The general regression equation has been calculated using the data characteristic of various types of water.Translated by Vladimir A. Puchkin.     

Vertical distribution of the hydrogen sulphide sources in the Black Sea

Data on the diffusion coefficientK _zand the concentration of H₂S in the Black Sea are used to compute the depth distribution of the vertical flux, and the intensity of the H₂S sources and sinks. On average, the total production of H₂S in the Black Sea reaches 37×10⁶ t/year. The main bulk of H₂S is produced not at the bottom, but in the layer of 450–1300 m. Destruction of H₂S prevails above the 400 m layer. Dissolved oxygen penetrating the H₂S zone can oxidize only half of the hydrogen sulphide produced in the sea.Translated by Mikhail M. Trufanov.     

Vertical transport in a stratified oceanic medium

On the basis of the analysis of the data of hydrological observations, we construct an empirical dependence of the vertical heat fluxQ on the temperature gradientT _z . The integral advective-diffusive transfer is taken into account by the method of Kolesnikov. We obtain a generalized dependence of the coefficient of vertical thermal diffusivityK _z onT _z and note that the well-known Rossby-Montgomery relation and Kolesnikov formula are its special cases. We analyze the possibility of application of this dependence to linear and nonlinear problems of thermal conductivity and diffusion in the ocean. Translated by Peter V. Malyshev and Dmitry V. Malyshev     

Mixed and mixing layer depths simulated by an OGCM

The global distributions of the mixed layer depth h _D , representing the depth of uniform density, and the mixing layer depth h _K , representing the depth of active turbulent mixing, were simulated using an ocean general circulation model (OGCM), and compared with each other, as well as with the mixed layer depth from the climatological data h _D *. The comparison between h _D and h _D * suggested that the threshold density difference Δ σ _θ should decrease from 0.09 kg m⁻³ to 0.02 kg m⁻³ with increasing latitude for consistent comparison between two mixed layer depths, due to the different nature of density profiles depending on latitude. The comparison between h _D and h _K revealed that h _K is deeper than h _D in the region where strong subsurface shear is present, such as the equatorial ocean and the region of the western boundary current. On the other hand, h _K is shallower than h _D in the high latitude ocean during convective cooling and early restratification.     

Observation and measurement of the bottom boundary layer flow in the prebreaking zone of shoaling waves

In this paper, the characteristics of the bottom boundary layer flow induced by nonlinear, asymmetric shoaling waves, propagating over a smooth bed of 1/15 uniform slope, is experimentally investigated. Flow visualization technique with thin-layered fluorescent dye was first used to observe the variation of the flow structure, and a laser Doppler velocimeter was then employed to measure the horizontal velocity, U.The bottom boundary layer flow is found to be laminar except within a small region near the breaking point. The vertical distribution of the phase-averaged velocity U at each phase is non-uniform, which is directly affected by the mean velocity, . The magnitude of increases from zero at the bottom to a local positive maximum at about z/δ2.02.5 (where z is the height above the sloping bottom and δ is the Stokes layer thickness), then decreases gradually to zero at z/δ6.07.0 approximately, and finally becomes negative as z/δ increases further. Moreover, as waves propagate towards shallower water, the rate of increase in the maximum onshore oscillating velocity component is greater than that of the offshore counterpart except near the breaking point. The free stream velocities in the profiles of the maximum onshore and offshore oscillating velocity components, and are found to appear at z/δ≥6.0. This implies that, if the Stokes layer thickness is used as a length scale, the non-dimensionalized boundary layer thickness remains constant in the pre-breaking zone. Although is greater than and the asymmetry of the maximum free stream velocities (i.e. ) increases with decrease of water depth, a universal similar profile can be established by plotting z/δ versus ( ) or ( ). The final non-dimensional profile is symmetric and unique for the distributions of the maximum onshore and offshore oscillating velocity components within the bottom boundary layer, which are induced by nonlinear, asymmetric shoaling waves crossing the pre-breaking zone.     

Formation of the vertical salinity profile in the Black Sea

The vertical distribution of salts brought by the Bosphorus undercurrent is numerically evaluated. By multiplying the average vertical salinity gradient by the diffusion coefficient,K _z , and the cross-section of the sea at the appropriate depth, we can determine the total vertical salt flux,Q(z). The derivative ofQ with respect toz depicts the salt source intensity distribution over depth. The highest intensity, Q/z, matches the 200 m depth level, i.e. the shelf edge. Below 1500 m, Q/z equals merely 0.1% of the value observed at a depth of 200 m. Above 37 m, salts are noted to sink, which corresponds to their outflow with the Bosphorus current. The distribution of Q/z and the respective values of mineral phosphorus and hydrogen sulphide are matched up.Translated by Vladimir A. Puchkin.     

边界层高度和动量粗糙度的不确定性对台风“天鸽”(1713)模拟的影响

基于中国气象局区域台风数值预报系统(CMA-TYM),通过一系列敏感性试验,分析研究了台风“天鸽”(1713)生命过程中不同时间阶段的移动路径及强度的数值模拟对模式参数化方案中边界层高度(h)和动量粗糙度(z₀)的敏感性。试验结果表明,使用不同参数化方案计算的h在“天鸽”(1713)初期的热带风暴阶段对热带气旋的移动路径有较明显的影响,在台风成熟后对热带气旋的移动路径影响不显著 ;台风中心附近最大10 m风速对h的变化不敏感,而最低海平面气压对h的变化却非常敏感。同时发现在台风发展初期阶段,边界层过薄或过厚都不利于台风强度的发展加强。这表明边界层高度h在热带气旋数值模拟和预报中是非常重要的,尤其是在台风发展的初期阶段。动量粗糙度z0的变化对台风“天鸽”的影响主要体现在台风增强阶段,台风中心附近10 m风速最大值在台风增强阶段对z₀敏感,尤其是在台风发展的初期阶段     

Effect of bottom slope in northeastern North Pacific on deep-water upwelling and overturning circulation

Hydrographic data show that the meridional deep current at 47°N is weak and southward in northeastern North Pacific; the strong northward current expected for an upwelling in a flat-bottom ocean is absent. This may imply that the eastward-rising bottom slope in the Northeast Pacific Basin contributes to the overturning circulation. After analysis of observational data, we examine the bottom-slope effect using models in which deep water enters the lower deep layer, upwells to the upper deep layer, and exits laterally. The analytical model is based on geostrophic hydrostatic balance, Sverdrup relation, and vertical advection–diffusion balance of density, and incorporates a small bottom slope and an eastward-increasing upwelling. Due to the sloping bottom, current in the lower deep layer intensifies bottomward, and the intensification is weaker for larger vertical eddy diffusivity (K _V), weaker stratification, and smaller eastward increase in upwelling. Varying the value of K _V changes the vertical structure and direction of the current; the current is more barotropic and flows further eastward as K _V increases. The eastward current is reproduced with the numerical model that incorporates the realistic bottom-slope gradient and includes boundary currents. The interior current flows eastward primarily, runs up the bottom slope, and produces an upwelling. The eastward current has a realistic volume transport that is similar to the net inflow, unlike the large northward current for a flat bottom. The upwelling water in the upper deep layer flows southward and then westward in the southern region, although it may partly upwell further into the intermediate layer.     

Eddy viscosities and velocities in combined wave-current flows

The eddy viscosities for the steady and the periodic components of combined wave-current flows have been studied quantitatively from the presently available experimental data. It has been found that inside the boundary interaction layer [z < δ] the eddy viscosity ε_c for the steady flow is increased in the presence of waves while outside the boundary interaction layer [z >δ] it is affected little by the wave motion, and that the eddy viscosity ε_w for the wave motion in the boundary layer is independent of the current strength U^*. On the other hand, a new eddy viscosity model is presented to give a good prediction of the velocity distributions of the waves and currents in comparison with experimental data.     

Chemical environment for red tides due toChattonella antiqua

Environmental parameters that affect the growth ofChattonella antiqua were monitored throughout the outbreak period of this species around the Ie-shima Islands, the Seto Inland Sea, in the summer of 1987 (20 July–13 August). Averaged cell concentration ofC. antiqua over the water column (21 m) was below 10 cells· ml^–1 on 20 July, gradually increased to reach the maximum of 250 cells·ml^–1 on 7 August, and then rapidly decreased to the value of 30 cells·ml^–1 on 13 August.Thermal stratifications were prominent from 20 July to 3 August and were destroyed after 4 August. Temperature and salinity were optimum for the growth ofC. antiqua throughout the survey period.At the bloom initiation period (20–21 July), concentrations of N- and P-nutrients (S _N andS _P ) were high throughout the water column. From 22 July to 3 August, whenC. antiqua increased its populations,S _N andS _P at the depth of 0–5m were low but those at the depth of 10–20m kept a high value. After 4 August,S _n andS _P at the depth of 10–20m decreased rapidly due to wind mixing coupled with the nutrient uptake byC. antiqua. When the populations ofC. antiqua reached the maximum (7–9 August), N-nutrients were depleted throughout the water column but P-nutrients were not. Concentrations of vitamin B₁₂ were almost in the same range as those of the previous years and were optimum for the growth ofC. antiqua.GP- value (growth potential of the seawater with respect to nitrogen and phos-phorus) was higher than 0.6 even at the surface layer (0–5 m) at the bloom-initiation period. During the bloom development period (22 July–3 August), GP at the surface layer (0–5m) was low (<0.2), but GP at the depth of 10–20m kept a rather high value (>0.4).In situ growth rates ofC. antiqua at the depth of 0 and 5m estimated from bottle experiments coincided well with the values expected from GP. A high value of GP at the surface layer in the initiation period and a shallow GP-cline in the development period, combined with the ability of diurnal vertical migration seemed to be at least one reason that natural populations ofC. antiqua grew at a rather high rate and formed red tides in the summer of 1987.     

Microstructure measurements in natural waters: Methodology and applications

Modern approaches to microstructure data processing, including wavelet denoising, are discussed. The wavelet procedure is applied to small-scale shear signals before estimating the dissipation rate ε and to the temperature/density profiles used to calculate Thorpe scales. Microstructure data obtained on the Mediterranean shelf of Catalonia are used to illustrate various approaches to the Thorpe displacement calculations. It is suggested that the Weibull probability function is an appropriate model for the Thorpe scale distribution. Microstructure measurements from the upper layer of the Boadella reservoir (Catalonia, Spain) support this finding.A new analytical approximation for the 1D Panchev–Kesich spectrum is deduced and the results of ε computation are compared with spectral fitting by the widely used Nasmyth spectrum. Applying the Kraichnan spectral model to compute ε from temperature spectra in the convective-viscous sub-range is examined as an alternative to the Batchelor spectrum. Microstructure measurements taken in Lake Banyoles (Catalonia, Spain) and in the North Atlantic were used for spectral calculations.Statistical analysis of eddy K_b and thermal K_θ diffusivities measured on a shallow shelf of the Black Sea shows the importance of process-orientated domain averaging of the diffusivities in obtaining good correspondence between K_b and K_θ in active turbulent regions. In weakly turbulent, stratified interior layers, the averaged K_b and K_θ differ significantly, which may point to the inapplicability of isotropic formulae used for ε and temperature dissipation χ_θ estimates, as well as to a dependence of the mixing efficiency γ on the Richardson number or in some cases on regions of fossil turbulence.     

Kinetics of nitrogen- or phosphorus-limited growth and effects of growth conditions on nutrient uptake inChattonella antiqua

Chattonella antiqua was grown in a nitrogen- or phosphorus-limited semicontinuous culture system. Using the cells in steady growth state, the relationship between growth rate and cell quota and effects of growth conditions on nitrate, ammonium and phosphate uptake were examined. Under nitrogen-limited conditions, growth rate as a function of nitrogen cell quota followed the empirical Droop equation and the uptake of nitrate and ammonium was not significantly affected by growth rate. Similarly, under phosphorus-limited conditions, the growth rate as a function of phosphorus cell quota also followed the Droop equation and phosphate uptake was not significantly affected by growth rate.Combining the results obtained in the present study with those from previous studies on nutrient uptake, half saturation constants for growth (K _g ) were calculated for nitrate, ammonium and phosphate. Comparisons ofK _g with nutrient concentrations in the Seto Inland Sea in summer, where red tides ofC. antiqua often occur, suggest that phosphate is one of the controlling factors for the population ofC. antiqua.     

A new expression for the production rate of sea water droplets on the sea surface

A new set of empirical formulas for the production rate and the number concentration of sea-water droplets on the sea surface are proposed, synthesizing past observation data of sea-salt particles in the sea and water droplets in wind-wave tanks. A new levelz _c is introduced as the effective wind-sea surface where seawater droplets are produced. The new formulas are expressed in linear functions in logarithmic scales ofu*²/v _p , a parameter to describe overall conditions of airsea boundary processes, whereu _* is the friction velocity of air,v the kinematic viscosity of air and _p the peak angular frequency of wind-wave part of wave spectra. A model of coexistence of spray droplets and suspended particles near the sea surface is proposed. As for the independent parameter, a comparison between the uses ofu*²/v _p and ofu _* ³ which was the traditional way of parameterization excluding wave measure, shows that the advantage of usingu*²/v _p is statistically significant with a confidence limit 89% in F-test.     

Circulation and mixing at the New England shelfbreak front: Results of purposeful tracer experiments

We present the results of six dye tracer experiments that measured the mixing and circulation at the shelfbreak front on the New England Shelf. The last three were conducted during the New England Shelfbreak Productivity Experiment (NESPEX) with concurrent isopycnal float deployments. The results are consistent with the Chapman and Lentz [Chapman, D.C., and Lentz, S.J. (1994). Trapping of a coastal density front by the bottom boundary layer. Journal of Physical Oceanography, 24, 1465–1479.] model prediction of the separation and upwelling along the shelfbreak front of bottom boundary layer (BBL) water forced by an Ekman buoyancy flux, but show considerable variability. Cross-shelf velocities at the detachment point are 2–3 × 10⁻² m/s. But seaward, over the slope region, dye tagged water was sheared from the main patch into small filaments that upwelled along the front with cross-shelf speeds up to 0.1 m/s. Cross-shelf diffusion was of order 10 m²/s in the mixed bottom layer and 1 m²/s in the interior along the front. Within the stratified front, the mean vertical diffusivity was K_z  4 × 10⁻⁶ m²/s. The dispersion of shelfwater in the slope region is effected by turbulent flow with advective speeds exceeding the small scale diffusive mixing. The mean flux of the detached BBL water is sufficient to account for the net loss of shelf water during its transit from Cape Cod to Cape Hatteras.     

Transverse mixing and surface heat exchange in the Waikato River: A comparison of two models

Abstract

Two mathematical models of different complexity were used to study transverse dispersion and surface heat transfer in the lower Waikato River. A simple analytical streamtube model (HPLUME) gave adequate temperature predictions in a reach where the channel was fairly regular but performed poorly where there were extensive shallows. In the latter reach, a two‐dimensional numerical model (SYSTEM21) gave good temperature and flow predictions once properly calibrated. Model calibration proved to be difficult in the Waikato River because the natural river temperature varied significantly along the channel. A search method was developed to estimate both the transverse dispersion and surface heat exchange coefficients from measured plume temperature profiles based on the observation that transverse variations in natural temperature were small. This method was used to calibrate SYSTEM21 in two separate reaches. Coefficient estimates were sensitive to measurement errors and slight departures in homogeneity of natural temperature and it would be desirable to corroborate the estimates of E_z using dye tests. In the upper reach, the average transverse dispersion coefficient was E_z/hu? = 1 which is high but within the range of published values. A sharp bend and buoyant spreading contributed to the high E_z value. E_z increased with river flow because both h and u? increase with flow. In the lower reach, E_z/hu? = 0.1 which is lower than expected but islands may have affected the model calibration. The surface heat exchange coefficient averaged K = 84–167 W m^?2 s^?1 °C^?1 which falls within the range of published values. No significant relationship could be found between surface heat exchange coefficient and meteorological variables.     

Prediction of wave pressures on smooth impermeable seawalls

Simple prediction methods are proposed to estimate the wave induced pressures on smooth impermeable seawalls. Based on the physics of the wave structure interaction, the sloped seawall is divided into a total of five zones (zones 1, 2 and 3 during run-up (corresponding pressures are called as positive pressures) and zones 4 and 5 during run-down (corresponding pressures are called negative pressures)) (Fig. 1). Zone 1 (0<z<d−H_i/2), where the wave pressure is governed by the partial reflection and phase shift; Zone 2 (d−H_i/2<z<d), where the effect of wave breaking and turbulence is significant; Zone 3 (d<z<Run-up height), where the pressure is induced by the run-up water; Zone 4 (Run-down<z<d), where the wave pressure is caused by the run-down effect and Zone 5 (0<z<d-Run down), where the negative wave pressures are due to partial reflection and phase shift effects. Here d is the water depth at the toe of the seawall, H_i is the incident wave height and z is the vertical elevation with toe of the seawall as origin and z is positive upward. For wave pressure prediction in zones 1 and 5, the empirical formula proposed by Ahrens et al. (1993) to estimate the wave reflection and Sutherland and Donoghue's recommendations (1998) for the estimation of phase shift of the waves caused by the sloped structures are used. Multiple regression analysis is carried out on the measured pressure data and empirical formulas are proposed for zones 2, 3 and 4. The recommendations of Van der Meer and Breteler (1990) and Schüttrumpf et al. (1994) for the prediction of wave run-down are used for pressure prediction at zone 4. Comparison of the proposed prediction formulas with the experimental results reveal that the prediction methods are good enough for practical purposes. The present study also shows a strong relation between wave reflection, wave run-up, wave run-down and phase shift of waves on wave pressures on the seawalls.     

Computations of self-propulsion free to sink and trim and of motions in head waves of the KRISO Container Ship (KCS) model

Two computations of the KCS model with motions are presented. Self-propulsion in model scale free to sink and trim are studied with the rotating discretized propeller from the Hamburg Model Basin (HSVA) at Fr = 0.26. This case is particularly complex to simulate due to the close proximity of the propeller to the rudder. The second case involves pitch and heave in regular head waves. Computations were performed with CFDShip-Iowa version 4.5, a RANS/DES CFD code designed for ship hydrodynamics. The self-propulsion computations were carried out following the procedure described in Carrica et al. [1], in which a speed controller is used to find the propeller rotational speed that results in the specified ship velocity. The rate of revolutions n, sinkage, trim, thrust and torque coefficients K_T, K_Q and resistance coefficient CT_(SP) are thus obtained. Comparisons between CFD and EFD show that the rate of revolutions n, thrust and torque coefficients K_T and K_Q have higher prediction accuracies than sinkage and trim. For the simulation of pitch and heave in head waves, the geometry includes KCS hull and rudder under three conditions with two Froude numbers and three wave length and amplitude combinations. 0th and 1st harmonic amplitudes and 1st harmonic phase are computed for total resistance coefficient C_T, heave motion z and pitch angle θ. Comparisons between CFD and EFD show that pitch and heave are much better predicted than the resistance. In both cases comparisons with simulations by other authors presented at the G2010 CFD Workshop [2] using different CFD methodologies are included.     

The hydrophobic fraction of organic matter in the Krka River Estuary

Surface active substances have been studied in the Krka River Estuary using an a.c. polarography method and various fractionation procedures. Total surface substances (SA_T) were separated into (particulate) heterodispersed and dissolved fractions by filtration. The dissolved fraction was separated into hydrophobic basic and neutral components, hydrophobic acid and hydrophilic components by sorption on to XAD-8 resin.The content and the composition of SA_T varied seasonally, depending on the biological activities taking place in the Krka River Estuary. The accumulation of SA_T, with hydrophobic properties, occurred at the fresh/saline water interface along the vertical profile. Daily changes in the form of the surface active substances along the vertical profile were observed. The formation of the heterodispersed fraction, especially observable in the upper freshwater layer and at the halocline, resulted in the change of physico-chemical properties of surface active substances. Organic matter appeared mainly in a dissolved form in the deeper saline water layer, showing negligible oscillations.The electrochemical characterization of the fractionated SA_T indicated that strong hydrophobicity of the surface active substances under natural pH conditions could be attributed to the significant contribution of the heterodispersed fraction of the organic material, which is the most reactive part of the surface active substances in the Krka River Estuary.The dissolved fraction of the surface active substances, prevailing in the deeper saline water layer, exhibited either hydrophilic or conditionally hydrophobic properties, depending on pH and ionic strength of the medium.     

Observation of wave-enhanced turbulence in the near-surface layer of the ocean during TOGA COARE

Dissipation rate statistics in the near-surface layer of the ocean were obtained during the month-long COARE Enhanced Monitoring cruise with a microstructure sensor system mounted on the bow of the research vessel. The vibration contamination was cancelled with the Wiener filter. The experimental technique provides an effective separation between surface waves and turbulence, using the difference in spatial scales of the energy-containing surface waves and small-scale turbulence. The data are interpreted in the coordinate system fixed to the ocean surface. Under moderate and high wind-speed conditions, we observed the average dissipation rate of the turbulent kinetic energy in the upper few meters of the ocean to be 3–20 times larger than the logarithmic layer prediction. The Craig and Banner (J. Phys. Oceanogr. 24 (1994) 2546) model of wave-enhanced turbulence with the surface roughness length from the water side z₀ parameterized according to the Terray et al. (J. Phys. Oceanogr. 26 (1996) 792) formula z₀=cH_s provides a reasonable fit to the experimental dissipation profile, where z is the depth (defined here as the distance to the ocean surface), c≈0.6, and H_s is the significant wave height. In the wave-stirred layer, however, the average dissipation profile deviates from the model (supposedly because of extensive removing of the bubble-disturbed areas close to the ocean surface). Though the scatter of individual experimental dissipation rates (10-min averages) is significant, their statistics are consistent with the Kolmogorov's concept of intermittent turbulence and with previous studies of turbulence in the upper ocean mixed layer.