Canopy-wake dynamics and wind sheltering effects on Earth surface fluxes

The atmospheric boundary layer adjustment at the abrupt transition from a canopy (forest) to a flat surface (land or water) is investigated in a wind tunnel experiment. Detailed measurements examining the effect of canopy turbulence on flow separation, reduced surface shear stress and wake recovery are compared to data for the classical case of a solid backward-facing step. Results provide new insights into the interpretation for flux estimation by eddy-covariance and flux gradient methods and for the assessment of surface boundary conditions in turbulence models of the atmospheric boundary layer in complex landscapes and over water bodies affected by canopy wakes. The wind tunnel results indicate that the wake of a forest canopy strongly affects surface momentum flux within a distance of 35–100 times the step or canopy height, and mean turbulence quantities require distances of at least 100 times the canopy height to adjust to the new surface. The near-surface mixing length in the wake exhibits characteristic length scales of canopy flows at the canopy edge, of the flow separation in the near wake and adjusts to surface layer scaling in the far wake. Components of the momentum budget are examined individually to determine the impact of the canopy wake. The results demonstrate why a constant flux layer does not form until far downwind in the wake. An empirical model for surface shear stress distribution from a forest canopy to a clearing or lake is proposed.     

CO2 fluxes near a forest edge: a numerical study

In contrast with recent advances on the dynamics of the flow at a forest edge, few studies have considered its role on scalar transport and, in particular, on CO2 transfer. The present study addresses the influence of the abrupt roughness change on forest atmosphere CO2 exchange and contrasts the concentration and flux fields against those of a uniform forested surface. We use an atmospheric boundary layer two-equation closure model that accounts for the flow dynamics and vertical divergence of CO2 sources/sinks within a plant canopy. This paper characterizes the spatial variation of CO2 fluxes as a function of both sources/sinks distribution and the vertical structure of the canopy. Results suggest that the ground source plays a major role in the formation of wave-like vertical CO2 flux behavior downwind of a forest edge, despite the fact that the contribution of foliage sources/sinks changes monotonously. Such a variation is caused by scalar advection in the trunk space and reveals itself as a decrease or increase in vertical fluxes over the forest relative to carbon dioxide exchange of the underlying forest. The effect was more pronounced in model forests where the leaf area is concentrated in the upper part of the canopy. These results can be useful both for interpretation of existing measurements of net ecosystem exchange of CO2 (NEE) from flux towers in limited fetch conditions and in planning future CO2 transport experiments.     

The structure of turbulence near a tall forest edge: the backward-facing step flow analogy revisited

Flow disturbances near tall forest edges are receiving significant attention in diverse disciplines including ecology, forest management, meteorology, and fluid mechanics. Current theories suggest that near a forest edge, when the flow originates from a forest into a large clearing, the flow retains its forest canopy turbulence structure at the exit point. Here, we propose that this framework is not sufficiently general for dense forested edges and suggest that the flow shares several attributes with backward-facing step (BFS) flow. Similar analogies, such as rotor-like circulations, have been proposed by a number of investigators, though the consequences of such circulations on the primary terms in the mean momentum balance at the forest clearing edge have rarely been studied in the field. Using an array of three triaxial sonic anemometers positioned to measure horizontal and vertical gradients of the velocity statistics near a forest edge, we show that the flow structure is more consistent with an intermittent recirculation pattern, rather than a continuous rotor, whose genesis resembles the BFS flow. We also show that the lateral velocity variance, v'2, is the moment that adjusts most slowly with downwind distance as the flow exits from the forest into the clearing. Surprisingly, the longitudinal and vertical velocity variances (u'2 and w'2) at the forest edge were comparable in magnitude to their respective values at the center of a large grass-covered forest clearing, suggesting rapid adjustment at the edge. Discussions on how the forest edge modifies the spectra and co-spectra of momentum fluxes, effective mixing length, and static pressure are also presented.     

Advection and resulting CO2 exchange uncertainty in a tall forest in central Germany

