Trade-offs and responsiveness of the single-layer urban canopy parametrization in WRF: an offline evaluation using the MOSCEM optimization algorithm and field observations

For an increasing number of applications, mesoscale modelling systems now aim to better represent urban areas. The complexity of processes resolved by urban parametrization schemes varies with the application. The concept of fitness-for-purpose is therefore critical for both the choice of parametrizations and the way in which the scheme should be evaluated. A systematic and objective model response analysis procedure (Multiobjective Shuffled Complex Evolution Metropolis (MOSCEM) algorithm) is used to assess the fitness of the single-layer urban canopy parametrization implemented in the Weather Research and Forecasting (WRF) model. The scheme is evaluated regarding its ability to simulate observed surface energy fluxes and the sensitivity to input parameters. Recent amendments are described, focussing on features which improve its applicability to numerical weather prediction, such as a reduced and physically more meaningful list of input parameters. The study shows a high sensitivity of the scheme to parameters characterizing roof properties in contrast to a low response to road-related ones. Problems in partitioning of energy between turbulent sensible and latent heat fluxes are also emphasized. Some initial guidelines to prioritize efforts to obtain urban land-cover class characteristics in WRF are provided. Copyright © 2010 Royal Meteorological Society and Crown Copyright.

Developing a reliable strategy to infer the effective soil hydraulic properties from field evaporation experiments for agro-hydrological models

The Richards equation has been widely used for simulating soil water movement. However, the take-up of agro-hydrological models using the basic theory of soil water flow for optimizing irrigation, fertilizer and pesticide practices is still low. This is partly due to the difficulties in obtaining accurate values for soil hydraulic properties at a field scale. Here, we use an inverse technique to deduce the effective soil hydraulic properties, based on measuring the changes in the distribution of soil water with depth in a fallow field over a long period, subject to natural rainfall and evaporation using a robust micro Genetic Algorithm. A new optimized function was constructed from the soil water contents at different depths, and the soil water at field capacity. The deduced soil water retention curve was approximately parallel but higher than that derived from published pedo-tranfer functions for a given soil pressure head. The water contents calculated from the deduced soil hydraulic properties were in good agreement with the measured values. The reliability of the deduced soil hydraulic properties was tested in reproducing data measured from an independent experiment on the same soil cropped with leek. The calculation of root water uptake took account for both soil water potential and root density distribution. Results show that the predictions of soil water contents at various depths agree fairly well with the measurements, indicating that the inverse analysis is an effective and reliable approach to estimate soil hydraulic properties, and thus permits the simulation of soil water dynamics in both cropped and fallow soils in the field accurately. (C) 2009 Elsevier B.V. All rights reserved.

Nature of vegetation and building morphology characteristics across a city: influence on shadow patterns and mean radiant temperatures in London

Vegetation and building morphology characteristics are investigated at 19 sites on a north-south LiDAR transect across the megacity of London. Local maxima of mean building height and building plan area density at the city centre are evident. Surprisingly, the mean vegetation height (zv3) is also found to be highest in the city centre. From the LiDAR data various morphological parameters are derived as well as shadow patterns. Continuous images of the effects of buildings and of buildings plus vegetationon sky view factor (Ψ) are derived. A general reduction of Ψ is found, indicating the importance of including vegetation when deriving Ψ in urban areas. The contribution of vegetation to the shadowing at ground level is higher during summer than in autumn. Using these 3D data the influence on urban climate and mean radiant temperature (T mrt ) is calculated with SOLWEIG. The results from these simulations highlight that vegetation can be most effective at reducing heat stress within dense urban environments in summer. The daytime average T mrt is found to be lowest in the densest urban environments due to shadowing; foremost from buildings but also from trees. It is clearly shown that this method could be used to quantify the influence of vegetation on T mrt within the urban environment. The results presented in this paper highlight a number of possible climate sensitive planning practices for urban areas at the local scale (i.e. 102- 5 × 103 m).

An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.

