Convective Atmospheric Boundary Layer Research Articles

Modification of wind turbine wakes by large-scale, convective atmospheric boundary layer structures

In this study, we consider the impact of large-scale, convective structures in an unstable atmospheric boundary layer on wind turbine wakes. Simulation data from a high-fidelity large-eddy simulation (LES) of the AWAKEN wind farm site matching unstable atmospheric conditions were analyzed, and both turbine performance and wake behavior were affected based on their location relative to the convective structures. Turbines located in updraft regions of the flow experienced lower inflow velocity and generated less power, but their wakes were observed to recover faster and saw greater turbulent kinetic energy mixing higher in the boundary layer. The opposite effect was found for turbines in the downdraft regions of the convective structures. A simplified model of this wake behavior was also developed based on a two-dimensional k–ε Reynolds-Averaged Navier–Stokes formulation. This simplified model included the effects of vertical transport, but could be efficiently solved as a parabolic system, and was found to capture similar wake modifications observed in the high-fidelity LES computations.

Quantification and assessment of the atmospheric boundary layer height measured during the AWAKEN experiment by a scanning LiDAR

The atmospheric boundary layer (ABL) height plays a key role in many atmospheric processes as one of the dominant flow length scales. However, a systematic quantification of the ABL height over the entire range of scales (i.e., with periods ranging from one minute to one year) is still lacking in literature. In this work, the ABL height is quantified based on high-resolution measurements collected by a scanning pulsed Doppler LiDAR during the recent American WAKE experimeNt (AWAKEN) campaign. The high availability of ABL height estimates (≈2200 collected over one year and each of them based on 10-min averaged statistics) allows to robustly assess five different ABL height models, i.e., one for convective thermal conditions and four for stable conditions. Thermal condition is quantified by a stability parameter spanning three orders of magnitude and probed by near-ground 3D sonic anemometry. The free-atmosphere stability, quantified by the Brunt–Väisälä frequency, is both calculated from simultaneous radiosonde measurements and obtained from the best fit of two of the chosen ABL height models. Good agreement is found between the data and three of the chosen models, quantified by mean absolute errors on the ABL height between 281 and 585 m. Furthermore, the seasonal variability of the convective ABL height model parameters (−15% to +23% with respect to the year baseline) agrees with the variability of buoyancy-generated turbulence caused by the variation in solar radiation throughout the year.

A Variant of the Local Similarity Theory and Approximations of Vertical Profiles of Turbulent Moments of the Atmospheric Convective Boundary Layer

The approximation of the turbulent moments of the atmospheric convective layer is based on a variant of the local similarity theory using the concepts of the semi-empirical theory of Prandtl turbulence. In the proposed variant of the local similarity theory, the second moment of vertical velocity and the “spectral” Prandtl mixing length are used as basic parameters. This approach allows us to extend Prandtl’s theory to turbulent moments of vertical velocity and buoyancy and additionally offer more than ten new approximations. The comparison of the proposed approximation with other variants of the theory of local similarity is considered. It is shown that the selected basic parameters significantly improve the agreement between the local similarity approximations and experimental data. The approximations are consistent with observations in the turbulent convective layer of the atmosphere, the upper boundary of which nearly corresponds to the lower boundary of the temperature inversion. Analytical approximations of local similarity can find applications in the construction of high-order moment closures in the vortex of resolving numerical turbulence models, as well as in the construction of “mass-flux” parametrization.

Large-scale semi-organized rolls in a sheared convective turbulence: Mean-field simulations