Potential losses by advection were estimated at Hainich Forest, Thuringia, Germany, where the tower is located at a gentle slope. Three approaches were used: (1) comparing nighttime eddy covariance fluxes to an independent value of total ecosystem respiration by bottom-up modeling of the underlying processes, (2) direct measurements of a horizontal CO2 gradient and horizontal wind speed at 2 m height in order to calculate horizontal advection, and (3) direct measurements of a vertical CO2 gradient and a three-dimensional wind profile in order to calculate vertical advection. In the first approach, nighttime eddy covariance measurements were compared to independent values of total ecosystem respiration by means of bottom-up modeling of the underlying biological processes. Turbulent fluxes and storage term were normalized to the fluxes calculated by the bottom-up model. Below a u(*) threshold of 0.6 m/s the normalized turbulent fluxes decreased with decreasing u(*), but the flux to the storage increased only up to values less than 20% of the modeled flux at low turbulence. Horizontal advection was measured by a horizontal CO2 gradient over a distance of 130 m combined with horizontal wind speed measurements. Horizontal advection occurred at most of the evenings independently of friction velocity above the canopy. Nevertheless, horizontal advection was higher when u(*) was low. The peaks of horizontal advection correlated with changes in temperature. A full mass balance including turbulent fluxes, storage, and horizontal and vertical advection resulted in an increase of spikes and scatter but seemed to generally improve the results from the flux measurements. The comparison of flux data with independent bottom-up modeling results as well as the direct measurements resulted in strong indications that katabatic flows along the hill slope during evening and night reduces the measured apparent ecosystem respiration rate. In addition, anabatic flows may occur during the morning. We conclude that direct measurements of horizontal and vertical advection are highly necessary at sites located even on gentle hill slopes.     

Observations of the cross-wind velocity variance in the stable boundary layer

Nine tower datasets over grassland, brush rangeland, snow covered plain, the ocean, three different pine forests, an aspen forest and an urban site, are used to document the scale-dependence of the cross-wind velocity variance in the stable boundary layer. The turbulence velocity variance scales with the surface momentum flux, as reported in previous studies. Such scaling removes the stability dependence of the variance at a given site, and also removes most of the differences between sites. The scaling is more effective with use of a record-dependent averaging time for defining the turbulent fluctuations. The variable averaging time is the timescale associated with the gap region in the heat flux multiresolution cospectra. On scales larger than turbulence and less than a few hours (mesoscale), variations in the cross-wind velocity variance at a given site are not related to local variables such as the friction velocity. Possible exceptions include suppression of turbulence and mesoscale motions in well-defined drainage flows and enhancement of turbulence and mesoscale motions in stronger winds downstream of a ridge. Larger mesoscale variance is associated with complex terrain and forested sites compared to the more homogeneous sites in flat terrain with short or no vegetation. These differences between sites are related to the absence of a gap region in the velocity spectra at the complex terrain and forested sites. The observed probability distribution functions of the total variance and the mesoscale variance are documented for different averaging times, stability classes and site characteristics.     

Large eddy simulation of flow and scalar transport in a vegetated channel

Predicting flow and mass transport in vegetated regions has a broad range of applications in ecology and engineering practice. This paper presents large eddy simulation (LES) of turbulent flow and scalar transport within a fully developed open-channel with submerged vegetation. To properly represent the scalar transport, an additional diffusivity was introduced within the canopy to account for the contribution of stem wakes, which were not resolved by the LES, to turbulent diffusion. The LES produced good agreement with the velocity and concentration fields measured in a flume experiment. The simulation revealed a secondary flow distributed symmetrically about the channel centerline, which differed significantly from the circulation in a bare channel. The secondary circulation accelerated the vertical spread of the plume both within and above the canopy layer. Quadrant analysis was used to identify the form and shape of canopy-scale turbulent structures within and above the vegetation canopy. Within the canopy, sweep events contributed more to momentum transfer than ejection events, whereas the opposite occurred above the canopy. The coherent structures were similar to those observed in terrestrial canopies, but smaller in scale due to the constraint of the water surface.     

