Quantitative knowledge of the surface energy balance is essential for the prediction of weather and climate. However, a multitude of studies from around the world indicate that the turbulent heat fluxes are generally underestimated using eddy-covariance measurements, and hence, the energy balance is not closed. This energy-balance-closure problem, which has been heavily covered in the literature for

We apply well-established flux–gradient relationships to deduce the aerodynamic and radiative properties of a winter wheat crop, using a neglected 1971 dataset (hourly averages), only recently resurrected as part of a historical review of precision CO2 measurements in Australia (Pearman et al. in Hist Rec Aust Sci 28:111–125, 2017). The aerodynamic roughness length (seasonal variation between 0.07

A computational steady Reynolds-averaged Navier–Stokes model is developed for a neutrally-stratified unidirectional flow in the atmospheric boundary layer (ABL) over generic two-dimensional cosine-shaped hills. Simulations are performed for a hill of a constant height and three different hill slopes. The hill surface is treated as smooth and aerodynamically rough. Atmospheric-boundary-layer inflow

We use large-eddy simulations with an immersed boundary method to study the performance of wind turbines and wind farms in hilly terrain. First, we analyze the performance of wind turbines in the vicinity of a two-dimensional hill. For turbines that are significantly taller than the hill, the performance improves as the flow speeds up over the hill. For turbines that have approximately the same or

All notation is as defined in the original article.

A synthetic-turbulence and temperature-fluctuation-generation method is developed and embedded in large-eddy simulations to investigate the effects of weak stable stratification (i.e. Richardson number \(Ri\le 1\)) on turbulence and dispersion following a simulated rural-to-urban transition. The modelling approach is validated by comparing predictions of mean velocity, turbulent stresses, and point-source

There are two models for the surface layer in the convective boundary layer (CBL). First is the standard Monin–Obukhov similarity theory, and second is the McNaughton–Laubach model (Laubach and McNaughton 2009, Boundary-Layer Meteorology 133:219–252; hereafter MNL model) developed based on the complex dynamical system approach to address characteristics of the unstable surface layer. The fundamental

We study the convective boundary layer (CBL) through low-order topological properties of updrafts and downdrafts, that is, based solely on the sign of the vertical velocity. The geometric representation of the CBL as a pair of two-dimensional cubical complexes, one each for updrafts and downdrafts, is exemplarily obtained from two simulations of the CBL, a realistic daily cycle and an idealized quasi-steady

A comparison between the spectral microphysics cloud scheme (MiStra) and the parametrized cloud scheme (PaStra) is presented. The main feature of MiStra consists of the treatment of aerosol particles, cloud droplets, and drizzle particles in a joint two-dimensional particle size distribution, whereas PaStra consists of a two-moment scheme for cloud droplets combined with a one-moment scheme for drizzle

The genesis and evolution of a sea-breeze front associated with a shallow convection episode at Fuerteventura Island (Canaries archipelago) are investigated using the Weather Research and Forecasting (WRF) numerical model. Three local and two non-local planetary boundary layer (PBL) parametrization schemes are used, and results of the numerical simulations are compared with observations from the local

Large-eddy simulations are widely used to study flows in the atmospheric boundary layer. As atmospheric boundary-layer flows of different atmospheric stratification have very different flow characteristics on different length scales, well-resolved simulations of these flows require very different meshes. The Parallelized Large-Eddy Simulation Model combined with a realizable dynamic subgrid model is

Recent research indicates that non-orthogonal sonic anemometers underestimate vertical wind velocity and consequently eddy-covariance fluxes of mass and energy. Whether this is a general problem among all non-orthogonal sonic anemometers, including those calibrated for flow-shadowing effects, is unknown. To investigate this, we test two sonic anemometer designs, orthogonal (3Vx-probe, Applied Technologies

Widespread and persistent summer multi-day episodes of dense wildfire smoke affected western Canada in 2017 and 2018. These events often occurred under otherwise clear-sky, anticyclonic weather conditions and can have significant impacts on surface temperatures, air quality, and surface radiation and energy budgets. Based on upward-pointing lidar observations, vertical temperature soundings and numerical

A method is proposed to estimate the threshold wind speed for dust emissions as a function of soil moisture in arid regions. This method, which is applicable at the local scale, employs a model of the surface heat budget to estimate the spatial distribution of thermal inertia-derived soil moisture (TISM) and an analytical footprint model to estimate dust source areas. It incorporates readily available

The present laboratory study is focused on the role of convective rolls in the enhancement of the heat flux from the sea and triggering of the process of rapid intensification of tropical cyclones. The appearance of coherent convective structures such as thermals and rolls is registered by different optical techniques and temperature measurements. Two-dimensional velocity fields are used for the investigation

A semi-empirical approach based on surface-renewal theory for estimating the friction velocity is tested for measurements taken in the inertial sublayer. For unstable cases, the input requirements are the mean wind speed and the high-frequency trace (10 or 20 Hz) of the air or sonic temperature. The method has been extended to traces of water vapour (H2O) and carbon dioxide (CO2) concentrations. For

A number of studies have indicated that non-stationarity of the wind field over the sea significantly disrupts the equilibrium between the wind, stress, and wave fields. However, no studies have systematically examined the impact of non-stationarity using a large dataset. Here, we examine the effect of non-stationarity of the wind field on the stress by analyzing six years of turbulent flux data from

A central-moments-based lattice Boltzmann model for large-eddy simulation of neutrally-stratified turbulent flows is described. Through comparative simulations of the airflow within and above a homogeneous plant canopy, the performance of the model is evaluated with respect to a conventional large-eddy-simulation model based on the incompressible Navier–Stokes equations. Simulated turbulence statistics

The atmospheric boundary layer undergoes transitions between stable and convective states. Over land, in undisturbed conditions, these transitions occur daily in the morning and late afternoon or early evening. Though less well studied and presenting more challenges than the fully stable and fully convective states, such transitions have been the subject of growing interest over the last few decades

Eddy-covariance fluxes serve as an essential benchmark for Earth system models and remote sensing data. However, two challenges prevent model-data intercomparisons from fully utilizing eddy-covariance fluxes. The first challenge stems from the differing and variable spatial representativeness of the eddy-covariance measurements, or footprint bias and transience. The second originates from the phenomenon

The scalar front generated by the horizontal self advection of a dipolar vortex through a modest scalar gradient is investigated. This physical scenario is an idealization of the emergence of strong temperature ramps in the stable atmospheric boundary layer. The proposed mechanism is discussed and a two-dimensional analogy is studied in depth using direct numerical simulation. More specifically, the

We investigate the temporal evolution of the planetary boundary-layer (PBL) height over the basin of Athens, Greece, during a 6-year period (2011–2016), using data from a Raman lidar system. The range-corrected lidar signals are selected around local noon (1200 UTC) and midnight (0000 UTC), for a total of 332 cases: 165 days and 167 nights. In this dataset, the extended-Kalman filtering technique is

Knowledge of the effects of small-scale fluctuations in temperature on light transmission in the atmosphere is necessary for the calibration of remote sensing instruments as well as for the understanding of turbulent heat transport in the atmospheric boundary layer. Recent developments in small unmanned aircraft systems (sUAS) have allowed for direct, spatial in situ estimation of temperature in the

While the increasing availability of computational power is enabling finer grid resolutions in numerical-weather-prediction models, representing land–atmosphere exchange processes remains challenging. This partially results from the fact that land-surface heterogeneity exists at all spatial scales, and its variability does not necessarily ‘average’ out with decreasing size. The work presented here

A comprehensive investigation is carried out to establish best practice guidelines for the modelling of statistically steady-state non-neutral urban boundary layers (UBL) using large-eddy simulation (LES). These steady-state simulations enable targeted studies under realistic non-neutral conditions without the complications associated with the inherently transient nature of the UBL. An extensive set

Cold-air outbreaks (CAO) lead to intense air–sea interactions, the appropriate representation of which are fundamental for climate modelling and numerical weather forecasting. We analyze a CAO event with low-level wind speeds of approximately 25 m s\(^{-1}\) observed in the north-western Mediterranean Sea. The marine atmospheric boundary layer (MABL) was sampled with an aircraft equipped for turbulence

The growth of the daytime convective atmospheric boundary layer, which plays a pivotal role in the vertical mixing and dispersal of water vapour and pollutants, is modulated by cloud radiative effects. Assessment of this cloud effect is sparse in any geographical region and non-existent over tropical coastal regions. We investigate the effect of clouds on the diurnal evolution of boundary-layer height

We assess the performance of turbulence closures of varying degrees of sophistication in the prediction of the mean flow and the thermal fields in a neutrally-stratified Ekman layer. The Reynolds stresses that appear in the Reynolds-averaged momentum equations are determined using both eddy-viscosity and complete differential Reynolds-stress-transport closures. The results unexpectedly show that the

Abstract Field experiments were undertaken in the summer of 2015 in Manchester, UK, to investigate the dispersion behaviour and infiltration into buildings of gas-phase pollutants over horizontal distances of 1–5 km. Inert cyclic perfluorocarbon tracers were released for 15 min at either one or three release points and samples taken in locations indoors and outdoors up to 2 km downwind. Background

Abstract High-frequency measurements are available at five heights within and above a row-gap trellised vineyard located on a 7\(^{\circ }\) slope in the Southern Okanagan Valley, British Columbia, Canada. During a 3-week campaign in July 2016, approximately 17% of the nocturnal conditions exhibit drainage flow along the local slope. Drainage conditions are characterized by temperature inversions beginning

Abstract Four cases of an overlying inversion imposed on a stable boundary layer are investigated, extending the earlier work of Hancock and Hayden (Boundary-Layer Meteorol 168:29–57, 2018), where no inversion was imposed. The inversion is imposed to one or other of two depths within the layer: midway or deep. Four cases of changed surface condition are also investigated, and it is seen that the surface

Abstract The Hong Kong Observatory conducted six flights in the atmospheric boundary layer of five tropical cyclones: tropical storm Jebi (1309), typhoon Kalmaegi (1415), severe tropical storm Linfa (1510), typhoon Mujigae (1522), and severe typhoon Nida (1604). Three-dimensional wind data with a 20-Hz sampling rate are available for a height range of 500–670 m, with the mean wind speed from these

Abstract We investigate the flow characteristics around step-up street canyons with various building aspect ratios (ratio of along-canyon building length to street-canyon width, and upwind building height to downwind building height) using a computational fluid dynamics (CFD) model. Simulated results are validated against experimental wind-tunnel results, with the CFD simulations conducted under the

Abstract We propose a new phenomenological model to represent the impact of wind-waves on the dissipation of turbulence kinetic energy near the sea surface. In this model, the momentum flux at a given height results from the averaged contribution of eddies attached to the sea surface whose sizes are related to the surface geometry. This yields a coupling between long wind-waves and turbulence at heights

Abstract Two-wavelength scintillometer systems can provide much needed measurements of area-averaged sensible and latent heat fluxes. However, these devices rarely have been deployed on canopy-covered complex terrain, and never in the circumpolar boreal biome, where large-scale fluxes are essential to hydroclimate modellers. We present a comparison of fluxes measured above a boreal-forested valley

Abstract The eddy-covariance (EC) technique is used to determine mass and energy fluxes between the Earth’s surface and the lower atmosphere at high temporal resolution. Despite the frequent and successful use of the EC technique at terrestrial sites, its application over water surfaces is rare. We present one season of EC measurements conducted on the Rappbode Reservoir, Germany’s largest drinking

Abstract Large-eddy simulation is used to examine the afternoon and early evening transition of the atmospheric boundary layer (ABL) over heterogeneous surfaces with respect to the decay of the turbulence kinetic energy (TKE). The heterogeneous surface sensible heat flux is prescribed based on the inverse transform of the Fourier spectrum whose slope is controlled to behave as \( \kappa^{0} \), \(

Abstract A resistance model for organized rough walls of varying geometry in fully-developed turbulent boundary layers is proposed. For a given flow, the shear stress on rough walls varies with the surface geometry. The model normalizes the roughness shear stress by its theoretical maximum. Based on an analysis of measurable drag on individual roughness elements, the expression for the normalized roughness

We investigate the effects of condensation and liquid water loading on the stably stratified surface layer, with an eye towards understanding the influence of turbulent mixing on fog formation. Direct numerical simulations of dry and moist open-channel flows are conducted, where in both a constant cooling rate is applied at the ground to mimic longwave radiative cooling. Depending on the cooling rate

We propose a single generalized function for dimensionless wind and temperature gradients in stable and unstable conditions based on Monin–Obukhov similarity theory. The proposed function may be matched with the currently accepted universal functions for the dimensionless wind-speed and temperature gradients using the empirical coefficients a and b to yield differences within ± 14% for both wind and

Considering that the characteristics of the heat balance between an isolated tree and the urban atmosphere have not yet been sufficiently clarified, we quantify both the sensible heat flux (\( H_{T} \)) and latent heat flux (\( lE_{T} \)) to and from an isolated Z. Serrata in the Japanese summer. To estimate the whole-tree transpiration rate (ET) and \( lE_{T} \) values, we apply a previously developed

When canopy flows are horizontally averaged to obtain mean profiles, the averaging operation can be defined either as an intrinsic average, normalized by the variable fluid volume, or as a superficial average, normalized by the total volume including solid canopy elements. Properties of spatial averages have been explored extensively in the context of flow through plant canopies, albeit with the assumption

Conventional in situ observations of meteorological variables are restricted to a limited number of levels near the surface, with the lowest observation often made around 1-m height. This can result in missed observations of both shallow fog, and the initial growth stage of thicker fog layers. At the same time, numerical experiments have demonstrated the need for high vertical grid resolution in the

The reduction in visibility that accompanies fog events presents a hazard to human safety and navigation. However, accurate fog prediction remains elusive, with numerical methods often unable to capture the conditions of fog formation, and observational methods having high false-alarm rates in order to obtain high hit rates of prediction. In this work, 5 years of observations from the Cabauw Experimental

The scaled standard deviations of temperature and humidity are investigated in complex terrain. The study area is a steep Alpine valley, with six measurement sites of different slope, orientation and roughness (i-Box experimental site, Inn Valley, Austria). Examined here are several assumptions forming the basis of Monin-Obukhov similarity theory (MOST), including constant turbulence fluxes with height

The performance of an atmospheric single-column model (SCM) is studied systematically for stably-stratified conditions. To this end, 11 years (2005-2015) of daily SCM simulations were compared to observations from the Cabauw observatory, The Netherlands. Each individual clear-sky night was classified in terms of the ambient geostrophic wind speed with a [Formula: see text] bin-width. Nights with overcast

Scalar dispersion from ground-level sources in arrays of buildings is investigated using wind-tunnel measurements and large-eddy simulation (LES). An array of uniform-height buildings of equal dimensions and an array with an additional single tall building (wind tunnel) or a periodically repeated tall building (LES) are considered. The buildings in the array are aligned and form long streets. The sensitivity

The correct simulation of the atmospheric boundary layer (ABL) is crucial for reliable weather forecasts in truly complex terrain. However, common assumptions for model parametrizations are only valid for horizontally homogeneous and flat terrain. Here, we evaluate the turbulence parametrization of the numerical weather prediction model COSMO with a horizontal grid spacing of Δ x = 1.1 km for the Inn

Simultaneous particle-image velocimetry and laser-induced fluorescence combined with large-eddy simulations are used to investigate the flow and pollutant dispersion behaviour in a rural-to-urban roughness transition. The urban roughness is characterized by an array of cubical obstacles in an aligned arrangement. A plane fence is added one obstacle height h upstream of the urban roughness elements

A generalized form of a recently developed minimum dissipation model for subfilter turbulent fluxes is proposed and implemented in the simulation of thermally stratified atmospheric boundary-layer flows. Compared with the original model, the generalized model includes the contribution of buoyant forces, in addition to shear, to the production or suppression of turbulence, with a number of desirable

Nocturnal boundary-layer phenomena in regions of complex topography are extremely diverse and respond to a multiplicity of forcing factors, acting primarily at the mesoscale and microscale. The interaction between different physical processes, e.g., drainage promoted by near-surface cooling and ambient flow over topography in a statically stable environment, may give rise to special flow patterns,

Carbon dioxide, latent and sensible energy fluxes were measured by means of the eddy covariance method above a mountain meadow situated on a steep slope in the Stubai Valley/Austria, based on the hypothesis that, due to the low canopy height, measurements can be made in the shallow equilibrium layer where the wind field exhibits characteristics akin to level terrain. In order to test the validity of

We present a proof-of-concept for the adaptive mesh refinement method applied to atmospheric boundary-layer simulations. Such a method may form an attractive alternative to static grids for studies on atmospheric flows that have a high degree of scale separation in space and/or time. Examples include the diurnal cycle and a convective boundary layer capped by a strong inversion. For such cases, large-eddy

One-dimensional urban models embedded in mesoscale numerical models may place several grid points within the urban canopy. This requires an accurate parametrization for shear stresses (i.e. vertical momentum fluxes) including the dispersive stress and momentum sinks at these points. We used a case study with a packing density of 33% and checked rigorously the vertical variation of spatially-averaged

Wind-tunnel experiments were carried out on fully-rough boundary layers with large roughness ( δ / h ≈ 10 , where h is the height of the roughness elements and δ is the boundary-layer thickness). Twelve different surface conditions were created by using LEGO™ bricks of uniform height. Six cases are tested for a fixed plan solidity ( λ P ) with variations in frontal density ( λ F ), while the other

The simulation of horizontally homogeneous boundary layers that have characteristics of weakly and moderately stable atmospheric flow is investigated, where the well-established wind engineering practice of using 'flow generators' to provide a deep boundary layer is employed. Primary attention is given to the flow above the surface layer, in the absence of an overlying inversion, as assessed from first-

The Monin-Obukhov similarity theory-based wind speed and potential temperature profiles are inherently coupled to each other. We have developed hybrid approaches to disentangle them, and as a direct consequence, the estimation of Obukhov length (and associated turbulent fluxes) from either wind-speed or temperature measurements becomes an effortless task. Additionally, our approaches give rise to two