夜间稳定边界层中小尺度地形激发的形式阻力和波动阻力

通过求解含有摩擦耗散的线性化大气动力学方程组,得到了在夜间稳定大气边界层中小尺度地形产生的波动阻力和形式阻力的解析解.结果表明边界层中的稳定度、风速和湍流状态、边界层厚度、上部残余层中的稳定度和风速以及地形高度和坡度,都会影响波动阻力和形式阻力的大小,应在数值模式的参数化方案中给予考虑.分析还表明,当地形坡度减到一定程度时,形式阻力可以忽略不计.     

Some effects of the air-water interface on gravity waves

Abstract

New rates of decay are presented for temporally-attenuated gravity waves in deep water, allowance being made for the energy dissipated in the Stokes interfacial boundary layer in the air. This decay-rate, involving air drag, may then be used to deduce a new “free-surface” boundary condition for the problem of the mass transport velocity due to progressive waves; for shallow-water waves, two specific velocity profiles are calculated, and indicate large differences in comparison with the corresponding profiles of Longuet-Higgins (1953) for a vacuum-water interface.     

Prediction of the drag law for air–sea momentum exchange

In this paper, we use the inertial coupling relation as a similarity model for the air–sea boundary layer, to predict the 10-m drag coefficient. Excellent agreement with the commonly used statistical relationship of Garratt (1992) is found for a fully developed growing wind wave sea with a constant inertial drag coefficient, K_I = 1.5 × 10^–3. This suggests that the inertial coupling model can be used to realistically predict the 10 m drag coefficient under more general wind wave conditions.Acknowledgements The paper was completed while JATB was a Fellow at the Hanse-Wissenschaftskolleg in Delmenhorst, Germany, in July and August 2004. The comments of two anonymous reviewers are gratefully acknowledged.     

Charnock dynamics: a model for the velocity structure in the wave boundary layer of the air–sea interface

We present a unified model of the air–sea boundary layer, which takes account of the air–sea momentum exchange across the sea surface. The recognition of the importance of the velocity shears in the water (which comprise a frictional shear and the Stokes shear due to the wave motion) in determining the sea surface roughness is a distinctive feature of the analysis, which leads to a prediction of the Charnock constant (α) in terms of two independent parameters, namely the wave age and the ratio of the Stokes shear to the Eulerian shear in the water. This expression is used to interpret the large observational variability of the Charnock constant. The 10-m drag coefficient can also be expressed using similar reasoning, and the introduction of a relation in which the ratio of the frictional shear in the water to the frictional shear in the air decreases with the friction velocity yields predictive relations for the variation of the 10-m drag coefficient at very high wind speeds both in the open ocean and in wind–wave tanks. The physical interpretation of this relation is that the production of spray essentially returns momentum from the ocean to the atmosphere, and this process becomes progressively more important as the wind speed increases.     

On the variability of the Charnock constant and the functional dependence of the drag coefficient on wind speed

A model for the air–sea interface, based on the coupled pair of similarity relations for “aerodynamically” rough flow in both fluids, is presented, which is applied to fetch-limited and high wind speed conditions which occur, for example, in hurricanes. It is shown that the specification of the maximum 10-m drag coefficient and the 10-m wind speed and the peak wave speed at which it occurs are sufficient to uniquely determine the drag law, which asymptotes at low wind speeds to a Charnock constant similar to that for the fully developed wind wave sea and is almost independent of the peak wave speed at the maximum in drag coefficient. A feature of the drag law is that it is of Charnock form, almost independent of the wave age, consistent with the transfer of momentum to the wave spectrum being due to the smaller rather than the dominant wavelengths. The analysis is also applied to a variable sea state in which either the surface wind or the surface Stokes drift vary, but the peak wave speed is kept constant. The corresponding variability in the Charnock constant is in general accord with observations.     

On the variability of the Charnock constant and the functional dependence of the drag coefficient on wind speed: Part II-Observations

