Cross‐shore variations of near‐surface wind velocity and atmospheric turbulence at the land‐sea boundary during CASP

Abstract

Airborne measurements of mean wind velocity and turbulence in the atmospheric boundary layer under wintertime conditions of cold offshore advection suggest that at a height of 50 m the mean wind speed increases with offshore distance by roughly 20% over a horizontal scale of order 10 km. Similarly, the vertical gust velocity and turbulent kinetic energy decay on scales of order 3.5 km by factors of 1.5 and 3.2, respectively. The scale of cross‐shore variations in the vertical fluxes of heat and downwind momentum is also 10 km, and the momentum flux is found to be roughly constant to 300 m, whereas the heat flux decreases with height. The stability parameter, z/L (where z = 50 m and L is the local Monin‐Obukhov length), is generally small over land but may reach order one over the warm ocean. The magnitude and horizontal length scales associated with the offshore variations in wind speed and turbulence are reasonably consistent with model results for a simple roughness change, but a more sophisticated model is required to interpret the combined effects of surface roughness and heat flux contrasts between land and sea.

Comparisons between aircraft and profile‐adjusted surface measurements of wind speed indicate that Doppler biases of 1–2 m s^?1 in the aircraft data caused by surface motions must be accounted for. In addition, the wind direction measurements of the Minimet anemometer buoy deployed in CASP are found to be in error by 25 ± 5°, possibly due to a misalignment of the anemometer vane. The vertical fluxes of heat and momentum show reasonably good agreement with surface estimates based on the Minimet data.     

Correction Of Eddy-Covariance Measurements Incorporating Both Advective Effects And Density Fluxes

Equations are presented to correct eddy-covariancemeasurements for both fluctuations in density andnon-zero mean advection, induced by convergence ordivergence of flow, and spatial source/sinkinhomogeneity, under steady-state and transientconditions. This correction collapses to theWebb–Pearman–Leuning expression ifthe mean vertical velocity is zero, and formally addsthe Webb–Pearman–Leuning expression to the correctionssuggested by Lee for conditions ofnon-zero vertical velocity and source/sink and meanscalar horizontal homogeneity. The equation requiresmeasurement of the mean vertical gradients of thescalar concentration of interest (air temperature,humidity, CO₂) as well as an accurateestimation of the mean vertical velocity, in additionto the vertical eddy covariance of the scalar. Simplemethods for the approximation of sensor tilt andcomplex terrain flow angle are presented, to allowestimation of non-zero mean vertical velocities. Theequations are applied to data from a maize crop and aforest to give examples of when the correction issignificant. In addition, a term for thethermodynamic expansion energy associated with watervapour flux is derived, which implies that the sonictemperature derived sensible heat flux will accuratelyinclude this contribution.     

稻麦轮作农田近地层湍流通量计算方案对比研究

稻麦轮作农田是我国典型农田类型,其模拟效果对我国农田气候模拟具有重要参考价值。气候中尺度模拟结果对近地层通量极为敏感,选择合适的通量计算方案对模拟效果至关重要。因此,对比分析稻麦轮作农田下不同的通量计算方案具有重要意义。本文选取了8种具有代表性的近地层湍流通量计算方案,采用寿县国家气候观象台实测资料对比分析了各方案在稻麦轮作农田的计算特征和差异。结果表明,在不同稳定度和不同风速情况下,各方案的误差特征各异。本文基于归一化标准差综合评价了各方案的准确度,总体而言,所有方案的动量通量总体平均归一化平均差为0.536,其中SS14(Sharan和Srivastava,2014年)方案最大为0.575,SS20(Sharan和Srivastava等,2020年)方案最小为0.517;所有方案的感热通量总体平均归一化标准差为0.638,其中GLGS20(Gryanik等,2020年)方案最大为0.871,SS14方案最小为0.476。此外,本研究还给出了稻麦轮作农田不同稳定层结和不同风速情况下各通量计算方案的误差特征。本文的研究结论,可为准确计算近地层湍流通量提供支撑。     

