Groundwater contamination usually originates from surface contamination. Contaminants then move downward through the vadose zone and finally reach the groundwater table. To date, however, analytical solutions of multi-species reactive transport are limited to transport only in the saturated zone. The motivation of this work is to utilize analytical solutions, which were previously derived for single-phase transport, to describe the reactive transport of multiple volatile contaminants in the unsaturated zone. A mathematical model is derived for describing transport with phase partitioning of sequentially reactive species in the vadose zone with constant flow velocity. Linear reaction kinetics and linear equilibrium partitioning between vapor, liquid, and solid phases are assumed in this model.