Based on a mean-field theory of a non-rotating turbulent convection [T. Elperin et al., Phys. Rev. E 66, 066305, (2002)], we perform mean-field simulations (MFS) of sheared convection that takes into account an effect of modification of the turbulent heat flux by the non-uniform large-scale motions. This effect is caused by the production of additional essentially anisotropic velocity fluctuations generated by tangling of the mean-velocity gradients by small-scale turbulent motions due to the influence of the inertial forces during the lifetime of turbulent eddies. These anisotropic velocity fluctuations contribute to the turbulent heat flux. As the result of this effect, there is an excitation of large-scale convective-shear instability, which causes the formation of large-scale semi-organized structures in the form of rolls. The lifetimes and spatial scales of these structures are much larger compared to the turbulent scales. By means of MFS performed for stress-free and no-slip vertical boundary conditions, we determine the spatial and temporal characteristics of these structures. Our study demonstrates that the modification of the turbulent heat flux by non-uniform flows leads to a strong reduction of the critical effective Rayleigh number (based on the eddy viscosity and turbulent temperature diffusivity) required for the formation of the large-scale rolls. During the nonlinear stage of the convective-shear instability, there is a transition from a two-layer vertical structure with two rolls in the vertical direction before the system reaches steady-state to a one-layer vertical structure with one roll after the system reaches steady state. This effect is observed for all effective Rayleigh numbers. We find that inside the convective rolls, the spatial distribution of the mean potential temperature includes regions with a positive vertical gradient of the potential temperature caused by the mean heat flux of the convective rolls. This study might be useful for understanding the origin of large-scale rolls observed in atmospheric convective boundary layers, as well as in numerical simulations and laboratory experiments.

A mechanism for coastal fog genesis at evening transition

AbstractTransitional changes in the atmospheric boundary layer (ABL) are known to facilitate the onset of terrestrial fog, which is defined as a condition with near‐surface visibility <1 km due to airborne water droplets. In particular, the evening transition from a daytime convective ABL to a night‐time stable ABL provides favorable conditions for fog. This article describes a local fog event observed during the evening transition at a Canadian islet in the north Atlantic known as Sable Island during the “Fog and Turbulence Interactions in the Marine Atmosphere (Fatima)” field campaign. The comprehensive set of data collected using a myriad of instruments covering a wide range of scales allowed identification of a novel mechanism underlying this fog event. Therein an ocean–land discontinuity created a flow regime consisting of several stacked boundary layers, interplay of which produced a thin low‐level cloud that then diffused downward to the surface, causing visibility reduction. This mechanism offers useful insights on the role of boundary layers, stratification, and turbulence in fog genesis over oceanic islands.

Spatially distributed atmospheric boundary layer properties in Houston – A value-added observational dataset

In 2022, Houston, TX became a nexus for field campaigns aiming to further our understanding of the feedbacks between convective clouds, aerosols and atmospheric boundary layer (ABL) properties. Houston’s proximity to the Gulf of Mexico and Galveston Bay motivated the collection of spatially distributed observations to disentangle coastal and urban processes. This paper presents a value-added ABL dataset derived from observations collected by eight research teams over 46 days between 2 June - 18 September 2022. The dataset spans 14 sites distributed within a ~80-km radius around Houston. Measurements from three types of instruments are analyzed to objectively provide estimates of nine ABL parameters, both thermodynamic (potential temperature, and relative humidity profiles and thermodynamic ABL depth) and dynamic (horizontal wind speed and direction, mean vertical velocity, updraft and downdraft speed profiles, and dynamical ABL depth). Contextual information about cloud occurrence is also provided. The dataset is prepared on a uniform time-height grid of 1 h and 30 m resolution to facilitate its use as a benchmark for forthcoming numerical simulations and the fundamental study of atmospheric processes.

Characterizing turbulence structures in convective and neutral atmospheric boundary layers via Koopman mode decomposition and unsupervised clustering

The atmospheric boundary layer (ABL) is a highly turbulent geophysical flow, which has chaotic and often too complex dynamics to unravel from limited data. Characterizing coherent turbulence structures in complex ABL flows under various atmospheric regimes is not systematically well established yet. This study aims to bridge this gap using large eddy simulations (LESs), Koopman theory, and unsupervised classification techniques. To this end, eight LESs of different convective, neutral, and unsteady ABLs are conducted. As the ratio of buoyancy to shear production increases, the turbulence structures change from roll vortices to convective cells. The quadrant analysis indicated that as this ratio increases, the sweep and ejection events decrease, and inward/outward interactions increase. The Koopman mode decomposition (KMD) is then used to characterize their turbulence structures. Our results showed that KMD can reveal non-trivial modes of highly turbulent ABL flows (e.g., transverse to the mean flow direction) and can reconstruct the primary dynamics of ABLs even under unsteady conditions with only ∼5% of the modes. We attributed the detected modes to the imposed pressure gradient (shear), Coriolis (inertial oscillations), and buoyancy (convection) forces by conducting novel timescale and quadrant analyses. We then applied the convolutional neural network combined with the K-means clustering to group the Koopman modes. This approach is displacement and rotation invariant, which allows efficiently reducing the number of modes that describe the overall ABL dynamics. Our results provide new insights into the dynamics of ABLs and present a systematic data-driven method to characterize their complex spatiotemporal patterns.

