DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences

Abstract

We explore the coupling of surface and subsurface flows on fully unstructured meshes that conform to complex soil structures. To accommodate the distorted meshes that inevitably result from explicit representation of complex soil structures, we leverage the structure of the Mimetic Finite Difference (MFD) spatial discretization scheme to couple surface and subsurface flows. The MFD method achieves second-order accuracy and maintains local mass conservation on distorted meshes. We couple the diffusion wave approximation for surface flows to the Richards equation for subsurface flow, ensuring continuity of both pressure and flux between the surface and subsurface. The MFD method is particularly convenient for this coupling because it uses face-based constraints in the subsurface system that can be expressed as face-pressure unknowns. Those unknowns are coincident with surface cell-based unknowns, thus allowing the discrete surface system to be directly substituted into the subsurface system and solved implicitly as a global system. Robust representation of the transition between wet and dry surface conditions requires upwinding of the relative permeability and is facilitated by globalization in the nonlinear solver. The approach and its implementation in the Advanced Terrestrial Simulator (ATS) are evaluated by comparison to previously published benchmarks. Using runoff from soils with patchymore » groundcover (duff) as an example, we show that the new method converges significantly faster in mesh convergence tests than the commonly used two-point flux approximation.« less

Authors:
; ; ; ; ; ; ;
Publication Date:
Research Org.:
Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1644243
Alternate Identifier(s):
OSTI ID: 1649323
Grant/Contract Number:  
LDRD-8872; LDRD-20150397DR; LDRD201200068DR; AC05-00OR22725
Resource Type:
Published Article
Journal Name:
Advances in Water Resources
Additional Journal Information:
Journal Name: Advances in Water Resources Journal Volume: 144 Journal Issue: C; Journal ID: ISSN 0309-1708
Publisher:
Elsevier
Country of Publication:
United Kingdom
Language:
English
Subject:
54 ENVIRONMENTAL SCIENCES; Mimetic finite differences; Integrated hydrology; Complex stratigraphy

Citation Formats

Coon, Ethan T., Moulton, J. David, Kikinzon, Evgeny, Berndt, Markus, Manzini, Gianmarco, Garimella, Rao, Lipnikov, Konstantin, and Painter, Scott L. Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences. United Kingdom: N. p., 2020. Web. doi:10.1016/j.advwatres.2020.103701.
Coon, Ethan T., Moulton, J. David, Kikinzon, Evgeny, Berndt, Markus, Manzini, Gianmarco, Garimella, Rao, Lipnikov, Konstantin, & Painter, Scott L. Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences. United Kingdom. https://doi.org/10.1016/j.advwatres.2020.103701
Coon, Ethan T., Moulton, J. David, Kikinzon, Evgeny, Berndt, Markus, Manzini, Gianmarco, Garimella, Rao, Lipnikov, Konstantin, and Painter, Scott L. Thu . "Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences". United Kingdom. https://doi.org/10.1016/j.advwatres.2020.103701.
@article{osti_1644243,
title = {Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences},
author = {Coon, Ethan T. and Moulton, J. David and Kikinzon, Evgeny and Berndt, Markus and Manzini, Gianmarco and Garimella, Rao and Lipnikov, Konstantin and Painter, Scott L.},
abstractNote = {We explore the coupling of surface and subsurface flows on fully unstructured meshes that conform to complex soil structures. To accommodate the distorted meshes that inevitably result from explicit representation of complex soil structures, we leverage the structure of the Mimetic Finite Difference (MFD) spatial discretization scheme to couple surface and subsurface flows. The MFD method achieves second-order accuracy and maintains local mass conservation on distorted meshes. We couple the diffusion wave approximation for surface flows to the Richards equation for subsurface flow, ensuring continuity of both pressure and flux between the surface and subsurface. The MFD method is particularly convenient for this coupling because it uses face-based constraints in the subsurface system that can be expressed as face-pressure unknowns. Those unknowns are coincident with surface cell-based unknowns, thus allowing the discrete surface system to be directly substituted into the subsurface system and solved implicitly as a global system. Robust representation of the transition between wet and dry surface conditions requires upwinding of the relative permeability and is facilitated by globalization in the nonlinear solver. The approach and its implementation in the Advanced Terrestrial Simulator (ATS) are evaluated by comparison to previously published benchmarks. Using runoff from soils with patchy groundcover (duff) as an example, we show that the new method converges significantly faster in mesh convergence tests than the commonly used two-point flux approximation.},
doi = {10.1016/j.advwatres.2020.103701},
journal = {Advances in Water Resources},
number = C,
volume = 144,
place = {United Kingdom},
year = {Thu Oct 01 00:00:00 EDT 2020},
month = {Thu Oct 01 00:00:00 EDT 2020}
}

