Boundary layer dynamics over London, UK, as observed using Doppler lidar during REPARTEE-II

Urban boundary layers (UBLs) can be highly complex due to the heterogeneous roughness and heating of the surface, particularly at night. Due to a general lack of observations, it is not clear whether canonical models of boundary layer mixing are appropriate in modelling air quality in urban areas. This paper reports Doppler lidar observations of turbulence profiles in the centre of London, UK, as part of the second REPARTEE campaign in autumn 2007. Lidar-measured standard deviation of vertical velocity averaged over 30 min intervals generally compared well with in situ sonic anemometer measurements at 190 m on the BT telecommunications Tower. During calm, nocturnal periods, the lidar underestimated turbulent mixing due mainly to limited sampling rate. Mixing height derived from the turbulence, and aerosol layer height from the backscatter profiles, showed similar diurnal cycles ranging from c. 300 to 800 m, increasing to c. 200 to 850 m under clear skies. The aerosol layer height was sometimes significantly different to the mixing height, particularly at night under clear skies. For convective and neutral cases, the scaled turbulence profiles resembled canonical results; this was less clear for the stable case. Lidar observations clearly showed enhanced mixing beneath stratocumulus clouds reaching down on occasion to approximately half daytime boundary layer depth. On one occasion the nocturnal turbulent structure was consistent with a nocturnal jet, suggesting a stable layer. Given the general agreement between observations and canonical turbulence profiles, mixing timescales were calculated for passive scalars released at street level to reach the BT Tower using existing models of turbulent mixing. It was estimated to take c. 10 min to diffuse up to 190 m, rising to between 20 and 50 min at night, depending on stability. Determination of mixing timescales is important when comparing to physico-chemical processes acting on pollutant species measured simultaneously at both the ground and at the BT Tower during the campaign. From the 3 week autumnal data-set there is evidence for occasional stable layers in central London, effectively decoupling surface emissions from air aloft.

Kinetics of layer hopping in a diblock copolymer lamellar phase

In the ordered state, symmetric diblock copolymers self-assemble into an anisotropic lamellar morphology. The equilibrium thickness of the lamellae is the result of a delicate balance between enthalpic and entropic energies, which can be tuned by controlling the temperature. Here we devise a simple yet powerful method of detecting tiny changes in the lamellar thickness using optical microscopy. From such measurements we characterize the enthalpic interaction as well as the kinetics of molecules as they hop from one layer to the next in order to adjust the lamellar thickness in response to a temperature jump. The resolution of the measurements facilitate a direct comparison to predictions from self-consistent field theory.

Multi-layer radial basis function networks: an extension to the radial basis function

This paper presents the initial research carried out into a new neural network called the multilayer radial basis function network (MRBF). The network extends the radial basis function (RBF) in a similar way to that in which the multilayer perceptron extends the perceptron. It is hoped that by connecting RBFs together in a layered fashion, an equivalent increase in ability can be gained, as is gained from using MLPs instead of single perceptrons. The results of a practical comparison between individual RBFs and MRBF's are also given.

Spatial and temporal effects of free-air CO2 enrichment (POPFACE) on leaf growth, cell expansion and cell production in a closed canopy of poplar

Leaf expansion in the fast-growing tree,Populus × euramericana was stimulated by elevated [CO2] in a closed-canopy forest plantation, exposed using a free air CO2 enrichment technique enabling long-term experimentation in field conditions. The effects of elevated [CO2] over time were characterized and related to the leaf plastochron index (LPI), and showed that leaf expansion was stimulated at very early (LPI, 0–3) and late (LPI, 6–8) stages in development. Early and late effects of elevated [CO2] were largely the result of increased cell expansion and increased cell production, respectively. Spatial effects of elevated [CO2] were also marked and increased final leaf size resulted from an effect on leaf area, but not leaf length, demonstrating changed leaf shape in response to [CO2]. Leaves exhibited a basipetal gradient of leaf development, investigated by defining seven interveinal areas, with growth ceasing first at the leaf tip. Interestingly, and in contrast to other reports, no spatial differences in epidermal cell size were apparent across the lamina, whereas a clear basipetal gradient in cell production rate was found. These data suggest that the rate and timing of cell production was more important in determining leaf shape, given the constant cell size across the leaf lamina. The effect of elevated [CO2] imposed on this developmental gradient suggested that leaf cell production continued longer in elevated [CO2] and that basal increases in cell production rate were also more important than altered cell expansion for increased final leaf size and altered leaf shape in elevated [CO2].

Long-term acclimation of leaf production, development, longevity and quality following 3 yr exposure to free-air CO2 enrichment during canopy closure in Populus