Lagrangian Stochastic Modeling of Stratified Atmospheric Boundary Layer

A single-column turbulence model for stratified atmospheric boundary layer (ABL), which solves the transport equations of turbulence probability density function (PDF) using a Lagrangian stochastic modeling (LSM) approach, is proposed in this study. This study adopts previously developed stochastic differential equations (SDEs) for particle velocity and temperature and extends the LSM to simulate inhomogeneous turbulence. The proposed LSM is tested for its ability to fully simulate statistics of inhomogeneous stratified turbulence. In the model, particles evolve by SDEs, and turbulence statistics are calculated by averaging the properties of particles. The model provides a full representation of turbulence PDF and simulates turbulent transport without any modeling assumption. The model performance is evaluated against large-eddy simulation (LES) results in the simulations of convective and stable ABL cases. For the convective ABL, LSM realistically simulates the entrainment process with the temperature and heat flux profiles that closely match with LES. The joint PDF simulated by LSM reproduces a curved and highly skewed shape, and some distinct features, like the asymmetric distribution of vertical velocity and the separation of the PDF in the entrainment zone, are simulated. LSM also reproduces the entrainment enhancement by wind shear in the simulation of sheared convective ABL. The LSM simulation of stable ABL predicts realistic turbulence intensity and mean field profiles, where Gaussian-like PDFs are simulated both in LSM and LES.

Saltation‐Induced Dust Emission of Dust Devils in the Convective Boundary Layer—An LES Study on the Meter Scale

AbstractDust devils are vertically oriented, columnar vortices that form within the atmospheric convective boundary layer (CBL) of dry regions. They are able to lift a sufficient amount of soil particles including dust to become visible and are considered as a potentially important dust source for the atmosphere. Mineral dust, a key component of atmospheric aerosols, influences the climate by affecting the radiation budget and cloud formation. Current estimates of the contribution of dust devils to the global, regional, and local dust release vary considerably from less than 1% to more than 50%. To address this uncertainty, we perform the highest resolved large‐eddy simulation (LES) study on dust emission in the CBL to date, using the PALM model system and the saltation‐based Air Force Weather Agency (AFWA) dust emission scheme. Our results show that under desert‐like conditions, dust devils are responsible for an average of 5% of regional dust emissions, with temporary maxima of up to 15%. This contrasts with previous measurement‐based (>35%) and LES‐based estimates (∼0.1%). Local emissions of dust devils (up to 10 mg m−2 s−1) are 1–3 orders of magnitude higher than the emission in the surroundings. This makes dust devils important for air quality and visibility. Additionally, our study reveals previously unknown large‐scale convective dust emission patterns. These patterns are tied to the CBL's cellular flow structure and are the main cause of dust release. Contrary to other studies, our findings clarify the important role of saltation‐induced dust emission.

Identification of Convective Boundary Layer Depth over the Great Lakes Region Using Aircraft Observations: A Comparison of Various Methods

Abstract This study evaluates the methods of identifying the height zi of 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 of zi is that above which the vertical velocity turbulent fluctuations are weak or absent. Thus, we use zi from the turbulence method as the “reference value” to which zi from 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−1 for soundings over land and 0.011 K m−1 for soundings over lake provided the estimates of zi that are most consistent with the turbulence method. The Rib threshold-based method, commonly used in numerical simulation studies, underestimated zi. Analyzing the methods’ performance on the averaging window zavg we recommend using zavg = 20 or 50 m for zi estimations 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 depth zi of the convective atmospheric boundary layer (CBL) strongly influences precipitation rates during lake-effect snowstorms (LES). However, various zi approximation 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 how zi from common estimation methods compare with “reference” zi derived 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 reference zi. 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.

