Long-term seasonal trends of nitrogen, phosphorus, and suspended sediment load from the non-tidal Susquehanna River Basin to Chesapeake Bay

Legacy sediment as a potential source of orthophosphate: Preliminary conceptual and geochemical models for the Susquehanna River, Chesapeake Bay watershed, USA

Nutrient pollution from agriculture and urban areas plus acid mine drainage (AMD) from legacy coal mines are primary causes of water-quality impairment in the Susquehanna River, which is the predominant source of freshwater and nutrients entering the Chesapeake Bay. Recent increases in the delivery of dissolved orthophosphate (PO4) from the river to the bay may be linked to long-term increases in pH, decreased acidity of precipitation, and decreased acidity, iron, and aluminum loading from widespread AMD. Since the 1950s, baseline pH increased from ~6.5 to ~8 in the West Branch and "North Branch" of the Susquehanna River, which drain bituminous and anthracite coalfields of Pennsylvania. A current baseline pH of ~8 and daily maxima exceeding 9 have been documented along the lower Susquehanna River. In response to improved river quality, bioavailable PO4 now may be released into solution from legacy sediment that has filled major impoundments in lower reaches of the river. At typical pH (5-8) of natural water, aqueous PO4 species tend to be adsorbed by hydrous iron, aluminum, and manganese oxides that coat soil and sediment particles; however, PO4 may be substantially desorbed at pH >8. We created a geochemical model that simulates equilibrium aqueous/solid distributions of PO4 as pH and other solution characteristics change. Considering current conditions in the lower Susquehanna River, the model demonstrates potential for extensive release of adsorbed PO4 at pH >8. Empirical data from laboratory experiments corroborate model results. The transfer of PO4 into the water column may increase algae growth, which removes CO2 and drives pH to higher values, facilitating additional PO4 release and exacerbating the potential for harmful algal blooms. Thus, legacy sediment is a currently unquantified source of PO4 that warrants consideration by resource managers and programs collaborating to reduce phosphorus loads to the bay and similar settings worldwide.

Effects of point and nonpoint source controls on total phosphorus load trends across the Chesapeake Bay watershed, USA

Reduction of total phosphorus (TP) loads has long been a management focus of Chesapeake Bay restoration, but riverine monitoring stations have shown mixed temporal trends. To better understand the regional patterns and drivers of TP trends across the Bay watershed, we compiled and analyzed TP load data from 90 non-tidal network stations using clustering and random forest (RF) approaches. These stations were categorized into two distinct clusters of short-term (2013-2020) TP load trends, i.e. monotonic increase (n = 35) and monotonic decline (n = 55). RF models were developed to identify likely regional drivers of TP trend clusters. Reductions in point sources and agricultural nonpoint sources (i.e. fertilizer) both contributed to water-quality improvement in our period of analysis, thereby demonstrating the effectiveness of nutrient management and the importance of continuing such efforts. In addition, declining TP trends have a larger chance to occur in carbonate areas but a smaller chance in Coastal Plain areas, with the latter likely reflecting the effect of legacy P. To provide spatially explicit information, TP trend clusters were predicted for the entire watershed at the scale of river segments, which are more directly relevant to watershed planning. Among the 975 river segments, 544 (56%) and 431 (44%) were classified as 'monotonic increase' and 'monotonic decrease', respectively. Furthermore, these predicted TP trend clusters were paired with our previously published total nitrogen (TN) trend clusters, showing that TP and TN both declined in 185 segments (19%) and neither declined in 337 segments (35%). Broadly speaking, large-scale nutrient reduction efforts are underway in many regions to curb eutrophication. Water-quality responses and drivers may differ among systems, but our work provides important new evidence on the effectiveness of management efforts toward controlling point and nonpoint sources.

What drove the nonlinear hypoxia response to nutrient loading in Chesapeake Bay during the 20th century?

Previous data analysis showed that the large expansion of hypoxia in Chesapeake Bay between 1950s and 1980s was correlated to the increased riverine nutrient loading, but the physical and biogeochemical processes driving this hypoxia response need to be better understood. Using a validated coupled hydrodynamic-biogeochemical model, we conducted a hindcast simulation of dissolved oxygen during the 40-year period (1950-1989) when the nutrient loading doubled. The model reproduced the observed decline in O₂ concentration at monitoring stations and the expansion of the hypoxic volume. The peak summer hypoxic volume expanded from ∼5 km³ during 1950-1969 to ∼10 km³ during 1970-1989. To discern how different physical and biochemical processes regulated dissolved O₂, we examined O₂ budget in a fixed control volume of the bottom water most susceptible to hypoxia. The increased water column respiration was found to be the dominant driver of the hypoxia expansion. Further analysis showed a nonlinear response to the nutrient loading. The accumulative hypoxia volume days per unit of nitrate load showed an abrupt (∼50 %) jump around 1968. The summer mean hypoxic volume increased with the winter-spring nutrient load, but it was 1.3 km³ (about 30 %) higher in 1968-1989 than in 1950-1967 at the same nutrient load. This upward shift in hypoxia was caused by the upward shift in the relationship between the water column respiration and winter-spring nutrient load. Hypoxia suppressed nitrification and denitrification processes in the sediment, amplifying nutrient recycling by 15 % and water column respiration by 12 %. Our modeling analysis demonstrated a feedback mechanism for driving the nonlinear hypoxia response to nutrient loading.

