Showing papers in "Boundary-Layer Meteorology in 2014"

Bridging the Transition from Mesoscale to Microscale Turbulence in Numerical Weather Prediction Models

Abstract: With a focus towards developing multiscale capabilities in numerical weather prediction models, the specific problem of the transition from the mesoscale to the microscale is investigated. For that purpose, idealized one-way nested mesoscale to large-eddy simulation (LES) experiments were carried out using the Weather Research and Forecasting model framework. It is demonstrated that switching from one-dimensional turbulent diffusion in the mesoscale model to three-dimensional LES mixing does not necessarily result in an instantaneous development of turbulence in the LES domain. On the contrary, very large fetches are needed for the natural transition to turbulence to occur. The computational burden imposed by these long fetches necessitates the development of methods to accelerate the generation of turbulence on a nested LES domain forced by a smooth mesoscale inflow. To that end, four new methods based upon finite amplitude perturbations of the potential temperature field along the LES inflow boundaries are developed, and investigated under convective conditions. Each method accelerated the development of turbulence within the LES domain, with two of the methods resulting in a rapid generation of production and inertial range energy content associated to microscales that is consistent with non-nested simulations using periodic boundary conditions. The cell perturbation approach, the simplest and most efficient of the best performing methods, was investigated further under neutral and stable conditions. Successful results were obtained in all the regimes, where satisfactory agreement of mean velocity, variances and turbulent fluxes, as well as velocity and temperature spectra, was achieved with reference non-nested simulations. In contrast, the non-perturbed LES solution exhibited important energy deficits associated to a delayed establishment of fully-developed turbulence. The cell perturbation method has negligible computational cost, significantly accelerates the generation of realistic turbulence, and requires minimal parameter tuning, with the necessary information relatable to mean inflow conditions provided by the mesoscale solution.

A Multi-layer Radiation Model for Urban Neighbourhoods with Trees

Abstract: A neighbourhood-scale multi-layer urban canopy model of shortwave and longwave radiation exchange that explicitly includes the radiative effects of tall vegetation (trees) is presented. Tree foliage is permitted both between and above buildings, and mutual shading, emission and reflection between buildings and trees are included. The basic geometry is a two-dimensional canyon with leaf area density profiles and probabilistic variation of building height. Furthermore, the model accounts for three-dimensional path lengths through the foliage. Ray tracing determines the receipt of direct shortwave irradiance by building and foliage elements. View factors for longwave and shortwave diffuse radiation exchange are computed once at the start of the simulation using a Monte Carlo ray tracing approach; for subsequent model timesteps, matrix inversion rapidly solves infinite reflections and interception of emitted longwave between all elements. The model is designed to simulate any combination of shortwave and longwave radiation frequency bands, and to be portable to any neighbourhood-scale urban canopy geometry based on the urban canyon. Additionally, the model is sufficiently flexible to represent forest and forest-clearing scenarios. Model sensitivity tests demonstrate the model is robust and computationally feasible, and highlight the importance of vertical resolution to the performance of urban canopy radiation models. Full model evaluation is limited by the paucity of within-canyon radiation measurements in urban neighbourhoods with trees. Where appropriate model components are tested against analytic relations and results from an independent urban radiation transfer model. Furthermore, system response tests demonstrate the ability of the model to realistically distribute shortwave radiation among urban elements as a function of built form, solar angle and tree foliage height, density and clumping. Separate modelling of photosynthetically-active and near-infrared shortwave bands is shown to be important in some cases. Increased canyon height-to-width ratio and/or tree cover diminishes the net longwave radiation loss of individual canyon elements (e.g., floor, walls), but, notably, has little effect on the net longwave loss of the whole urban canopy. When combined with parametrizations for the impacts of trees on airflow and hydrological processes in the urban surface layer, the new radiation model extends the applicability of urban canopy models and permits more robust assessment of trees as tools to manage urban climate, air quality, human comfort and building energy loads.

Evaluation of the Weather Research and Forecasting Mesoscale Model for GABLS3: Impact of Boundary-Layer Schemes, Boundary Conditions and Spin-Up

Abstract: We evaluated the performance of the three-dimensional Weather Research and Forecasting (WRF) mesoscale model, specifically the performance of the planetary boundary-layer (PBL) parametrizations. For this purpose, Cabauw tower observations were used, with the study extending beyond the third GEWEX Atmospheric Boundary-Layer Study (GABLS3) one-dimensional model intercomparison. The WRF model (version 3.4.1) contains 12 different PBL parametrizations, most of which have been only partially evaluated. The GABLS3 case offers a clear opportunity to evaluate model performance, focusing on time series of near-surface weather variables, radiation and surface flux budgets, vertical structure and the nighttime inertial oscillation. The model results revealed substantial differences between the PBL schemes. Generally, non-local schemes tend to produce higher temperatures and higher wind speeds than local schemes, in particular, for nighttime. The WRF model underestimates the 2-m temperature during daytime (about $$2$$ K) and substantially underestimates it at night (about $$4$$ K), in contrast to the previous studies where modelled 2-m temperature was overestimated. Considering the 10-m wind speed, during the night turbulent kinetic energy based schemes tend to produce lower wind speeds than other schemes. In all simulations the sensible and latent heat fluxes were well reproduced. For the net radiation and the soil heat flux we found good agreement with daytime observations but underestimations at night. Concerning the vertical profiles, the selected non-local PBL schemes underestimate the PBL depth and the low-level jet altitude at night by about 50 m, although with the correct wind speed. The latter contradicts most previous studies and can be attributed to the revised stability function in the Yonsei University PBL scheme. The local, turbulent kinetic energy based PBL schemes estimated the low-level jet altitude and strength more accurately. Compared to the observations, all model simulations show a similar structure for the potential temperature, with a consistent cold bias ( $$\approx $$ 2 K) in the upper PBL. In addition to the sensitivity to the PBL schemes, we studied the sensitivity to technical features such as horizontal resolution and domain size. We found a substantial difference in the model performance for a range of 12, 18 and 24 h spin-up times, longer spin-up time decreased the modelled wind speed bias, but it strengthened the negative temperature bias. The sensitivity of the model to the vertical resolution of the input and boundary conditions on the model performance is confirmed, and its influence appeared most significant for the non-local PBL parametrizations.