Nocturnal cooling below a forest canopy: Model and evaluation

Nocturnal cooling of air within a forest canopy and the resulting temperature profile may drive local thermally driven motions, such as drainage flows, which are believed to impact measurements of ecosystem–atmosphere exchange. To model such flows, it is necessary to accurately predict the rate of cooling. Cooling occurs primarily due to radiative heat loss. However, much of the radiative loss occurs at the surface of canopy elements (leaves, branches, and boles of trees), while radiative divergence in the canopy air space is small due to high transmissivity of air. Furthermore, sensible heat exchange between the canopy elements and the air space is slow relative to radiative fluxes. Therefore, canopy elements initially cool much more quickly than the canopy air space after the switch from radiative gain during the day to radiative loss during the night. Thus in modeling air cooling within a canopy, it is not appropriate to neglect the storage change of heat in the canopy elements or even to assume equal rates of cooling of the canopy air and canopy elements. Here a simple parameterization of radiatively driven cooling of air within the canopy is presented, which accounts implicitly for radiative cooling of the canopy volume, heat storage in the canopy elements, and heat transfer between the canopy elements and the air. Simulations using this parameterization are compared to temperature data from the Morgan–Monroe State Forest (IN, USA) FLUXNET site. While the model does not perfectly reproduce the measured rates of cooling, particularly near the top of the canopy, the simulated cooling rates are of the correct order of magnitude.

Privatisation, empowerment and accountability: What are the policy implications for establishing effective farmer organisations?

Weak institutional development and information flows have constrained the extent to which the small-holder farming sector in developing countries can significantly drive growth and poverty reduction. Thisis despite widely implemented economic liberalisation policies focussing on market efficiency. Farmerorganisations are viewed as a potential means of addressing public and private institutional failure but thishas frequently been limited by inequalities in access to power and information. This article investigatestwo issues that have received little research attention to date: what role downward accountability plays inenabling farmer organisations to improve services and markets, and what influences the extent to whichdownward accountability is achieved. Kenya Tea Development Agency (KTDA), one of the largest farmerorganisations in the world (>400,000 farmers) is examined alongside wider literature. Mixed methodswere used including key informant interviews, and eight months of participant observation followedby a questionnaire survey. The article concludes that without effective downward accountability farmerorganisations can become characterised by institutions and mechanisms that favour elites, restrictedweak coordination and regulation, and manipulated information flows. This in turn reduces individuals’incentives to invest. If farmer organisations are to realise their potential as a means of enabling the small-holder sector to significantly contribute to economic growth and poverty reduction, policy and researchneeds to address key factors which influence accountability including: how to ensure initial processes information of farmer organisations establish appropriate structures and rules; strong state regulation toenhance corporate accountability; transparent information provision regarding actions of farmer organi-sation leaders; and the role independent non-government organisations can play. Consequently attentionneeds to focus on developing means of legitimising rights, building poor people’s capacity to challengeexclusion, and moving from rights to obligations regarding information provision.

Effective interventions: 25 years outside the mainstream

