Background: This study assessed the removal of organic material and nutrients from full-scale subsurface flow (SSF) constructed wetlands (CWs) followed by anaerobic stabilization ponds under environmental conditions. Methods: The effluents were distributed evenly in 12 reed beds. Samples were taken twice monthly for a total of 6 months from several points in the wetland. Biochemical oxygen demand (BOD), chemical oxygen demand (COD), total suspended solids (TSS), and nutrient removal from the system and the longitudinal effect of the reed beds for removal of pollutions were determined. A full-scale model of flow, BOD, and nutrients in SSF in the CWs is presented. Results: The flow rate and concentrations of parameters indicated that removal of organic matter and nutrients in the cold months decreased rather than in the hot months, as expected. The removal efficiency for BOD, COD, and TSS and the strongest biological interactions showed no uniform trends. The beds showed the highest removal rates in the first few meters of bed. The hybrid Monod-Plug flow regime and the Stover-Kincannon models showed the best fit for the kinetics of the processes. Umax in the Stover-Kincannon model was 3.64 mg/l.d for nitrogen and 0.24 mg/l.d for phosphorus. These values are very low, which indicates lower consumption and inefficiency of the system for removing nitrogen and phosphorus. Conclusion: It can be concluded that the SSF in CWs are able to treat average wastewater as effectively as common mechanical systems at lower cost.