Groundwater contributions of flow and nitrogen in a headwater agricultural watershed

Nonpoint sources of nitrogen (N) and other nutrients are a major source of water pollution within the Chesapeake Bay watershed and other basins around the world. Human activities associated with agricultural practices can account for a large percentage of N loadings delivered to streams and rivers. This work aims to improve understanding of N transport from groundwater to surface waters, quantifying the principal hydrological processes driving water and N fluxes into and out of a headwater agricultural stream reach. The study site is a 175-m stream reach in a heavily cultivated 40-ha watershed in east-central Pennsylvania. This subwatershed is underlain by fractured shale bedrock, and receives most of its baseflow from groundwater, either by diffuse matrix discharge through the streambed or by localized discharge through riparian seeps. Samples of stream, seep, and shallow groundwater were collected approximately monthly under steady hydrologic conditions in 2017. Calculated matrix flow from hydraulic head and conductivity measurements paired with differential stream gauging was used to solve for the riparian seep flux using a mass balance approach. Riparian seep fluxes ranged from 45 to 217 m³/d, transporting 0.6–4.2 kg N d⁻¹ of nitrate-N from the fractured bedrock aquifer to the stream. Hydrochemical data suggest that the stream is mainly disconnected from the underlying aquifer and that seeps supply essentially all water and N to the system. Seeps are likely sourced with N in nearby agricultural fields and accelerated through the system with shorter residence times than shallow groundwater. Water isotope data reinforced this notion. This study underscores the importance of agriculture as a source of N to ground and surface waters. Identifying source areas that are causing groundwater enrichment of N and seep areas where N discharges to streams is beneficial for developing N pollution mitigation strategies and implementing management practices that aim to reduce nutrient loads to the Chesapeake Bay.     

基于改进SWAT模型的灌溉—施肥—耕作对乌梁素海流域营养物负荷及作物产量的影响

不合理的灌溉、施肥和耕作是导致乌梁素海流域农业面源污染的主要根源,乌梁素海作为我国北方地区重要的生态安全屏障,多年来面临着湖泊水环境污染、水生态退化等问题,科学开展湖泊水环境综合治理首先要解决流域内农业面源污染问题.研究通过修改土壤水平衡、溶质平衡、地下水平衡和作物生长等模块对SWAT 2012原始版本进行改进,并采用改进的SWAT模型构建了乌梁素海流域分布式水文模型,利用实测径流、硝态氮与总磷排放量、地下水埋深以及作物产量校正和验证模型.基于现状情景,以玉米、葵花和小麦3种主要作物为研究对象,设置了削减灌水量、施肥量及调整耕作方式3种农田管理情景.基于改进SWAT模型不同情景的模拟结果,计算分析各管理情景下的硝态氮与总磷负荷及对各作物产量的影响.结果表明,改进SWAT模型具有良好的模拟效果.不同作物削减5%夏灌水量增产最多达8.41%～10.32%,削减10%秋浇水量均明显减少硝态氮和总磷负荷.不同作物营养物负荷均随着氮磷施肥削减比例的增大呈现逐渐降低的趋势,但下降曲线逐渐趋于平缓;各作物产量随氮磷施肥削减比例的增加呈先增加后减少的趋势,其中玉米、小麦氮磷施肥削减比例达20%时产量开...     

Validation and Application of the Second Generation Three Dimensional Hydrodynamic Model of Chesapeake Bay

The validation and subsequent application of the current three-dimensional numerical hydrodynamic model of Chesapeake Bay is presented. The numerical model solves conservation equations for water mass, momentum, salinity, and heat on a boundary-fitted grid in the horizontal plane and a Cartesian z-grid in the vertical. A generalized ADI finite difference scheme is employed in conjunction with mode splitting technique, solving external and the internal modes. The 10-year boundary conditions including tide, slinity, temperature, wind, heat exchange coefficient, river and non-point source flows were constructed. Model validation was accomplished by demonstrating the model's ability to reproduce observed data over time scales ranging from tidal to seasonal periods. The major parameters compared include tidal elevation, intra-tidal and residual velocities, salinity, temperature, stratification, and flux calculated through the Bay mouth.After validation, the model was applied to simulate bay hydrodynamics for the 10 years of 1985–94. These results were used to drive the three-dimensional water quality model of Chesapeake Bay, which is discussed in a companion paper.     