Extension of the Averaging Time in Eddy-Covariance Measurements and Its Effect on the Energy Balance Closure

Abstract: The modified ogive analysis and the block ensemble average were employed to investigate the impact of the averaging time extension on the energy balance closure over six land-use types. The modified ogive analysis, which requires a steady-state condition, can extend the averaging time up to a few hours and suggests that an averaging time of 30 min is still overall sufficient for eddy-covariance measurements over low vegetation. The block ensemble average, which does not require a steady-state condition, can extend the averaging time to several days. However, it can improve the energy balance closure for some sites during specific periods, when secondary circulations exist in the vicinity of the sensor. These near-surface secondary circulations mainly transport sensible heat, and when near-ground warm air is transported upward, the sensible heat flux observed by the block ensemble average will increase at longer averaging times. These findings suggest an alternative energy balance correction for a ground-based eddy-covariance measurement, in which the attribution of the residual depends on the ratio of sensible heat flux to the buoyancy flux. The fraction of the residual attributed to the sensible heat flux by this energy balance correction is larger than in the energy balance correction that preserves the Bowen ratio.

Global Intermittency and Collapsing Turbulence in the Stratified Planetary Boundary Layer

Abstract: Direct numerical simulation of the turbulent Ekman layer over a smooth wall is used to investigate bulk properties of a planetary boundary layer under stable stratification. Our simplified configuration depends on two non-dimensional parameters: a Richardson number characterizing the stratification and a Reynolds number characterizing the turbulence scale separation. This simplified configuration is sufficient to reproduce global intermittency, a turbulence collapse, and the decoupling of the surface from the outer region of the boundary layer. Global intermittency appears even in the absence of local perturbations at the surface; the only requirement is that large-scale structures several times wider than the boundary-layer height have enough space to develop. Analysis of the mean velocity, turbulence kinetic energy, and external intermittency is used to investigate the large-scale structures and corresponding differences between stably stratified Ekman flow and channel flow. Both configurations show a similar transition to the turbulence collapse, overshoot of turbulence kinetic energy, and spectral properties. Differences in the outer region resulting from the rotation of the system lead, however, to the generation of enstrophy in the non-turbulent patches of the Ekman flow. The coefficient of the stability correction function from Monin–Obukhov similarity theory is estimated as $$\beta \approx 5.7$$ in agreement with atmospheric observations, theoretical considerations, and results from stably stratified channel flows. Our results demonstrate the applicability of this set-up to atmospheric problems despite the intermediate Reynolds number achieved in our simulations.

The Third GABLS Intercomparison Case for Evaluation Studies of Boundary-Layer Models. Part B: Results and Process Understanding

Abstract: We describe and analyze the results of the third global energy and water cycle experiment atmospheric boundary layer Study intercomparison and evaluation study for single-column models. Each of the nineteen participating models was operated with its own physics package, including land-surface, radiation and turbulent mixing schemes, for a full diurnal cycle selected from the Cabauw observatory archive. By carefully prescribing the temporal evolution of the forcings on the vertical column, the models could be evaluated against observations. We focus on the gross features of the stable boundary layer (SBL), such as the onset of evening momentum decoupling, the 2-m minimum temperature, the evolution of the inertial oscillation and the morning transition. New process diagrams are introduced to interpret the variety of model results and the relative importance of processes in the SBL; the diagrams include the results of a number of sensitivity runs performed with one of the models. The models are characterized in terms of thermal coupling to the soil, longwave radiation and turbulent mixing. It is shown that differences in longwave radiation schemes among the models have only a small effect on the simulations; however, there are significant variations in downward radiation due to different boundary-layer profiles of temperature and humidity. The differences in modelled thermal coupling to the land surface are large and explain most of the variations in 2-m air temperature and longwave incoming radiation among models. Models with strong turbulent mixing overestimate the boundary-layer height, underestimate the wind speed at 200 m, and give a relatively large downward sensible heat flux. The result is that 2-m air temperature is relatively insensitive to turbulent mixing intensity. Evening transition times spread 1.5 h around the observed time of transition, with later transitions for models with coarse resolution. Time of onset in the morning transition spreads 2 h around the observed transition time. With this case, the morning transition appeared to be difficult to study, no relation could be found between the studied processes, and the variation in the time of the morning transition among the models.

Turbulent Transport of Momentum and Scalars Above an Urban Canopy

