Discrete Spikes Research Articles

Incorporation of geology with rock physics enables subsalt poststack inversion: A case study in the Gulf of Mexico

One of the main goals of seismic inversion is to obtain high-resolution relative and absolute impedance for reservoir properties prediction. We aim to study whether the results from seismic inversion of subsalt data are sufficiently robust for reliable reservoir characterization. Approximately [Formula: see text] of poststack, wide-azimuth, anisotropic (vertical transverse isotropic) wave-equation migration seismic data from 50 Outer Continental Shelf blocks in the Green Canyon area of the Gulf of Mexico were inverted in this study. A total of four subsalt wells and four subsalt seismic interpreted horizons were used in the inversion process, and one of the wells was used for a blind test. Our poststack inversion method used an iterative discrete spike inversion method, based on the combination of space-adaptive wavelet processing to invert for relative acoustic impedance. Next, the dips were estimated from seismic data and converted to a horizon-like layer sequence field that was used as one of the inputs into the low-frequency model. The background model was generated by incorporating the well velocities, seismic velocity, seismic interpreted horizons, and the previously derived layer sequence field in the low-frequency model. Then, the relative acoustic impedance volume was scaled by adding the low-frequency model to match the calculated acoustic impedance logs from the wells for absolute acoustic impedance. Finally, the geological information and rock physics data were incorporated into the reservoir properties assessment for sand/shale prediction in two main target reservoirs in the Miocene and Wilcox formations. Overall, the poststack inversion results and the sand/shale prediction showed good ties at the well locations. This was clearly demonstrated in the blind test well. Hence, incorporating rock physics and geology enables poststack inversion in subsalt areas.

Reconstructing Stimuli from the Spike Times of Leaky Integrate and Fire Neurons

Reconstructing stimuli from the spike trains of neurons is an important approach for understanding the neural code. One of the difficulties associated with this task is that signals which are varying continuously in time are encoded into sequences of discrete events or spikes. An important problem is to determine how much information about the continuously varying stimulus can be extracted from the time-points at which spikes were observed, especially if these time-points are subject to some sort of randomness. For the special case of spike trains generated by leaky integrate and fire neurons, noise can be introduced by allowing variations in the threshold every time a spike is released. A simple decoding algorithm previously derived for the noiseless case can be extended to the stochastic case, but turns out to be biased. Here, we review a solution to this problem, by presenting a simple yet efficient algorithm which greatly reduces the bias, and therefore leads to better decoding performance in the stochastic case.

The Gross Temporal Correlation of Nearby Neuron Activity in the Macaque Postcentral Somatosensory Cortex Representing Orofacial Structures: With Special Reference to Numerical Methods for Analysis

Evaluation of the temporal correlation between two spike trains is an important issue when studying the dynamics of local neural networks. In the postcentral somatosensory cortex (SI), the temporal profile of neural activity during sustained stimulation is quite variable on the order of 100 msec, reflecting various thalamocortical inputs as well as various stages of the integrating process that might take place within local neural networks. Therefore, a numerical method is needed to evaluate the difference in the gross temporal profile over several hundreds of msec or more. To address this requirement, we tested a method using the dot product between a pair of spike trains. We applied the method to our own data, i. e., 88 pairs of SI neurons recorded simultaneously in one conscious monkey. Each discrete spike train was first smoothed by convolving with a Gaussian function, then the normalized dot product (similarity index) was taken between a pair of spike trains. We next conducted a simulation study using random spike trains (surrogate data set), in which the number of trials and the number of spikes of both units in each trial were set strictly to those of the original data set. The simulation was repeated 100 times for each pair, and the mean and standard deviation of the similarity index were taken. Finally, the difference in the similarity indexes was determined between the original and surrogate data set (z-score). The method successfully quantifies the gross temporal difference of activity in each pair of SI neurons.

Computing the inverse of the neurophysiological spike-response transform

AbstractConsider the transform from a discrete neuronal spike train to a continuous neurophysiological response such as postsynaptic membrane voltage or muscle contraction. Often, we can record only the response. Here we therefore ask about the inverse of this transform: given the response, how can we estimate from it the spike train that produced it? In previous work, we have developed a kernel-based "decoding" method for system identification of the forward transform. Using several variant approaches based on this method, here we demonstrate how to identify the spike train, in the minimal case from just the response, with synthetic as well as real synaptic and neuromuscular data.

Computing the inverse of the neurophysiological spike-response transform