The role of coherent turbulent structures in explaining scalar dissimilarity within the canopy sublayer

Scalar similarity is widely assumed in models and interpretation of micro-meteorological measurements. However, in the air space within and just above the canopy (the so-called canopy sublayer, CSL) scalar similarity is generally violated. The scalar dissimilarity has been mainly attributed to differences in the distribution of sources and sinks throughout the canopy. Since large-scale coherent structures in the CSL (e.g. double roller and sweep/ejection) arise from the instabilities generated by the interaction between the mean flow and the canopy, they may encode key dynamical features about the production term responsible for the source–sink dissimilarity of scalars. Therefore, it is reasonable to assume that the geometric attributes of coherent structures are tightly coupled to the onset and the vertical extent of scalar dissimilarity within the CSL. Large-eddy simulation (LES) runs were used to investigate the role of coherent structures in explaining scalar dissimilarity among three scalars (potential air temperature, water vapour and $\text{ CO }_2$ concentration) within the CSL under near-neutral conditions for horizontally uniform but vertically varying vegetation leaf area density. It was shown that coherent structures, when identified from the first mode of a novel proper orthogonal decomposition (POD) approach, were able to capture some features of the scalar dissimilarity in the original LES field. This skill was quantified by calculating scalar–scalar correlation coefficients and turbulent Schmidt numbers of the original field and the coherent structures, respectively. However, coherent structures tend to magnify the magnitude of scalar–scalar correlation, particularly in cases where this correlation is already strong. The ability of coherent structures to describe more complex features such as the scalar sweep-ejection cycle was also explored. It was shown that the first mode of the POD does not capture the relative importance of sweeps to ejections in the original LES field. However, the superposition of few secondary coherent structures, derived from higher order POD modes, largely diminish the discrepancies between the original field and the POD expansion.     

Stratified flow interactions with a suspended canopy

Field observations of the interactions between a stratified flow and a canopy suspended from the free surface above a solid boundary are described and analysed. Data were recorded in and around the canopy formed by a large long-line mussel farm. The canopy causes a partial blockage of the water flow, reducing velocities in the upper water column. Deceleration of the approaching flow results in a deepening of isopycnals upstream of the canopy. Energy considerations show that the potential for an approaching stratified flow to be diverted beneath a porous canopy is indicated by a densimetric Froude number. Strong stratification or low-velocities inhibit vertical diversion beneath the canopy, instead favouring a horizontal diversion around the sides. The effect on vertical mixing is also considered with a shear layer generated beneath the canopy and turbulence generated from drag within the canopy. In the observations, stratification is shown to be of sufficient strength to limit the effectiveness of the first mixing process, while the turbulence within the canopy is likely to enhance vertical exchange. Velocity and temperature microstructure measurements are used to investigate the effect of the canopy on turbulent dissipation and show that dissipation is enhanced within the canopy.     

Turbulence closure models and their application in RAMS

Turbulence closures are fundamental for modelling the atmospheric diffusion in numerical codes and the resulting eddy diffusivities are key parameters in describing the transport and dispersion in the boundary layer. In this work, four turbulence closure schemes have been applied for reproducing a neutral flow over schematic complex terrain using the meteorological model RAMS. Two of the closures, a one-equation (E-l) and a two-equations (E-) model, have been implemented in RAMS in alternative to the ones originally available. In these cases, an analytical method based on the similarity theory for the atmospheric surface layer and boundary layer is adopted to calculate the empirical constants of the turbulence closures. Some examples of numerical studies performed to simulate the flow and turbulence over a 3-D hill in wind-tunnel experiment in neutral stratification are presented and discussed. An intercomparison of simulations related to different closures is considered by analysing the main features of the flow over the hill and by comparing calculated vertical profiles of turbulent kinetic energy with measured ones.     

