scholarly journals Closure Schemes for Stably Stratified Atmospheric Flows without Turbulence Cutoff

2016 ◽  
Vol 73 (12) ◽  
pp. 4817-4832 ◽  
Author(s):  
Dan Li ◽  
Gabriel G. Katul ◽  
Sergej S. Zilitinkevich

Abstract Two recently proposed turbulence closure schemes are compared against the conventional Mellor–Yamada (MY) model for stably stratified atmospheric flows. The Energy- and Flux-Budget (EFB) approach solves the budgets of turbulent momentum and heat fluxes and turbulent kinetic and potential energies. The Cospectral Budget (CSB) approach is formulated in wavenumber space and integrated across all turbulent scales to obtain flow variables in physical space. Unlike the MY model, which is subject to a “critical gradient Richardson number,” both EFB and CSB models allow turbulence to exist at any gradient Richardson number and predict a saturation of flux Richardson number () at sufficiently large . The CSB approach further predicts the value of and reveals a unique expression linking the Rotta and von Kármán constants. Hence, all constants in the CSB model are nontunable and stability independent. All models agree that the dimensionless sensible heat flux decays with increasing . However, the decay rate and subsequent cutoff in the MY model appear abrupt. The MY model further exhibits an abrupt cutoff in the turbulent stress normalized by vertical velocity variance, while the CSB and EFB models display increasing trends. The EFB model produces a rapid increase in the ratio of turbulent potential energy and vertical velocity variance as is approached, suggesting a strong self-preservation mechanism. Vertical anisotropy in the turbulent kinetic energy is parameterized in different ways in MY and EFB, but this consideration is not required in CSB. Differences between EFB and CSB model predictions originate from how the vertical anisotropy is specified in the EFB model.

2021 ◽  
Author(s):  
Dino Zardi

<p>Sloping terrains of any inclination favour the development, under the daily cycle of day time surface heating and night time cooling, of thermally-driven organised flows, displaying peculiar boundary layer structures, and eventually triggering the development of atmospheric convection.</p><p>The ubiquitous occurrence over the Earth of variously tilted surfaces - from gently sloping plains to steep cliffs, or valley and basin sidewalls – makes the understanding of such flows of utmost importance in view of the appropriate forecasting of the associated boundary layer transport processes. Also, they display a highly conceptual relevance, as they represent a prototypal situations for many other thermally driven-flows over complex terrain.   </p><p>An appropriate surface-layer scaling for slope wind is derived extending the classical analysis for flat horizontal terrain situations to the cover inclines. In the former, momentum and heat fluxes at the surface are two independent quantities, and vertical profiles of velocity and temperature can only be connected to them by means  of similiarity relationships, as fluxes are nearly invariant with height.</p><p>Instead, equations governing slope winds show that the mean wind and temperature profiles are closely connected to the flux structure normal to the slope, as this is not constant. Also, surface values of momentum flux and sensible heat flux are shown to be proportional to each other.</p><p>Based on the above relationships, suitable expressions are derived for the slope-normal profiles of velocity and temperature, both in the viscous sublayer and in the fully turbulent surface layer, as well as for the appropriate scaling factors in the two regions.</p>


2014 ◽  
Vol 760 ◽  
pp. 494-539 ◽  
Author(s):  
Stimit K. Shah ◽  
Elie Bou-Zeid

AbstractDirect numerical simulations of stably stratified Ekman layers are conducted to study the effect of increasing static stability on turbulence dynamics and modelling in wall-bounded flows at three moderate Reynolds numbers. The flow field is analysed by examining the mean profiles of wind speed, potential temperature and momentum flux, as well as streamwise velocity and temperature spectra. The maximum stabilizing buoyancy flux that a flow can sustain while remaining fully turbulent is found to depend on the Reynolds number. The flows with the highest Reynolds number display a relatively well-developed inertial range and logarithmic layer, and are found to bear similarities to much higher-Reynolds-number flows like the ones encountered in the atmospheric boundary layer. In particular, the near-wall mean profiles follow the Monin–Obukhov similarity theory. However, several flow features, such as the critical Richardson number and the stress–strain alignment, are found to maintain significant dependence on the Reynolds number. The budgets of turbulence kinetic energy (TKE), vertical velocity variance, momentum and buoyancy fluxes, and temperature variance are analysed. The results indicate that the effect of stability on turbulence is first directly manifested in the vertical velocity variance budget, and results in damping of vertical motions. This then leads to a reduction in the downward transport of horizontal momentum components towards the surface, and consequently to a decrease in the shear production term in the TKE budget: changes in the vertical profile of TKE shear production with increasing Richardson number are significant and have a larger impact on TKE than direct buoyancy destruction. The reduction in vertical velocity variance also results in significant drops in the production terms in the momentum flux, buoyancy flux and temperature variance budgets. Various assumptions and parameters related to low-order turbulence closures are investigated. The results suggest that the vertical velocity variance is a more appropriate parameter than the full TKE on which to base eddy-diffusivity and viscosity models.