Nutrient limitation of phytoplankton in three tributaries of Chesapeake Bay: Detecting responses following nutrient reductions

Many coastal ecosystems suffer from eutrophication, algal blooms, and dead zones due to excessive anthropogenic inputs of nitrogen (N) and phosphorus (P). This has led to regional restoration efforts that focus on managing watershed loads of N and P. In Chesapeake Bay, the largest estuary in the United States, dual nutrient reductions of N and P have been pursued since the 1980s. However, it remains unclear whether nutrient limitation - an indicator of restriction of algal growth by supplies of N and P - has changed in the tributaries of Chesapeake Bay following decades of reduction efforts. Toward that end, we analyzed historical data from nutrient-addition bioassay experiments and data from the Chesapeake Bay long-term water-quality monitoring program for six stations in three tidal tributaries (i.e., Patuxent, Potomac, and Choptank Rivers). Classification and regression tree (CART) models were developed using concurrent collections of water-quality parameters for each bioassay monitoring location during 1990-2003, which satisfactorily predicted the bioassay-based measures of nutrient limitation (classification accuracy = 96%). Predictions from the CART models using water-quality monitoring data showed enhanced nutrient limitation over the period of 1985-2020 at four of the six stations, including the downstream station in each of these three tributaries. These results indicate detectable, long-term water-quality improvements in the tidal tributaries. Overall, this research provides a new analytical tool for detecting signs of ecosystem recovery following nutrient reductions. More broadly, the approach can be adapted to other waterbodies with long-term bioassays and water-quality data sets to detect ecosystem recovery.

Land use, hydrology, and climate influence water quality of China's largest river

Influences of multiple environmental factors on water quality patterns is less studied in large rivers. Landscape analysis, multiple statistical methods, and the water quality index (WQI) were used to detect water quality patterns and influencing factors in China's largest river, the Yangtze River. Compared with the dry season, the wet season had significantly higher total phosphorus (TP), chemical oxygen demand (COD), total suspended solids (TSS), and turbidity (TUR). The WQI indicated "Moderate" and "Good" water quality in the wet and dry seasons, respectively. Compared with other sites, the upper reach sites that immediately downstream of the Three Gorges Dam had lower TP, TN, TSS and TUR in both seasons, and had lower and higher water temperature in the wet and dry seasons, respectively. Water quality patterns were mainly driven by heterogeneity in land use (i.e., wetland, cropland, and urban land), hydrology (i.e., water flow, water level), and climate (i.e., rainfall, air temperature). Water quality in the wet season was primarily driven by land use while the joint effect of land use and hydrology primarily drove in the dry season. Decision-makers and regulators of large river basin management may need to develop programs that consider influences from both human and natural drivers for water quality conservation.

Nitrogen reductions have decreased hypoxia in the Chesapeake Bay: Evidence from empirical and numerical modeling

Seasonal hypoxia is a characteristic feature of the Chesapeake Bay due to anthropogenic nutrient input from agriculture and urbanization throughout the watershed. Although coordinated management efforts since 1985 have reduced nutrient inputs to the Bay, oxygen concentrations at depth in the summer still frequently fail to meet water quality standards that have been set to protect critical estuarine living resources. To quantify the impact of watershed nitrogen reductions on Bay hypoxia during a recent period including both average discharge and extremely wet years (2016-2019), this study employed both statistical and three-dimensional (3-D) numerical modeling analyses. Numerical model results suggest that if the nitrogen reductions since 1985 had not occurred, annual hypoxic volumes (O₂ < 3 mg L^-1) would have been ~50-120% greater during the average discharge years of 2016-2017 and ~20-50% greater during the wet years of 2018-2019. The effect was even greater for O₂ < 1 mg L^-1, where annual volumes would have been ~80-280% greater in 2016-2017 and ~30-100% greater in 2018-2019. These results were supported by statistical analysis of empirical data, though the magnitude of improvement due to nitrogen reductions was greater in the numerical modeling results than in the statistical analysis. This discrepancy is largely accounted for by warming in the Bay that has exacerbated hypoxia and offset roughly 6-34% of the improvement from nitrogen reductions. Although these results may reassure policymakers and stakeholders that their efforts to reduce hypoxia have improved ecosystem health in the Bay, they also indicate that greater reductions are needed to counteract the ever-increasing impacts of climate change.

Nutrient limitation of phytoplankton in Chesapeake Bay: Development of an empirical approach for water-quality management