Nocturnal accumulation of CO2 underneath a tropical forest canopy along a topographical gradient

Flux measurements of carbon dioxide and water vapor above tropical rain forests are often difficult to interpret because the terrain is usually complex. This complexity induces heterogeneity in the surface but also affects lateral movement of carbon dioxide (CO2) not readily detected by the eddy covariance systems. This study describes such variability using measurements of CO2 along vertical profiles and along a toposequence in a tropical rain forest near Manaus, Brazil. Seasonal and diurnal variation was recorded, with atmospheric CO2 concentration maxima around dawn, generally higher CO2 build-up in the dry season and stronger daytime CO2 drawdown in the wet season. This variation was reflected all along the toposequence, but the slope and valley bottom accumulated clearly more CO2 than the plateaus, depending on atmospheric stability. Particularly during stable nights, accumulation was along lines of equal altitude, suggesting that large amounts of CO2 are stored in the valleys of the landscape. Flushing of this store only occurs during mid-morning, when stored CO2 may well be partly transported back to the plateaus. It is clear that, for proper interpretation of tower fluxes in such complex and actively respiring terrain, the horizontal variability of storage needs to be taken into account not only during the night but also during the mornings.     

Using nocturnal cold air drainage flow to monitor ecosystem processes in complex terrain.

This paper presents initial investigations of a new approach to monitor ecosystem processes in complex terrain on large scales. Metabolic processes in mountainous ecosystems are poorly represented in current ecosystem monitoring campaigns because the methods used for monitoring metabolism at the ecosystem scale (e.g., eddy covariance) require flat study sites. Our goal was to investigate the potential for using nocturnal down-valley winds (cold air drainage) for monitoring ecosystem processes in mountainous terrain from two perspectives: measurements of the isotopic composition of ecosystem-respired CO2 (delta13C(ER)) and estimates of fluxes of CO2 transported in the drainage flow. To test if this approach is plausible, we monitored the wind patterns, CO2 concentrations, and the carbon isotopic composition of the air as it exited the base of a young (approximately 40 yr-old) and an old (>450 yr-old) steeply sided Douglas-fir watershed. Nocturnal cold air drainage within these watersheds was strong, deep, and occurred on more than 80% of summer nights. The depth of cold air drainage rapidly increased to tower height or greater when the net radiation at the top of the tower approached zero. The carbon isotope composition of CO2 in the drainage system holds promise as an indicator of variation in basin-scale physiological processes. Although there was little vertical variation in CO2 concentration at any point in time, we found that the range of CO2 concentration over a single evening was sufficient to estimate delta 13C(ER) from Keeling plot analyses. The seasonal variation in delta 13C(ER) followed expected trends: during the summer dry season delta 13C(ER) became less negative (more enriched in 13C), but once rain returned in the fall, delta 13C(ER) decreased. However, we found no correlation between recent weather (e.g., vapor pressure deficit) and delta 13C(ER) either concurrently or with up to a one-week lag. Preliminary estimates suggest that the nocturnal CO2 flux advecting past the 28-m tower is a rather small fraction (<20%) of the watershed-scale respiration. This study demonstrates that monitoring the isotopic composition and CO2 concentration of cold air drainage at the base of a watershed provides a new tool for quantifying ecosystem metabolism in mountainous ecosystems on the basin scale.     

Turbulence characteristics within sparse and dense canopies

Boundary layer interactions with canopies control various environmental processes. In the case of dense and homogeneous canopies, the so-called mixing layer analogy is most generally used. When the canopy becomes sparser, a transition occurs between the mixing layer and the boundary layer perturbed by interactions between element wakes. This transition has still to be fully understood and characterized. The experimental work presented here deals with the effect of the canopy density on the flow turbulence and involves an artificial canopy placed in a fully developed turbulent boundary layer. One and two-component velocity measurements are performed, both within and above the canopy. The influence of the spacing between canopy elements is studied. Longitudinal velocity statistical moments and Reynolds stresses are calculated and compared to literature data. For spacings greater than the canopy height, evidences of this transition are found in the evolution of the skewness factor, shear length scale and mixing length.     

