Comparisons of Planetary Boundary Layer Height from Ceilometer with ARM Radiosonde Data

Ceilometer measurements of aerosol backscatter profiles have been widely used to provide continuous PBLHT estimations. To investigate the robustness of ceilometer-estimated PBLHT under different atmospheric conditions, we compared ceilometer- and radiosonde-estimated PBLHTs using long term U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) ceilometer and balloon-borne sounding data at three ARM fixed-location atmospheric observatories and from three ARM mobile observatories deployed around the world for various field campaigns, which cover from Tropics to Polar regions and over both ocean and land surfaces. Statistical comparisons of ceilometer-estimated PBLHTs from the Vaisala CL31 ceilometer data with radiosonde-estimated PBLHTs from the ARM PBLHT-SONDE Value-added Product (VAP) are performed under different atmospheric conditions including stable and unstable atmospheric boundary layer, low-level cloud-free, and cloudy conditions at these ARM observatories. Under unstable atmospheric boundary layer conditions, good comparisons are found between ceilometer- and radiosonde-estimated PBLHTs at ARM low- and mid-latitude land observatories. However, it is still challenging to obtain reliable PBLHT estimations over ocean surfaces even using radiosonde data. Under stable atmospheric boundary layer conditions, ceilometer- and radiosonde-estimated PBLHTs have weak correlations. Among different PBLHT estimations utilizing the Heffter, the Liu-Liang, and the bulk Richardson number methods in the ARM PBLHT-SONDE VAP, ceilometer-estimated PBLHTs have better comparisons with the Liu-Liang method under unstable and with the bulk Richardson number method under stable atmospheric boundary layer conditions.

Emergence of optimal disturbances in a stratified turbulent shear flow under the stochastic forcing

Large-scale inclined organized structures in stably stratified turbulent shear flows were revealed in the numerical simulation and indirectly confirmed by the field measurements in the stable atmospheric boundary layer. Spatial scales and forms of these structures coincide with those of the optimal disturbances of a simplified linear model. In this paper, we clarify the relation between the organized structures and the optimal disturbances, analyzing a time series of turbulent fields obtained by the RANS model with eddy viscosity/diffusivity and stochastic forcing generating the small-scale turbulence.

The Influence of Slopes of Isolated Three-Dimensional Valleys on Near-Surface Turbulence

Motivated by a limited understanding of how valleys affect near-surface turbulence, characterizations of neutrally stable atmospheric-boundary-layer flows over isolated valleys are presented. In particular, the influence of the slopes of the three-dimensional ridges that form the idealized valleys are investigated. Flows over three distinct symmetric valley geometries were modelled in a large boundary-layer wind tunnel. For each valley geometry, the high-resolution measurements from the crests of each of the ridges and the midpoint between them are compared with an undisturbed moderately rough classed boundary-layer flow over flat terrain with homogeneous surface roughness. Flow separation originates above the crests of the first ridges of all geometries and generates recirculation zones. These are characterized by slope-dependent increases in three-dimensional near-surface turbulence when compared with the attached flows further upstream. The recirculation zones longitudinally extend to roughly half the valley width. Above the crests of the second ridges, the longitudinal velocity component decreases and turbulence intensity increases when compared with the flows above the crests of the first ridges. Results also exhibit significant increases of turbulence above the inner-valley regions of all geometries.

On the u☆-U Relationship in the Stable Atmospheric Boundary Layer over Arctic Sea Ice

A relationship between the friction velocity u☆ and mean wind speed U in a stable atmospheric boundary layer (ABL) over Arctic sea ice was considered. To that aim, the observations collected during the Surface Heat Budget of the Arctic Ocean (SHEBA) experiment were used. The observations showed the so-called “hockey-stick” shape of the u☆−U relationship, which consists of a slow increase of u☆ with increasing wind speed for U<Utr and a more rapid almost linear increase of u☆ for U>Utr, where Utr is the wind speed of transition between the two regimes. Such a relationship is most pronounced at the highest observational levels, namely at 9 and 14 m, and is also sharper when the air-surface temperature difference exceeds its average values for stable conditions. It is shown that the Monin–Obukhov similarity theory (MOST) reproduces the observed u☆−U relationship rather well. This suggests that at least for the SHEBA dataset, there is no contradiction between MOST and the “hockey-stick” shape of the u☆−U relationship. However, the SHEBA data, as well as the single-column simulations show that for cases with strong stability, u☆ significantly decreases with height due to the shallowness of the ABL. It was shown that when u☆ was assumed independent of height, the value of the normalized drag coefficient, i.e., of the so-called stability correction function for momentum, calculated using observations at a certain level, can be significantly underestimated. To overcome this, the decrease of u☆ with height was taken into account in the framework of MOST using local scaling instead of the scaling with surface fluxes. Using such an extended MOST brought the estimates of the normalized drag coefficient closer to the Businger–Dyer relation.

Spatiotemporal Dynamics of the Kinetic Energy in the Atmospheric Boundary Layer from Minisodar Measurements

