Flux contribution of coherent structures and its implications for the exchange of energy and matter in a tall spruce canopy

The flux contribution of coherent structures to the total exchange of energy and matter is investigated in a spruce canopy of moderate density in heterogeneous, complex terrain. The study deploys two methods of analysis to estimate the coherent exchange: conditional averages in combination with wavelet analysis, and quadrant analysis. The data were obtained by high-frequency single-point measurements using sonic anemometers and gas analysers at five observation heights above and within the canopy and subcanopy, and represent a period of up to 2.5 months. The study mainly addresses the momentum transfer and exchange of sensible heat throughout the roughness sublayer, while results are provided for the exchange of carbon dioxide and water vapour above the canopy. The magnitude of the flux contribution of coherent structures largely depends on the method of analysis, and it is demonstrated that these differences are attributed to differences in the sampling strategy between the two methods. Despite the differences, relational properties such as sweep and ejection ratios and the variation of the flux contribution with height were in agreement for both methods. The sweep phase of coherent structures is the dominant process close to and within the canopy, whereas the ejections gain importance with increasing distance to the canopy. The efficiency of the coherent exchange in transporting scalars exceeds that for momentum by a factor of two. The occurrence of coherent structures results in a flux error less than 4% for the eddy-covariance method. Based on the physical processes identified from the analysis of the ejection and sweep phases along the vertical profile in the roughness sublayer, a classification scheme for the identification of exchange regimes is developed. This scheme allows one to estimate the region of the canopy participating in the exchange of energy and matter with the above-canopy air under varying environmental conditions.     

Concurrency of Coherent Structures and Conditionally Sampled Daytime Sub-canopy Respiration

We investigated an alternative means for quantifying daytime ecosystem respiration from eddy-covariance data in three forests with different canopy architecture. Our hypothesis was that the turbulent transport by coherent structures is the main pathway for carrying detectable sub-canopy respiration signals through the canopy. The study extends previously published work by incorporating state-of-the-art wavelet decomposition techniques for the detection of coherent structures. Further, we investigated spatial and temporal variability of the respiration signal and coherent exchange at multiple heights, for three mature forest sites with varying canopy and terrain properties for one summer month. A connection between the coherent structures and identified sub-canopy respiration signal was clearly determined. Although not always visible in signals collected above the canopy, certain cases showed a clear link between conditionally sampled respiration events and coherent structures. The dominant time scales of the coherent structure ejection phase (20?C30 s), relative timing of maximum coincidence between respiration events and the coherent structure ejection phase (at approximately ?10 s from detection) and vertical transport upward through the canopy were shown to be consistent in time, across measurement heights and across the different forest sites. Best results were observed for an open canopy pine site. We conclude that the presented method is likely to be applicable at more open rather than dense (closed) canopies. The results provided a confirmation of the connection between below- and above-canopy scalar time series, and may help the development or refinement of direct methods for the determination of component fluxes from observations above the canopy.     

Coherent Structures at a Forest Edge: Properties, Coupling and Impact of Secondary Circulations

Little is known about the influence of coherent structures on the exchange process, mainly in the case of forest edges. Thus, in the framework of the ExchanGE processes in mountainous Regions (EGER) project, measurements of atmospheric turbulence were taken at different heights between a forest and an adjacent clear cutting using sonic anemometers and high-frequency optical gas analyzers. From these turbulence data, dominant coherent structures were extracted using an already existing wavelet methodology, which was developed for homogeneous forest canopies. The aim of this study is to highlight differences in properties of coherent structures between a forest and a clear cutting. Distinct features of coherent exchange at the forest edge are presented and a careful investigation of vertical and horizontal coupling by coherent structures around the surface heterogeneity is made. Within the forest, coherent structures are less frequent but possess larger time scales, indicating that only the largest coherent motions can penetrate through the forest canopy. At the forest edge, there is no crown layer that can hinder the vertical exchange of coherent structures, because these exhibit similar time scales at all heights. In contradiction to that, no improved vertical coupling was detected at the forest edge. This is mainly because the structures captured by the applied routine contribute less to total turbulent fluxes at the edge than within the forest. Thus, coherent structures with time scales between 10 and 40 s are not the dominant exchange mechanism at the forest edge. With respect to the horizontal direction, a consistent picture of coherent transport could be derived: along the forest edge there is mainly good coupling by coherent structures, whereas perpendicular to the forest edge there is mainly decoupling. Finally, it was found that there is a systematic modulation of coherent structures directly at the forest edge: strong ejection motions appear in all time series during the daytime, whereas strong sweeps dominate at night. An effect of wind direction relative to the forest edge is excluded. Consequently, it is hypothesized that this might be an indication of a quasi-stationary secondary circulation above the clear cutting that develops due to differences in surface temperature and roughness. Such circulations might be a relevant turbulent transport mechanism for ecosystem-atmosphere exchange in heterogeneous landscapes.     

