Simulation of water and solute transport with MACRO model in Cecil loamy sand soil
Hasan Merdun A C and Virgil L. Quisenberry BA Department of Agricultural Engineering, Faculty of Agriculture, Kahramanmaraş Sütçü İmam University, Kahramanmaraş 46060, Turkey.
B Department of Crop and Soil Environmental Sciences, Clemson University, Clemson, SC 29634, USA.
C Corresponding author. Email: panter@ksu.edu.tr or hmerdun@postaci.com
Australian Journal of Soil Research 42(8) 939-951 https://doi.org/10.1071/SR03126
Submitted: 13 August 2003 Accepted: 11 August 2004 Published: 14 December 2004
Abstract
Modelling preferential flow increases understanding of flow and transport processes in the unsaturated (vadose) zone and, hence, helps prevention of groundwater contamination. A dual-porosity model, MACRO, was evaluated for long-term drainage flow and short-term chloride-tagged water flow simulations in well-structured Cecil loamy sand soil. Water flow in micropores is calculated by the Richards’ equation, and simple gravity flow is assumed in the macropores. Solute transport in the micropores is calculated by the convection–dispersion equation (CDE), and the dispersion and diffusion in the CDE is neglected for the solute transport in the macropores. Based on the statistical criteria, the model accurately simulated drainage flow with depth and time. The average values of 3 statistical parameters (coefficient of residual mass, model efficiency, correlation coefficient) for drainage flow of different plots and times were 0.0057, 0.972, and 0.987, respectively. Similarly, the average values of the 3 statistical parameters in the same order for water flow and chloride transport of 4 plots were 0.097, 0.628, and 0.915, and 0.167, 0.938, and 0.982, respectively. The model simulated long-term drainage flow better than the short-term applied chloride-tagged water. The percentage recovery of measured water and chloride 2 h after the application ceased was 83 and 63 in the 1.05-m-deep profile of plot 1. This was a strong indication of preferential flow in this soil. Two-domain flow concept was required for acceptable simulation of this type of flow. In the 2-domain flow, boundary tension, boundary hydraulic conductivity, and effective diffusion path-length were the 3 most important parameters controlling flow and transport between the 2 domains. The effective diffusion path-length represented the structural development with depth in Cecil loamy sand soil. The relationships between the variability in flow and transport characteristics and fundamental soil properties and, hence, the associated key model parameters suggest that pedotransfer functions can be developed for the estimation of dual-porosity model parameters that control preferential flow.
Additional keywords: macropore flow, solute transport, modeling, macro, 2-domain, structured soil.
Andreu L,
Moreno F,
Jarvis NJ, Vachaud G
(1994) Application of the model MACRO to water movement and salt leaching in drained and irrigated marsh soils, Marismas, Spain. Agricultural Water Management 25, 71–88.
| Crossref | GoogleScholarGoogle Scholar |
Armstrong A,
Aden K,
Amraoui N,
Diekkruger B,
Jarvis NJ,
Mouvet C,
Nicholls P, Wittwer C
(2000) Comparison of the performance of the pesticide-leaching models on a cracking clay soil: results using the Brimstone Farm dataset. Agricultural Water Management 44, 85–104.
| Crossref | GoogleScholarGoogle Scholar |
Beven K, Germann P
(1981) Water flow in soil macropores II. A combined flow model. Journal of Soil Science 32, 15–29.
Booltink HWG,
Hatano R, Bouma J
(1993) Measurement and simulation of bypass flow in a structured clay soil: A physico-morphological approach. Journal of Hydrology 148, 149–168.
| Crossref | GoogleScholarGoogle Scholar |
Brooks RH, Corey AT
(1964) Hydraulic properties of porous media. Hydrology papers No. 3. Colorado State University, Fort Collins, CO, USA.
Bruce RR, Dane JH, Quisenberry VL, Powell NL, Thomas AW
(1983)
Durner, W (1992). Predicting the unsaturated hydraulic conductivity using multiporosity water retention curves. ‘Indirect methods for estimating the hydraulic properties of unsaturated soils’. pp. 185–202. (University of California: Riverside, CA)
Durner W
(1994) Hydraulic conductivity estimation for soils with heterogeneous pore structure. Water Resources Research 30, 211–223.
| Crossref | GoogleScholarGoogle Scholar |
Flury M
(1996) Experimental evidence of transport of pesticides through field soils: A review. Journal of Environmental Quality 25, 558–566.
van Genuchten MTh, Dalton FN
(1986) Models for simulating salt movement in aggregated field soils. Geoderma 38, 165–183.
| Crossref | GoogleScholarGoogle Scholar |
van Genuchten MTh, Wierenga CR
(1976) Movement of solutes in soil. ‘Soil chemistry-physico-chemical models, Part B’. (Ed. GH Bolt)
pp. 349–386. (Elsevier: New York)
Gerke HH, van Genuchten MT
(1993a) A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media. Water Resources Research 29, 305–319.
| Crossref | GoogleScholarGoogle Scholar |
Gerke HH, van Genuchten MT
(1993b) Evaluation of a first-order water transfer term for variably saturated dual-porosity flow models. Water Resources Research 29, 1225–1238.
| Crossref | GoogleScholarGoogle Scholar |
Germann P, Beven K
(1985) Kinematic wave approximation to infiltration into soils with sorbing macropores. Water Resources Research 21, 990–996.
Gottesburen B,
Aden K,
Barlund I,
Brown C,
Dust M,
Gorlitz G,
Jarvis NJ,
Rekolainen S, Schafer H
(2000) Comparison of the pesticide leaching models: Results using the Weiherbach dataset. Agricultural Water Management 44, 153–181.
| Crossref | GoogleScholarGoogle Scholar |
Hatfield MW
(1988) Water and anion movement in a typic Hapludult. PhD thesis, Clemson University, USA.
Jabro JD,
Jemison JM,
Fox RH, Fritton DD
(1994) Predicting bromide leaching under field conditions using SLIM and MACRO. Journal of Soil Science 157, 215–223.
Jarvis NJ
(1998) Modeling the impact of preferential flow on nonpoint source pollution. ‘Physical non-equilibrium in soils: modeling and application’. (Eds HH Selim, L Ma)
pp. 195–221. (Ann Arbor Press: Chelsea, MI)
Jarvis NJ,
Bergstrom L, Dik PE
(1991b) Modeling water and solute transport in macroporous soil. II. Chloride breakthrough under non-steady flow. Journal of Soil Science 42, 71–81.
Jarvis NJ,
Johnsson PE,
Dik PE, Messing I
(1991a) Modeling water and solute transport in macroporous soil. I. Model description and sensitivity analysis. Journal of Soil Science 42, 59–70.
Jarvis NJ, Larsson MH
(1998) The MACRO model (Version 4.1). Technical description, http://130.238.110.134:80/bgf/Macrohtm/macro.htm
Jarvis NJ,
Stahli M,
Bergstrom LF, Johnson H
(1994) Simulation of dichlorprop and bentazon leaching in soils of contrasting texture using the MACRO model. Journal of Environmental Science and Health A29, 1255–1277.
Jarvis NJ,
Villholth KG, Ulen B
(1999) Modeling particle mobilization and leaching in macroporous soil. European Journal of Soil Science 50, 621–632.
| Crossref | GoogleScholarGoogle Scholar |
Larsson MH, Jarvis NJ
(1999a) A dual-porosity model to quantify macropore flow effects on nitrate leaching. Journal of Environmental Quality 28, 1298–1307.
Larsson MH, Jarvis NJ
(1999b) Evaluation of a dual-porosity model to predict field-scale solute transport in a macroporous soil. Journal of Hydrology 215, 153–171.
| Crossref | GoogleScholarGoogle Scholar |
Larsson MH,
Jarvis NJ,
Torstensson G, Kasteel R
(1999) Quantifying the impact of preferential flow on solute transport to tile drains in a sandy field soil. Journal of Hydrology 215, 116–134.
| Crossref | GoogleScholarGoogle Scholar |
Loague K, Green RE
(1991) Statistical and graphical methods for evaluating solute transport models: overview and application. Journal of Contaminant Hydrology 7, 51–73.
| Crossref | GoogleScholarGoogle Scholar |
Mohanty BP,
Bowman RS,
Hendricks JMH, van Genuchten MTh
(1997) New piecewise-continuous functions for modeling preferential flow in an intermittent-flood irrigated field. Water Resources Research 33, 2049–2063.
| Crossref | GoogleScholarGoogle Scholar |
Mualem Y
(1976) A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resources Research 12, 513–522.
Nelson PA
(1990) Solute movement and root distribution as affected by soil structure. PhD thesis, Clemson University, USA.
Quisenberry VL, Philips RE
(1976) Percolation of surface applied water in the field. Soil Science Society of American Proceedings 40, 484–489.
Quisenberry VL,
Philips RE, Zeleznik JM
(1994) Spatial distribution of water and chloride macropore flow in a well-structured soil. Soil Science Society of America Journal 58, 1294–1300.
Quisenberry VL,
Smith BR,
Philips RE,
Scott HD, Nortcliff S
(1993) A soil classification system for describing water and chemical transport. Journal of Soil Science 156, 306–315.
Saxena RK,
Jarvis NJ, Bergstrom L
(1994) Interpreting non-steady state tracer breakthrough experiments in sand and clay soil using a dual-porosity model. Journal of Hydrology 162, 279–298.
| Crossref | GoogleScholarGoogle Scholar |
Skopp J,
Gardner WR, Tyler EJ
(1981) Solute movement in structured soils: Two-region model with small interaction. Soil Science Society of America Journal 45, 837–842.
Thomas GW, Swoboda AR
(1970) Anion exclusion effects on chloride movement in soils. Journal of Soil Science 110, 163–166.
Villholth KG, Jensen KH
(1998) Flow and transport processes in a macroporous subsurface-drained glacial till soil. II: Model analysis. Journal of Hydrology 207, 121–135.
| Crossref | GoogleScholarGoogle Scholar |
Zurmuhl T, Durner W
(1998) Determinations of parameters for bimodal hydraulic functions by inverse modeling. Soil Science Society of America Journal 62, 874–880.
