Subscribe Now Subscribe Today
Research Article

Phosphorus Transport Through a Saturated Soil Column: Comparison Between Physical Modeling and HYDRUS-3D Outputs

A.A. Naseri, Y. Hoseini, H. Moazed, F. Abbasi, H.M.V. Samani and S.A. Sakebi
Facebook Twitter Digg Reddit Linkedin StumbleUpon E-mail

The objective of the present study was to understand Phosphorus transport from soil columns at different water flux densities and to compare computer simulation results using HYDRUS-3D and the physical model against column-experiment data. Towards that end, two water flux densities (0.0075, 0.0065 m sec-1) and P concentration rates (33 mol cm-3) were used. The soil used has a silty clay- loam texture. The numerical model (HYDRUS-3D) successfully predicted P transport in the present experiment. Overall, the HYDRUS-3D model successfully simulated the water flow in the columns; however, it overestimated the final adsorbed PO4 concentrations in the soil. The present column experiment is useful for assessing relative behavior of P in increasing the movement of phosphorus ions out of the soil profile and into ground waters.

Related Articles in ASCI
Similar Articles in this Journal
Search in Google Scholar
View Citation
Report Citation

  How to cite this article:

A.A. Naseri, Y. Hoseini, H. Moazed, F. Abbasi, H.M.V. Samani and S.A. Sakebi, 2011. Phosphorus Transport Through a Saturated Soil Column: Comparison Between Physical Modeling and HYDRUS-3D Outputs. Journal of Applied Sciences, 11: 815-823.

DOI: 10.3923/jas.2011.815.823

Received: September 28, 2010; Accepted: January 18, 2011; Published: February 24, 2011


Phosphoros is a critical element in both natural and agricultural ecosystems throughout the world (Onweremadu, 2007) as its limited availability is often the main constraint for plant growth in highly weathered soils of the tropics (Bunemann et al., 2004). Phosphorus deficiency problems are common in well-weathered Oxisols and Ultisols because of strong acidic reactions and an abundance of Al and Fe ions (Saleque et al., 2004). A substantial problem to human health is contamination of ground waters through leaching of phosphorus. Usage of fertilizers and leakage of phosphorus through the soil has been a concern throughout the world. A considerable number of laboratory experiments have been conducted to enable modeling of contaminant transport to improve the simulation capability of models. HYDRUS-3D is a software package for simulating water, heat and solute transport in three-dimensional, variably saturated media (Simunek et al., 2005). Many researchers have used Hydrus to predict solute transport in the soil column (Pang et al., 2000; Ventrella et al., 2000). Phosphorus (P) in solutions is attached to positively charged sites on the surface of soil particles and is strongly adsorbed to the soil, forming stable bonds with the soil’s functional group. The soil adsorbs P molecules found in the soil solution; phosphorus then interacts with soil particles in its exchangeable form. Although phosphate is strongly adsorbed to soil particles, many recent studies have shown that phosphate can easily be found in subsurface fresh water (Enright and Madramootoo, 2003; Simard, 2005). The reason behind such results is the under-sorption and desorption process that is controlled by the equilibrium constant. As the concentration of phosphate increases in the soil solution, the phosphate is adsorbed to soil while, when the concentration of phosphate in the solution is depleted, the phosphate that has been adsorbed in the soil dissolves back into solution. Also, the macropore flow path is important to load phosphate through the soil (Stamm et al., 1998; Hooda et al., 1999). In order to improve predictions of phosphate transport in soil, computer modeling can be useful for simulating P distribution under laboratory and field conditions. Phillips (2006) investigated the accuracy of the HYDRUS-2D model in simulating the transport of reactive chemicals like potassium in undisturbed soil columns. The results showed that the HYDRUS model was successful in simulating water flow but it was not very conclusive with regard to reactive solute transport. Ben-Gal and Dudley (2003) simulated P movement and P distribution patterns under continuous application of P fertilizer from a subsurface, point-source emitter applied HYDRUS-2D. Both simulated and experimental results showed that an increase in P uptake by plants and greater P mobility, was observed under continuous applications of water and fertilizer as opposed to intermittent applications. However, Ben-Gal and Dudley (2003) concluded that the Langmuir isotherm, provided within the HYDRUS-2D model, is inadequate for describing P sorption. It assumes instantaneous equilibrium between the solid and solution interfaces. Rassam and Cook (2002) modeled water flow and solute transport in soil with a high concentration of acid sulfatteof tile drained fields using the HYDRUS-2D model. They concluded that numerical modeling as employed by the HYDRUS-2D model provided an effective tool in understanding the mechanisms involved in such a complex system, hence is a valuable management tool.

