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.

Compositional control of the electrical transport properties in the new series of thiospinels Ga1-xGexV4S8-delta (0 <=x<=1)

A new series of non-stoichiometric sulfides Ga1−xGexV4S8−δ (0≤x≤1; δ≤0.23) has been synthesized at high temperatures by heating stoichiometric mixtures of the elements in sealed quartz tubes. The samples have been characterized by powder X-ray diffraction, SQUID magnetometry and electrical transport-property measurements. Structural analysis reveals that a solid solution is formed throughout this composition range, whilst thermogravimetric data reveal sulfur deficiency of up to 2.9% in the quaternary phases. Magnetic measurements suggest that the ferromagnetic behavior of the end-member phase GaV4S8 is retained at x≤0.7; samples in this composition range showing a marked increase in magnetization at low temperatures. By contrast Ga0.25Ge0.75V4S8−δ appears to undergo antiferromagnetic ordering at ca. 15 K. All materials with x≠1 are n-type semiconductors whose resistivity falls by almost six orders of magnitude with decreasing Ga content, whilst the end-member phase GeV4S8−δ is a p-type semiconductor. The results demonstrate that the physical properties are determined principally by the degree of electron filling of narrow-band states arising from intracluster V–V interactions.

Structure and electrical transport properties of the ordered skutterudites MGe1.5S1.5 (M = Co, Rh, Ir)

High-resolution powder neutron diffraction data collected for the skutterudites MGe1.5S1.5 (M=Co, Rh, Ir) reveal that these materials adopt an ordered skutterudite structure (space group R3¯), in which the anions are ordered in layers perpendicular to the [111] direction. In this ordered structure, the anions form two-crystallographically distinct four-membered rings, with stoichiometry Ge2S2, in which the Ge and S atoms are trans to each other. The transport properties of these materials, which are p-type semiconductors, are discussed in the light of the structural results. The effect of iron substitution in CoGe1.5S1.5 has been investigated. While doping of CoGe1.5S1.5 has a marked effect on both the electrical resistivity and the Seebeck coefficient, these ternary skutterudites exhibit significantly higher electrical resistivities than their binary counterparts.

Time-lapse geophysical investigations over known archaeological features using electrical resistivity imaging and earth resistance

Electrical methods of geophysical survey are known to produce results that are hard to predict at different times of the year, and under differing weather conditions. This is a problem which can lead to misinterpretation of archaeological features under investigation. The dynamic relationship between a ‘natural’ soil matrix and an archaeological feature is a complex one, which greatly affects the success of the feature’s detection when using active electrical methods of geophysical survey. This study has monitored the gradual variation of measured resistivity over a selection of study areas. By targeting difficult to find, and often ‘missing’ electrical anomalies of known archaeological features, this study has increased the understanding of both the detection and interpretation capabilities of such geophysical surveys. A 16 month time-lapse study over 4 archaeological features has taken place to investigate the aforementioned detection problem across different soils and environments. In addition to the commonly used Twin-Probe earth resistance survey, electrical resistivity imaging (ERI) and quadrature electro-magnetic induction (EMI) were also utilised to explore the problem. Statistical analyses have provided a novel interpretation, which has yielded new insights into how the detection of archaeological features is influenced by the relationship between the target feature and the surrounding ‘natural’ soils. The study has highlighted both the complexity and previous misconceptions around the predictability of the electrical methods. The analysis has confirmed that each site provides an individual and nuanced situation, the variation clearly relating to the composition of the soils (particularly pore size) and the local weather history. The wide range of reasons behind survey success at each specific study site has been revealed. The outcomes have shown that a simplistic model of seasonality is not universally applicable to the electrical detection of archaeological features. This has led to the development of a method for quantifying survey success, enabling a deeper understanding of the unique way in which each site is affected by the interaction of local environmental and geological conditions.

The oscillatory bioelectrical signal from plants explained by a simulated electrical model and tested using Lempel-Ziv complexity

