Assessing the shear-sheltering theory applied to low-level jets in the nocturnal stable boundary layer

This paper investigates the existence of shear sheltering on turbulence data over a quasi-ideal experimental site in Oklahoma, USA. Originally developed for engineering flows, the shear-sheltering theory is predicated upon the idea of low-level jets blocking large eddies aloft, preventing them from propagating to the surface. In this scenario, suppression of low-frequency turbulence energy and reduction of surface fluxes would be expected. Results from the Oklahoma experiment show instead an enhancement of surface turbulence intensity and of the relative contribution of large scales to total (co)variances for low-level jet cases with strong shear, thus suggesting the absence of shear sheltering at the site. The results underline the complexity of surface-atmosphere interactions in nocturnal stable conditions. Atmospheric modeling of exchange using various scenarios of surface characteristics, flow regimes, and low-level jet properties is suggested to further assess the potential applicability of the shear-sheltering theory to atmospheric flows.     

Coherent eddies and turbulence in vegetation canopies: The mixing-layer analogy

This paper argues that the active turbulence and coherent motions near the top of a vegetation canopy are patterned on a plane mixing layer, because of instabilities associated with the characteristic strong inflection in the mean velocity profile. Mixing-layer turbulence, formed around the inflectional mean velocity profile which develops between two coflowing streams of different velocities, differs in several ways from turbulence in a surface layer. Through these differences, the mixing-layer analogy provides an explanation for many of the observed distinctive features of canopy turbulence. These include: (a) ratios between components of the Reynolds stress tensor; (b) the ratio K _H/K _M of the eddy diffusivities for heat and momentum; (c) the relative roles of ejections and sweeps; (d) the behaviour of the turbulent energy balance, particularly the major role of turbulent transport; and (e) the behaviour of the turbulent length scales of the active coherent motions (the dominant eddies responsible for vertical transfer near the top of the canopy). It is predicted that these length scales are controlled by the shear length scale % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamitamaaBa% aaleaacaWGtbaabeaakiabg2da9iaadwfacaGGOaGaamiAaiaacMca% caGGVaGabmyvayaafaGaaiikaiaadIgacaGGPaaaaa!3FD0!\[L_S = U(h)/U'(h)\] (where h is canopy height, U(z) is mean velocity as a function of height z, and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmyvayaafa% Gaeyypa0JaaeizaiaadwfacaGGVaGaaeizaiaadQhaaaa!3C32!\[U' = {\rm{d}}U/{\rm{d}}z\]). In particular, the streamwise spacing of the dominant canopy eddies is _x = mL _s, with m = 8.1. These predictions are tested against many sets of field and wind-tunnel data. We propose a picture of canopy turbulence in which eddies associated with inflectional instabilities are modulated by larger-scale, inactive turbulence, which is quasi-horizontal on the scale of the canopy.     

The Effect of the Vertical Source Distribution on Scalar Statistics within and above a Forest Canopy

Little is known about in-canopy processes that may alter forest–atmosphere exchanges of trace gases and aerosols. To improve our understanding of in-canopy mixing, we use large-eddy simulation to study the effect of scalar source/sink distributions on scalar concentration moments, fluxes, and correlation coefficients within and above an ideal forest canopy. Scalars are emitted from: (1) the ground, (2) the canopy, and (3) both the ground and the canopy; a scalar is also deposited onto the canopy. All scalar concentration moments, fluxes, and correlation coefficients are affected by the source location/distribution, as is the scalar segregation intensity. We conclude that vertical source/sink distribution has a profound impact on scalar concentration profiles, fluxes, correlation coefficient, and scalar segregation.     

Impact of Eddy Characteristics on Turbulent Heat and Momentum Fluxes in the Urban Roughness Sublayer

