Sparse polynomial approximation in positive order Sobolev spaces with bounded mixed derivatives and applications to elliptic problems with random loading

In the present paper we study the approximation of functions with bounded mixed derivatives by sparse tensor product polynomials in positive order tensor product Sobolev spaces. We introduce a new sparse polynomial approximation operator which exhibits optimal convergence properties in L2 and tensorized View the MathML source simultaneously on a standard k-dimensional cube. In the special case k=2 the suggested approximation operator is also optimal in L2 and tensorized H1 (without essential boundary conditions). This allows to construct an optimal sparse p-version FEM with sparse piecewise continuous polynomial splines, reducing the number of unknowns from O(p2), needed for the full tensor product computation, to View the MathML source, required for the suggested sparse technique, preserving the same optimal convergence rate in terms of p. We apply this result to an elliptic differential equation and an elliptic integral equation with random loading and compute the covariances of the solutions with View the MathML source unknowns. Several numerical examples support the theoretical estimates.

Carbon stable isotope analysis of cereal remains as a way to reconstruct water availability: preliminary results

Reconstructing past water availability, both as rainfall and irrigation, is important to answer questions about the way society reacts to climate and its changes and the role of irrigation in the development of social complexity. Carbon stable isotope analysis of archaeobotanical remains is a potentially valuable method for reconstructing water availability. To further define the relationship between water availability and plant carbon isotope composition and to set up baseline values for the Southern Levant, grains of experimentally grown barley and sorghum were studied. The cereal crops were grown at three stations under five different irrigation regimes in Jordan. Results indicate that a positive but weak relationship exists between irrigation regime and total water input of barley grains, but no relationship was found for sorghum. The relationship for barley is site-specific and inter-annual variation was present at Deir ‘Alla, but not at Ramtha and Khirbet as-Samra.

An exploration of the equivalent weights particle filter

Particle filters are fully non-linear data assimilation techniques that aim to represent the probability distribution of the model state given the observations (the posterior) by a number of particles. In high-dimensional geophysical applications the number of particles required by the sequential importance resampling (SIR) particle filter in order to capture the high probability region of the posterior, is too large to make them usable. However particle filters can be formulated using proposal densities, which gives greater freedom in how particles are sampled and allows for a much smaller number of particles. Here a particle filter is presented which uses the proposal density to ensure that all particles end up in the high probability region of the posterior probability density function. This gives rise to the possibility of non-linear data assimilation in large dimensional systems. The particle filter formulation is compared to the optimal proposal density particle filter and the implicit particle filter, both of which also utilise a proposal density. We show that when observations are available every time step, both schemes will be degenerate when the number of independent observations is large, unlike the new scheme. The sensitivity of the new scheme to its parameter values is explored theoretically and demonstrated using the Lorenz (1963) model.

Quantifying errors due to frequency changes and target location uncertainty for radar refractivity retrievals

Radar refractivity retrievals can capture near-surface humidity changes, but noisy phase changes of the ground clutter returns limit the accuracy for both klystron- and magnetron-based systems. Observations with a C-band (5.6 cm) magnetron weather radar indicate that the correction for phase changes introduced by local oscillator frequency changes leads to refractivity errors no larger than 0.25 N units: equivalent to a relative humidity change of only 0.25% at 20°C. Requested stable local oscillator (STALO) frequency changes were accurate to 0.002 ppm based on laboratory measurements. More serious are the random phase change errors introduced when targets are not at the range-gate center and there are changes in the transmitter frequency (ΔfTx) or the refractivity (ΔN). Observations at C band with a 2-μs pulse show an additional 66° of phase change noise for a ΔfTx of 190 kHz (34 ppm); this allows the effect due to ΔN to be predicted. Even at S band with klystron transmitters, significant phase change noise should occur when a large ΔN develops relative to the reference period [e.g., ~55° when ΔN = 60 for the Next Generation Weather Radar (NEXRAD) radars]. At shorter wavelengths (e.g., C and X band) and with magnetron transmitters in particular, refractivity retrievals relative to an earlier reference period are even more difficult, and operational retrievals may be restricted to changes over shorter (e.g., hourly) periods of time. Target location errors can be reduced by using a shorter pulse or identified by a new technique making alternate measurements at two closely spaced frequencies, which could even be achieved with a dual–pulse repetition frequency (PRF) operation of a magnetron transmitter.

On the spectra and pseudospectra of a class of non-self-adjoint random matrices and operators

