This paper presents the development of a numerical model for reactive transport of multicomponent chemicals in unsaturated soils. The model has been developed based on a coupled thermal, hydraulic, chemical, and mechanical (THCM) formulation, and extended by the inclusion of geochemical processes under mixed equilibrium and kinetically controlled reactions in–between the solid, aqueous, and gas phases in soil. This has been achieved by coupling the transport model, COMPASS, with the geochemical model, PHREEQC. Key coupling between the geochemical modelling and the flow of chemicals has been established via the inclusion of porosity modification from mineral precipitation–dissolution reactions and the consequential effects on flow processes. Verification of the developed model is addressed via a series of benchmark simulations with a focus on testing the coupling between the transport model and geochemical model. Good results have been achieved for the verification of the theoretical and numerical implementation of the new developments in the model. A simulation is presented to demonstrate the effects of mineral reactions on porosity evolution and chemical diffusion in a low porosity soil. The model developed is an advanced tool for studying the hydrogeochemical processes in unsaturated soils under variable THCM conditions.