Eddy-covariance observations above the densely built-up Centre of Nanjing were made from December 2011 to August 2012. Separate eddy-covariance systems installed at two levels on a 36-m tower located on a rooftop were operated simultaneously, and observations grouped into two sectors (A, B) according to the prevalent wind directions. For sector A, where the nearby buildings are all below the lower measurement level, the sensible heat and momentum fluxes are generally greater at the upper level. For sector B, where several high-rise buildings are located upwind, the sensible heat and momentum fluxes at the upper level are close to those at the lower level. The analysis shows that the turbulent eddy characteristics differ between the two wind sectors, leading to a different behaviour of turbulent exchange between the two levels. A hypothesis is proposed that addresses the vertical variation of turbulent fluxes in the urban roughness sublayer (RSL). For sector A, the buildings block the flow, change the trajectory of scalars, and distort the footprint of scalar fluxes; this ‘blocking effect’ is believed to lead to a smaller sensible heat flux above the canopy layer. Such an effect should decrease with height in the RSL, explaining the increase of the observed turbulent heat flux with height. In addition, the presence of non-uniform building heights adversely affects turbulence organization around the canopy top, and likely elevates the inflection point of the mean flow to a higher elevation close to the upper measurement level, where larger shear results in a larger momentum flux. For sector B, wake effects from the nearby high-rise buildings strongly reduce turbulence organization at higher elevations, leading to similar sensible heat and momentum fluxes at both measurement levels.     

Vertical Heat Transfer in the Lower Atmosphere over the Arctic Ocean During Clear-sky Periods

A diagnostic study of heat transfer within the lower atmosphere and between the atmosphere and the surface of the Arctic Ocean snow/ice pack during clear-sky conditions is conducted using data from the Surface Heat Budget of the Arctic Ocean (SHEBA) field experiment. Surface heat budgets computed for four cloudy and four clear periods show that, while the net turbulent heat fluxes at the surface are small during the cloudy periods, during the clear-sky periods they are a considerable source of surface heating, balancing significant portions of the conductive heat fluxes from within the snow/ice pack. Analysis of the dynamics and thermodynamics of the lower atmosphere during the clear-sky periods reveals that a considerable portion of the heat lost to the surface by turbulent heat fluxes is balanced by locally strong heating near the atmospheric boundary-layer (ABL) top due to the interaction of subsiding motions with the strong overlying temperature inversions surmounting the ABL. This heat is then entrained into the ABL and transported to the surface by turbulent mixing, maintained by a combination of vertical wind shear and wave-turbulence interactions. The frequency of stable, clear-sky periods, particularly during the winter, combined with these results, suggests that the downward transfer of heat through the lower atmosphere and into the surface represents an important component of the heat budgets of the lower atmosphere and snow/ice pack over the annual cycle     

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

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

Spectral Structure of Small-Scale Turbulent and Mesoscale Fluxes in the Atmospheric Boundary Layer over a Thermally Inhomogeneous Land Surface

Spectral analysis was performed on aircraft observations of a convective boundary layer (CBL) that developed over a thermally inhomogeneous, well-marked mesoscale land surface. The observations, part of the GAME-Siberia experiment, were recorded between April and June 2000 over the Lena River near Yakutsk City. A special integral parameter termed the ‘reduced depth of the CBL’ was used to scale the height of the mixed layer with variable depth. Analysis of wavelet cospectra and spectra facilitated the separation of fluxes and other variables into small-scale turbulent fluctuations (with scales less than the reduced depth of the CBL, approximately 2 km) and mesoscale fluctuations (up to 20 km). This separation approach allows for independent exploration of the scales. Analyses showed that vertical distributions obeyed different laws for small-scale fluxes and mesoscale fluxes (of sensible heat, water vapour, momentum and carbon dioxide) and for other variables (wind speed and air temperature fluctuations, coherence and degree of anisotropy). Vertical profiles of small-scale turbulent fluxes showed a strong decay that differed from generally accepted similarity models for the CBL. Vertical profiles of mesoscale fluxes and other variables clearly showed sharp inflections at the same relative (with respect to the reduced depth of the CBL) height of approximately 0.55 in the CBL. Conventional similarity models for sensible heat fluxes describe both small-scale turbulent and mesoscale flows. The present results suggest that mesoscale motions that reach up to the relative level of 0.55 could be initiated by thermal surface heterogeneity. Entrainment between the upper part of the CBL and the free atmosphere may cause mesoscale motions in that region of the CBL.     

