Integral transforms have been widely used for deriving analytical solutions for solute transport systems. Often, analytical solutions can only be written in closed form in frequency domains and numerical inverse-transforms have to be involved to obtain semi-analytical solutions in the time domain. For this reason, previously published closed form solutions are restricted either to a small number of species or to the same retardation assumption. In this paper, we applied the solution scheme proposed by Bauer et al. in the time domain. Using available analytical solutions of a single species transport with first-order decay without coupling with its parent species concentration as fundamental solutions, a daughter species concentration can be expressed as a linear function of those fundamental solutions. The implementation of the solution scheme is straight forward and exact analytical solutions are derived for one- and three-dimensional transport systems.