Some observations of turbulence and turbulent transport within and above plant canopies

Observations have been made of the structure of turbulence and turbulent exchange within plant canopy layers. A new three-dimensional anemometer was used to measure the eddy fluxes of heat and momentum, and the related cospectra, within and above a corn crop and above a red pine forest. Measured values of momentum and heat fluxes, at each height within the corn canopy, were relatively constant proportions of the flux above the canopy, for the period of a day's observation. Extensive regions obeying a –5/3 power relation were found. Isotropy was found above the forest at high frequencies while above and within the corn crop, the ratios of the lateral and vertical spectral densities to the longitudinal component were less than the expected value in the – 5/3 region. In all situations, the vertical velocity spectra were more peaked than a universal curve, particularly a vertical velocity spectrum from above the forest. It is suggested that the additional variance results from the mixing caused by the individual roughness elements. As expected, the spectra could not be normalized using the height above the soil surface to calculate a non-dimensional frequency, but scaling heights were estimated by matching the frequencies of the peak of each curve with that of the universal curve. Cospectra of uw and wT within the corn canopy were of similar shape and frequency regime, and were basically similar in shape to cospectra above the crop. All of the cospectra were more sharply peaked than universal cospectral curves.     

A Second-Order Closure for Neutrally Stratified Vegetative Canopy Flows

An existing second-order closure model is modified to include the effects on mean and turbulent motions of form and viscous drag in vegetative canopies. The additional physical mechanisms represented by the closure are viscous and pressure drag on canopy elements, their role in momentum absorption, in the creation of fine scale turbulent eddies and in enhancing the total viscous dissipation in the canopy airspace. Viscous dissipation is split into a standard 'isotropic contribution associated with the spectral eddy cascade and a foliage contribution associated with work against pressure and viscous drag on the foliage. Changes in the turbulent time scale that result from these mechanisms are included in the standard parameterisations of third moments and of the eddy cascade contribution to dissipation. The model is tested against a wind- tunnel 'wheat canopy, a corn canopy and a eucalypt canopy, a height range from 50 mm to 12.6 m. Model results show that the parameterisations of foliage interaction used in the closure are sufficiently robust to reproduce second-moment profiles within and above vegetative canopies to a high degree of accuracy without resorting to 'tuning of the model constants. The model also shows the natural emergence of two length scales, one associated with the familiar eddy cascade isotropic contribution to total dissipation and the other associated with the length scales of the canopy elements.     

The conductance of a maize crop and the underlying soil to ozone under various environmental conditions

Flux measurements of ozone and water vapour employing the eddy correlation technique were used to determine the surface conductance and canopy conductance to ozone. In the surface conductance to ozone, all surfaces at which ozone is destroyed and the transport process to these surfaces are included. The canopy conductance to ozone represents the ozone uptake of transpiring plant parts. The surface conductance to ozone of the maize crop and the underlying soil was generally larger than the canopy conductance to ozone. This means that beside the uptake by stomata, there was another important ozone sink. Under wet soil surface conditions, the surface conductance and the canopy conductance to ozone coincided. This indicates that the resistance of wet soil and the remaining plant parts (cuticle) to ozone was much larger than the stomatal or soil resistance. On the other hand, under dry soil conditions the conductances differ, largely caused by a variation in the transport process to the soil. The transport of ozone to soil increased with increasing friction velocity (u _*) and decreased with increasing atmospheric stability, leaf area index (LAI) or crop height (h). These effects for midday (unstable) conditions were parameterized with an in-crop aerodynamic resistance,r _inc in a very straightforward way;r _inc=13.9 LAIh/u _*+67 (cc.=0.77). If the ozone flux in air pollution models is described with a simple resistance model (Big Leaf model), the extra destruction at the soil should be modelled using an in-crop aerodynamic resistance. For these measurements the ozone flux to the soil was 0–65% of the total ozone flux measured above the crop. Under wet soil conditions, this was less than 20%; under dry soil conditions, this was 30–65%.     

Characteristics of canopy turbulence over a deciduous forest on complex terrain