Consider the transform from a discrete neuronal spike train to a continuous neurophysiological response such as postsynaptic membrane voltage or muscle contraction. Often, we can record only the response. Here we therefore ask about the inverse of this transform: given the response, how can we estimate from it the spike train that produced it? In previous work, we have developed a kernel-based "decoding" method for system identification of the forward transform. Using several variant approaches based on this method, here we demonstrate how to identify the spike train, in the minimal case from just the response, with synthetic as well as real synaptic and neuromuscular data.

Computing the inverse of the neurophysiological spike-response transform

Consider the transform from a discrete neuronal spike train to a continuous neurophysiological response such as postsynaptic membrane voltage or muscle contraction. Often, we can record only the response. Here we therefore ask about the inverse of this transform: given the response, how can we estimate from it the spike train that produced it? In previous work, we have developed a kernel-based decoding method for system identification of the forward transform. Using several variant approaches based on this method, here we demonstrate how to identify the spike train, in the minimal case from just the response, with synthetic as well as real synaptic and neuromuscular data.

How precise is the timing of action potentials?

How precise is the timing of action potentials?

Broadband IR supercontinuum generation using single crystal sapphire fibers

In this paper, an investigation on broadband IR supercontinuum generation in single crystal sapphire fibers is presented. It is experimentally demonstrated that broadband IR supercontinuum spectrum (up to 3.2microm) can be achieved by launching ultra-short femtosecond laser pulses into single crystal sapphire fiber with a dimension 115microm in diameter and 5cm in length, which covers both the near IR spectral region and the lower end of the mid-IR spectral range. Furthermore, the mechanism of supercontinuum generation in single crystal sapphire fibers is briefly addressed. When the fiber length is shorter than the dispersion length, the self-phase modulation dominates the broadening effect. In this case, the broad supercontinuum spectrum with a smooth profile can be obtained. However, when the fiber length is longer than the dispersion length, the soliton-related dynamics accompanied by the self-phase modulation dominates the broadening effect. There are discrete spikes in the spectrum (corresponding to different order solitons). The above assumption of supercontinuum generation mechanism is quantitatively modeled by the computer simulation program and verified by the experimental results. Thus, one can adjust the spectral profile by properly choosing the length of the sapphire fibers. The broad IR spectral nature of this supercontinuum source can be very useful in a variety of applications such as broadband LADAR, remote sensing, and multi-spectrum free space communications.

Dexmedetomidine Does Not Reduce Epileptiform Discharges in Adults With Epilepsy

There are limited data on the effect of dexmedetomidine on epileptiform electroencephalogram (EEG). The aim of this study was to investigate if dexmedetomidine will abolish epileptiform discharges in patients with medically refractory seizure disorders who were candidates for surgery to resect foci of epileptic activity. With approval from the Institutional Review Board and written informed consent, we enrolled 5 patients with medically intractable seizures who were undergoing continuous video/EEG monitoring. EEG and hemodynamic values were recorded from 15 minutes before, during, and for 60 minutes after a 60-minute dexmedetomidine infusion. Epileptiform discharges were counted for each 15-minute epoch during the study. Two of the 5 patients had a discrete spike focus in each hemisphere. Thus, we analyzed the activity of 7 distinct foci. Epileptiform activity did not decrease in any individual focus during dexmedetomidine infusion. Although dexmedetomidine did not have a statistically significant effect on interictal epileptiform activity for the group as a whole, the activity of 4 foci increased during dexmedetomidine infusion. Dexmedetomidine did not reduce seizure focus activity and thus may be a suitable anesthetic adjunct during seizure surgery.

Energetic electron precipitation during sawtooth injections

Abstract. We present simultaneous riometer observations of cosmic noise absorption in the nightside and dawn-noon sectors during sawtooth particle injections during 18 April 2002. Energetic electron precipitation (>30 keV) is a feature of magnetospheric substorms and cosmic radio noise absorption acts as a proxy for qualitatively measuring this precipitation. This event provides an opportunity to compare the absorption that accompanies periodic electron injections with the accepted paradigm of substorm-related absorption. We consider whether the absorption is consistent with the premise that these injections are quasi-periodic substorms and study the effects of sustained activity on the level of precipitation. Four consecutive electron injection events have been identified from the LANL (Los Alamos National Laboratory) geosynchronous data; the first two showing that additional activity can occur within the 2–4 h sawtooth periodicity. The first three events have accompanying absorption on the nightside that demonstrate good agreement with the expected pattern of substorm-absorption: discrete spike events with poleward motion at the onset followed by equatorward moving structures and more diffuse absorption, correlated with optical observations. Dayside absorption is linked to gradient-curvature drifting electrons observed at geostationary orbit and it is shown that low fluxes can lead to a lack of absorption as precipitation is suppressed; precipitation begins when the drifting electron flux surpasses some critical level following continuous injections of electrons from the magnetotail. In addition it is shown that the apparent motion of absorption determined from an azimuthal chain of riometers exhibits an acceleration that may be indicative of an energisation of the drifting electron population.

