Simulation of Water Infiltration Process into a Porous Unsaturated Soil: Application on Tangier’s Bay Region-Morocco

• Published in: Geotechnical and geological engineering Vol. 38; no. 1; pp. 353 - 365
• Main Authors: Belcaid, Aicha, Benaicha, Mouhcine, Le Palec, Georges, Draoui, Abdeslam
• Format: Journal Article
• Language: English
• Published: Cham Springer International Publishing 2020; Springer Nature B.V
• Subjects: Aeration zone; Civil Engineering; Computer simulation; Computer software; Earth and Environmental Science; Earth Sciences; Geotechnical Engineering & Applied Earth Sciences; Groundwater; Hydrogeology; Infiltration; Mathematical models; Moisture content; Original Paper; Physics; Pressure head; Saturation; Simulation; Soil; Soil physics; Soil pollution; Soil porosity; Soil properties; Soil structure; Soil water; Terrestrial Pollution; Unsaturated soils; Waste Management/Waste Technology; Water content; Water infiltration; Water pollution; Water properties; Hydrus 2D/3D; Unsaturated zone; Hydraulic conductivity; Pressure head; Water content
• ISSN: 0960-3182; 1573-1529
• DOI: 10.1007/s10706-019-01024-7

The study is focused on numerical simulation of the infiltration process of polluted water originating from a wadi into the unsaturated soil of a bay. This study addresses the linkage between unsaturated zone and groundwater and it discusses the ability of the unsaturated soil to receive the water dried at once by the wadi and the bay. The simulation requires a number of trials, something that is experimentally expensive and sometimes impossible. The infiltration process is governed by the Richard’s equation. To solve this later, the knowledge of soil water properties such as pressure head, water content and hydraulic conductivity is needed. Hence, the water retention and unsaturated hydraulic conductivity curves near saturation have been established for a range of water content going from h = (− 25 m; 25 m) to h = (− 175 m; 175 m). These curves are obtained using the simulating computer program Hydrus 2D/3D. Numerical results have been calibrated with the Musy and Soutter (Physique du sol (soil physics), Université de Technologie Compiègne, Lausanne, 1991) type relationships. It has been explained that the infiltration process from the saturated zones to the unsaturated zone is relatively slow depending on the soil structure.