Subsurface tile drainage (TD) is a dominant agriculture water management practice in the United States (US) to enhance crop production in poorly drained soils. Assessments of field‐level or watershed‐level (<50 km2) hydrologic impacts of TD are becoming common; however, a major gap exists in our understanding of regional (>105 km2) impacts of TD on hydrology. The National Water Model (NWM) is a distributed 1‐km resolution hydrological model designed to provide accurate streamflow forecasts at 2.7 million reaches across the US. The current NWM lacks TD representation which adds considerable uncertainty to streamflow forecasts in heavily tile‐drained areas. In this study, we quantify the performance of the NWM with a newly incorporated tile‐drainage scheme over the heavily tile‐drained Midwestern US. Employing a TD scheme enhanced the uncalibrated NWM performance by about 20–50% of the fully calibrated NWM (
Several studies have focused on the importance of river bathymetry (channel geometry) in hydrodynamic routing along individual reaches. However, its effect on other watershed processes such as infiltration and surface water (SW)‐groundwater (GW) interactions has not been explored across large river networks. Surface and sbsurface processes are interdependent, therefore, errors due to inaccurate representation of one watershed process can cascade across other hydraulic or hydrologic processes. This study hypothesizes that accurate bathymetric representation is not only essential for simulating channel hydrodynamics but also affects subsurface processes by impacting SW‐GW interactions. Moreover, quantifying the effect of bathymetry on surface and subsurface hydrological processes across a river network can facilitate an improved understanding of how bathymetric characteristics affect these processes across large spatial domains. The study tests this hypothesis by developing physically based distributed models capable of bidirectional coupling (SW‐GW) with four configurations with progressively reduced levels of bathymetric representation. A comparison of hydrologic and hydrodynamic outputs shows that changes in channel geometry across the four configurations has a considerable effect on infiltration, lateral seepage, and location of water table across the entire river network. For example, when using bathymetry with inaccurate channel conveyance capacity but accurate channel depth, peak lateral seepage rate exhibited 58% error. The results from this study provide insights into the level of bathymetric detail required for accurately simulating flooding‐related physical processes while also highlighting potential issues with ignoring bathymetry across lower order streams such as spurious backwater flow, inaccurate water table elevations, and incorrect inundation extents.
more » « less- Award ID(s):
- 1706612
- NSF-PAR ID:
- 10382366
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 58
- Issue:
- 11
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract Calib ). The calibrated NWM with tile drainage (CalibTD ) showed enhanced accuracy with higher event hit rates and lower false alarm rates thanCalib .CalibTD showed better performance in high‐flow estimations as TD increased streamflow peaks (14%), volume (2.3%), and baseflow (11%). Regional water balance analysis indicated that TD significantly reduced surface runoff (−7% to −29%), groundwater recharge (−43% to −50%), evapotranspiration (−7% to −13%), and soil moisture content (−2% to −3%). However, TD significantly increased soil profile lateral flow (27.7%) along with infiltration and soil water storage potential. Overall, our findings highlight the importance of incorporating the TD process into the operational configuration of the NWM. -
more » « less
Abstract A three‐dimensional unstructured‐grid hydrodynamic and water quality model (Semi‐implicit Cross‐scale Hydroscience Intergrated System Model‐Integrated Compartment Model) is applied successfully for Chesapeake Bay. The model is validated with observations of salinity, chlorophyll‐a, dissolved oxygen, nutrients, and phytoplankton productions from the year 1991 to 1995 for the mainstem and some major tributaries, based on multiple model skill scores. Model experiments are conducted to test the importance of having (1) an accurate representation of bathymetry to correctly predict hypoxia and other processes and (2) a high‐resolution model grid for tributaries to correctly simulate water quality variables. Comparison with the model experiment results with bathymetry smoothing indicates that bathymetry smoothing, as commonly used for many systems, changes the stratification and lateral circulation pattern, resulting in more salt intrusion into shallow water regions, and an increase in the freshwater age. Consequently, a model with bathymetry smoothing can lead to an unrealistic prediction of the distribution of hypoxia and phytoplankton production. Local grid refinement shows significant improvement of model simulations on local stratification and water quality variables. Overall, the use of high‐resolution unstructured grid model leads to a faithful representation of the complex geometry, and thus a seamless cross‐scale capability for simulating water quality processes in the Bay including tributaries and tidal creeks.
-
more » « less
Abstract Intermittent streams currently constitute >50% of the global river network, and the number of intermittent streams is expected to increase due to changes in land use and climate. Surface flows are known to expand and contract within the headwater channel network due to changes in the water table driven by climate, often changing seasonally. However, the underlying causes of disconnections and reconnections throughout the stream network remain poorly understood and may reflect subsurface flow capacity. We assess how 3D subsurface flowpaths control local surface flows at Gibson Jack Creek in the Rocky Mountains, Idaho, USA. Water table dynamics, hydraulic gradients, and hyporheic exchange were monitored along a 200‐m section of the stream throughout the seasonal recession in WY2018. Shallow lateral hillslope‐riparian‐stream connectivity was more frequent in transects spanning perennially flowing stream reaches than intermittent reaches. During low‐flow periods, larger losing vertical hydraulic gradients were observed in paired piezometers in intermittent reaches than in adjacent perennial reaches. Contrary to dominant conceptual models, longitudinal measurements of hydrologic exchange in both intermittent and perennial reaches were seasonally variable except for one perennial reach that showed consistent significant gains. Observed drying dynamics, as well as subsurface pathways, were highly variable even over short distances (30 m). Flow probability and subsurface flow capacity at upstream locations can be assessed with an outlet hydrograph and upstream flow measurements. Accurate characterization of subsurface storage, discharge, and connection is critical to understanding the drivers of drying cycles in intermittent streams and their likely responses to future change.
-
This study attempts to integrate a Surface Water (SW) model Soil and Water Assessment Tool (SWAT) with an existing steady-state, single layer, unconfined heterogeneous aquifer Analytic Element Method (AEM) based Ground Water (GW) model, named Bluebird AEM engine, for a comprehensive assessment of SW and GW resources and its management. The main reason for integrating SWAT with the GW model is that the SWAT model does not simulate the distribution and dynamics of GW levels and recharge rates. To overcome this issue, often the SWAT model is coupled with the numerical GW model (either using MODFLOW or FEFLOW), wherein the spatial and temporal patterns of the interactions are better captured and assessed. However, the major drawback in integrating the two models (SWAT with—MODFLOW/FEM) is its conversion from Hydrological Response Unit’s (HRU)/sub-basins to grid/elements. To couple them, a spatial translation system is necessary to move the inputs and outputs back and forth between the two models due to the difference in discretization. Hence, for effective coupling of SW and GW models, it may be desirable to have both models with a similar spatial discretization and reduce the need for rigorous numerical techniques for solving the PDEs. The objective of this paper is to test the proof of concept of integrating a distributed hydrologic model with an AEM model at the same spatial units, primarily focused on surface water and groundwater interaction with a shallow unconfined aquifer. Analytic Element Method (AEM) based GW models seem to be ideal for coupling with SWAT due to their innate character to consider the HRU, sub-basin, River, and lake boundaries as individual analytic elements directly without the need for any further discretization or modeling units. This study explores the spatio-temporal patterns of groundwater (GW) discharge rates to a river system in a moist-sub humid region with SWAT-AEM applied to the San Jacinto River basin (SJRB) in Texas. The SW-GW interactions are explored throughout the watershed from 2000–2017 using the integrated SWAT-AEM model, which is tested against stream flow and GW levels. The integrated SWAT-AEM model results show good improvement in predicting the stream flow (R2 = 0.65–0.80) and GW levels as compared to the standalone SWAT model. Further, the integrated model predicted the low flows better compared to the standalone SWAT model, thus accounting for the SW-GW interactions. Almost 80% of the stream network experiences an increase in groundwater discharge rate between 2000 and 2017 with an annual average GW discharge rate of 1853 Mm3/year. The result from the study seems promising for potential applications of SWAT-AEM coupling in regions with considerable SW-GW interactions.more » « less
-
more » « less
Abstract High‐resolution topography reveals that floodplains along meandering rivers in Indiana commonly contain intermittently flowing channel networks. We investigated how the presence of floodplain channels affects lateral surface‐water connectivity between a river and floodplain (specifically exchange flux and timescales of transport) as a function of flow stage in a low‐gradient river‐floodplain system. We constructed a two‐dimensional, surface‐water hydrodynamic model using Hydrologic Engineering Center's River Analysis System (HEC‐RAS) 2D along 32 km of floodplain (56 km along the river) of the East Fork White River near Seymour, Indiana, USA, using lidar elevation data and surveyed river bathymetry. The model was calibrated using land‐cover specific roughness to elevation‐discharge data from a U.S. Geological Survey gage and validated against high‐water marks, an aerial photo showing the spatial extent of floodplain inundation, and measured flow velocities. Using the model results, we analyzed the flow in the river, spatial patterns of inundation, flow pathways, river‐floodplain exchange, and water residence time on the floodplain. Our results highlight that bankfull flow is an oversimplified concept for explaining river‐floodplain connectivity because some stream banks are overtopped and major low‐lying floodplain channels are inundated roughly 19 days per year. As flow increased, inundation of floodplain channels at higher elevations dissected the floodplain, until the floodplain channels became fully inundated. Additionally, we found that river‐floodplain exchange was driven by bank height or channel orientation depending on flow conditions. We propose a conceptual model of river‐floodplain connectivity dynamics and developed metrics to analyze quantitatively complex river‐floodplain systems.
