Search NASA⌕ Search

SEARCH · Search NASA

Results for “Roughness sublayer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Probing eddy size and its effective mixing length in stably stratified roughness sublayer flows

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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The Turbulent Pressure Spectrum Within the Roughness Sublayer of a Subarctic Forest Canopy

The turbulent static pressure spectrum E pp (k x ) as a function of longitudinal wavenumber k x in the roughness sublayer of forested canopies is of interest to a plethora of problems such as pressure transport in the turbulent kinetic energy budget, pressure pumping from snow or forest floor, and coupling between flow within and above canopies. Long term static pressure measurements above a sub-arctic forested canopy for near-neutral conditions during the winter and spring were collected and analyzed for three snow cover conditions: trees and ground covered with snow, trees are snow free but the ground is covered with snow, and snow free cover. In all three cases, it is shown that E pp (k x ) obeys the attached eddy hypothesis at low wavenumbers (1/δ < k x < 1/z) —with E pp (k x ) ∝ u$^4_*$k$^{-1}_x$ and Kolmogorov scaling in the inertial subrange at higher wavenumbers—with E pp (k x ) ∝ ϵ 4/3 k x -7/3 , where u * is the friction velocity at the canopy top, is the mean turbulent kinetic energy dissipation rate, z is the distance from the snow top, and is the boundary layer depth. The implications of these two scaling laws to the normalized root-mean squared pressure C p = o p /u$^2_*$ and its newly proposed logarithmic scaling with normalized wall-normal distance z/δ are discussed for snow covered and snow free vegetation conditions. The work here also shows that k$^{-1}_x$ the in the E pp (k x ) appears more extensive and robust than its longitudinal velocity counterpart.

