To evaluate the treatment performance of vertical flow sand filters, the HYDRUS wetland module was used to simulate treatment in an experimental set-up. The laboratory filters were intermittently dosed with artificial greywater at a hydraulic loading rate of 0.032 m³ m⁻² day⁻¹ and an organic loading rate of 0.014 kg BOD5 m⁻² day⁻¹. The hydraulic properties of the filter were characterised, as were inflow and outflow concentrations of chemical oxygen demand (COD), biochemical oxygen demand (BOD), ammonia, nitrate and total nitrogen. The inverse simulation function of the HYDRUS software was used to calibrate the water flow model. The observed effect of water flowing faster along the column wall was included in the inverse simulations. The biokinetic model was calibrated by fitting heterotrophic biomass growth to measurements of potential respiration rate. Emphasis was put on simulating outflow concentrations of organic pollutants. The simulations were conducted using three models of varying degree of calibration effort and output accuracy. The effluent concentration was 245 mg COD L⁻¹ for the laboratory filters, 134 mg COD mg L⁻¹ for the model excluding wall flow effects and 338 mg COD mg L⁻¹ for the model including wall flow effects.