An analytical expression for the 10 m drag law in terms of the 10 m wind speed at the maximum in the 10 m drag coefficient, and the Charnock constant is presented, which is based on the results obtained from a model of the air-sea interface derived in Bye et al. (2010). This drag law is almost independent of wave age and over the mid-range of wind speeds (5?17 ms^?1) is very similar to the drag law based on observed data presented in Foreman and Emeis (2010). The linear fit of the observed data which incorporates a constant into the traditional definition of the drag coefficient is shown to arise to first-order as a consequence of the momentum exchange across the air-sea boundary layer brought about by wave generation and spray production which are explicitly represented in the theoretical model.     

利用瑞利激光雷达观测北京地区上平流层地形重力波活动

本文利用中国科学院空间科学与应用研究中心的瑞利激光雷达首次观测到了平流层地形重力波活动的现象,并结合美国国家环境预报中心(NCEP)的全球预报系统(GFS)的风场数据分析了该地形重力波的基本参数.与惯性重力波相比较,地形重力波的密度扰动没有下传的相位,在同一高度上,其扰动相位保持不变.北京空间科学与应用研究中心瑞利激光雷达自2012年开始观测实验以来,已经观测到多起地形重力波活动事件.本文以2013年11月11日的观测数据为例,研究北京上空的地形重力波活动,并结合GFS风场数据分析了北京上平流层地形重力波的波长、传播方向、传播速度等参量.通过分析得到在2013年11月11日北京上空存在一列传播方向为北偏西52.4°,水平波长为5.5km,平均垂直波长约为6.0km的地形重力波.     

Airside velocity measurements over the wind-sheared water surface using particle image velocimetry

An experimental investigation of the airflow structure in the near surface region over the wind-sheared air–water interface is reported. The two-dimensional velocity fields in a plane perpendicular to the water surface were measured using particle image velocimetry (PIV) technique over a wind speed range from 1.5 to 4.4 m s⁻¹. The results show a reduction in the mean velocity magnitudes and the tangential stresses when gravity waves appear on the surface. An enhanced vorticity layer was observed immediately above the water surface that extended to a height of approximately 2 cm. The vorticity was enhanced by an order of magnitude, and the energy dissipation rate was enhanced by a factor of 7 in this layer at all wind speeds. The vertical profiles of Reynolds stress, energy production, and dissipation indicate the contribution of surface waves in the enhanced transfer of momentum and energy between the two fluids. The results in this study show that the flow dynamics in a layer immediately adjacent to the water surface whose thickness is of the order of the significant wave height is significantly different from that at greater heights. Thus, it is concluded that the quantitative investigation of the flow in the immediate vicinity of the interface is vital for an improved understanding of the heat, mass, and momentum exchange between air and water. The present study demonstrates that PIV is an effective technique to accurately measure the velocity fields in this region.     

Vertical structure of the stable boundary layer detected by RASS-SODAR and in-situ measurements in SABLES 2006 field campaign

Data from the SABLES 2006 field campaign are used in order to analyse some of the main processes present along the nocturnal periods: surface-based inversions, low level jets, katabatic winds, wave-like motions, pressure perturbations, etc. These processes have an important influence on the vertical structure (both thermal and dynamical) of the atmospheric boundary layer, and can be better described with the synergetic combination of RASS-SODAR data and in-situ measurements (such as sonic anemometer data and high-resolution pressure series from microbarometers). It is shown how the different air masses and their evolution are easily identified when pressure and RASS-SODAR wind and temperature data are presented together. Likewise, periodic pressure fluctuations observed in the surface array of microbarometers reveal the existence of gravity wave motions whose propagation is better understood after locating the wave ducting layers with the help of RASS-SODAR average wind ant temperature profiles.     

Observations and modeling of wave-supported sediment gravity flows on the Po prodelta and comparison to prior observations from the Eel shelf