Variant of the Local Similarity Theory and Approximations of Vertical Profiles of Turbulent Moments of the Atmospheric Convective Boundary Layer

Variant of the Local Similarity Theory and Approximations of Vertical Profiles of Turbulent Moments of the Atmospheric Convective Boundary Layer

Stability Dependence of the Turbulent Dissipation Rate in the Convective Atmospheric Boundary Layer

AbstractTurbulent dissipation rate (ɛ) is a crucial parameter in turbulence theory, and an essential component of higher‐order planetary boundary layer schemes for numerical weather prediction and climate models. It is most often modeled diagnostically based on the dissipation scaling ɛ ∝ e3/2/L, where e and L are the turbulence kinetic energy (TKE) and the size of the largest turbulent eddies, respectively. Utilizing three‐month‐long vertically‐extended observations accompanied by high resolution large‐eddy simulations, scaling‐based ɛ‐models are evaluated, focusing on their stability dependence under daytime convective conditions. The analysis uncovers biases in the parameterized ɛ profiles that cannot be corrected through tuning of model constants. The biases are attributed to the limited and even opposing stability dependence of the modeled dissipation length. Close examination reveals violation of the dissipation scaling by the inclusion of TKE associated with organized convection. A self‐similar dissipation length is obtained when only the isotropic component of TKE is considered.

Numerical study of aerodynamic performance of horizontal axis dual-rotor wind turbine under atmospheric boundary layers

The use of horizontal axis dual-rotor wind turbine (DRWT) is a new strategy to enhance the capture rate of wind energy and increase the performance of wind farms. An actuator line model (ALM) and large eddy simulation(LES) are introduced to investigate the aerodynamic performance of DRWT, and its effects on downstream turbines under convective atmospheric boundary layer (CBL) and neutral atmospheric boundary layer (NBL). The results conclude that the dominant vibration frequencies of the power production for the front rotor of DRWT are nearly the same as those of the single-rotor wind turbine (SRWT) under both NBL and CBL, but the vibration amplitude is slightly higher. The strength and dominant frequencies of yaw moment (Myaw) for the front rotor of DRWT are almost the same as those of the SRWT in both NBL and CBL flows, while the results of blade-root out-of-plane bending moment (Moop) are different. There are obvious differences in wake development and wake meandering between the DRWT and the SRWT. For three turbines cases, the total power production is increased respectively by 3.3% and 3% under NBL and CBL at a tandem spacing of 5D when the DRWT is placed in the first row, while the results increase to 5.5% and 4.4% at a tandem spacing of 9D. The stability of Myaw and Moop of second-row turbine located 5D downstream behind the DRWT under both NBL and CBL and Moop of all turbines located behind the DRWT in three tandem spacings (5D, 7D, 9D) under CBL are deteriorated compared with those located behind the SRWT.

Analyses of Spatial Correlation and Coherence in ABL Flow with a Fleet of UAS

The spatial structures of turbulent flow in the atmospheric boundary layer (ABL) are complex and diverse. Multi-point spatial correlation measurements can help improve our understanding of these structures and their statistics. In this context, we investigate Taylor’s hypothesis and the statistics of spatial structures on the microscale. For the first time, simultaneous horizontally distributed wind measurements with a fleet of 20 quadrotor UAS (unmanned aerial systems) are realized. The measurements were taken at different heights and under different atmospheric conditions at the boundary layer field site in Falkenberg of the German Meteorological Service (DWD). A horizontal flight pattern has been specifically developed, consisting of measurements distributed along and lateral to the mean flow direction with separation distances of 5ldots 205 m. The validity of Taylor’s hypothesis is studied by examining the cross-correlations of longitudinally distributed UAS and comparing them with the autocorrelations of single UAS. To assess the similarity of flow structures on different scales, the lateral and longitudinal coherence of the streamwise velocity component is examined. Two modeling approaches for the decay of coherence are compared. The experimental results are in good agreement with the model approaches for neutral atmospheric conditions, whereas in stable and convective ABL, the exponential approaches are not unconditionally valid. The validation results and the agreement with the literature on coherence in the ABL underline the potential of the UAS fleet for the purpose of spatial turbulence measurements.

A Wavenumber–Frequency Spectrum Model for Sheared Convective Atmospheric Boundary Layer Flows

