An official website of the United States government
Stably stratified roughness sublayer flows are ubiquitous yet remain difficult to represent in models and to interpret using field experiments. Here, continuous high‐frequency potential temperature profiles from the forest floor up to 6.5 times the canopy height observed with distributed temperature sensing (DTS) are used to link eddy topology to roughness sublayer stability correction functions and coupling between air layers within and above the canopy. The experiments are conducted at two forest stands classified as hydrodynamically sparse and dense. Near‐continuous profiles of eddy sizes (length scales) and effective mixing lengths for heat are derived from the observed profiles using a novel conditional sampling approach. The approach utilizes potential temperature isoline fluctuations from a statically stable background state. The transport of potential temperature by an observed eddy is assumed to be conserved (adiabatic movement) and we assume that irreversible heat exchange between the eddy and the surrounding background occurs along the (vertical) periphery of the eddy. This assumption is analogous to Prandtl's mixing‐length concept, where momentum is transported rapidly vertically and then equilibrated with the local mean velocity gradient. A distinct dependence of the derived length scales on background stratification, height above ground, and canopy characteristics emerges from the observed profiles. Implications of these findings for (1) the failure of Monin–Obukhov similarity in the roughness sublayer and (2) above‐canopy flow coupling to the forest floor are examined. The findings have practical applications in terms of analysing similar DTS data sets with the proposed approach, modelling roughness sublayer flows, and interpreting nocturnal eddy covariance measurements above tall forested canopies.
more »
« less
Effects of Vegetation and Topography on the Boundary Layer Structure above the Amazon Forest
Abstract Observational data from two field campaigns in the Amazon forest were used to study the vertical structure of turbulence above the forest. The analysis was performed using the reduced turbulent kinetic energy (TKE) budget and its associated two-dimensional phase space. Results revealed the existence of two regions within the roughness sublayer in which the TKE budget cannot be explained by the canonical flat-terrain TKE budgets in the canopy roughness sublayer or in the lower portion of the convective ABL. Data analysis also suggested that deviations from horizontal homogeneity have a large contribution to the TKE budget. Results from LES of a model canopy over idealized topography presented similar features, leading to the conclusion that flow distortions caused by topography are responsible for the observed features in the TKE budget. These results support the conclusion that the boundary layer above the Amazon forest is strongly impacted by the gentle topography underneath.
more »
« less
- Award ID(s):
- 1644375
- PAR ID:
- 10181529
- Date Published:
- Journal Name:
- Journal of the Atmospheric Sciences
- Volume:
- 77
- Issue:
- 8
- ISSN:
- 0022-4928
- Page Range / eLocation ID:
- 2941 to 2957
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Dimensional analysis suggests that the dissipation length scale ( $$\ell _{{\it\epsilon}}=u_{\star }^{3}/{\it\epsilon}$$ ) is the appropriate scale for the shear-production range of the second-order streamwise structure function in neutrally stratified turbulent shear flows near solid boundaries, including smooth- and rough-wall boundary layers and shear layers above canopies (e.g. crops, forests and cities). These flows have two major characteristics in common: (i) a single velocity scale, i.e. the friction velocity ( $$u_{\star }$$ ) and (ii) the presence of large eddies that scale with an external length scale much larger than the local integral length scale. No assumptions are made about the local integral scale, which is shown to be proportional to $$\ell _{{\it\epsilon}}$$ for the scaling analysis to be consistent with Kolmogorov’s result for the inertial subrange. Here $${\it\epsilon}$$ is the rate of dissipation of turbulent kinetic energy (TKE) that represents the rate of energy cascade in the inertial subrange. The scaling yields a log-law dependence of the second-order streamwise structure function on ( $$r/\ell _{{\it\epsilon}}$$ ), where $$r$$ is the streamwise spatial separation. This scaling law is confirmed by large-eddy simulation (LES) results in the roughness sublayer above a model canopy, where the imbalance between local production and dissipation of TKE is much greater than in the inertial layer of wall turbulence and the local integral scale is affected by two external length scales. Parameters estimated for the log-law dependence on ( $$r/\ell _{{\it\epsilon}}$$ ) are in reasonable agreement with those reported for the inertial layer of wall turbulence. This leads to two important conclusions. Firstly, the validity of the $$\ell _{{\it\epsilon}}$$ -scaling is extended to shear flows with a much greater imbalance between production and dissipation, indicating possible universality of the shear-production range in flows near solid boundaries. Secondly, from a modelling perspective, $$\ell _{{\it\epsilon}}$$ is the appropriate scale to characterize turbulence in shear flows with multiple externally imposed length scales.more » « less
-
Abstract The Noah‐MP land surface model (LSM) relies on the Monin‐Obukhov (M‐O) Similarity Theory (MOST) to calculate land‐atmosphere exchanges of water, energy, and momentum fluxes. However, MOST flux‐profile relationships neglect canopy‐induced turbulence in the roughness sublayer (RSL) and parameterize within‐canopy turbulence in an ad hoc manner. We implement a new physics scheme (M‐O‐RSL) into Noah‐MP that explicitly parameterizes turbulence in RSL. We compare Noah‐MP simulations employing the M‐O‐RSL scheme (M‐O‐RSL simulations) and the default M‐O scheme (M‐O simulations) against observations obtained from 647 Snow Telemetry (SNOTEL) stations and two AmeriFlux stations in the western United States. M‐O‐RSL simulations of snow water equivalent (SWE) outperform M‐O simulations over 64% and 69% of SNOTEL sites in terms of root‐mean‐square‐error (RMSE) and correlation, respectively. The largest improvements in skill for M‐O‐RSL occur over closed shrubland sites, and the largest degradations in skill occur over deciduous broadleaf forest sites. Differences between M‐O and M‐O‐RSL simulated snowpack are primarily attributable to differences in aerodynamic conductance for heat underneath the canopy top, which modulates sensible heat flux. Differences between M‐O and M‐O‐RSL within‐canopy and below‐canopy sensible heat fluxes affect the amount of heat transported into snowpack and hence change snowmelt when temperatures are close to or above the melting point. The surface energy budget analysis over two AmeriFlux stations shows that differences between M‐O and M‐O‐RSL simulations can be smaller than other model biases (e.g., surface albedo). We intend for the M‐O‐RSL physics scheme to improve performance and uncertainty estimates in weather and hydrological applications that rely on Noah‐MP.more » « less
-
The significance of air flow within dense canopies situated on hilly terrain is not in dispute given its relevance to a plethora of applications in meteorology, wind energy, air pollution, atmospheric chemistry and ecology. While the mathematical description of such flows is complex, progress has proceeded through an interplay between experiments, mathematical modelling, and more recently large‐eddy simulations (LESs). In this contribution, LES is used to investigate the topography‐induced changes in the flow field and how these changes propagate to scalar transport within the canopy. The LES runs are conducted for a neutral atmospheric boundary layer above a tall dense forested canopy situated on a train of two‐dimensional sinusoidal hills. The foliage distribution is specified using leaf area density measurements collected in an Amazon rain forest. A series of LES runs with increasing hill amplitude are conducted to disturb the flow from its flat‐terrain state. The LES runs successfully reproduce the recirculation region and the flow separation on the lee‐side of the hill within the canopy region in agreement with prior laboratory and LES studies. Simulation results show that air parcels released inside the canopy have two preferential pathways to escape the canopy region: a “local” pathway similar to that encountered in flat terrain and an “advective” pathway near the flow‐separation region. Further analysis shows that the preferential escape location over the flow‐separation region leads to a “chimney”‐like effect that becomes amplified for air parcel releases near the forest floor. The work here demonstrates that shear‐layer turbulence is the main mechanism exporting air parcels out the canopy for both pathways. However, compared to flat terrain, the mean updraught at the flow separation induced by topography significantly shortens the in‐canopy residence time for air parcels released in the lower canopy, thus enhancing the export fraction of reactive gases.more » « less
-
Abstract Integrated quadrant analysis is a novel technique to identify and to characterize the trajectory and strength of turbulent coherent structures in the atmospheric surface layer. By integrating the three-dimensional velocity field characterized by traditional quadrant analysis with respect to time, the trajectory history of individual coherent structures can be preserved with Eulerian turbulence measurements. We develop a method to identify the ejection phase of coherent structures based on turbulence kinetic energy (TKE). Identifying coherent structures within a time series using TKE performs better than identifying them with the streamwise and vertical velocity components because some coherent structures are dominated by the cross-stream velocity component as they pass the sensor. By combining this identification method with the integrated quadrant analysis, one can animate or plot the trajectory of individual coherent structures from high-frequency velocity measurements. This procedure links a coherent ejection with the subsequent sweep and quiescent period in time to visualize and quantify the strength and the duration of a coherent structure. We develop and verify the method of integrated quadrant analysis with data from two field studies: the Eclipse Boundary Layer Experiment (EBLE) in Corvallis, Oregon in August 2017 (grass field) and the Vertical Cherry Array Experiment (VACE) in Linden, California in November 2019 (cherry orchard). The combined TKE identification method and integrated quadrant analysis are promising additions to conditional sampling techniques and coherent structure characterization because the identify coherent structures and couple the sweep and ejection components in space. In an orchard (VACE), integrated quadrant analysis verifies each coherent structure is dominated by a sweep. Conversely, above the roughness sublayer (EBLE), each coherent structure is dominated by an ejection.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1175/JAS-D-20-0063.1
