On the Relevance of $$\ln (z_{0} /z_{0T})= kB^{-1}$$

The assumption that the roughness Reynolds number \(( Re_{*})\) can be used as a basis for quantifying the boundary-layer property \({ kB}^{-1} (= \ln (z_{0}/z_{0T}))\) as in some modern numerical models is questioned. While \({ Re}_{*}\) is a useful property in studies of pipe flow, it appears to have only marginal applicability in the case of treeless terrain, as studied in the two experimental situations presented here. For both the daytime and night-time cases there appears to be little correlation between \({ kB}^{-1}\) and \({ Re}_{*}\). For daytime, the present studies indicate that the assumption \({ kB}^{-1} \approx 2\) is acceptable, while for night-time, the scatter involved in relating \({ kB}^{-1}\) to \({ Re}_{*}\) suggests there is little reason to assume a direct relationship. However, while the scatter affecting all of the night-time results is large, there remains a significant correlation between the heat and momentum fluxes upon which an alternative methodology for describing bulk air–surface exchange at night could be constructed. The friction coefficient (\(C_{f}\)) and the turbulent Stanton number \(({ St}_{*})\) are discussed as possible alternatives for describing bulk properties of the air layer adjacent to the surface. While describing the surface roughness in terms of the friction coefficient provides an attractive simplification relative to the conventional methodologies based on roughness length and stability considerations, use of the Stanton number shares many of uncertainties that affect \({ kB}^{-1}\). The transitions at dawn and dusk remain demanding situations to address.     

Water-Channel Estimation of Eulerian and Lagrangian Time Scales of the Turbulence in Idealized Two-Dimensional Urban Canopies

Lagrangian and Eulerian statistics are obtained from a water-channel experiment of an idealized two-dimensional urban canopy flow in neutral conditions. The objective is to quantify the Eulerian \((T^{\mathrm{E}})\) and Lagrangian \((T^{\mathrm{L}})\) time scales of the turbulence above the canopy layer as well as to investigate their dependence on the aspect ratio of the canopy, AR, as the latter is the ratio of the width (W) to the height (H) of the canyon. Experiments are also conducted for the case of flat terrain, which can be thought of as equivalent to a classical one-directional shear flow. The values found for the Eulerian time scales on flat terrain are in agreement with previous numerical results found in the literature. It is found that both the streamwise and vertical components of the Lagrangian time scale, \(T_\mathrm{u}^\mathrm{L} \) and \(T_\mathrm{w}^\mathrm{L} \), follow Raupach’s linear law within the constant-flux layer. The same holds true for \(T_\mathrm{w}^\mathrm{L} \) in both the canopies analyzed \((AR= 1\) and \(AR= 2\)) and also for \(T_\mathrm{u}^\mathrm{L} \) when \(AR = 1\). In contrast, for \(AR = 2\), \(T_\mathrm{u}^\mathrm{L} \) follows Raupach’s law only above \(z=2H\). Below that level, \(T_\mathrm{u}^\mathrm{L} \) is nearly constant with height, showing at \(z=H\) a value approximately one order of magnitude greater than that found for \(AR = 1\). It is shown that the assumption usually adopted for flat terrain, that \(\beta =T^{\mathrm{L}}/T^{\mathrm{E}}\) is proportional to the inverse of the turbulence intensity, also holds true even for the canopy flow in the constant-flux layer. In particular, \(\gamma /i_\mathrm{u} \) fits well \(\beta _\mathrm{u} =T_\mathrm{u}^\mathrm{L} /T_\mathrm{u}^\mathrm{E} \) in both the configurations by choosing \(\gamma \) to be 0.35 (here, \(i_\mathrm{u} =\sigma _\mathrm{u} / \bar{u} \), where \(\bar{u} \) and \(\sigma _\mathrm{u} \) are the mean and the root-mean-square of the streamwise velocity component, respectively). On the other hand, \(\beta _\mathrm{w} =T_\mathrm{w}^\mathrm{L} /T_\mathrm{w}^\mathrm{E} \) follows approximately \(\gamma /i_\mathrm{w} =0.65/\left( {\sigma _\mathrm{w} /\bar{u} } \right) \) for \(z > 2H\), irrespective of the AR value. The second main objective is to estimate other parameters of interest in dispersion studies, such as the eddy diffusivity of momentum \((K_\mathrm{{T}})\) and the Kolmogorov constant \((C_0)\). It is found that \(C_0\) depends appreciably on the velocity component both for the flat terrain and canopy flow, even though for the latter case it is insensitive to AR values. In all the three experimental configurations analyzed here, \(K_\mathrm{{T}}\) shows an overall linear growth with height in agreement with the linear trend predicted by Prandtl’s theory.     