1999 ◽  
Author(s):  
Namhyo Kim ◽  
David L. Rhode

Abstract A streamline curvature law of the wall is analytically derived to include very strong curved-channel wall curvature effects through a novel perturbation analysis. The new law allows improved analysis of such flows, and it provides the basis for improved wall function boundary conditions for their computation (CFD) over a wider range of y+, even for very strong curvature cases. The unique derivation is based on the Boussinesq eddy viscosity and curvature-corrected mixing length concepts, which is a linear function of the gradient Richardson number. For the first time, to include more complete curved flow physics, local streamline curvature effects in the gradient Richardson number are kept. To overcome the mathematical difficulty of keeping all of these local streamline curvature terms, a novel perturbation solution approach is successfully developed. This novel perturbation technique allows a closed-form analytical solution to many similar non-linear problems which previously required more complicated techniques. Qualitative and quantitative comparisons with measurements and previous curvature laws of the wall obtained by different approaches reveal that the new law shows improved representation of the wall curvature effects for all of the four test cases.


2021 ◽  
Vol 22 (10) ◽  
pp. 2547-2564
Author(s):  
Georg Lackner ◽  
Daniel F. Nadeau ◽  
Florent Domine ◽  
Annie-Claude Parent ◽  
Gonzalo Leonardini ◽  
...  

AbstractRising temperatures in the southern Arctic region are leading to shrub expansion and permafrost degradation. The objective of this study is to analyze the surface energy budget (SEB) of a subarctic shrub tundra site that is subject to these changes, on the east coast of Hudson Bay in eastern Canada. We focus on the turbulent heat fluxes, as they have been poorly quantified in this region. This study is based on data collected by a flux tower using the eddy covariance approach and focused on snow-free periods. Furthermore, we compare our results with those from six Fluxnet sites in the Arctic region and analyze the performance of two land surface models, SVS and ISBA, in simulating soil moisture and turbulent heat fluxes. We found that 23% of the net radiation was converted into latent heat flux at our site, 35% was used for sensible heat flux, and about 15% for ground heat flux. These results were surprising considering our site was by far the wettest site among those studied, and most of the net radiation at the other Arctic sites was consumed by the latent heat flux. We attribute this behavior to the high hydraulic conductivity of the soil (littoral and intertidal sediments), typical of what is found in the coastal regions of the eastern Canadian Arctic. Land surface models overestimated the surface water content of those soils but were able to accurately simulate the turbulent heat flux, particularly the sensible heat flux and, to a lesser extent, the latent heat flux.


2004 ◽  
Vol 43 (12) ◽  
pp. 1887-1899 ◽  
Author(s):  
Daniel J. Leathers ◽  
Daniel Graybeal ◽  
Thomas Mote ◽  
Andrew Grundstein ◽  
David Robinson

Abstract A one-dimensional snowpack model, a unique airmass identification scheme, and surface weather observations are used to investigate large ablation events in the central Appalachian Mountains of North America. Data from cooperative observing stations are used to identify large ablation events within a 1° latitude × 1° longitude grid box that covers the majority of the Lycoming Creek basin in northern Pennsylvania. All 1-day ablation events greater than or equal to 7.6 cm (3 in.) are identified for the period of 1950 through 2001. Seventy-one events are identified, and these days are matched with a daily airmass type derived using the Spatial Synoptic Classification technique. Average meteorological characteristics on ablation days of each airmass type are calculated in an effort to understand the diverse meteorological influences that led to the large ablation events. A one-dimensional mass and energy balance snowpack model (“SNTHERM”) is used to calculate surface/atmosphere energy fluxes responsible for ablation under each airmass type. Results indicate that large ablation events take place under diverse airmass/synoptic conditions in the central Appalachians. Five airmass types account for the 71 large ablation events over the 52-yr period. Forty-three of the events occurred under “moist” airmass types and 28 under “dry” airmass conditions. Large ablation events under dry airmass types are driven primarily by daytime net radiation receipt, especially net solar radiation. These events tend to occur early and late in the snow cover season when solar radiation receipt is highest and are characterized by relatively clear skies, warm daytime temperatures, and low dewpoint temperatures. Moist airmass types are characterized by cloudy, windy conditions with higher dewpoint temperatures and often with liquid precipitation. During these events sensible heat flux is most often the dominant energy flux to the snowpack during ablation episodes. However, in many cases there is also a significant input of energy to the snowpack associated with condensation. Combinations of high sensible and latent heat fluxes often result in extreme ablation episodes, similar to those witnessed in this area in January 1996.


