长江口湍流剖面的观测与分析

小尺度湍流过程对河口物质输运与能量交换至关重要。受传统观测方法的限制, 河口浅水区域的剖面观测资料至今较为匮乏, 进而限制了湍流过程的研究。为此, 采用新型5波束声学多普勒流速剖面仪(Nortek Signature 1 000 kHz AD2CP)在长江口开展湍流剖面观测, 获取高频、低噪的高质量湍流剖面数据, 并与声学多普勒点式流速仪(acoustic doppler velocimeters, ADV)同步观测的数据进行对比。结果表明, 通过AD2CP与ADV获得的近底部边界层摩阻流速u_*、拖曳系数C_d、雷诺应力S_R等特征参数基本一致, 底摩擦与波浪能量为河口区域湍动能的主要输入源。湍流垂向结构存在显著的非局地平衡, 即温盐等斜压作用引起的浮力通量、对流项以及强波浪作用影响的脉动压力做功、黏性输运等因素可能为长江口湍流非局地平衡的主因。     

海洋白浪寿命的定义及测量结果

通过对国内外白浪研究和应用的分析,首次提出了有效白浪寿命的定义,给出了计算白浪寿命的公式及测量方法和结果,并报告了以此方法在渤海实测的结果,得到了白浪寿命ＴＬ与海面风速Ｕ１０的关系为ＴＬ＝０．２６Ｕ１０以及白浪寿命概率分布近于瑞利分布等。     

A note on a critical wind speed for air–sea boundary processes

Wind and wind-generated waves were measured in a wind-wave tank. A clear transition was found in the relation between the wind speed U ₁₀ and the wind friction velocity u _* near u _* = 0.2 m/s, where U ₁₀ is the wind speed at 10 m height extrapolated from the measured wind profile in a logarithmic layer, and u _* = 0.2 m/s corresponds roughly to U ₁₀ = 8 m/s in the present measurement. Quite a similar transition was found in the relation between the spectral density of high frequency wind waves and u _*. These results suggest the existence of the critical wind speed for air–sea boundary processes, which was proposed by Munk (J Marine Res 6:203–218, 1947) more than half a century ago. His original idea of the critical wind speed was based on the discontinuities in such phenomena as white caps, wind stress, and evaporation, which commonly appear at a wind speed near 7 m/s. On the basis of the results of our present study and those of earlier studies, we discuss the phenomena which are relevant to the critical wind speed for the air–sea boundary processes. The conclusion is that the critical wind speed exists and it is attributed to the start of wave breaking rather than the Kelvin–Helmholtz instability, but the air–sea boundary processes are not discontinuous at a particular wind speed; because of the stochastic nature of breaking waves, the changes occur over a range of wind speeds. Detailed discussions are presented on the dynamical processes associated with the critical wind speed such as wind-induced change of sea surface roughness and high frequency wave spectrum. Future studies are required, however, to clarify the dynamical processes quantitatively. In particular, there is a need to further examine the gradual change of breaking patterns of wind waves with the increase of wind speed, and the associated change of the structure of the wind over wind waves, such as separation of the airflow at the crest of wind waves, the turbulent stress, and wave-induced stress. Studies on the dynamical structure of the high frequency wave spectrum are also needed.     

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.     

Comparison of Ocean Surface Wind Stress Computed with Different Parameterization Functions of the Drag Coefficient

Depending on the choice of reference wind speed, the quantitative and qualitative properties of the drag coefficient may vary. On the ocean surface, surface waves are the physical roughness at the air-sea interface, and they play an important role in controlling the air-sea exchange processes. The degree of dynamic influence of surface waves scales with wavelength. Drag coefficient computed with the reference wind speed at an elevation proportional to the wavelength (for example, U _λ/2) is fundamentally different from the drag coefficient computed with the wind speed at fixed 10 m elevation (U ₁₀). A comparison has been carried out to quantify the difference in wind stress computation using several different parameterization functions of the drag coefficient. The result indicates that the wind stress computed from U ₁₀ input using a drag coefficient referenced to U _λ/2 is more accurate than that computed with drag coefficient functions referenced to U ₁₀.     

