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.     

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.     

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.     

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.     

A closure approach with multi-scale mixing length,and its application to wind and shear stress profiles within the atmospheric boundary layer

A method to determine flow specific first-order closure for the turbulent flux of momentum in the atmospheric boundary layer (ABL) is presented. This is based on the premise that eddy viscosity is a flow rather than a fluid property, and the physically more realistic assumption that the transfer of momentum and other scalar quantities in a turbulent flow takes place by a large, but finite number of length scales, than the often used single length scale, the ‘mixing length’. The resulting eddy viscosity is flow specific and when applied to the study of the ABL, yields the vertical profiles of shear stress and mean wind velocity in good agreement with observations. The method may be extended to other types of turbulent flows, however it should be recognized that each type of flow may yield a different eddy viscosity profile. Using the derived eddy viscosity the paper presents simple analytical solutions of the ABL equations to determine observationally consistent wind speed and shear stress profiles in the ABL for a variety of practical applications including air pollution modelling.     

我国中亚热带人工林水热通量对干旱差异响应的模拟研究

季节性干旱现象在我国中亚热带地区时有发生,为了研究该区域大气-生态系统之间的相互作用关系及其碳水收支状况,2002年起在江西省千烟洲(26.7°N,115.1°E)人工林生态系统建立了通量观测塔。2003年7月该人工林生态系统遭遇了历史上少有的高温少雨天气,本研究应用基于生理生态学过程的EALCO(Ecological Assimilation of Land and Climate Observation)模型及2003和2004年通量观测数据对该生态系统的水热通量进行了模拟,同时分析了干旱胁迫对它们产生的影响。结果显示,模型能够很好的模拟该生态系统的能量通量的日变化,净辐射、显热和潜热通量模拟值与实测值相关系数的平方(R2)及标准差分别为0.99和8.05 W.m-2;0.81和41.02 W.m-2;0.90和31.49 W.m-2,模型可以解释87%的日蒸散量的变化。从模拟结果看,2003年7月下旬(发生较严重干旱胁迫)较2004年同期(干旱程度轻)相比,冠层及土壤水势下降约2倍,植物蒸腾的日变化形式改变,根系吸水滞后冠层蒸腾的时间缩短约半小时,冠层导度下降40%～60%。模拟与观测结果均表明,200...     

An introduction to flux measurements in difficult conditions

The origins of the aerodynamic techniques now widely used at sites around the world to measure continuous biosphere-atmosphere exchange of carbon and energy are briefly reviewed. A survey of the current state of this approach concludes that the technique often fails when standard analysis routines are applied to data from single towers in complex flows. In the daytime, problems are signaled by failure to close the surface energy balance because turbulent energy fluxes are routinely underestimated. Complex flows are more prevalent at night when they lead to failure to measure all the respired CO2. At such times, the aerodynamic methodology is commonly supplemented by biological models. A set of papers from a workshop on "Flux Measurements in Difficult Conditions" held at the National Center for Atmospheric Research (NCAR) in January 2006 are introduced next. Two papers review the causes and magnitude of these flow-based problems. Four papers describe intensive field experiments that detail the mechanisms that cause problematic complex flows. These experiments show, inter alia, that the technique of replacing nighttime eddy flux measurements by biological models can also be systematically biased. Finally, two model studies are used both to illustrate the physics behind these complex flows and to motivate an approach to systematic correction of single-tower results.     

Performance analysis of WRF and LES in describing the evolution and structure of the planetary boundary layer

We implemented the Weather Research and Forecast (WRF) model and WRF Large-Eddy Simulation (WRF–LES), focusing on calculations for the planetary boundary layer (PBL), and compared the results against a data set of a well-documented campaign, in the Houston–Galveston area, Texas, in summer 2006. A methodology using WRF in a mesoscale and LES was implemented to assess the performance of the model in simulating the evolution and structure of the PBL over Houston during the Vertical Mixing Experiment. Also, the WRF model in a real case mode was examined to explore potential differences between the results of each simulation approach. We analyzed both WRF results for key meteorological parameters like wind speed, wind direction and potential temperature, and compared the model results against the observations. The reasonably good agreement of LES results forced with observed surface fluxes provides confidence that LES describes turbulence quantities such as turbulent kinetic energy correctly and warrants further turbulence structure analysis. The LES results indicate a weak but noticeable nighttime turbulent kinetic energy which was produced by wind shear in Houston’s planetary boundary layer and which may likely be related to intermittent turbulence. This is supported by observations made at the University of Houston Moody Tower air quality station when intermittent peaks of carbon monoxide occurred in the evening, although the variability in wind conditions was very little.     

Flows through forest canopies in complex terrain