Test of Scintillometer Saturation Correction Methods Using Field Experimental Data

Saturation of large aperture scintillometer (LAS) signals can result in sensible heat flux measurements that are biased low. A field study with LASs of different aperture sizes and path lengths was performed to investigate the onset of, and corrections for, signal saturation. Saturation already occurs at \({C_n^2 \approx 0.074 D^{5/3} \lambda^{1/3} L^{-8/3}}\), where \({C_n^2}\) is the structure parameter of the refractive index, D is the aperture size, λ is the wavelength, L is the transect length, which is smaller than theoretically derived saturation limits. At a transect length of 1 km, a height of 2.5 m, and aperture ≈0.15 m the correction factor exceeds 5% already at \({C_n^2=2\times 10^{-12}{\rm m}^{-2/3}}\), which will affect many practical applications of scintillometry. The Clifford correction method, which only depends on \({C_n^2}\) and the transect geometry, provides good saturation corrections over the range of conditions observed in our study. The saturation correction proposed by Ochs and Hill results in correction factors that are too small in large saturation regimes. An inner length scale dependence of the saturation correction factor was not observed. Thus for practical applications the Clifford correction method should be applied.     

Seasonal Variations in Drag Coefficient over a Sastrugi-Covered Snowfield in Coastal East Antarctica

The surface of windy Antarctic snowfields is subject to drifting snow, which leads to the formation of sastrugi. In turn, sastrugi contribute to the drag exerted by the snow surface on the atmosphere and hence influence drifting snow. Although the surface drag over rough sastrugi fields has been estimated for individual locations in Antarctica, its variation over time and with respect to drifting snow has received little attention. Using year-round data from a meteorological mast, seasonal variations in the neutral drag coefficient at a height of 10 m \((C_{{ DN}10})\) in coastal Adelie Land are presented and discussed in light of the formation and behaviour of sastrugi based on observed aeolian erosion patterns. The measurements revealed high \(C_{{ DN}10} \) values \((\ge \) 2 \(\times \) 10\(^{-3})\) and limited drifting snow (35% of the time) in summer (December–February) versus lower \(C_{{ DN}10} \) values \((\approx \) 1.5 \(\times \) \(10^{-3})\) associated with more frequent drifting snow (70% of the time) in winter (March–November). Without the seasonal distinction, there was no clear dependence of \(C_{{ DN}10} \) on friction velocity or wind direction, but observations revealed a general increase in \(C_{{ DN}10} \) with rising air temperature. The main hypothesis defended here is that higher temperatures increase snow cohesion and the development of sastrugi just after snow deposition while inhibiting the sastrugi streamlining process by raising the erosion threshold. This increases the contribution of the sastrugi form drag to the total surface drag in summer when winds are lighter and more variable. The analysis also showed that, in the absence of erosion, single snowfall events can reduce \(C_{{ DN}10} \) to \(1\,\times \,10^{-3}\) due to the burying of pre-existing microrelief under newly deposited snow. The results suggest that polar atmospheric models should account for spatial and temporal variations in snow surface roughness through a dynamic representation of the sastrugi form drag.     

Estimating Sediment Mass Fluxes on Surfaces Sheltered by Live Vegetation

We present a simple model based on already existing and widely used equations for estimating particle mass fluxes on surfaces sheltered by live vegetation. Wind-tunnel measurements of vertical profiles of mass flux in three different dense live plant canopies, and as a function of the spatially averaged skin friction velocity \({u_{\tau }}'\), provide the baseline set of data. For the bare-sand surface, the total mass flux Q shows the typical \(b({u_\tau }' - {u_{\tau t}}')^{3 }\) increase with increasing skin friction velocity \({u_{\tau }}'\), where b is a constant and \({u_{\tau t}}'\) is the threshold at the onset of particle erosion. Similar relations, however, with different values for b and \({u_{\tau t}}'\) compared to the bare-sand surface were found for experiments with 5.25 and 24.5 plants \(\hbox {m}^{-2}\) and can be explained by the spatial variations of \(u_{\tau }\) for the canopy cases. Based on the resulting parameters b and \({u_{\tau t}}'\), which are found to be functions of the roughness density \(\lambda \), we present a final simple relation \(Q(\lambda ,\, {u_{\tau }}')\) used for estimating the total mass flux for surfaces sheltered by live vegetation.     

Katabatic Flow: A Closed-Form Solution with Spatially-Varying Eddy Diffusivities