In this paper we develop and apply methods for the spectral analysis of non-selfadjoint tridiagonal infinite and finite random matrices, and for the spectral analysis of analogous deterministic matrices which are pseudo-ergodic in the sense of E. B. Davies (Commun. Math. Phys. 216 (2001), 687–704). As a major application to illustrate our methods we focus on the “hopping sign model” introduced by J. Feinberg and A. Zee (Phys. Rev. E 59 (1999), 6433–6443), in which the main objects of study are random tridiagonal matrices which have zeros on the main diagonal and random ±1’s as the other entries. We explore the relationship between spectral sets in the finite and infinite matrix cases, and between the semi-infinite and bi-infinite matrix cases, for example showing that the numerical range and p-norm ε - pseudospectra (ε > 0, p ∈ [1,∞] ) of the random finite matrices converge almost surely to their infinite matrix counterparts, and that the finite matrix spectra are contained in the infinite matrix spectrum Σ. We also propose a sequence of inclusion sets for Σ which we show is convergent to Σ, with the nth element of the sequence computable by calculating smallest singular values of (large numbers of) n×n matrices. We propose similar convergent approximations for the 2-norm ε -pseudospectra of the infinite random matrices, these approximations sandwiching the infinite matrix pseudospectra from above and below.

Network optimization for enhanced resilience of urban heat island measurements

The urban heat island is a well-known phenomenon that impacts a wide variety of city operations. With greater availability of cheap meteorological sensors, it is possible to measure the spatial patterns of urban atmospheric characteristics with greater resolution. To develop robust and resilient networks, recognizing sensors may malfunction, it is important to know when measurement points are providing additional information and also the minimum number of sensors needed to provide spatial information for particular applications. Here we consider the example of temperature data, and the urban heat island, through analysis of a network of sensors in the Tokyo metropolitan area (Extended METROS). The effect of reducing observation points from an existing meteorological measurement network is considered, using random sampling and sampling with clustering. The results indicated the sampling with hierarchical clustering can yield similar temperature patterns with up to a 30% reduction in measurement sites in Tokyo. The methods presented have broader utility in evaluating the robustness and resilience of existing urban temperature networks and in how networks can be enhanced by new mobile and open data sources.

A scalable expressive ensemble learning using Random Prism: a MapReduce approach

The induction of classification rules from previously unseen examples is one of the most important data mining tasks in science as well as commercial applications. In order to reduce the influence of noise in the data, ensemble learners are often applied. However, most ensemble learners are based on decision tree classifiers which are affected by noise. The Random Prism classifier has recently been proposed as an alternative to the popular Random Forests classifier, which is based on decision trees. Random Prism is based on the Prism family of algorithms, which is more robust to noise. However, like most ensemble classification approaches, Random Prism also does not scale well on large training data. This paper presents a thorough discussion of Random Prism and a recently proposed parallel version of it called Parallel Random Prism. Parallel Random Prism is based on the MapReduce programming paradigm. The paper provides, for the first time, novel theoretical analysis of the proposed technique and in-depth experimental study that show that Parallel Random Prism scales well on a large number of training examples, a large number of data features and a large number of processors. Expressiveness of decision rules that our technique produces makes it a natural choice for Big Data applications where informed decision making increases the user’s trust in the system.

Effects of water availability on free amino acids, sugars, and acrylamide-forming potential in potato

Irrigation is used frequently in potato cultivation to maximize yield, but water availability may also affect the composition of the crop, with implications for processing properties and food safety. Five varieties of potatoes, including drought-tolerant and -sensitive types, which had been grown with and without irrigation, were analyzed to show the effect of water supply on concentrations of free asparagine, other free amino acids, and sugars and on the acrylamide-forming potential of the tubers. Two varieties were also analyzed under more severe drought stress in a glasshouse. Water availability had profound effects on tuber free amino acid and sugar concentrations, and it was concluded that potato farmers should irrigate only if necessary to maintain the health and yield of the crop, because irrigation may increase the acrylamide-forming potential of potatoes. Even mild drought stress caused significant changes in composition, but these differed from those caused by more extreme drought stress. Free proline concentration, for example, increased in the field-grown potatoes of one variety from 7.02 mmol/kg with irrigation to 104.58 mmol/kg without irrigation, whereas free asparagine concentration was not affected significantly in the field but almost doubled from 132.03 to 242.26 mmol/kg in response to more severe drought stress in the glasshouse. Furthermore, the different genotypes were affected in dissimilar fashion by the same treatment, indicating that there is no single, unifying potato tuber drought stress response.

Endophytic bacterial community composition in wheat (Triticum aestivum) is determined by plant tissue type, developmental stage and soil nutrient availability

Aims: To understand effects of tissue type, growth stage and soil fertilisers on bacterial endophyte communities of winter wheat (Triticum aestivum cv. Hereward). Methods: Endophytes were isolated from wheat grown under six fertiliser conditions in the long term Broadbalk Experiment at Rothamsted Research, UK. Samples were taken in May and July from root and leaf tissues. Results: Root and leaf communities differed in abundance and composition of endophytes. Endophytes were most abundant in roots and the Proteobacteria were most prevalent. In contrast, Firmicutes and Actinobacteria, the Gram positive phyla, were most prevalent in the leaves. Both fertiliser treatment and sample time influenced abundance and relative proportions of each phylum and genus in the endosphere. A higher density of endophytes was found in the Nil input treatment plants. Conclusions: Robust isolation techniques and stringent controls are critical for accurate recovery of endophytes. The plant tissue type, plant growth stage, and soil fertiliser treatment all contribute to the composition of the endophytic bacterial community in wheat. These results should help facilitate targeted development of endophytes for beneficial applications in agriculture.

