Single- and dual-porosity modelling of flow in reclaimed mine soil cores with embedded lignitic fragments

Lignitic mine soils represent a typical two-scale dual-porosity medium consisting of a technogenic mixture of overburden sediments that include lignitic components as dust and as porous fragments embedded within a mostly coarse-textured matrix. Flow and transport processes in such soils are not sufficiently understood to predict the course of soil reclamation or of mine drainage. The objective of this contribution is to identify the most appropriate conceptual model for describing small-scale heterogeneity effects on flow on the basis of the physical structure of the system. Multistep flow experiments on soil cores are analyzed using either mobile–immobile or mobile–mobile type 1D dual-porosity models, and a 3D numerical model that considers a local-scale distribution of fragments. Simulations are compared with time series' of upward infiltration and matric potential heads measured at two depths using miniature tensiometers. The 3D and the 1D dual-permeability models yielded comparable results as long as pressure heads are in local equilibrium; however, could describe either the upward infiltration or the matric potential curves but not both at the same time. The mobile–immobile type dual-porosity model failed to describe the data. A simultaneous match with pressure heads and upward infiltration data could only be obtained with the 1D dual-permeability model (i.e., mobile–mobile) by assuming an additional restriction of the inter-domain water transfer. These results indicate that for unsaturated flow conditions at higher matric potential heads (i.e., here >− 40 hPa), water in a restricted part of the fragment domain must be more mobile as compared to water in the sandy matrix domain. Closer inspections of the pore system and first neutron radiographic imaging support the hypothesis that a more continuous pore region exists at these pressure heads in the vicinity of the lignitic fragments possibly formed by fragment contacts and a lignitic dust interface-region between the two domains. The results suggest that the small-scale structure is too complex as to be represented by weighted contributions of individual components alone.