Hydrus (software)

Last updated
Hydrus
Developer(s) PC Progress
Stable release
3.02
Operating system Windows 10, Windows 8, Windows 7, Windows Vista, Windows XP
Type Hydrological modelling
License Public domain software (Hydrus-1D)
Proprietary (HYDRUS 2D/3D)
Website http://www.pc-progress.com/en/default.aspx?hydrus-3d

Hydrus is a suite of Windows-based modeling software that can be used for analysis of water flow, heat and solute transport in variably saturated porous media (e.g., soils). HYDRUS suite of software is supported by an interactive graphics-based interface for data-preprocessing, discretization of the soil profile, and graphic presentation of the results. While HYDRUS-1D simulates water flow, solute and heat transport in one-dimension, and is a public domain software, HYDRUS 2D/3D extends the simulation capabilities to the second and third dimensions, and is distributed commercially.

Contents

History

HYDRUS 1D

MainPg 04.gif
History of Hydrus-1D Development Hydrus1d history.tiff
History of Hydrus-1D Development

HYDRUS-1D traces its roots to the early work of van Genuchten [1] [2] and his SUMATRA and WORM models, as well as later work by Vogel (1987) [3] and Kool and van Genuchten (1989) [4] and their SWMI and HYDRUS models, respectively. While Hermitian cubic finite element numerical schemes were used in SUMATRA and linear finite elements in WORM and the older HYDRUS code for solution of both the water flow and solute transport equations, SWMI used finite differences to solve the flow equation.

Various features of these four early models were combined first in the DOS-based SWMI_ST model (Šimůnek et al., 1993), and later in the Windows-based HYDRUS-1D simulator (Šimůnek et al., 1998). [5] After releasing versions 1 (for 16-bit Windows 3.1) and 2 (for 32-bit Windows 95), the next two major updates (versions 3 and 4) were released in 2005 and 2008. These last two versions included additional modules applicable to more complex biogeochemical reactions than the standard HYDRUS modules.

While the standard modules of HYDRUS-1D can simulate the transport of solutes that are either fully independent or involved in the sequential first-order degradation chains, the two new modules can consider mutual interactions between multiple solutes, such as cation exchange and precipitation/dissolution.

Version 3 included the UNSATCHEM module (Suarez and Šimůnek, 1997 [6] ) for simulating carbon dioxide transport as well as the multi-component transport of major ions. The UNSATCHEM major ion module was recently included also in version 2 of HYDRUS (2D/3D) (Šimůnek et al., 2011 [7] ). Version 4 of HYDRUS-1D includes now not only the UNSATCHEM module, but also the HP1 program (Jacques and Šimůnek, 2005 [8] ), which resulted from coupling HYDRUS-1D with the biogeochemical program PHREEQC. [9]

HYDRUS 2D/3D

The current HYDRUS (2D/3D) suite of software and their predecessors have a long history. The origin of these models can be traced back to the early work of Dr. Shlomo Neuman and collaborators (e.g., Neuman, 1972 [10] ) who developed their UNSAT model at the Hydraulic Engineering Laboratory of Technion – Israel Institute of Technology, in Haifa, Israel, long before the introduction of personal computers. UNSAT was a finite element model simulating water flow in two-dimensional variably-saturated domains as described with the Richards equation. The model additionally considered root water uptake as well as a range of pertinent boundary conditions required to ensure wide applicability of the model. UNSAT was later modified by Davis and Neuman (1983) [11] at the University of Arizona, Tucson, such that the model could be run on personal computers.

This last version of UNSAT formed the basis of the SWMII model developed by Vogel (1987) during his stay at Wageningen University, the Netherlands. SWMII significantly extended the capabilities and ease of use of UNSAT. The code simulated variably-saturated water flow in two-dimensional transport domains, implemented the van Genuchten soil hydraulic functions (van Genuchten, 1980 [12] ) and modifications thereof, considered root water uptake by taking advantage of some of the features of the SWATRE model (Feddes et al., 1978 [13] ), and included scaling factors to enable simulations of flow in heterogeneous soils. The code also allowed the flow region to be composed of nonuniform soils having an arbitrary degree of local anisotropy. SWMII was a direct predecessor of the SWMS_2D model (Šimůnek et al., 1992 [14] ) developed later at US Salinity Laboratory.