地形和风速影响下的海气相互作用大涡模拟研究

当水流经过海洋地形时,水流的不稳定性会引起垂向混合并伴随大量湍流过程。针对传统海气耦合模式缺少在湍流尺度上讨论海洋地形与风速对海气相互作用影响的问题,使用并行大涡模拟海气耦合模式(the parallelized large eddy simulation model, PALM)在5 m/s的背景风场下,引入理想立方体地形,对比有无地形的影响;设置地形边长为L,高为3L (其中大气部分高L), L与水深H之比为L/H=1/2;然后保持地形条件不变。设置5、10和15 m/s三种风速,讨论风速对小尺度海气相互作用的影响。研究表明:地形在大气部分减弱顺风向速度,增强侧风向速度,影响0～5L的高度区域,而对垂向作用较小;无地形条件下湍流垂向涡黏系数K_m在-0.3L时,水深达到最大值0.024 m²/s,有地形条件下K_m在-0.8L时,达到最大值为0.16 m²/s,地形的存在使得上层海洋混合加强, K_m最大值增加1个数量级。随风速增大海洋和大气中的净热通量、淡水通量和浮力通量都相应...     

Sediment distribution and transport across the continental shelf and slope under idealized wind forcing

Resuspension, transport, and deposition of sediments over the continental shelf and slope are complex processes and there is still a need to understand the underlying spatial and temporal dynamical scales. As a step towards this goal, a two-dimensional slice model (zero gradients in the alongshore direction) based on the primitive flow equations and a range of sediment classes has been developed. The circulation is forced from rest by upwelling or downwelling winds, which are spatially uniform. Results are presented for a range of wind speeds and sediment settling speeds. Upwelling flows carry fine sediments (low settling speeds) far offshore within the surface Ekman layer, and significant deposition eventually occurs beyond the shelf break. However, coarser sediments quickly settle out of the deeper onshore component of the circulation, which can lead to accumulation of bottom sediments within the coastal zone. Downwelling flows are more effective at transporting coarse sediments off the shelf. However, strong vertical mixing at the shelf break ensures that some material is also carried into the surface Ekman layer and returned onshore. The concentrations and settling fluxes of coarse sediments decrease offshore and increase with depth under both upwelling and downwelling conditions, consistent with trends observed in sediment trap data. However, finer sediments decrease with depth (upwelling) or reach a maximum around the depth of the shelf break (downwelling). It is shown that under uniform wind conditions, suspended sediment concentrations and settling fluxes decay offshore over a length scale of order τ_s/ρf|w_s|, where τ_s is the wind stress, ρ the water density, f the Coriolis parameter, and w_s is the sediment settling velocity. This scaling applies to both upwelling and downwelling conditions, provided offshore transport is dominated by wind-driven advection, rather than horizontal diffusion.     

Estimation of the leeway drift of small craft

Small craft (<6·4 m) leeway is determined as a function of the wind speed in the range of 5–20 knots (3·6–10·3 m/sec). Leeway is calculated relative to the surface current by measurement of the separation distance of the small craft from a dyed patch of surface water at sea, using time-sequenced aerial photography. Leeway increases linearly with wind speed for small craft equipped with or without a sea anchor in the wind range studied. Leeway for small craft without sea anchor can be calculated from the equation U_L = 0.07 U_W + 0.04 where U_W is the wind speed at 2 m elevation. Leeway for small craft drifted off the be calculated from the equation U_LD = 0·05 U_W − 0·12. The small craft drifted off the downwind direction in about 80% of the experiments. The drift angle is variable and difficult to predict.     

A global wave parameter database for geophysical applications. Part 1: Wave-current–turbulence interaction parameters for the open ocean based on traditional parameterizations

