Coupling water flow and solute transport into a physically-based surface–subsurface hydrological model

• Published in: Advances in water resources Vol. 34; no. 1; pp. 128 - 136
• Main Authors: Weill, S., Mazzia, A., Putti, M., Paniconi, C.
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier Ltd 2011; Elsevier
• Subjects: Advection–dispersion equation; Channels; Earth Sciences; Earth, ocean, space; Environmental Engineering; Environmental Sciences; Exact sciences and technology; Flow–transport coupling; Freshwater; Global Changes; Hydrology; Hydrology. Hydrogeology; Integrated hydrologic modeling; Joining; Land; Marine; Mathematical analysis; Mathematical models; Old water/new water paradox; Parametrization; Richards equation; Sciences of the Universe; Surface–subsurface interactions; Transport; Water flow; Surface–subsurface interactions; Flow–transport coupling; Advection–dispersion equation; Richards equation; Integrated hydrologic modeling; Old water/new water paradox; experimental studies; algorithms; mass balance; Surface-subsurface interactions; storage; Advection-dispersion equation; coupling; parametrization; water resources; boundary conditions; performances; flow; advection; models; solutes; runoff; Flow-transport coupling; transport; finite element analysis; subsurface; channels; standard samples; drainage; dispersion; strategy
• ISSN: 0309-1708; 1872-9657
• DOI: 10.1016/j.advwatres.2010.10.001

► CATHY surface–subsurface flow model extended to solute transport phenomena. ► Model features boundary condition-based coupling and time-splitting method. ► Suitable for field, hillslope, and catchment scale simulations. ► Model addresses atmosphere/land surface/river/soil/groundwater interactions. ► Applications include contamination, hydrograph separation, water balance estimation. A distributed-parameter physically-based solute transport model using a novel approach to describe surface–subsurface interactions is coupled to an existing flow model. In the integrated model the same surface routing and mass transport equations are used for both hillslope and channel processes, but with different parametrizations for these two cases. For the subsurface an advanced time-splitting procedure is used to solve the advection–dispersion equation for transport and a standard finite element scheme is used to solve Richards equation for flow. The surface–subsurface interactions are resolved using a mass balance-based surface boundary condition switching algorithm that partitions water and solute into actual fluxes across the land surface and changes in water and mass storage. The time stepping strategy allows the different time scales that characterize surface and subsurface water and solute dynamics to be efficiently and accurately captured. The model features and performance are demonstrated in a series of numerical experiments of hillslope drainage and runoff generation.