光滑海面上风、温、湿廓线的层结订正模式及海气整体交换系数

详细介绍了光滑面标量粗糙度ＺＴ,Ｚｑ与风速粗糙度Ｚ０的相似表达式,论述了把Ｍｏｎｉｎ－Ｏｂｕｋｈｏｖ相似理论推广到光滑面上湍流气层的合理性,从而得到光滑面风、温、湿层结订正廓线与粗糙面廓线相统一的形式。总结了用该模式处理的三个海上梯度观测资料的计算结果,揭示出了微风时通过光滑海面的海气通量及整体交换系数受层结影响远大于风速影响的特征。并给出光滑界面上不同层结的整体交换系数随风速变化的拟合公式。     

Subfilter-scale Fluxes over a Surface Roughness Transition. Part I: Measured Fluxes and Energy Transfer Rates

A wind-tunnel experiment was designed and carried out to study the effect of a surface roughness transition on subfilter-scale (SFS) physics in a turbulent boundary layer. Specifically, subfilter-scale stresses are evaluated that require parameterizations and are key to improving the accuracy of large-eddy simulations of the atmospheric boundary layer. The surface transition considered in this study consists of a sharp change from a rough, wire-mesh covered surface to a smooth surface. The resulting magnitude jump in aerodynamic roughnesses, M = ln(z ₀₁/z ₀₂), where z ₀₁ and z ₀₂ are the upwind and downwind aerodynamic surface roughnesses respectively, is similar to that of past experimental studies in the atmospheric boundary layer. The two-dimensional velocity fields used in this study are measured using particle image velocimetry and are acquired at several positions downwind of the roughness transition as well as over a homogeneous smooth surface. Results show that the SFS stress, resolved strain rate and SFS transfer rate of resolved kinetic energy are dependent on the position within the boundary layer relative to the surface roughness transition. A mismatch is found in the downwind trend of the SFS stress and resolved strain rate with distance from the transition. This difference of behaviour may not be captured by some eddy-viscosity type models that parameterize the SFS stress tensor as proportional to the resolved strain rate tensor. These results can be used as a benchmark to test the ability of existing and new SFS models to capture the spatial variability SFS physics associated with surface roughness heterogeneities.     

On theoretical wind speed and temperature profiles over the sea with applications to data from sable island,Nova Scotia

Abstract

In October 1985, the Boundary‐Layer Research Division of the Atmospheric Environment Service conducted an experiment on Sable Island, Nova Scotia, where 10‐m wind measurements were made at a number of locations. Wind data were also collected at 4 levels on one of the 10‐m masts and at 6 levels on a 26‐m mast, both located on the South Beach. Other data used in the present study consisted of air temperature measured at 9 m and sea temperature measured at the beach.

The theory for wind speed and temperature profiles over the sea is reviewed. A method of deriving over‐sea profile parameters (u_*, θ_*, Z₀, L) from wind data at one level and the air‐sea temperature difference is described. The method is limited to applications either over homogeneous open ocean or, provided measurements are taken above the internal boundary layer generated by the change of roughness at the coastline, over a flat beach (without coastal orography). The heights at which the method is applied must be within the surface layer which must not have any discontinuities in wind speed or temperature in the vertical, such as are often associated with inversion layers. An application to data collected at beach sites in onshore flow during the October 1985 experiment is illustrated.

Once the above parameters are obtained, theoretical wind profiles may be computed and compared with observed profiles. In order to make a proper comparison it is essential to account for internal boundary layers generated at the shoreline by the step‐change in surface roughness. Only the data measured above the internal boundary layer are representative of over‐sea conditions and may, therefore, be used for verifying the theoretical profiles. The agreement between calculated and measured data is generally very good. One complication, however, is a slight upstream‐blockage effect due to a 7‐m high dune located about 140 m downwind of the 26‐m mast. Estimates of the magnitude of this effect partially account for small discrepancies in the results at the 26‐ and 10‐m mast locations.

An estimation of the most probable errors in the calculated parameters, based on assumed measurement errors, is included in the computer program. Results suggest that small measurement errors can explain the above discrepancies.     

Note on a Parametric Relation for Separating Flow over a Rough Hill

A parametrization method used to account for the effects of flow separation and wall roughness on the lower boundary condition for turbulent boundary layers is investigated against direct numerical simulation and laser Doppler anemometry data. The numerical simulation represents flow over a smooth, flat surface with a prescribed external adverse pressure gradient. The water-channel experiments cover flow over smooth and rough hills for two specified Reynolds numbers. Global optimization algorithms based on four different direct search methods are used to assess the parametrization function, C, in terms of local mean velocity profiles and the parametrization parameters u _* (friction velocity), ∂ _x p (local pressure gradient), z ₀ (effective roughness) and d (zero-plane displacement). The study investigates regions of attached and reversed flows, and forty-two velocity profiles are compared with the proposed expression for the function C, including two profiles that satisfy the solution of Stratford.     