Abstract: Turbulent transport of momentum and scalars over an urban canopy is investigated using the quadrant analysis technique. High-frequency measurements are available at three levels above the urban canopy (47, 140 and 280 m). The characteristics of coherent ejection–sweep motions (flux contributions and time fractions) at the three levels are analyzed, particularly focusing on the difference between ejections and sweeps, the dissimilarity between momentum and scalars, and the dissimilarity between the different scalars (i.e., temperature, water vapour and $$\hbox {CO}_{2})$$ . It is found that ejections dominate momentum and scalar transfer at all three levels under unstable conditions, while sweeps are the dominant eddy motions for transporting momentum and scalars in the urban roughness sublayer under neutral and stable conditions. The flux contributions and time fractions of ejections and sweeps can be adequately captured by assuming a Gaussian joint probability density function for flow variables. However, the inequality of flux contributions from ejections and sweeps is more accurately reproduced by the third-order cumulant expansion method (CEM). The incomplete cumulant expansion method (ICEM) also works well except for $$\hbox {CO}_{2}$$ at 47 m where the skewness of $$\hbox {CO}_{2}$$ fluctuations is significantly larger than that for vertical velocity. The dissimilarity between momentum and scalar transfers is linked to the dissimilarity in the characteristics of ejection–sweep motions and is further quantified by measures of transport efficiencies. Atmospheric stability is the controlling factor for the transport efficiencies of momentum and heat, and fitted functions from the literature describe their behaviour fairly accurately. However, transport efficiencies of water vapour and $$\hbox {CO}_{2}$$ are less affected by the atmospheric stability. The dissimilarity among the three scalars examined in this study is linked to the active role of temperature and to the surface heterogeneity effect.

Determining Boundary-Layer Height from Aircraft Measurements

Abstract: The height of the atmospheric boundary layer (ABL) is an important variable in both observational studies and model simulations. The most commonly used measurement for obtaining ABL height is a rawinsonde profile. Mesoscale or regional scale models use a bulk Richardson number based on profiles of the forecast variables. Here we evaluate the limitations of several frequently-used approaches for defining ABL height from a single profile, and identify the optimal threshold value for each method if profiles are the only available measurements. Aircraft measurements from five field projects are used, representing a variety of ABL conditions including stable, convective, and cloud-topped boundary layers over different underlying surfaces. ABL heights detected from these methods were validated against the ‘true’ value determined from aircraft soundings, where ABL height is defined as the top of the layer with significant turbulence. A detection rate was defined to denote how often the ABL height was correctly diagnosed with a particular method. The results suggest that the temperature gradient method provides the most reasonable estimates, although the detection rate and suitable detection criteria vary for different types of ABL. The Richardson number method, on the other hand, is in most cases inadequate or inferior to the other methods that were tried. The optimal range of the detection criteria is given for all ABL types examined in this study.

Wind-Tunnel Simulation of the Wake of a Large Wind Turbine in a Stable Boundary Layer: Part 2, the Wake Flow

Abstract: Measurements have been made in the wake of a model wind turbine in both a neutral and a stable atmospheric boundary layer, in the EnFlo stratified-flow wind tunnel, between 05 and 10 rotor diameters from the turbine, as part of an investigation of wakes in offshore winds In the stable case the velocity deficit decreased more slowly than in the neutral case, partly because the boundary-layer turbulence levels are lower and the consequentially reduced level of mixing, an ‘indirect’ effect of stratification A correlation for velocity deficit showed the effect of stratification to be the same over the whole of the measured extent, following a polynomial form from about five diameters After about this distance (for the present stratification) the vertical growth of the wake became almost completely suppressed, though with an increased lateral growth; the wake in effect became ‘squashed’, with peaks of quantities occurring at a lower height, a ‘direct’ effect of stratification Generally, the Reynolds stresses were lower in magnitude, though the effect of stratification was larger in the streamwise fluctuation than on the vertical fluctuations The vertical heat flux did not change much from the undisturbed level in the first part of the wake, but became much larger in the later part, from about five diameters onwards, and exceeded the surface level at a point above hub height

Large-Eddy Simulation of Pollutant Removal from a Three-Dimensional Street Canyon

Abstract: Large-eddy simulations were conducted to investigate the mechanism of pollutant removal from a three-dimensional street canyon. Five block configurations with aspect ratios (building height to length) of 1, 2, 4, 8 and \(\infty \) were used to create an urban-like array. A pollutant was released from a ground-level line source at the centre of the target canyon floor. For smaller aspect ratios, the relative contribution of the turbulent mass flux to net mass flux at the roof level, which was spatially averaged along the roof-level ventilation area, was closer to unity, indicating that turbulent motions mainly affected pollutant removal from the top of the canyon. As aspect ratio increased, the relative contribution became smaller, owing to strong upwind motions. However, the relative contribution again reached near unity for the infinite aspect ratio (i.e. a two-dimensional street canyon) because of lowered lateral flow convergence. At least 75 % of total emissions from the three-dimensional street canyon were attributable to turbulent motions. Pollutant removal by turbulent motions was related to the coherent structures of low-momentum fluid above the canyons. Though the coherent structure size of the low-momentum fluid differed, the positions of low-momentum fluid largely corresponded to instantaneous high concentrations of pollutant above the target canyon, irrespective of canyon geometry.

The Third GABLS Intercomparison Case for Evaluation Studies of Boundary-Layer Models. Part A: Case Selection and Set-Up

Abstract: We describe a novel methodology on the selection and composition of a single-case observational dataset from the comprehensive measurement program at the Cabauw observatory field site located in the Netherlands. The case can be regarded as the basis of the third case study conducted within the framework of the GEWEX (Global Energy and Water Exchange) Atmospheric Boundary-Layer Study (GABLS) and is meant to be used for the evaluation of single-column models. The ideal case is supposed to cover a period of at least 24 h with clear skies, moderate near-surface winds and a stable stratification during nighttime. From the multi-year data archive with Cabauw observations data for 1–2 July 2006 were found to best match the requirements, and were consequently selected for analysis. The dates contains a 24-h period with a nearly constant geostrophic wind of $$\approx $$ 7 m $$\hbox {s}^{-1}$$ , and a considerable wind shear in the vertical. It is also characterized by a pronounced low-level jet related to an inertial oscillation that developed around sunset when the atmosphere had decoupled from the surface. Detailed initial conditions, surface conditions and dynamical forcings are derived on the basis of local observations and the outcome of a conceptual and a three-dimensional atmospheric model. It appears that a very precise prescription of the forcings is a prerequisite to enable a meaningful evaluation of models against observational data.

