A Gaussian-mixture ensemble transform filter

We generalize the popular ensemble Kalman filter to an ensemble transform filter, in which the prior distribution can take the form of a Gaussian mixture or a Gaussian kernel density estimator. The design of the filter is based on a continuous formulation of the Bayesian filter analysis step. We call the new filter algorithm the ensemble Gaussian-mixture filter (EGMF). The EGMF is implemented for three simple test problems (Brownian dynamics in one dimension, Langevin dynamics in two dimensions and the three-dimensional Lorenz-63 model). It is demonstrated that the EGMF is capable of tracking systems with non-Gaussian uni- and multimodal ensemble distributions. Copyright © 2011 Royal Meteorological Society

Compliant Multi-electrode Micro-channel Array for Recording Afferent Nerve Activity

We have fabricated a compliant neural interface to record afferent nerve activity. Stretchable gold electrodes were evaporated on a polydimethylsiloxane (PDMS) substrate and were encapsulated using photo-patternable PDMS. The built-in microstructure of the gold film on PDMS allows the electrodes to twist and flex repeatedly, without loss of electrical conductivity. PDMS microchannels (5mm long, 100μm wide, 100μm deep) were then plasma bonded irreversibly on top of the electrode array to define five parallel-conduit implants. The soft gold microelectrodes have a low impedance of ~200kΩ at the 1kHz frequency range. Teased nerves from the L6 dorsal root of an anaesthetized Sprague Dawley rat were threaded through the microchannels. Acute tripolar recordings of cutaneous activity are demonstrated, from multiple nerve rootlets simultaneously. Confinement of the axons within narrow microchannels allows for reliable recordings of low amplitude afferents. This electrode technology promises exciting applications in neuroprosthetic devices including bladder fullness monitors and peripheral nervous system implants.

Voltage-gated sodium (NaV) channel blockade by plant cannabinoids does not confer anticonvulsant effects per se

Cannabidiol (CBD) is a non-psychoactive, well-tolerated, anticonvulsant plant cannabinoid, although its mechanism(s) of seizure suppression remains unknown. Here, we investigate the effect of CBD and the structurally similar cannabinoid, cannabigerol (CBG), on voltage-gated Na+ (NaV) channels, a common anti-epileptic drug target. CBG’s anticonvulsant potential was also assessed in vivo. CBD effects on NaV channels were investigated using patch-clamp recordings from rat CA1 hippocampal neurons in brain slices, human SH-SY5Y (neuroblastoma) cells and mouse cortical neurons in culture. CBG effects were also assessed in SH-SY5Y cells and mouse cortical neurons. CBD and CBG effects on veratridine-stimulated human recombinant NaV1.1, 1.2 or 1.5 channels were assessed using a membrane potential-sensitive fluorescent dye high-throughput assay. The effect of CBG on pentyleneterazole-induced (PTZ) seizures was assessed in rat. CBD (10M) blocked NaV currents in SH-SY5Y cells, mouse cortical neurons and recombinant cell lines, and affected spike parameters in rat CA1 neurons; CBD also significantly decreased membrane resistance. CBG blocked NaV to a similar degree to CBD in both SH-SY5Y and mouse recordings, but had no effect (50-200mg/kg) on PTZ-induced seizures in rat. CBD and CBG are NaV channel blockers at micromolar concentrations in human and murine neurons and recombinant cells. In contrast to previous reports investigating CBD, CBG had no effect upon PTZ-induced seizures in rat, indicating that NaV blockade per se does not correlate with anticonvulsant effects.

Fast identification algorithms for Gaussian process model

A class identification algorithms is introduced for Gaussian process(GP)models.The fundamental approach is to propose a new kernel function which leads to a covariance matrix with low rank,a property that is consequently exploited for computational efficiency for both model parameter estimation and model predictions.The objective of either maximizing the marginal likelihood or the Kullback–Leibler (K–L) divergence between the estimated output probability density function(pdf)and the true pdf has been used as respective cost functions.For each cost function,an efficient coordinate descent algorithm is proposed to estimate the kernel parameters using a one dimensional derivative free search, and noise variance using a fast gradient descent algorithm. Numerical examples are included to demonstrate the effectiveness of the new identification approaches.

Inhibition of the cardiac Na+ channel Nav1.5 by carbon monoxide

