Numerical Solution of the Advection-Dispersion Equation: Application to the Agricultural Drainage

Authors
C. Chavez 1
C. Fuentes 2
F. Brambila 3
A. Castañeda 4
1 Universidad Autonoma de Queretaro, Facultad de Ingenieria, Cerro de las Campanas s/n, 76010, Santiago de Querétaro, Qro., México.
2 Universidad Autónoma de Querétaro, Facultad de Ingeniería, Cerro de las Campanas s/n, 76010, Santiago de Querétaro, Qro., México
3 Universidad Nacional Autónoma de México, Facultad de Ciencias, Ciudad Universitaria, 07300, México D.F.
4 Universidad Autónoma de Querétaro. Facultad de Ingeniería, Cerro de las Campanas s/n, 76010, Santiago de Querétaro, Qro., México.
Abstract
Subsurface drainage systems are used to control the depth of the water table and to reduce or prevent soil salinity. Water flow in these systems is described by the Boussinesq Equation, and the Advection-Dispersion Equation coupled with the Boussinesq Equation is used to study the solute transport. The objective of this study was to propose a finite difference solution of the Advection-Dispersion Equation using a lineal radiation condition in the drains. The equations’ parameters were estimated from a methodology based on the granulometric curve and inverse problems. The algorithm needs the water flow values, which were calculated with the Boussinesq Equation, where a fractal radiation condition and variable drainable porosity were applied. To evaluate the solution descriptive capacity, a laboratory drainage experiment was used. In the experiment, the pH, temperature, and electric conductivity of drainage water were measured to find the salt’s concentration. The salts concentration evolution was reproduced using the finite difference solution of the Advection-Dispersion Equation, and the dispersivity parameter was found by inverse modelling. The numerical solution was used to simulate the leaching of saline soil. The result showed that this solution could be used as a new tool for the design of agricultural drainage systems, enabling the optimal development of crops according to their water needs and the degree of tolerance to salinity.

Keywords