Impacts of Mixing Processes in Nocturnal Atmospheric Boundary Layer on Urban Ozone Concentrations

Abstract: A number of open questions remain regarding the role of low-level jets (LLJs) and nocturnal mixing processes in the buildup of tropospheric ozone. The prevalence of southerly winds and LLJs in the U.S. Southern Great Plains during summer makes this region an ideal site for investigating the structure of the nocturnal boundary layer and its impacts on urban air quality. Ozone \((\mathrm{O}_{3})\) and nitrogen oxide concentrations measured at regulatory monitoring sites in the Oklahoma City (OKC) area and simulations with the Weather Research and Forecasting with Chemistry (WRF/Chem) model were analyzed to show how the nocturnal LLJ moderates boundary-layer mixing processes and air quality. Datasets collected during the Joint Urban 2003 campaign, which took place in July 2003 in OKC, provided detailed information about nocturnal boundary-layer structure and dynamics. In general, \({\mathrm{O}_{3}}\) time series show the expected behavior that urban \({\mathrm{O}_{3}}\) concentrations decrease at night due to nitrogen oxide titration reactions, but elevated \({\mathrm{O}_{3}}\) concentrations and secondary \({\mathrm{O}_{3}}\) peaks are also seen quite frequently after sunset. LLJs developed on most nights during the study period and were associated with strong vertical wind shear, which affected the boundary-layer stability and structure. Near-surface \({\mathrm{O}_{3}}\) concentrations are higher during less stable nights when active mixing persists throughout the night. The WRF/Chem model results agree well with the observations and further demonstrate the role of LLJs in moderating nocturnal mixing processes and air quality. The highest nocturnal \({\mathrm{O}_{3}}\) concentrations are linked to a strong LLJ that promotes both nocturnal long-range transport and persistent downward mixing of \({\mathrm{O}_{3}}\) from the residual layer to the surface.

Moving Beyond Monin–Obukhov Similarity Theory in Modelling Wind-Speed Profiles in the Lower Atmospheric Boundary Layer under Stable Stratification

Abstract: Monin–Obukhov similarity theory (MOST) is commonly used to model the wind-speed profile at altitudes relevant to wind-power production (e.g. 10–200 m). Though reasonably accurate for unstable to weakly stable stratification, this approach becomes less accurate under increasingly stable stratification, largely due to the constant-flux surface layer assumed by MOST becoming shallower than the altitude range of interest. Furthermore, above the surface layer, the Coriolis force has a considerable influence on the wind-speed profile (in particular in the formation of low-level jets) that cannot be modelled using similarity theory. Our goal is to compare the accuracy of alternative extrapolation models that are more physically appropriate above the surface layer. Using data from the 213-m Cabauw meteorological tower in the Netherlands between July 2007 and June 2008, it is shown that MOST is accurate only at low altitudes and low stability, and breaks down at high altitudes and high stability. Local similarity is generally more accurate than MOST across all altitudes and stabilities, though the model requires turbulent flux data at multiple altitudes that is generally impractical. In contrast, a two-layer MOST–Ekman model is found to be comparable to the other models at low stability ranges and considerably more accurate in the high stability range, while requiring only a measure of surface stability and the geostrophic wind.

Modelling Small-Scale Drifting Snow with a Lagrangian Stochastic Model Based on Large-Eddy Simulations

Abstract: Observations of drifting snow on small scales have shown that, in spite of nearly steady winds, the snow mass flux can strongly fluctuate in time and space. Most drifting snow models, however, are not able to describe drifting snow accurately over short time periods or on small spatial scales as they rely on mean flow fields and assume equilibrium saltation. In an attempt to gain understanding of the temporal and spatial variability of drifting snow on small scales, we propose to use a model combination of flow fields from large-eddy simulations (LES) and a Lagrangian stochastic model to calculate snow particle trajectories and so infer snow mass fluxes. Model results show that, if particle aerodynamic entrainment is driven by the shear stress retrieved from the LES, we can obtain a snow mass flux varying in space and time. The obtained fluctuating snow mass flux is qualitatively compared to field and wind-tunnel measurements. The comparison shows that the model results capture the intermittent behaviour of observed drifting snow mass flux yet differences between modelled turbulent structures and those likely to be found in the field complicate quantitative comparisons. Results of a model experiment show that the surface shear-stress distribution and its influence on aerodynamic entrainment appear to be key factors in explaining the intermittency of drifting snow.

Differences in Atmospheric Boundary-Layer Characteristics Between Wet and Dry Episodes of the Indian Summer Monsoon

Abstract: The differences and similarities in atmospheric boundary-layer (ABL) characteristics, in particular the ABL height, evolution and wind field, between two contrasting episodes of the Indian summer monsoon have been studied using measurements from wind profilers and an instrumented 50-m tower at Gadanki in India. The observed differences are discussed in light of various forcing mechanisms, in particular the effect of soil moisture on the surface energy balance and ABL. The differences in ABL height, its evolution and the wind field between episodes are quite pronounced. Wet episodes not only have a shallower ABL but also the growth is delayed by 1–4 h when compared with that for dry episodes. Abundant soil moisture during the wet episodes (a factor of two greater than during the dry episodes) reduces the buoyancy flux, and thereby not only limits the ABL height but also delays the commencement of ABL growth. The low-level jet (LLJ) is stronger during the dry episodes and has a larger diurnal range than during the wet episodes. The highest occurrence and magnitude of LLJ apparent at a height of 1.5 km during early morning hours shift progressively with height and time till the afternoon, following ABL evolution. The weaker LLJ during the wet episodes is attributed to its southward migration from its mean position (15 $$^{\circ }$$ N). Larger signal-to-noise ratio and spectral width values are observed during the early night to midnight, compared to noon-time, when the ABL is buoyantly turbulent.