Estimating the components of the sensible heat budget of a tall forest canopy in complex terrain

Ultrasonic wind measurements, sonic temperature and air temperature data at two heights in the advection experiment MORE II were used to establish a complete budget of sensible heat including vertical advection, horizontal advection and horizontal turbulent flux divergence. MORE II took place at the long-term Carbo-Europe IP site in Tharandt, Germany. During the growing period of 2003 three additional towers were established to measure all relevant parameters for an estimation of advective fluxes, primarily of CO₂. Additionally, in relation to other advection experiments, a calculation of the horizontal turbulent flux divergence is proposed and the relation of this flux to atmospheric stability and friction velocity is discussed. In order to obtain a complete budget, different scaling heights for horizontal advection and horizontal turbulent flux divergence are tested. It is shown that neglecting advective fluxes may lead to incorrect results. If advective fluxes are taken into account, the sensible heat budget based upon vertical turbulent flux and storage change only, is reduced by approximately 30%. Additional consideration of horizontal turbulent flux divergence would in turn add 5–10% to this sum (i.e., the sum of vertical turbulent flux plus storage change plus horizontal and vertical advection). In comparison with available energy horizontal advection is important at night whilst horizontal turbulent flux divergence is rather insignificant. Obviously, advective fluxes typically improve poor nighttime energy budget closure and might change ecosystem respiration fluxes considerably.     

Parametric relations for the atmospheric boundary layer

Some parametric relations for the atmospheric planetary boundary layer (PBL) are suggested for possible use in the various atmospheric circulation and air quality models, as well as in other applications. These are for parameterizing the mean wind and temperature profiles, the vertical fluxes of momentum, heat and moisture, the variances of velocity fluctuations and length and time scales in the PBL. The parametric relations for the PBL height, the vertical velocity at the top of the PBL and the total energy dissipation in the PBL are also discussed. Experimental and/or theoretical bases for the various parametric relation are given. Some of the suggested parameterizations should be considered as tentative, until they are properly validated.     

黄土高原复杂地形受中尺度运动影响的稳定边界层湍流特征

利用兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University,简称SACOL)2008年12月观测资料,研究了稳定边界层湍流特征.使用涡动相关资料研究湍流通量时,定义湍流的平均时间τ内的中尺度运动是造成湍流统计量变化范围大的主要原因,稳定情形? τ取几十秒至几分钟.对梯度理查森数大于0.3的强稳定情形的湍流尺度分解(MRD)谱分析表明,感热通量在112.4～449.9 s存在谱隙,尺度大于谱隙的中尺度运动造成了通量观测资料离散性大,甚至有支配性影响.动量通量的谱隙在112.4～224.9 s之间.弱风时,中尺度运动的影响更大,垂直风速标准差以0.1的比率随中尺度风速变化;垂直风速标准差同广义风速表现出很好的相关性,并随着广义风速消失而消失.三维风速标准差与摩擦速度呈很好的线性关系,垂直、水平、横风风速的无量纲标准差分别为1.35、2.54、2.21.对湍流动能的研究发现,在梯度理查森数大于0.3的条件下,仍然存在连续的湍流.以湍动能为依据,分析了湍流的平稳时间长度,其长度随稳定度变化而变化,2008年12月7～11日从133.5 s变化到856.2 s,湍流平稳时间长度反映了中尺度运动的发生频率.     

雷暴大风过程中对流层中低层动量通量和动能通量输送特征研究