The flow in a turbulent boundary layer upstream of a change in surface roughness

Modification of a turbulent flow upstream of a change in surface roughness has been studied by means of a stream function-vorticity model.A flow reduction is found upstream of a step change in surface roughness when a fluid flows from a smooth onto a rough surface. Above that layer and above the region of flow reduction downstream of a smooth-rough transition, a flow acceleration is observed. Similar flow modification can be seen at a rough-smooth transition with the exception that flow reduction and flow acceleration are reversed. Within a fetch of –500 < x/z ₀< + 500 (z ₀ is the maximum roughness length, the roughness transition is located at x/z ₀ = 0), flow reduction (flow acceleration) upstream of a roughness transition is one order of magnitude smaller than the flow reduction (flow acceleration) downstream of a smooth-rough (rough-smooth) transition. The flow acceleration (flow reduction) above that layer is two orders of magnitude.The internal boundary layer (IBL) for horizontal mean velocity extends to roughly 300z ₀ upstream of a roughness transition, whereas the IBL for turbulent shear stress as well as the distortion of flow equilibrium extend almost twice as far. For the friction velocity, an undershooting (overshooting) with respect to upstream equilibrium is predicted which precedes overshooting (undershooting) over new equilibrium just behind a roughness transition.The flow modification over a finite fetch of modified roughness is weaker than over a corresponding fetch downstream of a single step change in roughness and the flow stays closer to upstream equilibrium. Even in front of the first roughness change of a finite fetch of modified roughness, a distortion of flow equilibrium due to the second, downwind roughness change can be observed.     

Measurements of spray at an air-water interface

The spray content in the surface boundary layer above an air—water interface was determined by a series of measurements at various feteches and wind speeds in a laboratory facility. The droplet flux density N(z) can be described in terms of the scaling flux density N_* and von Karman constant K throguh the equation, N(z)/N_* = −(1/K) ln(z/z_0d) where z is height above the mean water level and z_0d is the droplet boundary layer thickness. N_* is given by a unique relationship in terms of the roughness Reynolds number u_*σ/ν where σ is the root-mean-square surface displacement. Spray inception occurred for u_* 0.3. The dominant mode of spray generation in the present and most other laboratory tests, as well as in available field data, appears to be bubble bursting.     

Regional roughness of the landes forest and surface shear stress under neutral conditions

Mean wind velocity profiles were measured by means of radio-windsondes over the Landes region in southwestern France, which consists primarily of pine forests with scattered villages and clearings with various crops. Analysis of neutral profiles indicated the existence of a logarithmic layer between approximately z – d ₀ = 67(±18)z ₀ and 128(+-32)z ₀ (z is the height above the ground, z ₀ the surface roughness and d ₀ the displacement height). The upper limit can also be given as z – d ₀ = 0.33 (±0.18)h, where h is the height of the bottom of the inversion. The profiles showed that the surface roughness of this terrain is around 1.2 m and the displacement height 6.0 m. Shear stresses derived from the profiles were in good agreement with those obtained just above the forest canopy at a nearby location with the eddy correlation method by a team from the Institute of Hydrology (Wallingford, England).     

Wind structure near a city centre

Wind speed and temperature profiles to a height of 8 m were recorded for 30-, 60-, and 90-min averaging times over a striated snow surface at the geographic South Pole during the austral winter of 1975. A gradient Richardson number was calculated for each averaging time to determine conditions of neutral stability under which the logarithmic wind law would hold. A log-linear regression technique was used to determine values of aerodynamic roughness height (Z ₀) for those profile averages recorded in conditions of neutral stability. A plot of Z ₀ as a function of average wind direction revealed a variation in Z ₀ of almost three orders of magnitude, from 0.01 to 7 cm, over 120 deg of wind direction. A simple model is presented to justify the fact that aerodynamic roughness is a function of wind direction and erosion history.     

A simple unified theory for flow in the canopy and roughness sublayer

