skip to main content

Title: On the Relationship between the TKE Dissipation Rate and the Temperature Structure Function Parameter in the Convective Boundary Layer
Abstract Under stably stratified conditions, the dissipation rate ε of turbulence kinetic energy (TKE) is related to the structure function parameter for temperature , through the buoyancy frequency and the so-called mixing efficiency. A similar relationship does not exist for convective turbulence. In this paper, we propose an analytical expression relating ε and in the convective boundary layer (CBL), by taking into account the effects of nonlocal heat transport under convective conditions using the Deardorff countergradient model. Measurements using unmanned aerial vehicles (UAVs) equipped with high-frequency response sensors to measure velocity and temperature fluctuations obtained during the two field campaigns conducted at Shigaraki MU observatory in June 2016 and 2017 are used to test this relationship between ε and in the CBL. The selection of CBL cases for analysis was aided by auxiliary measurements from additional sensors (mainly radars), and these are described. Comparison with earlier results in the literature suggests that the proposed relationship works, if the countergradient term γ D in the Deardorff model, which is proportional to the ratio of the variances of potential temperature θ and vertical velocity w , is evaluated from in situ (airplane and UAV) observational data, but fails if evaluated from large-eddy simulation (LES) results. This appears to be caused by the tendency of the variance of θ in the upper part of the CBL and at the bottom of the entrainment zone to be underestimated by LES relative to in situ measurements from UAVs and aircraft. We discuss this anomaly and explore reasons for it.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Journal of the Atmospheric Sciences
Page Range / eLocation ID:
2311 to 2326
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    This study evaluates the methods of identifying the heightziof the top of the convective boundary layer (CBL) during winter (December and January) over the Great Lakes and nearby land areas using observations taken by the University of Wyoming King Air research aircraft during the Lake-Induced Convection Experiment (1997/98) and Ontario Winter Lake-effect Systems (2013/14) field campaigns. Since CBLs facilitate vertical mixing near the surface, the most direct measurement ofziis that above which the vertical velocity turbulent fluctuations are weak or absent. Thus, we usezifrom the turbulence method as the “reference value” to whichzifrom other methods, based on bulk Richardson number (Rib), liquid water content, and vertical gradients of potential temperature, relative humidity, and water vapor mixing ratio, are compared. The potential temperature gradient method using a threshold value of 0.015 K m−1for soundings over land and 0.011 K m−1for soundings over lake provided the estimates ofzithat are most consistent with the turbulence method. The Ribthreshold-based method, commonly used in numerical simulation studies, underestimatedzi. Analyzing the methods’ performance on the averaging windowzavgwe recommend usingzavg= 20 or 50 m forziestimations for lake-effect boundary layers. The present dataset consists of both cloudy and cloud-free boundary layers, some having decoupled boundary layers above the inversion top. Because cases of decoupled boundary layers appear to be formed by nearby synoptic storms, we recommend use of the more general term, elevated mixed layers.

    Significance Statement

    The depthziof the convective atmospheric boundary layer (CBL) strongly influences precipitation rates during lake-effect snowstorms (LES). However, variousziapproximation methods produce significantly different results. This study utilizes extensive concurrently collected observations by project aircraft during two LES field studies [Lake-Induced Convection Experiment (Lake-ICE) and OWLeS] to assess howzifrom common estimation methods compare with “reference”ziderived from turbulent fluctuations, a direct measure of CBL mixing. For soundings taken both over land and lake; with cloudy or cloud-free conditions, potential temperature gradient (PTG) methods provided the best agreement with the referencezi. A method commonly employed in numerical simulations performed relatively poorly. Interestingly, the PTG method worked equally well for “coupled” and elevated decoupled CBLs, commonly associated with nearby cyclones.

    more » « less
  2. null (Ed.)
    Turbulence parameters in the lower troposphere (up to ~4.5 km) are estimated from measurements of high-resolution and fast-response cold-wire temperature and Pitot tube velocity from sensors onboard DataHawk Unmanned Aerial Vehicles (UAVs) operated at the Shigaraki Middle and Upper atmosphere (MU) Observatory during two ShUREX (Shigaraki UAV Radar Experiment) campaigns in 2016 and 2017. The practical processing methods used for estimating turbulence kinetic energy dissipation rate ε and temperature structure function parameter C T 2 from one-dimensional wind and temperature frequency spectra are first described in detail. Both are based on the identification of inertial (−5/3) subranges in respective spectra. Using a formulation relating ε and C T 2 valid for Kolmogorov turbulence in steady state, the flux Richardson number R f and the mixing efficiency χ m are then estimated. The statistical analysis confirms the variability of R f and χ m around ~ 0.13 − 0.14 and ~ 0.16 − 0.17 , respectively, values close to the canonical values found from some earlier experimental and theoretical studies of both the atmosphere and the oceans. The relevance of the interpretation of the inertial subranges in terms of Kolmogorov turbulence is confirmed by assessing the consistency of additional parameters, the Ozmidov length scale L O , the buoyancy Reynolds number R e b , and the gradient Richardson number Ri. Finally, a case study is presented showing altitude differences between the peaks of N 2 , C T 2 and ε , suggesting turbulent stirring at the margin of a stable temperature gradient sheet. The possible contribution of this sheet and layer structure on clear air radar backscattering mechanisms is examined. 
    more » « less
  3. Abstract

    The hydrostatic equilibrium addresses the approximate balance between the positive force of the vertical pressure gradient and the negative gravity force and has been widely assumed for atmospheric applications. The hydrostatic imbalance of the mean atmospheric state for the acceleration of vertical motions in the vertical momentum balance is investigated using tower, the global positioning system radiosonde, and Doppler lidar and radar observations throughout the diurnally varying atmospheric boundary layer (ABL) under clear-sky conditions. Because of the negligibly small mean vertical velocity, the acceleration of vertical motions is dominated by vertical variations of vertical turbulent velocity variances. The imbalance is found to be mainly due to the vertical turbulent transport of changing air density as a result of thermal expansion/contraction in response to air temperature changes following surface temperature changes. In contrast, any pressure change associated with air temperature changes is small, and the positive vertical pressure-gradient force is strongly influenced by its background value. The vertical variation of the turbulent velocity variance from its vertical increase in the lower convective boundary layer (CBL) to its vertical decrease in the upper CBL is observed to be associated with the sign change of the imbalance from positive to negative due to the vertical decrease of the positive vertical pressure-gradient force and the relative increase of the negative gravity force as a result of the decreasing upward transport of the low-density air. The imbalance is reduced significantly at night but does not steadily approach zero. Understanding the development of hydrostatic imbalance has important implications for understanding large-scale atmosphere, especially for cloud development.

    Significance Statement

    It is well known that the hydrostatic imbalance between the positive pressure-gradient force due to the vertical decrease of atmospheric pressure and the negative gravity forces in the vertical momentum balance equation has important impacts on the vertical acceleration of atmospheric vertical motions. Vertical motions for mass, momentum, and energy transfers contribute significantly to changing atmospheric dynamics and thermodynamics. This study investigates the often-assumed hydrostatic equilibrium and investigate how the hydrostatic imbalance is developed using field observations in the atmospheric boundary layer under clear-sky conditions. The results reveal that hydrostatic imbalance can develop from the large-eddy turbulent transfer of changing air density in response to the surface diabatic heating/cooling. The overwhelming turbulence in response to large-scale thermal forcing and mechanical work of the vast Earth surface contributes to the hydrostatic imbalance on large spatial and temporal scales in numerical weather forecast and climate models.

    more » « less
  4. Abstract

    Since the 1950s, a countergradient flux term has been added to some K-profile-based first-order PBL schemes, allowing them to simulate the slightly statically stable upper part of the convective boundary layer (CBL) observed in a limited number of aircraft soundings. There is, however, substantial uncertainty in inferring detailed CBL structure, particularly the level of neutral stability (zn), from such a limited number of soundings. In this study, composite profiles of potential temperature are derived from multiyear early afternoon radiosonde data over Beijing, China. The CBLs become slightly stable above zn ~ 0.31–0.33zi, where zi is the CBL depth. These composite profiles are used to evaluate two K-profile PBL schemes, the Yonsei University (YSU) and Shin–Hong (SH) schemes, and to optimize the latter through parameter calibration. In one-dimensional simulations using the WRF Model, YSU simulates a stable CBL above zn ~ 0.24zi, while default SH simulates a thick superadiabatic lower CBL with zn ~ 0.45zi. Experiments with the analytic solution of a K-profile PBL model show that adjusting the countergradient flux profile leads to significant changes in the thermal structure of CBL, informing the calibration of SH. The SH scheme replaces the countergradient heat flux term in its predecessor YSU scheme with a three-layer nonlocal heating profile, with fnl specifying the peak value and z*SL specifying the height of this peak value. Increasing fnl to 1.1 lowers zn, but to too low a value, while simultaneously increasing z*SL to 0.4 leads to a more appropriate zn ~ 0.36zi. The calibrated SH scheme performs better than YSU and default SH for real CBLs.

    more » « less
  5. Abstract Unsteadiness and horizontal heterogeneities frequently characterize atmospheric motions, especially within convective storms, which are frequently studied using large-eddy simulations (LES). The models of near-surface turbulence employed by atmospheric LES, however, predominantly assume statistically steady and horizontally homogeneous conditions (known as the equilibrium approach). The primary objective of this work is to investigate the potential consequences of such unrealistic assumptions in simulations of tornadoes. Cloud Model 1 (CM1) LES runs are performed using three approaches to model near-surface turbulence: the “semi-slip” boundary condition (which is the most commonly used equilibrium approach), a recently proposed nonequilibrium approach that accounts for some of the effects of turbulence memory, and a nonequilibrium approach based on thin boundary layer equations (TBLE) originally proposed by the engineering community for smooth-wall boundary layer applications. To be adopted for atmospheric applications, the TBLE approach is modified to account for the surface roughness. The implementation of TBLE into CM1 is evaluated using LES results of an idealized, neutral atmospheric boundary layer. LES runs are then performed for an idealized tornado characterized by rapid evolution, strongly curved air parcel trajectories, and substantial horizontal heterogeneities. The semi-slip boundary condition, by design, always yields a surface shear stress opposite the horizontal wind at the lowest LES grid level. The nonequilibrium approaches of modeling near-surface turbulence allow for a range of surface-shear-stress directions and enhance the resolved turbulence and wind gusts. The TBLE approach even occasionally permits kinetic energy backscatter from unresolved to resolved scales. Significance Statement The traditional approach of modeling the near-surface turbulence is not suitable for a tornado characterized by rapid evolution, strongly curved air parcel trajectories, and substantial horizontal heterogeneities. To understand the influence of statistically unsteady and horizontally heterogeneous near-surface conditions on tornadoes, this work adopts a fairly sophisticated approach from the engineering community and implements it into a widely used atmospheric model with necessary modifications. Compared to the traditional approach, the newly implemented approach produces more turbulent near-surface winds, more flexible surface-drag directions, and stronger wind gusts. These findings suggest a simulated tornado is very sensitive to the modeling approach of near-surface turbulence. 
    more » « less