2014年5月31日北京发生一次雷暴大风过程。以雷达资料同化结果为初始场,对此次过程进行高分辨率数值模拟。采用非静力平衡和非地转平衡的经向动量方程和质量权重动能方程,利用模拟资料,对雷暴大风过程中经向动量和质量权重动能进行收支分析,以此来研究雷暴过程中对流层中低层动量通量和动能通量输送特征,讨论地面大风的可能成因。分析结果表明,在对流层中低层,经向动量通量散度是影响经向动量局地变化的主要强迫项。雷暴系统后部的入流把中低层的经向动量倾斜向下输送,系统前部对流云区中低层的下沉气流也向下输送经向动量。这两支下传动量通量先后与近地面经向动量的水平通量汇合,向系统前沿输送经向动量。在北京西北部地形阻挡作用下,经向动量通量在系统前端近地面辐合,促进那里的经向动量局地增长,有利于增强那里的南风。质量动能收支的特征与经向动量收支类似,在近地面层质量动能的局地变化主要是由质量动能通量散度引起的。系统后部入流把中层质量动能向下传输到近地面层,然后与近地面质量动能的水平通量汇合,向系统前沿输送质量动能。相对来说,近地面层经向动量和质量动能的水平通量比下传通量更重要,这主要与低层较强的东南急流有关。     

Energy- and Flux-Budget Turbulence Closure Model for Stably Stratified Flows. Part II: The Role of Internal Gravity Waves

We advance our prior energy- and flux-budget (EFB) turbulence closure model for stably stratified atmospheric flow and extend it to account for an additional vertical flux of momentum and additional productions of turbulent kinetic energy (TKE), turbulent potential energy (TPE) and turbulent flux of potential temperature due to large-scale internal gravity waves (IGW). For the stationary, homogeneous regime, the first version of the EFB model disregarding large-scale IGW yielded universal dependencies of the flux Richardson number, turbulent Prandtl number, energy ratios, and normalised vertical fluxes of momentum and heat on the gradient Richardson number, Ri. Due to the large-scale IGW, these dependencies lose their universality. The maximal value of the flux Richardson number (universal constant ≈0.2–0.25 in the no-IGW regime) becomes strongly variable. In the vertically homogeneous stratification, it increases with increasing wave energy and can even exceed 1. For heterogeneous stratification, when internal gravity waves propagate towards stronger stratification, the maximal flux Richardson number decreases with increasing wave energy, reaches zero and then becomes negative. In other words, the vertical flux of potential temperature becomes counter-gradient. Internal gravity waves also reduce the anisotropy of turbulence: in contrast to the mean wind shear, which generates only horizontal TKE, internal gravity waves generate both horizontal and vertical TKE. Internal gravity waves also increase the share of TPE in the turbulent total energy (TTE = TKE + TPE). A well-known effect of internal gravity waves is their direct contribution to the vertical transport of momentum. Depending on the direction (downward or upward), internal gravity waves either strengthen or weaken the total vertical flux of momentum. Predictions from the proposed model are consistent with available data from atmospheric and laboratory experiments, direct numerical simulations and large-eddy simulations.     

Measurement of the Sensible Eddy Heat Flux Based on Spatial Averaging of Continuous Ground-Based Observations

Using the standard eddy-covariance (EC) method to quantify mass and energy exchange at a single location usually results in an underestimation of vertical eddy fluxes at the surface. In order to better understand the reasons for this underestimation, an experimental set-up is presented that is based on spatial averaging of air temperature data from a network of ground-based sensors over agricultural land. For eight days during the 34-day observational period in May and June 2007, additional contributions to the sensible heat flux of more than 50Wm⁻² were measured in the lower surface layer by applying the spatial EC method as opposed to the standard temporal EC method. Smaller but still significant additional sensible heat fluxes were detected for four more days. The additional energy is probably transported in organised convective structures resulting in a mean vertical wind velocity unequal to zero at the tower location. The results show that convective transport contributes significantly to the surface energy budget for measurement heights as low as 2–3 m. Since these structures may be quasi-stationary, they can hardly be captured by a single-location measurement. The spatial EC set-up presented here is capable of quantifying contributions to the sensible heat flux from structures up to the scale of our spatial sensor network, which covered an area 3.5 × 3.5 km. For future experiments aiming at closing the energy balance, the spatial EC method should be employed to measure both the sensible and latent heat fluxes. Experimental determination of the horizontal advection of sensible and latent heat should also be considered, since such transport must occur due to convergence and divergence related to convection.     