The mean flow profile within and above a tall canopy is well known to violate the standard boundary-layer flux–gradient relationships. Here we present a theory for the flow profile that is comprised of a canopy model coupled to a modified surface-layer model. The coupling between the two components and the modifications to the surface-layer profiles are formulated through the mixing layer analogy for the flow at a canopy top. This analogy provides an additional length scale—the vorticity thickness—upon which the flow just above the canopy, within the so-called roughness sublayer, depends. A natural form for the vertical profiles within the roughness sublayer follows that overcomes problems with many earlier forms in the literature. Predictions of the mean flow profiles are shown to match observations over a range of canopy types and stabilities. The unified theory predicts that key parameters, such as the displacement height and roughness length, have a significant dependence on the boundary-layer stability. Assuming one of these parameters a priori leads to the incorrect variation with stability of the others and incorrect predictions of the mean wind speed profile. The roughness sublayer has a greater impact on the mean wind speed in stable than unstable conditions. The presence of a roughness sublayer also allows the surface to exert a greater drag on the boundary layer for an equivalent value of the near-surface wind speed than would otherwise occur. This characteristic would alter predictions of the evolution of the boundary layer and surface states if included within numerical weather prediction models.     

Flux–Profile Relationships Over a Fetch Limited Beech Forest

The influence of an internal boundary layer and a roughness sublayer on flux–profile relationships for momentum and sensible heat have been investigated for a closed beech forest canopy with limited fetch conditions. The influence was quantified by derivation of local scaling functions for sensible heat flux and momentum (_h and _m) and analysed as a function of atmospheric stability and fetch. For heat, the influences of the roughness sublayer and the internal boundary layer were in agreement with previous studies. For momentum, the strong vertical gradient of the flow just above the canopy top for some wind sectors led to an increase in _m, a feature that has not previously been observed. For a fetch of 500 m over the beech forest during neutral atmospheric conditions, there is no height range at the site where profiles can be expected to be logarithmic with respect to the local surface. The different influence of the roughness sublayer on _h and _m is reflected in the aerodynamic resistance for the site. The aerodynamic resistance for sensible heat is considerably smaller than the corresponding value for momentum.     

Total energy equation method for calculating hurricane intensity

Summary Surface pressure reduction in hurricanes is calculated by applying the total energy equation (TEE) to ideal isentropic upflow in a vertical tube. The pressure reduction at the base of the tube, called the intensity, is calculated for three upflow processes: reversible upflow of air approaching equilibrium with the sea at the sea level pressure outside the tube; irreversible upflow of air approaching equilibrium with the sea at the sea-level pressure outside the tube; and upflow of air approaching equilibrium with the sea at the reduced surface pressure inside the tube. The sensitivity of intensity to the type of upflow process and to the sea surface temperature is investigated. Intensities calculated with the TEE are shown to be consistent with observations and to be close to intensities calculated with more complex methods. The TEE method is simple and can help understand the basic mechanism responsible for surface pressure reduction and for energy production. The method is used to show that approximately 20% of the heat taken from the sea during a hurricane is converted to mechanical energy. Received December 4, 2000/Revised March 29, 2001     

A theory for the scalar roughness and the scalar transfer coefficients over snow and sea ice

Although the bulk aerodynamic transfer coefficients for sensible (C _H ) and latent (C _E ) heat over snow and sea ice surfaces are necessary for accurately modeling the surface energy budget, they have been measured rarely. This paper, therefore, presents a theoretical model that predicts neutral-stability values of C _H and C _E as functions of the wind speed and a surface roughness parameter. The crux of the model is establishing the interfacial sublayer profiles of the scalars, temperature and water vapor, over aerodynamically smooth and rough surfaces on the basis of a surface-renewal model in which turbulent eddies continually scour the surface, transferring scalar contaminants across the interface by molecular diffusion. Matching these interfacial sublayer profiles with the semi-logarithmic inertial sublayer profiles yields the roughness lengths for temperature and water vapor. When coupled with a model for the drag coefficient over snow and sea ice based on actual measurements, these roughness lengths lead to the transfer coefficients. C _E is always a few percent larger than CH. Both decrease monotonically with increasing wind speed for speeds above 1 m s^–1, and both increase at all wind speeds as the surface gets rougher. Both, nevertheless, are almost always between 1.0 × 10^–3 and 1.5 × 10^–3.     

A study of the roughness effects of multiple windbreaks

A wind-tunnel study of the roughness effect of a network of uniformly spaced long parallel windbreaks has been attempted. The roughness parameterZ ₀ in the well-known logarithmic profile has been obtained as a function of the ratio of the spacingL to heightH of the barriers. In one part of the study, a naturally developing boundary layer was used and the windbreaks were not spread on the entire wind-tunnel floor. Spacing ratiosL/H of 1 to 20 were treated. In the other part of the study, spires were used to thicken the boundary layers artificially and the windbreaks were spread throughout the wind-tunnel floor. In this case, spacing ratiosL/H of 1 to 30 were examined.Values of the roughness parameter obtained from the above two cases have been compared to some (though very few) previously published pertinent values.     