Canopy turbulence plays an important role in mass and energy exchanges at the canopy-atmosphere interface. Despite extensive studies on canopy turbulence over a flat terrain, less attention has been given to canopy turbulence in a complex terrain. The purpose of this study is to scrutinize characteristics of canopy turbulence in roughness sublayer over a hilly forest terrain. We investigated basic turbulence statistics, conditionally sampled statistics, and turbulence spectrum in terms of different atmospheric stabilities, wind direction and vertical structures of momentum fluxes. Similarly to canopy turbulence over a homogeneous terrain, turbulence statistics showed coherent structure. Both quadrant and spectrum analysis corroborated the role of intermittent and energetic eddies with length scale of the order of canopy height, regardless of wind direction except for shift of peak in vertical wind spectrum to relatively high frequency in the down-valley wind. However, the magnitude of the momentum correlation coefficient in a neutral condition was smaller than typical value over a flat terrain. Further scrutiny manifested that, in the up-valley flow, temperature skewness was larger and the contribution of ejection to both momentum and heat fluxes was larger compared to the downvalley flow, indicating that thermal instability and weaker wind shear in up-valley flow asymmetrically affect turbulent transport within the canopy.     

Eddy correlation measurements of CO2, latent heat,and sensible heat fluxes over a crop surface

Eddy correlation equipment was used to measure mass and energy fluxes over a soybean crop. A rapid response CO₂ sensor, a drag anemometer, a Lyman-alpha hygrometer and a fine wire thermocouple were used to sense the fluctuating quantities.Diurnal fluxes of sensible heat, latent heat and CO₂ were calculated from these data. Energy budget closure was obtained by summing the sensible and latent heat fluxes determined by eddy correlation which balanced the sum of net radiation and soil heat flux. Peak daytime CO₂ fluxes were near 1.0 mg m^–2 (ground area) s^–1.The eddy correlation technique was also employed in this study to measure nocturnal CO₂ fluxes caused by respiration from plants, soil, and roots. These CO₂ fluxes ranged from - 0.1 to - 0.25 mg m^–2s^–1.From the data collected over mature soybeans, a relationship between CO₂ flux and photosynthetically active radiation (PAR) was developed. The crop did not appear to be light-saturated at PAR flux densities < 1800 Ei m^–2 s^–1. The light compensation point was found to be about 160 Ei m^–2 s^–1.Published as Paper No. 7402, Journal Series, Nebraska Agricultural Experiment Station. The work reported here was conducted under Nebraska Agricultural Experiment Station Project 27-003 and Regional Research Project 11–33.Post-doctoral Research Associate, Professor and Professor, respectively. Center for Agricultural Meteorology and Climatology, Institute of Agriculture and Natural Resources, University of Nebraska, Lincoln, NE 68583-0728.     

Influence of Nocturnal Low-level Jets on Eddy-covariance Fluxes over a Tall Forest Canopy

Observations of low-level jets (LLJs) at the Howland AmeriFlux site in the USA and the jet’s impact on nocturnal turbulent exchange and scalar fluxes over a tall forest canopy are discussed. Low-frequency motions and turbulent bursts characterize moderately strong LLJs, whereas low-frequency motions are suppressed during periods with strong LLJs and enhanced shear. An analysis based on the shear-sheltering hypothesis seeks to elucidate the effect of LLJs on flux measurements. In the absence of shear sheltering, large eddies penetrate the roughness sublayer causing enhanced mixing while during periods with shear sheltering, mixing is reduced. In the absence of the latter, ‘upside-down’ eddies are primarily responsible for the enhanced velocity variances, scalar and momentum fluxes. The integral length scales over the canopy are greater than the canopy height. The variance spectra and cospectra from the wavelet analysis indicate that large eddies (spatial scale greater than the low-level jet height) interact with active canopy-scale turbulence, contributing to counter-gradient scalar fluxes.     

Heat and momentum transport in an oak forest canopy

With the use of a sonic anemometer, vertical heat and momentum fluxes were measured at three different levels in an oak forest canopy. A quadrant analysis of the resulting data shows that approximately half of the transport occurs in extreme events lasting about 5 to 10% of the time. The partition of transport into momentum sweeps, bursts and interactions shows good agreement with existing data. The heat flux is analysed by observing the fluxes during the different momentum events and considering concurrent momentum and heat flux intensities by means of conditional probabilities. While low intensity (normal) events show similar probability distributions throughout the canopy, different structures appear at the three measurement heights for high intensity (extreme) events that can tentatively be explained by taking the temperature profile into account. This supports the idea that these events are coherent motions with scales comparable to the canopy height.     

The Validity of Similarity Theory in the Roughness Sublayer Above Forests