Aslan, Toprak [Finnish Meteorological Inst., Helsi↗
DOE OSTI2023

Turbulence regimes in the nocturnal roughness sublayer: Interaction with deep convection and tree mortality in the Amazon

We investigated the influence of seasonality and proximity to the forest canopy on nocturnal turbulence regimes in the roughness sublayer of a Central Amazon forest. Since convective systems of different scales are common in this region, we also analyzed the effect of extreme wind gusts (propagated from convective downdrafts) on the organization of the turbulence regimes, and their potential to cause the mortality of canopy trees. Our data include high-frequency winds, temperature and ozone concentration at different heights during the dry and wet seasons of 2014. In addition, we used critical wind-speed data derived from a tree-winching experiment and a modeling study conducted in the same study site. Two different turbulence regimes were identified at three heights above the canopy: a weakly stable (WS) and a very stable regime (VS). The threshold wind speeds that mark the transition between turbulence regimes were larger during the dry season and increased as a function of the height above the canopy. The turbulent fluxes of sensible heat and momentum during the WS accounted for 88% of the entire nighttime flux. Downdrafts occurred only in the WS and favored a fully coupled state of wind flow along the canopy profile. The destructive potential of winds was four times higher than on nights without downdrafts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Implementation and Evaluation of a Unified Turbulence Parameterization Throughout the Canopy and Roughness Sublayer in Noah‐MP Snow Simulations

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.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Merging flux-variance with surface renewal methods in the roughness sublayer and the atmospheric surface layer

In this article, Two micrometeorological methods utilizing high-frequency sampled air temperature were tested against eddy covariance (EC) sensible heat flux (𝐻) measurements at three sites representing agricultural, agro-forestry, and forestry systems. The two methods cover conventional and newly proposed forms of the flux-variance (FV) and surface renewal (SR) schemes of differing complexities. The sites represent measurements in surface, roughness, and roughness to surface transitional layers. Regression analyzes against EC show that the most reliable FV and SR forms estimate 𝐻 with slopes within ±10% from unity and coefficient of determination 𝑅2 > 0.9 across all the three sites. The best performance of both FV and SR was found at the agricultural site with measurements well within the surface layer, while the worst was found for the tall forest with measurements within the roughness sublayer where its thickness needed to be additionally estimated. The main variable driving 𝐻 in FV is the temperature variance, whereas in SR, it is the geometry of ramp-like structures. Since these structures are also responsible for most of the temperature variance, a novel FV-SR approach emerging from combining the methods is proposed and evaluated against EC measurements and conventional FV and SR schemes. The proposed FV-SR approach requiring only a single fast response thermocouple is potentially independent of calibration and ameliorates some of the theoretical objections that arise when combining ramp statistics with similarity arguments. The combination of methods also provides new insights into the contribution of coherent structures to the temperature variance and its dependence on atmospheric stratification. Other potential utility of the new method is to include it in multi-tool assessments of surface energy fluxes, since a convergence or divergence of the results has a high diagnostic value.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2006

Aeolian Sediment Transport through Large Patches of Roughness in the Atmospheric Inertial Sublayer

Roughness influences the flux of wind driven sand transport. In this paper, we report on sediment transport measurements for four different surface roughness configurations composed of the same size (solid) roughness elements in the atmospheric inertial sublayer (ISL). Results of these tests indicate that sediment transport rates through patches of roughness in the atmospheric inertial sublayer are to a large extent controlled and scale proportionally with the roughness density (lambda = n b h/S, where n is number of elements of breadth b and height h in area S) of the surface. However, element size apparently increases the magnitude of the reduction beyond that attributable to lambda. A sediment transport model that incorporates the effect of shear stress partitioning appears to predict reasonably well the effect of roughness on sand transport in the cases where the roughness elements are less than or equal to 0.10 m in height. However, when the dimensions of the roughness itself are equivalent to or are greater than the range of saltation lengths (vertical and horizontal), additional interactions of the elements with the saltation cloud appear to reduce the transport efficiency.

Gillies, John A.↗
DOE OSTI2020

Effects of Vegetation and Topography on the Boundary Layer Structure above the Amazon Forest

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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Turbulent Kinetic Energy Budgets over Gentle Topography Covered by Forests

Large-eddy simulations of flow over a “horizontally” uniform model forest are used to investigate the effects of gentle topography on the turbulent kinetic energy (TKE) budget within the canopy roughness sublayer. Despite significant differences between simulations using idealized sinusoidal topography and real topography of the Amazon forest, results indicate that the effects of topography are located predominantly in the upper canopy and above, and are mostly caused by mean advection of TKE. The “horizontally” averaged TKE budget from idealized and real gentle topographies are almost identical to that for flat terrain, including a clear inertial layer above the roughness sublayer in which shear production is balanced by local dissipation. At topography crests, where observational towers are usually located, mean vertical advection of TKE can be as important as horizontal advection. Here, we propose the use on an approximate TKE balance equation to estimate mean advection from single tower measurements, and introduce a new advection index that can be used as a proxy to quantify the importance of the topography on the TKE budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Machine Learning for Improving Surface-Layer-Flux Estimates

Abstract Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Topography-Induced TKE Budget Behavior Over an Amazon Forest

Data from three different heights (35, 50 and 81 m) of one of the ATTO Project towers in the Amazon forest were used to calculate the TKE (turbulence kinetic energy) budget and some other statistics within the RSL (roughness sublayer). The statistical analyses were carried out for unstable and stable cases. The vertical transport and the vertical advection terms do not explain the imbalances found in the TKE budget, highlighting the impor- tance of horizontal terms in complex terrain. Indeed, patterns independent of the time of day were found for vertical transport, horizontal and vertical velocity skewness, the mean vertical component velocity, and momentum and vertical TKE flux. Here, the highest values of mean vertical velocity were observed from the direction where there is a valley near the tower, which may indicate a topographical effect. Monin-Obukhov Similarity Theory (MOST) predictions for the mean velocity gradient, strictly applicable only in the inertial sublayer (ISL), gave reasonable predictions at 35 and 50 m, but dimensionless dissipation rates displayed large scatter. At 81 m, there were clear deviations from MOST for the two functions, disclosing the effect of topography. In addition, a neutral LES (large-eddy simu- lation) study showed good agreement with tower data for the same general wind direction.

ATTO project↗
DOE OSTI2025

The role of thermal stratification on the co‐spectral properties of momentum transport above an Amazonian forest

The influence of thermal stratification on the turbulent kinetic energy balance has been widely studied; however, its influence on the turbulent stress remains less explored in the presence of tall vegetated canopies and less ideal micrometeorological conditions. Here, the impact of thermal stratification on turbulent momentum flux is considered in the roughness sublayer (RSL) and the atmospheric surface layer (ASL) using the Amazon Tall Tower Observatory (ATTO) in Brazil. A scalewise co‐spectral budget (CSB) model is developed using standard closure schemes for the pressure–velocity decorrelation. The CSB revealed that the co‐spectrum ${F}_{wu}\left({k}_x\right)$ between longitudinal (u') and vertical (w') velocity fluctuations is impacted by the energy spectrum of the vertical velocity ${E}_{ww}\left({k}_x\right)$ and the much less studied longitudinal heat‐flux co‐spectrum ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$, where ${\theta}_{\mathrm{v}}^{\prime }$ are temperature fluctuations and ${k}_x$ is the longitudinal wavenumber. Under stable, very stable, and dynamic–convective conditions, the scaling exponent ${F}_{wu}\left({k}_x\right)$ in for the inertial subrange (ISR) scales is dominated by ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$ instead of ${E}_{ww}\left({k}_x\right)$. A near ${k}_x^{-7/3}$scaling in ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$ robust to large variations in thermal stratification is found, whereas the Kolmogorov ISR scaling for ${E}_{ww}\left({k}_x\right)\sim {k}_x^{-5/3}$ is not found. The scale‐dependent decorrelation time between u' and w' is dominated by ${\epsilon}^{-1/3}{k}_x^{-2/3}$ in the ISR, but is nearly constant for eddies larger than the vertical velocity integral scale, regardless of stability. Implications of these findings for generalized stability correction functions that are based on the turbulent stress budget instead of the turbulent kinetic energy budget are discussed.

canopy turbulence↗
NASA NTRS2020

Combining MODIS LAI with ICESat-Based Canopy Heights Improves Spaceborne Estimates of Vegetation Roughness Length for Momentum

Most land-surface models require parameterization of vertical wind profiles within the atmospheric boundary layer. For vegetated surfaces, it is common to assume a logarithmic profile in the surface layer, which includes estimates of vegetation roughness length for momentum (z0) and zero-plane displacement height (d0). This study finds that remotely-sensed forest canopy heights improve estimates of aerodynamic roughness length for momentum using a previously-developed representation of the roughness sublayer (Raupach 1992; Jasinski et al. 2005). Resulting roughness products consist of two datasets: 1) 14 years of 8-day snapshots of the global land surface at a nominal spatial resolution of 500-meters for users who wish to retain full temporal resolution and interannual variability; and 2) multiyear averages of the 8-day snapshots, here referred to as "climatologies" of roughness, which retain underlying seasonality. Both products are suitable for use in data assimilation and reanalyses such as the National Climate Assessment Land Data Assimilation System (NCA-LDAS), for which these products were initially developed.

