Modeling of flow and BOD fate in horizontal subsurface flow constructed wetlands

Abstract A numerical modeling effort of flow and Biochemical Oxygen Demand (BOD) fate in horizontal subsurface flow (HSF) constructed wetlands (CWs), under Mediterranean conditions, is presented. The Visual MODFLOW family computer code, based on the finite difference method, was used for the numerical simulation of flow in five pilot-scale HSF CW units containing various vegetation and porous material types, and operating under various temperature and hydraulic residence time (HRT) conditions. BOD fate and transport was simulated using the MT3DMS computer code. Experimental data from these CW units were used in estimating the BOD removal coefficient, required in the above modeling effort. Values of this removal coefficient were found to depend on temperature and HRT. Relations are presented for the prediction of this coefficient. The model was calibrated and verified using the measured BOD concentrations, and was used in test runs comparing performance under various vegetation, porous media size, temperature and HRT conditions. Model performance was found acceptable, indicating its usefulness in the simulation of BOD fate in CWs and in the design of these facilities.

[1]  Jan Vymazal,et al.  The use of sub-surface constructed wetlands for wastewater treatment in the Czech Republic: 10 years experience , 2002 .

[2]  C. Zheng,et al.  Applied contaminant transport modeling , 2002 .

[3]  Vassilios A. Tsihrintzis,et al.  Effect of temperature, HRT, vegetation and porous media on removal efficiency of pilot-scale horizontal subsurface flow constructed wetlands , 2007 .

[4]  Vassilios A Tsihrintzis,et al.  Performance and modeling of a vertical flow constructed wetland–maturation pond system , 2011, Journal of environmental science and health. Part A, Toxic/hazardous substances & environmental engineering.

[5]  Günter Langergraber,et al.  Recent developments in numerical modelling of subsurface flow constructed wetlands. , 2009, The Science of the total environment.

[6]  G Tchobanoglous,et al.  Time‐Dependent Retardation Model for Chemical Oxygen Demand Removal in a Subsurface‐Flow Constructed Wetland for Winery Wastewater Treatment , 2001, Water environment research : a research publication of the Water Environment Federation.

[7]  J. N. Petersen,et al.  Modeling Microbial Transport and Biodegradation in a Dual-porosity System , 1999 .

[8]  J L G Kumar,et al.  A review on numerous modeling approaches for effective, economical and ecological treatment wetlands. , 2011, Journal of environmental management.

[9]  Robert H. Kadlec,et al.  The inadequacy of first-order treatment wetland models , 2000 .

[10]  Günter Langergraber,et al.  Modeling of Processes in Subsurface Flow Constructed Wetlands: A Review , 2008 .

[11]  P. Domenico,et al.  Physical and chemical hydrogeology , 1990 .

[12]  David Steer,et al.  Efficiency of small constructed wetlands for subsurface treatment of single-family domestic effluent , 2002 .

[13]  C. Zheng.  A Modular Three-Dimensional Transport Model for Simulation of Advection, Dispersion and Chemical Reaction of Contaminants in Groundwater Systems , 1990 .

[14]  Christos S. Akratos,et al.  An artificial neural network model and design equations for BOD and COD removal prediction in horizontal subsurface flow constructed wetlands , 2008 .

[15]  John F. Brady,et al.  Dispersion in fixed beds , 1985, Journal of Fluid Mechanics.

[16]  Vassilios A. Tsihrintzis,et al.  Simulation of hydrodynamics and nitrogen transformation processes in HSF constructed wetlands and porous media using the advection–dispersion-reaction equation with linear sink-source terms , 2011 .

[17]  Hydrodynamic and boundary-layer dispersion in bidisperse porous media , 1999 .

[18]  Sherwood C. Reed,et al.  Subsurface flow wetlands : a performance evaluation , 1995 .

[19]  V. Tsihrintzis,et al.  Determination of Forchheimer equation coefficients a and b , 2007 .

[20]  Alexandros Stefanakis,et al.  Effect of wastewater step-feeding on removal efficiency of pilot-scale horizontal subsurface flow constructed wetlands , 2011 .

[21]  Hanadi S. Rifai,et al.  Comparison of biodegradation kinetics with an instantaneous reaction model for groundwater , 1990 .

[22]  J. Bear.  Hydraulics of Groundwater , 1979 .

[23]  Chris C. Tanner,et al.  Effect of loading rate and planting on treatment of dairy farm wastewaters in constructed wetlands—I. Removal of oxygen demand, suspended solids and faecal coliforms , 1995 .

[24]  Guenter Langergraber.  Numerical modelling: a tool for better constructed wetland design? , 2011, Water science and technology : a journal of the International Association on Water Pollution Research.

[25]  Stefano Marsili-Libelli,et al.  Identification of dynamic models for horizontal subsurface constructed wetlands , 2005 .

[26]  Margaret Greenway,et al.  Constructed wetlands in Queensland: Performance efficiency and nutrient bioaccumulation , 1999 .

[27]  Robert H. Kadlec,et al.  Deterministic and stochastic aspects of constructed wetland performance and design , 1997 .

[28]  E. A. Sudicky,et al.  Contaminant transport in fractured porous media: Analytical solution for a single fracture , 1981 .

[29]  D McNevin,et al.  Alternative analysis of BOD removal in subsurface flow constructed wetlands employing Monod kinetics. , 2001, Water research.

[30]  D. M. Griffin,et al.  The effect of temperature on biochemical oxygen demand removal in a subsurface flow wetland , 1999 .

[31]  Peter A Vanrolleghem,et al.  Model-based design of horizontal subsurface flow constructed treatment wetlands: a review. , 2004, Water research.

[32]  Karlheinz Spitz,et al.  A Practical Guide to Groundwater and Solute Transport Modeling , 1996 .

[33]  Arlen W. Harbaugh,et al.  User's documentation for MODFLOW-96, an update to the U.S. Geological Survey modular finite-difference ground-water flow model , 1996 .

[34]  Jian Huang,et al.  Nitrogen removal in constructed wetlands employed to treat domestic wastewater. , 2000 .

[35]  D Uhlmann,et al.  Bod removal rates of waste stabilization ponds as a function of loading, retention time, temperature and hydraulic flow pattern , 1979 .

[36]  V. Tsihrintzis,et al.  Constructed wetlands for wastewater and activated sludge treatment in north Greece: a review. , 2010, Water science and technology : a journal of the International Association on Water Pollution Research.