25 years ago when the Durham conferences were in full swing, I presented results of investigations on language and behaviour in autism. I tentatively proposed that early language in autism might tell us about the cognitive skills of people with ASD and the behaviour might lead to greater understanding of which brain systems might be affected. In this presentation, I will update these topics and present a summary of other work I have been involved with in attempting to improve the lives of people with autism and their families. Data on three people with autism at the early stages of speech development showed an unusual pattern of learning colour and number names early. One possibility was that this skill represented a sign of weak central coherence – they only attended to one dimension. Colleagues of mine were equally puzzled so we tried to find out if my results could be replicated – they were not (see Schafer, Williams & Smith, 2014). Instead we found this pattern was also seen in Down Syndrome, but that early vocabulary in autism was associated with low Colorado Meaningfulness at least in comprehension. The Colorado Meaningfulness of a word is a measure of how many words can be associated with it and often involve extensive use of context. Our data suggest that the number of contexts in which a particular word can appear has a role in determining vocabulary in ASD which is consistent with the weak central coherence theory of autism. In the course of this work I also came across a group of young people with autism who appeared to have a written vocabulary but not a spoken one. It seems possible that print might be a medium of communication when speech is not. Repetitive behaviour in autism remains a mystery. We can use functional analysis to determine why the behaviour occurs, but a worryingly large percentage of behaviours are described as being internally driven or sensory reinforced. What does that mean in terms of brain activity – could it be system analogous to epilepsy, where brain activity becomes inappropriately synchronised? At the moment I cannot claim to have solved this problem, but if sensation is a driver then sensory interventions should make a difference. Data from a recent study will be presented to suggest that for some individuals this is the case. Social behaviour remains the key however, and it remains to be seen whether it is possible for social behaviour to be aided. One route that has potential is direct teaching of skills through drama and working with others who do not have social difficulties of the same type. The picture is complicated by changes in social skills with age and experience, but the failure of people with ASD to interact when in settings of social contact is little researched.

Effective separation of Am(III) and Eu(III) from HNO3 Q1 Q2 solutions using CyMe4-BTPhen-functionalized silica-coated magnetic nanoparticles

It has been shown that CyMe4-BTPhen-functionalized silica-coated maghemite (c-Fe2O3) magnetic nanoparticles (MNPs) are capable of quantitative separation of Am(III) from Eu(III) from HNO3 solutions. These MNPs also show a small but significant selectivity for Am(III) over Cm(III) with a separation factor of around 2 in 4 M HNO3. The water molecule in the cavity of the BTPhen may also play an important part in the selectivity.

Bis amide-aromatic-ureasd - highly effective hydro- and organogelator systems

A series of hydro- and organo-supergelators have been synthesised via coupling of simple bis aromaticureas via alkyl amide linkages. These bis amide-aromatic-ureas exhibited reduced critical gelator concentrations, improved gelator stability, mechanical and dye removal properties for potential use in water purification, in comparison to related bis aromatic-ureas. Systematic structure studies via variation of the bis amide-aromatic-urea linker length as well as functionalization of the terminal aromatic moieties have enabled control over the gel properties.

A review of urban roughness sublayer turbulence

It is becoming increasingly important that we can understand and model flow processes in urban areas. Applications such as weather forecasting, air quality and sustainable urban development rely on accurate modelling of the interface between an urban surface and the atmosphere above. This review gives an overview of current understanding of turbulence generated by an urban surface up to a few building heights, the layer called the roughness sublayer (RSL). High quality datasets are also identified which can be used in the development of suitable parameterisations of the urban RSL. Datasets derived from physical and numerical modelling, and full-scale observations in urban areas now exist across a range of urban-type morphologies (e.g. street canyons, cubes, idealised and realistic building layouts). Results show that the urban RSL depth falls within 2 – 5 times mean building height and is not easily related to morphology. Systematic perturbations away from uniform layouts (e.g. varying building heights) have a significant impact on RSL structure and depth. Considerable fetch is required to develop an overlying inertial sublayer, where turbulence is more homogeneous, and some authors have suggested that the “patchiness” of urban areas may prevent inertial sublayers from developing at all. Turbulence statistics suggest similarities between vegetation and urban canopies but key differences are emerging. There is no consensus as to suitable scaling variables, e.g. friction velocity above canopy vs. square root of maximum Reynolds stress, mean vs. maximum building height. The review includes a summary of existing modelling practices and highlights research priorities.

Aerosol microphysics simulations of the Mt.~Pinatubo eruption with the UM-UKCA composition-climate model