高风速相干结构对通量输送影响的实验研究

切变湍流的相干结构是湍流研究中的重大发现,它表明湍流运动并非完全随机,其中具有可检测的有序结构.本文通过处理南京浦口地区大气边界层观测数据,来分析不稳定层结中高风速相干结构特征.本次观测项目包括对场地中央的气象铁塔上2 m和40 m高度上超声风速仪的脉动速度、温度测量以及风廓线雷达对边界层风速廓线的测量.对超声水平风速时间序列数据进行小波变换 (时间尺度400 s),通过阈值来识别这种高风速相干结构.与多普勒风廓线雷达测量结果对比后发现,这种方法确定的相干结构符合常规的认识,具有较长的时间尺度和较大的垂直尺度 (接近边界层厚度).分析三天相干结构特性得到无量纲空间间隔约为6,即每隔6个边界层厚度的水平位置出现一个高速相干结构.通过与垂直风速小波系数的比较,发现高风速相干结构与向下垂直风速之间有较好相关,这与湍流中 “阵风” 现象的研究结论相似.使用四象限分析方法分类得到两种动量通量输送为负的运动:较小水平风速的上扬 (ejection) 运动 (简称为上扬运动) 和较大水平风速的下扫 (sweep) 运动 (简称为下扫运动),这两种运动在整个湍流活动中处于主导地位.高风速相干结构通过促进下扫运动和抑制上扬运动来影响动量通量的输送.     

峨眉山地区近地层微气象特征研究

基于2019年12月至2020年11月峨眉山站梯度塔资料、辐射观测资料和地表通量资料,采用涡动相关法对峨眉山地区近地层的地表通量和蒸散发量的变化进行分析,并估算了零平面位移、空气动力粗糙度、空气热力粗糙度、动量通量输送系数和感热通量输送系数等重要的空气动力学和热力学参数.研究表明:近地面风速呈现高层高、低层低的特征,且...     

On the flow structure and turbulence during sweep and ejection events in a wind-tunnel model canopy

Particle image velocimetry (PIV) data obtained in a wind-tunnel model of a canopy boundary layer is used to examine the characteristics of mean flow and turbulence. The vector spacing varies between 1.7 and 2.5 times the Kolmogorov scales. Conditional sampling based on quadrants, i.e. based on the signs of velocity fluctuations, reveals fundamental differences in flow structure, especially between sweep and ejection events, which dominate the flow. During sweeps, the downward flow generates a narrow, highly turbulent, shear layer containing multiple small-scale vortices just below canopy height. During ejections, the upward flow expands this shear layer and the associated small-scale flow structures to a broad region located above the canopy. Consequently, during sweeps the turbulent kinetic energy (TKE), Reynolds stresses, as well as production and dissipation rates, have distinct narrow peaks just below canopy height, whereas during ejections these variables have broad maxima well above the canopy. Three methods to estimate the dissipation rate are compared, including spectral fits, measured subgrid-scale (SGS) energy fluxes at different scales, and direct measurements of slightly underresolved instantaneous velocity gradients. The SGS energy flux is 40–60% of the gradient-based (direct) estimates for filter sizes inside the inertial range, while decreasing with scale, as expected, within the dissipation range. The spectral fits are within 5–30% of the direct estimates. The spectral fits exceed the direct estimates near canopy height, but are lower well above and below canopy height. The dissipation rate below canopy height increases with velocity magnitude, i.e. it has the highest values during sweep and quadrant 1 events, and is significantly lower during ejection and quadrant 3 events. Well above the canopy, ejections are the most dissipative. Turbulent transport during sweep events acts as a source below the narrow shear layer within the canopy and as a sink above it. Transport during ejection events is a source only well above the canopy. The residual term in the TKE transport equation, representing mostly the effect of pressure–velocity correlations, is substantial only within the canopy, and is dominated by sweeps.     