Understanding the temporal and spatial roles of nutrient limitation on phytoplankton growth is necessary for developing successful management strategies. Chesapeake Bay has well-documented seasonal and spatial variations in nutrient limitation, but it remains unknown whether these patterns of nutrient limitation have changed in response to nutrient management efforts. We analyzed historical data from nutrient bioassay experiments (1992-2002) and data from long-term, fixed-site water-quality monitoring program (1990-2017) to develop empirical approaches for predicting nutrient limitation in the surface waters of the mainstem Bay. Results from classification and regression trees (CART) matched the seasonal and spatial patterns of bioassay-based nutrient limitation in the 1992-2002 period much better than two simpler, non-statistical approaches. An ensemble approach of three selected CART models satisfactorily reproduced the bioassay-based results (classification rate = 99%). This empirical approach can be used to characterize nutrient limitation from long-term water-quality monitoring data on much broader geographic and temporal scales than would be feasible using bioassays, providing a new tool for informing water-quality management. Results from our application of the approach to 21 tidal monitoring stations for the period of 2007-2017 showed modest changes in nutrient limitation patterns, with expanded areas of nitrogen-limitation and contracted areas of nutrient saturation (i.e., not limited by nitrogen or phosphorus). These changes imply that long-term reductions in nitrogen load have led to expanded areas with nutrient-limited phytoplankton growth in the Bay, reflecting long-term water-quality improvements in the context of nutrient enrichment. However, nutrient limitation patterns remain unchanged in the majority of the mainstem, suggesting that nutrient loads should be further reduced to achieve a less nutrient-saturated ecosystem.

A simulation-based approach to assess the power of trend detection in high- and low-frequency water quality records

To provide more precise understanding of water quality changes, continuous sampling is being used more in surface water quality monitoring networks. However, it remains unclear how much improvement continuous monitoring provides over spot sampling, in identifying water quality changes over time. This study aims (1) to assess our ability to detect trends using water quality data of both high and low frequencies and (2) to assess the value of using high-frequency data as a surrogate to help detect trends in other constituents. Statistical regression models were used to identify temporal trends and then to assess the trend detection power of high-frequency (15 min) and low-frequency (monthly) data for turbidity and electrical conductivity (EC) data collected across Victoria, Australia. In addition, we developed surrogate models to simulate five sediment and nutrients constituents from runoff, turbidity and EC. A simulation-based statistical approach was then used to the compare the power to detect trends between the low- and high-frequency water quality records. Results show that high-frequency sampling shows clear benefits in trend detection power for turbidity, EC, as well as simulated sediment and nutrients, especially over short data periods. For detecting a 1% annual trend with 5 years of data, up to 97% and 94% improvements on the trend detection probability are offered by high-frequency data compared with monthly data, for turbidity and EC, respectively. Our results highlight the benefits of upgrading monitoring networks with wider application of high-frequency sampling.

Trends in water quality in a subtropical Australian river-estuary system: Responses to damming, climate variability and wastewater discharges

The Logan-Albert estuary in southeast Queensland, Australia, has high biodiversity and supports multiple economic and recreational services. Elevated nutrient and sediment loads have been a longstanding management issue for the estuary. We investigated the spatial and seasonal patterns of nutrients and turbidity along the Logan-Albert estuary and assessed the effects of a recently constructed upstream dam. Nutrient concentrations and turbidity levels were analysed using 15 years of monitoring data from 19 water quality sites throughout the estuary. We hypothesised that the construction of Wyaralong Dam would act as a nutrient and sediment sink which may have positive effects on downstream water quality. Long-term trends of water quality constituents were evaluated using a non-parametric seasonal Mann-Kendall test and the effect of upstream impoundment was assessed with a Before-After Control-Impact (BACI) test. Nutrient concentrations and turbidity levels declined significantly with time in the upper Logan estuary and, to a lesser extent, in the lower Albert estuary. The general improvement of water quality in the upper Logan estuary was attributed to construction of the Wyaralong Dam. Significant decreases in concentrations of total phosphorus (TP) and oxidised nitrogen (NO_x-N) along the lower Albert were principally attributed to wetter conditions over the 15-year dataset, which diluted point-source loads from a nearby wastewater treatment plant (WWTP). Our results show that estuarine water quality changes can be highly dynamic with interactions amongst climate and management practices that necessitate long-term monitoring programs with good spatial coverage.

Factors driving nutrient trends in streams of the Chesapeake Bay watershed

Despite decades of effort toward reducing nitrogen and phosphorus flux to Chesapeake Bay, water-quality and ecological responses in surface waters have been mixed. Recent research, however, provides useful insight into multiple factors complicating the understanding of nutrient trends in bay tributaries, which we review in this paper, as we approach a 2025 total maximum daily load (TMDL) management deadline. Improvements in water quality in many streams are attributable to management actions that reduced point sources and atmospheric nitrogen deposition and to changes in climate. Nutrient reductions expected from management actions, however, have not been fully realized in watershed streams. Nitrogen from urban nonpoint sources has declined, although water-quality responses to urbanization in individual streams vary depending on predevelopment land use. Evolving agriculture, the largest watershed source of nutrients, has likely contributed to local nutrient trends but has not affected substantial changes in flux to the bay. Changing average nitrogen yields from farmland underlain by carbonate rocks, however, may suggest future trends in other areas under similar management, climatic, or other influences, although drivers of these changes remain unclear. Regardless of upstream trends, phosphorus flux to the bay from its largest tributary has increased due to sediment infill in the Conowingo Reservoir. In general, recent research emphasizes the utility of input reductions over attempts to manage nutrient fate and transport at limiting nutrients in surface waters. Ongoing research opportunities include evaluating effects of climate change and conservation practices over time and space and developing tools to disentangle and evaluate multiple influences on regional water quality.

