Non-testability of equal weights spatial dependence

We show that any invariant test for spatial autocorrelation in a spatial error or spatial lag model with equal weights matrix has power equal to size. This result holds under the assumption of an elliptical distribution. Under Gaussianity, we also show that any test whose power is larger than its size for at least one point in the parameter space must be biased.

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.

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.

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.

Effects of forage source and extruded linseed supplementation on methane emissions from growing dairy cattle of differing body weights

Changes in diet carbohydrate amount and type (i.e., starch vs. fiber) and dietary oil supplements can affect ruminant methane emissions. Our objectives were to measure methane emissions, whole-tract digestibility, and energy and nitrogen utilization from growing dairy cattle at 2 body weight (BW) ranges, fed diets containing either high maize silage (MS) or high grass silage (GS), without or with supplemental oil from extruded linseed (ELS). Four Holstein-Friesian heifers aged 13 mo (BW range from start to finish of 382 to 526 kg) were used in experiment 1, whereas 4 lighter heifers aged 12 mo (BW range from start to finish of 292 to 419 kg) were used in experiment 2. Diets were fed as total mixed rations with forage dry matter (DM) containing high MS or high GS and concentrates in proportions (forage:concentrate, DM basis) of either 75:25 (experiment 1) or 60:40 (experiment 2), respectively. Diets were supplemented without or with ELS (Lintec[AU1: Add manufacturer name and location.]; 260 g of oil/ kg of DM) at 6% of ration DM. Each experiment was a 4 × 4 Latin square design with 33-d periods, with measurements during d 29 to 33 while animals were housed in respiration chambers. Heifers fed MS at a heavier BW (experiment 1) emitted 20% less methane per unit of DM intake (yield) compared with GS (21.4 vs. 26.6, respectively). However, when repeated with heifers of a lower BW (experiment 2), methane yield did not differ between the 2 diets (26.6 g/kg of DM intake). Differences in heifer BW had no overall effect on methane emissions, except when expressed as grams per kilogram of digestible organic matter (OMD) intake (32.4 vs. 36.6, heavy vs. light heifers). Heavier heifers fed MS in experiment 1 had a greater DM intake (9.4 kg/d) and lower OMD (755 g/kg), but no difference in N utilization (31% of N intake) compared with heifers fed GS (7.9 kg/d and 799 g/kg, respectively). Tissue energy retention was nearly double for heifers fed MS compared with GS in experiment 1 (15 vs. 8% of energy intake, respectively). Heifers fed MS in experiment 2 had similar DM intake (7.2 kg/d) and retention of energy (5% of intake energy) and N (28% of N intake), compared with GS-fed heifers, but OMD was lower (741 vs. 765 g/kg, respectively). No effect of ELS was noted on any of the variables measured, irrespective of animal BW, and this was likely due to the relatively low amount of supplemental oil provided. Differences in heifer BW did not markedly influence dietary effects on methane emissions. Differences in methane yield were attributable to differences in dietary starch and fiber composition associated with forage type and source.

Implicit equal-weights particle filter

Filter degeneracy is the main obstacle for the implementation of particle filter in non-linear high-dimensional models. A new scheme, the implicit equal-weights particle filter (IEWPF), is introduced. In this scheme samples are drawn implicitly from proposal densities with a different covariance for each particle, such that all particle weights are equal by construction. We test and explore the properties of the new scheme using a 1,000-dimensional simple linear model, and the 1,000-dimensional non-linear Lorenz96 model, and compare the performance of the scheme to a Local Ensemble Kalman Filter. The experiments show that the new scheme can easily be implemented in high-dimensional systems and is never degenerate, with good convergence properties in both systems.

Utilisation by upland rice of plant residue- and fertiliser-phosphorus in two tropical acid soils

A dual isotopic technique was used to assess the effects of soil type, and residues of Gliricidia sepium, without and with added fertiliser-P on the utilisation of P. Upland rice (Oryza sativa) was grown for 70 days in two tropical acid soils of different P sorbing capacity and P status. Uniformly P-32-labelled soils were treated with inorganic fertiliser-P tagged with P-33, Gliricidia sepium residue applied at planting and 3 weeks earlier, and in a combination of fertiliser-P and Gliricidia applied at and 3 weeks before planting. There were significant responses of shoot and root weights, and total P uptake to Gliricidia- and/or fertiliser-P addition in the Ultisol (low P status) but not the Oxisol (high P status), suggesting that P in the latter soil was not yield limiting, despite the high standard P requirement. Similarly, incorporation of Gliricidia three weeks before planting further increased shoot weight only in the Ultisol. There were generally higher proportions, quantities and percent utilisations of the Gliricidia- P and fertiliser-P in the Ultisol than in the Oxisol. Gliricidia significantly increased the utilisation of fertiliser-P only in the Ultisol. However, early application of Gliricidia increased Gliricidia- P but not fertiliser-P utilisation in the Ultisol. Added fertiliser-P did not influence Gliricidia- P utilisation.

Influence of sodium chloride on seed germination and seedling root growth of cotton (Gossypium hirsutum L.)