A wind tunnel study of air flow in waving wheat: Single-point velocity statistics

We analyse single-point velocity statistics obtained in a wind tunnel within and above a model of a waving wheat crop, consisting of nylon stalks 47 mm high and 0.25 mm wide in a square array with frontal area index 0.47. The variability of turbulence measurements in the wind tunnel is illustrated by using a set of 71 vertical traverses made in different locations, all in the horizontally-homogeneous (above-canopy) part of the boundary layer. Ensemble-averaged profiles of the statistical moments up to the fourth order and profiles of Eulerian length scales are presented and discussed. They are consistent with other similar experiments and reveal the existence of large-scale turbulent coherent structures in the flow. The drag coefficient in this canopy as well as in other reported experiments is shown to exhibit a characteristic height-dependency, for which we propose an interpretation. The velocity spectra are analysed in detail; within and just above the canopy, a scaling based on fixed length and velocity scales (canopy height and mean horizontal wind speed at canopy top) is proposed. Examination of the turbulent kinetic energy and shear stress budgets confirms the role of turbulent transport in the region around the canopy top, and indicates that pressure transport may be significant in both cases. The results obtained here show that near the top of the canopy, the turbulence properties are more reminiscent of a plane mixing layer than a wall boundary layer.     

Coherent Turbulent Structures Across a Vegetation Discontinuity

The study of turbulent flow across a vegetation discontinuity is of significant interest as such landscape features are common, and as there is no available theory to describe this regime adequately. We have simulated the three-dimensional dynamics of the airflow across a discontinuity between a forest (with a leaf area index of 4) and a clearing surface using large-eddy simulation. The properties of the bulk flow, as well as the large-scale coherent turbulent structures across the forest-to-clearing transition and the clearing-to-forest transition, are systematically explored. The vertical transport of the bulk flow upstream of the leading edge gives rise to the enhanced gust zone around the canopy top, while the transport downstream of the trailing edge leads to the formation of a recirculation zone above the clearing surface. The large-scale coherent structures across the two transitions exhibit both similarities with and differences from those upstream of the corresponding transition. For example, the ejection motion is dominant over the sweep motion in most of the region 1?<?z/h < 2 (h is the canopy height) immediately downstream of the trailing edge, much as in the forested area upstream. Also, the streamwise vortex pair, which has previously been observed within the canopy sublayer and the atmospheric boundary layer, is consistently found across both transitions. However, the inflection observed both in the mean streamwise velocity, as well as in the vertical profiles of the coherent structures in the forested area, disappears gradually across the forest-to-clearing transition. The coherence of the turbulence, quantified by the percentage of the total turbulence kinetic energy that the coherent structures capture from the flow, decreases sharply immediately downstream of the trailing edge of the forest and increases downstream of the leading edge of the forest. The effects of the ratio of the forest/clearing lengths under a given streamwise periodicity on flow statistics and coherent turbulent structures are presented as well.     

Observations of Coherent Turbulence Structures in the Near-Neutral Atmospheric Boundary Layer