The Nieuwstadt closed-form solution for the stationary Ekman layer is generalized for katabatic flows within the conceptual framework of the Prandtl model. The proposed solution is valid for spatially-varying eddy viscosity and diffusivity (O’Brien type) and constant Prandtl number (Pr). Variations in the velocity and buoyancy profiles are discussed as a function of the dimensionless model parameters \(z_0 \equiv \hat{z}_0 \hat{N}^2 Pr \sin {(\alpha )} |\hat{b}_\mathrm{s} |^{-1}\) and \(\lambda \equiv \hat{u}_{\mathrm{ref}}\hat{N} \sqrt{Pr} |\hat{b}_\mathrm{s} |^{-1}\), where \(\hat{z}_0\) is the hydrodynamic roughness length, \(\hat{N}\) is the Brunt-Väisälä frequency, \(\alpha \) is the surface sloping angle, \(\hat{b}_\mathrm{s}\) is the imposed surface buoyancy, and \(\hat{u}_{\mathrm{ref}}\) is a reference velocity scale used to define eddy diffusivities. Velocity and buoyancy profiles show significant variations in both phase and amplitude of extrema with respect to the classic constant \(\textit{K}\) model and with respect to a recent approximate analytic solution based on the Wentzel-Kramers-Brillouin theory. Near-wall regions are characterized by relatively stronger surface momentum and buoyancy gradients, whose magnitude is proportional to \(z_0\) and to \(\lambda \). In addition, slope-parallel momentum and buoyancy fluxes are reduced, the low-level jet is further displaced toward the wall, and its peak velocity depends on both \(z_0\) and \(\lambda \).     

Edge-to-Stem Variability in Wet-Canopy Evaporation From an Urban Tree Row

Evaporation from wet-canopy (\(E_\mathrm{C}\)) and stem (\(E_\mathrm{S}\)) surfaces during rainfall represents a significant portion of municipal-to-global scale hydrologic cycles. For urban ecosystems, \(E_\mathrm{C}\) and \(E_\mathrm{S}\) dynamics play valuable roles in stormwater management. Despite this, canopy-interception loss studies typically ignore crown-scale variability in \(E_\mathrm{C}\) and assume (with few indirect data) that \(E_\mathrm{S}\) is generally \({<}2\%\) of total wet-canopy evaporation. We test these common assumptions for the first time with a spatially-distributed network of in-canopy meteorological monitoring and 45 surface temperature sensors in an urban Pinus elliottii tree row to estimate \(E_\mathrm{C}\) and \(E_\mathrm{S}\) under the assumption that crown surfaces behave as “wet bulbs”. From December 2015 through July 2016, 33 saturated crown periods (195 h of 5-min observations) were isolated from storms for determination of 5-min evaporation rates ranging from negligible to 0.67 \(\hbox {mm h}^{-1}\). Mean \(E_\mathrm{S}\) (0.10 \(\hbox {mm h}^{-1}\)) was significantly lower (\(p < 0.01\)) than mean \(E_\mathrm{C}\) (0.16 \(\hbox {mm h}^{-1}\)). But, \(E_\mathrm{S}\) values often equalled \(E_\mathrm{C}\) and, when scaled to trunk area using terrestrial lidar, accounted for 8–13% (inter-quartile range) of total wet-crown evaporation (\(E_\mathrm{S}+E_\mathrm{C}\) scaled to surface area). \(E_\mathrm{S}\) contributions to total wet-crown evaporation maximized at 33%, showing a general underestimate (by 2–17 times) of this quantity in the literature. Moreover, results suggest wet-crown evaporation from urban tree rows can be adequately estimated by simply assuming saturated tree surfaces behave as wet bulbs, avoiding problematic assumptions associated with other physically-based methods.     

Velocity Adjustment and Passive Scalar Diffusion in and Above an Urban Canopy in Response to Various Approach Flows

We used wind-tunnel experiments to investigate velocity-field adjustment and scalar diffusion behaviour in and above urban canopies located downwind of various roughness elements. Staggered arrays of rectangular blocks of various heights H and plan area ratios λ_p were used to model the urban canopies. The velocity field in the roughness sublayer (height \({z \lesssim 2H}\)) reached equilibrium at distances proportional to \({\sqrt{L_{\rm c}H}}\) where L _c is the canopy-drag length scale determined as a function of λ_p and the block side length L. A distance of about \({20\sqrt{L_{\rm c}H}}\) was required for adjustment at z = H/2 (in the canopy), and a distance of about \({10\sqrt{L_{\rm c}H}}\) was required at z = 2H (near the top of the roughness sublayer). Diffusion experiments from a ground emission source revealed that differences in upwind roughness conditions had negligible effects on the plume growth near the source (up to a few multiples of L from the source) if the source was located at a fetch F larger than about \({10\sqrt{L_{\rm c}H}}\) from the upwind edge of the canopy. However, at locations farther downwind (more than several multiples of L from the source), upwind conditions had considerable effects on the plume growth. For a representative urban canopy, it was shown that a much larger fetch than required for velocity-field adjustment in the roughness sublayer was necessary to eliminate the effects of upwind conditions on plume widths at 24L downwind from the source.     