Borak, Jordan S.↗
NASA NTRS2005

Bulk Surface Momentum Parameters for Satellite-Derived Vegetation Fields

The bulk aerodynamic parameters associated with the absorption of surface momentum by vegetated landscapes are theoretically estimated within the context of Raupach's roughness sublayer formulation. The parameters include the bulk plant drag coefficient, maximum u*/U(sub h), sheltering coefficient, and canopy area density at onset of sheltering. Parameters are estimated for the four principal IGBP land cover classes within the U.S. Southern Great Plains: evergreen needleleaf forests, grasslands, croplands, and open shrublands. The estimation approach applies the Method of Moments to roughness data from several international field experiments and other published sources. The results provide the necessary land surface parameters for satellite-based estimation of momentum aerodynamic roughness length and zero-plane displacement height for seasonally variable vegetation fields employed in most terrestrial and atmospheric simulation models used today. Construction of sample displacement and roughness maps over the Southern United States using MODIS land products demonstrates the potential of this approach for regional to global applications.

Jasinski, Michael F.↗
DOE OSTI2024

The vertical-velocity skewness in the atmospheric boundary layer without buoyancy and Coriolis effects

One of the main features of near-neutral atmospheric boundary layer (ABL) turbulence is the positive vertical velocity skewness $Sk_w$ above the roughness sublayer or the buffer region in smooth-walls. The $Sk_w$ variations are receiving renewed interest in many climate-related parameterizations of the ABL given their significance to cloud formation and to testing sub-grid schemes for Large Eddy Simulations (LES). The vertical variations of $Sk_w$ are explored here using wind tunnel and flume experiments collected above smooth, rough, and permeable-walls in the absence of buoyancy and Coriolis effects. These laboratory experiments form a necessary starting point to probe the canonical structure of $Sk_w$ as they deal with a key limiting case (i.e., near-neutral conditions). Diagnostic models based on cumulant expansions, realizability constraints, and constant mass flux approach routinely employed in the convective boundary layer as well as prognostic models based on third-order budgets are used to explain variations in $Sk_w$ for the idealized laboratory conditions. The failure of flux-gradient relations to model $Sk_w$ from the gradients of the vertical velocity variance σ$_w^2$ are explained and corrections based on models of energy transport offered. Novel links between the diagnostic and prognostic models are also featured, especially for the inertial term in the third-order budget of the vertical velocity fluctuation. The co-spectral properties of w′/σ w vs w′ 2 /σ$_w^2$ are also presented for the first time to assess the dominant scales governing $Sk_w$ in the inner and outer layers, where w′ is the fluctuating vertical velocity and σ w is the vertical velocity standard deviation.>

