Turbulence in the atmospheric boundary layer (ABL) is one of the most challenging and important phenomena in geophysical fluid dynamics. The ABL is the lowest kilometer or so of the atmosphere, directly influenced by the Earth’s surface through friction, heating, and moisture exchange. Turbulent motions in this layer govern the vertical transport of momentum, heat, water vapor, and pollutants, thus shaping weather, climate, and air quality. Recent years have witnessed remarkable progress in measuring, simulating, and conceptualizing ABL turbulence, driven by advances in observational technology, computational power, and theoretical frameworks. This article reviews these developments, highlights key discoveries, and discusses their implications for weather prediction, climate modeling, and environmental management.

Historical Perspective

The systematic study of atmospheric turbulence began in the early twentieth century, when Osborne Reynolds introduced the concept of Reynolds decomposition and the Reynolds-averaged Navier–Stokes (RANS) equations. His work for the first time provided a mathematical framework to separate mean flow from turbulent fluctuations, but the closure problem—how to represent turbulent fluxes in terms of mean quantities—remained intractable. Early empirical approaches relied on mixing-length theory, akin to molecular diffusion, but they could not capture the complexity of turbulent flows in a stratified atmosphere.

A major breakthrough came in the 1940s when Andrey Kolmogorov proposed his theory of turbulence, describing the energy cascade from large to small scales. Kolmogorov’s −5/3 law for the energy spectrum and his universal similarity hypotheses became cornerstones of turbulence theory. However, the atmospheric boundary layer posed additional challenges: it is strongly influenced by buoyancy (stability), surface roughness, and the Earth’s rotation. In the 1950s and 1960s, Monin and Obukhov developed similarity theory, which describes the vertical profiles of wind, temperature, and turbulent fluxes in the surface layer using a set of dimensionless stability parameters. This framework remains widely used in observational analysis and model parameterizations today.

Despite these theoretical advances, early field experiments were limited by instrument capabilities. Slow-response anemometers and thermometers could only capture large-scale fluctuations, and spatial coverage was sparse. The turbulence community relied on a handful of dedicated sites, such as the Kansas and Minnesota experiments in the United States, which provided the datasets that validated and refined Monin–Obukhov similarity theory. Yet the highly intermittent and three-dimensional nature of ABL turbulence demanded new approaches.

Recent Advances in Observational Technology

Remote Sensing

The last two decades have seen a revolution in observing turbulent motions. Coherent Doppler lidar systems can now measure wind velocity profiles with high temporal and spatial resolution, revealing details of turbulent eddies from meter to kilometer scales. Scanning lidars can map the three-dimensional structure of the boundary layer, detecting coherent features such as roll vortices, convective plumes, and shear-induced turbulence. Similarly, frequency-modulated continuous-wave (FMCW) radars and sodars provide complementary information on turbulence intensity and boundary layer height. These instruments can operate continuously and over large areas, enabling studies of turbulence in remote or complex terrain where in situ towers are impractical.

In Situ Instruments

While remote sensing offers unparalleled spatial coverage, in situ measurements remain essential for high-frequency flux measurements. Sonic anemometers now resolve all three wind components and temperature fluctuations at 10–20 Hz, allowing direct computation of turbulent fluxes via eddy covariance. The advent of fast-response gas analyzers for CO₂, H₂O, and other trace gases has transformed the study of land–atmosphere exchange. Unmanned aerial vehicles (UAVs) equipped with lightweight turbulence probes can now sample the boundary layer from near the surface up to several hundred meters, filling a critical gap between tower and airborne measurements. Research aircraft like the NSF/NCAR C-130 and the NASA Global Hawk carry sophisticated turbulence and flux instrumentation, probing large-scale coherent structures and spatial heterogeneity.

Computational Advances: Large Eddy Simulation and Direct Numerical Simulation

Parallel to observational progress, computational fluid dynamics has transformed our ability to simulate ABL turbulence. Large Eddy Simulation (LES) solves the filtered Navier–Stokes equations, resolving the energy-containing turbulent eddies while parameterizing the smaller, subgrid-scale motions. Modern LES codes can simulate domains of tens of kilometers with grid spacings of a few meters, capturing the multi-scale nature of the boundary layer—from surface-layer streaks to deep convective cells. LES has been used to study stably stratified boundary layers, cloud-topped boundary layers, and urban canopies, yielding insights that are impossible from observations alone due to incomplete coverage.

Direct Numerical Simulation (DNS), which resolves all scales down to the viscous dissipation range, remains limited to relatively low Reynolds numbers and small domains. Recent efforts have extended DNS to slightly unstable boundary layers and shear-driven turbulence, providing a benchmark for subgrid models in LES. The synergy between LES, DNS, and observations now forms a powerful triangle for understanding turbulence physics and developing improved parameterizations.

Key Discoveries and Theoretical Advances

Coherent Structures and Their Role in Transport

One of the most significant insights from high-resolution simulations and field campaigns is the ubiquity of coherent structures. In the convective boundary layer, rising thermals organize into large-scale roll vortices or spoke-like patterns. In the neutral and stable boundary layer, low-speed streaks and hairpin vortices dominate momentum transfer. These structures are not random; they exhibit preferred scales, orientations, and lifetimes. Their presence means that turbulent transport is highly intermittent and organized. For instance, a single gust of wind from a coherent downdraft can transport as much momentum in a few seconds as mean shear does over minutes. This has profound implications for parameterizing turbulent fluxes in large-scale models, which traditionally assume a continuous, isotropic mixing process.

Atmospheric Stability and Turbulence Regimes

