Transport of moisture and solutes in the unsaturated zone by preferential flow

Author(s):  
Frank Stagnitti ◽  
J.-Yves Parlange ◽  
Tammo S. Steenhuis ◽  
Jan Boll ◽  
Bruce Pivetz ◽  
...  
2017 ◽  
Vol 21 (1) ◽  
pp. 459-471 ◽  
Author(s):  
Mostaquimur Rahman ◽  
Rafael Rosolem

Abstract. Modelling and monitoring of hydrological processes in the unsaturated zone of chalk, a porous medium with fractures, is important to optimize water resource assessment and management practices in the United Kingdom (UK). However, incorporating the processes governing water movement through a chalk unsaturated zone in a numerical model is complicated mainly due to the fractured nature of chalk that creates high-velocity preferential flow paths in the subsurface. In general, flow through a chalk unsaturated zone is simulated using the dual-porosity concept, which often involves calibration of a relatively large number of model parameters, potentially undermining applications to large regions. In this study, a simplified parameterization, namely the Bulk Conductivity (BC) model, is proposed for simulating hydrology in a chalk unsaturated zone. This new parameterization introduces only two additional parameters (namely the macroporosity factor and the soil wetness threshold parameter for fracture flow activation) and uses the saturated hydraulic conductivity from the chalk matrix. The BC model is implemented in the Joint UK Land Environment Simulator (JULES) and applied to a study area encompassing the Kennet catchment in the southern UK. This parameterization is further calibrated at the point scale using soil moisture profile observations. The performance of the calibrated BC model in JULES is assessed and compared against the performance of both the default JULES parameterization and the uncalibrated version of the BC model implemented in JULES. Finally, the model performance at the catchment scale is evaluated against independent data sets (e.g. runoff and latent heat flux). The results demonstrate that the inclusion of the BC model in JULES improves simulated land surface mass and energy fluxes over the chalk-dominated Kennet catchment. Therefore, the simple approach described in this study may be used to incorporate the flow processes through a chalk unsaturated zone in large-scale land surface modelling applications.


2017 ◽  
Vol 65 (3) ◽  
pp. 287-296 ◽  
Author(s):  
Nicola Pastore ◽  
Claudia Cherubini ◽  
Concetta I. Giasi

AbstractWhen dealing with groundwater resources, a better knowledge of the hydrological processes governing flow in the unsaturated zone would improve the assessment of the natural aquifer recharge and its vulnerability to contamination. In North West Europe groundwater from unconfined chalk aquifers constitutes a major water resource, therefore the need for a good hydrological understanding of the chalk unsaturated zone is essential, as it is the main control for aquifer recharge. In the North Paris Basin, much of the recharge must pass through a regional chalk bed that is composed of a porous matrix with embedded fractures. The case study regards the role of the thick unsaturated zone of the Cretaceous chalk aquifer in Picardy (North of France) that controls the hydraulic response to rainfall. In order to describe the flow rate that reaches the water table, the kinematic diffusion theory has been applied that treats the unsaturated water flow equation as a wave equation composed of diffusive and gravitational components. The kinematic diffusion model has proved to be a convenient method to study groundwater recharge processes in that it was able to provide a satisfactory fitting both for rising and falling periods of water table fluctuation. It has also proved to give an answer to the question whether unsaturated flow can be described using the theory of kinematic waves. The answer to the question depends principally on the status of soil moisture. For higher values of hydraulic Peclet number (increasing saturation), the pressure wave velocities dominate and the preferential flow paths is provided by the shallow fractures in the vadose zone. With decreasing values of hydraulic Peclet number (increasing water tension), rapid wave velocities are mostly due to the diffusion of the flow wave. Diffusive phenomena are provided by matrix and fracture-matrix interaction.The use of a kinematic wave in this context constitutes a good simplified approach especially in cases when there is a lack of information concerning the hydraulic properties of the fractures/macropores close to saturation.