美国罗德河流域地下水中铝形成动力学

In the Rhode River estuarine/watershed area, increased acidity of precipitation from atmospheric acidic deposition has deleterious effects on the freshwater ecosystems. One of the characteristics of an acidic watershed is the mobilization of aluminum from soils to aquatic environment. Increased concentrations of aluminum in surface waters are toxic to living organisms. Detection of long-term changes of acidity and elevated Al concentrations in surface and soil waters is a central issue. In this present paper, the dynamics ofAl speciation in the shallow ground waters from Rhode River watershed 110 (la) were investigated. This research provides a unique, regional analysis of important controls on the transport of Al speciation through the rural watershed of the Chesapeake Bay.     

Predicting dissolved oxygen in the Chesapeake Bay: applications and implications

Eutrophic depletion of dissolved oxygen (DO) and its consequences for ecosystem dynamics have been a central theme of research, assessment and management policies for several decades in the Chesapeake Bay. Ongoing forecast efforts predict the extent of the summer hypoxic/anoxic area due to nutrient loads from the watershed. However, these models neither predict DO levels nor address the intricate interactions among various ecological processes. The prediction of spatially explicit DO levels in the Chesapeake Bay can eventually lead to a reliable depiction of the comprehensive ecological structure and functioning, and can also allow the quantification of the role of nutrient reduction strategies in water quality management. In this paper, we describe a three dimensional empirical model to predict DO levels in the Chesapeake Bay as a function of water temperature, salinity and dissolved nutrient concentrations (TDN and TDP). The residual analysis shows that predicted DO values compare well with observations. Nash–Sutcliffe efficiency (NSE) and root mean square error-observations standard deviation ratio (RSR) are used to evaluate the performance of the empirical model; the scores demonstrate the usability of model predictions (NSE, surface layer = 0.82–0.86; middle layer = 0.65–0.82; bottom layer = 0.70–0.82; RSR surface layer = 0.37–0.44; middle layer = 0.43–0.58 and bottom layer = 0.43–0.54). The predicted DO values and other physical outputs from downscaling of regional weather and climate predictions, or forecasts from hydrodynamic models, can be used to forecast various ecological components. Such forecasts would be useful for both recreational and commercial users of the Chesapeake Bay.     

Modelling atmospheric and hydrologic processes for assessment of meadow restoration impact on flow and sediment in a sparsely gauged California watershed

The restoration of meadowland using the pond and plug technique of gully elimination was performed in a 9‐mile segment along Last Chance Creek, Feather River Basin, California, in order to rehabilitate floodplain functions such as mitigating floods, retaining groundwater, and reducing sediment yield associated with bank erosion and to significantly alter the hydrologic regime. However, because the atmospheric and hydrological conditions have evolved over the restoration period, it was difficult to obtain a comprehensible evaluation of the impact of restoration activities by means of field measurements. In this paper, a new use of physically based models for environmental assessment is described. The atmospheric conditions over the sparsely gauged Last Chance Creek watershed (which does not have any precipitation or weather stations) during the combined historical critical dry and wet period (1982–1993) were reconstructed over the whole watershed using the atmospheric fifth‐generation mesoscale model driven with the US National Center for Atmospheric Research and US National Center for Environmental Prediction reanalysis data. Using the downscaled atmospheric data as its input, the watershed environmental hydrology (WEHY) model was applied to this watershed. All physical parameters of the WEHY model were derived from the existing geographic information system and satellite‐driven data sets. By comparing the prerestoration and postrestoration simulation results under the identical atmospheric conditions, a more complete environmental assessment of the restoration project was made. Model results indicate that the flood peak may be reduced by 10–20% during the wet year and the baseflow may be enhanced by 10–20% during the following dry seasons (summer to fall) in the postrestoration condition. The model results also showed that the hydrologic impact of the land management associated with the restoration mitigates bank erosion and sediment discharge during winter storm events. Copyright © 2013 John Wiley & Sons, Ltd.     

