Coupling groundwater and land surface processes: Idealized simulations to identify effects of terrain and subsurface heterogeneity on land surface energy fluxes

[1] This work investigates the role of terrain and subsurface heterogeneity on the interactions between groundwater dynamics and land surface energy fluxes using idealized simulations. A three-dimensional variably saturated groundwater code (ParFlow) coupled to a land surface model (Common Land Model) is used to account for both vertical and lateral water and pressure movement. This creates a fully integrated approach, coupling overland and subsurface flow while having an explicit representation of the water table and all land surface processes forced by atmospheric data. Because the water table is explicitly represented in these simulations, regions with stronger interaction between water table depth and the land surface energy balance (known as critical zones) can be identified. This study uses simple terrain and geologic configurations to demonstrate the importance of lateral surface and subsurface flows in determining land surface heat and moisture fluxes. Strong correlations are found between the land surface fluxes and water table depth across all cases, including terrain shape, subsurface heterogeneity, vegetation type, and climatological region. Results show that different land forms and subsurface heterogeneities produce very different water table dynamics and land surface flux responses to atmospheric forcing. Subsurface formation and properties have the greatest effect on the coupling between the water table and surface heat and moisture fluxes. Changes in landform and land surface slope also have an effect on these interactions by influencing the fraction of rainfall contributing to overland flow versus infiltration. This directly affects the extent of the critical zone with highest coupling strength along the hillside. Vegetative land cover, as seen in these simulations, has a large effect on the energy balance at the land surface but a small effect on streamflow and water table dynamics and thus a limited impact on the land surface-subsurface interactions. Although climate forcing has a direct effect on water table dynamics and feedbacks to the land surface, in this study it does not overcome that of subsurface heterogeneity and terrain.

[1]  Zong-Liang Yang,et al.  Preliminary study of spin‐up processes in land surface models with the first stage data of Project for Intercomparison of Land Surface Parameterization Schemes Phase 1(a) , 1995 .

[2]  J. Tóth.  A Theoretical Analysis of Groundwater Flow in Small Drainage Basins , 1963 .

[3]  R. Freeze,et al.  Blueprint for a physically-based, digitally-simulated hydrologic response model , 1969 .

[4]  Clemens Simmer,et al.  The Influence of Hydrologic Modeling on the Predicted Local Weather: Two-Way Coupling of a Mesoscale Weather Prediction Model and a Land Surface Hydrologic Model , 2002 .

[5]  Elfatih A. B. Eltahir,et al.  Representation of Water Table Dynamics in a Land Surface Scheme. Part I: Model Development , 2005 .

[6]  J. Vanderkwaak.  Numerical simulation of flow and chemical transport in integrated surface-subsurface hydrologic systems , 1999 .

[7]  K. Loague,et al.  Physics‐based hydrologic response simulation: platinum bridge, 1958 Edsel, or useful tool , 2004 .

[8]  R. Maxwell,et al.  Interdependence of groundwater dynamics and land-energy feedbacks under climate change , 2008 .

[9]  Zhenghui Xie,et al.  A new parameterization for surface and groundwater interactions and its impact on water budgets with the variable infiltration capacity (VIC) land surface model , 2003 .

[10]  Benjamin B. Mirus,et al.  First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach , 2009 .

[11]  D. Entekhabi,et al.  Groundwater-surface water interaction and the climatic spatial patterns of hillslope hydrological response , 1999 .

[12]  Aldo Fiori,et al.  Numerical analyses of subsurface flow in a steep hillslope under rainfall: The role of the spatial heterogeneity of the formation hydraulic properties , 2007 .

[13]  Dara Entekhabi,et al.  Hillslope and Climatic Controls on Hydrologic Fluxes , 1995 .

[14]  R. Maxwell,et al.  The groundwater land-surface atmosphere connection: Soil moisture effects on the atmospheric boundary layer in fully-coupled simulations , 2007 .

[15]  A. Pitman,et al.  Simulating the diurnal temperature range: Results from Phase I(a) of the Project for Intercomparison of Landsurface Parameterisation Schemes (PILPS ) , 1995 .

[16]  Xi Chen,et al.  Groundwater influences on soil moisture and surface evaporation , 2004 .

[17]  Ying Fan,et al.  Incorporating water table dynamics in climate modeling: 1. Water table observations and equilibrium water table simulations , 2007 .

[18]  G. Salvucci,et al.  Equilibrium analysis of groundwater–vadose zone interactions and the resulting spatial distribution of hydrologic fluxes across a Canadian Prairie , 1999 .

[19]  S. Ashby,et al.  A parallel multigrid preconditioned conjugate gradient algorithm for groundwater flow simulations , 1996 .

[20]  R. Maxwell,et al.  Capturing the influence of groundwater dynamics on land surface processes using an integrated, distributed watershed model , 2008 .

[21]  R. Maxwell,et al.  Integrated surface-groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model , 2006 .

[22]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[23]  Reed M. Maxwell,et al.  Development of a Coupled Land Surface and Groundwater Model , 2005 .

[24]  Fotini Katopodes Chow,et al.  Effects of soil moisture heterogeneity on boundary layer flow with coupled groundwater, land-surface, and mesoscale atmospheric modeling , 2006 .

[25]  R. Maxwell.  Coupled Surface–Subsurface Modeling across a Range of Temporal and Spatial Scales , 2009 .

[26]  R. Allan Freeze,et al.  Theoretical analysis of regional groundwater flow: 1. Analytical and numerical solutions to the mathematical model , 1966 .

[27]  L. A. Richards.  Capillary conduction of liquids through porous mediums , 1931 .

[28]  Zong-Liang Yang,et al.  Improving land‐surface model hydrology: Is an explicit aquifer model better than a deeper soil profile? , 2007 .

[29]  Peter Bayer,et al.  The Influence of Rain Sensible Heat and Subsurface Energy Transport on the Energy Balance at the Land Surface , 2009 .

[30]  Reed M. Maxwell,et al.  Quantifying the effects of three-dimensional subsurface heterogeneity on Hortonian runoff processes using a coupled numerical, stochastic approach , 2008 .

[31]  Jim E. Jones,et al.  Approved for Public Release; Further Dissemination Unlimited Newton-krylov-multigrid Solvers for Large-scale, Highly Heterogeneous, Variably Saturated Flow Problems , 2022 .

[32]  R. Allan Freeze,et al.  Role of subsurface flow in generating surface runoff: 1. Base flow contributions to channel flow , 1972 .

[33]  Marcel G. Schaap,et al.  Database-related accuracy and uncertainty of pedotransfer functions , 1998 .

[34]  Benjamin B. Mirus,et al.  Physics‐based hydrologic‐response simulation: foundation for hydroecology and hydrogeomorphology , 2006 .

[35]  S. Manabe,et al.  SIMULATED CLIMATOLOGY OF A GENERAL CIRCULATION MODEL WITH A HYDROLOGIC CYCLE , 1965 .

[36]  P. Huyakorn,et al.  A fully coupled physically-based spatially-distributed model for evaluating surface/subsurface flow , 2004 .

[37]  Clint Dawson,et al.  A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow , 2008 .

[38]  Pedro Viterbo,et al.  The land surface‐atmosphere interaction: A review based on observational and global modeling perspectives , 1996 .