Observations of the Temperature and Humidity Structure Parameter Over Heterogeneous Terrain by Airborne Measurements During the LITFASS-2003 Campaign

We used numerical simulations to investigate the general relationship between urban morphology and the intensity of wind gusts in built-up areas at the pedestrian level. The simulated urban boundary layer developed over a 19.2 km (length) \(\times \) 4.8 km (width) \(\times \) 1.0 km (height) simulation domain, with 2-m resolution in all directions, to explicitly resolve the detailed shapes of buildings and the flow at the pedestrian level. This complex computation was accomplished using the lattice Boltzmann method and by implementing a large-eddy simulation model. To generalize the results, a new parameter that expresses the intensity of gusts (the gust index, \({\tilde{U}}_{ max})\) was defined as the local maximum wind speed divided by the freestream velocity. In addition, this parameter was decomposed into the mean wind-speed ratio, \({\tilde{U}} \) and turbulent gust ratio, \({\tilde{U}}^{{\prime }}\) to evaluate the qualities of gusts. These parameters were useful for quantitatively comparing the gust intensities within urban canopies at different locations or even among different experiments. In addition, the entire horizontal domain was subdivided into homogeneous square patches, in which both the simulated gust parameters and the morphological characteristics of building geometries were averaged. This procedure masked the detailed structure of individual buildings but retained the bulk characteristics of the urban morphology. At the pedestrian level, the gust index decreased with increasing building cover. Compared to \({\tilde{U}} \), the quantity \({\tilde{U}}^{{\prime }}\) notably contributed to the index throughout the range of plan area index \((\lambda _p)\) values. The dependences of all normalized wind-speed ratios transiently changed at \(\lambda _p =~0.28\). In cases where \(\lambda _p < 0.28, {\tilde{U}} \) decreased with increasing \(\lambda _p \), although \({\tilde{U}}^{{\prime }}\) was almost constant. In cases where \(\lambda _p > 0.28, {\tilde{U}}\) was almost constant and \({\tilde{U}}^{{\prime }}\) decreased with increasing \(\lambda _p \). This was explained by the change in flow regimes within the building canyon. At a higher elevation above the canopy layer, \(\lambda _p \) becomes less relevant to normalized wind-speed ratios, and instead the aerodynamic roughness length became important.     

Observations of Atmospheric Methane and Carbon Dioxide Mixing Ratios: Tall-Tower or Mountain-Top Stations?

Mountain-top observations of greenhouse gas mixing ratios may be an alternative to tall-tower measurements for regional scale source and sink estimation. To investigate the equivalence or limitations of a mountain-top site as compared to a tall-tower site, we used the unique opportunity of comparing in situ measurements of methane (\(\hbox {CH}_{4}\)) and carbon dioxide (\(\hbox {CO}_{2}\)) mixing ratios at a mountain top (986 m above sea level, a.s.l.) with measurements from a nearby (distance 28.4 km) tall tower, sampled at almost the same elevation (1009 m a.s.l.). Special attention was given to, (i) how local wind statistics and greenhouse gas sources and sinks at the mountain top influence the observations, and (ii) whether mountain-top observations can be used as for those from a tall tower for constraining regional greenhouse gas emissions. Wind statistics at the mountain-top site are clearly more influenced by local flow systems than those at the tall-tower site. Differences in temporal patterns of the greenhouse gas mixing ratios observed at the two sites are mostly related to the influence of local sources and sinks at the mountain-top site. Major influences of local sources can be removed by applying a statistical filter (\(5{\mathrm{th}}\) percentile) or a filter that removes periods with unfavourable flow conditions. In the best case, the bias in mixing ratios between the mountain-top and the tall-tower sites after the application of the wind filter was \({-}0.0005\pm 0.0010\) ppm for methane (September, 0000–0400 UTC) and \(0.11\pm 0.18\) ppm for \(\hbox {CO}_{2}\) (February, 1200–1600 UTC). Temporal fluctuations of atmospheric \(\hbox {CH}_{4}\) and \(\hbox {CO}_{2}\) mixing ratios at both stations also showed good agreement (apart from \(\hbox {CO}_{2}\) during summertime) as determined by moving bi-weekly Pearson correlation coefficients (up to 0.96 for \(\hbox {CO}_{2}\) and 0.97 for \(\hbox {CH}_{4}\)). When only comparing mixing ratios minimally influenced by local sources (low bias and high correlation coefficients), our measurements indicate that mountain-top observations are comparable to tall-tower observations.     

