The geometry of random drift I. Stochastic distance and diffusion

1977 ◽  
Vol 9 (2) ◽  
pp. 238-249 ◽  
Author(s):  
Peter L. Antonelli ◽  
Curtis Strobeck

A stochastic distance measure is defined for a general diffusion process on a parameter space X. This distance is defined by where (gij) is the inverse of the covariance matrix of the diffusion equation. This permits the study of the geometry associated with a diffusion equation, since the matrix (gij) is the fundamental tensor of the Riemannian space (X, gij), and of a diffusion process in terms of Brownian motion. For the diffusion equation approximation to random drift with n alleles the covariance matrix is that of a multinomial distribution. The resulting stochastic distance is equal to twice the genetic distance as defined by Cavalli-Sforza and Edwards and is a generalization of the angular transformation of Fisher to n alleles. The geometry associated with the diffusion equation for random drift with n alleles is that of a part of an (n − 1)-sphere of radius two. We also show that the diffusion equation for random drift is not spherical Brownian motion, although it is approximated by it near the centroid of frequency space.

1977 ◽  
Vol 9 (02) ◽  
pp. 238-249 ◽  
Author(s):  
Peter L. Antonelli ◽  
Curtis Strobeck

A stochastic distance measure is defined for a general diffusion process on a parameter space X. This distance is defined by where (gij ) is the inverse of the covariance matrix of the diffusion equation. This permits the study of the geometry associated with a diffusion equation, since the matrix (gij ) is the fundamental tensor of the Riemannian space (X, gij ), and of a diffusion process in terms of Brownian motion. For the diffusion equation approximation to random drift with n alleles the covariance matrix is that of a multinomial distribution. The resulting stochastic distance is equal to twice the genetic distance as defined by Cavalli-Sforza and Edwards and is a generalization of the angular transformation of Fisher to n alleles. The geometry associated with the diffusion equation for random drift with n alleles is that of a part of an (n − 1)-sphere of radius two. We also show that the diffusion equation for random drift is not spherical Brownian motion, although it is approximated by it near the centroid of frequency space.


1978 ◽  
Vol 10 (03) ◽  
pp. 563-569 ◽  
Author(s):  
Peter L. Antonelli

In the present paper a spatially homogeneous distance measure of Edwards and Cavalli-Sforza type is derived for multiple allele random genetic drift with a (possibly vanishing) symmetric mutation field, using the technique of random time substitution. Since the mutation field is of gradient type in gene frequency space, the transformed process is proved to be Brownian motion relative to a new Riemannian geometry and a new time measure. The new geometry is conformally related to spherical geometry of the original process but is not of constant curvature, generally. A formula relating the stationary density of the old and new process is derived and the edge length formula for the new geometry on the n-simplex frequency space is given and analysed.


1978 ◽  
Vol 10 (3) ◽  
pp. 563-569 ◽  
Author(s):  
Peter L. Antonelli

In the present paper a spatially homogeneous distance measure of Edwards and Cavalli-Sforza type is derived for multiple allele random genetic drift with a (possibly vanishing) symmetric mutation field, using the technique of random time substitution. Since the mutation field is of gradient type in gene frequency space, the transformed process is proved to be Brownian motion relative to a new Riemannian geometry and a new time measure. The new geometry is conformally related to spherical geometry of the original process but is not of constant curvature, generally. A formula relating the stationary density of the old and new process is derived and the edge length formula for the new geometry on the n-simplex frequency space is given and analysed.


Author(s):  
ILYA V. TELYATNIKOV

We consider surface measures on the set of trajectories in a smooth compact Riemannian submanifold of Euclidean space generated by diffusion processes in the ambient space. A construction of surface measures on the path space of a smooth compact Riemannian submanifold of Euclidean space was introduced by Smolyanov and Weizsäcker for the case of the standard Brownian motion. The result presented in this paper extends the result of Smolyanov and Weizsäcker to the case when we consider measures generated by diffusion processes in the ambient space with nonidentical correlation operators. For every partition of the time interval, we consider the marginal distribution of the diffusion process in the ambient space under the condition that it visits the manifold at all times of the partition, when the mesh of the partition tends to zero. We prove the existence of some limit surface measures and the equivalence of the above measures to the distribution of some diffusion process on the manifold.


Author(s):  
Tianyu Ma ◽  
Vladimir S. Matveev ◽  
Ilya Pavlyukevich

AbstractWe show that geodesic random walks on a complete Finsler manifold of bounded geometry converge to a diffusion process which is, up to a drift, the Brownian motion corresponding to a Riemannian metric.