Incorporation of a high-roughness lower boundary into a mesoscale model for studies of dry deposition over complex terrain

For flow over natural surfaces, there exists a roughness sublayer within the atmospheric surface layer near the boundary. In this sublayer (typically 50z ₀ deep in unstable conditions), the Monin-Obukhov (M-O) flux profile relations for homogeneous surfaces cannot be applied. We have incorporated a modified form of the M-O stability functions (Garratt, 1978, 1980, 1983) in a mesoscale model to take account of this roughness sublayer and examined the diurnal variation of the boundary-layer wind and temperature profiles with and without these modifications. We have also investigated the effect of the modified M-O functions on the aerodynamic and laminar-sublayer resistances associated with the transfer of trace gases to vegetation. Our results show that when an observation height or the lowest level in a model is within the roughness sublayer, neglect of the flux-profile modifications leads to an underestimate of resistances by 7% at the most.     

Sensitivity of PBL model predictions to model design and uncertainties in environmental inputs

A simple time-dependent one-dimensional model of the planetary boundary layer (PBL) is described and used to examine the degree to which model design decisions affect model output variables. The model's sensitivity to changes in the environmental conditions is also explored. Averages of the surface fluxes, near-ground wind speeds and other PBL properties from 48 h simulations are compared to control runs. The model-calculated surface fluxes are most sensitive, in decreasing order of importance, to the vertical grid spacing, the form of closure between the surface temperature and the atmosphere, the use of vertical diffusivity smoothing, the choice of maximum time step and choice of turbulence closure scheme. These fluxes are relatively insensitive to mixing-length scaling or choice of implicit time step weighting factor. Sensitivity to changes in soil type exceeds any of the design criteria tested. The modeled fluxes are moderately sensitive to small variations in the horizontal pressure gradient, to unsteadiness in the geostrophic wind and to variations in surface roughness. They are relatively insensitive to uncertainties in local vertical velocities and small (25%) variations applied separately to soil thermal diffusivity or heat capacity. The sensitivity of the average PBL depth (Z _i ) to model and environmental changes are similar to those of surface fluxes except thatZ _i is more sensitive to changes in mixing length, albedo and imposed vertical velocity then are the surface fluxes.     

Mixed-layer heat advection and entrainment during the sea breeze

A model is developed to simulate the potential temperature and the height of the mixed layer under advection conditions. It includes analytic expressions for the effects of mixed-layer conditions upwind of the interface between two different surfaces on the development of the mixed layer downwind from the interface. Model performance is evaluated against tethersonde data obtained on two summer days during sea breeze flow in Vancouver, Canada. It is found that the mixed-layer height and temperature over the ocean has a small but noticeable effect on the development of the mixed layer observed 10 km inland from the coast. For these two clear days, the subsidence velocity at the inversion base capping the mixed layer is estimated to be about 30 mm s^–1 from late morning to late afternoon. When the effects of subsidence are included in the model, the mixed-layer height is considerably underpredicted, while the prediction for the mean potential temperature in the mixed layer is considerably improved. Good predictions for both height and temperature can be obtained when values for the heat entrainment ratio,c, 0.44 and 0.68 for these two days respectively for the period from 1000 to 1300 LAT, were used. These values are estimated using an equation including the additional effects on heat entrainment due to the mechanical mixing caused by wind shear at the top of the mixed layer and surface friction. The contribution of wind shear to entrainment was equal to, or greater than, that from buoyant convection resulting from the surface heat flux. Strong wind shear occurred near the top of the mixed layer between the lower level inland flow and the return flow aloft in the sea breeze circulation.Symbols c entrainment parameter for sensible heat - c _p specific heat of air at constant pressure, 1010 J kg^–1 K^–1 - d ₁ the thickness of velocity shear at the mixed-layer top, m - Q _H surface sensible heat flux, W m^–2 - u _m mean mixed-layer wind speed, m s^–1 - u _* friction velocity at the surface, m s^–1 - w subsidence velocity, m s^–1 - W subsidence warming,^oC s^–1 - w _e entrainment velocity, m s^–1 - w _* convection velocity in the mixed layer, m s^–1 - x downwind horizontal distance from the water-land interface, m - y dummy variable forx, m - Z height above the surface, m - Z _i height of capping inversion, m - Z _m mixed-layer depth, i.e.,Z _i–Z_s, m - Z _s height of the surface layer, m - lapse rate of potential temperature aboveZ _i, K m^–1 - potential temperature step atZ _i, K - u _h velocity step change at the mixed-layer top - _m mean mixed-layer potential temperature, K     