Clinical and electrophysiological characteristics of fasciculoventricular fibers in children

Introduction Fasciculoventricular (FV) fiber is a rare cause for ventricular preexcitation. It is usually described as an innocent bystander pathway. There is only limited data on FV fiber in children. Hence we evaluated the clinical and electrophysiological features of FV fiber in a group of 11 children. Methods and results Of 215 children with manifest preexcitation who had electrophysiological studies at the University of Miami, 11 (5.1%) had characteristics of FV fiber. FV fiber was not directly responsible for any arrhythmias in these children. Three children had supraventricular tachycardia due to associated left sided Kent fiber and FV fiber was identified after the ablation of Kent fiber. One child had associated hypertrophic cardiomyopathy, another child had atrial septal defect and a third child had ventricular septal defect. The electrophysiological testing of FV fiber revealed AH interval of 40–95 ms and H-delta interval of 15–40 ms. Mapping study showed that the FV fiber was located on the right side of the heart in all patients: right anteroseptal in 9, right midseptal in 1 and it could not be mapped well in 1 as the delta wave was intermittent. We identified a discrete FV fiber depolarization spike in 9 patients, with a local FV fiber depolarization to delta wave interval of 5–20 ms (11 + 4 ms). Conclusion FV fiber is an uncommon cause for ventricular preexcitation in children and it can be associated with other conduction abnormalities and heart defects. FV fibers are commonly located at the right anteroseptal region and are not usually involved in any tachycardia. Discrete local FV fiber potential can be identified during electrophysiological testing of these patients.

Analysis of various spike-trains from a digital spiking neuron

A digital spiking neuron (DSN), a quantized analog spiking neuron (QSN) and a discrete spike position map are introduced, and their relations are discussed. Spike position modulations for various spike-trains of DSN and QSN are considered. In addition, some applications are discussed.

Relationship of pore water freshening to accretionary processes in the Cascadia margin: Fluid sources and gas hydrate abundance

Drilling in the Cascadia accretionary complex enable us to evaluate the contribution of dehydration reactions and gas hydrate dissociation to pore water freshening. The observed freshening with depth and distance from the prism toe is consistent with enhanced conversion of smectite to illite, driven by increase in temperature and age of accreted sediments. Although they contain gas hydrate ‐as evidenced by discrete low chloride spikes‐ the westernmost sites drilled on Hydrate Ridge show no freshening trend with depth. Strontium data reveal that all the mélange samples contain deep fluids modified by reaction with the subducting oceanic crust. Thus we infer that, at the westernmost sites, accretion is too recent for the sediments to have undergone significant illitization. Our data demonstrate that a smooth decrease in dissolved chloride with depth cannot generally be used to infer the presence or to estimate the amount of gas hydrate in accretionary margins.

Epilepsia partialis continua—a clinical and electroencephalography study

Epilepsia partialis continua (EPC) is a rare type of localization-related motor epilepsy. Clinical spectrum, electroencephalography (EEG) characteristics and various prognostic factors in EPC were studied in 20 patients. Patients who fulfilled the criteria for EPC between the years 1985 and 1999 were included in this retrospective and prospective study. The mean age was 18 years (range 5 months–70 years). Eleven patients (55%) had Type 1 EPC and in the remaining nine (45%) patients there were features of Type 2 EPC. Among children Rasmussen’s encephalitis and viral encephalitis were the commonest cause for EPC. Encephalitis and vascular aetiology were frequently observed in adults. Tuberculous meningitis and tuberculomas occurred evenly in both the groups. The cause was unknown in two cases. Focal EEG abnormalities commonly consisted of discrete spikes, sharp waves (or) slow wave activity and periodic lateralized epileptiform discharges. The mean duration of follow up was 9.6 months with a range between 1 month and 4 years. Cognitive decline, motor deficits and pharmacoresistance to drugs were significantly seen among children with Type 2 EPC. Patients with Type 1 EPC had mild impairment of functional status with good response to treatment. The long-term prognosis depends upon the underlying cause.

