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 StatementIt 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
This content will become publicly available on March 27, 2026
Total energy conservation for explaining discrepancies between field observations and the traditional thermal energy balance
Abstract To address the inability of the traditional thermal energy balance to explain field observations in the atmospheric boundary layer (ABL), this study investigates total energy conservation for atmospheric applications. Total energy conservation serves as a unique constraint on variations among different energy forms, especially when these variations are interconnected, as observed in the non‐isothermal atmosphere. In contrast, the first law of thermodynamics for the derivation of the traditional thermal energy balance is a special case of total energy conservation for a closed system at rest with thermal energy variations only. By keeping both total and kinetic energy conservation, a generalized thermal energy balance is found to contain two additional terms related to airflows compared with the traditional balance. They represent the thermal energy contribution to non‐hydrostatic energy transfer via heat transfer by vertical airflows and dissipation heating via air viscosity due to airflow deformation. Because of the effective vertical heat transfer by airflows, the non‐hydrostatic energy transfer contributes to major stability‐dependent differences between the traditional and generalized thermal energy balances. The stability‐dependent bias of the traditional balance is consistent with disagreement between field observations and traditional theoretical expectations such as the well‐known observed surface thermal energy imbalance, among others.
more »
« less
- PAR ID:
- 10646713
- Publisher / Repository:
- Wiley
- Date Published:
- Journal Name:
- Quarterly Journal of the Royal Meteorological Society
- Volume:
- 151
- Issue:
- 771
- ISSN:
- 0035-9009
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract 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 between longitudinal () and vertical () velocity fluctuations is impacted by the energy spectrum of the vertical velocity and the much less studied longitudinal heat‐flux co‐spectrum , where are temperature fluctuations and is the longitudinal wavenumber. Under stable, very stable, and dynamic–convective conditions, the scaling exponent in for the inertial subrange (ISR) scales is dominated by instead of . A near scaling in robust to large variations in thermal stratification is found, whereas the Kolmogorov ISR scaling for is not found. The scale‐dependent decorrelation time between and is dominated by 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.more » « less
-
Abstract The water vapor transport associated with latent heat flux (LE) in the planetary boundary layer (PBL) is critical for the atmospheric hydrological cycle, radiation balance, and cloud formation. The spatiotemporal variability of LE and water vapor mixing ratio (rv) are poorly understood due to the scale‐dependent and nonlinear atmospheric transport responses to land surface heterogeneity. Here, airborne in situ measurements with the wavelet technique are utilized to investigate scale‐dependent relationships among LE, vertical velocity (w) variance (), andrvvariance () over a heterogeneous surface during the Chequamegon Heterogeneous Ecosystem Energy‐balance Study Enabled by a High‐density Extensive Array of Detectors 2019 (CHEESEHEAD19) field campaign. Our findings reveal distinct scale distributions of LE, , and at 100 m height, with a majority scale range of 120 m–4 km in LE, 32 m–2 km in , and 200 m–8 km in . The scales are classified into three scale ranges, the turbulent scale (8–200 m), large‐eddy scale (200 m–2 km), and mesoscale (2–8 km) to evaluate scale‐resolved LE contributed by and . The large‐eddy scale in PBL contributes over 70% of the monthly mean total LE with equal parts (50%) of contributions from and . The monthly temporal variations mainly come from the first two major contributing classified scales in LE, , and . These results confirm the dominant role of the large‐eddy scale in the PBL in the vertical moisture transport from the surface to the PBL, while the mesoscale is shown to contribute an additional ∼20%. This analysis complements published scale‐dependent LE variations, which lack detailed scale‐dependent vertical velocity and moisture information.more » « less
-
Abstract How large turbulent eddies influence non‐closure of the surface energy balance is an active research topic that cannot be uncovered by the mean continuity equation in isolation. It is demonstrated here that asymmetric turbulent flux transport of heat and water vapor by sweeps and ejections of large eddies under unstable atmospheric stability conditions reduce fluxes. Such asymmetry causes positive gradients in the third‐order moments in the turbulent flux budget equations, primarily attributed to substantially reduced flux contributions by sweeps and sustained large flux contributions by ejections. Small‐scale surface heterogeneity in heating generates ejecting eddies with larger air temperature variance than sweeping eddies, causing asymmetric flux transport in the atmospheric surface layer. Changes in asymmetry with increasing instability are congruent with observed increases in the surface energy balance non‐closure. To assess the contributions of asymmetric flux transport by large eddies to the non‐closure requires two eddy covariance systems on the tower to measure the gradients of the turbulent heat flux and other third‐order moments.more » « less
-
Abstract. Modelling the pressure in the Earth's interior is a common problem in Earth sciences. In this study we propose a method based on the conservation of the momentum of a fluid by using a hydrostatic scenario or a uniformly moving fluid to approximate the pressure. This results in a partial differential equation (PDE) that can be solved using classical numerical methods. In hydrostatic cases, the computed pressure is the lithostatic pressure. In non-hydrostatic cases, we show that this PDE-based approach better approximates the total pressure than the classical 1D depth-integrated approach. To illustrate the performance of this PDE-based formulation we present several hydrostatic and non-hydrostatic 2D models in which we compute the lithostatic pressure or an approximation of the total pressure, respectively. Moreover, we also present a 3D rift model that uses that approximated pressure as a time-dependent boundary condition to simulate far-field normal stresses. This model shows a high degree of non-cylindrical deformation, resulting from the stress boundary condition, that is accommodated by strike-slip shear zones. We compare the result of this numerical model with a traditional rift model employing free-slip boundary conditions to demonstrate the first-order implications of considering “open” boundary conditions in 3D thermo-mechanical rift models.more » « less
An official website of the United States government
