Subsurface constructed wetland (SCW) appears to be an economical and environmental-friendly practice to treat nitrogen-enriched (waste) water. Nevertheless, the removal mechanisms in SCW are complicated and rather time-consuming to conduct and assessment the efficiency of new experiments. This work mined data from literature and developed the machine learning models to elucidate the effect of influent inputs and predict ammonium removal rate (ARR) in SCW treatment. 755 sets and 11 attributes were applied in four modeled algorithms, including Random forest, Cubist, Support vector machines, and K-nearest neighbors. Six out of ten input features including ammonium (NH 4 ), total nitrogen (TN), hydraulic loading rate (HLR), the filter height (i.e., Height), aeration mode (i.e., Aeration), and types of inlet feeding (i.e., Feeding) have posed pronounced influences on the ARR. The Cubist algorithm appears the most optimal model showing the lowest RMSE i.e., 0.974 and the highest R 2 i.e., 0.957. The contribution of variables followed the order of NH 4 , HLR, TN, Aeration, Height and Feeding corresponding to 97, 93, 71, 49, 34, and 34%, respectively. The generalization ability to forecast ARR using testing data achieved the R 2 of 0.970 and the RMSE of 1.140 g/m 2 d, indicating that Cubist is a reliable tool for ARR prediction. User interface and web tool of final predictive model are provided to facilitate the application for designing and developing SCW system in real practice. • Four supervised ML algorithms were utilized to predict ARR by SCW. • Six over ten input variables were found to be the most pronounced parameters. • NH 4 -N (97%), HLR (93%), and TN (71%) contributed the most to model. • Cubist achieved the highest generalization with R 2 of 0.957 and RMSE of 0.974 g/m 2 d. • Strong correlation between the predicted and actual results by Cubist.