A mooring and tripod array was deployed from the fall of 2002 through the spring of 2003 on the Po prodelta to measure sediment transport processes associated with sediment delivered from the Po River. Observations on the prodelta revealed wave-supported gravity flows of high concentration mud suspensions that are dynamically and kinematically similar to those observed on the Eel shelf [Traykovski, P., Geyer, W.R., Irish, J.D., Lynch, J.F., 2000. The role of wave-induced density-driven fluid mud flows for cross-shelf transport on the Eel River continental shelf. Continental Shelf Research 20, 2113–2140]. Due to the dynamic similarity between the two sites, a simple one-dimensional (1D) across-shelf model with the appropriate bottom boundary condition was used to examine fluxes associated with this transport mechanism at both locations. To calculate the sediment concentrations associated with the wave-dominated and wave-current resuspension, a bottom boundary condition using a reference concentration was combined with an “active layer” formulation to limit the amount of sediment in suspension. Whereas the wave-supported gravity flow mechanism dominated the transport on the Eel shelf, on the Po prodelta flux due to this mechanism is equal in magnitude to transport due to wave resuspension and wind-forced mean currents in the cross-shore direction. Southward transport due to wave resuspension and wind forced mean currents move an order of magnitude more sediment along-shore than the down-slope flux associated wave-supported gravity flows.     

Generalized polarization relationships for acoustic gravity waves in the nonisothermal atmosphere with the wind

The nonlinear equation of the first order of the Riccati type has been obtained for the wave impedance of acoustic gravity waves in the nonisothermal atmosphere. The vertically nonuniform horizontal wind and the effect of viscosity on the horizontal components of the velocity field have been taken into account in the calculations. The boundary-value problem for the Riccati equation is defined by the boundary emission condition at high altitudes. Upon finding the wave impedance along with the generalized polarization relationship, all remaining disturbances of the atmospheric parameters related to acoustic gravity waves are found with the help of a simple integration. The results of using a developed formalism are illustrated by the numerical computation of acoustic gravity wave fields in the atmosphere with real vertical profiles of temperature and horizontal field velocity.     

Effects of the bottom boundary condition in numerical investigations of dense water cascading on a slope

The flow of dense water along continental slopes is considered. There is a large literature on the topic based on observations and laboratory experiments. In addition, there are many analytical and numerical studies of dense water flows. In particular, there is a sequence of numerical investigations using the dynamics of overflow mixing and entrainment (DOME) setup. In these papers, the sensitivity of the solutions to numerical parameters such as grid size and numerical viscosity coefficients and to the choices of methods and models is investigated. In earlier DOME studies, three different bottom boundary conditions and a range of vertical grid sizes are applied. In other parts of the literature on numerical studies of oceanic gravity currents, there are statements that appear to contradict choices made on bottom boundary conditions in some of the DOME papers. In the present study, we therefore address the effects of the bottom boundary condition and vertical resolution in numerical investigations of dense water cascading on a slope. The main finding of the present paper is that it is feasible to capture the bottom Ekman layer dynamics adequately and cost efficiently by using a terrain-following model system using a quadratic drag law with a drag coefficient computed to give near-bottom velocity profiles in agreement with the logarithmic law of the wall. Many studies of dense water flows are performed with a quadratic bottom drag law and a constant drag coefficient. It is shown that when using this bottom boundary condition, Ekman drainage will not be adequately represented. In other studies of gravity flow, a no-slip bottom boundary condition is applied. With no-slip and a very fine resolution near the seabed, the solutions are essentially equal to the solutions obtained with a quadratic drag law and a drag coefficient computed to produce velocity profiles matching the logarithmic law of the wall. However, with coarser resolution near the seabed, there may be a substantial artificial blocking effect when using no-slip.     

Differences between magnitudes of Taiwan earthquakes from catalogs of Taiwan and Beijing

IntroductionTaiwanlocatedinthecollisionboundalbetweentheEurasianandthePhilippineSeaplatesisoneofhighseismicityregionsintheworld.HundredsofearthquakeswithM25occurredperyearandmorethan40withM27since1900.Amongtheseevents,shalloweventswithdepthofseveraltensofkilometersandintermediate-deepeventswithdepthof100-200kinexistwhichrepresentsacharacterofthesubductionzone.ThemagnitudesofTaiwaneventslistedinthecatalogofChineseearthquakesaretakenfromsomehistoricaldocumentsandGutenbergandRichter'sworks(19…     

Measurements of Reynolds stress in a wind-driven lagoonal estuary

