Using MODFLOW to Model Riparian Wetland Shallow Groundwater and Nutrient Dynamics in an Appalachian Watershed

Abstract

Simulating shallow groundwater (SGW) flow dynamics and stream–SGW interactions using numerical modeling tools is necessary to develop a mechanistic understanding of water flow systems and improve confidence in water resource management practices. A three-dimensional (3D) SGW flow model was developed for a riparian wetland in a mixed forest and agricultural catchment in West Virginia (WV), Appalachia, USA, using a Modular 3D Groundwater Model (MODFLOW). The MODFLOW simulation was calibrated in steady (R2 = 0.98, ME = −0.21, and RMSE = 0.77), transient state (R2 = 0.97, ME = −0.41, and RMSE = 1.28) and validated (R2 = 0.97, ME = −0.28, and RMSE = 1.05) using observed SGW levels from thirteen nested piezometers under steady and transient states. An experimental MT3D transport scenario was developed to show the lateral transport of NO₃-N from the aquifer to stream cells. Relatively stable SGW head distribution was observed. In the downstream reach, SGW discharge varied from 948 m3/day to 907 m3/day in 2020, with creek seepage ranging from 802 m3/day to 790 m3/day. Similarly, SGW input to the stream ranged from 891 m3/day to 978 m3/day, while creek seepage ranged from 796 m3/day to 800 m3/day in 2021. In upstream reaches, losing stream conditions were observed in January, June, and September 2020 and January to April 2021, while gaining stream conditions prevailed during other months. Thus, an approximately monthly alternating gaining–losing stream condition was observed in the upstream area. An experimental MT3D transport scenario resulted in an advection–dispersion scenario, showing a cumulative loss of 947 g of NO3-N from SGW to the stream. Denitrification accounted for the cumulative loss of 1406 g of NO3-N from SGW, surpassing 639 g of nitrate from the SGW to the stream during the study period. Additionally, particle tracking using MODPATH indicated a long residence time for SGW nutrients, affirming the efficiency of nitrogen transformation through denitrification. This study is among the first to simulate hydrologic and nutrient interactions in riparian wetlands of a mixed land use catchment in the Appalachian region of the northeastern United States. The results better inform water resource management decisions and modeling efforts in the Appalachian region and similar physiographic regions globally.