2010年春季北京地区强沙尘暴过程的微气象学特征

利用北京大学校园地区PM10质量浓度观测资料、中国科学院大气物理研究所325m气象塔气象要素梯度和湍流观测资料,分析了北京地区2010年3月20～22日两次强沙尘暴过程微气象学要素和沙尘参量的时空演变以及湍流输送特征,为理解北京地区强沙尘暴天气沙尘输送规律和微气象学特征提供参考。结果表明:3月20～22日强沙尘暴过程前后不同高度温度先升后降,气压和相对湿度则相反。强沙尘暴来临时,高层风速先迅速增大,低层风速增加略有滞后,风切变明显加强,PM10浓度最大值和风速极大值出现时间较吻合。强沙尘暴过境时,不同高度向下的湍流动量输送、向上的湍流热量输送和湍流动能明显加强。与3月21日非沙尘暴日相比,强沙尘暴过程湍流动量通量增加,有利于沙尘粒子的水平和垂直输送过程;由于冷锋过境,水平热通量增大;垂直热通量因白天温度垂直梯度减小而减小,夜间因逆温层被破坏而增加;水平湍流动能对湍流动能占主要贡献,垂直湍流动能仅占水平湍流动能的10%～25%。     

A sensitive fast-response probe to measure turbulence and heat flux from any airplane

The theory, configuration, and accuracy of an inexpensive probe to measure turbulence from a small airplane are presented. The probe employs a nine-hole pressure-sphere design along with inprobe high-frequency pressure, temperature, and acceleration sensors. This sensor suite is specifically designed to extend mass, momentum and energy eddy-flux measurement to the higher frequencies characteristic of marine and nocturnal boundary layers. The probe is part of a mobile flux system, independent of the conveyance, which does not require a separate Inertial Navigation System.The new nine-port pressure sphere turbulence probe allows accurate turbulent velocity measurement with proper probe installation and appropriate computation technique for dynamic pressure. A thermistor in the central pressure port provides simultaneous temperature measurement, at a location symmetrical with respect to the flow, for accurate determination of true airspeed and heat flux. The probemounted temperature sensor gives heat fluxes with variance 5% of the mean in a weakly-turbulent marine boundary layer.Oak Ridge Associated Universities, assigned to NOAA/ATDD.     

Comparison of turbulence statistics within three boreal forest canopies

Three-dimensional sonic anemometers were used to measure velocities and temperatures within three natural boreal forest canopies. Vertical profiles of atmospheric turbulence statistics for a black spruce forest, a jack pine forest, and a trembling aspen forest, all located in southeastern Manitoba, were plotted and compared. The canopy structures were quite different, with total leaf-area indices of 2, 4 and 10, for the pine, aspen, and spruce forests, respectively.The profiles of the first and second moments differed among the canopies, where velocities decreased more rapidly in the top portions of the denser canopies. The velocity distributions were skewed and kurtotic within all canopies, and showed some differences among the canopies. Eulerian time scale profiles were generally similar among the canopies, and the vertical and streamwise time scale profiles were almost mirror images of each other. Eulerian length scale profiles showed some differences among canopies caused by differences in the velocity profiles. Ratios of vertical-to-horizontal time and length scales had a maximum in mid-canopy.Shear stress profiles were similar in the top parts of all canopies, and upward momentum fluxes were occasionally observed within the canopy trunk space. Countergradient heat fluxes were also observed sometimes. The countergradient fluxes and the skewed, kurtotic velocity distributions indicate the contribution of intermittent, large-scale eddies that are important for energy and mass transfer within canopies.     

洋面动量、感热和潜热通量计算的研究