Flux-gradient relationships based upon similarity theory have been reported to severely underestimate scalar fluxes in the roughness sublayer above forests, as compared to independent flux estimates (for example, eddy covariance or energy balance measurements). This paper presents the results of a unique three-month investigation into the validity of similarity theory in the roughness sublayer above forests. Eddy covariance and flux-gradient measurements of carbon dioxide (CO²) exchange were compared above a mixed deciduous forest at Camp Borden, Ontario, both before and after leaf senescence. The eddy covariance measurements used a Li-Cor infrared gas analyzer, and the flux-gradient (similarity theory) measurements featured a tunable diode laser Trace Gas Analysis System (TGAS). The TGAS resolved the CO² concentration difference to 300 parts per trillion by volume (ppt) based upon a half-hour sampling period. The measured enhancement factor (the ratio of independent flux estimates, in this case eddy covariance, to similarity theory fluxes) was smaller and occurred closer to the canopy than in most previous investigations of similarity theory. Very good agreement between the eddy covariance and similarity theory fluxes was found between 1.9 and 2.2 canopy heights (h^c), and the mean enhancement factors measured before and after leaf senescence were 1.10 plusmn; 0.06 and 1.24 ± 0.07, respectively. Larger discrepancies were measured closer to the canopy (1.2 to 1.4 h^c), and mean enhancement factors of 1.60 ± 0.10 and 1.82 ± 0.11 were measured before and after leaf senescence, respectively. Overall, the Borden results suggest that similarity theory can be used within the roughness sublayer with a greater confidence than previously has been believed.     

Estimation of Scalar Source/Sink Distributions in Plant Canopies Using Lagrangian Dispersion Analysis: Corrections for Atmospheric Stability and Comparison with a Multilayer Canopy Model

Source/sink distributions of heat, water vapour andCO₂ within a rice canopy were inferred using aninverse Lagrangian dispersion analysis and measuredmean profiles of temperature, specific humidity andCO₂ mixing ratio. Monin–Obukhov similarity theorywas used to account for the effects of atmosphericstability on _w(z), the standard deviation ofvertical velocity and _L(z), the Lagrangian timescale of the turbulence. Classical surface layer scaling was applied in the inertial sublayer (z > z_ruf)using the similarity parameter = (z - d)/L, where z is height above ground, d is the zero plane displacementheight for momentum, L is the Obukhov length,and z_ruf 2.3h_c, where h_c iscanopy height. A single length scale h_c, was usedfor the stability parameter 3 = h_c/L in the height range 0.25 < z/h_c < 2.5. This choice is justified by mixing layer theory, which shows that within the roughness sublayer there is one dominant turbulence length scaledetermined by the degree of inflection in the windprofile at the canopy top. In the absence of theoretical or experimental evidence for guidance,standard Monin–Obukhov similarity functions, with = h_c/L, were used to calculate the stabilitydependence of _w(z) and _L(z) in the roughness sublayer. For z/h_c < 0.25 the turbulence length and time scales are influenced by the presence of the lowersurface, and stability effects are minimal. With theseassumptions there was excellent agreement between eddycovariance flux measurements and deductions from theinverse Lagrangian analysis. Stability correctionswere particularly necessary for night time fluxes whenthe atmosphere was stably stratified.The inverse Lagrangian analysis provides a useful toolfor testing and refining multilayer canopy models usedto predict radiation absorption, energy partitioningand CO₂ exchanges within the canopy and at thesoil surface. Comparison of model predictions withsource strengths deduced from the inverse analysisgave good results. Observed discrepancies may be dueto incorrect specification of the turbulent timescales and vertical velocity fluctuations close to theground. Further investigation of turbulencecharacteristics within plant canopies is required toresolve these issues.     

Momentum- and Heat-Flux Parametrization at Dome C,Antarctica: A Sensitivity Study

An extensive meteorological observational dataset at Dome C, East Antarctic Plateau, enabled estimation of the sensitivity of surface momentum and sensible heat fluxes to aerodynamic roughness length and atmospheric stability in this region. Our study reveals that (1) because of the preferential orientation of snow micro-reliefs (sastrugi), the aerodynamic roughness length \(z_{0}\) varies by more than two orders of magnitude depending on the wind direction; consequently, estimating the turbulent fluxes with a realistic but constant \(z_{0}\) of 1 mm leads to a mean friction velocity bias of \(24\,\%\) in near-neutral conditions; (2) the dependence of the ratio of the roughness length for heat \(z_{0t}\) to \(z_{0}\) on the roughness Reynolds number is shown to be in reasonable agreement with previous models; (3) the wide range of atmospheric stability at Dome C makes the flux very sensitive to the choice of the stability functions; stability function models presumed to be suitable for stable conditions were evaluated and shown to generally underestimate the dimensionless vertical temperature gradient; as these models differ increasingly with increases in the stability parameter z / L, heat flux and friction velocity relative differences reached \(100\,\%\) when \(z/L > 1\); (4) the shallowness of the stable boundary layer is responsible for significant sensitivity to the height of the observed temperature and wind data used to estimate the fluxes. Consistent flux results were obtained with atmospheric measurements at heights up to 2 m. Our sensitivity study revealed the need to include a dynamical parametrization of roughness length over Antarctica in climate models and to develop new parametrizations of the surface fluxes in very stable conditions, accounting, for instance, for the divergence in both radiative and turbulent fluxes in the first few metres of the boundary layer.     