The Latitudes, Attitudes, and Platitudes of Watershed Phosphorus Management in North America

Phosphorus (P) plays a crucial role in agriculture as a primary fertilizer nutrient-and as a cause of the eutrophication of surface waters. Despite decades of efforts to keep P on agricultural fields and reduce losses to waterways, frequent algal blooms persist, triggering not only ecological disruption but also economic, social, and political consequences. We investigate historical and persistent factors affecting agricultural P mitigation in a transect of major watersheds across North America: Lake Winnipeg, Lake Erie, the Chesapeake Bay, and Lake Okeechobee/Everglades. These water bodies span 26 degrees of latitude, from the cold climate of central Canada to the subtropics of the southeastern United States. These water bodies and their associated watersheds have tracked trajectories of P mitigation that manifest remarkable similarities, and all have faced challenges in the application of science to agricultural management that continue to this day. An evolution of knowledge and experience in watershed P mitigation calls into question uniform solutions as well as efforts to transfer strategies from other arenas. As a result, there is a need to admit to shortcomings of past approaches, plotting a future for watershed P mitigation that accepts the sometimes two-sided nature of Hennig Brandt's "Devil's Element."

Evaluation of the total maximum allocated load of dissolved inorganic nitrogen using a watershed-coastal ocean coupled model

Due to the recent rapid increase in human activity and economic development, many coastal areas have recently experienced a high degree of land-based pollution. Evaluating the total maximum allocated load (TMAL) of dissolved inorganic nitrogen (DIN) nutrients and the remaining capacity is of importance for improving water quality. A considerable amount of nutrients derived from the coastal watershed can be found in wet seasons, which is non-negligible for the estimation of remaining capacity. Therefore, we use a watershed-coastal ocean coupled model combined with an optimization algorithm to tackle this issue. In contrast with previous studies, this study provides a method to estimate the spatiotemporal variations in TMALs and we then compare it to the current DIN nutrient load, including both point sources and non-point sources. Our results suggest that the TMAL of Daya Bay (DB), which is located in the northern part of the South China Sea, is about 7976 metric tons per year (t/yr) and ranges from 191 metric tons per month (t/month) to 1072 t/month. The increase of non-point source (NPS) DIN input also plays an important role in daily overload events during wet seasons. Moreover, the TMALs show an inverse exponential correlation with the water age, but only about 65% of the variance is explained. This suggests that the variations from the optimization algorithm and from local water function zoning plans are also important. According to our prediction of the DIN input, the TMAL of DB will soon be exhausted in the next several years. Consequently, prompt actions are necessary to consider the distribution of TMALs in urban developments and to decelerate the rapid growth of DIN input. Therefore, the results of this study will be helpful for both local pollution control and future urban planning.

Estuarine salinity recovery from an extreme precipitation event: Hurricane Harvey in Galveston Bay

With a warming climate and a more humid atmosphere, extreme precipitation events are projected to occur more frequently in future. Understanding how coastal systems respond to and recover from such acute events is of fundamental significance for environmental assessment and management. A hydrodynamic model was used to examine the estuarine responses in Galveston Bay to Hurricane Harvey, an extreme precipitation event with a return period of larger than 1000 years. The enormous freshwater input during Harvey caused long-lasting elevated water level, extraordinarily strong along-channel velocity, sharp decreases in salinity, and huge river plumes, all of which were well reproduced by the model. The salinity recovery time (T_R) was estimated as a measure of the system resiliency to stormwater input. Over the entire bay, the T_R had a mean of two months, but with great variability ranging from less than 10 days near the bay entrance to over three months in the inner part of Trinity Bay and the middle of East Bay. The spatially varying T_R was explained by different contributions of exchange flow and tidal pumping to salt flux. At the bay entrance, tidal pumping facilitated by the shelf current was the dominant mechanism for salt influx, while exchange flow and tidal pumping had a comparable contribution to salt influx to Trinity Bay. The spatial pattern of the T_R appears consistent with the changes in the phytoplankton community in the bay. A series of numerical experiments with different amounts of stormwater reveals a non-linear relationship between the bay-wide mean T_R and the amount of stormwater, with the rate of increase in T_R decreasing when stormwater input increases. The present approach using a hydrodynamic model will be able to provide a quick assessment of the environmental pressure from extreme events.

A REVIEW OF WATER QUALITY RESPONSES TO AIR TEMPERATURE AND PRECIPITATION CHANGES 2: NUTRIENTS, ALGAL BLOOMS, SEDIMENT, PATHOGENS

