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.     

Embedded large eddy simulation approach for pollutant dispersion around a model building in atmospheric boundary layer

In the present article, the potential of embedded large eddy simulation (ELES) approach to reliably predict pollutant dispersion around a model building in atmospheric boundary layer is assessed. The performance of ELES in comparison with large eddy simulation (LES) is evaluated in several ways. These include a number of qualitative and quantitative comparisons of time-averaged and instantaneous results with wind tunnel measurements supplemented by statistical data analyses using scatter plots and standard evaluation metrics. Results obtained by both LES and ELES approaches show very good agreement with the experiment. However, addition of turbulence to mean flow at Reynolds averaged Navier–Stokes (RANS)–LES interface in ELES approach not only increases the turbulence intensity, it also results in larger values of turbulent kinetic energy (TKE) as well as a shorter reattachment length in the wake region. Accordingly, higher levels of TKE predicted by ELES increase the local intensity of concentration leading to shorter plume shapes as compared with LES. In general, ELES shows better agreement with experiment on the surfaces of model building and also in the downstream wake region. In terms of computational costs, the CPU time required to obtain statistical values in ELES is about 49 % lower than that of LES and the number of iterations per time step is also reduced by 55 % as compared with LES.     

LES validation of urban flow,part I: flow statistics and frequency distributions

Essential prerequisites for a thorough model evaluation are the availability of problem-specific, quality-controlled reference data and the use of model-specific comparison methods. The work presented here is motivated by the striking lack of proportion between the increasing use of large-eddy simulation (LES) as a standard technique in micro-meteorology and wind engineering and the level of scrutiny that is commonly applied to assess the quality of results obtained. We propose and apply an in-depth, multi-level validation concept that is specifically targeted at the time-dependency of mechanically induced shear-layer turbulence. Near-surface isothermal turbulent flow in a densely built-up city serves as the test scenario for the approach. High-resolution LES data are evaluated based on a comprehensive database of boundary-layer wind-tunnel measurements. From an exploratory data analysis of mean flow and turbulence statistics, a high level of agreement between simulation and experiment is apparent. Inspecting frequency distributions of the underlying instantaneous data proves to be necessary for a more rigorous assessment of the overall prediction quality. From velocity histograms local accuracy limitations due to a comparatively coarse building representation as well as particular strengths of the model to capture complex urban flow features with sufficient accuracy are readily determined. However, the analysis shows that further crucial information about the physical validity of the LES needs to be obtained through the comparison of eddy statistics, which is focused on in part II. Compared with methods that rely on single figures of merit, the multi-level validation strategy presented here supports conclusions about the simulation quality and the model’s fitness for its intended range of application through a deeper understanding of the unsteady structure of the flow.     

LES validation of urban flow,part II: eddy statistics and flow structures

Time-dependent three-dimensional numerical simulations such as large-eddy simulation (LES) play an important role in fundamental research and practical applications in meteorology and wind engineering. Whether these simulations provide a sufficiently accurate picture of the time-dependent structure of the flow, however, is often not determined in enough detail. We propose an application-specific validation procedure for LES that focuses on the time dependent nature of mechanically induced shear-layer turbulence to derive information about strengths and limitations of the model. The validation procedure is tested for LES of turbulent flow in a complex city, for which reference data from wind-tunnel experiments are available. An initial comparison of mean flow statistics and frequency distributions was presented in part I. Part II focuses on comparing eddy statistics and flow structures. Analyses of integral time scales and auto-spectral energy densities show that the tested LES reproduces the temporal characteristics of energy-dominant and flux-carrying eddies accurately. Quadrant analysis of the vertical turbulent momentum flux reveals strong similarities between instantaneous ejection-sweep patterns in the LES and the laboratory flow, also showing comparable occurrence statistics of rare but strong flux events. A further comparison of wavelet-coefficient frequency distributions and associated high-order statistics reveals a strong agreement of location-dependent intermittency patterns induced by resolved eddies in the energy-production range. The validation concept enables wide-ranging conclusions to be drawn about the skill of turbulence-resolving simulations than the traditional approach of comparing only mean flow and turbulence statistics. Based on the accuracy levels determined, it can be stated that the tested LES is sufficiently accurate for its purpose of generating realistic urban wind fields that can be used to drive simpler dispersion models.     

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.     

Simulation of Ekman Boundary Layers by Large Eddy Model with Dynamic Mixed Subfilter Closure