Surface Heterogeneity Effects on Regional-Scale Fluxes in the Stable Boundary Layer: Aerodynamic Roughness Length Transitions

The effects of abrupt streamwise transitions of the aerodynamic roughness length ( $z_\mathrm{o}$ z o ) on the stable atmospheric boundary layer are evaluated using a series of large-eddy simulations based on the first Global Energy and Water Cycle Experiment Atmospheric Boundary Layer intercomparison study (GABLS1). Four $z_\mathrm{o}$ z o values spanning three orders of magnitude are used to create all possible binary distributions with each arranged into patches of characteristic length scales equal to roughly one-half, one, and two times the equivalent homogeneous boundary-layer height. The impact of the heterogeneity on mean profiles of wind speed and temperature, on surface fluxes of heat and momentum, and on internal boundary-layer dynamics are considered. It is found that $z_\mathrm{o}$ z o transitions do not significantly alter the functional relationship between the average surface fluxes and the mean profiles of wind speed and potential temperature. Although this suggests that bulk similarity theory is applicable for modelling the stable boundary layer over $z_\mathrm{o}$ z o heterogeneity, effective surface parameters must still be specified. Existing models that solve for effective roughness lengths of momentum and heat are evaluated and compared to values derived from the simulation data. The existing models are unable to accurately reproduce both the values of the effective aerodynamic roughness lengths and their trends as functions of patch length scale and stability. A new model for the effective aerodynamic roughness length is developed to exploit the benefits of the other models tested. It accurately accounts for the effects of the heterogeneity and stratification on the blending height and effective aerodynamic roughness length. The new model provides improved average surface fluxes when used with bulk similarity.     

Influence of foliar density and thermal stability on profiles of Reynolds stress and turbulence intensity in a deciduous forest

Observations were made of turbulence in an extensive deciduous forest on level terrain using a vertical array of seven three-dimensional sonic anemometer/thermometers within and above the canopy. Data were collected through the period of leaf fall and over a range of thermal stabilities. A bulk canopy drag coefficient was nearly independent of the density of the forest but decreased greatly with the onset of nocturnal stability. The depth of penetration of momentum into the forest increased with leaf fall but, again, was greatly curtailed by stable conditions. Turbulent velocities decreased with increasing depth in the forest but relative turbulence intensities increased to mid-canopy levels. Leaf density influenced turbulence levels but not as strongly as did thermal stability. Thermal effects were adequately described by the single parameter h/L, where h is the canopy height and L is the Monin-Obukhov length. The longitudinal and vertical velocity correlation coefficient was larger in magnitude than expected in the upper layers of the forest but decreased to a small value in the lowest layers where the Reynolds stress was small. The ratio _w /u ^*, where u ^* is the local friction velocity, reflected changes in the uw correlation, becoming smaller than usual in the upper canopy layers. It is believed that these effects result from the intermittent, spatially coherent structures that are responsible for a large fraction of the momentum flux to the forest.     

Stability Dependence of Height Scales and Effective Roughness Lengths of Momentum and Heat Transfer Over Roughness Changes

It is widely accepted that the correct formulation of an effective roughness length, defined as the area average of the roughness length in heterogeneous terrain, relies upon the appropriate de-termination of a height scale. At this height a meteorological quantity is approximately in equilibrium with local surface conditions and independent of horizontal position. This research note determines explicitly the different height scales from the perturbation solutions of flow velocity and temperature, as well as the fluxes of momentum and heat, in a stratified boundary layer. These solutions are derived from the asymptotic approximation theory and shown to capture major characteristics of momentum and heat transfer over heterogeneous terrain with changes of the underlying roughness lengths. The effective roughness lengths can then be computed by use of these height scales. The dependence of height scales and effective roughness lengths upon stratification is also discussed briefly.     

Impact of Eddy Characteristics on Turbulent Heat and Momentum Fluxes in the Urban Roughness Sublayer