We use a stratosphere–troposphere composition–climate model with interactive sulfur chemistry and aerosol microphysics, to investigate the effect of the 1991 Mount Pinatubo eruption on stratospheric aerosol properties. Satellite measurements indicate that shortly after the eruption, between 14 and 23 Tg of SO2 (7 to 11.5 Tg of sulfur) was present in the tropical stratosphere. Best estimates of the peak global stratospheric aerosol burden are in the range 19 to 26 Tg, or 3.7 to 6.7 Tg of sulfur assuming a composition of between 59 and 77 % H2SO4. In light of this large uncertainty range, we performed two main simulations with 10 and 20 Tg of SO2 injected into the tropical lower stratosphere. Simulated stratospheric aerosol properties through the 1991 to 1995 period are compared against a range of available satellite and in situ measurements. Stratospheric aerosol optical depth (sAOD) and effective radius from both simulations show good qualitative agreement with the observations, with the timing of peak sAOD and decay timescale matching well with the observations in the tropics and mid-latitudes. However, injecting 20 Tg gives a factor of 2 too high stratospheric aerosol mass burden compared to the satellite data, with consequent strong high biases in simulated sAOD and surface area density, with the 10 Tg injection in much better agreement. Our model cannot explain the large fraction of the injected sulfur that the satellite-derived SO2 and aerosol burdens indicate was removed within the first few months after the eruption. We suggest that either there is an additional alternative loss pathway for the SO2 not included in our model (e.g. via accommodation into ash or ice in the volcanic cloud) or that a larger proportion of the injected sulfur was removed via cross-tropopause transport than in our simulations. We also critically evaluate the simulated evolution of the particle size distribution, comparing in detail to balloon-borne optical particle counter (OPC) measurements from Laramie, Wyoming, USA (41° N). Overall, the model captures remarkably well the complex variations in particle concentration profiles across the different OPC size channels. However, for the 19 to 27 km injection height-range used here, both runs have a modest high bias in the lowermost stratosphere for the finest particles (radii less than 250 nm), and the decay timescale is longer in the model for these particles, with a much later return to background conditions. Also, whereas the 10 Tg run compared best to the satellite measurements, a significant low bias is apparent in the coarser size channels in the volcanically perturbed lower stratosphere. Overall, our results suggest that, with appropriate calibration, aerosol microphysics models are capable of capturing the observed variation in particle size distribution in the stratosphere across both volcanically perturbed and quiescent conditions. Furthermore, additional sensitivity simulations suggest that predictions with the models are robust to uncertainties in sub-grid particle formation and nucleation rates in the stratosphere.

Sunlit fractions on urban facets: impact of spatial resolution and approach

The extent of the surface area sunlit is critical for radiative energy exchanges and therefore for a wide range of applications that require urban land surface models (ULSM), ranging from human comfort to weather forecasting. Here a computational demanding shadow casting algorithm is used to assess the capability of a simple single-layer urban canopy model, which assumes an infinitely long rotating canyon (ILC), to reproduce sunlit areas on roof and roads over central London. Results indicate that the sunlit roads areas are well-represented but somewhat smaller using an ILC, while sunlit roofs areas are consistently larger, especially for dense urban areas. The largest deviations from real world sunlit areas are found for roofs during mornings and evenings. Indications that sunlit fractions on walls are overestimated using an ILC during mornings and evenings are found. The implications of these errors are dependent on the application targeted. For example, (independent of albedo) ULSMs used in numerical weather prediction applying ILC representation of the urban form will overestimate outgoing shortwave radiation from roofs due to the overestimation of sunlit fraction of the roofs. Complications of deriving height to width ratios from real world data are also discussed.

Simplified estimation of ray-path mirroring height for HF radiowaves reflected from the ionospheric F-region

From Milsom's equations, which describe the geometry of ray-path hops reflected from the ionospheric F-layer, algorithms for the simplified estimation of mirror-reflection height are developed. These allow for hop length and the effects of variations in underlying ionisation (via the ratio of the F2- and E-layer critical frequencies) and F2-layer peak height (via the M(3000)F2-factor). Separate algorithms are presented which are applicable to a range of signal frequencies about the FOT and to propagation at the MUF. The accuracies and complexities of the algorithms are compared with those inherent in the use of a procedure based on an equation developed by Shimazaki.