2018 ◽  
Vol 146 (12) ◽  
pp. 3949-3976 ◽  
Author(s):  
Herschel L. Mitchell ◽  
P. L. Houtekamer ◽  
Sylvain Heilliette

Abstract A column EnKF, based on the Canadian global EnKF and using the RTTOV radiative transfer (RT) model, is employed to investigate issues relating to the EnKF assimilation of Advanced Microwave Sounding Unit-A (AMSU-A) radiance measurements. Experiments are performed with large and small ensembles, with and without localization. Three different descriptions of background temperature error are considered: 1) using analytical vertical modes and hypothetical spectra, 2) using the vertical modes and spectrum of a covariance matrix obtained from the global EnKF after 2 weeks of cycling, and 3) using the vertical modes and spectrum of the static background error covariance matrix employed to initiate a global data assimilation cycle. It is found that the EnKF performs well in some of the experiments with background error description 1, and yields modest error reductions with background error description 3. However, the EnKF is virtually unable to reduce the background error (even when using a large ensemble) with background error description 2. To analyze these results, the different background error descriptions are viewed through the prism of the RT model by comparing the trace of the matrix , where is the RT model and is the background error covariance matrix. Indeed, this comparison is found to explain the difference in the results obtained, which relates to the degree to which deep modes are, or are not, present in the different background error covariances. The results suggest that, after 2 weeks of cycling, the global EnKF has virtually eliminated all background error structures that can be “seen” by the AMSU-A radiances.


Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. S101-S117 ◽  
Author(s):  
Alba Ordoñez ◽  
Walter Söllner ◽  
Tilman Klüver ◽  
Leiv J. Gelius

Several studies have shown the benefits of including multiple reflections together with primaries in the structural imaging of subsurface reflectors. However, to characterize the reflector properties, there is a need to compensate for propagation effects due to multiple scattering and to properly combine the information from primaries and all orders of multiples. From this perspective and based on the wave equation and Rayleigh’s reciprocity theorem, recent works have suggested computing the subsurface image from the Green’s function reflection response (or reflectivity) by inverting a Fredholm integral equation in the frequency-space domain. By following Claerbout’s imaging principle and assuming locally reacting media, the integral equation may be reduced to a trace-by-trace deconvolution imaging condition. For a complex overburden and considering that the structure of the subsurface is angle-dependent, this trace-by-trace deconvolution does not properly solve the Fredholm integral equation. We have inverted for the subsurface reflectivity by solving the matrix version of the Fredholm integral equation at every subsurface level, based on a multidimensional deconvolution of the receiver wavefields with the source wavefields. The total upgoing pressure and the total filtered downgoing vertical velocity were used as receiver and source wavefields, respectively. By selecting appropriate subsets of the inverted reflectivity matrix and by performing an inverse Fourier transform over the frequencies, the process allowed us to obtain wavefields corresponding to virtual sources and receivers located in the subsurface, at a given level. The method has been applied on two synthetic examples showing that the computed reflectivity wavefields are free of propagation effects from the overburden and thus are suited to extract information of the image point location in the angular and spatial domains. To get the computational cost down, our approach is target-oriented; i.e., the reflectivity may only be computed in the area of most interest.


2021 ◽  
pp. 2140011
Author(s):  
Tomás Caraballo ◽  
Tran Bao Ngoc ◽  
Tran Ngoc Thach ◽  
Nguyen Huy Tuan

This paper is concerned with the mathematical analysis of terminal value problems (TVP) for a stochastic nonclassical diffusion equation, where the source is assumed to be driven by classical and fractional Brownian motions (fBms). Our two problems are to study in the sense of well-posedness and ill-posedness meanings. Here, a TVP is a problem of determining the statistical properties of the initial data from the final time data. In the case [Formula: see text], where [Formula: see text] is the fractional order of a Laplace operator, we show that these are well-posed under certain assumptions. We state a definition of ill-posedness and obtain the ill-posedness results for the problems when [Formula: see text]. The major analysis tools in this paper are based on properties of stochastic integrals with respect to the fBm.


2013 ◽  
Vol 50 (1) ◽  
pp. 29-41
Author(s):  
Alexandra Chronopoulou ◽  
Georgios Fellouris

The problem of detecting an abrupt change in the distribution of an arbitrary, sequentially observed, continuous-path stochastic process is considered and the optimality of the CUSUM test is established with respect to a modified version of Lorden's criterion. We apply this result to the case that a random drift emerges in a fractional Brownian motion and we show that the CUSUM test optimizes Lorden's original criterion when a fractional Brownian motion with Hurst index H adopts a polynomial drift term with exponent H+1/2.


Sign in / Sign up

Export Citation Format

Share Document