Eddy-covariance observations above the densely built-up Centre of Nanjing were made from December 2011 to August 2012. Separate eddy-covariance systems installed at two levels on a 36-m tower located on a rooftop were operated simultaneously, and observations grouped into two sectors (A, B) according to the prevalent wind directions. For sector A, where the nearby buildings are all below the lower measurement level, the sensible heat and momentum fluxes are generally greater at the upper level. For sector B, where several high-rise buildings are located upwind, the sensible heat and momentum fluxes at the upper level are close to those at the lower level. The analysis shows that the turbulent eddy characteristics differ between the two wind sectors, leading to a different behaviour of turbulent exchange between the two levels. A hypothesis is proposed that addresses the vertical variation of turbulent fluxes in the urban roughness sublayer (RSL). For sector A, the buildings block the flow, change the trajectory of scalars, and distort the footprint of scalar fluxes; this ‘blocking effect’ is believed to lead to a smaller sensible heat flux above the canopy layer. Such an effect should decrease with height in the RSL, explaining the increase of the observed turbulent heat flux with height. In addition, the presence of non-uniform building heights adversely affects turbulence organization around the canopy top, and likely elevates the inflection point of the mean flow to a higher elevation close to the upper measurement level, where larger shear results in a larger momentum flux. For sector B, wake effects from the nearby high-rise buildings strongly reduce turbulence organization at higher elevations, leading to similar sensible heat and momentum fluxes at both measurement levels.     

Footprint prediction of scalar fluxes from analytical solutions of the diffusion equation

The use of analytical solutions of the diffusion equation for footprint prediction is explored. Quantitative information about the footprint, i.e., the upwind area most likely to affect a downwind flux measurement at a given height z, is essential when flux measurements from different platforms, particularly airborne ones, are compared. Analytical predictions are evaluated against numerical Lagrangian trajectory simulations which are detailed in a companion paper (Leclerc and Thurtell, 1990). For neutral stability, the structurally simple solutions proposed by Gash (1986) are shown to be capable of satisfactory approximation to numerical simulations over a wide range of heights, zero displacements and roughness lengths. Until more sophisticated practical solutions become available, it is suggested that apparent limitations in the validity of some assumptions underlying the Gash solutions for the case of very large surface roughness (forests) and tentative application of the solutions to cases of small thermal instability be dealt with by semi-empirical adjustment of the ratio of horizontal wind to friction velocity. An upper limit of validity of these solutions for z has yet to be established.     

Latent and sensible heat flux predictions from a uniform pine forest using surface renewal and flux variance methods

A surface renewal model that links organized eddy motion to the latent and sensible heat fluxes is tested with eddy correlation measurements carried out in a 13m tall uniform Loblolly pine plantation in Duke Forest, Durham, North Carolina. The surface renewal model is based on the occurance of ramp-like patterns in the scalar concentration measurements. To extract such ramp-like patterns from Eulerian scalar concentration measurements, a newly proposed time-frequency filtering scheme is developed and tested. The time-domain filtering is carried out using compactly-supported orthonormal wavelets in conjunction with the Universal Wavelet Thresholding approach of Donoho and Johnstone, while the frequency filtering is carried out by a band-pass sine filter centered around the ramp-occurrence frequency as proposed by other studies. The method was separately tested for heat and water vapour with good agreement between eddy correlation flux measurements and model predictions. The usefulness of the flux-variance method to predict sensible and latent heat fluxes is also considered. Our measurements suggest that the simple flux-variance method reproduces the measured heat and momentum fluxes despite the fact that the variances were measured within the roughness sublayer and not in the surface layer. Central to the predictions of water vapour fluxes using the flux-variance approach is the similarity between heat and water vapour transport by the turbulent air flow. This assumption is also investigated for this uniform forest terrain.     

Landscape Roughness Parameters for Sherwood Forest – Experimental Results

The aim of this work is to present experimentally evaluated effective roughnesses (zoe) of a partly forested landscape. Although the ratio of boundary-layer height to obstacle size was only of the order of 50, there still seemed to exist a height range of 75–200 m where surface-layer similarity was approximately valid. Attempts were made to use conventional wind profile analysis to evaluate zoe, but the small height range and the large number of variables initially led to unacceptable uncertainties. Fixing the displacement height zd, rather than fitting it, reduced the data scatter to an acceptable level. The profile-derived roughness lengths zop obtained in this way were in good agreement with previous work, and with an alternative roughness length estimate zof for which flux-derived profile parameters u* and * were used. This implies that the profile-derived roughnesses were consistent with the measured surface-layer momentum flux. Comparison of both roughness estimates also yielded an improved estimate of the displacement height. Besides this, the authors tested a landscape roughness evaluation method which makes use of the gustiness parameter Tu = u/U in the surface layer. The results obtained by this method were in fair agreement with the profile-derived data. In previous work, the gustiness method was advocated because it could be used at relatively low levels, perhaps even within the roughness sub-layer. At the present measuring site, this was not the case as the gustiness method was only valid in an approximate way, and for a limited height range.     