Response of cotton (Gossypium hirsutum L. cv. NIAB-78) to salinity, in terms of seed germination, seedling root growth and root Na+ and K+ content was determined in a laboratory experiment. Cotton seeds were exposed to increasing salinity levels using germination water with Sodium chloride concentrations of 0, 50, 100, 150 and 200 mM, to provide different degrees of salt stress. Germinated seeds were counted and roots were harvested at 24, 48, 72 and 96 h after the start of the experiment. It appeared that seed germination was only slightly affected by an increase in salinity (in most cases the differences between treatment were non-significant), whereas root length, root growth rate, root fresh and dry weights were severely affected, generally highly significant differences in these variables were found for comparisons involving most combinations of salinity levels, in particular with increased incubation period. K+ contents decreased with increasing salinity levels, although differences in K+ content were only significant when comparing the control and the 4 salinity levels. Na+ content of the roots increased with increasing levels of NaCl in the germination water, suggesting an exchange of K+ for Na+. The ratio K+/Na+ strongly decreased with rising levels of salinity from around 4.5 for the control to similar to 1 at 200 mM NaCl.

Temperature-responsive properties and drug solubilization capacity of amphiphilic copolymers based on N-vinylpyrrolidone and vinyl propyl ether

A series of amphiphilic copolymers were synthesized by free-radical copolymerization of N-vinylpyrrolidone (NVP) with vinyl propyl ether (VPE), and the structure of the copolymers was characterized by elemental analysis and gel permeation chromatography. The reactivity of VPE in copolymerization was found to be significantly lower than the reactivity of NVP, which resulted in a decrease of copolymers’ yields and molecular weights with higher content of VPE in the feed mixture. An investigation of the behavior of the copolymers in aqueous solutions at different temperatures by dynamic light scattering revealed the presence of lower critical solution temperature, which depending on the content of VPE ranged within 23−38 °C. Aqueous solutions of these copolymers were studied by fluorescent spectroscopy with pyrene as a polarity probe to reveal the formation of hydrophobic domains. The copolymers were found to be useful for enhancing the solubility of riboflavin in water.

A new linear initialization in SOM for biomolecular data

In the past decade, the amount of data in biological field has become larger and larger; Bio-techniques for analysis of biological data have been developed and new tools have been introduced. Several computational methods are based on unsupervised neural network algorithms that are widely used for multiple purposes including clustering and visualization, i.e. the Self Organizing Maps (SOM). Unfortunately, even though this method is unsupervised, the performances in terms of quality of result and learning speed are strongly dependent from the neuron weights initialization. In this paper we present a new initialization technique based on a totally connected undirected graph, that report relations among some intersting features of data input. Result of experimental tests, where the proposed algorithm is compared to the original initialization techniques, shows that our technique assures faster learning and better performance in terms of quantization error.

Redefinition of the kilogram, ampere, kelvin and mole: a proposed approach to implementing CIPM recommendation 1 (CI-2005)

The International System of Units (SI) is founded on seven base units, the metre, kilogram, second, ampere, kelvin, mole and candela corresponding to the seven base quantities of length, mass, time, electric current, thermodynamic temperature, amount of substance and luminous intensity. At its 94th meeting in October 2005, the International Committee for Weights and Measures (CIPM) adopted a recommendation on preparative steps towards redefining the kilogram, ampere, kelvin and mole so that these units are linked to exactly known values of fundamental constants. We propose here that these four base units should be given new definitions linking them to exactly defined values of the Planck constant h, elementary charge e, Boltzmann constant k and Avogadro constant NA, respectively. This would mean that six of the seven base units of the SI would be defined in terms of true invariants of nature. In addition, not only would these four fundamental constants have exactly defined values but also the uncertainties of many of the other fundamental constants of physics would be either eliminated or appreciably reduced. In this paper we present the background and discuss the merits of these proposed changes, and we also present possible wordings for the four new definitions. We also suggest a novel way to define the entire SI explicitly using such definitions without making any distinction between base units and derived units. We list a number of key points that should be addressed when the new definitions are adopted by the General Conference on Weights and Measures (CGPM), possibly by the 24th CGPM in 2011, and we discuss the implications of these changes for other aspects of metrology.

An update on the kilogram

We need a new definition for the kilogram! The present definition was sanctioned by the first General Conference on Weights and Measures (CGPM, Conférence Général des Poids et Mesures) in 1889, with a minor revision to the words in 1901, and remains unchanged after 116 years. It is the only base unit of the International System of Units (the SI) that is still defined in terms of a prototype artifact, the International Prototype of the Kilogram (IPK), which is kept in a safe at the International Bureau of Weights and Measures (the BIPM, Bureau International des Poids et Mesures) in Sèvres, near Paris.

Definitions of the units radian, neper, bel and decibel

The definition of coherent derived units in the International System of Units (SI) is reviewed, and the important role of the equations defining physical quantities is emphasized in obtaining coherent derived units. In the case of the dimensionless quantity plane angle, the choice between alternative definitions is considered, leading to a corresponding choice between alternative definitions of the coherent derived unit - the radian, degree or revolution. In this case the General Conference on Weights and Measures (CGPM) has chosen to adopt the definition that leads to the radian as the coherent derived unit in the SI. In the case of the quantity logarithmic decay (or gain), also sometimes called decrement, and sometimes called level, a similar choice of defining equation exists, leading to a corresponding choice for the coherent derived unit - the neper or the bel. In this case the CGPM has not yet made a choice. We argue that for the quantity logarithmic decay the most logical choice of defining equation is linked to that of the radian, and is that which leads to the neper as the corresponding coherent derived unit. This should not prevent us from using the bel and decibel as units of logarithmic decay. However, it is an important part of the SI to establish in a formal sense the equations defining physical quantities, and the corresponding coherent derived units.