Turbulence structures of high Reynolds number flow in the near-neutral atmospheric boundary layer (ABL) are investigated based on observations at Shionomisaki and Shigaraki, Japan. A Doppler sodar measured the vertical profiles of winds in the ABL. Using the integral wavelet transform for the time series of surface wind data, the pattern of a descending high-speed structure with large vertical extent (from the surface to more than 200-m level) is depicted from the Doppler sodar data. Essentially this structure is a specific type of coherent structure that has been previously shown in experiments on turbulent boundary-layer flows. Large-scale high-speed structures in the ABL are extracted using a long time scale (240 s) for the wavelet transform. The non-dimensional interval of time between structures is evaluated as 3.0–6.2 in most cases. These structures make a large contribution to downward momentum transfer in the surface layer. Quadrant analyses of the turbulent motion measured by the sonic anemometer (20-m height) suggest that the sweep motion (high-speed downward motion) plays a substantial role in the downward momentum transfer. In general, the contribution of sweep motions to the momentum flux is nearly equal to that of ejection motions (low-speed upward motions). This contribution of sweep motions is related to the large-scale high-speed structures.     

Divergence of turbulent fluxes in the surface layer: case of a coastal city

This study quantifies the processes that take place in the layer between the mean building height and the measurement level of an energy balance micrometeorological tower located in the dense old core of a coastal European city. The contributions of storage, vertical advection, horizontal advection and radiative divergence for heat are evaluated with the available measurements and with a three-dimensional, high-resolution meteorological simulation that had been evaluated against observations. The study focused on a summer period characterized by sea-breeze flows that affect the city. In this specific configuration, it appears that the horizontal advection is the dominant term. During the afternoon when the sea breeze is well established, correction of the sensible heat flux with horizontal heat advection increases the measured sensible heat flux up to 100 W m⁻². For latent heat flux, the horizontal moisture advection converted to equivalent latent heat flux suggests a decrease of 50 W m⁻². The simulation reproduces well the temporal evolution and magnitude of these terms.     

Turbulent Structures in a Pine Forest with a Deep and Sparse Trunk Space: Stand and Edge Regions

Forested landscapes often exhibit large spatial variability in vertical and horizontal foliage distributions. This variability may affect canopy-atmosphere exchanges through its action on the development of turbulent structures. Here we investigate in neutral stratification the turbulent structures encountered in a maritime pine forest characterized by a high, dense foliated layer associated with a deep and sparse trunk space. Both stand and edge regions are considered. In situ measurements and the results of large-eddy simulations are used and analyzed together. In stand conditions, far from the edge, canopy-top structures appear strongly damped by the dense crown layer. Turbulent wind fluctuations within the trunk space, where the momentum flux vanishes, are closely related to these canopy-top structures through pressure diffusion. Consequently, autocorrelation and spectral analyses are not quite appropriate to characterize the vertical scale of coherent structures in this type of canopy, as pressure diffusion enhances the actual scale of structures. At frequencies higher than those associated with canopy-top structures, wind fluctuations related to wake structures developing behind tree stems are observed within the trunk space. They manifest themselves in wind velocity spectra as secondary peaks in the inertial subrange region, confirming the hypothesis of spectral short-cuts in vegetation canopies. In the edge region specific turbulent structures develop just below the crown layer, in addition to canopy-top structures. They are generated by the wind shear induced by the sub-canopy wind jet that forms at the edge. These structures provide a momentum exchange mechanism similar to that observed at the canopy top but in the opposite direction and with a lower magnitude. They may develop as in plane mixing-layer flows, with some perturbations induced by canopy-top structures. Wake structures are also observed within the trunk space in the edge region.     

Variability of Sub-Canopy Flow,Temperature, and Horizontal Advection in Moderately Complex Terrain