整体传输公式是一种常用的计算洋面海气通量的方法。在作气候平均计算时,存在两种平均方法:方法一F_1=AB及方法二F_2=AB。显然,方法一是合理的,而方法二是对方法一的近似。本文利用COADS资料针对不同的平均时间计算分析了方法二对方法一的近似程度,结果表明:利用方法二计算长期气候平均,平滑时间超过5天则对计算结果必须进行订正;对感热和潜热通量,月平均方法二在20％的精度内可代替方法一,而旬平均方法二则可在10％的精度内代替方法一;动量通量月平均和旬平均精度分别为40％和20％。因此方法二在计算月平均动量通量时必须作订正,或直接用方法一进行计算,本文同时计算比较了1972年(E1 Nino年)和1975年(非E1 Nino年)太平洋洋面的各项能量及总能量收支,发现1972年洋面能量收支年变化幅度比1975年小得多,其差别主要决定于感热和潜热的变化。     

Turbulence Measurements above a Pine Forest

Eddy fluxes of momentum, sensible and latent heat, and turbulence spectra measured over the Thetford Forest during 10 days in the Spring of 1973 are described. The measured total heat flux (H + λE) for 122 20-min periods agreed closely on average with independent estimates from an energy balance method. There was evidence that the energy balance data gave small systematic overestimates of available energy during the hours before noon, compensated by slight underestimates for the remainder of the day. A comparison of measured wind speeds and friction velocities in neutral stability confirmed the validity of the aerodynamic method for estimating momentum fluxes at heights of a few roughness lengths above the canopy. In stable conditions the log-linear wind profile $U=(u_{?t}/k)({? ln}((z-d)/z_{0})+←pha(z-d-z_{0})/L)$ with α = 3.4 ± 0.4 provided a good fit to the data. Spectra in unstable conditions were generally more sharply peaked than those measured by other workers over smoother terrain: differences were less marked in the case of vertical velocity in stable conditions. Temperature spectra in these stable conditions showed high energy at relatively low wavenumbers, and wT cospectra showed a cospectral gap; both of these results were associated with an intermittent sawtooth structure in the temperature fluctuations.     

Towards Closing the Surface Energy Budget of a Mid-latitude Grassland

Observations for May and August, 2005, from a long-term grassland meteorological station situated in central Netherlands were used to evaluate the closure of the surface energy budget. We compute all possible enthalpy changes, such as the grass cover heat storage, dew water heat storage, air mass heat storage and the photosynthesis energy flux, over an averaging time interval. In addition, the soil heat flux was estimated using a harmonic analysis technique to obtain a more accurate assessment of the surface soil heat flux. By doing so, a closure of 96% was obtained. The harmonic analysis technique appears to improve closure by 9%, the photosynthesis for 3% and the rest of the storage terms for a 3% improvement of the energy budget closure. For calm nights (friction velocity u _* < 0.1 m s⁻¹) when the eddy covariance technique is unreliable for measurement of the vertical turbulent fluxes, the inclusion of a scheme that calculates dew fluxes improves the energy budget closure significantly.     

鄱阳湖地区复杂地表条件下一次强降水过程的近地面边界层特征

使用鄱阳湖北部70 m气象塔湍流和梯度观测数据,分析了2011年6月6日夜间一次暴雨过程中近地面边界层特征.结果表明,此次过程是在高空低槽和西南急流的天气背景下,受鄱阳湖复杂地表影响产生的局地性强降水.强降水发生前受东南暖平流影响,近地面边界层中水汽累积,不稳定性增加;强降水过程中,近地层感热、潜热通量迅速增加,同时,近地面层湍流动量通量下传和水平输送增加,鄱阳湖的水汽输送加强降水强度.另外,强降水过程中,近地面湍流动能迅速增大并达到最大值,而平均动能的增大发生在强降水结束后,表明地表作用明显,近地面边界层的湍流场为暴雨提供动力条件.尺度分析表明,强降水前,中尺度动量通量占主要地位,降水过程中湍流动量通量显著加强.     

Air-sea interaction including a shallow and coastal zone