Many numerical models have been tested, such as NCSWAP, LEACHM (Hutson and Wagenet, 1992; RZWQM Team, 1995), CHAIN-2D (Simunek and Van Genuchten, 1994), CHAIN _IR (Zhang, 1997) and HYDRUS-2D (Simcnek et al., 1999b). As a result of increasing concerns over groundwater pollution with phosphorus, the objectives of the present study were defined as follows:

Study P transport from soil columns with water flux at different P concentrations
Compare computer simulation results using numerical and physical models with column-experiment data


The research was conducted during the cropping season between 2009 and 2010 at Omidie, southeastern Iran. The study site is located at latitude 31°-53'-14" N and 48°-40'-34" E, with an elevation of 55 meters (handheld Global Positioning System – GPS) (Garmin Ltd, Kansas, USA). The predominant parent material underlying the area from which most soil was formed is coastal plain sand (Benin formation) of the Miocene-Oligocene geological era. These soils are referred to as silty clay loam and are characteristically acidic, have low cation exchange capacity, low base saturation and low fertility and are challenged by multiple nutrient deficiencies (Oti, 2002). The area has a semi-arid tropical climate (Igwe and Stahr, 2004) with two distinct seasons, namely the wet and dry seasons. Rainfall distribution is bimodal with peaks during the months of July and September. Temperatures are high and change only slightly during the year.

Some physical and chemical properties are illustrated in Table 1.

Soil columns: The columns had an ID of 10 and 100 cm long. They were made of transparent polyvinyl chloride (PVC). The physical model used in the experiments (designed and built in the Water Science Faculty of Shahid Chamran University, Ahvaz) consists of a PVC pipe 100 cm tall. The model had three parts - input, output and porous media- which were placed on a foot stool; water and contaminants were flowed from the bottom of the column. To prevent soil particles from falling off the top of the column, a Plexiglas plate was placed on the soil column. Water, along with contaminants, was flowed from the source through a pipe into the porous media. To prevent air from entering the system, a control valve was placed along this path. Sampling tubes were placed on the body model in flowing lengths 5, 12, 25, 50 and 75 cm from the bottom of the model. Drainage pipes were placed 32, 55 and 82 cm from the bottom of the model. The schema used in the physical model is shown in Fig. 1.

Water and phosphate application: The soil columns were saturated by adding water from the bottom of the column to reach saturated conditions and phosphate solution was applied at steady rate. The soil columns at different depths were monitored to collect the leachate. The water was applied at two different constant rates: 0.0075, 0.006 m sec-1 with a P concentration of 33 mol cm-3. The P concentration was calculated with a scanning spectrophotometer. BTC curves were extracted at different depths and are shown in Fig. 4 to 15. At the end of experiment at each depth, the soil was changed to execute experiments at another depth.

Soil hydraulic properties: The van Genuchten-Mualem model (Van Genuchten, 1980) was used to describe the water retention and conductivity (K(h)) functions for both matrix and macropore domains (Mualem, 1976):


Fig. 1: Schema of physical model

Table 1: The number of physical and chemical properties of soil used in the present study
SOC: Soil organic carbon, Alsat: Aluminium saturation, Ca: Calcium, Al: Aluminium CEC: Cation exchange capacity





Where :
(h) = The soil water content at the matrix head, (h)
h = The soil matrix head (L)
= The saturated water content (L3 L-3)
KS = The saturated hydraulic conductivity (LT-1)
he = The air-entry potential (L-1)
Kr = The relative hydraulic conductivity (-)
Se = The relative water saturation (-) m, n and α are the fitting parameters of the retention curve
= The residual water content (L3 L-3)
K(h) = Unsaturated hydraulic conductivity (LT-1)The values of α, m and n are obtained by fitting Eq. 2 to the soil water retention data using the RETC model (Van Genuchten et al., 1991)

