ABSTRACT The problem of solute transport through the unsaturated zone of heterogeneous porous formations is discussed, focusing on mechanistic, stochastic vadose-zone transport models based on two different approaches to modeling flow in these formations, namely the independent vertical columns and the stochastic continuum approaches. It is shown here that the resulting transport models differ in their predictions with respect to (i) how solute spreading evolves with time; and (ii) how water saturation of the formation affects solute spreading. Predictions of transport models based on the first approach are restricted to solute spreading in the longitudinal direction only. Predictions of transport models based on the second approach are in qualitative agreement with the results of numerical simulations of flow and transport in partially saturated, heterogeneous formations, and are consistent with results obtained by the theory of transport by groundwater flow. The need for additional field-scale characterizations of the spatial variability of the soil properties relevant to vadose-zone transport, and for additional, carefully designed field-scale transport experiments and simulations, to validate further existing stochastic transport models, is emphasized.