Very-Large-Scale Motions in the Atmospheric Boundary Layer Educed by Snapshot Proper Orthogonal Decomposition

Abstract: Large-eddy simulations of the atmospheric boundary layer (ABL) under a wide range of stabilities are conducted to educe very-large-scale motions and then to study their dynamics and how they are influenced by buoyancy. Preliminary flow visualizations suggest that smaller-scale motions that resemble hairpins are embedded in much larger scale streamwise meandering rolls. Using simulations that represent more than 150 h of physical time, many snapshots in the $$xy$$ -, $$yz$$ - and $$xz$$ -planes are then collected to perform snapshot proper orthogonal decomposition and further investigate the large structures. These analyses confirm that large streamwise rolls that share several features with the very-large-scale motions observed in laboratory studies arise as the dominant modes under most stabilities, but the effect of the surface kinematic buoyancy flux on the energy content of these dominant modes is very significant. The first two modes in the $$yz$$ -plane in the neutral case contain up to 3 % of the total turbulent kinetic energy; they also have a vertical tilt angle in the $$yz$$ -plane of about 0 to 30 $$^\circ $$ due to the turning effect associated with the Coriolis force. Unstable cases also feature streamwise rolls, but in the convective ABL they are strengthened by rising plumes in between them, with two to four rolls spanning the whole domain in the first few modes; the Coriolis effect is much weaker in the unstable ABL. These rolls are no longer the dominant modes under stable conditions where the first mode is observed to contain sheet-like motions with high turbulent kinetic energy. Using these proper orthogonal decomposition modes, we are also able to extract the vertical velocity fields corresponding to individual modes and then to correlate them with the horizontal velocity or temperature fields to obtain the momentum and heat flux carried by individual modes. Structurally, the fluxes are explained by the topology of their corresponding modes. However, the fraction of the fluxes produced by the modes is invariably smaller than the fraction of energy they contain, particularly under stable conditions where the first modes are found to perform weak counter-gradient fluxes.

Multi-Scale Sensible Heat Fluxes in the Suburban Environment from Large-Aperture Scintillometry and Eddy Covariance

Abstract: Sensible heat fluxes ( $$Q_{H}$$ ) are determined using scintillometry and eddy covariance over a suburban area. Two large-aperture scintillometers provide spatially integrated fluxes across path lengths of 2.8 and 5.5 km over Swindon, UK. The shorter scintillometer path spans newly built residential areas and has an approximate source area of 2–4 $$\text {km}^{2}$$ , whilst the long path extends from the rural outskirts to the town centre and has a source area of around 5–10 $$\text {km}^{2}$$ . These large-scale heat fluxes are compared with local-scale eddy-covariance measurements. Clear seasonal trends are revealed by the long duration of this dataset and variability in monthly $$Q_{H}$$ is related to the meteorological conditions. At shorter time scales the response of $$Q_{H}$$ to solar radiation often gives rise to close agreement between the measurements, but during times of rapidly changing cloud cover spatial differences in the net radiation ( $$Q^{*}$$ ) coincide with greater differences between heat fluxes. For clear days $$Q_{H}$$ lags $$Q^{*}$$ , thus the ratio of $$Q_{H}$$ to $$Q^{*}$$ increases throughout the day. In summer the observed energy partitioning is related to the vegetation fraction through use of a footprint model. The results demonstrate the value of scintillometry for integrating surface heterogeneity and offer improved understanding of the influence of anthropogenic materials on surface-atmosphere interactions.

Nesting Turbulence in an Offshore Convective Boundary Layer Using Large-Eddy Simulations

Abstract: The applicability of the one-way nesting technique for numerical simulations of the heterogeneous atmospheric boundary layer using the large-eddy simulation (LES) framework of the Weather Research and Forecasting model is investigated. The focus of this study is on LES of offshore convective boundary layers. Simulations were carried out using two subgrid-scale models (linear and non-linear) with two different closures [diagnostic and prognostic subgrid-scale turbulent kinetic energy (TKE) equations]. We found that the non-linear backscatter and anisotropy model with a prognostic subgrid-scale TKE equation is capable of providing similar results when performing one-way nested LES to a stand-alone domain having the same grid resolution but using periodic lateral boundary conditions. A good agreement is obtained in terms of velocity shear and turbulent fluxes, while velocity variances are overestimated. A streamwise fetch of 14 km is needed following each domain transition in order for the solution to reach quasi-stationary results and for the velocity spectra to generate proper energy content at high wavelengths, however, a pile-up of energy is observed at the low-wavelength portion of the spectrum on the first nested domain. The inclusion of a second nest with higher resolution allows the solution to reach effective grid spacing well within the Kolmogorov inertial subrange of turbulence and develop an appropriate energy cascade that eliminates most of the pile-up of energy at low wavelengths. Consequently, the overestimation of velocity variances is substantially reduced and a considerably better agreement with respect to the stand-alone domain results is achieved.

Large-Eddy Simulation of Mesoscale Circulations Forced by Inhomogeneous Urban Heat Island

Abstract: The large-eddy simulation mode of the Weather Research and Forecasting model is employed to simulate the planetary boundary-layer characteristics and mesoscale circulations forced by an ideal urban heat island (UHI). In our simulations, the horizontal heterogeneity of the UHI intensity distribution in urban areas is considered and idealized as a cosine function. Results indicate that the UHI heating rate and the UHI intensity heterogeneity affect directly the spatial distribution of the wind field; a stronger UHI intensity produces a maximum horizontal wind speed closer to the urban centre. The strong advection of warm air from the urban area to the rural area in the upper part of the planetary boundary-layer causes a more stable atmospheric stratification over both the urban and rural areas. The mesoscale sensible heat flux caused by the UHI circulation increases with UHI intensity but vanishes when the background wind speed is sufficiently high $$(>$$ 3.0 $$\mathrm{{m\,s}}^{-1})$$ .

