  In the present study we examine non-Gaussian spreading of solutes subject to
advection, dispersion and kinetic sorption (adsorption/desorption). We start
considering the behavior of a single particle and apply a random walk to
describe advection/dispersion plus a Markov chain to describe kinetic sorption.
We show in a rigorous way that this model leads to a set of differential
equations. For this combination of stochastic processes such a derivation is
new. Then, to illustrate the mechanism that leads to non-Gaussian spreading we
analyze this set of equations at first leaving out the Gaussian dispersion term
(microdispersion). The set of equations now transforms to the telegrapher's
equation. Characteristic for this system is a longitudinal spreading, that
becomes Gaussian only in the long-time limit. We refer to this as kinetics
induced spreading. When the microdispersion process is included back again, the
characteristics of the telegraph equations are still present. Now two spreading
phenomena are active, the Gaussian microdispersive spreading plus the kinetics
induced non-Gaussian spreading. In the long run the latter becomes Gaussian as
well. Another non-Gaussian feature shows itself in the 2D situation. Here, the
lateral spread and the longitudinal displacement are no longer independent, as
should be the case for a 2D Gaussian spreading process. In a displacing plume
this interdependence is displayed as a `tailing' effect. We also analyze
marginal and conditional moments, which confirm this result. With respect to
effective properties (velocity and dispersion) we conclude that effective
parameters can be defined properly only for large times (asymptotic times). In
the two-dimensional case it appears that the transverse spreading depends on
the longitudinal coordinate. This results in `cigar-shaped' contours.