Sub-lethal carbon monoxide (CO) exposure is frequently associated with myocardial arrhythmias and our recent studies have demonstrated that these may be attributable to modulation of cardiac Na+ channels, causing an increase in the late current and an inhibition of the peak current. Using a recombinant expression system, we demonstrate that CO inhibits peak human Nav1.5 current amplitude without activation of the late Na+ current observed in native tissue. Inhibition was associated with a hyperpolarizing shift in the steady-state inactivation properties of the channels and was unaffected by modification of channel gating induced by anemone toxin (rATX-II). Systematic pharmacological assessment indicated that no recognised CO-sensitive intracellular signalling pathways appeared to mediate CO inhibition of Nav1.5. Inhibition was, however, markedly suppressed by inhibition of nitric oxide (NO) formation, but NO donors did not mimic or occlude channel inhibition by CO, indicating that NO alone did not account for the actions of CO. Exposure of cells to dithiothreitol immediately before CO exposure also dramatically reduced the magnitude of current inhibition. Similarly, L-cysteine and N-ethylmaleimide significantly attenuated the inhibition caused by CO. In the presence of DTT and the NO inhibitor L-NAME, the ability of CO to inhibit Nav1.5 was almost fully prevented. Our data indicate that inhibition of peak Na+ current (which can lead to Brugada-syndrome like arrhythmias) occurs via a mechanism distinct from induction of the late current, requires NO formation and is dependent on channel redox state.

Observation impact in data assimilation: the effect of non-Gaussian observation error.

Data assimilation methods which avoid the assumption of Gaussian error statistics are being developed for geoscience applications. We investigate how the relaxation of the Gaussian assumption affects the impact observations have within the assimilation process. The effect of non-Gaussian observation error (described by the likelihood) is compared to previously published work studying the effect of a non-Gaussian prior. The observation impact is measured in three ways: the sensitivity of the analysis to the observations, the mutual information, and the relative entropy. These three measures have all been studied in the case of Gaussian data assimilation and, in this case, have a known analytical form. It is shown that the analysis sensitivity can also be derived analytically when at least one of the prior or likelihood is Gaussian. This derivation shows an interesting asymmetry in the relationship between analysis sensitivity and analysis error covariance when the two different sources of non-Gaussian structure are considered (likelihood vs. prior). This is illustrated for a simple scalar case and used to infer the effect of the non-Gaussian structure on mutual information and relative entropy, which are more natural choices of metric in non-Gaussian data assimilation. It is concluded that approximating non-Gaussian error distributions as Gaussian can give significantly erroneous estimates of observation impact. The degree of the error depends not only on the nature of the non-Gaussian structure, but also on the metric used to measure the observation impact and the source of the non-Gaussian structure.

Gaussian anamorphosis in the analysis step of the EnKF: a joint state-variable/observation approach

The analysis step of the (ensemble) Kalman filter is optimal when (1) the distribution of the background is Gaussian, (2) state variables and observations are related via a linear operator, and (3) the observational error is of additive nature and has Gaussian distribution. When these conditions are largely violated, a pre-processing step known as Gaussian anamorphosis (GA) can be applied. The objective of this procedure is to obtain state variables and observations that better fulfil the Gaussianity conditions in some sense. In this work we analyse GA from a joint perspective, paying attention to the effects of transformations in the joint state variable/observation space. First, we study transformations for state variables and observations that are independent from each other. Then, we introduce a targeted joint transformation with the objective to obtain joint Gaussianity in the transformed space. We focus primarily in the univariate case, and briefly comment on the multivariate one. A key point of this paper is that, when (1)-(3) are violated, using the analysis step of the EnKF will not recover the exact posterior density in spite of any transformations one may perform. These transformations, however, provide approximations of different quality to the Bayesian solution of the problem. Using an example in which the Bayesian posterior can be analytically computed, we assess the quality of the analysis distributions generated after applying the EnKF analysis step in conjunction with different GA options. The value of the targeted joint transformation is particularly clear for the case when the prior is Gaussian, the marginal density for the observations is close to Gaussian, and the likelihood is a Gaussian mixture.

Calibration of structure in a distributed forecasting model for a semiarid flash flood: dynamic surface storage and channel roughness