We examine the space–time structure of the wind and temperature fields, as well as that of the resulting spatial temperature gradients and horizontal advection of sensible heat, in the sub-canopy of a forest with a dense overstorey in moderately complex terrain. Data were collected from a sensor network consisting of ten stations and subject to orthogonal decomposition using the multiresolution basis set and stochastic analyses including two-point correlations, dimensional structure functions, and various other bulk measures for space and time variability. Despite some similarities, fundamental differences were found in the space–time structure of the motions dominating the variability of the sub-canopy wind and temperature fields. The dominating motions occupy similar spatial, but different temporal, scales. A conceptual space–time diagram was constructed based on the stochastic analysis that includes the important end members of the spatial and temporal scales of the observed motions of both variables. Short-lived and small-scale motions govern the variability of the wind, while the diurnal temperature oscillation driven by the surface radiative transfer is the main determinant of the variability in the temperature signal, which occupies much larger time scales. This scale mismatch renders Taylor’s hypothesis for sub-canopy flow invalid and aggravates the computation of meaningful estimates of horizontal advective fluxes without dense spatial information. It may further explain the ambiguous and inconclusive results reported in numerous energy and mass balance and advection studies evaluating the hypothesis that accounting for budget components other than the change in storage term and the vertical turbulent flux improves the budget closure when turbulent diffusion is suppressed in plant canopies. Estimates of spatial temperature gradients and advective fluxes were sensitive to the network geometry and the spatial interpolation method. The assumption of linear spatial temperature gradients was not supported by the results, and leads to increased spatial and temporal variability of inferred spatial gradients and advection estimates. A method is proposed to estimate the appropriate minimum network size of wind and temperature sensors suitable for an evaluation of energy and mass balances by reducing spatial and temporal variability of the spatially sampled signals, which was estimated to be on the order of 200 m at the study site.     

Organised Motion in a Tall Spruce Canopy: Temporal Scales,Structure Spacing and Terrain Effects

This study investigates the organised motion near the canopy-atmosphere interface of a moderately dense spruce forest in heterogeneous, complex terrain. Wind direction is used to assess differences in topography and surface properties. Observations were obtained at several heights above and within the canopy using sonic anemometers and fast-response gas analysers over the course of several weeks. Analysed variables include the three-dimensional wind vector, the sonic temperature, and the concentration of carbon dioxide. Wavelet analysis was used to extract the organised motion from time series and to derive its temporal scales. Spectral Fourier analysis was deployed to compute power spectra and phase spectra. Profiles of temporal scales of ramp-like coherent structures in the vertical and longitudinal wind components showed a reversed variation with height and were of similar size within the canopy. Temporal scales of scalar fields were comparable to those of the longitudinal wind component suggesting that the lateral scalar transport dominates. The existence of a – 1 power law in the longitudinal power spectra was confirmed for a few cases only, with a majority showing a clear 5/3 decay. The variation of effective scales of organised motion in the longitudinal velocity and temperature were found to vary with atmospheric stability, suggesting that both Kelvin-Helmholtz instabilities and attached eddies dominate the flow with increasing convectional forcing. The canopy mixing-layer analogy was observed to be applicable for ramp-like coherent structures in the vertical wind component for selected wind directions only. Departures from the prediction of m = Λ _w L _s ⁻¹ = 8–10 (where Λ _w is the streamwise spacing of coherent structures in the vertical wind w and L _s is a canopy shear length scale) were caused by smaller shear length scales associated with large-scale changes in the terrain as well as the vertical structure of the canopy. The occurrence of linear gravity waves was related to a rise in local topography and can therefore be referred to as mountain-type gravity waves. Temporal scales of wave motion and ramp-like coherent structures were observed to be comparable.     

Nocturnal Intermittent Coupling Between the Interior of a Pine Forest and the Air Above It

A 1-year set of measurements of CO₂ and energy turbulent fluxes above and within a 25-m pine forest in southern Brazil is analyzed. The study focuses on the coupling state between two levels and its impact on flux determination by the eddy-covariance method. The turbulent series are split in their typical temporal scales using the multiresolution decomposition, a method that allows proper identification of the time scales of the turbulent events. Initially, four case studies are presented: a continually turbulent, a continually calm, a calm then turbulent, and an intermittent night. During transitions from calm to turbulent, large scalar fluxes of opposing signs occur at both levels, suggesting the transference of air accumulated in the canopy during the stagnant period both upwards and downwards. Average fluxes are shown for the entire period as a function of turbulence intensity and a canopy Richardson number, used as an indicator of the canopy coupling state. Above the canopy, CO₂ and sensible heat fluxes decrease in magnitude both at the neutral and at the very stable limit, while below the canopy they increase monotonically with the canopy Richardson number. Latent heat fluxes decrease at both levels as the canopy air becomes more stable. The average temporal scales of the turbulent fluxes at both levels approach each other in neutral conditions, indicating that the levels are coupled in that case. Average CO₂ fluxes during turbulent periods that succeed very calm ones are appreciably larger than the overall average above the canopy and smaller than the average or negative within the canopy, indicating that the transfer of air accumulated during calm portions at later turbulent intervals affects the flux average. The implications of this process for mean flux determination are discussed.     