Theoretical analysis of boundary layer turbulence has suggested a feasibility of sufficiently accurate turbulence resolving simulations at relatively coarse meshes. However, large eddy simulation (LES) codes, which employ traditional eddy-viscosity turbulence closures, fail to provide adequate turbulence statistics at coarse meshes especially within a surface layer. Manual tuning of parameters in these turbulence closures may correct low order turbulence statistics but severely harms spectra of turbulence kinetic energy (TKE). For more than decade, engineering LES codes successfully employ dynamic turbulence closures. A dynamic Smagorinsky turbulence closure (DSM) has been already tried in environmental LES. The DSM is able to provide adequate turbulence statistics at coarse meshes but it is not completely consistent with the LES equations. This paper investigates applicability of an advanced dynamic mixed turbulence closure (DMM) to simulations of Ekman boundary layers of high Reynolds number flows. The DMM differs from the DSM by explicit calculation of the Leonard term in the turbulence stress tensor. The Horizontal Array Turbulence Study (HATS) field program has revealed that the Leonard term is indeed an important component of the real turbulence stress tensor. This paper presents validation of a new LES code LESNIC. The study shows that the LES code with the DMM provides rather accurate low order turbulence statistics and the TKE spectra at very coarse meshes. These coarse LES maintain more energetic small scale fluctuations of velocity especially within the surface layer. This is critically important for success of simulations. Accurate representation of higher order turbulence statistics, however, requires essentially better LES resolution. The study also shows that LES of the Ekman boundary layer cannot be directly compared with conventionally neutral atmospheric boundary layers. The depth of the boundary layer is an important scaling parameter for turbulence statistics.     

Numerical simulation of water mixing and renewals in the Barcelona harbour area: the winter season

In the present paper, we use numerical simulation to investigate currents, mixing and water renewal in Barcelona harbour under typical conditions of wind forcing for the winter season. This site is of particular importance due to the interplay between touristic and commercial activities, requiring detailed and high-definition studies of water quality within the harbour. We use Large Eddy Simulation (LES) which directly resolves the anisotropic and energetic large scales of motion and parametrizes the small, dissipative, ones. Small-scale turbulence is modelled by the anisotropic Smagorinsky model (ASM) to be employed in presence of large cell anisotropy. The complexity of the harbour is modelled using a combination of curvilinear, structured, non-staggered grid and the immersed boundary method. Boundary conditions for wind and currents at the inlets of the port are obtained from in-situ measurements. Analysis of the numerical results is carried out based on both instantaneous and time-averaged velocity fields. First- and second-order statistics, such as turbulent kinetic energy and horizontal and vertical eddy viscosities, are calculated and their spatial distribution is discussed. The study shows the presence of intense current in the narrow and elongated part of the harbour together with sub-surface along-shore elongated rolling structures (with a time scale of a few hours), and they contribute to the vertical water mixing. Time-averaged velocity field reveals intense upwelling and downwelling zones along the walls of the harbour. The analysis of second-order statistics shows strong inhomogeneity of turbulent kinetic energy and horizontal and vertical eddy viscosities in the horizontal plane, with larger values in the regions characterized by stronger currents. The water renewal within the port is quantified for particular sub-domain regions, showing that the complexity of the harbour is such that certain in-harbour basins have a water renewal of over five days, including the yacht marina area. The LES solution compares favourably with available current-meter data. The LES solution is also compared with a RANS solution obtained in literature for the same site under the same forcing conditions, the comparison demonstrating a large sensitivity of properties to model resolution and frictional parametrization.     

A Large-Eddy Simulation Study of the Convective Boundary Layer over Philadelphia during the 1999 Summer NE-OPS Campaign

This study presents a large-eddy simulation (LES) study of the convective boundary layer on August 1, 1999 over Philadelphia, PA during a summer ozone episode. The study is an evaluation of the Colorado State University's Regional Atmospheric Modeling System Version 4.3 (RAMS4.3) with the LES option using Northeast Oxidant and Particulate Study (NE-OPS) data. Simulations were performed with different imposed sensible heat fluxes at the ground surface. The model was initialized with the atmospheric sounding data collected at Philadelphia at 1230 UTC and model integrations continued till 2130 UTC. The resulting mean profiles of temperature and humidity obtained from the LES model were compared with atmospheric soundings, tethered balloon and aircraft data collected during the NE-OPS 1999 field campaign. Also the model-derived vertical profiles of virtual temperature were compared with NE-OPS Radio Acoustic Sounder System (RASS) data while the humidity profiles were compared with NE-OPS lidar data. The comparison of the radiosonde data with the LES model predictions suggests that the growth of the mixing layer is reasonably well simulated by the model. Overall, the agreement of temperature predictions of the LES model with the radiosonde observations is good. The model appears to underestimate humidity values for the case of higher imposed sensible heat flux. However, the humidity values in the mixing layer agree quite well with radiosonde observations for the case of lower imposed sensible heat flux. The model-predicted temperature and humidity profiles are in reasonable agreement with the tethered balloon data except for some small overestimation of temperature at lower layers and some underestimation of humidity values. However, the humidity profiles as simulated by the model agree quite well with the tethered balloon data for the case of lower imposed sensible heat flux. The model-predicted virtual temperature profile is also in better agreement with RASS data for the case of lower imposed sensible heat flux. The model-predicted temperature profile further agrees quite well with aircraft data for the case of lower imposed heat flux. However, the relative humidity values predicted by the model are lower compared with the aircraft data. The model-predicted humidity profiles are only in partial agreement with the lidar data. The results of this study suggest that the explicitly resolved energetic eddies seem to provide the correct forcing necessary to produce good agreement with observations for the case of an imposed sensible heat flux of 0.1 K m s^–1 at the surface.     

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.     

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.     