Acoustic Doppler current profilers (ADCPs) have been used to measure Reynolds stresses in tidally dominated environments where wave action was minimal. In this paper, we examine observations from a microtidal estuary where the effects of wind stress and surface waves dominate the velocity variance. Reynolds stress measurements in this setting require a technique for addressing surface gravity wave contamination. We present here a method of reducing the effect of wave motion on Reynolds stresses by subtracting coincident observations along the axis of the ADCP beam. Linear wave theory is used to account for the attenuation of wave orbital velocities with depth. Using this method, Reynolds stress values are brought in line with those predicted by drag laws at the surface and bottom. The apparent Reynolds stress that is removed by the along-axis subtraction is shown to be largely due to the interaction of a slight tilt (1°) in the ADCP and the wave orbital velocity. During periods of stronger wind and waves, there is evidence of enhanced near-surface turbulence and momentum flux, presumably due to breaking waves. During these events, our calculated Reynolds stress magnitudes still appear reasonable, although the directions are suspect. We develop a diagnostic technique that clearly demarcates this region when it occurs. Coincident density profile measurements are used with the ADCP data to compute gradient Richardson numbers throughout the water column. Enhanced Reynolds stresses appear to correspond to Richardson numbers less than one. Responsible editor: Alejandro Souza     

The wave-induced boundary layer under long internal waves

The boundary layer formed under the footprint of an internal solitary wave is studied by numerical simulation for waves of depression in a two-layer model of the density stratification. The inviscid outer flow, in the perspective of boundary-layer theory, is based on an exact solution for the long wave-phase speed, yielding a family of fully nonlinear solitary wave solutions of the extended Korteweg–de Vries equation. The wave-induced boundary layer corresponding to this outer flow is then studied by means of simulation employing the Reynolds-averaged Navier–Stokes (RANS) formulation coupled with a turbulence closure model validated for wall-bounded flows. Boundary-layer characteristics are computed for an extensive range of environmental conditions and wave amplitudes. Boundary-layer transition, identified by monitoring the eddy viscosity, is correlated in terms of a boundary-layer Reynolds number. The frictional drag is evaluated for laminar, transitional, and turbulent cases, and correlations are presented for the friction coefficient plus relevant measures of the boundary-layer thickness.     

A Model Study of the Strong and Weak Wind,Stably Stratified Nocturnal Boundary Layer: Influence of Gentle Slopes

With the exception of intermittency and waves, a brief review of the observed and modeled mean structure of the nocturnal boundary layer (NBL) is presented. The effect of gentle slopes on strong and weak wind NBL was investigated here using a one-dimensional model, with a simple correction term to account for the slope effects, identical to the one used by Brost and Wyngaard (1978). The study indicates that the wind profiles, temperature profiles and surface layer turbulence characteristics are extremely sensitive to the imposed geostrophic wind when small slopes are present especially for light winds. This is due to the complex interaction between the buoyancy driven slope flow and the imposed geostrophic wind that in turn influence the shear generation of turbulence. Finally, the current issues in the modeling of weak wind boundary layer are discussed.     

Gibt es eine „kritische Windgeschwindigkeit“ für die Grenzfläche Wasser — Luft?

Zusammenfassung Aus Windprofilmessungen in feststehenden H?hen über den Wellenk?mmen lassen sich keine zuverl?ssigen Schlüsse über die effektive Schubkraft des Windes an der Meeresoberfl?che ziehen. Der Grund liegt im wesentlichen in der verschiedenartigen vertikalen Windgeschwindigkeitsverteilung über Wellenbergen und Wellent?lern und in der Tatsache, da? die Windprofilmessungen im wesentlichen nur das Windprofil über den Bergen widerspiegeln. Der Versuch, die Windprofilmessungen zu „reduzieren“, ist mit zu gro?en Fehlern behaftet, zeigt aber qualitativ eine systematische Erh?hung der aus Windmessungen berechneten Reibungsfaktoren an. Damit kann der Anschlu? an die aus Windstaubeobachtungen ermittelten Reibungsfaktoren gewonnen werden, so da? der vermeintliche „Sprung“ von einer „glatten“ zu einer „rauhen“ Grenzfl?che verschwindet, und eine „kritische Windgeschwindigkeit“, die diesen Sprung angeben soll, gar nicht existiert.