This paper demonstrates the oscillatory characteristics of electrical signals acquired from two ornamental plant types (Epipremnum pinnatum and Philodendron scandens - Family Araceae), using a noninvasive acquisition system. The electrical signal was recorded using Ag/AgCl superficial electrodes inside a Faraday cage. The presence of the oscillatory electric generator was shown using a classical power spectral density. The Lempel and Ziv complexity measurement showed that the plant signal was not noise despite its nonlinear behavior. The oscillatory characteristics of the signal were explained using a simulated electrical model that establishes that for a frequency range from 5 to 15 Hz, the oscillatory characteristic is higher than for other frequency ranges. All results show that non-invasive electrical plant signals can be acquired with improvement of signal-to-noise ratio using a Faraday cage, and a simple electrical model is able to explain the electrical signal being generated. (C) 2010 Elsevier B.V. All rights reserved.

Evaluation of the pH and electrical conductivity in milk from jersey cows during the first month of lactation

In order to assess the influence of the colostrum period on pH and, electrical conductivity, we collected 418 milk samples from 127 Jersey cows. The samples were collected from healthy udders that did not present any bacterial growth in the microbiological examination. They were divided into eight groups as follows < 1/2 day; 1/2 and 1 degrees day; 2 degrees day; 3 degrees day; 4 degrees and 5 degrees day; 6 degrees and 7 degrees day; 8 degrees to 15 degrees day; 16 degrees to 30 degrees days of lactation. The samples were collected before milking and the following analyses were conducted: pH, electrical conductivity. In the first 24 hours of lactation, there was an reduction in electrical conductivity value, associated with an increase in pH value. We observed that transition of secretion from colostrum to milk, occurs during the first week of lactation; from 6(rd) day of lactation for pH value and 3(th) day for electrical conductivity value. We recommend the use the following figures as normal ranges for the first 24 hours of lactation (colostrum period): pH <= 6,51 and electrical conductivity <= 6,33 mS/cm; while for the interval between 2(nd) and 7(th) days of lactation (transition from colostrum to milk) we suggest the use of the values as normal ranges: pH <= 6,66 and electrical conductivity <= 5,93 mS/cm.

Use of electrical resistivity to identify collapsible soils in Brazil

Three soil profiles in Ilha Solteira, Brazil were investigated to establish their potential for collapsible behavior. The soil profiles were identified using terrain evaluation techniques and simple laboratory tests. Geophysical surveys were undertaken as they are quick and relatively cheap. The results were correlated with trial pit descriptions and cone and standard penetration tests. The study has shown that electrical resistivity is a useful tool for the preliminary identification of horizons of collapsible soils, before more expensive intrusive and laboratory work is undertaken.

Neuronal Expression of Cd36, Cd44, and Cd83 Antigen Transcripts Maps to Distinct and Specific Murine Brain Circuits

Cells recruited by the innate immune response rely on surface-expressed molecules in order to receive signals from the local environment and to perform phagocytosis, cell adhesion, and others processes linked to host defense. Hundreds of surface antigens designated through a cluster of differentiation (CD) number have been used to identify particular populations of leukocytes. Surprisingly, we verified that the genes that encode Cd36 and Cd83 are constitutively expressed in specific neuronal cells. For instance, Cd36 mRNA is expressed in some regions related to circuitry involved in pheromone responses and reproductive behavior. Cd44 expression, reanalyzed and detailed here, is associated with the laminar formation and midline thalamic nuclei in addition to striatum, extended amygdala, and a few hypothalamic, cortical, and hippocampal regions. A systemic immune challenge was able to increase Cd44 expression quickly in the area postrema and motor nucleus of the vagus but not in regions presenting expressive constitutive expression. In contrast to Cd36 and Cd44, Cd83 message was widely distributed from the olfactory bulb to the brain stem reticular formation, sparing the striatopallidum, olivary region, and cerebellum. Its pattern of expression nevertheless remained strongly associated with hypothalamic, thalamic, and hindbrain nuclei. Unlike the other transcripts, Cd83 mRNA was rapidly modulated by restraint stress. Our results indicate that these molecules might play a role in specific neural circuits and present functions other than those attributed to leukocyte biology. The data also suggest that these surface proteins, or their associated mRNA, could be used to label neurons in specific circuits/regions. J. Comp. Neurol. 517:906-924, 2009. (C) 2009 Wiley-Liss, Inc.