Recent progress on boundary layer flow within and above tall forest canopies in complex terrain is reviewed from the perspective of developing methods to interpret carbon dioxide fluxes from tower measurements in real terrain. Two examples of complex terrain are considered in detail: a forest edge, which exemplifies nonuniform forests, and hilly terrain, which can lead to drainage currents at night. Dynamical arguments show that, when boundary layer winds approach a forest edge, the mean wind adjusts on a length scale of approximately 3L(c), where L(c) is the canopy drag length scale, which depends inversely on the leaf area density of the forest. Over a further distance that also scales on L(c), turbulence in the flow adjusts, and the mixing and transport in the canopy approaches the homogeneous limit. Even low hills change the neutral flow within and above the forest canopy substantially. When the canopy is tall, pressure gradients drive flow up both the upwind and downwind slopes of the hill, leading to an ejection of air out of the top of the canopy just downwind of the crest. This flow at the crest can then advect scalar out of the top of the forest, leading to large variations in the flux of scalar across the hill. At night, when the air near the ground cools and becomes stably stratified, turbulence within the canopy can collapse, even when the flow above the canopy remains turbulent. This leads to a decoupling of the air motions within the canopy from those above. The air above the canopy can then continue to pass up and over the hill, as it does in the neutral case, but at the same time, air within the canopy drains down the hill slopes as drainage currents. These analyses will help us understand when flux towers are reliably measuring the net ecosystem exchange and suggest ways of correcting the flux tower data in more complex situations.     

Turbulence and turbulent flux events in a small estuary

Relatively little systematic research has been conducted on the turbulence characteristics of small estuaries. In the present study, detailed measurements were conducted in a small subtropical estuary with a focus on turbulent flux events. Acoustic Doppler velocimeters were installed in the mid-estuary at fixed locations and sampled simultaneously and continuously for 50 h. A turbulent flux event analysis was performed for the entire data sets extending the technique of Narasimha et al. (Phil Trans R Soc Ser A 365:841–858, 2007) to the unsteady open channel flow motion and to turbulent sub-events. Turbulent bursting events were defined in terms of the instantaneous turbulent flux. The data showed close results for all ADV units. The very-large majority of turbulent events lasted between 0.04 and 0.3 s with an average of 1 to 4 turbulent events observed per second. A number of turbulent bursting events consisted of consecutive turbulent sub-events, with between 1 and 3 sub-events per main event on average. For all ADV systems, the number of events, event duration and event amplitude showed some tidal trends, with basic differences between high- and low-water periods. A comparison between the present estuary data and the atmospheric boundary layer results of Narasimha et al. (Phil Trans R Soc Ser A 365:841–858, 2007) showed a number of similarities and demonstrated the significance of turbulent events in environmental flows. A burstiness index of 0.85 was found for the present data.     

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.     

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.     

Simulation of exchange flow between open water and floating vegetation using a modified RNG <Emphasis Type="Italic">k</Emphasis>-<Emphasis Type="Italic">ε</Emphasis> turbulence model

The paper focuses on the numerical simulation of the exchange flow between open water and floating vegetation, which plays an important role in maintaining the ecological balance by transporting nutrient matter. The simulation was conducted using a new solver developed upon OpenFOAM. A modified RNG k-ε turbulence model, which is expected to model both the high- and low-Reynolds number flows correctly, was used to determine the eddy viscosity. Several particular terms were added into the momentum equations and turbulence model equations to model the effects of vegetation and buoyancy. Among these terms, the term for the effect of vegetation in the ε-equation was re-modelled. The model was validated by properly predicting the profiles of mean velocity and turbulent kinetic energy for flows through suspended canopies. The density flow between open and vegetated water was simulated with the same conditions as those of the experiment conducted by Zhang and Nepf. The predicted results agreed well with the experimental data and provided more detailed information of such exchange flow. The convection between the root layer and the layer beneath the roots, which was not observed in the experiment, was observed in the numerical simulation.     

Large-eddy simulation of low-level jet-like flow in a canopy

The computational method of Large-Eddy Simulations has been used to study the weak, neutrally stable drainage flow within tree canopies. The computational results show that a secondary velocity maximum that resembles a jet is formed within the canopy under the nocturnal flow conditions. This jet-like flow is important in the analysis and measurements of the net ecosystem-atmosphere exchange (NEE) for carbon dioxide (CO₂). A uniformly distributed, plane source was placed within the canopy in order to simulate the nocturnal production of CO₂. The NEE is calculated as the sum of the integration of the rate of change of the concentration of CO₂ over the computational domain, the vertical turbulent flux measured directly by eddy-covariance (EC) method, and the advection terms, which are not taken into account in the EC method. Numerical results of the velocity and concentration fields, within and above the canopy, are presented and their impact on the CO₂ transport is investigated in detail. The computational results show that 15–20% of NEE is drained out by the advection process under the canopy. The results also show that the turbulent fluctuations in the lateral direction are also significant and may result in 2–5% CO₂ transport.     

Design and Testing of a Turbulence Probe for Harsh Flows

