Correlations of control variables in variational data assimilation

Variational data assimilation systems for numerical weather prediction rely on a transformation of model variables to a set of control variables that are assumed to be uncorrelated. Most implementations of this transformation are based on the assumption that the balanced part of the flow can be represented by the vorticity. However, this assumption is likely to break down in dynamical regimes characterized by low Burger number. It has recently been proposed that a variable transformation based on potential vorticity should lead to control variables that are uncorrelated over a wider range of regimes. In this paper we test the assumption that a transform based on vorticity and one based on potential vorticity produce an uncorrelated set of control variables. Using a shallow-water model we calculate the correlations between the transformed variables in the different methods. We show that the control variables resulting from a vorticity-based transformation may retain large correlations in some dynamical regimes, whereas a potential vorticity based transformation successfully produces a set of uncorrelated control variables. Calculations of spatial correlations show that the benefit of the potential vorticity transformation is linked to its ability to capture more accurately the balanced component of the flow.

A policy critique of Stansted Airport's expansion to 25 million passengers per annum (MPPA)

In this paper we undertake a preliminary assessment of the regional planning and development implications of BAA Stansted Airport’s planning permission to grow to 25 million passengers per annum (mppa) by 2010. Our concern is not simply to consider the overall growth of the airport on the airport site itself but the nature and type of growth both on- and off-site. In this document we focus on the submitted planning permission documents and test them. The methodology we employed was to draw on published and unpublished numerical estimates of the airport’s growth – particularly including estimates produced by the airport owner, BAA, and their economic and planning consultants DTZ Pieda - and critically, and systematically analyse their figures. We adopted this approach because unless the figures which were employed in the initial calculations were correct then all of the subsequent projections which flow from them - and the polices which could then be based on them – could be flawed. The analysis is divided into two parts – firstly, are the growth forecasts correct?; and secondly, what do these forecasts actually mean in developmental terms? In effect, what we have done is to produce a critique of the existing body of evidence by questioning underpinning assumptions and then draw some preliminary conclusions for the region based on this analysis. A major focus of this report has been analyse the figures involved in the planning application to expand Stansted to 25mppa. Ironically, one of our key findings, that the local impact of Stansted’s proposed expansion in employment terms might well be less than was originally thought, might make it easier to gain the acceptance of the relevant local authorities involved to allow the development to take place. Our main overall findings are that the BAA projections over-estimate the local employment impact of the airport’s proposed growth and under-estimate its potential regional ‘transportation’ employment effect. These two findings are, of course, related to each other in important ways, and we also feel that they have potentially significant medium and long-term economic, competitiveness and planning policy implications for the East of England region

Belowground biomass functions and expansion factors in high elevation Norway spruce

Biomass allocation to above- and belowground compartments in trees is thought to be affected by growth conditions. To assess the strength of such influences, we sampled six Norway spruce forest stands growing at higher altitudes. Within these stands, we randomly selected a total of 77 Norway spruce trees and measured volume and biomass of stem, above- and belowground stump and all roots over 0.5 cm diameter. A comparison of our observations with models parameterised for lower altitudes shows that models developed for specific conditions may be applicable to other locations. Using our observations, we developed biomass functions (BF) and biomass conversion and expansion factors (BCEF) linking belowground biomass to stem parameters. While both BF and BCEF are accurate in belowground biomass predictions, using BCEF appears more promising as such factors can be readily used with existing forest inventory data to obtain estimates of belowground biomass stock. As an example, we show how BF and BCEF developed for individual trees can be used to estimate belowground biomass at the stand level. In combination with existing aboveground models, our observations can be used to quantify total standing biomass of high altitude Norway spruce stands.

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].

Metal cluster expansion by addition of low valent group 14 reagents: III. Reaction of bis[bis(trimethylsilyl)methyl]tin with M3(CO)12 or its derivatives (M = Fe, Ru, or Os): the crystal and molecular structures of M3(μ-SnR2)2(CO)10(R = CH(SiMe3)2; M = Ru or Os)

The stannylene [SnR2] (R = CH(SiMe3)2) reacts in different ways with the three dodecacarbonyls of the iron triad: [Fe3(CO)12] gives [Fe2(CO)8(μ-SnR2)], [Ru3(CO)12] gives the planar pentametallic cluster [Ru3(CO)10(μ-SnR2)2], for which a full structural analysis is reported, while [Os3(CO)12] fails to react. Different products are also obtained from three nitrile derivatives: [Fe3-(CO)11(MeCN)] gives [Fe2(CO)6(μ-SnR2)2], which has a structure significantly different from that of known Fe2Sn2 clusters, [Ru3(CO)10(MeCN)2] gives the pentametallic cluster described above, while [Os3(CO)10(MeCN)2] gives the isostructural osmium analogue, which shows the unusual feature of a CO group bridging two osmium atoms.