Research in the past decade has greatly refined our understanding of how static stability modulates turbulence. The surface layer follows Monin–Obukhov similarity, but above the surface, the stably stratified boundary layer is often characterized by intermittent turbulence, low-level jets, and gravity waves. Observations from the SHEBA project in the Arctic and the FLOSSII experiment in the U.S. Great Plains have revealed that turbulence in the stable boundary layer can collapse and suddenly burst, triggered by shear instabilities or topographically induced waves. These phenomena are poorly represented in operational models. The concept of the “critical Richardson number” — below which turbulence is sustained — has been replaced by a more nuanced view involving hysteresis and subcritical transitions.

Surface Heterogeneity and Urban Effects

Real surfaces are rarely homogeneous; land use, vegetation, topography, and urbanization create strong spatial variations in roughness, albedo, and moisture availability. Recent field campaigns like LAFE (Land-Atmosphere Feedback Experiment) and experiments over cities (e.g., BUBBLE in Basel) have shown that surface heterogeneity generates secondary circulations that can dominate day-to-day turbulence intensity. In urban areas, the roughness of buildings and the heat island effect create deep turbulent layers with distinct scaling laws. These findings are now being incorporated into next-generation land surface models and urban parameterizations.

Intermittency and Non-stationarity

Another key discovery is that ABL turbulence is often non-stationary and intermittent, especially under stable conditions and during transitions (e.g., morning and evening). Traditional spectral and correlation methods assume stationarity, but new analytical tools like wavelet transforms, Hilbert–Huang transforms, and Lagrangian particle statistics are revealing rich intermittent dynamics. This has implications for flux averaging periods, gap-filling of eddy covariance data, and the interpretation of field measurements.

Implications for Weather and Climate Modeling

Accurate representation of ABL turbulence is essential for numerical weather prediction (NWP) and climate models. Turbulence controls surface fluxes of heat, moisture, and momentum, which feed back to the large-scale circulation. Errors in boundary layer schemes can lead to biases in near-surface temperature, wind, and precipitation, as well as in the timing of fog, cloud formation, and sea breeze circulations.

Recent advances have led to new parameterizations that incorporate the effects of coherent structures, stability dependence, and surface heterogeneity. The YSU (Yonsei University) and MYNN (Mellor–Yamada–Nakanishi–Niino) schemes are widely used in operational models, but they still rely on first-order closure or simple turbulent kinetic energy (TKE) budgets. Many research groups are now developing scale-aware parameterizations that adjust based on model resolution, as well as unified schemes that seamlessly handle convective and stable conditions. Some models are even moving toward explicit LES-like representations within the boundary layer, a technique known as “cloud-resolving” or “convection-permitting” modeling.

For climate projections, the correct handling of ABL turbulence is critical. For instance, the strength of the land–atmosphere coupling influences the severity of heatwaves and droughts. In the tropical marine boundary layer, shallow cumulus clouds—which depend on turbulent transport—have a large impact on the global energy budget. Improving parameterizations of these clouds remains a top priority for the climate modeling community.

Future Directions

Integration of Machine Learning

Machine learning (ML) is rapidly becoming a powerful tool for understanding and parameterizing turbulence. Neural networks, random forests, and deep learning models have been trained on LES and observational data to predict turbulent fluxes, eddy diffusivities, and subgrid-scale stress. ML-based closures can capture non-local effects and conditional averages that traditional algebraic closures miss. However, challenges remain: ML models must respect physical constraints (e.g., conservation laws, realizability) and generalize to out-of-sample conditions. Hybrid approaches that combine physics-based models with data-driven corrections may offer the best path forward.

Multi-Scale Campaigns and Global Networks

Future field experiments will emphasize integrated, multi-platform measurements that span from the surface to the top of the ABL and extend over large spatial domains. The upcoming Atmospheric Radiation Measurement (ARM) User Facility’s “Boundary Layer and Cloud” campaigns will combine lidar, radar, UAVs, and balloon soundings to capture turbulence across multiple scales. Global networks of flux towers, like FLUXNET, continue to grow and provide long-term datasets that can be used to test parameterizations under diverse climate regimes. Such data will also serve as training sets for ML models.

High-Resolution Global Modeling

With exascale computing on the horizon, global models with grid spacings of a few kilometers will become feasible. At these resolutions, many boundary layer features—including large convective rolls and shallow cumulus—will be partially resolved, reducing the reliance on parameterizations. However, the smallest scales (the inertial subrange and below) will still need stochastic or deterministic closure. Research into scale-adaptive and stochastic parameterizations is therefore a key frontier.

Extreme Environments and Climate Change

Understanding turbulence in a warming climate is crucial. The Arctic boundary layer is undergoing rapid changes—sea ice loss alters surface roughness and heat fluxes, leading to new patterns of turbulence and cloud formation. Similarly, future cities with higher temperatures and denser structures will experience modified urban boundary layers. Observations and models must be extended to these emerging regimes to ensure that predictions remain reliable.

Conclusion

The study of turbulence in the atmospheric boundary layer has advanced from a largely phenomenological discipline to a quantitative, process-oriented science. High-resolution observations from lidar, UAVs, and flux towers, coupled with sophisticated simulations using LES and DNS, have revealed the centrality of coherent structures, the importance of stability and surface heterogeneity, and the intermittent nature of turbulent transport. These insights are being woven into improved parameterizations for weather and climate models, while machine learning and exascale computing promise further breakthroughs. As the ABL continues to respond to global change, sustained investment in turbulence research will be essential for accurate environmental prediction and informed decision-making.

External resources for further reading: