Open AccessAn Improved Dual‐Permeability Model of Water Flow and Solute Transport in the Vadose Zone
Author(s) -
Larsbo Mats,
Roulier Stephanie,
Stenemo Fredrik,
Kasteel Roy,
Jarvis Nicholas
Publication year - 2005
Publication title -
vadose zone journal
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.036
H-Index - 81
ISSN - 1539-1663
DOI - 10.2136/vzj2004.0137
Subject(s) - vadose zone , hydraulic conductivity , macropore , richards equation , permeability (electromagnetism) , porosity , porous medium , mechanics , water flow , soil science , soil water , chemistry , geotechnical engineering , geology , physics , mesoporous material , biochemistry , membrane , catalysis
We introduce an improved, one‐dimensional, non‐steady‐state dual‐permeability model (MACRO 5.1). The model simulates water flow and solute transport in the vadose zone of structured soils by coupling a high‐conductivity–low porosity macropore domain to a low‐conductivity–high porosity domain representing the soil matrix. Mass exchange between the domains is approximated by first‐order expressions. The numerical solutions are briefly described, focusing on the dual‐permeability formulation. The solution method for water flow in macropores was verified by comparing simulation results with analytical solutions for a “kinematic wave”. The model was tested against high time‐resolution measurements of water flow and nonreactive (Cl − ) solute transport in transient microlysimeter experiments. The objective was to test the identifiability of four key model parameters determining the degree of preferential flow using the generalized likelihood uncertainty estimation (GLUE) procedure. The parameters were chosen either because they are difficult or impossible to measure directly or because they were considered sensitive on the basis of earlier experience with the model. The measurements, indicating strong preferential flow, were adequately reproduced by the model simulations (overall model efficiency = 0.62). The GLUE procedure conditioned the saturated matrix hydraulic conductivity, the macroporosity, and the mass exchange coefficient (diffusion pathlength), indicating that these parameters would be identifiable in inverse modeling approaches based on microlysimeter experiments. The conditioning of the kinematic exponent was poor, which was attributed primarily to correlation with the macroporosity.