Linear Approximations of the Second Turbulent Moments of the Atmospheric Convective Surface Layer in a Forced-Convection Sublayer

The Monin–Obukhov similarity theory for the convective surface layer distinguishes two limiting cases: a dynamic limit and a free-convection limit. The dynamic limit for the convective surface layer is defined as a flow with a logarithmic profile of wind and a zero buoyancy flux at the underlying surface. The free-convection limit is characterized by a zero wind speed and a positive buoyancy flux at the underlying surface. The limits of the generalized Monin–Obukhov similarity theory are able to describe the higher order turbulent moments. In this paper, it is assumed that the convective surface layer consists of two sublayers: the lower dynamic sublayer adjacent to the surface and the upper forced-convection sublayer. The turbulent moments can be approximated separately for each sublayer. Linear approximations are suggested for the turbulent moments of the vertical velocity and the potential temperature variance in the forced-convection sublayer. The first-order expansion terms of them correspond to the free-convection limits of the Monin–Obukhov theory under no-wind conditions. The second-order expansion terms describe profiles of the turbulent moments in under convective conditions with a moderate wind. A comparison between the proposed approximations and experimental data strongly suggests that the linear approximation is correct within a forced-convection range.     

卫星SAR图像的海洋大气边界层特征参数的反演

针对中国南海ERS-2／SAR三维对流涡旋图像,分别利用Young的协方差方法和惯性子域方法,反演了对流海洋大气边界层高度、海洋大气边界层垂直对流尺度速度、表面浮力通量、Obukhov长度、稳定性校正后的海面风速和海面风速稳定性校正因子。结果表明,2种反演方法的结果符合较好,在一定程度上说明利用协方差方法和惯性子域方法从卫星SAR图像反演海洋大气边界层参数的有效性。     

Preliminary results of measurements of atmospheric turbulence over the sea

We perform the experimental verification of the applicability of the theory of similarity to the wave boundary layer and the assessment of wave-induced perturbations of the air flow depending on various conditions of stratification of the atmosphere and the state of the sea. The measurements were carried out from a stationary platform located in the coastal part of the Black Sea. The experimental procedure is based on the simultaneous measurements of the profile and fluctuations of the wind speed at 5–6 levels in the 1.3–21-m layer, the elevations of the sea surface, the directions of waves and winds, and the mean gradients of temperature and humidity of air. The structure of the boundary layer in the region of measurements depends on the direction of the wind. For weak and moderate onshore winds (< 9 m/sec), the approximate balance is preserved between the production and dissipation of turbulent energy in the cases of unstable and neutral stratification. On the average, the estimates of friction velocity according to the profiles are higher than the dissipative estimates by 10% mainly due to the deficiency of dissipation near the surface. For the offshore wind, the structure of the boundary layer abruptly changes and is determined not by the local parameters but by strong turbulent eddies formed over the dry land. The intensity of low-frequency turbulent fluctuations and the gradient of wind velocity near the surface in the coastal zone are 1.5–2 times higher than for the open sea. __________ Translated from Morskoi Gidrofizicheskii Zhurnal, No. 3, pp. 42–61, May–June, 2007.     

Average velocity field of the air flow over the water surface in a laboratory modeling of storm and hurricane conditions in the ocean

Laboratory experiments on studying the structure of the turbulent air boundary layer over waves were carried out at the Wind-Wave Channel of the Institute of Applied Physics, Russian Academy of Sciences (IAP RAS), in conditions modeling the near-water boundary layer of the atmosphere under strong and hurricane winds and the equivalent wind velocities from 10 to 48 m/s at the standard height of 10 m. A modified technique of Particle Image Velocimetry (PIV) was used to obtain turbulent pulsation averaged velocity fields of the air flow over the water surface curved by a wave and average profiles of the wind velocity. The measurements showed that the logarithmic part of the velocity profile of the air flow in the channel was observed in the immediate vicinity from the water surface (at a distance of 30 mm) and could be detected only using remote methods (PIV). According to the measured velocity profiles, dependences of aerodynamic drag factors of the water surface on the wind velocity at a height of 10 m were retrieved; they were compared with results of contact measurements carried out earlier on the same setup. It is shown that they agree with an accuracy of up to 20%; at moderate and strong wind velocities the coincidence falls within the experimental accuracy.     

