Space-time duality for fractional diffusion

Zolotarev proved a duality result that relates stable densities with different indices. In this paper, we show how Zolotarev duality leads to some interesting results on fractional diffusion. Fractional diffusion equations employ fractional derivatives in place of the usual integer order derivatives. They govern scaling limits of random walk models, with power law jumps leading to fractional derivatives in space, and power law waiting times between the jumps leading to fractional derivatives in time. The limit process is a stable L\'evy motion that models the jumps, subordinated to an inverse stable process that models the waiting times. Using duality, we relate the density of a spectrally negative stable process with index $1<\alpha<2$ to the density of the hitting time of a stable subordinator with index $1/\alpha$, and thereby unify some recent results in the literature. These results also provide a concrete interpretation of Zolotarev duality in terms of the fractional diffusion model.Comment: 16 page

A study of atmospheric diffusion from the LANDSAT imagery

LANDSAT multispectral scanner data of the smoke plumes which originated in eastern Cabo Frio, Brazil and crossed over into the Atlantic Ocean, are analyzed to illustrate how high resolution LANDSAT imagery can aid meteorologists in evaluating specific air pollution events. The eleven LANDSAT images selected are for different months and years. The results show that diffusion is governed primarily by water and air temperature differences. With colder water, low level air is very stable and the vertical diffusion is minimal; but water warmer than the air induces vigorous diffusion. The applicability of three empirical methods for determining the horizontal eddy diffusivity coefficient in the Gaussian plume formula was evaluated with the estimated standard deviation of the crosswind distribution of material in the plume from the LANDSAT imagery. The vertical diffusion coefficient in stable conditions is estimated using Weinstock's formulation. These results form a data base for use in the development and validation of meso scale atmospheric diffusion models

Soft and hard wall in a stochastic reaction diffusion equation

We consider a stochastically perturbed reaction diffusion equation in a bounded interval, with boundary conditions imposing the two stable phases at the endpoints. We investigate the asymptotic behavior of the front separating the two stable phases, as the intensity of the noise vanishes and the size of the interval diverges. In particular, we prove that, in a suitable scaling limit, the front evolves according to a one-dimensional diffusion process with a non-linear drift accounting for a "soft" repulsion from the boundary. We finally show how a "hard" repulsion can be obtained by an extra diffusive scaling.Comment: 33 page

Stable L\'{e}vy diffusion and related model fitting

A fractional advection-dispersion equation (fADE) has been advocated for heavy-tailed flows where the usual Brownian diffusion models fail. A stochastic differential equation (SDE) driven by a stable L\'{e}vy process gives a forward equation that matches the space-fractional advection-dispersion equation and thus gives the stochastic framework of particle tracking for heavy-tailed flows. For constant advection and dispersion coefficient functions, the solution to such SDE itself is a stable process and can be derived easily by least square parameter fitting from the observed flow concentration data. However, in a more generalized scenario, a closed form for the solution to a stable SDE may not exist. We propose a numerical method for solving/generating a stable SDE in a general set-up. The method incorporates a discretized finite volume scheme with the characteristic line to solve the fADE or the forward equation for the Markov process that solves the stable SDE. Then we use a numerical scheme to generate the solution to the governing SDE using the fADE solution. Also, often the functional form of the advection or dispersion coefficients are not known for a given plume concentration data to start with. We use a Levenberg--Marquardt (L-M) regularization method to estimate advection and dispersion coefficient function from the observed data (we present the case for a linear advection) and proceed with the SDE solution construction described above.Comment: Published at https://doi.org/10.15559/18-VMSTA106 in the Modern Stochastics: Theory and Applications (https://vmsta.org/) by VTeX (http://www.vtex.lt/