Understanding Diurnality and Inter-Seasonality of a Sub-tropical Urban Heat Island

We quantify the spatial and temporal aspects of the urban heat-island (UHI) effect for Kanpur, a major city in the humid sub-tropical monsoon climate of the Gangetic basin. Fixed station measurements are used to investigate the diurnality and inter-seasonality in the urban–rural differences in surface temperature (\({\Delta } T_\mathrm{s}\)) and air temperature (\({\Delta } T_\mathrm{c}\)) separately. The extent of the spatial variations of the nighttime \({\Delta } T_\mathrm{c}\) and \({\Delta } T_\mathrm{s}\) is investigated through mobile campaigns and satellite remote sensing respectively. Nighttime \({\Delta } T_\mathrm{c}\) values dominate during both the pre-monsoon (maximum of 3.6 \(^\circ \hbox {C}\)) and the monsoon (maximum of 2.0 \(^\circ \hbox {C}\)). However, the diurnality in \({\Delta } T_\mathrm{s}\) is different, with higher daytime values during the pre-monsoon, but very little diurnality during the monsoon. The nighttime \({\Delta } T_\mathrm{s}\) value is mainly associated with differences in the urban–rural incoming longwave radiative flux (\(r^{2}=0.33\) during the pre-monsoon; 0.65 during the monsoon), which, in turn, causes a difference in the outgoing longwave radiative flux. This difference may modulate the nighttime \({\Delta } T_\mathrm{c}\) value as suggested by significant correlations (\(r^{2}=0.68\) for the pre-monsoon; 0.50 for the monsoon). The magnitude of \({\Delta } T_\mathrm{c}\) may also be modulated by advection, as it is inversely related with the urban wind speed. A combination of in situ, remotely sensed, and model simulation data were used to show that the inter-seasonality in \({\Delta } T_\mathrm{s}\), and, to a lesser extent, in \({\Delta } T_\mathrm{c}\), may be related to the change in the land use of the rural site between the pre-monsoon and the monsoon periods. Results suggest that the degree of coupling of \({\Delta } T_\mathrm{s}\) and \({\Delta } T_\mathrm{c}\) may be a strong function of land use and land cover.     

Quantifying the Influence of Random Errors in Turbulence Measurements on Scalar Similarity in the Atmospheric Surface Layer

The influence of random errors in turbulence measurements on scalar similarity for temperature, water vapour, \(\hbox {CO}_{2}\), and \(\hbox {NH}_{3}\) is investigated using two eddy-covariance datasets collected over a lake and a cattle feedlot. Three measures of scalar similarity, namely, the similarity constant in the flux–variance relationship, the correlation coefficient between two scalars and the relative transport efficiency, are examined. The uncertainty in the similarity constant \(C_{s}\) in the flux–variance relationship resulting from random errors in turbulence measurements is quantified based on error propagation analyses and a Monte-Carlo sampling method, which yields a distribution instead of a single value for \(C_{s}\). For different scalars, the distributions of \(C_{s}\) are found to significantly overlap, implying that scalars are transported similarly under strongly unstable conditions. The random errors in the correlation coefficients between scalars and the relative transport efficiencies are also quantified through error propagation analyses, and they increase as the atmosphere departs from neutral conditions. Furthermore, the correlation coefficients between three scalars (water vapour, \(\hbox {CO}_{2}\), and \(\hbox {NH}_{3}\)) are statistically different from unity while the relative transport efficiencies are not, which highlights the difference between these two measures of scalar similarity. The results suggest that uncertainties in these measures of scalar similarity need to be quantified when using them to diagnose the existence of dissimilarity among different scalars.     

The Effect of Breaking Waves on $$\hbox {CO}_2$$ Air–Sea Fluxes in the Coastal Zone

The influence of wave-associated parameters controlling turbulent \(\hbox {CO}_2\) fluxes through the air–sea interface is investigated in a coastal region. A full year of high-quality data of direct estimates of air–sea \(\hbox {CO}_2\) fluxes based on eddy-covariance measurements is presented. The study area located in Todos Santos Bay, Baja California, Mexico, is a net sink of \(\hbox {CO}_2\) with a mean flux of \(-1.3\, \upmu \hbox {mol m}^{-2}\hbox {s}^{-1}\) (\(-41.6\hbox { mol m}^{-2}\hbox {yr}^{-1}\)). The results of a quantile-regression analysis computed between the \(\hbox {CO}_2\) flux and, (1) wind speed, (2) significant wave height, (3) wave steepness, and (4) water temperature, suggest that the significant wave height is the most correlated parameter with the magnitude of the flux but the behaviour of the relation varies along the probability distribution function, with the slopes of the regression lines presenting both positive and negative values. These results imply that the presence of surface waves in coastal areas is the key factor that promotes the increase of the flux from and into the ocean. Further analysis suggests that the local characteristics of the aqueous and atmospheric layers might determine the direction of the flux.     