2009 ◽  
Vol 13 (7) ◽  
pp. 987-998 ◽  
Author(s):  
Z. Gao ◽  
D. H. Lenschow ◽  
Z. He ◽  
M. Zhou

Abstract. In order to examine energy partitioning and CO2 exchange over a steppe prairie in Inner Mongolia, China, fluxes of moisture, heat and CO2 in the surface layer from June 2007 through June 2008 were calculated using the eddy covariance method. The study site was homogenous and approximately 1500 m×1500 m in size. Seasonal and diurnal variations in radiation components, energy components and CO2 fluxes are examined. Results show that all four radiation components changed seasonally, resulting in a seasonal variation in net radiation. The radiation components also changed diurnally. Winter surface albedo was higher than summer surface albedo because during winter the snow-covered surface increased the surface albedo. The seasonal variations in both sensible heat and CO2 fluxes were stronger than those of latent heat and soil heat fluxes. Sensible heat flux was the main consumer of available energy for the entire experimental period. The energy imbalance problem was encountered and the causes are analyzed.


2014 ◽  
Vol 11 (4) ◽  
pp. 5969-5995
Author(s):  
C. C. van Heerwaarden ◽  
A. J. Teuling

Abstract. This study investigates the difference in land–atmosphere interactions between grassland and forest during typical heat wave conditions in order to understand the controversial results of Teuling et al. (2010) (T10, hereafter), who have found the systematic occurrence of higher sensible heat fluxes over forest than over grassland during heat wave conditions. With a simple, but accurate coupled land–atmosphere model, we are able to reproduce the findings of T10 for both normal summer and heat wave conditions, and to carefully explore the sensitivity of the coupled land–atmosphere system to changes in incoming radiation and early-morning temperature. Our results emphasize the importance of fast processes during the onset of heat waves, since we are able to explain the results of T10 without having to take into account changes in soil moisture. In order to disentangle the contribution of differences in several static and dynamic properties between forest and grassland, we have performed an experiment in which new land use types are created that are equal to grassland, but with one of its properties replaced by that of forest. From these, we conclude that the closure of stomata in the presence of dry air is by far the most important process in creating the different behavior of grassland and forest during the onset of a heat wave. However, we conclude that for a full explanation of the results of T10 also the other properties (albedo, roughness and the ratio of minimum stomatal resistance to leaf-area index) play an important, but indirect role; their influences mainly consist of strengthening the feedback that leads to the closure of the stomata by providing more energy that can be converted into sensible heat. The model experiment also confirms that, in line with the larger sensible heat flux, higher atmospheric temperatures occur over forest.


2021 ◽  
pp. 079
Author(s):  
Ivan Bašták Ďurán ◽  
Pascal Marquet

Le schéma de turbulence Toucans est utilisé dans la configuration opérationnelle Alaro du modèle Aladin depuis début 2015. Son développement a été initié, guidé et en grande partie conçu par Jean-François Geleyn. Ce développement a commencé avec le prédécesseur du schéma Toucans, le schéma « pseudo-pronostique » en énergie cinétique turbulente, lui-même basé sur l'ancien schéma de turbulence de Louis, mais étendu dans Toucans à un schéma pronostique. Le schéma Toucans a pour objectif de traiter de manière cohérente les fonctions qui dépendent de la stabilité verticale de l'atmosphère, de l'influence de l'humidité et des échelles de longueur de la turbulence (de mélange et de dissipation). De plus, de nouvelles caractéristiques ont été ajoutées : une représentation améliorée pour les stratifications très stables (absence de nombre de Richardson critique), une meilleure représentation de l'anisotropie, un paramétrage unifié de la turbulence et des nuages par l'ajout d'une deuxième énergie turbulente pronostique et la paramétrisation des moments du troisième ordre. The Toucans turbulence scheme is a turbulence scheme that is used in the operational Alaro configuration of the Aladin model since early 2015. Its development was initiated, guided and to a large extend authored by Jean-François Geleyn. The development started with the predecessor of the Toucans scheme, the "pseudo-prognostic" turbulent kinetic energy scheme which itself was built on the "Louis" turbulence scheme, but extended to a prognostic scheme. The Toucans scheme aims for a consistent treatment of stability dependency functions, influence of moisture, and turbulence length scales. Additionally, new features were added to the turbulence scheme: improved representation of turbulence in very stable stratification (absence of critical gradient Richardson number), better representation of anisotropy, unified parameterization of turbulence and clouds via addition of second prognostic turbulence energy, and parameterization of third order moments.