Ocean surface mixing and drift are influenced by the mixed layer depth, buoyancy fluxes and currents below the mixed layer. Drift and mixing are also functions of the surface Stokes drift U_ss, volume Stokes transport T_S, a wave breaking height scale H_swg, and the flux of energy from waves to ocean turbulence Φ_oc. Here we describe a global database of these parameters, estimated from a well-validated numerical wave model, that uses traditional forms of the wave generation and dissipation parameterizations, and covers the years 2003–2007. Compared to previous studies, the present work has the advantage of being consistent with the known physical processes that regulate the wave field and the air–sea fluxes, and also consistent with a very large number of in situ and satellite observations of wave parameters. Consequently, some of our estimates differ significantly from previous estimates. In particular, we find that the mean global integral of Φ_oc is 68 TW, and the yearly mean value of T_S is typically 10–30% of the Ekman transport, except in well-defined regions where it can reach 60%. We also have refined our previous estimates of U_ss by using a better treatment of the high frequency part of the wave spectrum. In the open ocean, U_ss  0.013U₁₀, where U₁₀ is the wind speed at 10 m height.     

基于大涡模拟耦合模式的小尺度海气相互作用研究

大气和海洋是影响地球气候系统的两个重要因素，它们之间的相互作用是海洋和大气研究的重要课题，海气耦合模式则是研究海气相互作用的重要工具，而海气耦合模式重点考虑的参数是海气通量。针对传统的大尺度海气耦合模式缺少湍流尺度分析的问题，本文使用并行大涡模拟海气耦合模式（The Parallelized Large-Eddy Simulation Model，PALM），在小尺度上探究风速对海气通量及湍流动能收支（Turbulence Kinetic Energy Budget，TKE Budget）的影响，设置了5、10和15m/s三种地转风速度对大气边界层（Atmospheric Boundary Layer，ABL）和海洋混合层（Oceanic Mixed Layer，OML）进行海气耦合模拟。研究表明：海气通量的分布与风速大小密切相关，风速越大，净热通量和浮力通量相对越大，由于温度上升导致海水蒸发加剧，使得大气的淡水通量增大；海洋湍流动能收支各项在近海面处受风速影响较大，且随着深度加深而逐渐减弱。本研究初步展示了小尺度海气耦合模式在海气通量研究中的应用，对进行小尺度海气相互作用研究具有一定的意义。     

Modeling the eddy transport of momentum and heat: Comparison with direct measurements in free atmosphere

Direct measurements of eddy diffusivities for momentum K _m and heat K _h by Doppler radar and by a radio acoustic sounding system in the upper troposphere and lower stratosphere were used to examine the applicability of three Reynolds-averaged Navier-Stokes (RANS) schemes of stratified turbulence in the environment: the E — ? turbulence scheme modified for stratified flows, the algebraic two-parameter E — ? Reynolds-stress scheme, and the three-parameter \(E - \varepsilon - \overline {\theta ^2 } \) turbulence scheme. All turbulence parameters-the turbulent kinetic energy (E), the dissipation rate (?), and vertical profiles of potential temperature (atmospheric stability) and mean wind velocity-were derived from direct measurements for all three turbulence schemes. It is shown that the profile of the vertical diffusivity of momentum (K _m ) obtained from the three-parameter RANS turbulence scheme agrees well with its directly measured analog. The profile of K _m calculated by the two-parameter turbulence schemes fits measurements rather qualitatively.     

A numerical study of three-dimensional wave interaction with a square cylinder

In this study, a three-dimensional numerical model is used to study the wave interaction with a vertical rectangular pile. The model employs the large eddy simulation (LES) method to model the effect of small-scale turbulence. The velocity and vorticity fields around the pile are presented and discussed. The drag and inertial coefficients are calculated based on the numerical computation. The calculated coefficients are found to be in a reasonable range compared with the experimental data. Additional analyses are performed to assess the relative importance of drag and initial effects, which could be quantified by the force-related Keulegan and Carpenter (KC) number: KC_f=UT/(4πL). Here U is the maximum fluid particle velocity, T the wave period and L the length of structure aligned with the wave propagation direction. For small KC_f, the effective drag coefficient is proportional to 1/KC_f, provided the wavelength is much longer than the structural length. When wavelength is comparable to the structure dimension, the effective drag coefficient would be reduced significantly due the cancellation of forces, which has been demonstrated by numerical results.     

Adjustment of Wind Waves to Sudden Changes of Wind Speed