The impact of spike timing variability on the signal-encoding performance of neural spiking models.

It remains unclear whether the variability of neuronal spike trains in vivo arises due to biological noise sources or represents highly precise encoding of temporally varying synaptic input signals. Determining the variability of spike timing can provide fundamental insights into the nature of strategies used in the brain to represent and transmit information in the form of discrete spike trains. In this study, we employ a signal estimation paradigm to determine how variability in spike timing affects encoding of random time-varying signals. We assess this for two types of spiking models: an integrate-and-fire model with random threshold and a more biophysically realistic stochastic ion channel model. Using the coding fraction and mutual information as information-theoretic measures, we quantify the efficacy of optimal linear decoding of random inputs from the model outputs and study the relationship between efficacy and variability in the output spike train. Our findings suggest that variability does not necessarily hinder signal decoding for the biophysically plausible encoders examined and that the functional role of spiking variability depends intimately on the nature of the encoder and the signal processing task; variability can either enhance or impede decoding performance.

Effect of sequence length, sequence frequency, and data acquisition rate on the performance of a hadamard transform time-of-flight mass spectrometer

Various factors influencing the performance of a Hadamard transform time-of-flight mass spectrometer (HT-TOFMS) have been investigated. Using a nitrogen corona discharge to produce an ion stream of N 2 +, N 3 +, and N 4 +, it is found for spectra containing only N 4 + that the signal-to-noise ratio (SNR) closely approaches the value calculated from the ion background by assuming that the ion background follows a Poisson distribution. In contrast, for a more intense beam containing N 2 +, N 3 +, and N 4 +, the SNR is less than its theoretical value because of the appearance of discrete spikes in the mass spectrum caused by deviations in the actual modulation sequence from the ideal one. These spikes can be reduced, however, by decreasing the modulation voltage. Under these optimized conditions, the pseudo-random sequence length is varied to understand how it alters SNR, mass resolution, and scan speed. When the length of the pseudo-random sequence is doubled, the SNR increases by 2 while the time necessary to record a mass spectrum also doubles. Mass resolution can be varied between 500 and 1200 at m/ z = 609 as the sequence length, modulation speed (10 MHz, 25 MHz), and acquisition rate (up to 50 MHz) are changed. Scan speeds of 6000 passes per s can be obtained using a sequence containing 4095 elements modulated at 25 MHz. The capability to tailor the HT-TOFMS to increase the scan speed and resolution with a constant 50% duty cycle makes the technique extremely appealing as a mass analyzer for measuring rapid changes in the composition of an ion stream.

Information theoretic analysis of pulmonary stretch receptor spike trains.

Primary afferent neurons transduce physical, continuous stimuli into discrete spike trains. Investigators have long been interested in interpreting the meaning of the number or pattern of action potentials in attempts to decode the spike train back into stimulus parameters. Pulmonary stretch receptors (PSRs) are visceral mechanoreceptors that respond to deformation of the lungs and pulmonary tree. They provide the brain stem with feedback that is used by cardiorespiratory control circuits. In anesthetized, paralyzed, artificially ventilated rabbits, we recorded the action potential trains of individual PSRs while continuously manipulating ventilator rate and volume. We describe an information theoretic-based analytical method for evaluating continuous stimulus and spike train data that is of general applicability to any continuous, dynamic system. After adjusting spike times for conduction velocity, we used a sliding window to discretize the stimulus (average tracheal pressure) and response (number of spikes), and constructed co-occurrence matrices. We systematically varied the number of categories into which the stimulus and response were evenly divided at 26 different sliding window widths (5, 10, 20, 30,..., 230, 240, 250 ms). Using the probability distributions defined by the co-occurrence matrices, we estimated associated stimulus, response, joint, and conditional entropies, from which we calculated information transmitted as a fraction of the maximum possible, as well as encoding and decoding efficiencies. We found that, in general, information increases rapidly as the sliding window width increases from 5 to approximately 50 ms and then saturates as observation time increases. In addition, the information measures suggest that individual PSRs transmit more "when" than "what" type of information about the stimulus, based on the finding that the maximum information at a given window width was obtained when the stimulus was divided into just a few (usually <6) categories. Our results indicate that PSRs provide quite reliable information about tracheal pressure, with each PSR conveying about 31% of the maximum possible information about the dynamic stimulus, given our analytical parameters. When the stimulus and response are divided into more categories, slightly less information is transmitted, and this quantity also saturates as a function of observation time. We consider and discuss the importance of information contained in window widths on the time scales of an excitatory postsynaptic potential and Hering-Breuer reflex central delay.