Flux-gradient profiles for momentum and heat over an urban surface

Summary In this paper, we evaluate the applicability of flux-gradient relationships for momentum and heat for urban boundary layers within the Monin-Obukhov similarity (MOS) theory framework. Although the theory is widely used for smooth wall boundary layers, it is not known how well the theory works for urban layers. To address this problem, we measured the vertical profiles of wind velocity, air temperature, and fluxes of heat and momentum over a residential area and compared the results to theory. The measurements were done above an urban canopy whose mean height z_h is 7.3 m. 3-D sonic anemometers and fine wire thermocouples were installed at 4 heights in the region 1.5z_h < z < 4z_h. We found the following: (1) The non-dimensional horizontal wind speed has good agreement with the stratified logarithmic profile predicted using the semi-empirical Monin-Obukov similarity (MOS) function, when it was scaled by the surface friction velocity that is derived from the shear stress extrapolated to the roof-top level. (2) The scaled gradient of horizontal wind speed followed a conventional semi-empirical function for a flat surface at a level (z/z_h = 2.9), whereas, in the vicinity of the canopy height was larger than the commonly-used empirical relationship. (3) The potential temperature profile above the canopy shows dependency on the atmospheric stability and the scaled gradient of temperature is in good agreement with a conventional shear function for heat. In the case of heat, the dependency on height was not found. (4) The flux-gradient relationship for momentum and heat in the region 1.5z_h < z < 4z_h was rather similar to that for flat surfaces than that for vegetated canopies.     

Computational parameter estimation for a maize crop

During a whole growing season, the evolution of the displacement height, d, and roughness length, z ₀, of a maize crop has been estimated by a measurement programme. The results have been used to check different types of existing models to calculate these parameters from canopy characteristics only; a simple geometric model and two matching models have been investigated. A geometric model is based on geometric features of the surface only. After a simple modification, the geometric model gives good results for the displacement height as well as for the roughness length.A matching model, based on gradient-diffusion theory, yields good results for the displacement height. The roughness parameter, however, is overestimated by 17%. By a simple modification, the model results could be improved considerably.A matching model, based on a second-order closure procedure, yields excellent results for the displacement height and good results for the roughness length. But it appears that, when applying this model, the plant density index and plant area density distribution as a function of height must be well known.     

1998年南海季风爆发期间近海面层大气湍流结构和通量输送的观测研究

1998年5月14日至6月22日,在西沙永兴岛近海铁塔上进行了一次海-气通量观测试验,观测期包括了西南季风爆发前、爆发、爆发后风速加强等几个阶段。这次试验获得该年西南季风爆发的天气特征以及由涡度相关法、廓线法计算的动量、感热、潜热通量及湍流强度等一些统计量分布。分析结果表明,观测期湍流强度σu、σv、σw与平均风速之比为0.096、0.066、0.045;在近中性条件下（z/L≈0）,各相似函数基本为常数,σu/u＊≈3,σv/u＊≈2,σw/u＊≈1.25;在稳定条件下u、v、w三方向σ/u＊近似相等;在稳定与不稳定条件下σt/t＊随稳定度参数z/L的变化趋势相似,但符号相反。摩擦速度u＊随风速的变化接近于关系式u＊=0.029U10+0.006（U10为10m高度30min平均风速）,空气动力粗糙度长度z0变化在0.01～0.35mm之间。观测期动量通量变化在0.05～0.30N/m2之间,季风爆发期明显增大;曳力系数CD平均为1.12×10-3,它随平均风速变化可以表示为关系式103CD=0.003U210+0.020U10+0.836或103CD=0.056U10+0.732。感热通量由海洋输向大气,平均值为7.8W/m2,数值在0～15W/m2之间变化,季风爆发后明显增大;潜热通量数值一般变化于50～200W/m2,白天中午数值较高,夜间较低。季风爆发后也明显增大。观测期间鲍恩比（BowenRatio）日平均为0.05左右。