Multispecies contaminant transport in the Earth's subsurface is commonly modelled using advection-dispersion equations coupled via first-order reactions. Analytical and semi-analytical solutions for such problems are highly sought after but currently limited to either one species, homogeneous media, certain reaction networks, specific boundary conditions or a combination thereof. In this paper,...