Flash floods pose a significant danger for life and property. Unfortunately, in arid and semiarid environment the runoff generation shows a complex non-linear behavior with a strong spatial and temporal non-uniformity. As a result, the predictions made by physically-based simulations in semiarid areas are subject to great uncertainty, and a failure in the predictive behavior of existing models is common. Thus better descriptions of physical processes at the watershed scale need to be incorporated into the hydrological model structures. For example, terrain relief has been systematically considered static in flood modelling at the watershed scale. Here, we show that the integrated effect of small distributed relief variations originated through concurrent hydrological processes within a storm event was significant on the watershed scale hydrograph. We model these observations by introducing dynamic formulations of two relief-related parameters at diverse scales: maximum depression storage, and roughness coefficient in channels. In the final (a posteriori) model structure these parameters are allowed to be both time-constant or time-varying. The case under study is a convective storm in a semiarid Mediterranean watershed with ephemeral channels and high agricultural pressures (the Rambla del Albujón watershed; 556 km 2 ), which showed a complex multi-peak response. First, to obtain quasi-sensible simulations in the (a priori) model with time-constant relief-related parameters, a spatially distributed parameterization was strictly required. Second, a generalized likelihood uncertainty estimation (GLUE) inference applied to the improved model structure, and conditioned to observed nested hydrographs, showed that accounting for dynamic relief-related parameters led to improved simulations. The discussion is finally broadened by considering the use of the calibrated model both to analyze the sensitivity of the watershed to storm motion and to attempt the flood forecasting of a stratiform event with highly different behavior.

Sensitivity of a hydraulic model to changes in channel erosion during extreme flooding

Recent research into flood modelling has primarily concentrated on the simulation of inundation flow without considering the influences of channel morphology. River channels are often represented by a simplified geometry that is implicitly assumed to remain unchanged during flood simulations. However, field evidence demonstrates that significant morphological changes can occur during floods to mobilise the boundary sediments. Despite this, the effect of channel morphology on model results has been largely unexplored. To address this issue, the impact of channel cross-section geometry and channel long-profile variability on flood dynamics is examined using an ensemble of a 1D-2D hydraulic model (LISFLOOD-FP) of the 1:2102 year recurrence interval floods in Cockermouth, UK, within an uncertainty framework. A series of hypothetical scenarios of channel morphology were constructed based on a simple velocity based model of critical entrainment. A Monte-Carlo simulation framework was used to quantify the effects of channel morphology together with variations in the channel and floodplain roughness coefficients, grain size characteristics, and critical shear stress on measures of flood inundation. The results showed that the bed elevation modifications generated by the simplistic equations reflected a good approximation of the observed patterns of spatial erosion despite its overestimation of erosion depths. The effect of uncertainty on channel long-profile variability only affected the local flood dynamics and did not significantly affect the friction sensitivity and flood inundation mapping. The results imply that hydraulic models generally do not need to account for within event morphodynamic changes of the type and magnitude modelled, as these have a negligible impact that is smaller than other uncertainties, e.g. boundary conditions. Instead morphodynamic change needs to happen over a series of events to become large enough to change the hydrodynamics of floods in supply limited gravel-bed rivers like the one used in this research.

Channel equalization for indoor lighting communications networks

We consider indoors communications networks using modulated LEDs to transmit the information packets. A generic indoor channel equalization formulation is proposed assuming the existence of both line of sight and diffuse emitters. The proposed approach is of relevance to emergent indoors distributed sensing modalities for which various lighting based network communications protocols are considered.

Oxygen-dependent hydroxylation by Factor Inhibiting HIF (FIH) regulates the TRPV3 ion channel

Factor Inhibiting HIF (FIH) is an oxygen-dependent asparaginyl hydroxylase that regulates the hypoxia-inducible factors (HIFs). Several proteins containing ankyrin repeat domains have been characterised as substrates of FIH, although there is little evidence for a functional consequence of hydroxylation on these substrates. This study demonstrates that the transient receptor potential vanilloid 3 (TRPV3) channel is hydroxylated by FIH on asparagine 242 within the cytoplasmic ankyrin repeat domain. Hypoxia, FIH inhibitors and mutation of asparagine 242 all potentiated TRPV3-mediated current, without altering TRPV3 protein levels, indicating that oxygen-dependent hydroxylation inhibits TRPV3 activity. This novel mechanism of channel regulation by oxygendependent asparaginyl hydroxylation is likely to extend to other ion channels.