Abstract Parameterizing turbulence in the atmospheric boundary layer as a function of space and time is essential for weather and climate models. Here, we explore a model for wavenumber–frequency spectra based on a linear random advection approach to characterize sheared convective atmospheric boundary layer flows. Building on previous works, we obtain the wavenumber–frequency spectrum as a product of the wavenumber spectrum and a Gaussian frequency distribution, whose mean and variance are given by the mean advection and random sweeping velocities, respectively. The applicability of the model is tested with direct numerical simulation data in the mixed layer and the entrainment zone for the streamwise and vertical velocity components and buoyancy. To obtain a fully analytical model, we propose using a von Kármán wavenumber spectrum parameterized by the characteristic variances and integral length scales. These parameters are height dependent and vary considerably with the relative balance of buoyancy and shear forces. The introduced analytical model relies on fitting parameters obtained from numerical data in the relevant range of scales. The comparison of the von Kármán–based spectra for velocity and buoyancy to simulation results shows that the main features of the measured spectra are captured by the model.

Evolution and Features of Dust Devil‐Like Vortices in Turbulent Rayleigh‐Bénard Convection—An Experimental Study

AbstractWe present an experimental study simulating atmospheric dust devils in a controlled laboratory experiment. The experimental facility, called the “Barrel of Ilmenau” (www.ilmenauer-fass.de) represents a classical Rayleigh‐Bénard set‐up and is believed to model the phenomena in a convective atmospheric boundary layer fairly well. Our work complements and extends the numerical work of Giersch and Raasch (2021)https//doi.org/10.1029/2020jd034334 by experiments. Dust devils are thermal convective vortices with a vertical axis of rotation visualized by entrained soil particles. They evolve in the convective atmospheric boundary layer and are believed to substantially contribute to the aerosol transport into the atmosphere. Thus, their evolution, size, lifetime, and frequency of occurrence are of particular research interest. Extensive experimental studies have been conducted by field measurements and laboratory experiments so far. Beyond that, our study is the first attempt of Rayleigh‐Bénard convection (RBC) in air to investigate dust devil‐like vortices in a laboratory experiment. Up to now, this set‐up mimics the natural process of dust devil evolution as closest to reality. The flow measurement was carried out by particle tracking velocimetry using neutrally buoyant soap bubbles. We initially identified dust devil‐like vortices by eye from the Lagrangian velocity field, and in a later, more sophisticated analysis by a specific algorithm from the corresponding Eulerian velocity field. We analyzed their frequency of occurrence, observation time, and size. With our work, we could demonstrate that turbulent RBC is an appropriate model to mimic the natural process of the evolution of dust devils in the convective atmospheric boundary layer without artificial stimulation.

Large‐eddy simulation of soil moisture heterogeneity‐induced secondary circulation with ambient winds

AbstractLand surface heterogeneity in conjunction with ambient winds influences the convective atmospheric boundary layer by affecting the distribution of incoming solar radiation and forming secondary circulations. This study performed coupled large‐eddy simulation (ICON‐LEM) with a land surface model (TERRA‐ML) over a flat river corridor mimicked by soil moisture heterogeneity to investigate the impact of ambient winds on secondary circulations. The coupled model employed double‐periodic boundary conditions with a spatial scale of 4.8 km. All simulations used the same idealized initial atmospheric conditions with constant incident radiation of 700 W⋅m−2 and various ambient winds with different speeds (0 to 16 m⋅s−1) and directions (e.g., cross‐river, parallel‐river, and mixed). The atmospheric states are decomposed into ensemble‐averaged, mesoscale, and turbulence. The results show that the secondary circulation structure persists under the parallel‐river wind conditions independently of the wind speed but is destroyed when the cross‐river wind is stronger than 2 m⋅s−1. The soil moisture and wind speed determine the influence on the surface energy distribution independent of the wind direction. However, secondary circulations increase advection and dispersive heat flux while decreasing turbulent energy flux. The vertical profiles of the wind variance reflect the secondary circulation, and the maximum value of the mesoscale vertical wind variance indicates the secondary circulation strength. The secondary circulation strength positively scales with the Bowen ratio, stability parameter (−Zi/L), and thermal heterogeneity parameter under cross‐river wind and mixed wind conditions. The proposed similarity analyses and scaling approach provide a new quantitative perspective on the impact of the ambient wind under heteronomous soil moisture conditions on secondary circulation.