Seasonal and diurnal variations of coherent structures over a deciduous forest

Coherent structures in turbulent flow above a midlatitude deciduous forest are identified using a wavelet analysis technique. Coupling between motions above the canopy (z/h=1.5, whereh is canopy height) and within the canopy (z/h=0.6) are studied using composite velocity and temperature fields constructed from 85 hours of data. Data are classified into winter and summer cases, for both convective and stable conditions. Vertical velocity fluctuations are in phase at both observation levels. Horizontal motions associated with the structures within the canopy lead those above the canopy, and linear analysis indicates that the horizontal motions deep in the canopy should lead the vertical motions by 90°. On average, coherent structures are responsible for only about 40% of overall turbulent heat and momentum fluxes, much less than previously reported. However, our large data set reveals that this flux fraction comes from a wide distribution that includes much higher fractions in its upper extremes. The separation distanceL _s between adjacent coherent structures, 6–10h, is comparable to that obtained in previous observations over short canopies and in the laboratory. Changes in separation between the summer and winter (leafless) conditions are consistent withL _s being determined by a local horizontal wind shear scale.     

A Second-Order Closure for Neutrally Stratified Vegetative Canopy Flows

An existing second-order closure model is modified to include the effects on mean and turbulent motions of form and viscous drag in vegetative canopies. The additional physical mechanisms represented by the closure are viscous and pressure drag on canopy elements, their role in momentum absorption, in the creation of fine scale turbulent eddies and in enhancing the total viscous dissipation in the canopy airspace. Viscous dissipation is split into a standard 'isotropic contribution associated with the spectral eddy cascade and a foliage contribution associated with work against pressure and viscous drag on the foliage. Changes in the turbulent time scale that result from these mechanisms are included in the standard parameterisations of third moments and of the eddy cascade contribution to dissipation. The model is tested against a wind- tunnel 'wheat canopy, a corn canopy and a eucalypt canopy, a height range from 50 mm to 12.6 m. Model results show that the parameterisations of foliage interaction used in the closure are sufficiently robust to reproduce second-moment profiles within and above vegetative canopies to a high degree of accuracy without resorting to 'tuning of the model constants. The model also shows the natural emergence of two length scales, one associated with the familiar eddy cascade isotropic contribution to total dissipation and the other associated with the length scales of the canopy elements.     

A Case Study of Turbulence in the Nocturnal Boundary Layer During the Indian Summer Monsoon

Observations from the Cloud-Aerosol Interaction and Precipitation Enhancement Experiment-Integrated Ground Observation Campaign (CAIPEEX-IGOC) provide a rare opportunity to investigate nocturnal atmospheric surface-layer processes and surface-layer turbulent characteristics associated with the low-level jet (LLJ). Here, an observational case study of the nocturnal boundary layer is presented during the peak monsoon season over Peninsular India using data collected over a single night representative of the synoptic conditions of the Indian summer monsoon. Datasets based on Doppler lidar and eddy-covariance are used for this purpose. The LLJ is found to generate nocturnal turbulence by introducing mechanical shear at higher levels within the boundary layer. Sporadic and intermittent turbulent events observed during this period are closely associated with large eddies at the scale of the height of the jet nose. Flux densities in the stable boundary layer are observed to become non-local under the influence of the LLJ. Different turbulence regimes are identified, along with transitions between turbulent periods and intermittency. Wavelet analysis is used to elucidate the presence of large-scale eddies and associated intermittency during nocturnal periods in the surface layer. Although the LLJ is a regional-scale phenomenon it has far reaching consequences with regard to surface-atmosphere exchange processes.     