Synaptic vesicle glycoprotein 2A modulates vesicular release and calcium channel function at peripheral sympathetic synapses

Synaptic vesicle glycoprotein (SV)2A is a transmembrane protein found in secretory vesicles and is critical for Ca2+-dependent exocytosis in central neurons, although its mechanism of action remains uncertain. Previous studies have proposed, variously, a role of SV2 in the maintenance and formation of the readily releasable pool (RRP) or in the regulation of Ca2+ responsiveness of primed vesicles. Such previous studies have typically used genetic approaches to ablate SV2 levels; here, we used a strategy involving small interference RNA (siRNA) injection to knockdown solely presynaptic SV2A levels in rat superior cervical ganglion (SCG) neuron synapses. Moreover, we investigated the effects of SV2A knockdown on voltage-dependent Ca2+ channel (VDCC) function in SCG neurons. Thus, we extended the studies of SV2A mechanisms by investigating the effects on vesicular transmitter release and VDCC function in peripheral sympathetic neurons. We first demonstrated an siRNA-mediated SV2A knockdown. We showed that this SV2A knockdown markedly affected presynaptic function, causing an attenuated RRP size, increased paired-pulse depression and delayed RRP recovery after stimulus-dependent depletion. We further demonstrated that the SV2A–siRNA-mediated effects on vesicular release were accompanied by a reduction in VDCC current density in isolated SCG neurons. Together, our data showed that SV2A is required for correct transmitter release at sympathetic neurons. Mechanistically, we demonstrated that presynaptic SV2A: (i) acted to direct normal synaptic transmission by maintaining RRP size, (ii) had a facilitatory role in recovery from synaptic depression, and that (iii) SV2A deficits were associated with aberrant Ca2+ current density, which may contribute to the secretory phenotype in sympathetic peripheral neurons.

Estimation of Gaussian process regression model using probability distance measures

A new class of parameter estimation algorithms is introduced for Gaussian process regression (GPR) models. It is shown that the integration of the GPR model with probability distance measures of (i) the integrated square error and (ii) Kullback–Leibler (K–L) divergence are analytically tractable. An efficient coordinate descent algorithm is proposed to iteratively estimate the kernel width using golden section search which includes a fast gradient descent algorithm as an inner loop to estimate the noise variance. Numerical examples are included to demonstrate the effectiveness of the new identification approaches.

Gaussian processes autoencoder for dimensionality reduction

Learning low dimensional manifold from highly nonlinear data of high dimensionality has become increasingly important for discovering intrinsic representation that can be utilized for data visualization and preprocessing. The autoencoder is a powerful dimensionality reduction technique based on minimizing reconstruction error, and it has regained popularity because it has been efficiently used for greedy pretraining of deep neural networks. Compared to Neural Network (NN), the superiority of Gaussian Process (GP) has been shown in model inference, optimization and performance. GP has been successfully applied in nonlinear Dimensionality Reduction (DR) algorithms, such as Gaussian Process Latent Variable Model (GPLVM). In this paper we propose the Gaussian Processes Autoencoder Model (GPAM) for dimensionality reduction by extending the classic NN based autoencoder to GP based autoencoder. More interestingly, the novel model can also be viewed as back constrained GPLVM (BC-GPLVM) where the back constraint smooth function is represented by a GP. Experiments verify the performance of the newly proposed model.

On-line Gaussian mixture density estimator for adaptive minimum bit-error-rate beamforming receivers

We develop an on-line Gaussian mixture density estimator (OGMDE) in the complex-valued domain to facilitate adaptive minimum bit-error-rate (MBER) beamforming receiver for multiple antenna based space-division multiple access systems. Specifically, the novel OGMDE is proposed to adaptively model the probability density function of the beamformer’s output by tracking the incoming data sample by sample. With the aid of the proposed OGMDE, our adaptive beamformer is capable of updating the beamformer’s weights sample by sample to directly minimize the achievable bit error rate (BER). We show that this OGMDE based MBER beamformer outperforms the existing on-line MBER beamformer, known as the least BER beamformer, in terms of both the convergence speed and the achievable BER.