Biphasic kinetics of activation and signaling for PAR1 and PAR4 thrombin receptors in platelets.

Thrombin activates platelets in an ordered sequence of events that includes shape change, increase in cytoplasmic Ca(2+), activation of the alphaIIbbeta3 integrin, granule secretion, aggregation, and formation of a stable hemostatic plug. Activation of this process has also been implicated in the pathogenesis of atherosclerosis, stroke, and thrombosis. There are two identified thrombin-activated receptors on the surface of human platelets. PAR1 is a high-affinity thrombin receptor, and PAR4 is a low apparent affinity thrombin receptor of uncertain function. The goal of these studies is to determine the kinetics of thrombin activation of PAR1 and PAR4 and to relate the individual inputs from each receptor to platelet Ca(2+) signaling, secondary autocrine stimulation, and aggregation. Using a combination of PAR-specific peptide ligands and anti-PAR1 reagents, we separated the biphasic thrombin Ca(2+) response of platelets into two discrete components-a rapid spike response caused by PAR1, followed by a slower prolonged response from PAR4. Despite having a 20-70-fold slower rate of activation, PAR4 produces the majority of the integrated Ca(2+) signal that is sustained by the continuous presence of catalytically active thrombin. Surprisingly, PAR4 activation is much more effective than PAR1 activation in mounting secondary autocrine Ca(2+) signals from secreted ADP. The strong ADP response due to activated PAR4, however, requires prior activation of PAR1 as would normally occur during treatment of platelets with thrombin. Thus, the late signal generated by activated PAR4 is not redundant with the early signal from PAR1 and instead serves to greatly extend the high intracellular Ca(2+) levels that support the late phase of the platelet aggregation process.

In-situ Growth and Growth Kinetics of Epitaxial (100) CoSi2 Layer on (100) Si by Reactive Chemical Vapor Deposition

ABSTRACTUniform epitaxial CoSi2 layers have been grown in situ on a (100) Si substrate at temperatures near 600 °C by reactive chemical-vapor deposition of cyclopentadienyl dicarbonyl cobalt, Co(η5-C5H5)(CO)2. The growth kinetics of an epitaxial CoSi2 layer on a Si (100) substrate was investigated at temperatures ranging from 575 to 650 °C. In initial deposition stage, platelike discrete CoSi2 spikes were nucleated along the &lt;111&gt; directions in (100) Si substrate with a twinned structure. The discrete CoSi2 plates with both {111} and (100) planes grew into an epitaxial layer with a flat interface on (100) Si. For epitaxial CoSi2 growth on (100) Si, the activation energy of the parabolic growth was found to be 2.80 eV. The growth rate seems to be controlled by the diffusion of Co through the CoSi2 layer.

Complex effusive events at Kīlauea as documented by the GOES satellite and remote video cameras

GOES provides thermal data for all of the Hawaiian volcanoes once every 15 min. We show how volcanic radiance time series produced from this data stream can be used as a simple measure of effusive activity. Two types of radiance trends in these time series can be used to monitor effusive activity: (a) Gradual variations in radiance reveal steady flow-field extension and tube development. (b) Discrete spikes correlate with short bursts of activity, such as lava fountaining or lava-lake overflows. We are confident that any effusive event covering more than 10,000 m2 of ground in less than 60 min will be unambiguously detectable using this approach. We demonstrate this capability using GOES, video camera and ground-based observational data for the current eruption of Kīlauea volcano (Hawai'i). A GOES radiance time series was constructed from 3987 images between 19 June and 12 August 1997. This time series displayed 24 radiance spikes elevated more than two standard deviations above the mean; 19 of these are correlated with video-recorded short-burst effusive events. Less ambiguous events are interpreted, assessed and related to specific volcanic events by simultaneous use of permanently recording video camera data and ground-observer reports. The GOES radiance time series are automatically processed on data reception and made available in near-real-time, so such time series can contribute to three main monitoring functions: (a) automatically alerting major effusive events; (b) event confirmation and assessment; and (c) establishing effusive event chronology.