An experiment was conducted in a small wind-wave facility at the Ocean Engineering Laboratory, California, to address the following question: when the wind speed changes rapidly, how quickly and in what manner do the short wind waves respond? To answer this question we have produced a very rapid change in wind speed between U _low (4.6 m s^?1) and U _high (7.1 m s^?1). Water surface elevation and air turbulence were monitored up to a fetch of 5.5 m. The cycle of increasing and decreasing wind speed was repeated 20 times to assure statistical accuracy in the measurement by taking an ensemble mean. In this way, we were able to study in detail the processes by which the young laboratory wind waves adjust to wind speed perturbations. We found that the wind-wave response occurs over two time scales determined by local equilibrium adjustment and fetch adjustment, Δt ₁/T = O(10) and Δt ₂/T = O(100), respectively, in the current tank. The steady state is characterized by a constant non-dimensional wave height (H/gT ² or equivalently, the wave steepness for linear gravity waves) depending on wind speed. This equilibrium state was found in our non-steady experiments to apply at all fetches, even during the long transition to steady state, but only after a short initial relaxation Δt ₁/T of O(10) following a sudden change in wind speed. The complete transition to the new steady state takes much longer, Δt ₂/T of O(100) at the largest fetch, during which time energy propagates over the entire fetch along the rays (dx/dt = c _g) and grows under the influence of wind pumping. At the same time, frequency downshifts. Although the current study is limited in scale variations, we believe that the suggestion that the two adjustment time scales are related to local equilibrium adjustment and fetch adjustment is also applicable to the ocean.     

沿岸海域三维斜压场的数值模拟 Ⅱ、渤海湍流能量数值计算

研究建立渤海海域三维斜压场的湍流能量模型。对计算网格无法分辩的湍流运动，引入湍流能量求得垂向动量、质量及能量交换系统，从而形成封闭的湍流运动方程组。定量地计算了湍流能量、湍流速度及有风、无风状况下湍流能量的垂直分布等，探讨了湍流运动对潮波运动的影响。     

粗糙海面L 和C 双波段的代价函数多参量遥感反演分析

利用代价函数(cost function)方法,通过分析粗糙海面L和C双波段多极化遥感亮温对海表盐度、温度、风速和有效波高等参数的敏感性以及L和C双波段多极化的代价函数收敛特性,建立了反演海表盐度、温度、风速和有效波高等多参数的L和C双波段多极化代价函数模式。双波段遥感模式分析结果表明:(1)对于双参数的联合反演,L和C双波段垂直极化代价函数联合反演海表盐度和温度可以获得较好的反演结果。(2)L波段垂直极化和C波段水平极化代价函数联合反演海表盐度和风速较好。(3)对于三参数联合反演,L波段垂直极化和C波段的双极化联合反演盐度、温度和风速的精度较高。(4)L波段亮温对有效波高的敏感性较低(C波段经验模式不含有效波高),使得有效波高反演误差较大,L和C波段经验模式不适合反演有效波高参数。另外,为了定量分析L和C双波段代价函数的多参量遥感反演结果,采用加性噪音模拟亮温方法,对上述L和C双波段多极化模式的盐度、温度和风速等多参数联合反演误差进行了分析,均得出较好的结果。结论表明L和C双波段代价函数联合反演多参量可以明显提高参量反演精度,为粗糙海表面多参量的反演提供了新的方法和途径。     

A pair of Langmuir cells in two laboratory tanks (II) on generation mechanism

Vertical and cross-wind profiles of mean currents were measured systematically in vertical cross-sections of two wind-wave tanks with aspect ratios of order one to study the secondary flow in the tanks. A pair of Langmuir cells turned out to be driven by a close combination of the pressure gradient along the tank and the side-wall effects. That is, part of the adverse pressure gradient produced a parabolic cross-wind profile with the smallest downwind current at the centerline and the largest current along the two sidewalls. As a result, upwelling occurred in the center zone where the return flow was strongest, probably because of the entrainment action of the wind-driven current. In order to compensate for this upwelling, downwelling occurred along the two side-walls from the flow continuity. The resulting vertical circulation formed a pair of Langmuir circulations across the span and served to maintain the parabolic profile formed by the pressure gradient. A positive feedback mechanism is thus found between the primary and secondary circulations through upwelling of the return flow in the center zone. Vertical shears of the span-averaged downwind current measured in two tanks were found to be systematically different from each other. This difference seems to depend on the magnitude of the advective Reynolds stresses in the two tanks.     