2019 ◽  
Vol 23 (12) ◽  
pp. 5033-5058
Author(s):  
Guillaume Bigeard ◽  
Benoit Coudert ◽  
Jonas Chirouze ◽  
Salah Er-Raki ◽  
Gilles Boulet ◽  
...  

Abstract. The heterogeneity of Agroecosystems, in terms of hydric conditions, crop types and states, and meteorological forcing, is difficult to characterize precisely at the field scale over an agricultural landscape. This study aims to perform a sensitivity study with respect to the uncertain model inputs of two classical approaches used to map the evapotranspiration of agroecosystems: (1) a surface energy balance (SEB) model, the Two-Source Energy Balance (TSEB) model, forced with thermal infrared (TIR) data as a proxy for the crop hydric conditions, and (2) a soil–vegetation–atmosphere transfer (SVAT) model, the SEtHyS model, where hydric conditions are computed from a soil water budget. To this end, the models' skill was compared using a large and unique in situ database covering different crops and climate conditions, which was acquired over three experimental sites in southern France and Morocco. On average, the models provide 30 min estimations of latent heat flux (LE) with a RMSE of around 55 W m−2 for TSEB and 47 W m−2 for SEtHyS, and estimations of sensible heat flux (H) with a RMSE of around 29 W m−2 for TSEB and 38 W m−2 for SEtHyS. A sensitivity analysis based on realistic errors aimed to estimate the potential decrease in performance induced by the spatialization process. For the SVAT model, the multi-objective calibration iterative procedure (MCIP) is used to determine and test different sets of parameters. TSEB is run with only one set of parameters and provides acceptable performance for all crop stages apart from the early growing season (LAI < 0.2 m2 m−2) and when hydric stress occurs. An in-depth study on the Priestley–Taylor key parameter highlights its marked diurnal cycle and the need to adjust its value to improve flux partitioning between the sensible and latent heat fluxes (1.5 and 1.25 for France and Morocco, respectively). Optimal values of 1.8–2 were highlighted under cloudy conditions, which is of particular interest due to the emergence of low-altitude drone acquisition. Under developed vegetation (LAI > 0.8 m2 m−2) and unstressed conditions, using sets of parameters that only differentiate crop types is a valuable trade-off for SEtHyS. This study provides some scientific elements regarding the joint use of both approaches and TIR imagery, via the development of new data assimilation and calibration strategies.


2020 ◽  
Vol 66 (258) ◽  
pp. 543-555 ◽  
Author(s):  
Lindsey Nicholson ◽  
Ivana Stiperski

AbstractWe present the first direct comparison of turbulence conditions measured simultaneously over exposed ice and a 0.08 m thick supraglacial debris cover on Suldenferner, a small glacier in the Italian Alps. Surface roughness, sensible heat fluxes (~20–50 W m−2), latent heat fluxes (~2–10 W m−2), topology and scale of turbulence are similar over both glacier surface types during katabatic and synoptically disturbed conditions. Exceptions are sunny days when buoyant convection becomes significant over debris-covered ice (sensible heat flux ~ −100 W m−2; latent heat flux ~ −30 W m−2) and prevailing katabatic conditions are rapidly broken down even over this thin debris cover. The similarity in turbulent properties implies that both surface types can be treated the same in terms of boundary layer similarity theory. The differences in turbulence between the two surface types on this glacier are dominated by the radiative and thermal contrasts, thus during sunny days debris cover alters both the local surface turbulent energy fluxes and the glacier component of valley circulation. These variations under different flow conditions should be accounted for when distributing temperature fields for modeling applications over partially debris-covered glaciers.


Sign in / Sign up

Export Citation Format

Share Document