Impact of climate change on diffuse pollutant fluxes at the watershed scale

This study aims to assess watershed‐scale impacts of changing climate on sediment, phosphorus, nitrogen and pesticide (atrazine) fluxes over the 21st century at the watershed scale. In particular, changes in dissolved and particulate forms of water quality constituents in response to climate change are investigated. The hydrologic model Soil and Water Assessment Tool was calibrated and evaluated in a primarily agricultural watershed in the Midwestern United States to simulate hydrologic and water quality processes on a daily basis over the 2015–2099 time horizon. The model was then driven with 112 distinct statistically downscaled climate projections representing Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios (IPCC SRES) low, moderate and high greenhouse gas emission scenarios. Projected hydrologic and water quality responses were categorized according to the three IPCC SRES emission scenarios for summarizing and synthesizing results over early‐century (2015–2034), mid‐century (2045–2064) and late‐century (2080–2099) assessment. Results revealed clear warming trends in the study area, whereas small increases in precipitation were predicted. Streamflow, sediment and total nutrient loads did not differ noticeably between assessment periods. However, the proportion of dissolved to total nutrients increased significantly from early‐century to late‐century periods. With the exception of total atrazine in the mid‐century period, predicted pollutant loads for a given assessment period did not differ between emission pathways for a given assessment period. Changes in pollutant fluxes showed pronounced monthly variability. The projected increase in readily available forms of nutrients has important implications for the ecological health of water systems and management of drinking water supplies. Copyright © 2013 John Wiley & Sons, Ltd.     

The impact of submersed aquatic vegetation on the development of river mouth bars

This work is inspired by the sudden resurgence of the submersed aquatic vegetation (SAV) bed in the Chesapeake Bay (USA). Because the SAV bed occurs at the mouth of the Bay's main tributary (Susquehanna River), it plays a significant role in modulating sediment and nutrient inputs from the Susquehanna to the Bay. Previous model studies on the impact of submersed aquatic vegetation on the development of river mouth bars lacked a complete mechanistic understanding. This study takes advantage of new advances in 3D computational models that include explicit physical-sedimentological feedbacks to obtain this understanding. Specifically, we used Delft3D, a state-of-the-art hydrodynamic model that provides fine-scale computations of three-dimensional flow velocity and bed shear stress, which can be linked to sediment deposition and erosion. Vegetation is modeled using a parameterization of hydraulic roughness that depends on vegetation height, stem density, diameter, and drag coefficient. We evaluate the hydrodynamics, bed shear stresses, and sediment dynamics for different vegetation scenarios under conditions of low and high river discharge. Model runs vary the vegetation height, density, river discharge, and suspended-sediment concentration. Numerical results from the idealized model show that dense SAV on river mouth bars substantially diverts river discharge into adjacent channels and promotes sediment deposition at ridge margins, as well as upstream bar migration. Increasing vegetation height and density forms sandier bars closer to the river mouth and alteration of the bar shape. Thus, this study highlights the important role of SAV in shaping estuarine geomorphology, which is especially relevant for coastal management. © 2019 John Wiley & Sons, Ltd.     

Phytoplankton Kinetics in the Chesapeake Bay Eutrophication Model

The CE-Qual-ICM model computes phytoplankton biomass and production as a function of temperature, light, and nutrients. Biomass is computed as carbon while inorganic nitrogen, phosphorus, and silica are considered as nutrients. Model formulations for production, metabolism, predation, nutrient limitation, and light limitation are detailed. Methods of parameter determination and parameter values are presented. Results of model application to a ten-year period in Chesapeake Bay indicate the model provides reasonable representations of observed biomass, nutrient concentrations, and limiting factors. Computed primary production agrees with observed under light-limited conditions. Under strongly nutrient-limited conditions, computed product is less than observed. The production characteristics of the model are similar to behavior reported for several similar models. Process omitted from the model that may account for production shortfalls include variable algal stoichiometry, use of urea as nutrient, and vertical migration by phytoplankton.     

Application of a watershed runoff model to North-east Pond River,Newfoundland: to study water balance and hydrological characteristics owing to atmospheric change