In this paper we review the published, scientific literature addressing the response of nutrients, sediment, pathogens and cyanobacterial blooms to historical and potential future changes in air temperature and precipitation. The goal is to document how different attributes of water quality are sensitive to these drivers, to characterize future risk, to inform management responses and to identify research needs to fill gaps in our understanding. Results suggest that anticipated future changes present a risk of water quality and ecosystem degradation in many U.S. locations. Understanding responses is, however, complicated by inherent high spatial and temporal variability, interactions with land use and water management, and dependence on uncertain changes in hydrology in response to future climate. Effects on pollutant loading in different watershed settings generally correlate with projected changes in precipitation and runoff. In all regions, increased heavy precipitation events are likely to drive more episodic pollutant loading to water bodies. The risk of algal blooms could increase due to an expanded seasonal window of warm water temperatures and the potential for episodic increases in nutrient loading. Increased air and water temperatures are also likely to affect the survival of waterborne pathogens. Responding to these challenges requires understanding of vulnerabilities, and management strategies to reduce risk.

Chesapeake Bay's water quality condition has been recovering: Insights from a multimetric indicator assessment of thirty years of tidal monitoring data

To protect the aquatic living resources of Chesapeake Bay, the Chesapeake Bay Program partnership has developed guidance for state water quality standards, which include ambient water quality criteria to protect designated uses (DUs), and associated assessment procedures for dissolved oxygen (DO), water clarity/underwater bay grasses, and chlorophyll-a. For measuring progress toward meeting the respective states' water quality standards, a multimetric attainment indicator approach was developed to estimate combined standards attainment. We applied this approach to three decades of monitoring data of DO, water clarity/underwater bay grasses, and chlorophyll-a data on annually updated moving 3-year periods to track the progress in all 92 management segments of tidal waters in Chesapeake Bay. In 2014-2016, 40% of tidal water segment-DU-criterion combinations in the Bay (n = 291) are estimated to meet thresholds for attainment of their water quality criteria. This index score marks the best 3-year status in the entire record. Since 1985-1987, the indicator has followed a nonlinear trajectory, consistent with impacts from extreme weather events and subsequent recoveries. Over the period of record (1985-2016), the indicator exhibited a positive and statistically significant trend (p < 0.05), indicating that the Bay has been recovering since 1985. Patterns of attainment of individual DUs are variable, but improvements in open water DO, deep channel DO, and water clarity/submerged aquatic vegetation have combined to drive the improvement in the Baywide indicator in 2014-2016 relative to its long-term median. Finally, the improvement in estimated Baywide attainment was statistically linked to the decline of total nitrogen, indicating responsiveness of attainment status to the reduction of nutrient load through various management actions since at least the 1980s.

Agricultural conservation practices can help mitigate the impact of climate change

Agricultural conservation practices (CPs) are commonly implemented to reduce diffuse nutrient pollution. Climate change can complicate the development, implementation, and efficiency of agricultural CPs by altering hydrology, nutrient cycling, and erosion. This research quantifies the impact of climate change on hydrology, nutrient cycling, erosion, and the effectiveness of agricultural CP in the Susquehanna River Basin in the Chesapeake Bay Watershed, USA. We develop, calibrate, and test the Soil and Water Assessment Tool-Variable Source Area (SWAT-VSA) model and select four CPs; buffer strips, strip-cropping, no-till, and tile drainage, to test their effectiveness in reducing climate change impacts on water quality. We force the model with six downscaled global climate models (GCMs) for a historic period (1990-2014) and two future scenario periods (2041-2065 and 2075-2099) and quantify the impact of climate change on hydrology, nitrate-N (NO₃-N), total N (TN), dissolved phosphorus (DP), total phosphorus (TP), and sediment export with and without CPs. We also test prioritizing CP installation on the 30% of agricultural lands that generate the most runoff (e.g., critical source areas-CSAs). Compared against the historical baseline and with no CPs, the ensemble model predictions indicate that climate change results in annual increases in flow (4.5±7.3%), surface runoff (3.5±6.1%), sediment export (28.5±18.2%) and TN export (9.5±5.1%), but decreases in NO₃-N (12±12.8%), DP (14±11.5), and TP (2.5±7.4%) export. When agricultural CPs are simulated most do not appreciably change the water balance, however, tile drainage and strip-cropping decrease surface runoff, sediment export, and DP/TP, while buffer strips reduce N export. Installing CPs on CSAs results in nearly the same level of performance for most practices and most pollutants. These results suggest that climate change will influence the performance of agricultural CPs and that targeting agricultural CPs to CSAs can provide nearly the same level of water quality effects as more widespread adoption.

Declining ambient water phosphorus concentrations in Massachusetts' rivers from 1999 to 2013: Environmental protection works

Over the last century, nutrient concentrations in streams, rivers, lakes and ponds have increased substantially in the United States. Elevated phosphorus levels are a concern due to their ability to cause changes in freshwater ecosystems that are detrimental to humans and wildlife. In the present study, long-term trends in total phosphorus (TP) concentrations from 20 rivers in central Massachusetts from 1999 to 2013 were investigated. Kendall's correlation coefficients were used to demonstrate that 18 of the 20 rivers had significant reductions in TP concentrations (P < 0.05). A similar trend was found when flow-adjusted TP concentrations were analyzed. At the beginning of monitoring activities, the average TP concentration in 9 of the 20 rivers was greater than 0.05 mg/L and 6 of these 9 rivers contained TP concentrations greater than 0.1 mg/L; about fifteen years later, only 3 rivers contained TP greater than 0.05 mg/L and none had concentrations> 0.1 mg/L. TP decreases were greater in rivers with more anthropogenic inputs. Principal component analysis (PCA) revealed that the decline of TP in these Massachusetts streams is likely the result of advancements in wastewater treatment and implementation of effective non-point source management practices.

Watershed export of fine sediment, organic carbon, and chlorophyll-a to Chesapeake Bay: Spatial and temporal patterns in 1984-2016

Chesapeake Bay has long experienced nutrient enrichment and water clarity deterioration. This study provides new quantification of loads and yields for sediment (fine and coarse grained), organic carbon (total, dissolved, and particulate), and chlorophyll-a from the monitored nontidal Chesapeake Bay watershed (MNTCBW), all of which are expected to drive estuarine water clarity. We conducted an integrated analysis of nine major tributaries to the Bay to understand spatial and temporal export patterns over the last thirty years (1984-2016). In terms of spatial pattern, export of these constituents from the MNTCBW was strongly dominated (~90%) by the three largest tributaries (i.e., Susquehanna, Potomac, and James). Among the nine tributaries, the ranking of constituent export generally follows the order of their watershed sizes, with other factors such as land use and reservoir playing important roles in some exceptions. In terms of partitioning, suspended sediment (SS) export was dominated by fine-grained sediment (SS_fine) in all nine tributaries; overall, ~90% of the MNTCBW SS is SS_fine. Total organic carbon (TOC) export was dominated by dissolved organic carbon (DOC) in all tributaries except Potomac River; overall, ~60% of the MNTCBW TOC is DOC. A comparison with literature shows that the MNTCBW SS and TOC yields were ~80% and ~60% of the respective medians of worldwide watersheds. In terms of temporal pattern, flow-normalized yields from the MNTCBW show overall increases in SS (both long-term [1984-2016] and short-term [2004-2016]), SS_fine (long-term and short-term), TOC (long-term), and chlorophyll-a (short-term). The rises in SS, SS_fine, and TOC were largely driven by Susquehanna River where Conowingo Reservoir's trapping efficiency has greatly diminished in the last twenty years. Overall, these new results on the status and trends of sediment, organic carbon, and chlorophyll-a provide the foundation for building potential linkages between riverine inputs and estuarine water clarity patterns.

Synthesis of nutrient and sediment export patterns in the Chesapeake Bay watershed: Complex and non-stationary concentration-discharge relationships

Derived from river monitoring data, concentration-discharge (C-Q) relationships are useful indicators of riverine export dynamics. A top-down synthesis of C-Q patterns was conducted for suspended sediment (SS), total phosphorus (TP), and total nitrogen (TN) for nine major tributaries (15 monitoring sites) to Chesapeake Bay, which represent diverse characteristics in terms of land use, physiography, and hydrological settings. Model coefficients from the recently-developed Weighted Regressions on Time, Discharge, and Season (WRTDS) method were used to make informative interpretation of C-Q relationships. Unlike many previous C-Q studies that focused on stormflow conditions, this approach allows simultaneous examination of various discharge conditions within an uncertainty framework. This synthesis on WRTDS coefficients (i.e., the sensitivity of concentration to discharge) has offered new insights on the complexity of watershed function. Results show that watershed export has been dominated by mobilization patterns for SS and TP (particulate-dominated species) and chemostasis patterns for TN (dissolved-dominated species) under many river discharge conditions. Among nine possible modalities of low-flow vs. high-flow patterns, the three most frequent modalities are mobilization vs. mobilization (17 cases), chemostasis vs. mobilization (13 cases), and chemostasis vs. chemostasis (7 cases), representing 82% of all 45 watershed-constituent pairs. The general lack of dilution patterns may suggest that none of these constituents has been supply-limited in these watersheds. For many watershed-constituent combinations, results show clear temporal non-stationarity in C-Q relationships under selected time-invariant discharges, reflecting major changes in dominant watershed sources due to anthropogenic actions. These results highlight the potential pitfalls of assuming fixed C-Q relationships in the record. Overall, this work demonstrates the utility of WRTDS model coefficients for interpretation of river water-quality data and for generation of sensible hypotheses on dominant processes in different watersheds. The approach is readily adaptable to other river systems, where long-term discretely-sampled data are available, to decipher complex interactions between hydrological and biogeochemical processes.

Chesapeake Bay Dissolved Oxygen Criterion Attainment Deficit: Three Decades of Temporal and Spatial Patterns