The equivalent-weights particle filter in a high-dimensional system

In general, particle filters need large numbers of model runs in order to avoid filter degeneracy in high-dimensional systems. The recently proposed, fully nonlinear equivalent-weights particle filter overcomes this requirement by replacing the standard model transition density with two different proposal transition densities. The first proposal density is used to relax all particles towards the high-probability regions of state space as defined by the observations. The crucial second proposal density is then used to ensure that the majority of particles have equivalent weights at observation time. Here, the performance of the scheme in a high, 65 500 dimensional, simplified ocean model is explored. The success of the equivalent-weights particle filter in matching the true model state is shown using the mean of just 32 particles in twin experiments. It is of particular significance that this remains true even as the number and spatial variability of the observations are changed. The results from rank histograms are less easy to interpret and can be influenced considerably by the parameter values used. This article also explores the sensitivity of the performance of the scheme to the chosen parameter values and the effect of using different model error parameters in the truth compared with the ensemble model runs.

The effect of the equivalent-weights particle filter on dynamical balance in a primitive equation model

The disadvantage of the majority of data assimilation schemes is the assumption that the conditional probability density function of the state of the system given the observations [posterior probability density function (PDF)] is distributed either locally or globally as a Gaussian. The advantage, however, is that through various different mechanisms they ensure initial conditions that are predominantly in linear balance and therefore spurious gravity wave generation is suppressed. The equivalent-weights particle filter is a data assimilation scheme that allows for a representation of a potentially multimodal posterior PDF. It does this via proposal densities that lead to extra terms being added to the model equations and means the advantage of the traditional data assimilation schemes, in generating predominantly balanced initial conditions, is no longer guaranteed. This paper looks in detail at the impact the equivalent-weights particle filter has on dynamical balance and gravity wave generation in a primitive equation model. The primary conclusions are that (i) provided the model error covariance matrix imposes geostrophic balance, then each additional term required by the equivalent-weights particle filter is also geostrophically balanced; (ii) the relaxation term required to ensure the particles are in the locality of the observations has little effect on gravity waves and actually induces a reduction in gravity wave energy if sufficiently large; and (iii) the equivalent-weights term, which leads to the particles having equivalent significance in the posterior PDF, produces a change in gravity wave energy comparable to the stochastic model error. Thus, the scheme does not produce significant spurious gravity wave energy and so has potential for application in real high-dimensional geophysical applications.

A moment problem for random discrete measures

Let X be a locally compact Polish space. A random measure on X is a probability measure on the space of all (nonnegative) Radon measures on X. Denote by K(X) the cone of all Radon measures η on X which are of the form η =

Twin experiments with the equivalent weights particle filter and HadCM3

This paper investigates the use of a particle filter for data assimilation with a full scale coupled ocean–atmosphere general circulation model. Synthetic twin experiments are performed to assess the performance of the equivalent weights filter in such a high-dimensional system. Artificial 2-dimensional sea surface temperature fields are used as observational data every day. Results are presented for different values of the free parameters in the method. Measures of the performance of the filter are root mean square errors, trajectories of individual variables in the model and rank histograms. Filter degeneracy is not observed and the performance of the filter is shown to depend on the ability to keep maximum spread in the ensemble.

Heathland restoration on former agricultural land: effects of artificial acidification on the availability and uptake of toxic metal cations

Lowland heath is an internationally important habitat type that has greatly declined in abundance throughout Western Europe. In recent years this has led to a growing interest in the restoration of heathland on agricultural land. This generally requires the use of chemical treatments to return soil chemical conditions to those appropriate for the support of heathland ecosystems. However, the potential for negative impacts on the environment due to the potential of these treatments to increase the availability of trace metals via raised soil acidity requires investigation. A large-scale field study investigated the effect of two chemical treatments used in heathland restoration, elemental sulphur and ferrous sulphate, on soil acidity and whether it is possible to predict the effect of the treatments on availability of two potentially toxic cations (Al and Cd) in the soil along with their subsequent accumulation in the shoots of the grass Agrostis capillaris. Results showed that both treatments decreased soil pH, but that only elemental sulphur produced a pH similar to heathland soil. The availability of Al, measured by extraction with 1 M ammonium nitrate, could not be predicted by soil pH, depth in the soil and total Al concentration in the soil. By contrast, availability of Cd could be predicted from these three variables. Concentrations of both Al and Cd in the shoots of A. capillaris showed no significant relationship with the extractable concentration in the soil. Results are discussed in light of the possible environmental impacts of the chemical restoration techniques.