Water flow equation:
The Richard equation describes one-dimensional transient water flow in soil (Richards, 1931):

Where :
(h) = The volumetric water content (L3 L-3)
K(h) = The unsaturated hydraulic conductivity (LT-1)
z = The vertical coordination (L) taken positively upward
h = The matrix head (L)
t = The time (T)

Solute transport equation: Solute transport in porous media can be described with the one-dimensional convective-dispersive equation (CDE) (Simunek et al., 1999a):

Fig. 2: Defined model in HYDROUS-3D and sampling depths

Where :
S = The sorbed solute concentration (MM-1)
ρb = The soil bulk density (ML-3)
qw = The volumetric water Flux (LT-1)
C = The solute concentration in solution (ML-3)
D = The effective dispersion coefficient (L2 T-1)

The effective dispersion coefficient (D) is given by Bear (1972):

λL = The longitudinal dispersive (L),
D0 = The aqueous ionic or molecular diffusion coefficient of nitrate in water (L2 T-1),
τ = The tortuosity factor given by Millington and Quirk (1961):


By manipulating Eq. 7, we obtained the modified Convective-dispersion equation:

Fig. 3: P transport and its concentration in soil column

Table 2: T-test analysis tools for equality of data



R = The retardation factor given by the following equation

Kd = The partition coefficient of solute (Cm3 g-1)
V = The pore-water velocity (LT-1)

Numerical simulation: The water flow and solute transport equations with initial and boundary conditions were solved numerically with the HYDRUS-3D code (Simunek et al., 1999b). The HYDRUS-3D code is based on the Galerkin finite elements method for a space weighting scheme, while the time derivatives for the solute transport equation were approximated by a Crank-Nicholson finite differences scheme. Figure 2 is shown in the definition model in HYDRUS-3d and sampling depths. Figure 3 shows P transport and its concentration in the soil column. The Student’s t-test was carried out on the data shown in Table 2.


Phosphate transport in soil columns: The overlay of the BTC curve from experiment and HYDRUS-3D are shown in Fig. 4 to 15. From information on Fig. 4 and 5 it can be seen that in the 5 cm soil column depth the hydrus model accurately predicted the phosphorus movement in the soil. The t-test shows no significant difference between the results of the model and the soil column experiment.

Fig. 4: BTC curve from experimental and HYDRUS-3D first test at 5 cm depth

Fig. 5: BTC curve from experimental and HYDRUS-3D second test at 5 cm depth

Fig. 6: BTC curve from experimental and HYDRUS-3D first test at 12 cm depth

Fig. 7: BTC curve from experimental and HYDRUS-3D second test at 12 cm depth

Fig. 8: BTC curve from experimental and HYDRUS-3D first test at 25 cm depth

Fig. 9: BTC curve from experimental and HYDRUS-3D second test at 25 cm depth

Fig. 10: BTC curve from experimental and HYDRUS-3D first test at 50 cm depth

Fig. 11: BTC curve from experimental and HYDRUS-3D second test at 50 cm depth

Fig. 12: BTC curve from experimental and HYDRUS-3D first test at 75 cm depth

Fig. 13: BTC curve from experimental and HYDRUS-3D second test at 75 cm depth

Fig. 14: BTC curve from experimental and HYDRUS-3D first test at 100 cm depth

Fig. 15: BTC curve from experimental and HYDRUS-3D second test at 100 cm depth

Fig. 6 and 7, comparing the BTC curve at 12 cm depth, show that the result of the two methods in the second experiment has a higher retardation factor, which indicates that the phosphorus adsorption isotherm is more effective in the hydrus model. This can be been seen in Fig. 8 and 9. The results show that the hydrus results at 25 cm depth most closely match column outputs. As illustrated in Fig. 10 and 11, results of the hydrus model at 50 cm depth show faster phosphorus movement than was actually observed in the first experiment. This phenomenon is demonstrated in Fig. 12 and 13, which show results at 75 cm depth. Figure 14 and 15 show the BTC curve at 100 cm depth. These figures illustrate that the result of the hydrus model at 100 cm depth in the soil column does not show good similarity in both experiments. The predicted solute movement by the model is slower than in reality.