---

On the problem of the “critical wind velocity” at the air — sea boundary surface

Summary From the measurements of wind profiles at fixed heights above the wave crests no reliable conclusions can be drawn with regard to the wind's effective shearing force at the sea surface. This is essentially due to the diversified vertical distribution of the wind speed above the waves' crests and troughs as well as to the fact that when measuring wind profiles above the waves it is mainly the profiles above the wave crests that are obtained. Attempts to “reduce” wind profiles have proved to be subject to considerable errors; however, they show qualitatively a systematic increase of the frictional factors as computed from wind measurements. This links up with the frictional factors ascertained from observations of the wind effect so that the supposed “leap” from a “smooth” to a “rough” boundary surface and the “critical wind speed” corresponding to this leap do not exist at all.

     

Parameterization of gravity wave momentum deposition based on nonlinear wave interactions: basic formulation and sensitivity tests

We present a practical method for parameterization of gravity wave drag based on the Medvedev and Klaassen (1995. Journal of Geophysical Research 100, 25,841–25,853) theory of gravity wave spectral evolution and saturation. The only tuning necessary for the scheme involves assumptions about the nature of the source spectrum of subgrid-scale gravity waves, i.e. the wavenumbers, launch heights and amplitudes of the spectral components. In this paper we employ a column model with representative distributions of mean wind and temperature to examine the sensitivity of the parameterized wave drag to the source spectra. For the range of anticipated variability of source spectra in the troposphere the scheme produces plausible results consistent with observations and with theoretical estimates. Computationally, the scheme is as efficient as a multiple-wave Lindzen scheme, and suitable for use in general circulation models (GCMs) of the atmosphere.     

Gravity wave reflection: Case study based on rocket data

Since gravity waves significantly influence the atmosphere by transporting energy and momentum, it is important to study their wave spectrum and their energy dissipation rates. Besides that, knowledge about gravity wave sources and the propagation of the generated waves is essential. Originating in the lower atmosphere, gravity waves can move upwards; when the background wind field is equal to their phase speed a so-called critical layer is reached. Their breakdown and deposition of energy and momentum is possible. Another mechanism which can take place at critical layers is gravity wave reflection.In this paper, gravity waves which were observed by foil chaff measurements during the DYANA (DYnamics Adapted Network for the Atmosphere) campaign in 1990 in Biscarrosse (44°N, 1°W)—as reported by Wüst and Bittner [2006. Non-linear wave–wave interaction: case studies based on rocket data and first application to satellite data. Journal of Atmospheric and Solar-Terrestrial Physics 68, 959–976]—are investigated to look for gravity wave reflection processes. Following nonlinear theory, energy dissipation rates according to Weinstock [1980. Energy dissipation rates of turbulence in the stable free atmosphere. Journal of the Atmospheric Sciences 38, 880–883] are calculated from foil chaff cloud and falling sphere data and compared with the critical layer heights. Enhanced energy dissipation rates are found at those altitudes where the waves’ phase speed matches the zonal background wind speeds. Indication of gravity wave trapping is found between two altitudes of around 95 and 86 km.     

Simulation of Atmospheric Boundary Layer Characteristics during Indian Summer Monsoon using Observations from Monsoon Trough Boundary Layer Experiment at Jodhpur, India

The diurnal structure of the boundary layer during Indian summer monsoon period is studied using a one-dimensional meteorological boundary layer model and the observations collected from the Monsoon Trough Boundary Layer Experiment conducted in 1990 at Jodhpur, India. The model was initialized with the observed temperature profiles at 0530 LST on 17 July, 1990 at Jodhpur and was run for 26 hours. The study is carried out with a geostrophic wind speed of 9.5 m s⁻¹ corresponding to the strong wind simulation. The mean thermodynamic and wind structure simulated by the model are in good agreement with those observed from 30 m tower. The computed surface layer characteristics such as the surface fluxes, TKE and standard deviations of velocity components are found to be reasonably in good agreement with those based on turbulence measurements. The shear and buoyancy budget computed from the model are also compared with the turbulence measurements. The integrated cooling budget in the nocturnal boundary layer is examined.