The Roles of Low-level Jets in “21·7” Henan Extremely Persistent Heavy Rainfall Event

An extremely heavy rainfall event lasting from 17 to 22 July 2021 occurred in Henan Province of China, with accumulated precipitation of more than 1000 mm over a 6-day period that exceeded its mean annual precipitation. The present study examines the roles of persistent low-level jets(LLJs) in maintaining the precipitation using surface station observations and reanalysis datasets. The LLJs triggered strong ascending motions and carried moisture mainly from the outflow of Typhoon In-fa(2021). Th...     

Turbulence Characteristics of the Shear-Free Convective Boundary Layer Driven by Heterogeneous Surface Heating

Large-eddy simulations (LESs) are employed to investigate the turbulence characteristics in the shear-free convective boundary layer (CBL) driven by heterogeneous surface heating. The patterns of surface heating are arranged as a chessboard with two different surface heat fluxes in the neighbouring patches, and the heterogeneity scale Λ in four different cases is taken as 1.2, 2.5, 5.0 and 10.0 km, respectively. The results are compared with those for the homogeneous case. The impact of the heterogeneity scale on the domain-averaged CBL characteristics, such as the profiles of the potential temperature and the heat flux, is not significant. However, different turbulence characteristics are induced by different heterogeneous surface heating. The greatest turbulent kinetic energy (TKE) is produced in the case with the largest heterogeneity scale, whilst the TKE in the other heterogeneous cases is close to that for the homogeneous case. This result indicates that the TKE is not enhanced unless the scale of the heterogeneous surface heating is large enough. The potential temperature variance is enhanced more significantly by a larger surface heterogeneity scale. But this effect diminishes with increasing CBL height, which implies that the turbulent eddy structures are changed during the CBL development. Analyses show that there are two types of organized turbulent eddies: one relates to the thermal circulations induced by the heterogeneous surface heating, whilst the other identifies with the inherent turbulent eddies (large eddies) induced by the free convection. At the early stage of the CBL development, the dominant scale of the organized turbulent eddies is controlled by the scale of the surface heterogeneity. With time increasing, the original pattern breaks up, and the vertical velocity eventually displays horizontal structures similar to those for the homogeneous heating case. It is found that after this transition, the values of λ/z _i (λ is the dominant horizontal scale of the turbulent eddies, z _i is the boundary-layer height) ≈1.6, which is just the aspect ratio of large eddies in the CBL.     

Simplified expressions for adjusting higher-order turbulent statistics obtained from open path gas analyzers

When density fluctuations of scalars such as CO₂ are measured with open-path gas analyzers, the measured vertical turbulent flux must be adjusted to take into account fluctuations induced by ‘external effects’ such as temperature and water vapour. These adjustments are needed to separate the effects of surface fluxes responsible for ‘natural’ fluctuations in CO₂ concentration from these external effects. Analogous to vertical fluxes, simplified expressions for separating the ‘external effects’ from higher-order scalar density turbulence statistics are derived. The level of complexity in terms of input to these expressions are analogous to that of the Webb–Pearman–Leuning (WPL), and are shown to be consistent with the conservation of dry air. It is demonstrated that both higher-order turbulent moments such as the scalar variances, the mixed velocity-scalar covariances, and the two-scalar covariance require significant adjustments due to ‘external effects’. The impact of these adjustments on the turbulent CO₂ spectra, probability density function, and dimensionless similarity functions derived from flux-variance relationships are also discussed.     

Dispersion of a Point-Source Release of a Passive Scalar Through an Urban-Like Array for Different Wind Directions