A comparison between the HYDRUS-3D outputs and observed data from the soil column indicate that the model made sufficiently accurate predictions. Break-through curves are shown in Fig. 4 to 14. In all figures it is clear that the outputs of the hydrus model are close to the experiment results. The phosphorus leached through the soil profile and reached different depths in the soil column. Previous studies have indicated that although phosphate is strongly adsorbed to soil particles, its ions can move through soil columns. This phenomenon specifies the effect of macropores on the transport process. For this reason we must use light rates of P fertilizer and minimize irrigation water to reduce the loss of phosphorus. Additionally, it is obvious that the results of the experiment with field conditions may be different than with laboratory conditions and it is therefore necessary for this test to be done under field conditions. The difference between the model and experiments may be as a result of rapid movement of water through the soil without sufficient time for reaction between the soil and the solute (Hillel, 1998). As illustrated in Fig. 4 to 14, as the soil depth in the experiment increases, the difference between the model and soil column outputs also increases. This phenomenon may be related to the increase in error due to preferential flow paths and/or the choice of the incorrect adsorption isotherm.


The authors of this article would like to thank the Shahid Chamran University of Ahvaz, College of Graduate Studies and Research, who made this research possible through a doctoral thesis. Also, the authors would like to thank the Center of Excellence in Irrigation and Drainage Networks Management, Shahid Chamran University of Ahvaz, Iran for partial financial support of this research.

1:  Bear, J., 1972. Dynamics of Fluids in Porous Media. Dover Publicatons, Inc., New York, pp: 764.

2:  Bunemann, E.K., F. Steinebruner, P.C. Smithson, E. Frossard and A. Oberson, 2004. Phosphorus dynamics in a highly weathered soil as revealed by isotopic labelling techniques. Soil Sci. Soc. Am. J., 68: 1645-1655.
Direct Link  |  

3:  Enright, P. and C.A. Madramootoo, 2003. Phosphorus losses in surface runoff and subsurface drainage waters on two agricultural fields in Quebec. CSAE Annual Meeting Paper 03-111, McGill University.

4:  Hillel, D., 1988. Environmental Soil Physics. Academic Press, San Diego, CA, pp: 771.

5:  Hooda, P.S., M. Moynagh, I.F. Svoboda, A.C. Edwards, H.A. Anderson and G. Sym, 1999. Phosphorus loss in drainflow from intensively managed grassland soils. J. Environ. Qual., 28: 1235-1242.
Direct Link  |  

6:  Hutson, J.L. and R.J. Wagenet, 1992. LEACHM. Leaching Estimation and Chemistry Model: A Process Based Model of Water and Solute Movement, Transformation, Plant Uptake and Chemical Reactions in Unsaturated Zone. Version 3, Cornell University, Ithaca, New York.

7:  Igwe, C.A. and K. Stahr, 2004. Waterstable aggregates of flooded chemical properties. Aust. J. Soil Res., 42: 171-179.
Direct Link  |  

8:  Millington, R.J. and J.P. Quirk, 1961. Permeability of porous solids. Trans. Faraday Soc., 57: 1200-1207.
CrossRef  |  Direct Link  |  

9:  Mualem, Y., 1976. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resourc. Res., 12: 513-522.
CrossRef  |  Direct Link  |  

10:  Onweremadu, E.U., E.T. Eshett and G.E. Osuji, 2007. Temporal variability of selected heavy metals in automobile soil. Int. J. Environ. Sci. Technol., 4: 35-41.
Direct Link  |  

11:  Oti, N.N., 2002. Discriminate functions for classifying erosion degraded lands at Otamiri, Southeastern Nigeria. Agro-Sci., 3: 34-40.

12:  Pang, L., M.E. Close, J.P.C. Watt and K.W. Vincent, 2000. Simulation of picloram, atrazine, and simazine leaching through two new zealand soils and into groundwater using HYDRUS-2D. J. Contaminant Hydrol., 44: 19-46.
CrossRef  |  

