Three kinetic models, for predicting the removal of nitrogen and organics in vertical flow wetlands, have been developed and evaluated. These models were established by combining first-order, Monod and multiple Monod kinetics with continuous stirred-tank reactor (CSTR) flow pattern. Critical evaluations of these models using three statistical parameters, coefficient of determination, relative root mean square error and model efficiency, indicated that when the Monod/multiple Monod kinetics was combined with CSTR flow pattern it allowed close match between theoretical prediction and experiment data of nitrogen and organics removal. The kinetic coefficients (derived from Monod/multiple Monod kinetics) was found to increase with pollutant loading, indicating that the coefficients may vary based on different factors, such as influent pollutant concentration, hydraulic loading, and water depth. Overall, this study demonstrated the validity of combining Monod and multiple Monod kinetics with CSTR flow pattern for the modelling and design of vertical flow wetland systems. (C) 2010 Elsevier Ltd. All rights reserved.