The dispersion of a point-source release of a passive scalar in a regular array of cubical, urban-like, obstacles is investigated by means of direct numerical simulations. The simulations are conducted under conditions of neutral stability and fully rough turbulent flow, at a roughness Reynolds number of Re _τ  = 500. The Navier–Stokes and scalar equations are integrated assuming a constant rate release from a point source close to the ground within the array. We focus on short-range dispersion, when most of the material is still within the building canopy. Mean and fluctuating concentrations are computed for three different pressure gradient directions (0°, 30°, 45°). The results agree well with available experimental data measured in a water channel for a flow angle of 0°. Profiles of mean concentration and the three-dimensional structure of the dispersion pattern are compared for the different forcing angles. A number of processes affecting the plume structure are identified and discussed, including: (i) advection or channelling of scalar down ‘streets’, (ii) lateral dispersion by turbulent fluctuations and topological dispersion induced by dividing streamlines around buildings, (iii) skewing of the plume due to flow turning with height, (iv) detrainment by turbulent dispersion or mean recirculation, (v) entrainment and release of scalar in building wakes, giving rise to ‘secondary sources’, (vi) plume meandering due to unsteady turbulent fluctuations. Finally, results on relative concentration fluctuations are presented and compared with the literature for point source dispersion over flat terrain and urban arrays.     

Application of transilient turbulent theory to study interactions between the atmospheric boundary layer and forest canopies

The new Forest-Land-Atmosphere ModEl called FLAME is presented. The first-order, nonlocal turbulence closure called transilient turbulence theory (Stull, 1993) is applied to study the interactions between a forested land-surface and the atmospheric boundary layer (ABL). The transilient scheme is used for unequal vertical grid spacing and includes the effects of drag, wake turbulence, and interference to vertical mixing by plant elements. Radiation transfer within the vegetation and the equations for the energy balance at the leaf surface have been taken from Norman (1979). Among others, the model predicts profiles of air temperature, humidity and wind velocity within the ABL, sensible and latent heat fluxes from the soil and the vegetation, the stomata and aerodynamic resistances, as well as profiles of temperature and water content in the soil. Preliminary studies carried out for a cloud free day and idealized initial conditions are presented. The canopy height is 30 m within a vertical domain of 3 km. The model is able to capture some of the effects usually observed within and above forested areas, including the relative wind speed maximum in the trunk space and the counter gradient-fluxes in the lower part of the plant stand. Of special interest is the determination of the location and magnitude of the turbulent mixing between model layers, which permits one to identify the effects of large eddies transporting momentum and scalar quantities into the canopy. A comparison between model simulations and field measurements will be presented in a future paper.     

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

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

Characteristics of canopy turbulence over a deciduous forest on complex terrain

Canopy turbulence plays an important role in mass and energy exchanges at the canopy-atmosphere interface. Despite extensive studies on canopy turbulence over a flat terrain, less attention has been given to canopy turbulence in a complex terrain. The purpose of this study is to scrutinize characteristics of canopy turbulence in roughness sublayer over a hilly forest terrain. We investigated basic turbulence statistics, conditionally sampled statistics, and turbulence spectrum in terms of different atmospheric stabilities, wind direction and vertical structures of momentum fluxes. Similarly to canopy turbulence over a homogeneous terrain, turbulence statistics showed coherent structure. Both quadrant and spectrum analysis corroborated the role of intermittent and energetic eddies with length scale of the order of canopy height, regardless of wind direction except for shift of peak in vertical wind spectrum to relatively high frequency in the down-valley wind. However, the magnitude of the momentum correlation coefficient in a neutral condition was smaller than typical value over a flat terrain. Further scrutiny manifested that, in the up-valley flow, temperature skewness was larger and the contribution of ejection to both momentum and heat fluxes was larger compared to the downvalley flow, indicating that thermal instability and weaker wind shear in up-valley flow asymmetrically affect turbulent transport within the canopy.     

Dispersion of a Passive Scalar Within and Above an Urban Street Network

