scholarly journals Stochastic slowly adapting ionic currents may provide a decorrelation mechanism for neural oscillators by causing wander in the intrinsic period

2016 ◽  
Vol 116 (3) ◽  
pp. 1189-1198 ◽  
Author(s):  
Sharon E. Norman ◽  
Robert J. Butera ◽  
Carmen C. Canavier

Oscillatory neurons integrate their synaptic inputs in fundamentally different ways than normally quiescent neurons. We show that the oscillation period of invertebrate endogenous pacemaker neurons wanders, producing random fluctuations in the interspike intervals (ISI) on a time scale of seconds to minutes, which decorrelates pairs of neurons in hybrid circuits constructed using the dynamic clamp. The autocorrelation of the ISI sequence remained high for many ISIs, but the autocorrelation of the ΔISI series had on average a single nonzero value, which was negative at a lag of one interval. We reproduced these results using a simple integrate and fire (IF) model with a stochastic population of channels carrying an adaptation current with a stochastic component that was integrated with a slow time scale, suggesting that a similar population of channels underlies the observed wander in the period. Using autoregressive integrated moving average (ARIMA) models, we found that a single integrator and a single moving average with a negative coefficient could simulate both the experimental data and the IF model. Feeding white noise into an integrator with a slow time constant is sufficient to produce the autocorrelation structure of the ISI series. Moreover, the moving average clearly accounted for the autocorrelation structure of the ΔISI series and is biophysically implemented in the IF model using slow stochastic adaptation. The observed autocorrelation structure may be a neural signature of slow stochastic adaptation, and wander generated in this manner may be a general mechanism for limiting episodes of synchronized activity in the nervous system.

Genetics ◽  
2000 ◽  
Vol 154 (3) ◽  
pp. 1403-1417 ◽  
Author(s):  
David J Cutler

Abstract Rates of molecular evolution at some protein-encoding loci are more irregular than expected under a simple neutral model of molecular evolution. This pattern of excessive irregularity in protein substitutions is often called the “overdispersed molecular clock” and is characterized by an index of dispersion, R(T) > 1. Assuming infinite sites, no recombination model of the gene R(T) is given for a general stationary model of molecular evolution. R(T) is shown to be affected by only three things: fluctuations that occur on a very slow time scale, advantageous or deleterious mutations, and interactions between mutations. In the absence of interactions, advantageous mutations are shown to lower R(T); deleterious mutations are shown to raise it. Previously described models for the overdispersed molecular clock are analyzed in terms of this work as are a few very simple new models. A model of deleterious mutations is shown to be sufficient to explain the observed values of R(T). Our current best estimates of R(T) suggest that either most mutations are deleterious or some key population parameter changes on a very slow time scale. No other interpretations seem plausible. Finally, a comment is made on how R(T) might be used to distinguish selective sweeps from background selection.


Hydrology ◽  
2021 ◽  
Vol 8 (2) ◽  
pp. 79
Author(s):  
Tribeni C. Sharma ◽  
Umed S. Panu

On a global basis, there is trend that a majority of reservoirs are sized using a draft of 75% of the mean annual flow (0.75 MAF). The reservoir volumes based on the proposed drought magnitude (DM) method and the sequent peak algorithm (SPA) at 0.75 MAF draft were compared at the annual, monthly and weekly scales using the flow sequences of 25 Canadian rivers. In our assessment, the monthly scale is adequate for such analyses. The DM method, although capable of using flow data at any time scale, has been demonstrated using monthly standardized hydrological index (SHI) sequences. The moving average (MA) smoothing of the monthly SHI sequences formed the basis in the DM method for estimating the reservoir volume through the use of the extreme number theorem, and the hypothesis that drought magnitude is equal to the product of the drought intensity and drought length. The truncation level in the SHI sequences was found as SHIo [ = (0.75 ‒ 1) µo/σo], where µo and σo are the overall mean and standard deviation of the monthly flows. The DM-based estimates for the deficit volumes and the SPA-based reservoir volumes were found comparable within an error margin of ±18%.


1992 ◽  
Vol 36 (01) ◽  
pp. 1-16
Author(s):  
G. A. Athanassoulis ◽  
P. B. Vranas ◽  
T. H. Soukissian