Boundary layer flow↗
DOE OSTI2025

Scaling Exponents of Turbulent Static Pressure Structure Function in the Inertial Subrange

The measured variations in the turbulent static pressure structure function D pp (r) with scale $r$ in the roughness sublayer above a subarctic forest are empirically shown to exhibit exponents that are smaller than $r$ 4/3 predicted for the inertial subrange (ISR). Three hypotheses are offered to explain these deviations. The first is based on conventional intermittency correction to the averaged turbulent kinetic energy dissipation rate, the second is based on shearing introducing deviations from locally isotropic state that must be sensed by both velocity and pressure structure functions, and the third is based on large and inertial scale pressure interactions that persist at values of $r$ within the resolvable ISR. The third hypothesis is shown to yield superior results, which allows a new formulation for D pp (r) to be derived that accommodates such finite interactions.

structure function↗
NASA NTRS1987

A comparative study of some mathematical models of the mean wind structure and aerodynamic drag of plant canopies

A semianalytical method for describing the mean wind profile and shear stress within plant canopies and for estimating the roughness length and the displacement height is presented. This method incorporates density and vertical structure of the canopy and includes simple parameterizations of the roughness sublayer and shelter factor. Some of the wind profiles examined are consistent with first-order closure techniques while others are consistent with second-order closure techniques. Some profiles show a shearless region near the base of the canopy; however, none displays a secondary maximum there. Comparing several different analytical expressions for the canopy wind profile against observations suggests that one particular type of profile (an Airy function which is associated with the triangular foliage surface area density distribution) is superior to the others. Because of the numerical simplicity of the methods outlined, it is suggested that they may be profitably used in large-scale models of plant-atmosphere exchanges.

Massman, William↗
NASA NTRS1994

An evaluation of the regional acid deposition model surface module for ozone uptake at three sites in the San Joaquin Valley of California

Plants and soils act as major sinks for the destruction of tropospheric ozone, especially during daylight hours when plant stomata open and are thought to provide the dominant pathway for the uptake of ozone. The present study, part of the California Ozone Deposition Experiment, compares predictions of the regional acid deposition model ozone surface conductance module with surface conductance data derived from eddy covariance measurements of ozone flux taken at a grape, a cotton, and a grassland site in the San Joaquin Valley of California during the summer of 1991. Results indicate that the model (which was developed to provide long-term large-area estimates for the eastern United States) significantly overpredicts the surface conductance at all times of the day for at least two important types of plant cover of the San Joaquin Valley and that it incorrectly partitions the ozone flux between transpiring and nontranspiring components of the surface at the third site. Consequently, the model either overpredicts or inaccurately represents the observed deposition velocities. Other results indicate that the presence of dew does not reduce the rate of ozone deposition, contradicting to model assumptions, and that model assumptions involving the dependency of stomata upon environmental temperature are unnecessary. The effects of measurement errors and biases, arising from the presence of the roughness sublayer and possible photochemical reactions, are also discussed. A simpler model for ozone surface deposition (at least for the San Joaquin Valley) is proposed and evaluated.

Massman, W. J.↗