We consider multi-component reactive transport in heterogeneous porous media with uncertain hydraulic and chemical properties. This parametric uncertainty is quantified by treating relevant flow and transport parameters as random fields, which renders the governing equations stochastic. We adopt a stochastic Lagrangian framework to replace a three-dimensional advection–reaction transport equation with a one-dimensional equation for solute travel times. We derive approximate expressions for breakthrough curves and their temporal moments. To illustrate our general theory, we consider advective transport of dissolved species undergoing an irreversible bimolecular reaction.