The SWMS_2D model (Šimůnek et al., 1992 [15] ) considerably extended the capabilities of SWMII by including provisions for solute transport. Solute transport was described using the standard advection-dispersion equation that included linear sorption, first-order degradation in both the liquid and solid phases, and zero-order production in both phases. Several other numerical improvements were at the time also implemented in SWMS_2D. These included solution of the mixed form of the Richards equation as suggested by Celia et al. (1990), [16] thus providing excellent mass balances in the water flow calculations. While SWMII could simulate water flow in either two-dimensional vertical or horizontal planes, SWMS_2D extended the range of applications also to three-dimensional axisymmetrical flow domains around a vertical axis of symmetry. Examples are flow to a well, infiltration from a surface ring or tension disc infiltrometer, and infiltration from a surface or subsurface dripper.

History of HYDRUS 2D/3D Development Hydrus2d3dhistory.tiff
History of HYDRUS 2D/3D Development

The first major upgrade of SWMS_2D was released under the name CHAIN_2D (Šimůnek et al., 1994b [17] ). This model greatly expanded upon the capabilities of SWMS_2D by including, among other things, sequential first-order solute decay chains and heat transport. The temperature dependence of the soil hydraulic properties was included by considering the effects of temperature on surface tension, dynamic viscosity and the density of water. The heat transport equation in CHAIN_2D considered transport due to conduction and advection with flowing water. The solute transport equations considered advective-dispersive transport in the liquid phase, as well as diffusion in the gaseous phase. The transport equations also included provisions for nonlinear nonequilibrium reactions between the solid and liquid phases, linear equilibrium reactions between the liquid and gaseous phase, zero-order production and two first-order degradation reactions: one which was independent of other solutes, and one which provided the coupling between solutes involved in the sequential first-order decay reactions.

The SWMS_2D and CHAIN_2D models formed the bases of versions 1.0 (for 16-bit Windows 3.1) and 2.0 (for 32-bit Windows 95) of HYDRUS-2D (Šimůnek et al., 1999 [18] ). A unique feature of HYDRUS-2D was that it used a Microsoft Windows-based Graphics User Interface (GUI) to manage the input data required to run the program, as well as for nodal discretization and editing, parameter allocation, problem execution, and visualization of results. It could handle flow regions delineated by irregular boundaries, as well as three-dimensional regions exhibiting radial symmetry about the vertical axis. The code includes the MeshGen2D mesh generator, which was specifically designed for variably-saturated subsurface flow and transport problems. The mesh generator may be used for defining very general domain geometries, and for discretizing the transport domain into an unstructured finite element mesh. HYDRUS-2D has been recently fully replaced with HYDRUS (2D/3D) as described below.

The HYDRUS (2D/3D) (version 1) software package (Šimůnek et al., 2006; [19] Šejna and Šimůnek, 2007 [20] ) is an extension and replacement of HYDRUS-2D (version 2.0) and SWMS_3D (Šimůnek et al., 1995 [21] ). This software package is a complete rewrite of HYDRUS-2D and its extensions for two- and three-dimensional geometries. In addition to features and processes available in HYDRUS-2D and SWMS_3D, the new computational modules of HYDRUS (2D/3D) consider (a) water flow and solute transport in a dual-porosity system, thus allowing for preferential flow in fractures or macropores while storing water in the matrix, (b) root water uptake with compensation, (c) the spatial root distribution functions, (d) the soil hydraulic property models of Kosugi and Durner, (e) the transport of viruses, colloids, and/or bacteria using an attachment/detachment model, filtration theory, and blocking functions, (f) a constructed wetland module (only in 2D), (g) the new hysteresis model to eliminate pumping by keeping track of historical reversal points, and many other options.

Simulated processes

Both HYDRUS models may be used to simulate movement of water, heat, and multiple solutes in variably saturated media. Both programs use linear finite elements to numerically solve the Richards equation for saturated-unsaturated water flow and Fickian-based advection dispersion equations for both heat and solute transport. The flow equation also includes a sink term to account for water uptake by plant roots as a function of both water and salinity stress. The unsaturated soil hydraulic properties can be described using van Genuchten, Brooks and Corey, modified van Genuchten, Kosugi, and Durner type analytical functions. The heat transport equation considers conduction as well as advection with flowing water. The solute transport equations assume advective-dispersive transport in the liquid phase, and diffusion in the gaseous phase. The transport equations further include provisions for nonlinear and/or non-equilibrium reactions between the solid and liquid phases, linear equilibrium reactions between the liquid and gaseous phases, zero-order production, and two first-order degradation reactions: one which is independent of other solutes, and one which provides the coupling between solutes involved in sequential first order decay reactions. In addition, physical non-equilibrium solute transport can be accounted for by assuming a two-region, dual-porosity type formulation which partitions the liquid phase into mobile and immobile regions.