Turbulent Flux Measurements Above and Below the Overstory of a Boreal Aspen Forest

Turbulent flux measurements both above and beneath the canopy of a boreal aspen forest are described. Velocity skewness showed that, beneath the aspen canopy, turbulence was dominated by intermittent, downward penetrating gusts. Eulerian horizontal length scales calculated from integration of the autocorrelation function or spectral peaks were 9.0 and 1.4 times the mean aspen height of 21.5 m respectively. Above-canopy power spectral slopes for all velocity components followed the -2/3 power law, whereas beneath-canopy slopes were closer to -1 and showed a spectral short cut in the horizontal and vertical components. Cospectral patterns were similar both above and beneath the canopy. The Monin–Obukhov similarity function for the vertical wind velocity variance was a well-defined function of atmospheric stability, both above and beneath the canopy. Nocturnal flux underestimation and departures of this similarity function from that expected from Monin–Obukhov theory were a function of friction velocity. Energy balance closure greater than 80% was achieved at friction velocities greater than 0.30 and 0.10 m s-1, above and below the aspen canopy, respectively. Recalculating the latent heat flux using various averaging periods revealed a minimum of 15 min were required to capture 90% of the 30-min flux. Linear detrending reduced the flux at shorter averaging periods compared to block averaging. Lack of energy balance closure and erratic flux behaviour led to the recalculation of the latent and sensible heat fluxes using the ratio of net radiation to the sum of the energy balance terms.     

Observations of planview flux patterns within convective structures of the marine atmospheric surface layer

Air/sea flux variability on horizontal scales from 50 m to several km results, in part, from the presence of coherent convective structures within the atmospheric boundary layer. The horizontal distribution of fluxes within these convective updrafts and downdrafts is, therefore, central to studies of air/sea interaction and remote sensing of sea surface wind and wave fields. This study derives these flux patterns from observations of the Marine Atmospheric Surface Layer (MASL).Research aircraft flights through the MASL provide an optimal means for sampling large numbers of the above-mentioned coherent structures. The NCAR Electra flew numerous legs through the MASL at a height of 50 m during the 1987 stratocumulus phase of Project FIRE (First ISSCP (International Satellite Cloud Climatology Program) Regional Experiment).In situ measurements from these legs serve as the dataset for this paper. The data are processed in such a way as to retain only the turbulence fluctuations. Conditional sampling, based on the vertical velocity field, results in the isolation of convective updrafts and downdrafts. Compositing of the data for these two classes of convective drafts results in horizontal planviews of the vertical fluxes of buoyancy, absolute humidity, along-meanwind component of momentum, and vertical velocity. To ensure dynamical similarity, these horizontal planviews are oriented in a coordinate system aligned with the mean wind.     

Surface Temperature and Surface-Layer Turbulence in a Convective Boundary Layer