Metal cluster expansion by addition of low-valent group 4B reagents: crystal and molecular structure of [Os3SnH2(CO)10{CH(SiMe3)2}2]; a compound with hydrogen bridging osmium and tin

Cluster expansion of [Os3H2(CO)10] with [SnR2][R = CH(SiMe3)2] take place in high yield to give [Os3SnH2(CO)10R2], the first closed triosmium–main-group metal cluster to be structurally characterized; a novel feature is the presence of a hydrogen atom bridging the tin atom and one of the osmium atoms.

A simple column model to explore anticipated problems in variational assimilation of satellite observations

We investigate a simplified form of variational data assimilation in a fully nonlinear framework with the aim of extracting dynamical development information from a sequence of observations over time. Information on the vertical wind profile, w(z ), and profiles of temperature, T (z , t), and total water content, qt (z , t), as functions of height, z , and time, t, are converted to brightness temperatures at a single horizontal location by defining a two-dimensional (vertical and time) variational assimilation testbed. The profiles of T and qt are updated using a vertical advection scheme. A basic cloud scheme is used to obtain the fractional cloud amount and, when combined with the temperature field, this information is converted into a brightness temperature, using a simple radiative transfer scheme. It is shown that our model exhibits realistic behaviour with regard to the prediction of cloud, but the effects of nonlinearity become non-negligible in the variational data assimilation algorithm. A careful analysis of the application of the data assimilation scheme to this nonlinear problem is presented, the salient difficulties are highlighted, and suggestions for further developments are discussed.

Conditioning of incremental variational data assimilation, with application to the Met Office system

Implementations of incremental variational data assimilation require the iterative minimization of a series of linear least-squares cost functions. The accuracy and speed with which these linear minimization problems can be solved is determined by the condition number of the Hessian of the problem. In this study, we examine how different components of the assimilation system influence this condition number. Theoretical bounds on the condition number for a single parameter system are presented and used to predict how the condition number is affected by the observation distribution and accuracy and by the specified lengthscales in the background error covariance matrix. The theoretical results are verified in the Met Office variational data assimilation system, using both pseudo-observations and real data.

Conditioning and preconditioning of the variational data assimilation problem

Numerical weather prediction (NWP) centres use numerical models of the atmospheric flow to forecast future weather states from an estimate of the current state. Variational data assimilation (VAR) is used commonly to determine an optimal state estimate that miminizes the errors between observations of the dynamical system and model predictions of the flow. The rate of convergence of the VAR scheme and the sensitivity of the solution to errors in the data are dependent on the condition number of the Hessian of the variational least-squares objective function. The traditional formulation of VAR is ill-conditioned and hence leads to slow convergence and an inaccurate solution. In practice, operational NWP centres precondition the system via a control variable transform to reduce the condition number of the Hessian. In this paper we investigate the conditioning of VAR for a single, periodic, spatially-distributed state variable. We present theoretical bounds on the condition number of the original and preconditioned Hessians and hence demonstrate the improvement produced by the preconditioning. We also investigate theoretically the effect of observation position and error variance on the preconditioned system and show that the problem becomes more ill-conditioned with increasingly dense and accurate observations. Finally, we confirm the theoretical results in an operational setting by giving experimental results from the Met Office variational system.

Biomass functions and expansion factors in young Norway spruce (Picea abies [L.] Karst) trees

Roots, stems, branches and needles of 160 Norway spruce trees younger than 10 years were sampled in seven forest stands in central Slovakia in order to establish their biomassfunctions (BFs) and biomassexpansionfactors (BEFs). We tested three models for each biomass pool based on the stem base diameter, tree height and the two parameters combined. BEF values decreased for all spruce components with increasing height and diameter, which was most evident in very young trees under 1 m in height. In older trees, the values of BEFs did tend to stabilise at the height of 3–4 m. We subsequently used the BEFs to calculate dry biomass of the stands based on average stem base diameter and tree height. Total stand biomass grew with increasing age of the stands from about 1.0 Mg ha−1 at 1.5 years to 44.3 Mg ha−1 at 9.5 years. The proportion of stem and branch biomass was found to increase with age, while that of needles was fairly constant and the proportion of root biomass did decrease as the stands grew older.

'Agricultural poverty' and the expansion of artisanal mining in sub-Saharan Africa: Experiences from southwest Mali and southeast Ghana