1013号“鲇鱼”台风大气边界层低层动量通量逆梯度输送的观测特征

为了进一步揭示大气边界层低层动量通量逆梯度输运的观测事实,本文利用福建南部赤湖镇海边100m铁塔上观测得到的超声风温仪资料,采用诊断方法分析了2010年13号"鲇鱼"台风影响过程中近海面边界层大气中动量垂直通量逆梯度输送特征。结果表明:在逐小时时间尺度上台风近海面边界层大气中湍流动量垂直通量输送以沿梯度方向为主,但也存在一定比例的逆梯度方向输送现象,且越往高层该现象越明显,本文中其百分比最高不超过16%。并且这一现象出现的频率在台风环流的不同区域存在差别:在台风环流的内核区域最多,其次是登陆以后残余环流区,而在早期的台风外围环流中较少。湍流低频扰动的相干结构也是影响因素,当水平扰动与垂直扰动出现同位相变化时容易出现动量垂直通量逆梯度输送。大气动量通量逆梯度的水平和垂直平均空间尺度分别为258m和35m,平均时间尺度分别为123s和13s,均小于通常的顺梯度湍流低频扰动的时空尺度的1—3倍。     

沿海市区近地面大气层风、温、湿湍流特征的观测和分析

用三向风速仪和高灵敏温湿度脉动观测装置,测量沿海城市上空大气表面层风速、温度和湿度脉动。并用涡动相关法计算了大气表面层动量、热量和水汽湍流通量,给出了复杂地形上各物理量的谱特征,讨论了各湍流通量的特性。结果指出,风速、温度和湿度谱密度分布基本上都服从—5/3律。由于建筑物的影响,在f≥0.2Hz频段内拢动能量明显偏大。动量和热量输送主要由f<0.2Hz的扰动完成的,高频扰动的贡献可以忽略不计。     

Polar lows and tropical hurricanes: Their energy and sizes and a quantitative criterion for their generation

Similarity and dimension considerations applied to convection in a rotating fluid allows one to estimate the sizes and horizontal velocities of generated vortices. To do this, it is necessary to know the buoyancy flux in the fluid and the angular velocity of fluid rotation [1, 2]. The author’s preliminary efforts [3] have shown that the sizes, wind speeds, and total kinetic energy can thus be estimated correctly for tropical cyclones (TCs), as well as for polar lows (PLs) (which are often called explosive mesocyclones because they take just a few hours to develop). In this study, the sensible and latent heat fluxes for U = 33 m/s and the related buoyancy fluxes are estimated on the basis of climatology, bulk formulas, and the velocity scale of convection in a rotating fluid. In the tropics, at hurricane wind speeds U ≥ 33 m/s and climatological air humidity r = 80%, the total heat flux at the water surface temperature T _s ≥ 26°C becomes equal to or greater than 700 W/m². Due to the Clausius-Clapeyron equation, the latent heat flux to the atmosphere (the main part of the flux in the tropics) decreases substantially at lower values of T _s. Thus, an energy flux from the ocean to the atmosphere of 700 W/m² or greater should be regarded as the first necessary condition for TC genesis instead of the temperature T _s. Low static stability, which must be at least half its climatological value as estimated here, is another necessary condition [4]. In polar regions, total fluxes roughly twice those in the tropics are needed for the formation of explosive mesocyclones, PLs, which is explained by the much smaller role of latent heat, greater geostrophicity, and stronger static stability of the atmosphere there. Enthalpy fluxes and wind speeds are interrelated: the larger the flux is, the stronger the convection, the higher the concentration of angular momentum in an ascending convective air column, and the greater the azimuthal velocity in the vortex are, which in turn enhances the transfer of energy from the ocean. Considering the problem with the use of simple analytic relations makes it possible, for the first time, to find a numerical criterion for their generation. It is hoped that this material may be useful for educational purposes as well.