Neuromuscular electrical stimulation improves GLUT-4 and morphological characteristics of skeletal muscle in rats with heart failure

Aim: Changes in skeletal muscle morphology and metabolism are associated with limited functional capacity in heart failure, which can be attenuated by neuromuscular electrical stimulation (ES). The purpose of the present study was to analyse the effects of ES upon GLUT-4 protein content, fibre structure and vessel density of the skeletal muscle in a rat model of HF subsequent to myocardial infarction. Methods: Forty-four male Wistar rats were assigned to one of four groups: sham (S), sham submitted to ES (S+ES), heart failure (HF) and heart failure submitted to ES (HF+ES). The rats in the ES groups were submitted to ES of the left leg during 20 days (2.5 kHz, once a day, 30 min, duty cycle 50%- 15 s contraction/15 s rest). After this period, the left tibialis anterior muscle was collected from all the rats for analysis. Results: HF+ES rats showed lower values of lung congestion when compared with HF rats (P = 0.0001). Although muscle weight was lower in HF rats than in the S group, thus indicating hypotrophy, 20 days of ES led to their recovery (P < 0.0001). In both groups submitted to ES, there was an increase in muscle vessel density (P < 0.04). Additionally, heart failure determined a 49% reduction in GLUT-4 protein content (P < 0.03), which was recovered by ES (P < 0.01). Conclusion: In heart failure, ES improves morphological changes and raises GLUT-4 content in skeletal muscle.

Effect of neuromuscular electrical stimulation and isotonic exercises in flexor and extensor muscles of knee of hemiplegic patients

Aim. To verify the muscular force and resistance to the movement of the flexor and extensor muscles of the knee of patients with spasticity after treatment with neuromuscular electrical stimulation (NMES) and isotonic exercises. Patients and methods. The patients this study were divided into group I (NMES) and group 2 (isotonic exercises). Their muscular torque and resistance to the movement of the flexor and extensor knee muscles were measured by the isokinetic dynamometer and the degree of spasticity by the modified Ashworth scale before and after ten sessions. Results. Alterations in the scores of the modified Ashworth scale were not observed. An increase in the flexor torque in group 1 (p = 0.041) and in group 2 (p = 0.001) was verified. In the passive mode, group 1 presented a reduction of resistance to the flexion movement (p = 0.026), while in group 2, a reduction of resistance to both the flexion (p = 0,029) and extension movements (p = 0.019) was verified. Conclusions. The two therapeutical resources had their efficiency proven only for the increase of the force of the flexor muscles. The resistance to movement, the isotonic exercises were more effective because they promoted a reduction in the resistance of the flexor and extensor knee muscles.

Two-level network design with intermediate facilities: An application to electrical distribution systems

We consider the two-level network design problem with intermediate facilities. This problem consists of designing a minimum cost network respecting some requirements, usually described in terms of the network topology or in terms of a desired flow of commodities between source and destination vertices. Each selected link must receive one of two types of edge facilities and the connection of different edge facilities requires a costly and capacitated vertex facility. We propose a hybrid decomposition approach which heuristically obtains tentative solutions for the vertex facilities number and location and use these solutions to limit the computational burden of a branch-and-cut algorithm. We test our method on instances of the power system secondary distribution network design problem. The results show that the method is efficient both in terms of solution quality and computational times. (C) 2010 Elsevier Ltd. All rights reserved.

A Parallel Hardware Architecture for Scale and Rotation Invariant Feature Detection

This paper proposes a parallel hardware architecture for image feature detection based on the Scale Invariant Feature Transform algorithm and applied to the Simultaneous Localization And Mapping problem. The work also proposes specific hardware optimizations considered fundamental to embed such a robotic control system on-a-chip. The proposed architecture is completely stand-alone; it reads the input data directly from a CMOS image sensor and provides the results via a field-programmable gate array coupled to an embedded processor. The results may either be used directly in an on-chip application or accessed through an Ethernet connection. The system is able to detect features up to 30 frames per second (320 x 240 pixels) and has accuracy similar to a PC-based implementation. The achieved system performance is at least one order of magnitude better than a PC-based solution, a result achieved by investigating the impact of several hardware-orientated optimizations oil performance, area and accuracy.