Formation of a circulation in the quasiisopycnic model of the Black Sea taking into account the stochastic nature of the wind stress

We discuss the results of numerical experiments carried out within the framework of a new version of the quasiisopycnic model of the Black Sea [1, 2]. We take into account the stochasticity of the wind stress, vertical shift of the flow velocity in the upper quasihomogeneous layer, and diapycnic mixing. We generalize the procedure of “convective adjustment,” which enables us to consider the propagation and transformation of waters of the Sea of Marmara. It is shown that the consideration of these factors enables one to determine the structure and seasonal variability of hydrothermodynamic fields in the sea more correctly. Translated by Peter V. Malyshev and Dmitry V. Malyshev     

The wind field over wind-waves

The velocity fluctuations of wind over wind-waves in a wind tunnel are measured with a X-type hot-wire anemometer at some heights over the water surface.The observed vertical profiles of the wave-induced velocity fluctuations and the wave-induced Reynolds stress at the wave spectral peak frequency are different from those expected from the inviscid quasi-laminar model;i.e., the observed vertical profiles of the power spectral density of the wave-induced horizontal or vertical velocity fluctuations of wind have the minimum value at the height much heigher than the critical layer, and the value of the wave-induced Reynolds stress is negative at several heights over the water surface. From the comparison between the experimental results and the numerical solutions of a linear model of the turbulent shear flow over the wavy boundary, it is shown that the discrepancy described above can be attributed to the atmospheric turbulence.     

Model analysis of the currents and wind waves in the Vistula Lagoon of the Baltic Sea

This paper considers the results of the numerical simulation of the steady currents and waves field structures in the Vistula Lagoon under constant wind forcing. The currents?? structure is essentially 3-dimensional. The direction of the near-shore drift is determined by the wind and shore orientation, but the currents involve two layers in the deeper part of the basin: the upper layer downwind current and the upwind compensative current in the underlying layer. The wind waves depend upon the wind speed, fetch, and depth until the wind speed is less than 6?C8 m/s. As the wind increases, the fetch??s dependence disappears. The wind height depends only on the wind??s magnitude, and, under stormy conditions (when the wind??s speed exceeds 15 m/s), the waves are limited by the water depth alone.     

Suspension of sediment particles over a ripple due to turbulent convection under unsteady flow conditions

We analyzed the motions of small sediment particles over a sinusoidal ripple due to an unsteady turbulent boundary layer flow using Large Eddy Simulation. The motions of sediment particles are described in terms of the Lagrangian framework as it is helpful in studying the structure of sediment suspension in detail. Strong coherent vortical structures are well developed along the upslope of the ripple surface during the accelerating flow phase, which effectively drag the particles to the ripple crest. At the maximum flow rate and at the decelerating flow phase, a cloud of vortical structures is developed vertically in the lee area of the ripple. Sediment particles render strong dispersion in the vertical direction when they are captured by these turbulent vortices, causing convective sediment flux that cannot be explained by the mean flows. The convective sediment suspension is strongest at the time of flow deceleration, while over a flat bed at the time of flow reversal. This observation suggests that bed form effect should be considered in modeling convective sediment flux.     

Thermohydrodynamics of the ocean decay of nonlinear oscillations in a bounded basin of variable depth

We study the time decay of surges of a liquid in a round shallow-water basin of variable depth. The dependence of the logarithmic decrement of oscillations on the bottom topography and wind velocity is analyzed. The role of convective acceleration and bottom friction in the formation of both the level of vertical displacement of the surface of the basin and the velocity field of horizontal wave currents is estimated. __________ Translated from Morskoi Gidrofizicheskii Zhurnal, No. 2, pp. 3–11, March–April, 2006.     

Least Square Finite Element Method for Viscous Splitting of Unsteady Incompressible Navier–Stokes Equations

In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton's method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C~0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.     

Effect of the tangential component of a force field on convection in a rotating plane layer

Thermal convection in an inclined plane layer having boundaries with different temperatures and rotating around an axis perpendicular to its plane is studied experimentally. It is shown that the convection and heat transfer are determined by two different convective mechanisms—gravitational and thermovibrational ones; they manifest themselves in the threshold excitation of cellular convective structures with different sizes. The thermovibrational convection is caused by the period-averaged mass force arising as a result of tidal oscillations of a nonisothermal fluid with respect to the cavity under the action of the external force-field component tangential to the layer. The map of convection regimes on the plane of governing parameters—the gravitational Rayleigh number and its vibrational analog—has been constructed. It has been found that the thermovibrational convection can develop even in a layer heated from above. The role of the dimensionless velocity of rotation, which has a stabilizing effect on the excitation of both gravitational and thermovibrational convection, has been studied.     