An Improved Three-Dimensional Simulation of the Diurnally Varying Street-Canyon Flow

Abstract: The impact of diurnal variations of the heat fluxes from building and ground surfaces on the fluid flow and air temperature distribution in street canyons is numerically investigated using the PArallelized Large-eddy Simulation Model (PALM). Simulations are performed for a 3 by 5 array of buildings with canyon aspect ratio of one for two clear summer days that differ in atmospheric instability. A detailed building energy model with a three-dimensional raster-type geometry—Temperature of Urban Facets Indoor-Outdoor Building Energy Simulator (TUF-IOBES)—provides urban surface heat fluxes as thermal boundary conditions for PALM. In vertical cross-sections at the centre of the spanwise canyon the mechanical forcing and the horizontal streamwise thermal forcing at roof level outweigh the thermal forces from the heated surfaces inside the canyon in defining the general flow pattern throughout the day. This results in a dominant canyon vortex with a persistent speed, centered at a constant height. Compared to neutral simulations, non-uniform heating of the urban canyon surfaces significantly modifies the pressure field and turbulence statistics in street canyons. Strong horizontal pressure gradients were detected in streamwise and spanwise canyons throughout the day, and which motivate larger turbulent velocity fluctuations in the horizontal directions rather than in the vertical direction. Canyon-averaged turbulent kinetic energy in all non-neutral simulations exhibits a diurnal cycle following the insolation on the ground in both spanwise and streamwise canyons, and it is larger when the canopy bottom surface is paved with darker materials and the ground surface temperature is higher as a result. Compared to uniformly distributed thermal forcing on urban surfaces, the present analysis shows that realistic non-uniform thermal forcing can result in complex local airflow patterns, as evident, for example, from the location of the vortices in horizontal planes in the spanwise canyon. This study shows the importance of three-dimensional simulations with detailed thermal boundary conditions to explore the heat and mass transport in an urban area.

Evaluation of Two Energy Balance Closure Parametrizations

Abstract: A general lack of energy balance closure indicates that tower-based eddy-covariance (EC) measurements underestimate turbulent heat fluxes, which calls for robust correction schemes. Two parametrization approaches that can be found in the literature were tested using data from the Canadian Twin Otter research aircraft and from tower-based measurements of the German Terrestrial Environmental Observatories (TERENO) programme. Our analysis shows that the approach of Huang et al. (Boundary-Layer Meteorol 127:273–292, 2008), based on large-eddy simulation, is not applicable to typical near-surface flux measurements because it was developed for heights above the surface layer and over homogeneous terrain. The biggest shortcoming of this parametrization is that the grid resolution of the model was too coarse so that the surface layer, where EC measurements are usually made, is not properly resolved. The empirical approach of Panin and Bernhofer (Izvestiya Atmos Oceanic Phys 44:701–716, 2008) considers landscape-level roughness heterogeneities that induce secondary circulations and at least gives a qualitative estimate of the energy balance closure. However, it does not consider any feature of landscape-scale heterogeneity other than surface roughness, such as surface temperature, surface moisture or topography. The failures of both approaches might indicate that the influence of mesoscale structures is not a sufficient explanation for the energy balance closure problem. However, our analysis of different wind-direction sectors shows that the upwind landscape-scale heterogeneity indeed influences the energy balance closure determined from tower flux data. We also analyzed the aircraft measurements with respect to the partitioning of the “missing energy” between sensible and latent heat fluxes and we could confirm the assumption of scalar similarity only for Bowen ratios $$\approx $$ 1.

Advances in Air–Sea \hbox {CO}_2 Flux Measurement by Eddy Correlation

Abstract: Eddy-correlation measurements of the oceanic $$\hbox {CO}_2$$ flux are useful for the development and validation of air–sea gas exchange models and for analysis of the marine carbon cycle. Results from more than a decade of published work and from two recent field programs illustrate the principal interferences from water vapour and motion, demonstrating experimental approaches for improving measurement precision and accuracy. Water vapour cross-sensitivity is the greatest source of error for $$\hbox {CO}_2$$ flux measurements using infrared gas analyzers, often leading to a ten-fold bias in the measured $$\hbox {CO}_2$$ flux. Much of this error is not related to optical contamination, as previously supposed. While various correction schemes have been demonstrated, the use of an air dryer and closed-path analyzer is the most effective way to eliminate this interference. This approach also obviates density corrections described by Webb et al. (Q J R Meteorol 106:85–100, 1980). Signal lag and frequency response are a concern with closed-path systems, but periodic gas pulses at the inlet tip provide for precise determination of lag time and frequency attenuation. Flux attenuation corrections are shown to be $$<$$ 5 % for a cavity ring-down analyzer (CRDS) and dryer with a 60-m inlet line. The estimated flux detection limit for the CRDS analyzer and dryer is a factor of ten better than for IRGAs sampling moist air. While ship-motion interference is apparent with all analyzers tested in this study, decorrelation or regression methods are effective in removing most of this bias from IRGA measurements and may also be applicable to the CRDS.

One Year of Surface-Based Temperature Inversions at Dome C, Antarctica