The hydrological sensitivities to long-term climate change of a watershed in Eastern Canada were analysed using a deterministic watershed runoff model developed to simulate watershed acidification. This model was modified to study atmospheric change effects in the watershed. Water balance modelling techniques, modified for assessing climate effects, were developed and tested for a watershed using atmospheric change scenarios from both state of the art general circulation models and a series of hypothetical scenarios. The model computed daily surface, inter- and groundwater flows from the watershed. The moisture, infiltration and recharge rate are also computed in the soil reservoirs. The thirty years of simulated data can be used to evaluate the effects of climatic change on soil moisture, recharge rate and surface and subsurface flow systems. The interaction between surface and subsurface water is discussed in relation to climate change. These hydrological results raise the possibility of major environmental and socioeconomic difficulties and have significant implications for future water resource planning and management. © 1997 John Wiley & Sons, Ltd.     

Impacts of environmental stressors on the water resources of intensively managed hydrologic systems

Watersheds are complex systems due to their surface and subsurface spatially connected water fluxes and biochemical processes that shape Earth's critical zone. In intensively managed landscapes, the implementation of watershed management practices (WMPs) regulate their short‐term responses, whereas climate variability controls the long‐term processes. Understanding their responses to anthropogenic and natural stressors requires a holistic approach that takes into account their multiscale spatio‐temporal linkages. The objective of this study was to simulate the impacts of spatially and temporally varying WMPs and projected climate changes on the surface and groundwater resources in the Upper Sangamon River Basin (USRB), a watershed in central Illinois greatly impacted by agricultural and industrial operations. The physically based hydrologic model MIKE‐SHE was used to simulate the hydrologic responses of the basin to different WMPs and climatic conditions. The simulation of a WMP was varied spatially across the basin to determine the spectrum of responses and critical conditions. In general, the wetlands and forested riparian buffer scenarios were found to cause a reduction in the average streamflow, whereas crop rotation had varied responses depending on the location of implementation and the climate condition assumed. Reductions of up to 30% in the average streamflow were found for the forested riparian buffer under the ESM 2M climate projections, whereas an increase of up to 13% with the crop rotation schemes under CM3 climate was predicted. The model results showed that the installation of tile drains across the USRB increased the water table depth (from ground level) by up to 56%, making crop production possible. Groundwater level in USRB appeared to be more sensitive to future climatic conditions than to WMP implementation. The impacts of WMPs are determined to depend on the climate conditions under which they are applied. Investigating individual and combined stressors' effects over the critical zone at a watershed scale can lead to a more comprehensive analysis of the risk and trade‐offs in every managerial decision that will enable an efficient use of resources.     

Improved water quality in response to pollution control measures at Masan Bay, Korea

The total pollution load management system (TPLMS) was first applied in 2007 to the highly developed Masan Bay watershed, Korea. To evaluate the effect of TPLMS on water quality improvement, we analyzed the water qualities in rivers and bay during 2005-2010, targeting chemical oxygen demand (COD), suspended sediment (SS), total nitrogen (TN), and total phosphorus (TP) loads. Land-based pollutant loading all decreased during this period, with a significant reduction in COD and SS loads (p<0.01). The COD reduction in seawater, following the TPLMS implementation, was also significant (p<0.01). Time-lagged responses in COD and Chl-a supported an estimated seawater residence time of ~1 month. Land-based nutrient loads were also significantly reduced for TN (p<0.01) and TP (p<0.05), however, significant reductions were not observed in the bay, indicating potential alternative nutrient inputs from non-point sources into the bay system.     

Assessment of Impact of Storm on Point Source Pollutant Transport in Estuary by Dissolved Tracer Modeling

A continuously discharged dissolved conservative tracer was simulated with the Chesapeake Bay Estuary Model Package to study pollutant transport in the estuary in response to point source loads and the impact of the November, 1985 storm. A visualization technology is applied to show 3-dimensional concentration variations in a continuous daily time sequence. The differential responses of daily net transport during storms versus inter-storm periods can be observed from an MPEG movie. It may take 2–3 months for a tracer to travel from the fall-line to the mouth of a river during relatively dry seasons, only 2 weeks in some medium storms, and less than 5 days in a big storm. Plots of daily concentrations from eleven selected locations in the estuary provide quantitative information on the response of tracer concentration to flows. The magnitude and time of tracer peaks related with different weather events in these locations reflect the combined effects of flows from various directions to these locations. The lower tributaries (which are closer to the Bay mouth) are affected more than the upper tributaries by a source discharged at a mid-tributary. A storm can transport materials more effectively to the Bay and affect adjacent tributaries more severely.     