The effects of elevated CO2 on leaf development in three genotypes of Populus were investigated during canopy closure, following exposure to elevated CO2 over 3 yr using free-air enrichment.• Leaf quality was altered such that nitrogen concentration per unit d. wt (Nmass) declined on average by 22 and 13% for sun and shade leaves, respectively, in elevated CO2. There was little evidence that this was the result of ‘dilution’ following accumulation of nonstructural carbohydrates. Most likely, this was the result of increased leaf thickness. Specific leaf area declined in elevated CO2 on average by 29 and 5% for sun and shade leaves, respectively.• Autumnal senescence was delayed in elevated CO2 with a 10% increase in the number of days at which 50% leaf loss occurred in elevated as compared with ambient CO2.• These data suggest that changes in leaf quality may be predicted following long-term acclimation of fast-growing forest trees to elevated CO2, and that canopy longevity may increase, with important implications for forest productivity.

Advances in continuously profiling the thermodynamic state of the boundary layer: integration of measurements and methods

This paper describes advances in ground-based thermodynamic profiling of the lower troposphere through sensor synergy. The well-documented integrated profiling technique (IPT), which uses a microwave profiler, a cloud radar, and a ceilometer to simultaneously retrieve vertical profiles of temperature, humidity, and liquid water content (LWC) of nonprecipitating clouds, is further developed toward an enhanced performance in the boundary layer and lower troposphere. For a more accurate temperature profile, this is accomplished by including an elevation scanning measurement modus of the microwave profiler. Height-dependent RMS accuracies of temperature (humidity) ranging from 0.3 to 0.9 K (0.5–0.8 g m−3) in the boundary layer are derived from retrieval simulations and confirmed experimentally with measurements at distinct heights taken during the 2005 International Lindenberg Campaign for Assessment of Humidity and Cloud Profiling Systems and its Impact on High-Resolution Modeling (LAUNCH) of the German Weather Service. Temperature inversions, especially of the lower boundary layer, are captured in a very satisfactory way by using the elevation scanning mode. To improve the quality of liquid water content measurements in clouds the authors incorporate a sophisticated target classification scheme developed within the European cloud observing network CloudNet. It allows the detailed discrimination between different types of backscatterers detected by cloud radar and ceilometer. Finally, to allow IPT application also to drizzling cases, an LWC profiling method is integrated. This technique classifies the detected hydrometeors into three different size classes using certain thresholds determined by radar reflectivity and/or ceilometer extinction profiles. By inclusion into IPT, the retrieved profiles are made consistent with the measurements of the microwave profiler and an LWC a priori profile. Results of IPT application to 13 days of the LAUNCH campaign are analyzed, and the importance of integrated profiling for model evaluation is underlined.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Hydrogen-bonding-driven self-assembly of PEGylated organosilica nanoparticles with poly(acrylic acid) in aqueous solutions and in layer-by-layer deposition at solid surfaces

PEGylated organosilica nanoparticles have been synthesized through self-condensation of (3-mercaptopropyl)trimethoxysilane in dimethyl sulfoxide into thiolated nanoparticles with their subsequent reaction with methoxypoly(ethylene glycol) maleimide. The PEGylated nanoparticles showed excellent colloidal stability over a wide range of pH in contrast to the parent thiolated nanoparticles, which have a tendency to aggregate irreversibly under acidic conditions (pH < 3.0). Due to the presence of a poly(ethylene glycol)-based corona, the PEGylated nanoparticles are capable of forming hydrogen-bonded interpolymer complexes with poly(acrylic acid) in aqueous solutions under acidic conditions, resulting in larger aggregates. The use of hydrogen-bonding interactions allows more efficient attachment of the nanoparticles to surfaces. The alternating deposition of PEGylated nanoparticles and poly(acrylic acid) on silicon wafer surfaces in a layer-by-layer fashion leads to multilayered coatings. The self-assembly of PEGylated nanoparticles with poly(acrylic acid) in aqueous solutions and at solid surfaces was compared to the behavior of linear poly(ethylene glycol). The nanoparticle system creates thicker layers than the poly(ethylene glycol), and a thicker layer is obtained on a poly(acrylic acid) surface than on a silica surface, because of the effects of hydrogen bonding. Some implications of these hydrogen-bonding-driven interactions between PEGylated nanoparticles and poly(acrylic acid) for pharmaceutical formulations are discussed.

Determining the contribution of volcanic ash and boundary layer aerosol in backscatter lidar returns: a three‐component atmosphere approach