Previous laboratory and atmospheric experiments have shown that turbulence influences the surface temperature in a convective boundary layer. The main objective of this study is to examine land-atmosphere coupled heat transport mechanism for different stability conditions. High frequency infrared imagery and sonic anemometer measurements were obtained during the boundary layer late afternoon and sunset turbulence (BLLAST) experimental campaign. Temporal turbulence data in the surface-layer are then analyzed jointly with spatial surface-temperature imagery. The surface-temperature structures (identified using surface-temperature fluctuations) are strongly linked to atmospheric turbulence as manifested in several findings. The surface-temperature coherent structures move at an advection speed similar to the upper surface-layer or mixed-layer wind speed, with a decreasing trend with increase in stability. Also, with increasing instability the streamwise surface-temperature structure size decreases and the structures become more circular. The sequencing of surface- and air-temperature patterns is further examined through conditional averaging. Surface heating causes the initiation of warm ejection events followed by cold sweep events that result in surface cooling. The ejection events occur about 25 % of the time, but account for 60–70 % of the total sensible heat flux and cause fluctuations of up to 30 % in the ground heat flux. Cross-correlation analysis between air and surface temperature confirms the validity of a scalar footprint model.     

Early-Morning Flow Transition in a Valley in Low-Mountain Terrain Under Clear-Sky Conditions

We investigate the evolution of the early-morning boundary layer in a low-mountain valley in south-western Germany during COPS (convective and orographically induced precipitation study) in summer 2007. The term low-mountain refers to a mountainous region with a relief of gentle slopes and with an absolute altitude that remains under a specified height (usually 1,500 m a.s.l.). A subset of 23 fair weather days from the campaign was selected to study the transition of the boundary-layer flow in the early morning. The typical valley atmosphere in the morning hours was characterized by a stable temperature stratification and a pronounced valley wind system. During the reversal period—called the low wind period—of the valley wind system (duration of 1–2 h), the horizontal flow was very weak and the conditions for free convection were fulfilled close to the ground. Ground-based sodar observations of the vertical wind show enhanced values of upward motion, and the corresponding statistical properties differ from those observed under windless convective conditions over flat terrain. Large-eddy simulations of the boundary-layer transition in the valley were conducted, and statistical properties of the simulated flow agree with the observed quantities. Spatially coherent turbulence structures are present in the temporal as well as in the ensemble mean analysis. Thus, the complex orography induces coherent convective structures at predictable, specific locations during the early-morning low wind situations. These coherent updrafts, found in both the sodar observations and the simulation, lead to a flux counter to the gradient of the stably stratified valley atmosphere and reach up to the heights of the surrounding ridges. Furthermore, the energy balance in the surface layer during the low wind periods is closed. However, it becomes unclosed after the onset of the valley wind. The partition into the sensible and the latent heat fluxes indicates that missing flux components of sensible heat are the main reason for the unclosed energy balance in the considered situations. This result supports previously published investigations on the energy balance closure.     

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.     

昆明城市热岛效应的数值模拟研究

本文利用WRF（V3.9.1）模式中耦合Noah/SLUCM方案作为Control试验,研究了土地利用类型（Md04试验）、陆面过程（NoUCM试验）和湖泊（Nolake试验）对城市热岛强度及昆明城市气象要素水平、垂直的时空分布影响。主要结论如下:（1）四个试验城市热岛强度的平均日变化趋势相似,白天城市热岛强度较弱、夜间较强,在20时（北京时,下同）左右达到最大值。城市冠层（湖泊）对城市热岛有较明显的减（增）温,Control-NoUCM（Nolake）试验中,平均日最大差值为?0.79°C（+1.07°C）。（2）从能量平衡方程分析Control-Md04试验,感热（潜热）通量的差值为+46.18（?79.71）W m?2,潜热通量释放大于感热通量的绝对值。Control-NoUCM试验中,感热（潜热）通量的差值为?40.88（+29.60）W m?2;因NoUCM试验未考虑几何建筑物储热与遮挡,太阳辐射大部分被地表所吸收,导致感热通量偏大。（3）四种试验中,15（07）时边界层高度达到最大（小）值。NoUCM（Nolake）试验中城市边界层高度分别降低103 m（32 m）左右,而Md04试验中城市边界层高度增加102 m左右。（4）湖泊（滇池）对城市热岛环流影响的试验表明,湖泊上空垂直运动较弱,但水平方向湖陆风较大,这有利于向城市输送水汽,增加干空气湿度,使城市中空气的水汽含量增加,同时增大潜热能量释放,降低感热通量,减小了垂直温度梯度。