Energy- and flux-budget theory for surface layers in atmospheric convective turbulence

The energy- and flux-budget (EFB) theory developed previously for atmospheric stably stratified turbulence is extended to the surface layer in atmospheric convective turbulence. This theory is based on budget equations for turbulent energies and fluxes in the Boussinesq approximation. In the lower part of the surface layer in the atmospheric convective boundary layer, the rate of turbulence production of the turbulent kinetic energy (TKE) caused by the surface shear is much larger than that caused by the buoyancy, which results in three-dimensional turbulence of very complex nature. In the upper part of the surface layer, the rate of turbulence production of TKE due to the shear is much smaller than that caused by the buoyancy, which causes unusual strongly anisotropic buoyancy-driven turbulence. Considering the applications of the obtained results to the atmospheric convective boundary-layer turbulence, the theoretical relationships potentially useful in modeling applications have been derived. The developed EFB theory allows us to obtain a smooth transition between a stably stratified turbulence to a convective turbulence. The EFB theory for the surface layer in a convective turbulence provides an analytical expression for the entire surface layer including the transition range between the lower and upper parts of the surface layer, and it allows us to determine the vertical profiles for all turbulent characteristics, including TKE, the intensity of turbulent potential temperature fluctuations, the vertical turbulent fluxes of momentum and buoyancy (proportional to potential temperature), the integral turbulence scale, the turbulence anisotropy, the turbulent Prandtl number, and the flux Richardson number.

The Impacts of Wind Shear on Spatial Variation of the Meteorological Element Field in the Atmospheric Convective Boundary Layer Based on Large Eddy Simulation

As wind shear increases, the quasi-two-dimensional structure of flows becomes more significant in the convective boundary layer (CBL), indicating that wind shear plays an essential role in the variation of the field of atmospheric flow. Therefore, sensitive numerical experiments based on Large Eddy Simulation (LES) techniques were conducted to comprehensively investigate the effects of wind shear on the spatial variations in the velocity and potential temperature (θ) horizontal fields. Under the constant surface heat flux condition, the main findings are summarized. Firstly, in the CBL, the variances of the streamwise velocity (u), cross-stream velocity (v), and θ enhance as wind shear increases, whereas the variance of vertical velocity (w) is insensitive to wind shear. Secondly, in the CBL, with increasing wind shear, low-wavenumber Power Spectrum Densities (PSDs) of u, v, w, and θ increase significantly, suggesting that the increasing wind shear always enhances the large-scale motions of the atmosphere (i.e., low-wavenumber PSD). Therefore, it is more likely that some mesoscale weather processes will be triggered. Thirdly, generally, in the high-wavenumber range, with increasing wind shear, the PSDs of u, v, and θ increase slightly, whereas the PSD of w decreases slightly. This study provides a new perspective for understanding the role of wind shear in the spatial variations of the horizontal fields of meteorological elements under the same conditions of surface heat flux.

Impact of turbulence on aeolian particle entrainment: results from wind-tunnel experiments

Abstract. We hypothesize that large eddies play a major role in the entrainment of aeolian particles. To test this, wind-tunnel experiments are carried out to measure the particle entrainment rate for various sizes and flow conditions. Wind-tunnel flows are usually neutrally stratified with no large eddies, which are typically seen in convective atmospheric boundary layers. Here, a novel technique is applied, by deploying a piece of randomly fluttering cloth, to generate large eddies similar to convective eddies, here referred to as quasi-convective turbulence. The characteristics of quasi-convective turbulence are analyzed with respect to neutral turbulence in the Monin–Obukhov similarity framework, and the probability distributions of surface shear stress are examined. We show that for a given mean flow speed and in comparison with neutral flow conditions, quasi-convective turbulence increases the surface shear stress and alters its probability distribution and hence substantially enhances the entrainment of sand and dust particles. Our hypothesis is thus confirmed by the wind-tunnel experiments. We also explain why large eddies are important to aeolian entrainment and transport.