HYDRUS models may be used to analyze water and solute movement in unsaturated, partially saturated, or fully saturated homogeneous of layered media. The codes incorporates hysteresis by assuming that drying scanning curves are scaled from the main drying curve, and wetting scanning curves from the main wetting curve. Root water uptake can be simulated as a function of both water and salinity stress, and can be either compensated or uncompensated. The HYDRUS software packages additionally implement a Marquardt–Levenberg type parameter estimation technique for inverse estimation of soil hydraulic and/or solute transport and reaction parameters from measured transient or steady-state flow and/or transport data. The programs are for this purpose written in such a way that almost any application that can be run in a direct mode can equally well be run in an inverse mode, and thus for model calibration and parameter estimation.

The HYDRUS packages use a Microsoft Windows-based graphical user interface (GUI) to manage the input data required to run the program, as well as for nodal discretization and editing, parameter allocation, problem execution, and visualization of results. All spatially distributed parameters, such as those for various soil horizons, the root water uptake distribution, and the initial conditions for water, heat and solute movement, are specified in a graphical environment. The program offers graphs of the distributions of the pressure head, water content, water and solute fluxes, root water uptake, temperature and solute concentrations in the subsurface at pre-selected times. Also included is a small catalog of unsaturated soil hydraulic properties, as well as pedotransfer functions based on neural networks.

Both HYDRUS models also consider various provisions for simulating non-equilibrium flow and transport. The flow equation for the latter purpose can consider dual-porosity-type flow with a fraction of the water content being mobile, and a fraction immobile. The transport equations additionally were modified to allow consideration of kinetic attachment/detachment processes of solutes to the solid phase, and hence of solutes having a finite size. This attachment/detachment feature has been used by many recently to simulate the transport of viruses, colloids, and bacteria.

HYDRUS model further include modules for simulating carbon dioxide transport (only HYDRUS-1D) and major ion chemistry modules, adopted from the UNSATCHEM program. HYDRUS-1D can thus be used in applications evaluating overall salinity, the concentration of individual soluble cations, as well as of the Sodium Adsorption Ratio and the Exchangeable Sodium Percentage.

Applications

Both HYDRUS-1D and HYDRUS (2D/3D) has been used in hundreds, if not thousands of applications referenced in peer-reviewed journal articles and many technical reports. Both software packages are also used in classrooms of many universities in courses covering Soil Physics, Processes in the Vadose Zone, or Vadose Zone Hydrology. A selected list of hundreds of applications of both HYDRUS software packages are given at:

http://www.pc-progress.com/en/Default.aspx?h3d-references

http://www.pc-progress.com/en/Default.aspx?h1d-references

The website also provides many specific applications in the libraries of HYDRUS projects at:

http://www.pc-progress.com/en/Default.aspx?h1d-library

http://www.pc-progress.com/en/Default.aspx?h3d-applications

HYDRUS software also provides capabilities for simulating water flow and solute transport for specialized domains.

Constructed Wetland Module