The transport of a passive scalar from a continuous point-source release in an urban street network is studied using direct numerical simulation (DNS). Dispersion through the network is characterized by evaluating horizontal fluxes of scalar within and above the urban canopy and vertical exchange fluxes through the canopy top. The relative magnitude and balance of these fluxes are used to distinguish three different regions relative to the source location: a near-field region, a transition region and a far-field region. The partitioning of each of these fluxes into mean and turbulent parts is computed. It is shown that within the canopy the horizontal turbulent flux in the street network is small, whereas above the canopy it comprises a significant fraction of the total flux. Vertical fluxes through the canopy top are predominantly turbulent. The mean and turbulent fluxes are respectively parametrized in terms of an advection velocity and a detrainment velocity and the parametrization incorporated into a simple box-network model. The model treats the coupled dispersion problem within and above the street network in a unified way and predictions of mean concentrations compare well with the DNS data. This demonstrates the usefulness of the box-network approach for process studies and interpretation of results from more detailed numerical simulations.     

Low-Level Jets in the Moscow Region in Summer and Winter Observed with a Sodar Network

We evaluate the statistical properties of low-level jets (LLJs) observed by means of a network of Doppler sodars in the Moscow region, Russia. Continuous long-term measurements of the echo-signal intensity and wind-velocity profiles were carried out in July 2005 and in 2008–2010 synchronously in the centre of Moscow and at a rural site. The summertime nocturnal LLJs have a very clear diurnal cycle and exhibit features predicted by the Blackadar mechanism. In contrast, the long-lasting wintertime jets do not have any clear diurnal variability. The urban environment strongly influences LLJs in both seasons: above the city LLJs are higher, weaker and observed more rarely than at the rural site. In very cold periods (air temperature below −8°C) no LLJs were observed over the city, instead convection emerged in the urban boundary layer. The results are based on observations made in July 2005, January and December 2009, and January 2011.     

Large-Eddy Simulation of Coherent Turbulence Structures Associated with Scalar Ramps Over Plant Canopies

Large-eddy simulations were performed of a neutrally-stratified turbulent flow within and above an ideal, horizontally- and vertically-homogeneous plant canopy. Three simulations were performed for shear-driven flows in small and large computational domains, and a pressure-driven flow in a small domain, to enable the nature of canopy turbulence unaffected by external conditions to be captured. The simulations reproduced quite realistic canopy turbulence characteristics, including typical ramp structures appearing in time traces of the scalar concentration near the canopy top. Then, the spatial structure of the organised turbulence that caused the scalar ramps was examined using conditional sampling of three-dimensional instantaneous fields, triggered by the occurrence of ramp structures. A wavelet transform was used for the detection of ramp structures in the time traces. The ensemble-averaged results illustrate that the scalar ramps are associated with the microfrontal structure in the scalar, the ejection-sweep structure in the streamwise and vertical velocities, a laterally divergent flow just around the ramp-detection point, and a positive, vertically-coherent pressure perturbation. These vertical structures were consistent with previous measurements made in fields or wind tunnels. However, the most striking feature is that the horizontal slice of the same structure revealed a streamwise-elongated region of high-speed streamwise velocity impacting on another elongated region of low-speed velocity. These elongated structures resemble the so-called streak structures that are commonly observed in near-wall shear layers. Since elongated structures of essentially similar spatial scales were observed in all of the runs, these streak structures appear to be inherent in near-canopy turbulence. Presumably, strong wind shear formed just above the canopy is involved in their formation. By synthesis of the ensemble-averaged and instantaneous results, the following processes were inferred for the development of scalar microfronts and their associated flow structures: (1) a distinct scalar microfront develops where a coherent downdraft associated with a high-speed streak penetrates into the region of a low-speed streak; (2) a stagnation in flow between two streaks of different velocities builds up a vertically-coherent high-pressure region there; (3) the pressure gradients around the high-pressure region work to reduce the longitudinal variations in streamwise velocity and to enhance the laterally-divergent flow and lifted updrafts downstream of the microfront; (4) as the coherent mother downdraft impinges on the canopy, canopy-scale eddies are formed near the canopy top in a similar manner as observed in conventional mixing-layer turbulence.