Low dissolved oxygen (DO) conditions are a recurring issue in waters of Chesapeake Bay, with detrimental effects on aquatic living resources. The Chesapeake Bay Program partnership has developed criteria guidance supporting the definition of state water quality standards and associated assessment procedures for DO and other parameters, which provides a binary classification of attainment or impairment. Evaluating time series of these two outcomes alone, however, provides limited information on water quality change over time or space. Here we introduce an extension of the existing Chesapeake Bay water quality criterion assessment framework to quantify the amount of impairment shown by space-time exceedance of DO criterion ("attainment deficit") for a specific tidal management unit (i.e., segment). We demonstrate the usefulness of this extended framework by applying it to Bay segments for each 3-year assessment period between 1985 and 2016. In general, the attainment deficit for the most recent period assessed (i.e., 2014-2016) is considerably worse for deep channel (DC; n = 10) segments than open water (OW; n = 92) and deep water (DW; n = 18) segments. Most subgroups - classified by designated uses, salinity zones, or tidal systems - show better (or similar) attainment status in 2014-2016 than their initial status (1985-1987). Some significant temporal trends (p < 0.1) were detected, presenting evidence on the recovery for portions of Chesapeake Bay with respect to DO criterion attainment. Significant, improving trends were observed in seven OW segments, four DW segments, and one DC segment over the 30 3-year assessment periods (1985-2016). Likewise, significant, improving trends were observed in 15 OW, five DW, and four DC segments over the recent 15 assessment periods (2000-2016). Subgroups showed mixed trends, with the Patuxent, Nanticoke, and Choptank Rivers experiencing significant, improving short-term (2000-2016) trends while Elizabeth experiencing a significant, degrading short-term trend. The general lack of significantly improving trends across the Bay suggests that further actions will be necessary to achieve full attainment of DO criterion. Insights revealed in this work are critical for understanding the dynamics of the Bay ecosystem and for further assessing the effectiveness of management initiatives aimed toward Bay restoration.

Watershed 'Chemical Cocktails': Forming Novel Elemental Combinations in Anthropocene Fresh Waters

In the Anthropocene1, watershed chemical transport is increasingly dominated by novel combinations elements, which are hydrologically linked together as 'chemical cocktails.' Chemical cocktails are novel because human activities greatly enhance elemental concentrations and their probability for biogeochemical interactions and shared transport along hydrologic flowpaths. A new chemical cocktail approach advances our ability to: trace contaminant mixtures in watersheds, develop chemical proxies with high-resolution sensor data, and manage multiple water quality problems. We explore the following questions: (1) Can we classify elemental transport in watersheds as chemical cocktails using a new approach? (2) What is the role of climate and land use in enhancing the formation and transport of chemical cocktails in watersheds? To address these questions, we first analyze trends in concentrations of carbon, nutrients, metals, and salts in fresh waters over 100 years. Next, we explore how climate and land use enhance the probability of formation of chemical cocktails of carbon, nutrients, metals, and salts. Ultimately, we classify transport of chemical cocktails based on solubility, mobility, reactivity, and dominant phases: (1) sieved chemical cocktails (e.g., particulate forms of nutrients, metals and organic matter); (2) filtered chemical cocktails (e.g., dissolved organic matter and associated metal complexes); (3) chromatographic chemical cocktails (e.g., ions eluted from soil exchange sites); and (4) reactive chemical cocktails (e.g., limiting nutrients and redox sensitive elements). Typically, contaminants are regulated and managed one element at a time, even though combinations of elements interact to influence many water-quality problems such as toxicity to life, eutrophication, infrastructure and water treatment. A chemical cocktail approach significantly expands evaluations of water-quality signatures and impacts beyond single elements to mixtures. High-frequency sensor data (pH, specific conductance, turbidity, etc.) can serve as proxies for chemical cocktails and improve real-time analyses of water-quality violations, identify regulatory needs, and track water quality recovery following and extreme climate events. Ultimately, a watershed chemical cocktail approach is necessary for effectively co-managing groups of contaminants and provides a more holistic approach for studying, monitoring, and managing water quality in the Anthropocene.

A SWAT model validation of nested-scale contemporaneous stream flow, suspended sediment and nutrients from a multiple-land-use watershed of the central USA

There is an ongoing need to validate the accuracy of predictive model simulated pollutant yields, particularly from multiple-land-use (i.e. forested, agricultural, and urban) watersheds. However, there are seldom sufficient observed data sets available that supply requisite spatial and temporal resolution and coupled multi-parameter constituents for rigorous model performance assessment. Four years of hydroclimate and water quality data were used to validate SWAT model estimates of monthly stream flow, suspended sediment, total phosphorus, nitrate, nitrite, ammonium, and total inorganic nitrogen from 5 nested-scale gauging sites located in a multiple-land-use watershed of the central USA. The uncalibrated SWAT model satisfactorily simulated monthly stream flow with Nash-Sutcliffe efficiency (NSE) values ranging from 0.50 near the headwaters, to 0.75 near the watershed outlet. However, the uncalibrated model did not accurately simulate monthly sediment, total phosphorus, nitrate, nitrite, ammonium, and total inorganic nitrogen with NSE values<0.05. Calibrating the SWAT model to multiple gauging sites within the watershed improved estimates of monthly stream flow (NSE=0.83), sediment (NSE=0.78), total phosphorus (NSE=0.81), nitrate (NSE=0.90), and total inorganic nitrogen (NSE=0.86). However, NSE values were <-0.16 for nitrite and ammonium estimates. Additionally, model performance decreased for sediment, nitrate, and total inorganic nitrogen during the validation period with NSE values<0.62, 0.52, and 0.36, respectively. Results highlight the benefits of calibrating the SWAT model to multiple gauging sites and provide guidance to SWAT model (or similar models) users wishing to improve model performance at multiple scales.

