Nitrogen (N) contamination within agricultural‐karst landscapes and aquifers is widely reported; however, the complex hydrological pathways of karst make N fate difficult to ascertain. We developed a hydrologic and N numerical model for agricultural‐karst, including simulation of soil, epikarst, phreatic, and quick flow pathways as well as biochemical processes such as nitrification, mineralization, and denitrification. We tested the model on four years of nitrate (NO3−) data collected from a phreatic conduit and an overlying surface channel in the Cane Run watershed, Kentucky, USA. Model results indicate that slow to moderate flow pathways (phreatic and epikarst) dominate the N load and account for nearly 90% of downstream NO3−delivery. Further, quick flow pathways dilute NO3−concentrations relative to background aquifer levels. Net denitrification distributed across soil, epikarst, and phreatic water removes approximately 36% of the N inputs to the system at rates comparable to nonkarst systems. Evidence is provided by numerical modeling that NO3−accumulation via evapotranspiration in the soil followed by leaching through the epikarst acts as a control on spring NO3−concentration and loading. Compared to a fluvial‐dominated immature karst system, mature‐karst systems behave as natural detention basins for NO3−, temporarily delaying NO3−delivery to downstream waters and maintaining elevated NO3−concentrations for days to weeks after hydrologic activity ends. This study shows the efficacy of numerical modeling to elucidate complex pathways, processes, and timing of N in karst systems.
Rainfall runoff and leaching are the main driving forces that nitrogen, an important non‐point source (NPS) pollutant, enters streams, lakes, and groundwater. Hydrological and transport processes thus play a pivotal role in NPS nitrogen pollution. Existing hydro‐environmental models for nitrogen pollution often over‐simplify the within‐watershed processes. It is unclear how such simplification affects the pollution assessment regarding the formation and distribution of denitrification hot spots—which is important for the design of land‐based countermeasures. To study this problem, we developed a model, DHSVM‐N, and its variant, DHSVM‐N_alt. DHSVM‐N is developed by integrating nitrogen‐related processes of SWAT into a comprehensive process‐based hydrological model, the Distributed Hydrology Soil and Vegetation Model (DHSVM). DHSVM‐N includes detailed representations of nitrate transport process at a fine spatial resolution with good landscape connectivity to accommodate interactions between hydrological and biogeochemical processes along the flow travel pathways. Because of the lack of spatially distributed observational data for validation, a model‐to‐model comparison study is conducted. Through comparison studies on a representative catchment using SWAT, DHSVM‐N and DHSVM‐N_alt, we quantify the critical roles of hydrological processes and nitrate transport processes in modeling the denitrification process. That is, the capabilities to give reasonable soil moisture estimates and to account for essential processes that take place along flow pathways are keys to simulate denitrification hot spots and their spatial variation. Furthermore, DHSVM‐N results show that terrestrial denitrification from hotspots alone can reach as high as 36% of the annual stream nitrate export of the watershed.
more » « less- NSF-PAR ID:
- 10498294
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 60
- Issue:
- 4
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
Nonpoint source (NPS) pollution is a pressing issue worldwide, especially in the Chesapeake Bay, where sediment, nitrogen (N), and phosphorus (P) are the most critical water quality concerns. Despite significant efforts by federal, state, and local governments, the improvement in water quality has been limited. Investigating the spatial distribution of NPS hotspots can help understand NPS pollutant output and guide control measures. We hypothesize that as land cover changes from natural (e.g., forestland) and agricultural to suburban and ultra-urban, the distribution of NPS pollution source areas becomes increasingly spatially uniform. To test this hypothesis, we analyzed three real watersheds with varying land cover (Greensboro watershed for agriculture, Watts Branch watershed for suburban, and Watershed 263 for ultra-urban) and three synthetic watersheds developed based on the Watts Branch watershed, which ranged from forested and agricultural to ultra-urban but had the same soil, slope, and weather conditions. The Soil and Water Assessment Tool (SWAT) was selected as a phenomenological model for the analysis, and SWAT-CUP was used for model calibration and validation. The hydrologic responses of the three real and synthetic watersheds were simulated over ten years (1993–2002 or 2002–2011), and calibration and validation results indicated that SWAT could properly predict the export of runoff and three target NPS pollution constituents (sediment, total nitrogen, and total phosphorus). The results showed that the distribution of NPS pollutant outputs becomes increasingly uniform as land cover changes from agriculture to ultra-urban across watersheds. This research suggests that the spatial distribution of NPS pollution source areas is a function of the major land cover category of study watersheds, and control strategies should be adapted accordingly. If NPS pollution is distributed unevenly across a watershed, hotspot areas output a disproportionate amount of pollution and require more targeted and intensive control measures. Conversely, if the distribution of NPS pollution is more uniform across a watershed, the control strategies need to be more widespread and encompass a larger area.more » « less
-
Abstract. Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in two simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions.more » « less
-
more » « less
Abstract The compounding effects of anthropogenic legacies for environmental pollution are significant, but not well understood. Here, we show that centennial‐scale legacies of milldams and decadal‐scale legacies of road salt salinization interact in unexpected ways to produce hot spots of nitrogen (N) in riparian zones. Riparian groundwater and stream water concentrations upstream of two mid‐Atlantic (Pennsylvania and Delaware) milldams, 2.4 and 4 m tall, were sampled over a 2 year period. Clay and silt‐rich legacy sediments with low hydraulic conductivity, stagnant and poorly mixed hydrologic conditions, and persistent hypoxia in riparian sediments upstream of milldams produced a unique biogeochemical gradient with nitrate removal via denitrification at the upland riparian edge and ammonium‐N accumulation in near‐stream sediments and groundwaters. Riparian groundwater ammonium‐N concentrations upstream of the milldams ranged from 0.006 to 30.6 mgN L−1while soil‐bound values were 0.11–456 mg kg−1. We attribute the elevated ammonium concentrations to ammonification with suppression of nitrification and/or dissimilatory nitrate reduction to ammonium (DNRA). Sodium inputs to riparian groundwater (25–1,504 mg L−1) from road salts may further enhance DNRA and ammonium production and displace sorbed soil ammonium‐N into groundwaters. This study suggests that legacies of milldams and road salts may undercut the N buffering capacity of riparian zones and need to be considered in riparian buffer assessments, watershed management plans, and dam removal decisions. Given the widespread existence of dams and other barriers and the ubiquitous use of road salt, the potential for this synergistic N pollution is significant.
-
The Baltimore Ecosystem Study (BES) has established a network of long-term permanent biogeochemical study plots. These plots will provide long-term data on vegetation, soil and hydrologic processes in the key ecosystem types within the urban ecosystem. The current network of study plots includes eight forest plots, chosen to represent the range of forest conditions in the area, and four grass plots. These plots are complemented by a network of 200 less intensive study plots located across the Baltimore metropolitan area. Plots are currently instrumented with lysimeters (drainage and tension) to sample soil solution chemistry, time domain reflectometry probes to measure soil moisture, dataloggers to measure and record soil temperature and trace gas flux chambers to measure the flux of carbon dioxide, nitrous oxide and methane from soil to the atmosphere. Measurements of in situ nitrogen mineralization, nitrification and denitrification were made at approximately monthly intervals from Fall 1998 - Fall 2000. Detailed vegetation characterization (all layers) was done in summer 1998. Data from these plots has been published in Groffman et al. (2006, 2009) and Groffman and Pouyat (2009). In November of 1998 four rural, forested plots were established at Oregon Ridge Park in Baltimore County northeast of the Gwynns Falls Watershed. Oregon Ridge Park contains Pond Branch, the forested reference watershed for BES. Two of these four plots are located on the top of a slope; the other two are located midway up the slope. In June of 2010 measurements at the mid-slope sites on Pond Branch were discontinued. Monuments and equipment remain at the two plots. These plots were replaced with two lowland riparian plots; Oregon upper riparian and Oregon lower riparian. Each riparian sites has four 5 cm by 1-2.5 meter depth slotted wells laid perpendicular to the stream, four tension lysimeters at 10 cm depth, five time domain reflectometry probes, and four trace gas flux chambers in the two dominant microtopographic features of the riparian zones - high spots (hummocks) and low spots (hollows). Four urban, forested plots were established in November 1998, two at Leakin Park and two adjacent to Hillsdale Park in west Baltimore City in the Gwynns Falls. One of the plots in Hillsdale Park was abandoned in 2004 due to continued vandalism. In May 1999 two grass, lawn plots were established at McDonogh School in Baltimore County west of the city in the Gwynns Falls. One of these plots is an extremely low intensity management area (mowed once or twice a year) and one is in a low intensity management area (frequent mowing, no fertilizer or herbicide use). In 2009, the McDonogh plots were abandoned due to management changes at the school. Two grass lawn plots were established on the campus of the University of Maryland, Baltimore County (UMBC) in fall 2000. One of these plots is in a medium intensity management area (frequent mowing, moderate applications of fertilizer and herbicides) and one is in a high intensity management area (frequent mowing, high applications of fertilizer and herbicides). Literature Cited Bowden R, Steudler P, Melillo J and Aber J. 1990. Annual nitrous oxide fluxes from temperate forest soils in the northeastern United States. J. Geophys. Res.-Atmos. 95, 13997 14005. Driscoll CT, Fuller RD and Simone DM (1988) Longitudinal variations in trace metal concentrations in a northern forested ecosystem. J. Environ. Qual. 17: 101-107 Goldman, M. B., P. M. Groffman, R. V. Pouyat, M. J. McDonnell, and S. T. A. Pickett. 1995. CH4 uptake and N availability in forest soils along an urban to rural gradient. Soil Biology and Biochemistry 27:281-286. Groffman PM, Holland E, Myrold DD, Robertson GP and Zou X (1999) Denitrification. In: Robertson GP, Bledsoe CS, Coleman DC and Sollins P (Eds) Standard Soil Methods for Long Term Ecological Research. (pp 272-290). Oxford University Press, New York Groffman PM, Pouyat RV, Cadenasso ML, Zipperer WC, Szlavecz K, Yesilonis IC,. Band LE and Brush GS. 2006. Land use context and natural soil controls on plant community composition and soil nitrogen and carbon dynamics in urban and rural forests. Forest Ecology and Management 236:177-192. Groffman, P.M., C.O. Williams, R.V. Pouyat, L.E. Band and I.C. Yesilonis. 2009. Nitrate leaching and nitrous oxide flux in urban forests and grasslands. Journal of Environmental Quality 38:1848-1860. Groffman, P.M. and R.V. Pouyat. 2009. Methane uptake in urban forests and lawns. Environmental Science and Technology 43:5229-5235. DOI: 10.1021/es803720h. Holland EA, Boone R, Greenberg J, Groffman PM and Robertson GP (1999) Measurement of Soil CO2, N2O and CH4 exchange. In: Robertson GP, Bledsoe CS, Coleman DC and Sollins P (Eds) Standard Soil Methods for Long Term Ecological Research. (pp 258-271). Oxford University Press, New York Robertson GP, Wedin D, Groffman PM, Blair JM, Holland EA, Nadelhoffer KJ and. Harris D. 1999. Soil carbon and nitrogen availability: Nitrogen mineralization, nitrification and carbon turnover. In: Standard Soil Methods for Long Term Ecological Research (Robertson GP, Bledsoe CS, Coleman DC and Sollins P (Eds) Standard Soil Methods for Long Term Ecological Research. (pp 258-271). Oxford University Press, New York Savva, Y., K. Szlavecz, R. V. Pouyat, P. M. Groffman, and G. Heisler. 2010. Effects of land use and vegetation cover on soil temperature in an urban ecosystem. Soil Science Society of America Journal 74:469-480."more » « less