Towards a Validation of Scintillometer Measurements: The LITFASS-2009 Experiment

Scintillometry has been increasingly used over the last decade for the experimental determination of area-averaged turbulent fluxes at a horizontal scale of a few kilometres. Nevertheless, a number of assumptions in the scintillometer data processing and interpretation still call for a thorough evaluation, in particular over heterogeneous terrain. Moreover, a validation of the path-averaged structure parameters derived from scintillometer data (and forming the basis for the flux calculations) by independent measurements is still missing. To achieve this, the LITFASS-2009 field campaign has been performed around the Meteorological Observatory Lindenberg ?C Richard-A?mann-Observatory of the German Meteorological Service (DWD) in July 2009. The experiment combined tower-based in-situ turbulence measurements, field-scale laser scintillometers, long-range optical (large-aperture) and microwave scintillometers, and airborne turbulence measurements using an automatically operating unmanned aircraft. The paper describes the project design and strategy, and discusses first results. Daytime near-surface values of the temperature structure parameter, ${C_{T}^{2}}$ , over different types of farmland differ by more than one order of magnitude in their dependence on the type and status of the vegetation. Considerable spatial variability in ${C_{T}^{2}}$ was also found along the flight legs at heights between 50 and 100?m. However, it appeared difficult to separate the effects of heterogeneity from the temporal variability of the turbulence fields. Aircraft measurements and scintillometer data agreed in magnitude with respect to the temporal variation of the path-averaged ${C_{T}^{2}}$ values during the diurnal cycle. The decrease of ${C_{T}^{2}}$ with height found from the scintillometer measurements close to the surface and at 43?m under daytime convective conditions corresponds to free-convection scaling, whereas the aircraft measurements at 54 and 83?m suggest a different behaviour.     

Turbulence in Continental Stratocumulus,Part I: External Forcings and Turbulence Structures

Comprehensive, ground-based observations from the US Department of Energy Atmospheric Radiation Measurements program Southern Great Plains site are used to study the variability of turbulence forcings and cloud-scale turbulence structures in a continental stratocumulus cloud. The turbulence observations are made from an upward facing cloud (35 GHz) Doppler radar. Cloud base and liquid water path are characterized using a lidar at the surface and a microwave radiometer. The turbulence characterizations are compared and contrasted with those observed in marine stratocumulus clouds. During the 16-h observation period used in this study the cloud-base and cloud-top heights evolve with time and changes in liquid water path observed by the radiometer are consistent with variations in cloud depth. Unlike marine stratocumulus clouds, a diurnal cycle of cloud thickness and liquid water path is not observed. The observed surface latent, sensible, and virtual sensible heat fluxes and the radiative fluxes exhibit a diurnal cycle with values increasing from sunrise to afternoon and decreasing afterwards. During the night, the sensible heat, virtual sensible heat and the net radiative fluxes at the surface are slightly negative. Solar radiative heating prevails in the cloud layer during the day and strong radiative cooling exists at cloud top even during the day. Unlike marine stratocumulus, surface heating described by the convective velocity scale \(W_\mathrm{s}^{*}\) and cloud-top cooling described by \(W_\mathrm{r}^{*}\) are both important in driving the in-cloud turbulence during the day, whereas cloud-top cooling is the exclusive contributor during the night. The combined \(W_\mathrm{s}^{*}\) and \(W_\mathrm{r}^{*}\) (the total velocity scale \(W_\mathrm{t}^{*})\) provides a useful way to track the evolution of the turbulence structure in the cloud. The variance of the radar-measured radial velocity, which is related to resolved turbulence, follows the diurnal cycle and is consistent with the total velocity scale \(W_\mathrm{t}^{*}\) variations. It is higher during the day and lower during the night, which is contrary to that in marine stratocumulus. The \(W_\mathrm{t}^{*}\) values are lowest around sunset when the radiative cooling is also small due to upper-level clouds observed above the low-level stratus. The vertical distribution of the variance results from the surface heating during the day and cloud-top cooling during the night. The squared spectrum width, which is related to turbulence structures within the radar sampling volume (unresolved turbulence) also follows the diurnal cycle. Its vertical distribution indicates that the unresolved turbulence more closely relates to the processes near cloud top. Turbulence in the cloud requires about an hour to respond to the external forcings of surface heating and cloud-top radiative cooling. Positive skewness prevails during the day and negative skewness prevails at night with a sharp transition around sunset. Resolved turbulence dominates near cloud base whereas unresolved turbulence dominates near cloud top. The turbulence characteristics and variability defined in this study can be used to evaluate the time evolution of turbulence structures in large eddy simulation forced by surface and cloud-top radiative forcings.     

