Simulation of Heat Flow in a Synthetic Watershed: The Role of the Unsaturated Zone
Open Access
Abstract
Future climate forecasts suggest atmospheric warming, with expected effects on aquatic systems (e.g., cold-water fisheries). Here we apply a recently published and computationally efficient approach for simulating unsaturated/saturated heat transport with coupled flow (MODFLOW) and transport (MT3D-USGS) models via a synthetic three-dimensional (3D) representation of a temperate watershed. Key aspects needed for realistic representation at the watershed-scale include climate drivers, a layering scheme, consideration of surface-water groundwater interactions, and evaluation of transport parameters influencing heat flux. The unsaturated zone (UZ), which is typically neglected in heat transport simulations, is a primary focus of the analysis. Results from three model versions are compared—one that neglects UZ heat-transport processes and two that simulate heat transport through a (1) moderately-thick UZ and (2) a UZ of approximately double thickness. The watershed heat transport is evaluated in terms of temperature patterns and trends in the UZ, at the water table, below the water table (in the groundwater system), and along a stream network. Major findings are: (1) Climate forcing is the product of infiltration temperatures and infiltration rates; they combine into a single heat inflow forcing function. (2) The UZ acts as a low-pass filter on heat pulses migrating downward, markedly dampening the warming recharge signal. (3) The effect of warming on the watershed is also buffered by the mixing of temperatures at discharge points where shallow and deep flow converge. (4) The lateral extent of the riparian zone, defined as where the water table is near land surface (<1 m), plays an important role in determining the short-term dynamics of the stream baseflow response to heat forcing. Runoff generated from riparian areas is particularly important in periods when rejected infiltration during warm and wet periods generates extra runoff from low-lying areas to surface water.Keywords
Funding Information
- U.S. Geological Survey Land Change Science/Climate Research and Development
This publication has 24 references indexed in Scilit:
- Simulation of climate change effects on streamflow, groundwater, and stream temperature using GSFLOW and SNTEMP in the Black Earth Creek Watershed, Wisconsin2016
- MT3D-USGS version 1: A U.S. Geological Survey release of MT3DMS updated with new and expanded transport capabilities for use with MODFLOW2016
- Simulation of groundwater and surface-water resources of the Santa Rosa Plain watershed, Sonoma County, California2014
- Appraising options to reduce shallow groundwater tables and enhance flow conditions over regional scales in an irrigated alluvial aquifer systemJournal of Hydrology, 2013
- Simulation of climate-change effects on streamflow, lake water budgets, and stream temperature using GSFLOW and SNTEMP, Trout Lake Watershed, Wisconsin2013
- Modeling Variably Saturated Subsurface Solute Transport with MODFLOW‐UZF and MT3DMSGroundwater, 2012
- Comment on “Evaluating Interactions between Groundwater and Vadose Zone Using the HYDRUS‐Based Flow Package for MODFLOW” by Navin Kumar C. Twarakavi, Jirka Šimůnek, and Sophia SeoVadose Zone Journal, 2009
- Importance of Unsaturated Zone Flow for Simulating Recharge in a Humid ClimateGroundwater, 2008
- The Role of Groundwater in Delaying Lake AcidificationWater Resources Research, 1986
- An Experimental Investigation of Runoff Production in Permeable SoilsWater Resources Research, 1970