The contribution of advective fluxes to net ecosystem exchange in a high-elevation, subalpine forest

The eddy covariance technique, which is used in the determination of net ecosystem CO2 exchange (NEE), is subject to significant errors when advection that carries CO2 in the mean flow is ignored. We measured horizontal and vertical advective CO2 fluxes at the Niwot Ridge AmeriFlux site (Colorado, USA) using a measurement approach consisting of multiple towers. We observed relatively high rates of both horizontal (F(hadv)) and vertical (F(vadv)) advective fluxes at low surface friction velocities (u(*)) which were associated with downslope katabatic flows. We observed that F(hadv) was confined to a relatively thin layer (0-6 m thick) of subcanopy air that flowed beneath the eddy covariance sensors principally at night, carrying with it respired CO2 from the soil and lower parts of the canopy. The observed F(vadv) came from above the canopy and was presumably due to the convergence of drainage flows at the tower site. The magnitudes of both F(hadv) and F(vadv) were similar, of opposite sign, and increased with decreasing u(*), meaning that they most affected estimates of the total CO2 flux on calm nights with low wind speeds. The mathematical sign, temporal variation and dependence on u(*) of both F(hadv) and F(vadv) were determined by the unique terrain of the Niwot Ridge site. Therefore, the patterns we observed may not be broadly applicable to other sites. We evaluated the influence of advection on the cumulative annual and monthly estimates of the total CO2 flux (F(c)), which is often used as an estimate of NEE, over six years using the dependence of F(hadv) and F(vadv) on u(*). When the sum of F(hadv) and F(vadv) was used to correct monthly F(c), we observed values that were different from the monthly F(c) calculated using the traditional u(*)-filter correction by--16 to 20 g C x m(-2) x mo(-1); the mean percentage difference in monthly Fc for these two methods over the six-year period was 10%. When the sum of F(hadv) and F(vadv) was used to correct annual Fc, we observed a 65% difference compared to the traditional u(*)-filter approach. Thus, the errors to the local CO2 budget, when F(hadv) and F(vadv) are ignored, can become large when compounded in cumulative fashion over long time intervals. We conclude that the "micrometeorological" (using observations of F(hadv) and F(vadv)) and "biological" (using the u(*) filter and temperature vs. F(c) relationship) corrections differ on the basis of fundamental mechanistic grounds. The micrometeorological correction is based on aerodynamic mechanisms and shows no correlation to drivers of biological activity. Conversely, the biological correction is based on climatic responses of organisms and has no physical connection to aerodynamic processes. In those cases where they impose corrections of similar magnitude on the cumulative F(c) sum, the result is due to a serendipitous similarity in scale but has no clear mechanistic explanation.     

The assessment of environmentally sensitive forest road construction in Calabrian pine forest areas of Turkey

Forest road construction by bulldozers in Calabrian Pine (Pinus brutia Ten.) forests on mountainous terrain of Turkey causes considerable damage to the environment and the forest standing alongside the road. This situation obliges a study of environmentally sound road construction in Turkey. This study was carried out in 4 sample sites of Antalya Forest Directorate in steep (34-50% gradient) and very steep terrain (51-70% gradient) conditions with bulldozer and excavator machine and direct damages to forest during road construction was determined, including forest area losses and damages to downhill trees in mountainous areas. It was determined that in steep terrain when excavators were used, less forest area (22.16%) was destroyed compared to bulldozers and 26.54% less area in very steep terrain. The proportion of damage on trees where bulldozer worked was nearly twofold higher than excavator was used. The results of this research show that the environmentally sensitive techniques applied for the road construction projects are considerably superior to the traditional use of bulldozers on steep slopes. The environmentally sound forest road construction by use of excavator must be considered an appropriate and reliable solution for mountainous terrain where areas of sensitive forest ecosystems are to be opened up.     