The force of wind on the ground created by turbulent eddies is commonly used to describe the horizontal flux of material during wind erosion. Here we present the Murdoch Turbulence Probe, an instrument for use in both clean and eroding flows which uses pressure differences to measure the three components of wind velocity. Correlation techniques calculate the forces near the ground and turbulence statistics in nearly real time, including turbulent velocity fluctuations from less than 0.1 Hz to 200 Hz, mean flow velocities, Reynolds stresses as well as the integral length and time scales. In the portable wind-tunnel used by Agriculture Western Australia, turbulence statistics were recorded over stable surfaces and in blowing sand from the initiation of erosion up to the time the sand supply was exhausted. Estimates of the friction velocity derived from the turbulence probe were compared with estimates obtained from the wind speed profile measured with a rake of pitot and static tubes. The Murdoch Turbulence Probe appears to work well in sandblasting conditions. Relative turbulence intensities ranged from 0.11 to 0.2 and are in close agreement with values in the literature. The ratio of the turbulence to the friction velocity (3 to 3.2) is at the high end of the reported range. The Reynolds stress measurements agree closely with predictions of the threshold friction velocities of the sand and estimates from the wind speed profile with a von Kármán constant of about 0.3, lower than the commonly accepted value of 0.4. We suggest that the wind-tunnel profile represents the `outer layer' of the boundary-layer that may best be described by a `Wake Law' or `Defect Law'. At about 54 mm above the surface, the friction velocity decreases from 0.64 m/s to 0.39 m/s and the mean velocity increases from 9.6 m/s to 11.6 m/s as the supply of sand is depleted. In addition to the friction velocity, other scales may be necessary to characterise the overriding effect of the wind and in extending wind-tunnel results to the field.     

Turbulence structure and coherent vortices in open-channel flows with wind-induced water waves

When wind-induced water waves appear over the free-surface flows such as natural rivers and artificial channels, large amounts of oxygen gas and heat are transported toward the river bed through the interface between water and wind layers. In contrast, a bed region is a kind of turbulent boundary layer, in which turbulence generation and its transport is promoted by the production of bed shear stress. In particular, coherent hairpin vortices, together with strong ejection events toward the outer part of the layer, promote mass and momentum exchanges between the inner and outer layers. It is inferred that such a near-bed turbulence may be influenced significantly by these air–water interfacial fluctuations accompanied with free-surface velocity shear and wind-induced water waves. However, these wind effects on the wall-turbulence structure are less understood. To address these exciting and challenging topics, we conducted particle imagery velocimetry (PIV) measurements in open-channel flows combined with air flows, and furthermore the present measured data allows us to investigate the effects of air–water interactions on turbulence structure through the whole depth region.     

Analytical expression for predicting the reduced settling velocity of small particles in turbulence

The influence of turbulence on the settling velocity of small particles remains an inconclusive research subject. Both enhanced and retarded particle settling compared to quiescent settling have been reported in previous literature, and several theories have been proposed. Among the mechanisms that account for reduced settling velocities in turbulence, the loitering effect is an important one that generally exists in various conditions. This study focuses on an analytical prediction of reduced particle settling velocities due to the loitering effect. By considering the velocity autocorrelation function as a step function analogous to the free path theory, and using a modified integral time scale to capture the essence of the loitering effect, an analytical expression is derived for predicting the reduction in particle settling velocities in turbulent flows. Calculation results of the expression are then examined by comparing with the results of a random walk model, direct numerical simulations that conditionally captured the reduction of particle settling velocities due to the loitering effect, and several representative experiments. Major possible influencing factors on the prediction of the analytical expression and applicable conditions of the expression are then further discussed. The proposed analytical expression is shown to be suitable for predicting the reduced settling velocities of small particles with relatively weak inertia in turbulent flows and could provide a reasonable explanation for reported cases in which small to moderate reductions in particle settling velocities were observed.

     

我国中亚热带人工林水热通量对干旱差异响应的模拟研究

季节性干旱现象在我国中亚热带地区时有发生,为了研究该区域大气-生态系统之间的相互作用关系及其碳水收支状况,2002年起在江西省千烟洲（26.7°N,115.1°E）人工林生态系统建立了通量观测塔。2003年7月该人工林生态系统遭遇了历史上少有的高温少雨天气,本研究应用基于生理生态学过程的EALCO（Ecological Assimilation of Land and Climate Observation）模型及2003和2004年通量观测数据对该生态系统的水热通量进行了模拟,同时分析了干旱胁迫对它们产生的影响。结果显示,模型能够很好的模拟该生态系统的能量通量的日变化,净辐射、显热和潜热通量模拟值与实测值相关系数的平方（R2）及标准差分别为0.99和8.05 W.m-2;0.81和41.02 W.m-2;0.90和31.49 W.m-2,模型可以解释87%的日蒸散量的变化。从模拟结果看,2003年7月下旬（发生较严重干旱胁迫）较2004年同期（干旱程度轻）相比,冠层及土壤水势下降约2倍,植物蒸腾的日变化形式改变,根系吸水滞后冠层蒸腾的时间缩短约半小时,冠层导度下降40%～60%。模拟与观测结果均表明,2003年7月下旬每天正午的波文比大都介于1～2.2,而2004年同期正午的波文比则介于0.2～0.6。EALCO模型通过Ball模型将植物碳水过程耦合在一起,从而可以很好的模拟植物的气孔行为,进而准确的模拟植物水热过程对干旱的响应。土壤水分匮乏对冠层导度的限制是2003年干旱期间冠层潜热通量模拟值下降的根本原因。     

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.