We performed special experimental and theoretical research concerning the energy--mass exchange process in shallow waters, which took into account the influence of the basin depth on the evaporation and sensible heat exchange. Data was obtained from deep water basins, as well as from shallow waters, at an open sea, as well as in its coastal zone. A new parameterization model of evaporation and friction velocity from shallow water surfaces under different wind velocities was investigated. Results of models and measurements of the energy/mass exchange intensity of a small shallow lake (LITFASS-experiments) were compared. The validation of these models with the eddy-covariance measurements of the LITFASS-98 and LITFASS-2003 experiments showed good results for the wind sector, with good fetch conditions. Therefore, the models examined may be used for calculating the evaporation of lakes where a standard data set of wind velocity, air and water temperature, air moisture and the depth of the lake is available. It should be remembered that in a coastal zone the roughness of the water surface is transformed and therefore there are additional calculation difficulties for the intensity of the energy-mass exchange. Based on our experimental data of the basin depth influence on the water-atmosphere exchange, a new model for the calculation of the energy-mass exchange in a coastal zone was developed. Our new model in combination with the empirical dependence for the calculation of the energy-mass exchange in a coastal zone allows calculation of the momentum, heat and humidity fluxes values at different distances from a shore.     

Modelling Vegetation-Atmosphere Co2 Exchange By A Coupled Eulerian-Langrangian Approach

A Eulerian-Lagrangian canopy microclimate model wasdeveloped with the aim of discerning physical frombiophysical controls of CO₂ and H₂O fluxes. The model couples radiation attenuation with mass,energy, and momentum exchange at different canopylevels. A unique feature of the model is its abilityto combine higher order Eulerian closure approachesthat compute velocity statistics with Lagrangianscalar dispersion approaches within the canopy volume. Explicit accounting for within-canopy CO₂,H₂O, and heat storage is resolved by consideringnon-steadiness in mean scalar concentration andtemperature. A seven-day experiment was conducted inAugust 1998 to investigate whether the proposedmodel can reproduce temporal evolution of scalar(CO₂, H₂O and heat) fluxes, sources andsinks, and concentration profiles within and above auniform 15-year old pine forest. The modelreproduced well the measured depth-averaged canopy surfacetemperature, CO₂ and H₂O concentrationprofiles within the canopy volume, CO₂ storageflux, net radiation above the canopy, and heat andmass fluxes above the canopy, as well as the velocitystatistics near the canopy-atmosphere interface. Implications for scaling measured leaf-levelbiophysical functions to ecosystem scale are alsodiscussed.     

Increased heat fluxes near a forest edge

Summary ?Observations of sensible and latent heat flux above forest downwind of a forest edge show these fluxes to be larger than the available energy over the forest. The enhancement averages to 56 W m⁻², or 16% of the net radiation, at fetches less than 400 m, equivalent to fetch to height ratios less than 15. The enhancement of turbulent energy fluxes is explained by advection and increases with the difference in temperature and humidity of the air over the upwind area as compared to the forest. The relatively high temperature and humidity of the upwind air are not caused by high surface heat fluxes, but are explained by the relatively low aerodynamic roughness of the upwind surface. Although the heat fluxes over forest are enhanced, the momentum fluxes are almost adjusted to the underlying forest. The different behaviour of heat and momentum fluxes is explained by absorption of momentum by pressure gradients near the forest edge. It is concluded that fetch requirements to obtain accurate surface fluxes from atmospheric observations need to be more stringent for scalar fluxes as compared to momentum fluxes. Received November 23, 2001; accepted May 13, 2002     

A bulk model for the atmospheric planetary boundary layer

The integrated momentum and thermodynamic equations through the planetary boundary layer (PBL) are solved numerically to predict the mean changes of wind and potential temperature from which surface fluxes are computed using bulk transfer coefficients of momentum and heat. The second part of the study involves a formulation and testing of a PBL height model based on the turbulent energy budget equation where turbulent fluxes of wind and heat are considered as the source of energy. The model exhibits capability of predicting the PBL height development for both stable and unstable regimes of observed conditions. Results of the model agree favourably with those of Deardorff's (1974a) and Tennekes' (1973) models in convective conditions.Contribution number 396.