Abstract: In 2005 the Study of Stable Boundary Layer Environment at Dome C (STABLEDC) experimental campaign was conducted at the plateau station of Concordia at Dome C, Antarctica. Temperature profiles measured with a microwave radiometer were used to study the characteristics of surface-based temperature inversions over the course of a year. Statistics of temperature profiles for every month are discussed; the difference between daytime and nocturnal cases observed during the summer months disappears during winter. Surface-based temperature inversions occurred in 70 % of the time during summer, and almost all of the time during winter. During winter the occurrence of warming events leads to a decrease in the temperature difference between the top and the base of the inversion (i.e. the inversion strength). The inversion strength maxima ranged between \(3\,^{\circ }\mathrm{C}\) (December) and \(35\,^{\circ }\mathrm{C}\) (August) corresponding to gradients of 0.1 and \(0.3\,^{\circ }\mathrm{C}\, \mathrm{m}^{-1}\) , respectively. The average surface-based inversion height presents a daily cycle during the summer months with values up to 200 m in the morning hours, while it affects a layer always deeper than 100 m during the winter months. The relationships between inversion strength and the downward longwave radiative flux, absolute temperature, and wind speed are examined. The inversion strength decreases as the longwave radiation increases. A clear anti-correlation between inversion strength and near-surface temperature is evident throughout the year. During the winter, the largest inversion strength values were observed under low wind-speed conditions; in contrast, a clear dependence was not found during the summer.

Thermal stratification effects on flow over a generic urban canopy

Abstract: The influence of local surface heating and cooling on flow over urban-like roughness is investigated using large-eddy simulations. By adjusting the incoming or outgoing heat flux from the ground surface, various degrees of local thermal stratification, represented by a Richardson number $$(Ri_\tau )$$ , were attained. Drag and heat transfer coefficients, turbulence structure, integral length scales, and the strength of quadrant events that contribute to momentum and heat fluxes were obtained and are compared with locally stable, neutral and unstable flows. With increasing $$Ri_\tau $$ , or equivalently as the flow characteristics change from local thermal instability to stability, a gradual decline in the drag and heat transfer coefficients is observed. These values are found to be fairly independent of the type of thermal boundary condition (constant heat flux or constant temperature) and domain size. The maps of anisotropy invariants showed that for the values of $$Ri_\tau $$ considered, turbulence structures are almost the same in shape for neutral and unstable cases but differ slightly from those in the stable case. The degree of anisotropy is found to decrease as $$Ri_\tau $$ increases from $$-2$$ to 2.5. Compared to the neutral case, the integral length scales are shortened in the streamwise and vertical direction by ground cooling, but enhanced in the vertical direction with ground heating. Quadrant analysis showed that an increase in floor heating increases the strength of ejections above the canopy. However, the contributions of updrafts or downdrafts to the heat flux are found not to be significantly influenced by the type of local thermal stratification for the values of $$Ri_\tau $$ considered. From the octant analysis, the transport mechanisms of momentum and heat above the canopy are found to be very similar in both locally unstable and stable flows.

Predictability of the Meteorological Conditions Favourable to Radiative Fog Formation During the 2011 ParisFog Campaign

Abstract: Radiative fog formation is a complex phenomenon involving local physical and microphysical processes that take place when particular meteorological conditions occur. This study aims at quantifying the ability of a regional numerical weather model to analyze and forecast the conditions favourable to radiative fog formation at an instrumental site in the Paris area. Data from the ParisFog campaign have been used in order to quantify the meteorological conditions favorable to radiative fog formation (pre-fog conditions) by setting threshold values on the key meteorological variables driving this process: 2-m temperature tendency, 10-m wind speed, 2-m relative humidity and net infrared flux. Data from the ParisFog observation periods of November 2011 indicate that use of these thresholds leads to the detection of 87 % of cases in which radiative fog formation was observed. In order to evaluate the ability of a regional weather model to reproduce adequately these conditions, the same thresholds are applied to meteorological model fields in both analysis and forecast mode. It is shown that, with this simple methodology, the model detects 74 % of the meteorological conditions finally leading to observed radiative fog, and 48 % 2 days in advance. Finally, sensitivity tests are conducted in order to evaluate the impact of using larger time or space windows on the forecasting skills.

Detection of Entrainment Influences on Surface-Layer Measurements and Extension of Monin–Obukhov Similarity Theory

Abstract: We present a method to detect influences of boundary-layer processes on surface-layer measurements, using statistics and spectra of surface-layer variables only. We validated our detection method with boundary-layer measurements. Furthermore, we confirm that Monin–Obukhov similarity functions fit well to temperature-variance data obtained at two different homogeneous surfaces. However, we found that humidity variance measurements deviate from the universal functions above one of the two studied surfaces for days on which entrained air reached the surface layer. These results confirm that Monin–Obukhov similarity theory should be used with care in the analysis of surface-layer data. Finally, we propose the use of an extra term in flux-variance relations that depends on the entrainment ratio for humidity and on the boundary-layer height. If boundary-layer measurements are not available, we show how the entrainment ratio for humidity can be approximated from the skewness of the humidity distribution.

A Length Scale Defining Partially-Resolved Boundary-Layer Turbulence Simulations

Abstract: Numerical weather prediction (NWP) model forecasts at horizontal grid lengths in the range of 100 m to 1 km are now possible. Within this range of grid lengths, the convective boundary layer (CBL) is partially resolved and thus in the so-called ‘grey zone’. For simulations in the grey zone, numerical dissipation sources from both the advection scheme and the subgrid model are likely to be significant. Until now, these effects have not been incorporated fully into our understanding of the grey zone. In order to quantify these effects, a dissipation length scale is defined based on the second moment of the turbulent kinetic energy (TKE) spectrum. An ensemble of simulations of a CBL are performed using a large-eddy model across the grey-zone resolutions and for a range of subgrid model, advection scheme and vertical grid configurations. The dissipation length scale distinguishes the effects of the different model configurations in the grey zone. In the middle of the boundary layer, the resolved TKE is strongly controlled by the numerical dissipation. This leads to a similarity law for the resolved TKE in the grey zone using the dissipation length scale. A new definition of the grey zone emerges where the inversion depth and dissipation length scale are the same size. This contrasts with the typical definition using the horizontal grid length. At the inversion, however, the variation of the dissipation length scale with grid length is less predictable, reflecting significant challenges for modelling entrainment in the grey zone. The dissipation length scale is thus a simple diagnostic to aid both NWP and large-eddy modellers in understanding the grey zone.