A new approach for calculating the long-term statistics of sea waves is proposed. A rational long-term stochastic model is introduced which recognizes that the wave climate at a given site in the ocean consists of a random succession of individual sea states, each sea state possessing its own duration and intensity. This model treats the sea-surface elevation as a random function of a "fast" time variable, and the time history of the spectral characteristics of the successive sea states as a random function of a "slow" time variable. By developing an appropriate conceptual framework, it becomes possible to express various probabilistic characteristics of the sea-surface elevation, which are sensible only in the fast-time scale, in terms of the statistics of sea-states duration and intensity, which is meaningful only in the slow-time scale. As an example, we study the random quantity MU(T) = "number of maxima of the sea-surface elevation lying above the level u and occurring during a long-term time period [0,T]." Exploiting the proposed framework, it is shown that, under certain clearly defined assumptions, Mu(T) can be given the structure of a renewal-reward (cumulative) process, whose interarrival times correspond to the duration of successive sea states. Thus, using renewal theory, the complete characterization of the probability structure of MU(T) is obtained. As a consequence, the long-term probability distribution function of the individual wave height is rigorously defined and calculated. The relation of the present results with corresponding ones previously obtained is thoroughly discussed. The proposed model can be extended twofold: either by replacing some of the simplifying assumptions by more realistic ones, or by extending the model for treating the corresponding problems for ship and structures responses.


Author(s):  
Anindya Chatterjee ◽  
Joseph P. Cusumano

Abstract We present a new observer-based method for parameter estimation for nonlinear oscillatory mechanical systems where the unknown parameters appear linearly (they may each be multiplied by bounded and Lipschitz continuous but otherwise arbitrary, possibly nonlinear, functions of the oscillatory state variables and time). The oscillations in the system may be periodic, quasiperiodic or chaotic. The method is also applicable to systems where the parameters appear nonlinearly, provided a good initial estimate of the parameter is available. The observer requires measurements of displacements. It estimates velocities on a fast time scale, and the unknown parameters on a slow time scale. The fast and slow time scales are governed by a single small parameter ϵ. Using asymptotic methods including the method of averaging, it is shown that the observer’s estimates of the unknown parameters converge like e−ϵt where t is time, provided the system response is such that the coefficient-functions of the unknown parameters are not close to being linearly dependent. It is also shown that the method is robust in that small errors in the model cause small errors in the parameter estimates. A numerical example is provided to demonstrate the effectiveness of the method.


Radiocarbon ◽  
1989 ◽  
Vol 31 (2) ◽  
pp. 169-177 ◽  
Author(s):  
Ad Stolk ◽  
Koos Hogervorst ◽  
Henk Berendsen

Large numbers of 14C dates of the base and top of Holocene peat layers may be plotted in 14C histograms in order to establish statistically a chronology of periods of essentially clastic sedimentation and peat formation. Due to the non-linearity of the 14C time scale in terms of calendar years, clustering of 14C dates on random peat growth may occur. This seriously hampers the interpretation of histograms. A quantitative method and computer program were developed to correct the histograms for this effect. The correction factor that has to be applied depends on the calibration curve and the interval width of the correction parameter dy. For peat samples, an interval width of 100 14C yr and a calibration curve based on a 100-yr moving average seems to be a reasonable choice.


Plasma ◽  
2018 ◽  
Vol 1 (1) ◽  
pp. 126-143
Author(s):  
Jeffrey Parker ◽  
Lynda LoDestro ◽  
Alejandro Campos

One route to improved predictive modeling of magnetically confined fusion reactors is to couple transport solvers with direct numerical simulations (DNS) of turbulence, rather than with surrogate models. An additional challenge presented by coupling directly with DNS is the inherent fluctuations in the turbulence, which limit the convergence achievable in the transport solver. In this article, we investigate the performance of one numerical coupling method in the presence of turbulent fluctuations. To test a particular numerical coupling method for the transport solver, we use an autoregressive-moving-average model to generate stochastic fluctuations efficiently with statistical properties resembling those of a gyrokinetic simulation. These fluctuations are then added to a simple, solvable problem, and we examine the behavior of the coupling method. We find that monitoring the residual as a proxy for the error can be misleading. From a pragmatic point of view, this study aids us in the full problem of transport coupled to DNS by predicting the amount of averaging required to reduce the fluctuation error and obtain a specific level of accuracy.


2007 ◽  
Vol 40 (7) ◽  
pp. 1629-1634 ◽  
Author(s):  
Jonathan B. Dingwell ◽  
Domenic F. Napolitano ◽  
David Chelidze

2001 ◽  
Vol 09 (04) ◽  
pp. 291-314 ◽  
Author(s):  
STEFAN SCHUSTER ◽  
MARKO MARHL

The behavior of calcium oscillations near bifurcations is analyzed for three different models. For the model developed by Somogyi and Stucki [42], it is shown that the range of oscillations is bounded by supercritical and subcritical Hopf bifurcations. Near the latter, canard orbits arise, that is, quasi-harmonic oscillations with a very small amplitude grow very fast to become pulsed oscillations. The potential biological significance of this behavior is discussed. A time-scale analysis of this model is performed and an approximation formula for the oscillation period is derived. For two models that we presented earlier [30, 31], it is shown that a homoclinic bifurcation and an infinite period bifurcation, respectively, occur. These imply that the oscillation period can reach arbitrarily high values. This behavior is discussed in the light of frequency encoding, and the scaling laws of the oscillation period are given.


Sign in / Sign up

Export Citation Format

Share Document