Observations of Lake-Breeze Events During the Toronto 2015 Pan-American Games

Enhanced meteorological observations were made during the 2015 Pan and Parapan American Games in Toronto in order to measure the vertical and horizontal structure of lake-breeze events. Two scanning Doppler lidars (one fixed and one mobile), a C-band radar, and a network including 53 surface meteorological stations (mesonet) provided pressure, temperature, humidity, and wind speed and direction measurements over Lake Ontario and urban areas. These observations captured the full evolution (prior, during, and after) of 27 lake-breeze events (73% of observation days) in order to characterize the convective and dynamic processes driving lake breezes at the local scale and mesoscale. The dominant signal of a passing lake-breeze front (LBF) was an increase in dew-point temperature of \(2.3 \pm 0.3 \,^{\circ }\hbox {C}\), coinciding with a \(180^{\circ }\) shift in wind direction and a decrease in air temperature of \(2.1 \pm 0.2 \,^{\circ }\hbox {C}\). Doppler lidar observations over the lake detected lake breezes 1 hour (on average) before detection by radar and mesonet. On days with the synoptic flow in the offshore direction, the lidars observed wedge-shaped LBFs with shallow depths, which inhibited the radar’s ability to detect the lake breeze. The LBF’s ground speed and inland penetration distance were found to be well-correlated (\(r = 0.78\)), with larger inland penetration distances occurring on days with non-opposing (non-offshore) synoptic flow. The observed enhanced vertical motion \(({>} 1\hbox { m s}^{-1})\) at the LBF, observed by the lidar on 54% of lake-breeze days, was greater (at times \({>} 2.5\hbox { m s}^{-1}\)) than that observed in previous studies and longer-lasting over the lake than over land. The weaker and less pronounced lake-breeze structure over land is illustrated in two case studies highlighting the lifetime of the lake-breeze circulation and the impact of propagation distance on lake-breeze intensity.     

Experimental Investigation of Turbulent Momentum Transfer in a Neutral Boundary Layer over a Rough Surface

The turbulent characteristics of the neutral boundary layer developing over rough surfaces are not well predicted with operational weather-forecasting models. The problem is attributed to inadequate mixing-length models, to the anisotropy of the flow and to a lack of controlled experimental data against which to validate numerical studies. Therefore, in order to address directly the modelling difficulties for the development of a neutral boundary layer over rough surfaces, and to investigate the turbulent momentum transfer of such a layer, a set of hydraulic flume experiments were carried out. In the experiments, the mean and turbulent quantities were measured by a particle image velocimetry (PIV) technique. The measured velocity variances and fluxes \({(\overline{{u_{i}^{\prime}}{u_{j}^{\prime}}})}\) in longitudinal vertical planes allowed the vertical and longitudinal gradients (?/?z and ?/?x) of the mean and turbulent quantities (fluxes, variances and third-order moments) to be evaluated and the terms of the evolution equations for ?e/?t, \({\partial \overline{u^{\prime 2}}/\partial t}\), \({\partial \overline{w^{\prime 2}}/\partial t}\) and \({\partial \overline{{u^{\prime}}{w^{\prime}}}/\partial t}\) to be quantified, where e is the turbulent kinetic energy. The results show that the pressure-correlation terms allow the turbulent energy to be transferred equitably from \({\overline{{u^{\prime}}^{2}}}\) to \({\overline{{w^{\prime}}^{2}}}\). It appears that the repartition between the constitutive terms of the budget of e, \({\overline{{u^{\prime}}^{2}}}\), \({\overline{{w^{\prime}}^{2}}}\) and \({\overline{{u^{\prime}}{w^{\prime}}}}\) is not significantly affected by the development of the rough neutral boundary layer. For the whole evolution, the transfers of energy are governed by the same terms that are also very similar to the smooth-wall case. The PIV measurements also allowed the spatial integral scales to be computed directly and to be compared with the dissipative and mixing length scales, which were also computed from the data.     

Investigation of Turbulence Parametrization Schemes with Reference to the Atmospheric Boundary Layer Over the Aegean Sea During Etesian Winds