Large-eddy simulation of starting buoyant jets

A series of Large Eddy Simulations (LES) are performed to investigate the penetration of starting buoyant jets. The LES code is first validated by comparing simulation results with existing experimental data for both steady and starting pure jets and lazy plumes. The centerline decay and the growth rate of the velocity and concentration fields for steady jets and plumes, as well as the simulated transient penetration rate of a starting pure jet and a starting lazy plume, are found to compare well with the experiments. After validation, the LES code is used to study the penetration of starting buoyant jets with three different Reynolds numbers from 2000 to 3000, and with a wide range of buoyancy fluxes from pure jets to lazy plumes. The penetration rate is found to increase with an increasing buoyancy flux. It is also observed that, in the initial Period of Flow Development, the two penetrative mechanisms driven by the initial buoyancy and momentum fluxes are uncoupled; therefore the total penetration rate can be resolved as the linear addition of these two effects. A fitting equation is proposed to predict the penetration rate by combining the two independent mechanisms.     

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

季节性干旱现象在我国中亚热带地区时有发生,为了研究该区域大气-生态系统之间的相互作用关系及其碳水收支状况,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...     

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.     

Measurements and modeling of open-channel flows with finite semi-rigid vegetation patches

The hydrodynamics of flows through a finite length semi-rigid vegetation patch (VP) were investigated experimentally and numerically. Detailed measurements have been carried out to determine the spatial variation of velocity and turbulence profiles within the VP. The measurement results show that an intrusion region exists in which the peak Reynolds stress remains near the bed. The velocity profile is invariant within the downstream part of the VP while the Reynolds stress profile requires a longer distance to attain the spatially invariant state. Higher vegetation density leads to a shorter adjustment length of the transition region, and a higher turbulence level within the VP. The vegetation density used in the present study permits the passing through of water and causes the peak Reynolds stress and turbulence kinetic energy each the maximum at the downstream end of the patch. A 3D Reynolds-averaged Navier–Stokes model incorporating the Spalart–Allmaras turbulence closure was employed subsequently to replicate the flow development within the VP. The model reproduced transitional flow characteristics well and the results are in good agreement with the experimental data. Additional numerical experiments show that the adjustment length can be scaled by the water depth, mean velocity and maximum shear stress. Empirical equations of the adjustment lengths for mean velocity and Reynolds stress were derived with coefficients quantified from the numerical simulation results.     

A Fine-Scale <Emphasis Type="Italic">k</Emphasis>−ε Model for Atmospheric Flow over Heterogeneous Landscapes

A multi-purpose model for small-scale atmospheric flows over heterogeneous landscapes is being developed. The aim of this research is to build a tool able to predict the dynamical (wind, turbulence) and diffusive (gases, particles) fields over landscapes characterised by heterogeneous plant cover. In its present stage of development the model is based on the numerical integration of neutral atmospheric flow equations, using an energy-dissipation closure scheme and over a domain that may include vegetation layers. Three validation cases of the model are presented: (i) response of the airflow to a change in surface roughness; (ii) airflow within and above a horizontally homogeneous plant canopy; (iii) airflow over two complex forest-to-clearing and clearing-to-forest transitions. All simulations provide results in good agreement with the experimental data, except for turbulent kinetic energy just after a clearing-to-forest transition. This result is not surprising for a statistical k– model in a flow region characterised by strong distorsion and intermittent turbulence. However the overall good performance of the model is promising for environmental research at fine scales over heterogeneous landscapes.     

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

季节性干旱现象在我国中亚热带地区时有发生,为了研究该区域大气-生态系统之间的相互作用关系及其碳水收支状况,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年干旱期间冠层潜热通量模拟值下降的根本原因。     

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.     

Ozon in Waldökosystemen aus atmosphärenchemischer und pflanzenphysiologischer Sicht

The trend of rising ozone concentrations in forest ecosystems and the phytotoxicity of ozone demand a realistic risk assessment according to an internationally accepted and flux-based quality standard. Ozone fluxes within the canopy are influenced by chemical gas-phase reactions with nitrogen oxide and biogenic hydrocarbons and by surface deposition processes. Therefore, a differentiation of the ozone flux within the canopy is needed between stomatal uptake and other transport pathways. The Eddy Covariance technique is the method of choice for the determination of trace gas fluxes relevant for ozone chemistry. This method is also used for stomatal conductance measurements based on evapotranspiration fluxes and for emission measurements of biogenic hydrocarbons by PTR-MS. Although considerably research efforts were directed to canopy measurements in recent years, the underlying processes are not fully understood yet. Thus, major differences occur in the ratios of stomatal ozone uptake, non-stomatal deposition and gas-phase chemistry between different studies. Furthermore, the vertical concentration gradients within the canopy measured at several forest sites are rather inconsistent and the existing deposition models do rarely account for chemical transformation and detoxification processes. Only a simultaneous measurement of all photochemically relevant trace gases, plant physiological parameters at different sites and forest species over entire vegetation periods, and model parameterization according to the measurement results from the experimental sites will contribute to the clarification of the canopy processes and will ensure realistic risk assessments.     

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.