THE DIAGNOSTIC CALCULATION OF EQUATORIAL FLOWS IN THE EASTERN PART OF THE WESTERN PACIFIC OCEAN

On the basis of observational data of the eastern part of the West Pacific Ocean, a diagnostic calculation of equatorial flow for this region is performed by using the authors' model equations and computing scheme and methods. For the first cruise (January 3-March 4, 1979), the results show: (i) The primary driving force of the equatorial surface flows comes from the prevailing northeasterly wind field, with an average uniform wind speed Vw = 6.3 m/s. The steady westward wind produces divergent westward flows in the surface layers, causing an upwelling near the equator. The importance of the steady wind stress in the upper layer (120 m) decreases with depth and becomes insignificant at the level of z = 75 m or z = 100 m, (ii) The equatorial undercurrent is a strong eastward and equator-ward flow, with its eastward component of undercurrent larger than its meridional component. The core of the undercurrent is at the thermocline, and its maximum velocity is 88-90 cm/s at the level of z=200 m. The deeper f     

华北冷涡背景下青岛三次混合型对流天气环境场条件对比分析

利用探空观测资料,对比分析华北冷涡背景下青岛三次混合型对流天气过程环境场条件,揭示出现短时强降水、雷暴大风和冰雹天气时的水汽、稳定度和垂直风切变差异特征.分析结果表明:500 hPa上冷涡中心位于42°N的华北冷涡、850 hPa低涡系统和偏南风急流以及地面气旋是这三次混合型对流天气的影响系统;在这三次混合型对流过程中...     

Models of near-bed dynamics and sediment movement at the Iberian margin

A three dimensional time-dependent baroclinic hydrodynamic model, including sediment transport and incorporating a turbulence energy sub-model, is used in cross sectional form to examine sediment movement at the shelf edge off North West Iberia at 42°40.5’N where measurements were made as part of the OMEX-II programme. These calculations are complemented by a simpler, in essence time-independent model, which is used to examine the sensitivity of the sediment distribution over the slope (from a shelf-break source) to changes in the specified values of horizontal and vertical diffusion coefficients. The philosophy of the paper is to use idealized tidal, wind and wind wave forcing to examine changes in sediment distribution resulting from these processes. Calculations with the time-dependent and steady state models give insight into both the role of events and long-term effects. The steady state model focuses on the off-shelf region, whilst the time-dependent model considers on-shelf events.Tidal calculations showed that for the stratification used here the internal tide in the OMEX region was primarily confined to the shelf edge and ocean. A mean on-shelf sediment transport in the surface layer and off-shelf transport at the bed was found. Across-shelf circulations produced by up-welling/down-welling favourable winds gave rise to on-shelf/off-shelf currents in the bottom boundary layer with an opposite flow in the surface layer. In the case of an up-welling favourable wind, sediment suspension was at a maximum in the near coastal region, with sediment being advected off shore in the surface layer. With a down-welling favourable wind, surface sediment was advected towards the shore, but there was offshore transport at the bed. Near the shelf edge any upwelling flow had the tendency to return this sediment to the surface layer from whence it was transported on-shore. So in essence the sediment was trapped within an on-shelf circulation cell. Wind waves effects increased the total bed stress and hence the sediment concentration and its transport, although its pattern was determined by tidal and wind forcing.The time independent model with increased/decreased lateral diffusivity gave an enhanced/reduced horizontal sediment distribution for a given settling velocity. As the settling velocity increases, the down-slope movement of sediment is increased, with a reduction in the thickness of the near-bed sediment layer, but with little change in its horizontal extent.     

Modeling a neutrally stratified turbulent air flow over a horizontal rough surface