Investigation of the Planetary Boundary Layer in the Swiss Alps Using Remote Sensing and In Situ Measurements

Abstract: The development of the planetary boundary layer (PBL) has been studied in a complex terrain using various remote sensing and in situ techniques. The high-altitude research station at Jungfraujoch (3,580 m a.s.l.) in the Swiss Alps lies for most of the time in the free troposphere except when it is influenced by the PBL reaching the station, especially during the summer season. A ceilometer and a wind profiler were installed at Kleine Scheidegg, a mountain pass close to Jungfraujoch, located at an altitude of 2,061 m a.s.l. Data from the ceilometer were analyzed using two different algorithms, while the signal-to-noise ratio of the wind profiler was studied to compare the retrieved PBL heights. The retrieved values from the ceilometer and wind profiler agreed well during daytime and cloud-free conditions. The results were additionally compared with the PBL height estimated by the numerical weather prediction model COSMO-2, which showed a clear underestimation of the PBL height for most of the cases but occasionally also a slight overestimation especially around noon, when the PBL showed its maximum extent. Air parcels were transported upwards by slope winds towards Jungfraujoch when the PBL was higher than 2,800 m a.s.l. during cloud-free cases. This was confirmed by the in situ aerosol measurements at Jungfraujoch with a significant increase in particle number concentration, particle light absorption and scattering coefficients when PBL-influenced air masses reached the station in the afternoon hours. The continuous aerosol in situ measurements at Jungfraujoch were clearly influenced by the local PBL development but also by long-range transport phenomena such as Saharan dust or pollution from the south.

Intercomparison of Methods for the Simultaneous Estimation of Zero-Plane Displacement and Aerodynamic Roughness Length from Single-Level Eddy-Covariance Data

Abstract: We applied three approaches to estimate the zero-plane displacement $$d$$ through the aerodynamic measurement height $$z$$ (with $$z = z_{m}- d$$ and $$z_{m}$$ being the measurement height above the surface), and the aerodynamic roughness length $$z_{0}$$ , from single-level eddy covariance data. Two approaches (one iterative and one regression-based) were based on the universal function in the logarithmic wind profile and yielded an inherently simultaneous estimation of both $$d$$ and $$z_{0}$$ . The third approach was based on flux–variance similarity, where estimation of $$d$$ and consecutive estimation of $$z_{0}$$ are independent steps. Each approach was further divided into two methods differing either with respect to the solution technique (profile approaches) or with respect to the variable (variance of vertical wind and temperature, respectively). All methods were applied to measurements above a large, growing wheat field where a uniform canopy height and its frequent monitoring provided plausibility limits for the resulting estimates of time-variant $$d$$ and $$z_{0}$$ . After applying, for each approach, a specific data filtering that accounted for the range of conditions (e.g. stability) for which it is valid, five of the six methods were able to describe the temporal changes of roughness parameters associated with crop growth and harvest, and four of them agreed on $$d$$ to within 0.3 m most of the time. Application of the same methods to measurements with a more heterogeneous footprint consisting of fully-grown sugarbeet and a varying contribution of adjacent harvested fields exhibited a plausible dependence of the roughness parameters on the sugarbeet fraction. It also revealed that the methods producing the largest outliers can differ between site conditions and stability. We therefore conclude that when determining $$d$$ for canopies with unknown properties from single-level measurements, as is increasingly done, it is important to compare the results of a number of methods rather than rely on a single one. An ensemble average or median of the results, possibly after elimination of methods that produce outliers, can help to yield more robust estimates. The estimates of $$z_{0}$$ were almost exclusively physically plausible, although $$d$$ was considered unknown and estimated simultaneously with the methods and results described above.

Influence from Surrounding Land on the Turbulence Measurements Above a Lake

Abstract: Turbulence measurements taken at a Swedish lake are analyzed. Although the measurements took place over a relatively large lake with several km of undisturbed fetch, the turbulence structure was found to be highly influenced by the surrounding land during daytime. Variance spectra of both horizontal velocity and scalars during both unstable and stable stratification displayed a low frequency peak. The energy at lower frequencies showed a daily variation, increasing in the morning and decreasing in the afternoon. This behaviour is explained by spectral lag, where the low frequency energy due to large eddies that originate from the convective boundary layer above the surrounding land. When the air is advected over the lake the small eddies rapidly equilibrate with the new surface forcing. However, the large eddies remain for an appreciable distance and influence the turbulence in the developing lake boundary layer. The variances of the horizontal velocity and scalars are increased by these large eddies, while the turbulent fluxes are mainly unaffected. The drag coefficient, Stanton number and Dalton number used to parametrize the momentum flux, heat flux and latent heat flux respectively all compare well with current parametrizations developed for open sea conditions. The diurnal cycle of the partial pressure of methane, $$p\mathrm{CH}_{4}$$ , observed at this site is closely related to the diurnal cycle of the lake-air methane flux. An idealized two-dimensional model simulation of the boundary layer at a lake site indicates that the strong response of $$p\mathrm{CH}_{4}$$ to the surface methane flux is due to the shallow internal boundary layer that develops above the lake, allowing methane to accumulate in a relatively small volume.