Constructed wetlands (CWs) are engineered water treatment systems that optimize the treatment processes found in natural environments. CWs are popular systems which efficiently treat various types of polluted water and are therefore sustainable, environmentally friendly solutions. A large number of physical, chemical and biological processes are simultaneously active and mutually influence each other. HYDRUS offers two biokinetic model formulations: (a) the CW2D module (Langergraber and Šimůnek, 2005 [22] ), and/or the CW M1 (Constructed Wetland Model #1) biokinetic model (Langergraber et al., 2009b [23] ).

Related Research Articles

<span class="mw-page-title-main">Hydrogeology</span> Study of the distribution and movement of groundwater

Hydrogeology is the area of geology that deals with the distribution and movement of groundwater in the soil and rocks of the Earth's crust. The terms groundwater hydrology, geohydrology, and hydrogeology are often used interchangeably.

Soil moisture is the water content of the soil. It can be expressed in terms of volume or weight. Soil moisture measurement can be based on in situ probes or remote sensing methods.

Used in hydrogeology, the groundwater flow equation is the mathematical relationship which is used to describe the flow of groundwater through an aquifer. The transient flow of groundwater is described by a form of the diffusion equation, similar to that used in heat transfer to describe the flow of heat in a solid. The steady-state flow of groundwater is described by a form of the Laplace equation, which is a form of potential flow and has analogs in numerous fields.

<span class="mw-page-title-main">Water retention curve</span>

Water retention curve is the relationship between the water content, θ, and the soil water potential, ψ. This curve is characteristic for different types of soil, and is also called the soil moisture characteristic.

The United States Environmental Protection Agency (EPA) Storm Water Management Model (SWMM) is a dynamic rainfall–runoff–subsurface runoff simulation model used for single-event to long-term (continuous) simulation of the surface/subsurface hydrology quantity and quality from primarily urban/suburban areas.

<span class="mw-page-title-main">Hydrological transport model</span>

An hydrological transport model is a mathematical model used to simulate the flow of rivers, streams, groundwater movement or drainage front displacement, and calculate water quality parameters. These models generally came into use in the 1960s and 1970s when demand for numerical forecasting of water quality and drainage was driven by environmental legislation, and at a similar time widespread access to significant computer power became available. Much of the original model development took place in the United States and United Kingdom, but today these models are refined and used worldwide.

The Richards equation represents the movement of water in unsaturated soils, and is attributed to Lorenzo A. Richards who published the equation in 1931. It is a quasilinear partial differential equation; its analytical solution is often limited to specific initial and boundary conditions. Proof of the existence and uniqueness of solution was given only in 1983 by Alt and Luckhaus. The equation is based on Darcy-Buckingham law representing flow in porous media under variably saturated conditions, which is stated as

Groundwater models are computer models of groundwater flow systems, and are used by hydrologists and hydrogeologists. Groundwater models are used to simulate and predict aquifer conditions.

<span class="mw-page-title-main">Soil salinity control</span> Controlling the problem of soil salinity

Soil salinity control refers to controlling the process and progress of soil salinity to prevent soil degradation by salination and reclamation of already salty (saline) soils. Soil reclamation is also known as soil improvement, rehabilitation, remediation, recuperation, or amelioration.

FEFLOW is a computer program for simulating groundwater flow, mass transfer and heat transfer in porous media and fractured media. The program uses finite element analysis to solve the groundwater flow equation of both saturated and unsaturated conditions as well as mass and heat transport, including fluid density effects and chemical kinetics for multi-component reaction systems.

<span class="mw-page-title-main">HydroGeoSphere</span>

HydroGeoSphere (HGS) is a 3D control-volume finite element groundwater model, and is based on a rigorous conceptualization of the hydrologic system consisting of surface and subsurface flow regimes. The model is designed to take into account all key components of the hydrologic cycle. For each time step, the model solves surface and subsurface flow, solute and energy transport equations simultaneously, and provides a complete water and solute balance.

MIKE 11 is a computer program that simulates flow and water level, water quality and sediment transport in rivers, flood plains, irrigation canals, reservoirs and other inland water bodies. MIKE 11 is a 1-dimensional river model. It was developed by DHI.

SVFLUX is a finite element seepage analysis program developed by SoilVision Systems Ltd.. The software is designed to analyze both saturated and unsaturated flow through the ground through the solving of Richard's equation. The program is used in the fields of civil engineering and hydrology in order to analyze seepage and groundwater regional flow. The software is used for the calculation of flow rates, pore-water pressures, and pumping rates associated with regional groundwater flow. The software can be coupled with CHEMFLUX in order to calculate diffusion, advection, and decay rates or with SVHEAT in order to calculate thermal gradients and freeze/thaw fronts.

<span class="mw-page-title-main">SHETRAN</span>

SHETRAN is a hydrological modelling system for water flow, solute and sediment transport in river catchments. SHETRAN is a physically based, distributed model (PBDM) that can simulate the entire land phase of the hydrologic cycle including surface water flow and groundwater flow. The plan area of the catchment in SHETRAN is usually in the range of one to a few thousand square kilometres and the horizontal depth of the subsurface is usually less than 100m.

DPHM-RS is a semi-distributed hydrologic model developed at University of Alberta, Canada.

<span class="mw-page-title-main">TELEMAC</span> Open source software for computational fluid dynamics

In computational fluid dynamics, TELEMAC is short for the open TELEMAC-MASCARET system, or a suite of finite element computer program owned by the Laboratoire National d'Hydraulique et Environnement (LNHE), part of the R&D group of Électricité de France. After many years of commercial distribution, a Consortium was officially created in January 2010 to organize the open source distribution of the open TELEMAC-MASCARET system now available under GPLv3.

Geochemical modeling or theoretical geochemistry is the practice of using chemical thermodynamics, chemical kinetics, or both, to analyze the chemical reactions that affect geologic systems, commonly with the aid of a computer. It is used in high-temperature geochemistry to simulate reactions occurring deep in the Earth's interior, in magma, for instance, or to model low-temperature reactions in aqueous solutions near the Earth's surface, the subject of this article.

Reactive transport modeling in porous media refers to the creation of computer models integrating chemical reaction with transport of fluids through the Earth's crust. Such models predict the distribution in space and time of the chemical reactions that occur along a flowpath. Reactive transport modeling in general can refer to many other processes, including reactive flow of chemicals through tanks, reactors, or membranes; particles and species in the atmosphere; gases exiting a smokestack; and migrating magma.

<span class="mw-page-title-main">Finite water-content vadose zone flow method</span>

The finite water-content vadose zone flux method represents a one-dimensional alternative to the numerical solution of Richards' equation for simulating the movement of water in unsaturated soils. The finite water-content method solves the advection-like term of the Soil Moisture Velocity Equation, which is an ordinary differential equation alternative to the Richards partial differential equation. The Richards equation is difficult to approximate in general because it does not have a closed-form analytical solution except in a few cases. The finite water-content method, is perhaps the first generic replacement for the numerical solution of the Richards' equation. The finite water-content solution has several advantages over the Richards equation solution. First, as an ordinary differential equation it is explicit, guaranteed to converge and computationally inexpensive to solve. Second, using a finite volume solution methodology it is guaranteed to conserve mass. The finite water content method readily simulates sharp wetting fronts, something that the Richards solution struggles with. The main limiting assumption required to use the finite water-content method is that the soil be homogeneous in layers.

The Marine Unsaturated Model is a two-dimensional finite element model capable of simulating the migration of water and solutes in saturated-unsaturated porous media while accounting for the impact of solute concentration on water density and viscosity, as saltwater is heaving and more viscous than freshwater. The detailed formulation of the MARUN model is found in and. The model was used to investigate seepage flow in trenches and dams, the migration of brine following evaporation and, submarine groundwater discharge, and beach hydrodynamics to explain the persistence of some of the Exxon Valdez oil in Alaska beaches.

References

  1. van Genuchten, M. Th. 1978. Mass transport in saturated-unsaturated media: one-dimensional solutions. Research Rep. No. 78-WR-11. Water Resources Program. Princeton Univ. Princeton, NJ.
  2. van Genuchten, M. Th. 1987. A numerical model for water and solute movement in and below the root zone. Research Report No 121. U.S. Salinity laboratory, USDA, ARS, Riverside, California.
  3. Vogel, T. 1987. SWMII – Numerical model of two-dimensional flow in a variably saturated porous medium. Research Rep. No. 87. Dept. of Hydraulics and Catchment Hydrology. Agricultural Univ., Wageningen, The Netherlands.
  4. Kool, J. B., and M. Th. van Genuchten. 1989. HYDRUS, One-dimensional variably saturated flow and transport model, including hysteresis and root water uptake. Version 3.2. User Manual. Hydrogeologic Inc., Herndon, VA, 116 p.
  5. Šimůnek, J., M. Šejna, and M. Th. van Genuchten. 1998. The HYDRUS-1D software package for simulating the one-dimensional movement of water, heat, and multiple solutes in variably-saturated media. Version 1.0. IGWMC – TPS – 70, International Ground Water Modeling Center, Colorado School of Mines, Golden, Colorado, 186pp.
  6. Suarez, D. L., and J. Šimůnek. 1997. UNSATCHEM: Unsaturated water and solute transport model with equilibrium and kinetic chemistry. Soil Sci. Soc. Am. J., 61:1633–1646.
  7. Šimůnek, J., M. Th. van Genuchten, and M. Šejna. 2011. The HYDRUS Software Package for Simulating Two- and Three-Dimensional Movement of Water, Heat, and Multiple Solutes in Variably-Saturated Media. Technical Manual, Version 2.0, PC Progress, Prague, Czech Republic, pp. 258.
  8. Jacques, D., and J. Šimůnek. 2005. User Manual of the Multicomponent Variably-Saturated Flow and Transport Model HP1, Description, Verification and Examples. Version 1.0. SCK•CEN-BLG-998, Waste and Disposal, SCK•CEN, Mol, Belgium, 79 pp.
  9. Parkhurst, D. L., and C. A. J. Appelo. 1999. User’s guide to PHREEQC (Version 2): A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations. Water-Resour. Invest. Rep. 99-4259. USGS, Denver, CO.
  10. Neuman, S. P. 1972. Finite element computer programs for flow in saturated-unsaturated porous media. Second Annual Report, Project No. A10-SWC-77, Hydraulic Engineering Lab., Technion, Haifa, Israel.
  11. Davis, L. A., and S. P. Neuman. 1983. Documentation and user's guide: UNSAT2 – Variably saturated flow model. Final Rep., WWL/TM-1791-1, Water, Waste & Land, Inc., Ft. Collins, CO.
  12. van Genuchten, M. Th. 1980. A closedform equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J., 44:892898.
  13. Feddes, R. A., P. J. Kowalik, and H. Zaradny. 1978. Simulation of Field Water Use and Crop Yield. John Wiley & Sons, New York, NY.
  14. Šimůnek, J., T. Vogel and M. Th. van Genuchten. 1992. The SWMS_2D code for simulating water flow and solute transport in two-dimensional variably saturated media, Version 1.1, Research Report No. 126, U. S. Salinity Laboratory, USDA, ARS, Riverside, CA.
  15. Šimůnek, J., T. Vogel and M. Th. van Genuchten. 1992. The SWMS_2D code for simulating water flow and solute transport in two-dimensional variably saturated media, Version 1.1, Research Report No. 126, U. S. Salinity Laboratory, USDA, ARS, Riverside, CA.
  16. Celia, M. A., and E. T. Bououtas, R. L. Zarba. 1990. A general mass-conservative numerical solution for the unsaturated flow equation. Water Resour. Res., 26:1483-1496.
  17. Šimůnek, J., and M. Th. van Genuchten. 1994. The CHAIN_2D code for simulating two-dimensional movement of water flow, heat, and multiple solutes in variably-saturated porous media. Version 1.1, Research Report No 136, U.S. Salinity laboratory, USDA, ARS, Riverside, California.
  18. Šimůnek, J., M. Šejna, and M. Th. van Genuchten. 1999. The HYDRUS-2D software package for simulating two-dimensional movement of water, heat, and multiple solutes in variably saturated media. Version 2.0, IGWMC – TPS – 53, International Ground Water Modeling Center, Colorado School of Mines, Golden, Colorado, 251pp.
  19. Šimůnek, J., M. Th. van Genuchten, and M. Šejna. 2006. The HYDRUS Software Package for Simulating Two- and Three-Dimensional Movement of Water, Heat, and Multiple Solutes in Variably-Saturated Media, Technical Manual, Version 1.0, PC Progress, Prague, Czech Republic, pp. 241.
  20. Šejna, M., and J. Šimůnek. 2007. HYDRUS (2D/3D): Graphical User Interface for the HYDRUS Software Package Simulating Two- and Three-Dimensional Movement of Water, Heat, and Multiple Solutes in Variably-Saturated Media, published online at www.-pc-progress.cz, PC-Progress, Prague, Czech Republic.
  21. Šimůnek, J., K. Huang, and M. Th. van Genuchten. 1995. The SWMS_3D code for simulating water flow and solute transport in three-dimensional variably saturated media. Version 1.0, Research Report No. 139, U.S. Salinity Laboratory, USDA, ARS, Riverside, California, 155 pp.
  22. Langergraber, G., and J. Šimůnek, Modeling variably-saturated water flow and multi-component reactive transport in constructed wetlands, Vadose Zone J., 4(4), 924–938, 2005.
  23. Langergraber, G., D. Rousseau, J. García, and J. Mena, CWM1 – A general model to describe biokinetic processes in subsurface flow constructed wetlands, Water Sci. Technol., 59(9), 1687-1697, 2009.