Multivariate adaptive regression splines for estimating riverine constituent concentrations

Regression-based methods are commonly used for riverine constituent concentration/flux estimation, which is essential for guiding water quality protection practices and environmental decision making. This paper developed a multivariate adaptive regression splines model for estimating riverine constituent concentrations (MARS-EC). The process, interpretability and flexibility of the MARS-EC modelling approach, was demonstrated for total nitrogen in the Patuxent River, a major river input to Chesapeake Bay. Model accuracy and uncertainty of the MARS-EC approach was further analysed using nitrate plus nitrite datasets from eight tributary rivers to Chesapeake Bay. Results showed that the MARS-EC approach integrated the advantages of both parametric and nonparametric regression methods, and model accuracy was demonstrated to be superior to the traditionally used ESTIMATOR model. MARS-EC is flexible and allows consideration of auxiliary variables; the variables and interactions can be selected automatically. MARS-EC does not constrain concentration-predictor curves to be constant but rather is able to identify shifts in these curves from mathematical expressions and visual graphics. The MARS-EC approach provides an effective and complementary tool along with existing approaches for estimating riverine constituent concentrations.     

Seasonal and annual loads of hydrophobic organic contaminants from the Susquehanna River basin to the Chesapeake Bay

Water from the Susquehanna River was collected and analyzed for polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyl (PCB) congeners to estimate seasonal and annual riverine loads to the Chesapeake Bay. Temporal variations in the chemical loads resulted from the large changes in the water flow rates and in the particle-associated contaminant concentrations. Concentrations of PCBs and PAHs in river particles (ng/g) were twice as great as those in the northern Chesapeake Bay, indicating that the Susquehanna River is an important source of these contaminants to the bay. The river carries a majority of its hydrophobic organic contaminants (HOCs) in the particulate phase. During periods of high flow, large amounts of suspended particles in the river result in elevated HOC levels and increased loadings of these contaminants to the bay. From 1997 to 1998, 60% of the total annual HOC loading occurred in the early spring coincident with high river flows. The total PCB and PAH annual loadings from the Susquehanna River to the Chesapeake Bay were 76 and 3160 kg/year, respectively and 75% of the loaded organic contaminants were in the particulate phase. Principal component analysis of PAH and PCB congener patterns in the particles reveals that the river suspended particles were dominated by autochthonous production in the summer and by resuspended sediment and watershed erosion during the winter and early spring.     

Assessing impacts of agricultural water interventions in the Kothapally watershed,Southern India

The paper describes a hydrological model for agricultural water intervention in a community watershed at Kothapally in India, developed through integrated management and a consortium approach. The impacts of various soil and water management interventions in the watershed are compared to no‐intervention during a 30‐year simulation period by application of the calibrated and validated ARCSWAT 2005 (Version 2.1.4a) modelling tool. Kothapally receives, on average, 800 mm rainfall in the monsoon period. 72% of total rainfall is converted as evaporation and transpiration (ET), 20% is stored by groundwater aquifer, and 8% exported as outflow from the watershed boundary in current water interventions. ET, groundwater recharge and outflow under no‐intervention conditions are found to be 64, 9, and 19%, respectively. Check dams helped in storing water for groundwater recharge, which can be used for irrigation, as well minimising soil loss. In situ water management practices improved the infiltration capacity and water holding capacity of the soil, which resulted in increased water availability by 10–30% and better crop yields compared to no‐intervention. Water outflows from the developed watershed were more than halved compared to no‐intervention, indicating potentially large negative downstream impacts if these systems were to be implemented on a larger scale. On the other hand, in the watershed development program, sediment loads to the streams were less than one‐tenth. It can be concluded that the hydrological impacts of large‐scale implementation of agricultural water interventions are significant. They result in improved rain‐fed agriculture and improved productivity and livelihood of farmers in upland areas while also addressing the issues of poverty, equity, and gender in watersheds. There is a need for case‐specific studies of such hydrological impacts along with other impacts in terms of equity, gender, sustainability, and development at the mesoscale. Copyright © 2011 John Wiley & Sons, Ltd.     