Instantaneous transport of a passive scalar in a turbulent separated flow

The results of large-eddy simulations of flow and transient solute transport over a backward facing step and through a 180° bend are presented. The simulations are validated successfully in terms of hydrodynamics and tracer transport with experimental velocity data and measured residence time distribution curves confirming the accuracy of the method. The hydrodynamics are characterised by flow separation and subsequent recirculation in vertical and horizontal directions and the solute dispersion process is a direct response to the significant unsteadiness and turbulence in the flow. The turbulence in the system is analysed and quantified in terms of power density spectra and covariance of velocity fluctuations. The injection of an instantaneous passive tracer and its dispersion through the system is simulated. Large-eddy simulations enable the resolution of the instantaneous flow field and it is demonstrated that the instabilities of intermittent large-scale structures play a distinguished role in the solute transport. The advection and diffusion of the scalar is governed by the severe unsteadiness of the flow and this is visualised and quantified. The analysis of the scalar mass transport budget quantifies the mechanisms controlling the turbulent mixing and reveals that the mass flux is dominated by advection.     

Quantifying vertical mixing in estuaries

Estuarine turbulence is notable in that both the dissipation rate and the buoyancy frequency extend to much higher values than in other natural environments. The high dissipation rates lead to a distinct inertial subrange in the velocity and scalar spectra, which can be exploited for quantifying the turbulence quantities. However, high buoyancy frequencies lead to small Ozmidov scales, which require high sampling rates and small spatial aperture to resolve the turbulent fluxes. A set of observations in a highly stratified estuary demonstrate the effectiveness of a vessel-mounted turbulence array for resolving turbulent processes, and for relating the turbulence to the forcing by the Reynolds-averaged flow. The observations focus on the ebb, when most of the buoyancy flux occurs. Three stages of mixing are observed: (1) intermittent and localized but intense shear instability during the early ebb; (2) continuous and relatively homogeneous shear-induced mixing during the mid-ebb, and weakly stratified, boundary-layer mixing during the late ebb. The mixing efficiency as quantified by the flux Richardson number Rf was frequently observed to be higher than the canonical value of 0.15 from Osborn (J Phys Oceanogr 10:83–89, 1980). The high efficiency may be linked to the temporal–spatial evolution of shear instabilities.     

Modeling seed dispersal distances: implications for transgenic Pinus taeda.

Predicting forest-tree seed dispersal across a landscape is useful for estimating gene flow from genetically engineered (GE) or transgenic trees. The question of biocontainment has yet to be resolved, although field-trial permits for transgenic forest trees are on the rise. Most current field trials in the United States occur in the Southeast where Pinus taeda L., an indigenous species, is the major timber commodity. Seed dispersal distances were simulated using a model where the major determinants were: (1) forest canopy height at seed release, (2) terminal velocity of the seeds, (3) absolute seed release, and (4) turbulent-flow statistics, all of which were measured or determined within a P. taeda plantation established from seeds collected from wild forest-tree stands at the Duke Forest near Durham, North Carolina, USA. In plantations aged 16 and 25 years our model results showed that most of the seeds fell within local-neighborhood dispersal distances, with estimates ranging from 0.05 to 0.14 km from the source. A fraction of seeds was uplifted above the forest canopy and moved via the long-distance dispersal (LDD) process as far as 11.9-33.7 km. Out of 10(5) seeds produced per hectare per year, roughly 440 seeds were predicted to be uplifted by vertical eddies above the forest canopy and transported via LDD. Of these, 70 seeds/ha traveled distances in excess of 1 km from the source, a distance too great to serve as a biocontainment zone. The probability of LDD occurrence of transgenic conifer seeds at distances exceeding 1 km approached 100%.     

Modelling hydraulic jump using the bubbly two-phase flow method

