|
|
|
|||||||||||||||||
| JOURNAL HOME | HELP | CONTACT PUBLISHER | SUBSCRIBE | ARCHIVE | SEARCH | TABLE OF CONTENTS |
ORIGINAL RESEARCH |
a Institute for Soil Landscape Research, Leibniz-Centre for Agricultural Landscape Research (ZALF), Eberswalder Strasse 84, D-15374 Müncheberg, Germany
b Dep. of Hydraulics and Hydrology, Faculty of Civil Engineering, Czech Technical Univ., Prague, Czech Republic
c Faculty for Agricultural and Environmental Sciences, Soil Physics and Resources Protection, Univ. of Rostock, Justus-von-Liebig-Weg 6, D-18059 Rostock, Germany
Correspondence: * Corresponding author (hgerke{at}zalf.de).
Received for publication 10 February 2007. Preferential flow has been hypothesized as an important factor for chemical leaching from tile-drained agricultural fields with structured soils originating from glacial till sediments. Previous studies showed that one-dimensional single-porosity models (1D-SPM) failed and that one-dimensional dual-permeability models (1D-DPERM) were limited in explaining both Br leaching and residual Br distribution, although tile water outflow peaks could somehow be reproduced. The objective of this paper was to analyze the tile outflow and leaching patterns using a two-dimensional (2D)-DPERM and a standard 2D-SPM for comparison. Flow and transport were simulated in a 2D vertical cross-section of 5.9 m length and 2 m depth using previously tested parameters. Simulated drainage rates and Br-effluent concentrations were made comparable with collector data from a field experiment by weighing results for irrigated and nonirrigated plots according to their area fractions. The 2D-DPERM simulations for surface application of Br in dissolved form in both domains overestimated the observed initial outflow concentration peaks, in contrast to closer approximation of observations assuming Br application in the soil matrix domain only. The simulated 2D mass transfer rate distribution showed most intensive exchange between domains near the water table and in the topsoil. Results from the 2D-DPERM analyses suggest that conditions at the soil surface, near the water table, and of the field-scale mixing are significantly affecting leaching patterns, in addition to local nonequilibrium effects. Here, the description of preferential flow toward tile drain could be strongly improved with the 2D-DPERM compared with the 2D-SPM. Further improvements remain challenging with respect to DPERM numerical modeling and field experimentation, with special attention toward soil structure and soil surface conditions.
Abbreviations: 1D, one-dimensional 2D, two-dimensional 3D, three-dimensional BC, boundary condition Br, bromide CDE, convection–dispersion equation DPERM, dual-permeability model DPM, dual-porosity model ET, evapotranspiration MIM, mobile–immobile model PF, preferential flow SM, soil matrix SPM, single-porosity model
This article has been cited by other articles:
![]() |
J. Dusek, H. H. Gerke, and T. Vogel Surface Boundary Conditions in Two-Dimensional Dual-Permeability Modeling of Tile Drain Bromide Leaching Vadose Zone J., November 1, 2008; 7(4): 1287 - 1301. [Abstract] [Full Text] [PDF] |
||||
| JOURNAL HOME | HELP | CONTACT PUBLISHER | SUBSCRIBE | ARCHIVE | SEARCH | TABLE OF CONTENTS |