1. Abassi, F., Simunek, J., van Genuchten, M. Th., Feyen, J., Adamsen, F. J., Hunsaker, D. J., Strelkoff, T. S. and Shouse, P. 2003. Overland Flow and Solute Transport: Model Development and Field-data Analysis. J. Irrig. Drain. Eng., 129 (2): 71–81.
2. Boussinesq, J. 1904. Recherches The´oriques sur L’e´coulement des Nappes d’eau Infiltre´es Dans le Sol et Sur le De´bit des Sources. J. Math. Pure. Appl., 10: 5-78.
3. Carslaw, H. S. and Jaeger, J. C. 1959. Conduction of Heat in Solids. Oxford University Press, Oxford, 507 pp.
4. Chávez, C. 2010. Solución Mumérica de las Ecuaciones de Transferencia de Agua y Solutos en Riego y Drenaje. Dr. in Eng. Thesis, Universidad Autónoma de Querétaro México, 200 PP. (in Spanish)
5. Chávez, C., Fuentes, C., Zataráin, F. and Zavala, M. 2011. Finite Difference Solution of the Boussinesq Equation with Variable Drainable Porosity and Fractal Radiation Boundary Condition. Agrociencia, 45(8): 911-927.
6. CONAGUA. 2010. Statistics on Water in Mexico. México, D. F., PP. 17-117. (in Spanish)
7. Conte, S. D. and De Boor, C. 1980. Elementary Numerical Analysis: An Algorithmic Approach. McGraw-Hill, New York, 432 PP.
8. Elder, J. W. 1959. The Dispersion of Marked Fluid in Turbulent Shear Flow. J. Fluid Mech., 5: 544-560.
9. Fayer, M. J. 2000. UNSAT-H Version 3.0: Unsaturated Soil Water and Heat Flow Model. Theory, User Manual, and Examples. Pacific Northwest National Laboratory 13249, USA, 184 PP.
10. Feyen, J., Jacques, D., Timmerman, A. and Vanderborght, J. 1998. Modelling Water Flow and Solute Transport in Heterogeneous Soils: A Review of Recent Approaches. J. Agric. Eng. Res., 70: 231-256.
11. Fischer H. B. 1967. The Mechanics of Dispersion in Natural Streams. J. Hydraul. Div., Am. Soc. Civ. Eng., 93 (6): 187-216.
12. Freund, R. W. and Hoppe, R. W. 2007. Stoer/Bulirsch: Numerische Mathematik 1. Springer-Lehrbuch, Germany, 410 PP.
13. Fuentes, C. 1992. Approche Fractale des Transferts Hydriques Dans les Sols No-Saturés. Tesis de Doctorado, Universidad Joseph Fourier de Grenoble, Francia, 276 PP.
14. Fuentes, C., Brambila, F., Vauclin, M., Parlange, J-Y. and Haverkamp, R. 2001. Fractal Modeling of Hydraulic Conductivity in Non-Saturated soils. Hydraul. Eng. México, 16(2): 119-137.
15. Fuentes, C., Zavala, M. and Saucedo, H. 2009. Relationship between the Storage Coefficient and the Soil-water Retention Curve in Subsurface Agricultural Drainage Systems: Water Table Drawdown. J. Irrig. Drain. Eng., 135(3): 279-285.
16. Haverkamp, R., Leij, F. J., Fuentes, C., Sciortino, A. and Ross P. J. 2005. Soil Water Retention. I. Introduction of a Shape Index. Soil Sci. Soc. Am. J., 69: 1881-1890.
17. Holly, F. M. 1985. Dispersion in Rivers and Coastal Waters. 1. Physical Principles and Dispersion Equations. 3. In: “Developments in Hydraulic Engineering”, (Ed.): Novak, P.. Elsevier, New York, 1: 1-38.
18. Holly, F. M. 1975. Two Dimensional Mass Dispersion in Rivers. Hydrologic Papers, Colorado State University Press, Fort Collins, Colorado, 78 PP.
19. IMTA. 1998. Manual de Diseño e Instalación de Drenaje Parcelario en Zonas áridas y Semiáridas bajo Riego. México, 71 PP.
20. Jansson, P. E. and Karlberg, L. 2001. Coupled Heat and Mass Transfer Model for Soil-plant-atmosphere Systems. Royal Institute of Technology, Department of Civil and Environmental Engineering, Stockholm, 325 PP.
21. Marquardt, D. W. 1963. An Algorithm for Least-squares Estimation of Nonlinear Parameters. SIAM J. Appl. Math., 11: 431-441.
22. Mirabzadeh, M. and Mohammadi, K. 2006. A Dynamic Programming Solution to Solute Transport and Dispersion Equations in Groundwater. J. Agric. Sci. Technol., 8: 233-241.
23. Mousavi1, S. F., Mostafazadeh-Fard, B., Farkhondeh, A. and Feizi, M. 2009. Effects of Deficit Irrigation with Saline Water on Yield, Fruit Quality and Water Use Efficiency of Cantaloupe in an Arid Region. J. Agr. Sci. Tech., 11: 469-479.
24. Richards, L. A. 1931. Capillary Conduction of Liquids through Porous Mediums. Physics, 1: 318-333.
25. Ritzema, H. P., Satyanarayama, T. V., Raman, S. and Boonstra, J. 2008. Subsurface Drainage to Combat Waterlogging and Salinity in Irrigated Lands in India: Lessons Learned in Farmers’ Fields. Agric. Water Manage., 95(3): 179-189.
26. Ruiz-Cerda, E., Aldaco, N. R., Montemayor, T. J., Fortis, H. M., Olague, R. J. and Villagómez, G. J. 2007. Aprovechamiento y Mejoramiento de un Suelo Salino Mediante Vermicomposta. Tec. Pecu. Mex., 45(1): 19-24.
27. Rutherford, J. C. 1994. River Mixing. Wiley, New York, 347 PP.
28. Simunek, J. 2005. Models of Water Flow and Solute Transport in the Unsaturated Zone. “Encyclopedia of Hydrological Sciences”. (Ed.): Anderson, John Wiley & Sons, Ltd., Chichester, England, 1171-1180.
29. Simunek, J. and van Genuchten, M. Th. 1999. Using the HYDRUS-1D and HYDRUS-2D Codes for Estimating Unsaturated Soil Hydraulic and Solute Transport Parameters. In: “Characterization and Measurement of the Hydraulic Properties of Unsaturated Porous Media” (Eds.): van Genuchten, M. Th. Simunek, J. and "Sejna, M. University of California, Riverside, CA, PP. 1523–1536.
30. Simunek, J., Sejna, M. and van Genuchten, M. Th. 1998. The HYDRUS-1D Software Package for Simulating the Movement of Water, Heat, and Multiple Solutes in Variably Saturated Media, Version 2.0. United States Salinity Laboratory, USDA-ARS, Riverside, Calif.
31. Simunek, J., van Genuchten M.Th, Sejna, M., Toride, N. and Leij, F.J. 1999. The STANMOD Computer Software for Evaluating Solute Transport in Porous Media Using Analytical Solutions of Convection-Dispersion Equation, Versions 1.0 and 2.0. IGWMC – TPS-71, International Ground Water Modeling Center, Colorado School of Mines, Golden, 32 PP.
32. Singh, S., Ghosh, N. C., Pandey, R. P., Galkate, R. V., Thomas, T. and Jaiswal, R. K. 2009. Numerical Solution of 1D Boussinesq Equation for Water Table Fluctuation between Drains in Response to Recharge and ET in A Sloping Aquifer. Int. J. Eco. Econ. Stat., 14(9): 45-54.
33. Taylor, G. I. 1954. The Dispersion of Matter in Turbulent Flow through a Pipe. Proc. R. Soc. London, Ser. A. 223: 446–48.
34. van Dam, J. C., Huygen, J., Wesseling, J. G., Feddes, R. A., Kabat, P., van Valsum, P. E. V., Groenendijk, P. and van Diepen, C. A. 1997. Theory of SWAP, Version 2.0. Simulation of Water Flow, Solute Transport and Plant Growth in the Soil-water-atmosphere-plant Environment. Technical Document 45, Report 71, Department Water Resources, WAU, DLO Winand Staring Centre, Wageningen.
35. van Genuchten, M. 1980. A Closed-Form Equation for Predicting the Hydraulic Conductivity of the Unsaturated Soils. Soil Sci. Soc. Amer. J., 44: 892-898.
36. Verhoest, N., Pauwels, V., Troch, P. and de Troch, F. 2002. Analytical Solution for Transient Water Table Heights and Outflows from Inclined Ditch-drained Terrains. J. Irrig. Drain Eng., 128(6): 358–364.
37. Villareal, E. and Bello, S. 1964. The Concentration and Electrical Conductivity in Aqueous Solutions of Electrolytes. Rev. Mex. Fis., 13(2): 55-74.
38. Zataráin, F., Fuentes, C., Palacios, V. O. L., Mercado, E., Brambila, F. and Villanueva, N. 1998. Modelación del Transporte de Agua y Solutos en el Suelo Agrociencia, 32(4): 373-383. (in Spanish)
39. Zavala, M., Fuentes, C. and Saucedo, H. 2007. Nonlinear Radiation in the Boussinesq Equation of Agricultural Drainage. J. Hydrol., 332(3): 374-380.
40. Zerihun, D., Furman, A., Warrick, A. W. and Sánchez, C. A. 2005. Coupled Surface-Subsurface Solute Transport Model for Irrigation Borders and Basin. I. Model Development. J. Irri. Drain. Eng., 131(3): 396-406.
Journal of Agricultural Science and Technology
Volume 16, Issue 6
November and December 2014
Pages 1357-1388