Works referenced in this record:

Mathematical Modeling of Watershed Hydrology
journal, July 2002


Weighted Implicit Finite-Volume Model for Overland Flow
journal, September 1998


A comparison of two physics-based numerical models for simulating surface water–groundwater interactions
journal, April 2010


Hyperresolution global land surface modeling: Meeting a grand challenge for monitoring Earth's terrestrial water: OPINION
journal, May 2011

  • Wood, Eric F.; Roundy, Joshua K.; Troy, Tara J.
  • Water Resources Research, Vol. 47, Issue 5
  • DOI: 10.1029/2010WR010090

Support-Operator Finite-Difference Algorithms for General Elliptic Problems
journal, April 1995

  • Shashkov, Mikhail; Steinberg, Stanly
  • Journal of Computational Physics, Vol. 118, Issue 1
  • DOI: 10.1006/jcph.1995.1085

Application of the mixed hybrid finite element approximation in a groundwater flow model: Luxury or necessity?
journal, November 1994

  • Mosé, R.; Siegel, P.; Ackerer, P.
  • Water Resources Research, Vol. 30, Issue 11
  • DOI: 10.1029/94WR01786

Improved numerical solvers for implicit coupling of subsurface and overland flow
journal, December 2014


A general mass-conservative numerical solution for the unsaturated flow equation
journal, July 1990

  • Celia, Michael A.; Bouloutas, Efthimios T.; Zarba, Rebecca L.
  • Water Resources Research, Vol. 26, Issue 7
  • DOI: 10.1029/WR026i007p01483

Hydrologic conditions controlling runoff generation immediately after wildfire: HYDROLOGIC CONDITIONS CONTROLLING RUNOFF AFTER WILDFIRE
journal, March 2012

  • Ebel, Brian A.; Moody, John A.; Martin, Deborah A.
  • Water Resources Research, Vol. 48, Issue 3
  • DOI: 10.1029/2011WR011470

Modeling Coupled Surface-Subsurface Flow Processes: A Review
journal, May 2008


A Mathematical Model for the Coupled Simulation of Surface Runoff and Infiltration
journal, January 2000

  • Gandolfi, C.; Savi, F.
  • Journal of Agricultural Engineering Research, Vol. 75, Issue 1
  • DOI: 10.1006/jaer.1999.0484

Physically based modeling in catchment hydrology at 50: Survey and outlook: SURVEY OF PHYSICALLY BASED MODELING IN CATCHMENT HYDROLOGY
journal, September 2015

  • Paniconi, Claudio; Putti, Mario
  • Water Resources Research, Vol. 51, Issue 9
  • DOI: 10.1002/2015WR017780

HydroGeoSphere: A Fully Integrated, Physically Based Hydrological Model
journal, December 2011


Modeling of Conjunctive Two-Dimensional Surface-Three-Dimensional Subsurface Flows
journal, February 2002


Unsteady one‐dimensional flow over a plane: Partial equilibrium and recession hydrographs
journal, April 1980


Simulating Surface and Subsurface Water Balance Changes Due to Burn Severity
journal, January 2018

  • Atchley, Adam L.; Kinoshita, Alicia M.; Lopez, Sonya R.
  • Vadose Zone Journal, Vol. 17, Issue 1
  • DOI: 10.2136/vzj2018.05.0099

A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow
journal, May 2008


Global threats to human water security and river biodiversity
journal, September 2010

  • Vörösmarty, C. J.; McIntyre, P. B.; Gessner, M. O.
  • Nature, Vol. 467, Issue 7315
  • DOI: 10.1038/nature09440

Nonlinear Krylov acceleration applied to a discrete ordinates formulation of the k-eigenvalue problem
journal, April 2013

  • Calef, Matthew T.; Fichtl, Erin D.; Warsa, James S.
  • Journal of Computational Physics, Vol. 238
  • DOI: 10.1016/j.jcp.2012.12.024

Managing complexity in simulations of land surface and near-surface processes
journal, April 2016


First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach
journal, March 2009

  • Ebel, Brian A.; Mirus, Benjamin B.; Heppner, Christopher S.
  • Hydrological Processes, Vol. 23, Issue 13
  • DOI: 10.1002/hyp.7279

