HAL will be down for maintenance from Friday, June 10 at 4pm through Monday, June 13 at 9am. More information
Skip to Main content Skip to Navigation
Conference papers

Using complementarity framework to couple subsurface flow and seepage processes: a physically based basis to integrate hotspots reactivity at the hillslope scale

Jean-Raynald De Dreuzy 1 Jean Marçais 1 Jocelyne Erhel 2
2 FLUMINANCE - Fluid Flow Analysis, Description and Control from Image Sequences
IRMAR - Institut de Recherche Mathématique de Rennes, IRSTEA - Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture, Inria Rennes – Bretagne Atlantique
Abstract : Landscape structure and geological heterogeneities are major controls on subsurface flow dynamic. Particularly it has strong impact on saturated areas (potential hotspots) emergence, promoting seepage production. This control leads to highly non linear flow response where, for each hillslope location, two possible states can be distinguished: with or without seepage. Different algorithmic solutions have been proposed to model this process. A solution is to solve subsurface flow for an assumed water table position and then iterate on the water table position until convergence is met. Seepage areas and seepage values are then deduced from the locations where water table intersects surface. A second way to proceed is to explicitly couple groundwater equations and surface water equations with an exchange flux. Here we developed a novel approach using the complementarity framework to reconcile in a single system the two states potentially encountered (with or without seepage). The complementarity framework manages the current state and the possible transition between states thanks to a specifically devoted equation. This framework is applied to the 1D hillslope storage Boussinesq equations (Troch et al. 2003). Reformulating this complementarity system differently enables to partition effectively the local flux balance between storage variation and seepage. This differential algebraic equations (DAEs) system has the major benefit to be directly solvable with built-in ode libraries. Finally the system is regularized to enhance fast and efficient solving. This model is stable with fast spatial convergence. It respects the mass balance locally beyond the tolerance limits and shows limited sensitivity to the value of the regularization parameter. The model appears to be robust, able to solve complex realistic case with presence of landscape heterogeneity and real hydrologic forcing. This model will then be used as a physical basis to implement biogeochemistry reactivity.
Document type :
Conference papers
Complete list of metadata

Contributor : Isabelle Dubigeon Connect in order to contact the contributor
Submitted on : Tuesday, February 14, 2017 - 9:10:41 AM
Last modification on : Friday, May 20, 2022 - 9:04:52 AM


  • HAL Id : insu-01467053, version 1


Jean-Raynald De Dreuzy, Jean Marçais, Jocelyne Erhel. Using complementarity framework to couple subsurface flow and seepage processes: a physically based basis to integrate hotspots reactivity at the hillslope scale . American Geophysical Union Fall Meeting 2016, American Geophysical Union, Dec 2016, San Francisco, United States. pp.H41H-03. ⟨insu-01467053⟩



Record views