The role of land surface versus drainage network characteristics in controlling water quality and quantity in a small urban watershed

The processes that control run‐off quantity and quality in urban watersheds are complex and not well understood. Although impervious surface coverage has traditionally been used to examine altered hydrologic response in urban watersheds, several studies suggest that other elements of the urban landscape, particularly those associated with urban infrastructure and the drainage system, play an equally important role. The relative importance of impervious surfaces, stormwater ponds, expansion of the drainage network, and drainage network structures in controlling hydrologic response was examined in the subwatersheds of the Kromma Kill, an urban watershed located in Albany County, NY. In this study, geographic information systems was used to compute geospatial land surface and drainage network properties of 5 Kromma Kill subwatersheds. In these same subwatersheds, water quantity (rainfall and run‐off) and quality (macroinvertebrates, nitrate, total nitrogen, dissolved oxygen, total dissolved solids, and nonpurgable organic carbon) parameters were measured. Strong and significant correlations were identified between land surface and drainage network properties and field observations. Causal relationships were then tested using the Environmental Protection Agency's Stormwater Management Model. Field and model analyses suggest that whereas percent imperviousness is a dominant control on water quality, drainage density and slope are equally important. However, for water quantity, whereas imperviousness is positively correlated with increased run‐off volumes, drainage network properties and slope are the dominant controls on run‐off volumes. Results have important implications for stormwater management plans, especially those aimed at reducing the effective impervious surface coverage of urban watersheds. Reducing the percentage of effective imperviousness in a watershed is not a “one size fits all” solution and can help to meet some management objectives, such as reducing nitrogen concentrations and improving water quality, but may not serve as the most effective, and therefore economical, solution for every management objective including reducing run‐off volumes.     

Aerial flux of particulate hydrocarbons to the Chesapeake Bay Estuary

Bulk hydrocarbon deposition rates have been measured over a 15 month period at four stations in south-eastern Virginia surrounding the lower Chesapeake Bay. A nearly linear trend of atmospheric particulate deposition was recorded. Deposition rates at the urban station (195 μg m^?2 day^?1) were aproximately three times greater than those recorded for nonurban and coastal locations (mean value 69 μg m^?2 day^?1). The increased levels at the urban location were attributed to localized source inputs. Anthropogenic hydrocarbons accounted for approximately 50% of the total deposition at all stations. Significant biogenic inputs were indicated by the odd/even n-alkane distribution. A minimum flux to the water surface, based on mean nonurban deposition rates (24 mg yr^?1), indicated an annual particulate hydrocarbon flux of +275 metric tons. Little information is available for the comparison of additional source inputs; however, the data reported here indicate that the aerial deposition of hydrocarbons is of the same order of magnitude as the input from municipal wastewater facilities and accidental discharge and is a potentially significant source of hydrocarbon pollution to the Chesapeake Bay Estuary.     

Atmospheric deposition to watersheds in complex terrain

Single collection stations for wet or bulk deposition are generally inadequate to describe atmospheric inputs to watersheds in complex terrain. Atmospheric deposition is delivered by wet, dry and cloud deposition processes, and these processes are controlled by a wide range of landscape features, including canopy type and structure, topographic exposure, elevation and slope orientation. As a result, there can be a very high degree of spatial variability within a watershed, and a single sampling point, especially at low elevation, is unlikely to be representative. Atmospheric inputs at the watershed scale can be calculated from the whole watershed mass balance if the outputs and within-watershed sources and sinks are known with sufficient accuracy. Alternatively, indices of atmospheric deposition such as Pb accumulation in the forest floor and SO²⁻₄ flux in throughfall can be used to characterize patterns of total deposition, and these indices can be used to model deposition to the entire watershed based on known landscape features such as elevation and canopy type. © 1997 John Wiley & Sons, Ltd.