The neutrally stratified boundary layer over a smooth rough surface is consider. The turbulent flow is simulated using a finite-difference eddy-resolving model of the atmospheric boundary layer (ABL). The model includes different turbulence closure schemes and numerical approximations for advection components of the momentum balance equation. We investigate the quality of reproduction of spectral characteristics of the turbulent flow and the model’s capabilities to reproduce the observed profile of mean wind velocity near the rough surface. It is shown that the best result is obtained by coupling a numerical scheme of higher order of accuracy with a mixed closure scheme based on an adaptive estimation of the mixing length for subgrid-scale fluctuations. Here, we are able to reproduce the asymptotics of the fluctuation spectrum of the longitudinal component of wind velocity near the surface and within the boundary layer as well as the logarithmic profile of mean velocity near the surface.     

Thermochaline Convection in the Edge-Ice Zone in the Barents Sea to the East of Spitsbergen

The analysis of the thermochaline (double-diffusive) convection in the edge-ice zone of the Barents Sea to the east of Spitsbergen was performed by using the data of two hydrological surveys of the Norwegian Polar Institute in 1999 and 2000. The Turner's model, where the critical Rayleigh number is determined by using the empirical relationship between the Nusselt number and the Rayleigh number normalized to its critical value, was used for analysis. The results of calculation of thickness of the convective layer were compared with the data of supervision. The estimates of the vertical convective velocity were obtained.     

Parameterization of the depth of the upper quasi-homogeneous layer based on measurements in the North Atlantic (Section along 53° N)

The dependence of the variation in the depth of the upper mixed layer (MLD) on the governing parameters (the momentum flux, the buoyancy fluxes at the ocean surface, and the density gradient in the pycnocline) is considered. It is shown that, in the spring storm season, wind mixing dominates over convective mixing. In this case, the MLD is linearly correlated with the Ekman scale calculated from the friction velocity observed approximately 12 h before the measurement of the MLD.     

Internal flow structure of short wind waves

The minimum value of wind stress under which the flow velocity in short wind waves exceeds the phase speed is estimated by calculating the laminar boundary layer flow induced by the surface tangential stress with a dominant peak at the wave crest as observed in previous experiments. The minimum value of the wind stress is found to depend strongly on, the ratio of the flow velocity just below the boundary layer and the phase speed, but weakly onL, the wavelength. For wind waves previously studied (=0.5,L=10 cm), the excess flow appears when the air friction velocityu _* is larger than about 30 cm sec^–1. The present results confirm that the excess flow found in my previous experiments is associated with the local growth of a laminar boundary layer flow near the wave crest.     

OBS海洋环境信号分析与应用

为研究海洋环境信号在OBS(Ocean Bottom Seismograph)原始数据中的规律及应用,根据OBS原始数据的波形及频谱特征,将研究区划分为5个时间段,依次为旧涌浪阶段、风浪渐强阶段、风浪全盛阶段、风浪消退阶段和新涌浪阶段。结合海洋天气预报,认为上述现象是由偏南风风浪对海流的影响造成的。参考野外地震数据采集记录班报,得到各阶段的时长和距离,计算风浪渐强、全盛和消退阶段OBS附近海流的平均速度。结果表明:OBS原始资料中浅海海洋环境噪音增强的主要因素是风浪,且风浪引起的噪音信号的波形变化特征是渐进式的;OBS可用于接收某种特殊阶段(如台风、海啸等)的噪音信号,并根据噪音信号的波形特征、频谱变化规律和持续时间估算该阶段的海流速度变化。     

Spatial spectra and characteristic horizontal scales of temperature and velocity fluctuations in the convective boundary layer of the atmosphere

Using a high-resolution LES numerical model, we calculated the turbulent thermal convection for high ratios of horizontal and vertical sizes of the computational domain (26: 26: 1). The natural analog of the simulated process is a planetary boundary layer (PBL) of the atmosphere growing with height in the background of stably stratified overlying air layers over a horizontally homogeneous heated surface under a weak average wind. We obtained the spectral distributions of variances of fluctuations in potential temperature and velocity components in ranges corresponding to scales from a few tens of meters to a few tens of kilometers. We found energetically significant segments of the spectrum of large-scale fluctuations in the potential temperature for which the power dependences S ～ k ^?1/3 and S ～ k ^?4/3 are satisfied with good accuracy. We calculated the characteristic spatial scales of horizontal fluctuations in velocity and temperature. We obtained a dependence of these scales on the height of the growing convective PBL. We discuss the characteristic features of large-scale distributions in terms of the self-similarity of the growing boundary layer behavior.