Why do people engage in artisanal and small-scale mining (ASM) – labour-intensive mineral extraction and processing activity – across sub-Saharan Africa? This paper argues that ‘agricultural poverty’, or hardship induced by an over-dependency on farming for survival, has fuelled the recent rapid expansion of ASM operations throughout the region. The diminished viability of smallholder farming in an era of globalization and overreliance on rain-fed crop production restricted by seasonality has led hundreds of thousands of rural African families to ‘branch out’ into ASM, a move made to secure supplementary incomes. Experiences from Komana West in Southwest Mali and East Akim District in Southeast Ghana are drawn upon to illustrate how a movement into the ASM economy has impacted farm families, economically, in many rural stretches of sub-Saharan Africa.

Ex vivo expansion of limbal stem cells is affected by substrate properties

Limbal epithelial stem cells play a key role in the maintenance and regulation of the corneal surface. Damage or destruction of these cells results in vascularisation and corneal opacity. Subsequent limbal stem cell transplantation requires an ex vivo expansion step and preserving cells in an undifferentiated state remains vital. In this report we seek to control the phenotype of limbal epithelial stem cells by the novel application of compressed collagen substrates. We have characterised the mechanical and surface properties of conventional collagen gels using shear rheology and scanning electron microscopy. In doing so, we provide evidence to show that compressive load can improve the stiffness of collagen substrates. In addition Western blotting and immunohistochemistry display increased cytokeratin 3 (CK3) protein expression relating to limbal epithelial cell differentiation on stiff collagen substrates. Such gels with an elastic modulus of 2900 Pa supported a significantly higher number of cells than less stiff collagen gels (3 Pa). These findings have substantial influence in the development of ocular surface constructs or experimental models particularly in the fields of stem cell research, tissue engineering and regenerative medicine.

An item response analysis of the motor and behavioral subscales of the unified Huntington's disease rating scale in huntington disease gene expansion carriers

Although the Unified Huntington's Disease Rating Scale (UHDRS) is widely used in the assessment of Huntington disease (HD), the ability of individual items to discriminate individual differences in motor or behavioral manifestations has not been extensively studied in HD gene expansion carriers without a motor-defined clinical diagnosis (ie, prodromal-HD or prHD). To elucidate the relationship between scores on individual motor and behavioral UHDRS items and total score for each subscale, a nonparametric item response analysis was performed on retrospective data from 2 multicenter longitudinal studies. Motor and behavioral assessments were supplied for 737 prHD individuals with data from 2114 visits (PREDICT-HD) and 686 HD individuals with data from 1482 visits (REGISTRY). Option characteristic curves were generated for UHDRS subscale items in relation to their subscale score. In prHD, overall severity of motor signs was low, and participants had scores of 2 or above on very few items. In HD, motor items that assessed ocular pursuit, saccade initiation, finger tapping, tandem walking, and to a lesser extent, saccade velocity, dysarthria, tongue protrusion, pronation/supination, Luria, bradykinesia, choreas, gait, and balance on the retropulsion test were found to discriminate individual differences across a broad range of motor severity. In prHD, depressed mood, anxiety, and irritable behavior demonstrated good discriminative properties. In HD, depressed mood demonstrated a good relationship with the overall behavioral score. These data suggest that at least some UHDRS items appear to have utility across a broad range of severity, although many items demonstrate problematic features.

Wave scattering by an axisymmetric ice floe of varying thickness

The problem of water wave scattering by a circular ice floe, floating in fluid of finite depth, is formulated and solved numerically. Unlike previous investigations of such situations, here we allow the thickness of the floe (and the fluid depth) to vary axisymmetrically and also incorporate a realistic non-zero draught. A numerical approximation to the solution of this problem is obtained to an arbitrary degree of accuracy by combining a Rayleigh–Ritz approximation of the vertical motion with an appropriate variational principle. This numerical solution procedure builds upon the work of Bennets et al. (2007, J. Fluid Mech., 579, 413–443). As part of the numerical formulation, we utilize a Fourier cosine expansion of the azimuthal motion, resulting in a system of ordinary differential equations to solve in the radial coordinate for each azimuthal mode. The displayed results concentrate on the response of the floe rather than the scattered wave field and show that the effects of introducing the new features of varying floe thickness and a realistic draught are significant.

L1-regularisation for ill-posed problems in variational data assimilation

We consider four-dimensional variational data assimilation (4DVar) and show that it can be interpreted as Tikhonov or L2-regularisation, a widely used method for solving ill-posed inverse problems. It is known from image restoration and geophysical problems that an alternative regularisation, namely L1-norm regularisation, recovers sharp edges better than L2-norm regularisation. We apply this idea to 4DVar for problems where shocks and model error are present and give two examples which show that L1-norm regularisation performs much better than the standard L2-norm regularisation in 4DVar.