The spatial structure of the marine atmospheric boundary layer (MABL) over the Aegean Sea is investigated using the Weather Research and Forecasting (WRF) mesoscale model. Two ‘first-order’ non-local and five ‘1.5-order’ local planetary boundary-layer (PBL) parametrization schemes are used. The predictions from the WRF model are evaluated against airborne observations obtained by the UK Facility for Airborne Atmospheric Measurements BAe-14 research aircraft during the Aegean-GAME field campaign. Statistical analysis shows good agreement between measurements and simulations especially at low altitude. Despite the differences between the predicted and measured wind speeds, they reach an agreement index of 0.76. The simulated wind-speed fields close to the surface differ substantially among the schemes (maximum values range from 13 to \(18\hbox { m s}^{-1}\) at 150-m height), but the differences become marginal at higher levels. In contrast, all schemes show similar spatial variation patterns in potential temperature fields. A warmer (1–2 K) and drier (2–3\(\hbox { g kg}^{-1})\) layer than is observed, is predicted by almost all schemes under stable conditions (eastern Aegean Sea), whereas a cooler (up to 2 K) and moister (1–2\(\hbox { g kg}^{-1})\) layer is simulated under near-neutral to nearly unstable conditions (western Aegean Sea). Almost all schemes reproduce the vertical structure of the PBL and the shallow MABL (up to 300 m) well, including the low-level jet in the eastern Aegean Sea, with non-local schemes being closer to observations. The simulated PBL depths diverge (up to 500 m) due to the different criteria applied by the schemes for their calculation. Under stable conditions, the observed MABL depth corresponds to the height above the sea surface where the simulated eddy viscosity reaches a minimum; under neutral to slightly unstable conditions this is close to the top of the simulated entrainment layer. The observed sensible heat fluxes vary from ?40 to \(25\hbox { W m}^{-2}\), while the simulated fluxes range from ?40 to \(40\hbox { W m}^{-2}\); however, all of the schemes’ predictions are close to the observations under unstable conditions. Finally, all schemes overestimate the friction velocity, although the simulated range (from 0.2 to \(0.5\hbox { m s}^{-1})\) is narrower than that observed (from 0.1 to \(0.7\hbox { m s}^{-1})\).     

High-Resolution Vertical Profile Measurements for Carbon Dioxide and Water Vapour Concentrations Within and Above Crop Canopies

We present a portable elevator-based facility for measuring \(\hbox {CO}_{2}\), water vapour, temperature and wind-speed profiles between the soil surface and the atmospheric surface layer above crop canopies. The end of a tube connected to a closed-path gas analyzer is continuously moved up and down over the profile range (in our case, approximately 2 m) while concentrations are logged at a frequency of \(20 \hbox { s}^{-1}\). Using campaign measurements in winter wheat, winter barley and a catch crop mixture (spring 2015 to autumn 2016) during different stages of crop development and different times of the day, we demonstrate a simple approach to correct for time lags, and the resulting profiles of 30-min mean mole fractions of \(\hbox {CO}_{2}\) and \(\hbox {H}_{2}\hbox {O}\) over height increments of 0.025 m. The profiles clearly show the effects of soil respiration and photosynthetic carbon assimilation, varying both during the diurnal cycle and during the growing season. Profiles of temperature and wind speed are based on a ventilated finewire thermocouple and a hot-wire anemometer, respectively. Measurements over bare soil and a short plant canopy were analyzed in the framework of Monin–Obukhov similarity theory to check the validity of the measurements and raw-data-processing approach. Derived fluxes of \(\hbox {CO}_{2}\), latent and sensible heat and momentum show good agreement with eddy-covariance measurements.     

The Effects of Heat Advection on UK Weather and Climate Observations in the Vicinity of Small Urbanized Areas

Weather and climate networks traditionally follow rigorous siting guidelines, with individual stations located away from frost hollows, trees or urban areas. However, the diverse nature of the UK landscape suggests that the feasibility of siting stations that are truly representative of regional climate and free from distorting local effects is increasingly difficult. Whilst the urban heat island is a well-studied phenomenon and usually accounted for, the effect of warm urban air advected downwind is rarely considered, particularly at rural stations adjacent to urban areas. Until recently, urban heat advection (UHA) was viewed as an urban boundary-layer process through the formation of an urban plume that rises above the surface as it is advected. However, these dynamic UHA effects are shown to also have an impact on surface observations. Results show a significant difference in temperatures anomalies (\(p\,< \,0.001\)) between observations taken downwind of urban and rural areas. For example, urban heat advection from small urbanized areas (\(\sim \)1\(\,\hbox {km}^{2}\)) under low cloud cover and wind speeds of 2–3\(\,\hbox {m}\,\hbox {s}^{-1}\) is found to increase mean nocturnal air temperatures by 0.6\(\,^{\circ }\hbox {C}\) at a horizontal distance of 0.5 km. Fundamentally, these UHA results highlight the importance of careful interpretation of long-term temperature data taken near small urban areas.