Numerical simulation of phosphorus removal in horizontal subsurface flow constructed wetlands