The advection-dispersion equation that is used to model the solute transport in a porous medium is based on the premise that the fluctuating components of the flow velocity, hence the fluxes, due to a porous matrix can be assumed to obey a relationship similar to Fick's law. This introduces phenomenological coefficients which are dependent on the scale of the experiments. This book presents an approach, based on sound theories of stochastic calculus and differential equations, which removes this basic premise. This leads to a multiscale theory with scale independent coefficients. This book illustrates this outcome with available data at different scales, from experimental laboratory scales to regional scales.