Spatiotemporal dynamics of the atmospheric kinetic energy and its components caused by the ordered and turbulent motions of air masses are estimated from minisodar measurements of three velocity vector components and their variances within the lowest 5–200 m layer of the atmosphere, with a particular emphasis on the turbulent kinetic energy. The layered structure of the total atmospheric kinetic energy has been established. From the diurnal hourly dynamics of the altitude profiles of the turbulent kinetic energy (TKE) retrieved from minisodar data, four layers are established by the character of the altitude TKE dependence, namely, the near-ground layer, the surface layer, the layer with a linear TKE increase, and the transitive layer above. In the first layer, the most significant changes of the TKE were observed in the evening hours. In the second layer, no significant changes in the TKE values were observed. A linear increase in the TKE values with altitude was observed in the third layer. In the fourth layer, the TKE slightly increased with altitude and exhibited variations during the entire observation period. The altitudes of the upper boundaries of these layers depended on the time of day. The MKE values were much less than the corresponding TKE values, they did not exceed 50 m2/s2. From two to four MKE layers were distinguished based on the character of its altitude dependence. The two-layer structures were observed in the evening and at night (under conditions of the stable atmospheric boundary layer). In the morning and daytime, the four-layer MKE structures with intermediate layers of linear increase and subsequent decrease in the MKE values were observed. Our estimates demonstrated that the TKE contribution to the total atmospheric kinetic energy considerably (by a factor of 2.5–3) exceeded the corresponding MKE contribution.

Spatiotemporal Dynamics of the Kinetic Energy in the Atmospheric Boundary Layer from Minisodar Measurements

Spatiotemporal dynamics of the atmospheric kinetic energy and its components caused by the ordered and turbulent motions of air masses are estimated from minisodar measurements of three velocity vector components and their variances within the lowest 5&ndash;200 m layer of the atmosphere, with a particular emphasis on the turbulent kinetic energy. The layered structure of the total atmospheric kinetic energy has been established. From the diurnal hourly dynamics of the altitude profiles of the turbulent kinetic energy (TKE) retrieved from minisodar data, four layers are established by the character of the altitude TKE dependence, namely, the near-ground layer, the surface layer, the layer with a linear TKE increase, and the transitive layer above. In the first layer, the most significant changes of the ТКЕ were observed in the evening hours. In the second layer, no significant changes in the TKE values were observed. A linear increase in the TKE values with altitude was observed in the third layer. In the fourth layer, the TKE slightly increased with altitude and exhibited variations during the entire observation period. The altitudes of the upper boundaries of these layers depended on the time of day. The MKE values were much less than the corresponding TKE values, they did not exceed 50 m2/s2. From two to four MKE layers were distinguished based on the character of its altitude dependence. The two-layer structures were observed in the evening and at night (under conditions of the stable atmospheric boundary layer). In the morning and daytime, the four-layer MKE structures with intermediate layers of linear increase and subsequent decrease in the MKE values were observed. Our estimates demonstrated that the ТКЕ contribution to the total atmospheric kinetic energy considerably (by a factor of 2.5&ndash;3) exceeded the corresponding МКЕ contribution.

Explicit Algebraic Reynolds-stress Modelling of a Convective Atmospheric Boundary Layer Including Counter-Gradient Fluxes

In a recent study (Želi et al. in Bound Layer Meteorol 176:229–249, 2020), we have shown that the explicit algebraic Reynolds-stress (EARS) model, implemented in a single-column context, is able to capture the main features of a stable atmospheric boundary layer (ABL) for a range of stratification levels. We here extend the previous study and show that the same formulation and calibration of the EARS model also can be applied to a dry convective ABL. Five different simulations with moderate convective intensities are studied by prescribing surface heat flux and geostrophic forcing. The results of the EARS model are compared to large-eddy simulations of Salesky and Anderson (J Fluid Mech 856:135–168, 2018). It is shown that the EARS model performs well and is able to capture the counter-gradient heat flux in the upper part of the ABL due to the presence of the non-gradient term in the relation for vertical turbulent heat flux. The model predicts the full Reynolds-stress tensor and heat-flux vector and allows us to compare other important aspects of a convective ABL such as the profiles of vertical momentum variance. Together with the previous studies, we show that the EARS model is able to predict the essential features of the ABL. It also shows that the EARS model with the same model formulation and coefficients is applicable over a wide range of stable and moderately unstable stratifications.

Scale-Dependent Turbulent Dynamics and Phase-Space Behavior of the Stable Atmospheric Boundary Layer

The structure of turbulent dynamics in a stable atmospheric boundary layer was studied by means of a phase-space description. Data from the CASES-99 experiment, decomposed in local modes (with increasing time scale) using empirical mode decomposition, were analyzed in order to extract the proper time lag and the embedding dimension of the phase-space manifold, and subsequently to estimate their scale-dependent correlation dimension. Results show that the dynamics are low-dimensional and anisotropic for a large scale, where the flow is dominated by the bulk motion. Then, they become progressively more high-dimensional while transiting into the inertial sub-range. Finally, they reach three-dimensionality in the range of scales compatible with the center of the inertial sub-range, where the phase-space-filling turbulent fluctuations dominate the dynamics.