A solution of the lidar equation is discussed, that permits combining backscatter and depolarization measurements to quantitatively distinguish two different aerosol types with different depolarization properties. The method has been successfully applied to simultaneous observations of volcanic ash and boundary layer aerosol obtained in Exeter, United Kingdom, on 16 and 18 April 2010, permitting the contribution of the two aerosols to be quantified separately. First a subset of the atmospheric profiles is used where the two aerosol types belong to clearly distinguished layers, for the purpose of characterizing the ash in terms of lidar ratio and depolarization. These quantities are then used in a three‐component atmosphere solution scheme of the lidar equation applied to the full data set, in order to compute the optical properties of both aerosol types separately. On 16 April a thin ash layer, 100–400 m deep, is observed (average and maximum estimated ash optical depth: 0.11 and 0.2); it descends from ∼2800 to ∼1400 m altitude over a 6‐hour period. On 18 April a double ash layer, ∼400 m deep, is observed just above the morning boundary layer (average and maximum estimated ash optical depth: 0.19 and 0.27). In the afternoon the ash is entrained into the boundary layer, and the latter reaches a depth of ∼1800 m (average and maximum estimated ash optical depth: 0.1 and 0.15). An additional ash layer, with a very small optical depth, was observed on 18 April at an altitude of 3500–4000 m. By converting the lidar optical measurements using estimates of volcanic ash specific extinction, derived from other works, the observations seem to suggest approximate peak ash concentrations of ∼1500 and ∼1000 mg/m3,respectively, on the two observations dates.

Wind-direction effects on urban-type flows

Practically all extant work on flows over obstacle arrays, whether laboratory experiments or numerical modelling, is for cases where the oncoming wind is normal to salient faces of the obstacles. In the field, however, this is rarely the case. Here, simulations of flows at various directions over arrays of cubes representing typical urban canopy regions are presented and discussed. The computations are of both direct numerical simulation and large-eddy simulation type. Attention is concentrated on the differences in the mean flow within the canopy region arising from the different wind directions and the consequent effects on global properties such as the total surface drag, which can change very significantly—by up to a factor of three in some circumstances. It is shown that for a given Reynolds number the typical viscous forces are generally a rather larger fraction of the pressure forces (principally the drag) for non-normal than for normal wind directions and that, dependent on the surface morphology, the average flow direction deep within the canopy can be largely independent of the oncoming wind direction. Even for regular arrays of regular obstacles, a wind direction not normal to the obstacle faces can in general generate a lateral lift force (in the direction normal to the oncoming flow). The results demonstrate this and it is shown how computations in a finite domain with the oncoming flow generated by an appropriate forcing term (e.g. a pressure gradient) then lead inevitably to an oncoming wind direction aloft that is not aligned with the forcing term vector.

Wind-driven mixing below the oceanic mixed layer

This study describes the turbulent processes in the upper ocean boundary layer forced by a constant surface stress in the absence of the Coriolis force using large-eddy simulation. The boundary layer that develops has a two-layer structure, a well-mixed layer above a stratified shear layer. The depth of the mixed layer is approximately constant, whereas the depth of the shear layer increases with time. The turbulent momentum flux varies approximately linearly from the surface to the base of the shear layer. There is a maximum in the production of turbulence through shear at the base of the mixed layer. The magnitude of the shear production increases with time. The increase is mainly a result of the increase in the turbulent momentum flux at the base of the mixed layer due to the increase in the depth of the boundary layer. The length scale for the shear turbulence is the boundary layer depth. A simple scaling is proposed for the magnitude of the shear production that depends on the surface forcing and the average mixed layer current. The scaling can be interpreted in terms of the divergence of a mean kinetic energy flux. A simple bulk model of the boundary layer is developed to obtain equations describing the variation of the mixed layer and boundary layer depths with time. The model shows that the rate at which the boundary layer deepens does not depend on the stratification of the thermocline. The bulk model shows that the variation in the mixed layer depth is small as long as the surface buoyancy flux is small.

Dispersion of a point-source release of a passive scalar through an urban-like array for different wind directions

The dispersion of a point-source release of a passive scalar in a regular array of cubical, urban-like, obstacles is investigated by means of direct numerical simulations. The simulations are conducted under conditions of neutral stability and fully rough turbulent flow, at a roughness Reynolds number of Reτ = 500. The Navier–Stokes and scalar equations are integrated assuming a constant rate release from a point source close to the ground within the array. We focus on short-range dispersion, when most of the material is still within the building canopy. Mean and fluctuating concentrations are computed for three different pressure gradient directions (0◦ , 30◦ , 45◦). The results agree well with available experimental data measured in a water channel for a flow angle of 0◦ . Profiles of mean concentration and the three-dimensional structure of the dispersion pattern are compared for the different forcing angles. A number of processes affecting the plume structure are identified and discussed, including: (i) advection or channelling of scalar down ‘streets’, (ii) lateral dispersion by turbulent fluctuations and topological dispersion induced by dividing streamlines around buildings, (iii) skewing of the plume due to flow turning with height, (iv) detrainment by turbulent dispersion or mean recirculation, (v) entrainment and release of scalar in building wakes, giving rise to ‘secondary sources’, (vi) plume meandering due to unsteady turbulent fluctuations. Finally, results on relative concentration fluctuations are presented and compared with the literature for point source dispersion over flat terrain and urban arrays. Keywords Direct numerical simulation · Dispersion modelling · Urban array