2009 ◽  
Vol 13 (7) ◽  
pp. 1215-1233 ◽  
Author(s):  
T. Blume ◽  
E. Zehe ◽  
A. Bronstert

Abstract. Spatial patterns as well as temporal dynamics of soil moisture have a major influence on runoff generation. The investigation of these dynamics and patterns can thus yield valuable information on hydrological processes, especially in data scarce or previously ungauged catchments. The combination of spatially scarce but temporally high resolution soil moisture profiles with episodic and thus temporally scarce moisture profiles at additional locations provides information on spatial as well as temporal patterns of soil moisture at the hillslope transect scale. This approach is better suited to difficult terrain (dense forest, steep slopes) than geophysical techniques and at the same time less cost-intensive than a high resolution grid of continuously measuring sensors. Rainfall simulation experiments with dye tracers while continuously monitoring soil moisture response allows for visualization of flow processes in the unsaturated zone at these locations. Data was analyzed at different spacio-temporal scales using various graphical methods, such as space-time colour maps (for the event and plot scale) and binary indicator maps (for the long-term and hillslope scale). Annual dynamics of soil moisture and decimeter-scale variability were also investigated. The proposed approach proved to be successful in the investigation of flow processes in the unsaturated zone and showed the importance of preferential flow in the Malalcahuello Catchment, a data-scarce catchment in the Andes of Southern Chile. Fast response times of stream flow indicate that preferential flow observed at the plot scale might also be of importance at the hillslope or catchment scale. Flow patterns were highly variable in space but persistent in time. The most likely explanation for preferential flow in this catchment is a combination of hydrophobicity, small scale heterogeneity in rainfall due to redistribution in the canopy and strong gradients in unsaturated conductivities leading to self-reinforcing flow paths.


2014 ◽  
Vol 2 (1) ◽  
pp. 181-195 ◽  
Author(s):  
D. M. Krzeminska ◽  
T. A. Bogaard ◽  
T.-H. Debieche ◽  
F. Cervi ◽  
V. Marc ◽  
...  

Abstract. The unsaturated zone largely controls groundwater recharge by buffering precipitation while at the same time providing preferential flow paths for infiltration. The importance of preferential flow on landslide hydrology is recognised in the literature; however, its monitoring and quantification remain difficult. This paper presents a combined hydrological and hydrochemical analysis of small-scale sprinkling experiments. It aims at showing the potential of such experiments for studying the spatial differences in dominant hydrological processes within a landslide. This methodology was tested in the highly heterogeneous black marls of the Super-Sauze landslide. The tests were performed in three areas characterised by different displacement rates, surface morphology and local hydrological conditions. Special attention was paid to testing the potential of small-scale sprinkling experiments for identifying and characterising preferential flow patterns and dominant hydrological processes.


1997 ◽  
Vol 28 (4-5) ◽  
pp. 307-328 ◽  
Author(s):  
Nils-Otto Kitterød ◽  
E. Langsholt ◽  
W. K. Wong ◽  
L. Gottschalk

The spatial distribution of soil moisture defines preferential flow paths in the unsaturated zone. Hence, three dimensional (3D) estimates of soil moisture are of great importance to understand transport of contaminants as well as remediation processes in the unsaturated zone. In this study 3D estimates conditioned on spatially frequent observations of soil moisture, have been obtained by kriging. The observations were divided into subdomains consistent with the local stratigraphy and directional semivariogram analysis was applied. It was found difficult to clearly identify a 3D semivariogram function in this case, but from a georadar survey two semivariogram functions were derived, describing two different sedimentological units. By conditioning the estimates of soil moisture on the sedimentological architecture computed by indicator kriging, more accurate estimates were achieved. These improvements were quantified by a ‘jackknife’ cross validation procedure. Besides the practical aspects of finding the most important flow paths estimates of soil moisture are valuable when validating unsaturated flow models.


Sign in / Sign up

Export Citation Format

Share Document