13:  Richards, L.A., 1931. Capillary conduction of liquids through porous mediums. Physics, 1: 318-333.
CrossRef  |  Direct Link  |  

14:  RZWQM Team, 1995. RZWQM users manual. GPSR Technical Rep. No. 5, USDA-ARS, Great Plains Systems Res., Fort Collins, CO.

15:  Saleque, M.A., U.A. Naher, A. Islam, A.B.M.B.U. Pathan, A.T.M.S. Hossain and C.A. Meisner, 2004. Inorganic and organic phosphorus fertilizer effects on the phosphorus fractionation in wetland rice soils. Soil Sci. Soc. Am. J., 68: 1635-1644.
CrossRef  |  Direct Link  |  

16:  Simard, G., 2005. Monitoring and simulation of nutrient transport from agricultural fields. M.Sc. Thsies, Montreal, Quebec: McGill University, Department Biroresource Engineering.

17:  Simunek, J., M.Th. van Genuchten and M. Ejna, 2005. The HYDRUS-1D Software Package for Simulating One-Dimensional Movement of Water, Heat, and Multiple Solutes in Variably-Saturated Media. Version 3.0, Department of Environmental Sciences, University of California Riverside, USA., pp: 1-270.

18:  Simunek, J. and M.T. van Genuchten, 1994. The CHAIN-2D code for simulating the two-dimensional movement of water, heat and multiple solutes in variably-saturated porous media version 1.1. Research Report No. 136, US. Salinity Laboratory. USDA-ARS. Riverside, Califormia.

19:  Simunek, J., M. Sejna and M.T. van Genuchten, 1999. HYDRUS-2D/MESHGEN-2D, simulating water flow and solute transport in two-dimensional variably saturated media. US Salinity Laboratory, USDA/ARS, Riverside, California-distributed by International Ground Water Modeling Center, Colorado School of Mines, Golden, CO 80401, USA.

20:  Simunek, J., M. Seja and M.T. van Genuchten, 1999. The HYDRUS-2D Software Package for Simulating the Two-Dimensional Movement of Water, Heat and Multiple Solutes in Variably-Saturated Media. Version 2.0, USDA/ARS, California, USA.

21:  Stamm, C., H. Fluhler, R. Gachter, J. Leuenberger and H. Wunderli, 1998. Preferential transport of phosphorus in drained grassland soils. J. Environ. Qual., 27: 515-522.
Direct Link  |  

22:  Van Genuchten, M.T., 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J., 44: 892-898.
CrossRef  |  Direct Link  |  

23:  Van Genuchten, M.T., F.J. Leij and S.R. Yates, 1991. The RETC Code for Qualifying the Hydraulic Functions of Unsaturated Soils. USSL, USDA, Agricultural Research Services, California, USA.

24:  Ventrella, D., B.P. Mohanty, J. Imunek, N. Lasavio and M.T. van Genuchten, 2000. Water and chloride transport in a fine-textured soil: Field experimentand modeling. J. Soil Sci., 165: 624-631.
Direct Link  |  

25:  Zhang, R., 1997. CHAIN-IR: Irrigation simulations of water flow and solute transport with nitrogen transformation version 1.1. Report No. B-1961. Agricultural Experiment Station, University of Wyoming, Laramie, WY.

26:  Ben-Gal, A. and L.M. Dudley, 2003. Phosphorus availability under continuous point source irrigation. Soil Sci. Soc. Am. J., 67: 1449-1456.
Direct Link  |  

27:  Phillips, I.R., 2006. Modelling water and chemical transport in large undisturbed soil cores using HYDRUS-2D. Aust. J. Soil Res., 44: 27-34.
CrossRef  |  Direct Link  |  

28:  Rassam, D.W. and F.J. Cook, 2002. Numerical simulations of water flow and solute transport applied to acid sulfate soils. J. Irrig. Drain. Eng., 128: 107-115.
CrossRef  |  Direct Link  |  

©  2021 Science Alert. All Rights Reserved