Transport of an inert solute in a heterogeneous aquifer is governed by two mechanisms: advection by the random velocity field V(x) and pore-scale dispersion of coefficients Ddij. The velocity field is assumed to be stationary and of constant mean U and of correlation scale I much larger than the pore-scale d. It is assumed that Ddij=&agr;dijU are constant. The relative effect of the two mechanisms is quantified by the Peclet numbers Peij=UI/Ddij=I/&agr;dij, which as a rule are much larger than unity. The main aim of the study is to determine the impact of finite, though high, Pe on ⟨C⟩ and &sgr;C2, the concentration mean and variance, respectively. The solution, derived in the past, for Pe=∞ is reconsidered first. By assuming a normal X probability density function (p.d.f.), closed form solutions are obtained for ⟨C⟩ and &sgr;C2. Recasting the problem in an Eulerian framework leads to the same results if certain closure conditions are adopted. The concentration moments for a finite Pe are derived subsequently in a Lagrangean framework. The pore-scale dispersion is viewed as a Brownian motion type of displacement Xd of solute subparticles, of scale smaller than d, added to the advective displacements X. By adopting again a normal p.d.f. for the latter, explicit expressions for ⟨C⟩ and &sgr;C2 are obtained in terms of quadratures over the joint p.d.f. of advective two particles trajectories. While the influence of high Pe on ⟨C⟩ is generally small, it has a significant impact on &sgr;C2. Simple results are obtained for a small V0, for which trajectories are fully correlated. In particular, the concentration coefficient of variation at the center tends to a constant value for large time. Comparison of the present solution, obtained in terms of a quadrature, with the Monte Carlo simulations of Graham and McLaughlin <1989> shows a very good agreement.¿ 1997 American Geophysical Union