Formulation of the thermal internal boundary layer in a mesoscale model

Mesoscale models using a non-local K-scheme for parameterization of boundary-layer processes require an estimate of the planetary boundary layer (PBL) height z _i at all times. In this paper, two-dimensional sea-breeze experiments are carried out to evaluate three different formulations for the advective contribution in the z _i prognostic equation of Deardorff (1974).Poor representation of the thermal internal boundary layer in the sea breeze is obtained when z _i is advected by the wind at level z _i . However, significantly better results are produced if the mean PBL wind is used for the advecting velocity, or if z _i is determined simply by checking for the first sufficiently stable layer above the ground.A Lagrangian particle model is used to demonstrate the effect of each formulation on plume dispersion by the sea breeze.     

The influence of atmospheric stratification on the growth of water waves

The atmospheric surface layer over sea has a density stratification which varies with moisture content and air/sea temperature difference. This influences the growth of water waves. To study the effect quantitatively, the Reynolds equations are solved numerically. For given wind speed and surface roughness, wave growth is found to be more rapid in unstably stratified conditions than in stable conditions. This is due to an increase in turbulence, primarily caused by an increase of mixing length.Under the assumption of a Charnock relation between surface roughness and friction velocity, it is found that for large inverse wave age (u _*/c>0.07), the effect of stratification on wave growth is weell described by Monin-Obukhov scaling of the friction velocity. For smaller values ofu _*/c, Monin-Obukhov scaling overpredicts.The effect on duration-limited wave growth is studied with the third-generation WAM surface wave model driven by 10 m winds. Effects of stratification on the significant wave height are found to be of the order of 10%. The results are comparable to those of a recent reanalysis of field measurements, although the measured stratification effect is somewhat stronger. Implementation of a stratification-dependent growth in wave models is recommended, as it can lead to small but significant improvements in wave forecasts when accurate air and sea temperatures are available.     

基于S波段双偏振雷达和X波段相控阵雷达的超级单体观测分析

利用广州S波段双偏振雷达和X波段相控阵雷达资料,对2022年3月26日一次降雹超级单体风暴成熟阶段的雷达观测特征开展分析,结果表明:超级单体呈现出钩状回波、回波悬垂、中气旋、三体散射等经典结构特征。径向速度上观测到中低层辐合、高层辐散以及中气旋和反气旋共存的双涡旋结构,有助于超级单体的维持发展。偏振特征分析发现,超级单体低层出现了反射率因子(ZDR)弧,低层强回波区对应偏小的差分反射率(ZDR)、低的相关系数(CC)和大的差分相移率(KDP),符合融化的冰雹特征。中层观测到ZDR环、CC环和三体散射(TBSS)的偏振特征。高层强回波区对应低的ZDR、较高的CC和低的KDP,对应空中干的大冰雹。垂直方向上观测到ZDR柱和KDP柱,ZDR柱最大发展高度达到8 km。X波段相控阵雷达更快的扫描速度还精细监测到超级单体钩状回波和中气旋的形成演变过程,低层也观测到与S波段双偏振雷达类似的ZDR弧特征和融化中的冰雹特征,但是使用中要留意衰减造成的影响。     

Near-Wall Flow Development After A Step Change In Surface Roughness

An experimental study of the initial flow field downstream of a step change in surface roughness is presented. The roughness length of the downstream surface was approximately tenfold that of the upstream roughness and, unlike all previous studies, attention was concentrated on the roughness sublayer region beneath the inertial (log-law) region. The experiments were conducted at a boundary layer Reynolds number of about 6 × 10⁴ (based on layer thickness andfree-stream velocity) and around a longitudinal location where the (downstream) roughness length, z_o2, was about 1% of the boundary-layer thickness atthe roughness change point.The thickness of the roughness sublayer was found for the two roughness. It was observed that the vertical profiles of mean velocity and turbulence characteristics started to show similarity after about 160z₀₂ downstream of the roughness change. The presence of a shear stress overshoot is shown to depend strongly on the precise location (with respect to the roughness elements) at which the measurements are made and the thickness of the equilibrium layer is shown to be very sensitive to the way it is defined. It is demonstrated that the growing equilibrium layer has first to encompass the roughness sublayer before any thickness of inertial sublayer can be developed. It follows that, in somepractical cases, like flows across some urban environments, the latter(log-law) region may never exist at all.