Decadal-scale export of nitrogen, phosphorus, and sediment from the Susquehanna River basin, USA: Analysis and synthesis of temporal and spatial patterns

The export of nitrogen (N), phosphorus (P), and suspended sediment (SS) is a long-standing management concern for the Chesapeake Bay watershed, USA. Here we present a comprehensive evaluation of nutrient and sediment loads over the last three decades at multiple locations in the Susquehanna River basin (SRB), Chesapeake's largest tributary watershed. Sediment and nutrient riverine loadings, including both dissolved and particulate fractions, have generally declined at all sites upstream of Conowingo Dam (non-tidal SRB outlet). Period-of-record declines in riverine yield are generally smaller than those in source input, suggesting the possibility of legacy contributions. Consistent with other watershed studies, these results reinforce the importance of considering lag time between the implementation of management actions and achievement of river quality improvement. Whereas flow-normalized loadings for particulate species have increased recently below Conowingo Reservoir, those for upstream sites have declined, thus substantiating conclusions from prior studies about decreased reservoir trapping efficiency. In regard to streamflow effects, statistically significant log-linear relationships between annual streamflow and annual constituent load suggest the dominance of hydrological control on the inter-annual variability of constituent export. Concentration-discharge relationships revealed general chemostasis and mobilization effects for dissolved and particulate species, respectively, both suggesting transport-limitation conditions. In addition to affecting annual export rates, streamflow has also modulated the relative importance of dissolved and particulate fractions, as reflected by its negative correlations with dissolved P/total P, dissolved N/total N, particulate P/SS, and total N/total P ratios. For land-use effects, period-of-record median annual yields of N, P, and SS all correlate positively with the area fraction of non-forested land but negatively with that of forested land under all hydrological conditions. Overall, this work has informed understanding with respect to four major factors affecting constituent export (i.e., source input, reservoir modulation, streamflow, and land use) and demonstrated the value of long-term river monitoring.

Influence of Reservoir Infill on Coastal Deep Water Hypoxia

Ecological restoration of the Chesapeake through the Chesapeake Bay total maximum daily load (TMDL) requires the reduction of nitrogen, phosphorus, and sediment loads in the Chesapeake watershed because of the tidal water quality impairments and damage to living resources they cause. Within the Chesapeake watershed, the Conowingo Reservoir has been filling in with sediment for almost a century and is now in a state of near-full capacity called . The development of the Chesapeake TMDL in 2010 was with the assumption that the Conowingo Reservoir was still effectively trapping sediment and nutrients. This is now known not to be the case. In a TMDL, pollutant loads beyond the TMDL allocation, which are brought about by growth or other conditions, must be offset. Using the analysis tools of the Chesapeake TMDL for assessing the degree of water quality standard attainment, the estimated nutrient and sediment loads from a simulated dynamic equilibrium infill condition of the Conowingo Reservoir were determined. The influence on Chesapeake water quality by a large storm and scour event of January 1996 on the Susquehanna River was estimated, and the same storm and scour events were also evaluated in the more critical living resource period of June. An analysis was also made on the estimated influence of more moderate high flow events. The infill of the Conowingo reservoir had estimated impairments of water quality, primarily on deep-water and deep-channel dissolved oxygen, because of increased discharge and transport of organic and particulate inorganic nutrients from the Conowingo Reservoir.

Long-Term Changes in Sediment and Nutrient Delivery from Conowingo Dam to Chesapeake Bay: Effects of Reservoir Sedimentation

Reduction of suspended sediment (SS), total phosphorus (TP), and total nitrogen is an important focus for Chesapeake Bay watershed management. The Susquehanna River, the bay's largest tributary, has drawn attention because SS loads from behind Conowingo Dam (near the river's mouth) have been rising dramatically. To better understand these changes, we evaluated histories of concentration and loading (1986-2013) using data from sites above and below Conowingo Reservoir. First, observed concentration-discharge relationships show that SS and TP concentrations at the reservoir inlet have declined under most discharges in recent decades, but without corresponding declines at the outlet, implying recently diminished reservoir trapping. Second, best estimates of mass balance suggest decreasing net deposition of SS and TP in recent decades over a wide range of discharges, with cumulative mass generally dominated by the 75∼99.5th percentile of daily Conowingo discharges. Finally, stationary models that better accommodate effects of riverflow variability also support the conclusion of diminished trapping of SS and TP under a range of discharges that includes those well below the literature-reported scour threshold. Overall, these findings suggest that decreased net deposition of SS and TP has occurred at subscour levels of discharge, which has significant implications for the Chesapeake Bay ecosystem.