A framework for classifying and comparing distributed hillslope and catchment hydrologic models: DISTRIBUTED MODEL REVIEW
journal, May 2007

  • Kampf, Stephanie K.; Burges, Stephen J.
  • Water Resources Research, Vol. 43, Issue 5
  • DOI: 10.1029/2006WR005370

New preconditioning strategy for Jacobian-free solvers for variably saturated flows with Richards’ equation
journal, August 2016


Hydrologic-Response simulations for the R-5 catchment with a comprehensive physics-based model
journal, April 2001

  • VanderKwaak, Joel E.; Loague, Keith
  • Water Resources Research, Vol. 37, Issue 4
  • DOI: 10.1029/2000WR900272

Integrated surface–groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model
journal, July 2006


A Second-Order Accurate, Finite Volume-Based, Integrated Hydrologic Modeling (FIHM) Framework for Simulation of Surface and Subsurface Flow
journal, November 2009

  • Kumar, Mukesh; Duffy, Christopher J.; Salvage, Karen M.
  • Vadose Zone Journal, Vol. 8, Issue 4
  • DOI: 10.2136/vzj2009.0014

The role of bedrock topography on subsurface storm flow: THE ROLE OF BEDROCK TOPOGRAPHY ON SUBSURFACE STORM FLOW
journal, December 2002

  • Freer, Jim; McDonnell, J. J.; Beven, K. J.
  • Water Resources Research, Vol. 38, Issue 12
  • DOI: 10.1029/2001WR000872

Surface-subsurface model intercomparison: A first set of benchmark results to diagnose integrated hydrology and feedbacks
journal, February 2014

  • Maxwell, Reed M.; Putti, Mario; Meyerhoff, Steven
  • Water Resources Research, Vol. 50, Issue 2
  • DOI: 10.1002/2013WR013725

Design and Application of a Gradient-Weighted Moving Finite Element Code II: in Two Dimensions
journal, May 1998


A fully coupled physically-based spatially-distributed model for evaluating surface/subsurface flow
journal, April 2004


An Investigation of Numerical Grid Effects in Parameter Estimation
journal, November 2006


Design and Application of a Gradient-Weighted Moving Finite Element Code I: in One Dimension
journal, May 1998


A scalable coupled surface–subsurface flow model
journal, August 2015


A Wetland Simulation Module for the MODFLOW Ground Water Model
journal, September 1998


Mixed-hybrid finite element approximations of second-order elliptic boundary-value problems
journal, May 1977

  • Babuska, I.; Oden, J. T.; Lee, J. K.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 11, Issue 2
  • DOI: 10.1016/0045-7825(77)90058-5

Mimetic finite difference method
journal, January 2014

  • Lipnikov, Konstantin; Manzini, Gianmarco; Shashkov, Mikhail
  • Journal of Computational Physics, Vol. 257
  • DOI: 10.1016/j.jcp.2013.07.031

Connections between groundwater flow and transpiration partitioning
journal, July 2016


Hyper-resolution global hydrological modelling: what is next?: “Everywhere and locally relevant”
journal, December 2014

  • Bierkens, Marc F. P.; Bell, Victoria A.; Burek, Peter
  • Hydrological Processes, Vol. 29, Issue 2
  • DOI: 10.1002/hyp.10391

Analysis of coupling errors in a physically-based integrated surface water–groundwater model
journal, December 2012


Overland Flow on an Infiltrating Surface
journal, August 1971


A generalized Richards equation for surface/subsurface flow modelling
journal, March 2009


Newton–Krylov-multigrid solvers for large-scale, highly heterogeneous, variably saturated flow problems
journal, July 2001


The mimetic finite difference method for elliptic and parabolic problems with a staggered discretization of diffusion coefficient
journal, January 2016

  • Lipnikov, Konstantin; Manzini, Gianmarco; Moulton, J. David
  • Journal of Computational Physics, Vol. 305
  • DOI: 10.1016/j.jcp.2015.10.031

Comparative Study of Coupling Approaches for Surface Water and Subsurface Interactions
journal, May 2009


A Closed-form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils1
journal, January 1980


Fill and Spill Hillslope Runoff Representation With a Richards Equation‐Based Model
journal, November 2019

  • Camporese, M.; Paniconi, C.; Putti, M.
  • Water Resources Research, Vol. 55, Issue 11
  • DOI: 10.1029/2019WR025726

Relationships among some locally conservative discretization methods which handle discontinuous coefficients
journal, December 2004