A signi®cant source of uncertainty in multicomponent transport models is due to reaction parameter estimation error. In this paper, two codes are developed to compute reaction parameter sensitivity coecients, which provide a quantitative measure of the impact of this uncertainty. One of the codes solves the full coupled system of sensitivity equations, while the other applies the operator splitting approach to decouple the sensitivity equations for each chemical component. CPU time comparisons demonstrate that operator splitting signi®cantly reduces the simulation time for solving both the state and sensitivity equations. These comparisons also con®rm that solution of the sensitivity equations is more ecient than a traditional direct perturbation sensitivity analysis. Sensitivity calculations for an example transport problem with cobalt and nitrilotriacetate (NTA) are used to gain insight into the relative signi®cance of reaction processes and to rank individual reaction parameters in terms of importance. In this example, speciation reactions cause an indirect linking of the NTA biodegradation and the cobalt sorption reactions, making the NTA concentration nearly as sensitive to the sorption parameters as to the degradation parameters.