Hydraulic jumps have complex flow structures, characterised by strong turbulence and large air contents. It is difficult to numerically predict the flows. It is necessary to bolster the existing computer models to emphasise the gas phase in hydraulic jumps, and avoid the pitfall of treating the phenomenon as a single-phase water flow. This paper aims to improve predictions of hydraulic jumps as bubbly two-phase flow. We allow for airflow above the free surface and air mass entrained across it. We use the Reynolds-averaged Navier–Stokes equations to describe fluid motion, the volume of fluid method to track the interface, and the k–ε model for turbulence closure. A shear layer is shown to form between the bottom jet flow and the upper recirculation flow. The key to success in predicting the jet flow lies in formulating appropriate bottom boundary conditions. The majority of entrained air bubbles are advected downstream through the shear layer. Predictions of the recirculation region’s length and air volume fraction within the layer are validated by available measurements. The predictions show a linear growth of the shear layer. There is strong turbulence at the impingement, and the bulk of the turbulence kinetic energy is advected to the recirculation region via the shear layer. The predicted bottom-shear-stress distribution, with a peak value upstream of the toe of the jump and a decaying trend downstream, is realistic. This paper reveals a significant transient bottom shear stress associated with temporal fluctuations of mainly flow velocity in the jump. The prediction method discussed is useful for modelling hydraulic jumps and advancing the understanding of the complex flow phenomenon.     

城市化对北京夏季极端高温影响的数值研究

利用一个耦合了城市冠层模式（UCM）的区域数值模拟系统（WRF／NCAR）,引入由LandsatTM提取的京津冀区域30m分辨率下垫面GIS数据集代替美国USGS地表分类数据,对2009年6月24—25日出现在北京地区的一次超过40℃极端高温天气过程进行了高分辨率数值模拟,用以考察WRF／UCM系统对北京“城市热岛”及城市高温天气的模拟效果,并分析了城市化对北京地区城市高温和地表能量平衡的影响及其日变化特征。结果表明：采用精细化下垫面分类数据集后能更好地模拟出主要高温区的分布特征,并能较好再现夜间的“城市热岛”效应。城市化发展对近地层气温的影响主要表现在会促使城区及其下风向近郊区气温的升高,增幅可达0．5～2℃,这与城市热岛及其下游效应有关。城市下垫面的高粗糙度对近地层风速表现出明显的阻挡效应,表现在模拟的10m风场减弱明显。考虑了城市下垫面属性的敏感性试验更好地再现了城区温度的日变化,其模拟的日间最高温度与实际观测值更为接近,也较好地模拟出了城区具有较高最低温度的特征。通过城区与郊区能量平衡过程差异的分析表明,城市化可以显著改变能量平衡中各项所占的比重。地表对近地层大气热量输送过程的变化表明随着城市下垫面的日愈扩大,会显著增强白天地表对大气的向上感热输送,增大城区日间出现高温的可能。夜间,模式反映出地表能量收入来自土壤热通量的向上输送,同样由于城区的潜热通量小,收入的能量仍主要以感热形式加热大气,夜间城区具有较高的最低温度并表现出较强的热岛特征,主要与夜间感热加热的持续相关。     

Turbulent kinetic energy budgets in a model canopy: comparisons between LES and wind-tunnel experiments

A comparative study of turbulence in a wind-tunnel model canopy is performed, using Large eddy simulation (LES) and experimental data from PIV and hot-wire anemometry measurements. The model canopy is composed of thin cylindrical stalks. In the LES, these are represented using a plant-scale approach, while the scale-dependent Lagrangian dynamic model is used as subgrid-scale model. LES predictions of turbulence statistics and energy spectra are found to be in good agreement with the experimental data. Turbulent kinetic energy (TKE) budgets from the LES simulation are analyzed to provide more information absent in the measurements. Results confirm that sloshing motions at the low levels of the canopy are mainly driven by pressure fluctuations. A difference between the energy flux obtained from the energy spectrum and the SGS dissipation rate is observed, consistent with a spectral bypass mechanism in which the real spectral flux due to cascade is smaller than that implied by the energy-spectrum level, due to direct drain by the canopy.