Introduction

The Chesapeake Bay is an estuary which has undergone considerable water quality degradation from anthropogenic point and diffuse sources of N and P (Cooper and Brush 1993; Kemp et al. 2005). N- and P-rich sources from the atmosphere and land impair use of receiving waters for drinking and recreation, result in algal blooms and hypoxia (Officer et al. 1984; Nixon 1995; Glibert et al. 1995; Harding Jr. and Perry 1997; Follett and Follett 2001; Dagg and Breed 2003; Hagy et al. 2004; Kemp et al. 2005; Harding et al. 2019), and often result in losses of marshlands and submerged aquatic grasses (Orth and Moore 1983; Davis 1985; Deegan et al. 2012). In rural watersheds, most diffuse sources of N and P are dominated by agriculture (e.g., David et al. 1997; Lee et al. 2001; van Grieken et al. 2019), whereas in urbanized watersheds, N and P inputs are primarily from point sources (e.g., Groffman et al. 2003).

Anthropogenic sources of N and P can be controlled by a variety of mechanisms. Point sources in most developed countries are typically regulated by environmental agencies, and the costs of nutrient reductions are borne by the industry or municipality via the tax base. In the USA, point sources have discharge permits allocating specific pollution caps (National Pollution Discharge Elimination System (NPDES)) and are self-monitored, with data reported to a state or federal environmental agency. Likewise, control of agricultural diffuse sources in developed countries is usually encouraged by government agencies via adoption of agricultural best management practices (BMPs, Ches. Bay Found. 2019). These are typically not monitored, and estimated reductions are based on previous field trials. Agricultural BMPs are conservation practices, such as fertilizer management, drainage control structures, or winter cover crops, that are intended to reduce losses of N and P from agricultural fields to surface runoff or groundwater but may not increase cash crop yields. BMPs are largely voluntary and may have benefits for the public but not necessarily for farmers. As a result, BMPs are subsidized by agencies such as the US Dept. of Agriculture and state agencies with incentive payments and “cost-shares” to partially off-set costs to the farmer. Reductions in diffuse nutrient pollution are more challenging and are typically quantified by monitoring of atmospheric deposition or single land uses (e.g., Reckhow et al. 1980; Beaulac and Reckhow 1982; Clesceri et al. 1986; Clark et al. 2000; Schlesinger and Bernhardt 2014), by comparing watersheds with varying amounts of land uses (e.g., Jordan et al. 2003; Fisher et al. 2006a, 2010), or by hydrochemical modeling (e.g., Lee et al. 2001; Alexander et al. 2002).

Recognition of the symptoms of eutrophication in Chesapeake Bay and its causes led to considerable efforts to monitor and reduce the inputs of N and P (EPA 2010). A total maximum daily load (TMDL) for the Chesapeake was established by the US EPA in 2010 to address degraded water quality, with a goal to achieve by 2025, a 25%, 24%, and 20% reduction in N, P, and total suspended solids (TSS), respectively (EPA 2017). Efforts include a major monitoring program in the watershed and estuary (EPA 2018), development of water quality criteria (Harding Jr. et al. 2014), wastewater treatment upgrades (EPA 2016), and implementation of agricultural and urban BMPs (e.g., Hively et al. 2018).

Widespread improvements in water quality have yet to occur in Chesapeake Bay (Hirsch et al. 2010; Liu et al. 2017; Ator et al. 2019; Moyer and Blomquist 2020). However, parts of the Chesapeake have experienced improved water quality, and Gurbisz et al. (2016) and Lefcheck et al. (2018) provided evidence that a partial recovery of submerged grasses in Chesapeake Bay was due to nutrient reductions from watershed sources. In addition, large reductions in sewage N and P from the Blue Plains wastewater plant servicing Washington DC and discharging into the Potomac estuary resulted in recovery of submerged grasses in shallow tidal fresh to oligohaline areas of the estuary (Jaworski et al. 2007). Furthermore, Eshleman et al. (2013) provided evidence that diffuse source reductions in the non-tidal Potomac River are related to reductions in atmospheric N deposition resulting from the 1990 US Clean Air Act Acid Rain Program amendment and other regulatory programs. However, despite these signs of recovery in the tidal fresh Potomac, the tidal saline Potomac estuary shows little sign of change in algal biomass (measured as chlorophyll a) and bottom water dissolved oxygen (Bricker et al. 2014).

The Choptank is a tributary of Chesapeake Bay on the Delmarva Peninsula (Fig. 1). For the last 20 years, we have worked with farmers, wastewater treatment plant operators, government agencies, and water quality groups to encourage conservation efforts and to discern trends in water quality in the Choptank basin. Here, we evaluate whether the TMDL and other enhanced management of atmospheric deposition, point sources, and agriculture have led to improved water quality in streams and in the Choptank estuary. Our main hypothesis is that management efforts to reduce inputs from the three known nutrient sources in the Choptank watershed (atmospheric deposition, agriculture, waste water) have improved water quality in the Choptank estuary.

Fig. 1
figure 1

The Choptank Basin in the Mid-Atlantic region. The USGS gauge at Greensboro MD is shown, along with the gauged watersheds, the nearby forested watershed, the NADP precipitation monitoring station, the permitted wastewater treatment plants, and the EPA Bay Program estuarine monitoring stations EE2.1, ET5.2, and ET5.1

Study Site Description

The Choptank watershed and estuary are components of a coastal plain tributary draining into Chesapeake Bay (Fig. 1). The Choptank watershed lies primarily in the state of Maryland, with a portion in Delaware (Fig. 1). The area has a humid, temperate climate with an average annual precipitation of 1.12 m year−1. In the estuary, semi-diurnal tides with an amplitude of ~ 0.5 m in the lower estuary are amplified by occasional wind tides driven by on-shore and off-shore components of wind. Tides in the upper basin are amplified by the shape of the estuary to as much as 1 m and penetrate as far upstream as the USGS gauging station at Greensboro MD (Fig. 1).

Eutrophication of the Choptank estuary (Fisher et al. 2006a) is a microcosm of eutrophication in Chesapeake Bay as a whole. For example, agriculture and forests are the dominant land uses in both the Choptank and Chesapeake basins, and population density in the Choptank (59 km−2 in 2000, 52 km−2 in 2010; Fisher et al. 2006a) is similar to other Chesapeake basins (30–70 km−2; Carlozo et al. 2008). As for the Chesapeake (e.g., Kemp et al. 2005), the degradation of the Choptank estuary is well documented (e.g., Staver and Brinsfield 1996; Newell et al. 2004; Benitez and Fisher 2004; Fisher et al. 2006a, b; Sutton et al. 2009). The Choptank watershed and estuary are also focus areas for programs of the US Department of Agriculture (USDA) and US National Oceanographic and Atmospheric Administration (NOAA).

The Choptank basin is well monitored. There is a National Atmospheric Deposition Program (NADP) monitoring station at Wye Mills MD just west of the Choptank basin (nadp.slh.wisc.edu/NADP, station MD013, Fig. 1), and there is a USGS gauging station near Greensboro MD (01491000, 294 km2 in the upper Choptank, hereafter “USGS Greensboro station or watershed”), with a 50-year history of discharge and water quality (Fig. 1; Fisher et al. 1998, 2006a, 2010). Since 2003, we have monitored hydrology and baseflow water chemistry at 15 small watersheds (8–51 km2) dominated by agriculture and forest (Fig. 1; Table 1). These monitored watersheds represent 33% of the area of the Choptank watershed. We have also monitored a small (~ 2 km2), 100% forested watershed (Marshy Hope) in the nearby Nanticoke basin (Fig. 1). The forested watershed is part of the Chesapeake Forest Land managed by MD DNR and is selectively logged at ~ 30-year intervals. Monitoring details can be found in Fisher et al. (2006a, 2010) and Koskelo et al. (2012, 2018).

Table 1 Summary of Choptank watershed properties, TN and TP trends, average ± se of observed TN and TP concentrations, and N and P export coefficients for the watershed and for agricultural lands

There are three estuarine stations monitored since 1984 by the MD Dept. of Natural Resources (Fig. 1). ET5.1 is a tidal fresh station primarily influenced by diffuse agricultural inputs from the watershed: 60% of the watershed’s agricultural lands are upstream of ET5.1, but only 13% of the watershed’s wastewater discharges are upstream (Fig. 1; Table 2). ET5.2 is a mesohaline station near Cambridge MD located < 2 km downstream from the largest waste water plant in the basin. Station EE2.1 lies near the connection with Chesapeake Bay and is primarily influenced by tidal exchange with bay waters.

Table 2 Discharge values for the 11 municipal wastewater treatment plants in the Choptank basin in reference year 2017 after installation of tertiary treatment (Fig. 1)

Land use in the Choptank basin is representative of the Delmarva landscape. It is dominated by agriculture and forest (~ 60% and 35%, respectively; Fisher et al. 2006b, 2010), while the USGS Greensboro watershed is 49% agricultural and 46% forested. The primary agricultural products are corn, wheat, and soybeans for poultry feed (Staver and Brinsfield 2001). The watershed contains numerous concentrated animal feeding operations (CAFO, poultry, dairy, beef), which produce manure that is applied as organic fertilizer on crop fields (McCarty et al. 2008). Grain crops and animal production facilities dominate the agricultural land, while forested areas are largely restricted to poorly drained lands where agriculture is unprofitable or housing is impractical (Fisher et al. 2006b). Although the Choptank is primarily agricultural and rural, there are 11 towns serviced by sewage treatment facilities (Table 2). Municipal wastewater treatment plants (WWTPs) in urban areas represent 90–100% of all point source N and P inputs to the Choptank (CAST 2017). Because many urban areas in the Choptank were historically established as easily accessible ports (Cambridge, Easton, Oxford, Secretary), most of the WWTPs discharge into the deeper, stratified portions of the estuary (Fig. 1; Lee et al. 2001).

Basin inputs of N and P to the Choptank estuary in the 1990s were estimated by Lee et al. (2001). They reported that atmospheric deposition represented 6% of N and 2% of P inputs, that agriculture was the dominant source (90% N, 49% P), and that wastewater was responsible for 4% of N and 49% of P inputs. Here, we compare current inputs of N and P from these three sources to those measured in the 1990s.

Methods

Data Sources

We have integrated multiple data sources to evaluate long-term changes in inputs of N and P from atmospheric deposition, agriculture, and wastewater treatment plants for comparison with changes in estuarine water quality. The data utilized are either publicly available datasets or our own measurements of hydrology, N, and P from small watersheds. We have combined data from these sources to estimate inputs to the Choptank estuary for reference year 2017 for comparison with data reported by Lee et al. (2001) for reference year 1998. The N and P inputs for years 2000–2018 reported here were then used to estimate changes in inputs since 1998 for comparison with trends in estuarine water quality observed in monitoring programs between the years 1998 and 2017.

Atmospheric Deposition

Inputs of atmospheric inorganic N were accessed for the local NADP station (Fig. 1). To estimate atmospheric N deposition to the estuary for reference year 1998, we have not used the atmospheric TN input data reported by Lee et al. (2001). Those values were based on ~ 1.5 years of TN samples, and there is considerable interannual variation in N deposition. To avoid using a single year of N data for the reference years 1998 and 2017, we used the values predicted for 1998 and 2017 by the regression through all available data (Fig. 2). NADP does not report atmospheric P or organic N deposition. Atmospheric TP inputs for reference year 1998 were estimated from the 1.5 years of TP measurements reported by Lee et al. (2001), and to our knowledge, no additional data on atmospheric P deposition in the Choptank are available.

Fig. 2
figure 2

Trends in annual precipitation (PPT, mm year−1), volume-weighted mean DIN (NH4+NO3, mg N L−1), and atmospheric N deposition (kg N ha−1 year−1) in precipitation at NADP station MD013 just outside of the Choptank basin (Fig. 1). Although precipitation was increasing at + 0.9% year−1, DIN was decreasing at − 1.4% year−1, and atmospheric N deposition decreased at − 0.9% year−1

Monitored Watersheds

Continuous hydrology data are available from the non-tidal freshwater USGS Greensboro station since 1949 and base flow and storm flow hydrology is available from our 16 monitored watersheds since 2003 (Fig. 1; Fisher et al. 2006a, 2010; Koskelo et al. 2012, 2018; Fox et al. 2014, n.d.). Baseflow here is defined as the discharge in streams from inflowing groundwater after 3 days with no rain, and stormflow is defined as the additional discharge associated with water from precipitation events (direct rainfall, overland flow, shallow throughflow). The separation of base and storm flows in the Choptank basin is described in detail in Koskelo et al. (2012, 2018).

N and P concentration data are also available for the USGS Greensboro station and at our monitored watersheds (Fig. 1). USGS data were collected 1–2 times month−1 and accessed at the USGS Website (https://waterdata.usgs.gov/nwis/dv/?site_no=01491000&agency_cd=USGS&amp; referred_module=sw). At our forested and agricultural watersheds (Table 2), we have measured nutrient concentrations in baseflow monthly since 2003 to the present. Trends in concentration data were primarily evaluated using linear regression; however, some watersheds exhibited changing trends during 2003–2018, and we used a quadratic or other fit to test for such changes. In these cases, the most recent linear trend over a minimum of 5 years in the data is reported in Table 1. To evaluate the importance of these trends in the water chemistry of monitored areas equitably, the water quality trends were area weighted to account for differences in watershed land areas (2–294 km2).

During 2013 to 2019, we collected 6–10 volume-weighted mean nutrient concentrations in stormflows at selected watersheds annually. However, we do not yet have sufficient data to capture long-term trends. We have therefore estimated annual N and P export from our watersheds using baseflow data corrected by + 14% for TN and + 64% for TP according to relationships reported in Koskelo et al. (2018) for the Choptank. The 16-year baseflow data records from these watersheds corrected for storm flows enabled us to evaluate trends in N and P concentrations and fluxes from agricultural and forested areas within the Choptank Basin. “Methods” details on these data are provided in Fisher et al. (2006a, b, 2010), Koskelo et al. (2012, 2018), and Fox et al. (2014, n.d.).

Integration of Terrestrial Diffuse Sources

We combined our direct measurements and those of USGS to compute water yields (m year−1) and nutrient yield coefficients (kg N or P ha−1 year−1) for forested and agricultural land in the 17 monitored watersheds. Monthly baseflow measurements of TN and TP concentrations (mmol m−3) were aggregated by calendar year and combined with long-term mean water yields at the USGS gauging station (0.4 m year−1 = 0.4 m3 m−2 year−1) to compute average watershed N and P yields (Y, kg ha−1 year−1) after unit adjustments. We used the 100% forested reference watershed (Marshy Hope, Table 2) to estimate N and P export coefficients for forest land (Yfor, kg ha−1 year−1) within each monitored watershed. We used this simplified approach because forest and agriculture in the monitored watersheds represent 91–100% (average = 96%) of the areas of each monitored watershed. The forest contribution to the watershed export coefficient of each watershed was inserted into a simple mass balance:

$$ {Y}_{\mathrm{w}}\ast {A}_{\mathrm{w}}={Y}_{\mathrm{ag}}\ast {A}_{\mathrm{ag}}+{Y}_{\mathrm{for}}\ast {A}_{\mathrm{for}} $$
(1)

where Yw, Yag, and Yfor are the watershed, agricultural land, and forested land yield coefficients for N or P (kg ha−1 year−1) and Aw, Aag, and Afor are the areas of the watershed, agricultural land, and forested land (ha). Solving for the agricultural contribution (Yag) to Yw, we obtained:

$$ {Y}_{\mathrm{ag}}=\frac{\left({Y}_{\mathrm{w}}\ast {A}_{\mathrm{w}}-{Y}_{\mathrm{for}}\ast {A}_{\mathrm{for}}\right)}{A_{\mathrm{ag}}} $$
(2)

The forest export coefficients of 2.1 kg N and 0.15 kg P ha−1 year−1 (Table 1) compare well with those reported in a review of the biogeochemistry of relatively undisturbed forests in the USA (0.9 kg N and 0.1 kg P ha−1 year−1) by Clark et al. (2000); however, our forest export coefficients are slightly elevated, probably due to periodic harvest activity similar to forested areas within our agricultural watersheds.

The measured Yfor and average Yag coefficients of monitored areas (33% of land; Table 1) were used to extrapolate over unmonitored areas within the Choptank basin (67% of land). Using the monitored watersheds as samples of the unmonitored areas, terrestrial inputs from unmonitored land areas were obtained from land use information in unmonitored areas (Fisher et al. 2006b) combined with Yfor and Yag in monitored areas. Inputs of directly measured watersheds and estimated inputs from unmonitored areas were then combined as total terrestrial inputs.

Wastewater treatment plants

Discharge and chemistry data were accessed for the 11 wastewater treatment plants with NPDES permits within the Choptank Basin (see Fig. 1). We used the database of the Chesapeake Area Scenario Tool (CAST 2017) to obtain measurements of annual export of N and P reported by each facility. Using the urban collection areas of each wastewater plant from Lee et al. (2001) and the CAST wastewater data, we calculated N and P yield coefficients for wastewater inputs as we have for forested and agricultural lands in our watersheds (Table 3).

Table 3 Summary of inputs to Choptank Estuary

Estuarine water quality stations

There is a 34-year history of estuarine water quality testing by MD Dept. Natural Resources and the EPA Bay Program. Data on Secchi depth and N, P, and chlorophyll a concentrations at the three monitoring stations in the Choptank (Fig. 1) were obtained from the Chesapeake Information Management System (CIMS, https://www.chesapeakebay.net/what/downloads /cbp_water_quality_database_1984_present) to examine long-term trends in estuarine water quality.

GIS, Graphics, and Statistics

ET5.1 and ET5.2 station watershed boundaries were based on the Choptank Basin boundary and 12 digit watersheds (https://data.imap.maryland.gov/datasets/maryland-watersheds-12-digit-watersheds). Portions of the watershed boundaries were digitized along the midpoints between stream reaches using ArcMap v12.5. Curve fitting and statistical tests for water quality trends were done within SigmaPlot v12.5, and parametric or non-parametric tests were used as appropriate. The use of a linear vs. quadratic or other models for concentration trends was determined based on the Akaike Information Criterion (AIC, Akaike 1973):

$$ \mathrm{AIC}=n\ast \ln \left(\mathrm{SSerror}/n\right)+2\ast K $$
(3)

where n is sample size, SSerror is sum of squares error, and K is number of model parameters + 1. The model with p < 0.05 and the lowest AIC score was chosen. If models had similar AIC values (± 7), we chose the more complex model (Burnham and Anderson 2002).

Results

Nutrient Inputs

Atmospheric Deposition

Regional atmospheric deposition of dissolved inorganic N (DIN = NO3 + NO2 + NH4+) on the estuarine surface was the second largest N source in the Choptank (Table 3). DIN data from the nearby NADP station for 1983–2017 (Fig. 2) showed a long-term, volume-weighted mean ± se of 48 ± 2 μM (0.67 ± 0.03 mg N L−1), with a linear decrease in annual average DIN concentrations (r2 = 0.59, p < 0.001, n = 35) in precipitation of − 0.89 μmol L−1 year−1 (− 0.012 mg N L−1 year−1, equivalent to − 1.4% year−1; Fig. 2, lower panel). Over the same time period, annual precipitation depth averaged 1093 ± 29 mm year−1 and increased by ~ 9 mm year−1 (0.9% year−1, r2 = 0.24, p = 0.003, n = 35), partially cancelling the effects of decreasing DIN concentrations. The DIN flux from wet precipitation averaged 7.2 ± 0.2 kg N ha−1 year−1 over the period of record, and despite the opposing trends in concentrations and precipitation, the DIN flux decreased by − 0.075 kg N ha−1 year−2 (− 0.9% year−1, r2 = 0.35, p < 0.001, n = 35, Fig. 2, upper panel). Revised inputs to the estuary from atmospheric deposition were predicted from the regression in Fig. 2 to estimate deposition in both reference years. This resulted in N deposition of 7.4 and 6.0 kg N ha−1 year−1 for 1998 and 2017, respectively, a 20% decrease consistent with the entire NADP dataset (Fig. 2).

Terrestrial Diffuse Sources

Agriculture was the largest N and P source in the Choptank (Table 3). At the USGS Greensboro station (Fig. 1), annually averaged N and P concentrations steadily increased during 1964–2017 (TN: r2 = 0.39, p < 0.001, n = 37; TP: r2 = 0.61, p < 0.001, n = 48; Fig. 3). The 15 smaller agricultural watersheds (Table 1; Fig. 1) exhibited higher concentrations, but less consistent patterns in terms of N and P trends from 2003–2018. For baseflow TP, no agricultural watersheds exhibited increases in concentrations, 11 had no changes, and two decreased. For baseflow TN, one agricultural watershed showed increases in TN, 4 showed no changes, and 10 showed decreases (Table 1). Figure 4 shows the extremes of the baseflow concentration trends observed at the outlets of the small agricultural watersheds. At South Forge watershed there was a 2.8% year−1 decrease in TN after 2008 (Fig. 4, left). In contrast, TN at Spring watershed steadily increased by 1.2% year−1 during 2003–2018 (Fig. 4, right). Despite the decreases in TN draining from 10 of the 15 agricultural watersheds, the concentrations remained high (Figs. 3 and 4; Table 1).

Fig. 3
figure 3

Interannual trends in N and P concentrations at the USGS near Greensboro in the upper Choptank basin (Fig. 1)

Fig. 4
figure 4

Trends in monthly N and P concentrations at two of the agriculturally dominated watersheds in the upper Choptank basin (Fig. 1). Both South Forge and Spring were experimental watersheds with enhanced BMPs, and these examples illustrate the range of agricultural watershed responses to BMP additions

Area-weighted trends of all monitored watersheds were used to evaluate changes in water quality in the entire monitored area of the Choptank. This revealed that 53 and 55% of the monitored area showed increasing concentrations of N and P, respectively; 9 and 41% showed no changes in N and P, and 38 and 4% showed decreases in N and P concentrations. Overall, we observed increasing or stable concentrations of TN and TP in 62 and 96% of the monitored agricultural watersheds, respectively, providing little evidence of widespread improvements in water quality.

Yield coefficients for agricultural land in 2017 were estimated for each monitored watershed (Table 1; Eq. 2). Yag estimated for the individual watersheds varied over 13–55 kg N and 0.25–0.96 kg P ha−1 year−1, averaging 33 ± 4 kg N and 0.49 ± 0.05 kg P ha−1 year−1. These values for reference year 2017 (Table 1) represented an increase compared to the reference year 1998 compiled by Lee et al. (2001), who reported average export coefficients for agriculture of 20 kg N and 0.21 kg P ha−1 year−1.

Total watershed inputs of terrestrial diffuse sources to the Choptank estuary from agriculture and forests are shown in Table 3. There are no data available on urban diffuse sources, but this land use is only 6% of the land area of the Choptank basin. Terrestrial diffuse sources of N and P from forested and agricultural lands increased between the two reference years partially due to inclusion of storm flow inputs, which were not available in 1998, and partially due to observed changes in water quality for reference year 2017. Terrestrial N inputs (forest + agriculture) to the estuary increased by 65% from 2230 to 3680 Mg N year−1, primarily due to agricultural contributions (86–93% of the total). Terrestrial P inputs increased by > 300% from 16.8 to 61.9 Mg P year−1, also primarily due to agriculture (39–77%) and our inclusion of stormflow P in the estimates.

Terrestrial Point Sources

There have been large changes in the treatment of human sewage in the Choptank basin. In the 1980s, all WWTPs had at most secondary treatment, with large and variable discharges augmented by combined storm and wastewater lines (Fig. 5, top panel). Discharge volumes increased during the 1980s and 1990s and then slowly declined after 2000 due to installation of separate storm drains in the 1990s. Export of N and P from the wastewater plants was also large and variable during 1984–2001, but most WWTPs in the Choptank were upgraded to tertiary treatment after 2000 using wastewater taxes imposed by the state of MD (Chesapeake Bay Restoration Act or “flush tax”). The advanced tertiary treatment installed during 2000–2010 and declining discharge volumes resulted in exponentially decreasing discharges of N (78%) and P (91%) for all plants by 2017, including the Cambridge WWTP (Fig. 5, middle and lower panels).

Fig. 5
figure 5

Trends in discharges from all Choptank WWTPs and the Cambridge WWTP. Data source: CAST, Chesapeake Assessment Scenario Tool; BNR, biological nutrient reduction; ENR, enhanced nutrient reduction

Inputs are a function of both intensity (yield coefficient, kg ha−1 year−1) and area (ha). Normalizing point sources to their collection areas enables comparisons with the yield coefficients of other inputs. During 1988–2000, the data of Lee et al. (2001) resulted in wastewater N yield coefficients (24 kg N ha−1 year−1) somewhat higher than the average for agricultural lands (20 kg N ha−1 year−1); however, the wastewater P coefficients during 1988–2000 (6.9 kg P ha−1 year−1) were more than an order of magnitude higher than agricultural diffuse sources (0.21 kg P ha−1 year−1). In contrast, the data reported here for 2001–2018 after advanced tertiary treatment show greatly reduced coefficients of 5.7 kg N ha−1 year−1 (smaller than atmospheric N deposition) and 1.7 kg P ha−1 year−1, lower than in the 1990s, but still the highest P yield coefficient.

Summary of Inputs

Total estimated N and P inputs to the Choptank estuary increased between the two time periods (Table 3). During 1988–2000, atmospheric and terrestrial N inputs to the estuary were 2483 Mg N year−1, 88% of which were contributed by agriculture. The N inputs increased to 3872 Mg N year−1 in 2001–2017, a 53% increase, which was primarily due to increasing N from terrestrial diffuse sources (Table 1). Agriculture was responsible for 92% of the inputs calculated for reference year 2017, and the remaining inputs were < 5% each. Unfortunately, the decreases in atmospheric N deposition (− 19% from the 1998 estimates) and wastewater N discharges (− 76%) were overwhelmed by the increases in agricultural N.

P inputs to the Choptank estuary also increased by 20% from the earlier to later time period estimates (Table 3). In the earlier time period, P inputs were 58 Mg P year−1, primarily from wastewater plants (28 Mg year−1) and agriculture (23 Mg year−1). However, during 2001–2018, there were 91% decreases in wastewater P inputs, making current wastewater P the smallest input for which data are available. There are currently no measurements that we are aware of for atmospheric P deposition in the Choptank, but it is likely that atmospheric P input to the estuary is much smaller than current wastewater inputs, given the values measured by Lee et al. (2001; Table 3). Agricultural inputs increased from 1988 to 2000 (23 Mg P year−1) to 2001–2018 (53 Mg P year−1). Overall, estimated P inputs to the estuary have increased from 58 to 69 Mg P year−1 for the two time periods.

The Choptank Estuary

Estuarine station ET5.1 is the most upstream sampling location in the estuary (Fig. 1). ET5.1 had the highest average chlorophyll a concentrations, most shallow Secchi depth, and highest average N and P concentrations of any river sampling location, including the USGS Greensboro station (Table 4; Fig. 6). Annual average chlorophyll a concentrations were high during 1984–1998 (20–30 μg L−1), and then steadily declined during 1998–2005 to 5–12 μg L−1, largely remaining at these concentrations through 2018 (Fig. 6). Summer bottom dissolved O2 was close to atmospheric equilibrium (~ 7 mg O2 L−1) in 1984, and then steadily declined during most of the monitoring period to under-saturated values of 4–6 mg O2 L−1. Secchi depths were shallow (~ 0.3 m) and remained unchanged during the entire period.

Table 4 Summary of mean station properties in the Choptank River and estuary for 1984–2018
Fig. 6
figure 6

Trends in annually averaged surface chlorophyll a, July–Sept bottom-dissolved oxygen and annual average Secchi depth at the tidal fresh stations ET5.1 near the turbidity maximum in the Choptank estuary. Data source: MD DNR

Mesohaline estuarine station ET5.2 lies further downstream in the estuary near Cambridge MD (Fig. 1). During 1984–2005, this station exhibited degrading water quality, with increases in annual average chlorophyll a in surface waters and decreases in summer bottom dissolved O2 and Secchi depths (Fisher et al. 2006a, 2010; Fig. 7). However, after 2005 chlorophyll a concentrations stabilized and then slightly decreased. Furthermore, O2 concentrations in bottom waters started to increase, and Secchi depth increased. These recent changes in water quality occurred as our estimates of current N and P inputs to the estuary were increasing (Table 3) but corresponded closely in time with reductions in wastewater N and P inputs (Fig. 5), much of which occurred at the Cambridge wastewater treatment plant < 2 km upstream from estuarine station ET5.2.

Fig. 7
figure 7

Trends in annually averaged surface chlorophyll a, July–Sept bottom-dissolved oxygen and annual average Secchi depth at the inner mesohaline station ET5.2 near the chlorophyll maximum in the Choptank estuary. Data source: MD DNR

Station EE2.1 is the outermost mesohaline estuarine station in the Choptank near the main stem of Chesapeake Bay (Fig. 1). At EE2.1, chlorophyll a steadily increased from 5–10 to 10–20 μg L−1 during 1984–2018 (Fig. 8), gradually converging with chlorophyll a concentrations at station ET5.2 (Fig. 7). Summer bottom dissolved O2 was under-saturated at EE2.1 (4–7 mg O2 L−1), with no trend during 1984–2018, and was similar to the improving values at ET5.2 after 2015. Secchi depths at EE2.1 decreased from 1.5–2.0 m in the 1980s to 1.0–1.5 m in 2015–2018, in parallel with increasing concentrations of chlorophyll a.

Fig. 8
figure 8

Trends in annually averaged surface chlorophyll a, July–Sept bottom-dissolved oxygen and annual average Secchi depth at the outer mesohaline station EE2.1 in the Choptank estuary near the main axis of Chesapeake Bay. Data source: MD DNR

Discussion

Spatial and Temporal Patterns

We have shown that there were changes in nutrient inputs to the Choptank estuary between reference years 1998 and 2017. Atmospheric N deposition to the estuarine surface was the second largest N input to the estuary (Table 3) but decreased at − 0.075 kg N ha−1 year−1, or approximately − 0.9% N year−1, for a total reduction of 20% between 1998 and 2017 (Fig. 2). Decreases in atmospheric N deposition via precipitation are commonly observed in the Eastern USA (e.g., Eshleman et al. 2013; Lloret and Valiela 2016; Detenbeck et al. 2019) as a result of the Clean Air Act (Sullivan et al. 2018). In Europe, similar rates of decline in atmospheric N deposition have occurred due to legislation, but the current European deposition rates are much higher (10–50 kg N ha−1 year−1; Dirnbӧck et al. 2018) than in the USA.

Terrestrial diffuse sources, primarily agriculture, were the largest inputs of both N and P to the Choptank estuary. Unlike atmospheric deposition, terrestrial diffuse inputs of N and P increased between 1998 and 2017. The increases in agricultural N and P export coefficients are consistent with the increasing concentration trends at the USGS Greensboro station (Figs. 3 and 4) and also support the corrections based on Koskelo et al. (2018) that we applied for storm flows. Storm contributions to the estimates of N and P inputs in the earlier time period reported by Lee et al. (2001) were not available, and our inclusion of storm effects partially contributes to the differences between the two reference years. However, the changes reported in Table 3 (63 and 134% for N and P, respectively) are larger than the additions from stormflows (14% and 64%). There was large spatial variability in the trends within the monitored areas of the Choptank, and our smaller agricultural watersheds integrate less area (8–51 km2, Table 1) and illustrate variability at spatial scales smaller than the USGS Greensboro station (294 km2).

Increasing concentrations of N and P in agricultural areas have been commonly reported. For example, similar increases in baseflow N have been observed in the Mississippi River, despite application of agricultural conservation efforts, and reductions in dead zones in the Gulf of Mexico may take decades due to large stores of legacy N and P in soils and groundwater (van Meter et al. 2018; Stackpoole et al. 2019). Similarly, the River Seine in France has a history of increases in agricultural N inputs due to agricultural intensification since the 1950s (Billen et al. 2001; Romero et al. 2016).

Tertiary treatment of urban point sources in the Choptank resulted in large decreases in N (− 78%) and P (− 91%) in WWTP discharges to mesohaline areas (Fig. 5; Table 2). These large changes appear to have influenced water quality at the ET5.2 estuarine monitoring station. Prior to 2000, annual average surface chlorophyll a was increasing and summer bottom dissolved O2 was decreasing, clear symptoms of increasing eutrophication. As wastewater inputs declined due to upgrades to tertiary treatment after 2000 (Fig. 5), the water quality of the Choptank estuary at station ET5.2 improved (Fig. 7). This station is < 2 km downstream from the Cambridge MD WWTP, the largest in the Choptank. Improvements in estuarine water quality occurred at ET5.2 despite overall increases in N and P inputs to the estuary (Table 3).

Unlike estuarine monitoring station ET5.2, ET5.1, and EE2.1 showed decreasing water quality. Water quality at ET5.1, upstream of ET5.2 and primarily upstream of the wastewater inputs, showed a mixed pattern of water quality responses (decreasing chlorophyll a, no change in Secchi depth, and decreasing summer bottom dissolved oxygen to under-saturated concentrations; Fig. 6). Since water depth at ET5.1 is ~ 10 m and vertically unstratified (Berndt 1999), the shallow Secchi depth of 0.3 m (indicative of poor water transparency) reflects high turbidity and low light availability in most of the water column. The declining chlorophyll a and summer bottom dissolved O2 indicate increasing heterotrophy at this station. At the most saline estuarine station EE2.1, downstream of ET5.2 and strongly influenced by exchange with Chesapeake Bay, chlorophyll a increased as Secchi depth decreased, but with no trends in summer bottom O2, which were close to atmospheric equilibrium (Fig. 8). The data from this station is consistent with our input budgets showing increasing N and P moving from the Choptank watershed into the estuary (Table 3).

These three stations (Table 4; Figs. 6, 7 and 8) portray gradients of water quality in the Choptank estuary. At the freshwater end there are high chlorophyll a and nutrients indicative of local watershed nutrient sources (Table 4). Moving into the mesohaline stations, N and P concentrations are diluted by water from Chesapeake Bay with lower concentrations of N and P, and only station ET5.2, in close proximity to the Cambridge WWTP, shows signs of improving water quality. The nutrient gradients of N and P across these stations are indicative of net movements of N and P downstream through the Choptank estuary into the Chesapeake Bay.

Reductions in atmospheric deposition and wastewater discharge occurred a few years prior to improvements in estuarine water quality at ET5.2 (Figs. 2, 5, and 7). This sequence of events and the proximity of ET5.2 to the Cambridge WWTP (Table 2), suggest that reductions in atmospheric N deposition and wastewater N and P discharges are responsible for the improving estuarine water quality at ET5.2. This indicates that the investments in mitigating N emissions to the atmosphere and decreasing N and P discharges by upgrading to tertiary treatment at local wastewater plants were responsible for the improvements in estuarine water quality at ET5.2 (Fig. 7). However, we cannot eliminate the influence of augmented estuarine filter feeders such as zooplankton and oyster populations (Bricker et al. 2020, https://www.chesapeakeprogress.com/abundant-life/oysters) or other restoration efforts that occurred in the estuary. Some quantitative estimates on removal of N and P by filter feeders is available (e.g., Newell et al. 2004), but additional research would be useful. Regardless of the cause of the improving water quality at ET5.2, Millette et al. (2019) have proposed using rates of increases in chlorophyll a, such as the increases during 1985–2000 at ET5.2 (Fig. 7), to target areas for nutrient reduction, and the subsequent improvements in water quality which we observed at ET5.2 support this.

The proximity of the Cambridge WWTP to ET5.2 makes this station particularly sensitive to changes in inputs from municipal wastewater. It is likely that the changes at ET5.2 are a local response to reduced WWTP inputs. Furthermore, unmeasured sources downstream of ET5.2 or occasional Chesapeake Bay inputs (Sanford and Boicourt 1990) could be responsible for increasing chlorophyll a and decreasing Secchi at EE2.1. Most agricultural inputs are upstream of EE2.1, along with most of the drainage basin, but two small wastewater treatment plants are close to EE2.1 (Fig. 1). However, the major gradient of N and P concentrations is from freshwater to the Bay (Table 4). We do not have enough information to understand the differences between the interannual variations in water quality at stations EE2.1 and ET5.2, and this is clearly an important area for future research in the Choptank.

Similar effects of reductions in WWTP N and P have been observed throughout other eutrophic estuarine and coastal areas. Examples include the Seine River and estuary (France, Billen et al. 2001), Roskilde Fjord (Denmark, Staehr et al. 2017), Tampa Bay (USA, Johansson and Lewis III 1992), Gunston Cove (Potomac River USA, Jones 2020), the Bilbao estuary (Spain, Garcia-Barcina et al. 2006), and the Avon-Heathcote estuary (New Zealand, Barr et al. 2020; Zeldis et al. 2020). However, in the Stockholm inner archipelago of the Baltic Sea, reductions in WWTP P inputs were buffered by internal recycling from sediments (Walve et al. 2018). In a review by Boesch (2019), quick improvements in estuarine water quality usually resulted from mandated wastewater treatment plant upgrades, with smaller or no reductions from terrestrial diffuse inputs to eutrophic coastal areas worldwide.

The spatial and temporal characteristics of the three estuarine stations in the Choptank are consistent with the conceptual model of estuarine biogeochemical structure proposed by Fisher et al. (1988). Station ET5.1 lies in or near the turbidity maximum of the Choptank estuary where fresh and salt water first mix. The station is ~ 10 m deep, with no salinity stratification (Berndt 1999), quite turbid (indicated by the shallow Secchi depth in Fig. 6), with N and P concentrations higher than those at the USGS Greensboro station (Fig. 1; Table 4). This indicates additional sources of N and P in the watershed between the USGS Greensboro station and ET5.1, and we suggest that the source is likely agricultural areas between the two monitoring stations. Although declining chlorophyll a at ET5.1 can be viewed as indicative of improving water quality, it is more likely due to light limitation of phytoplankton growth under turbid, nutrient-rich conditions (Table 4; Fisher et al. 1999).

Estuarine station ET5.2 has much lower concentrations of N and P compared to tidal fresh station ET5.1 (Table 4). Between these stations, the Choptank River is sinuous with large areas of oligohaline wetlands (Fig. 1), and we suggest that these wetlands are an important landscape sink within the Choptank basin for upstream N, P, and sediment. Sixty percent of the agricultural diffuse sources are upstream of these wetlands (Table 4), which have been shown to have high rates of denitrification and burial of N and P both in the Choptank (Malone et al. 2003) and the Patuxent tributaries of Chesapeake Bay (Jordan et al. 2003; Boynton et al. 2008). This feature of the geomorphology of the estuarine tributaries of Chesapeake Bay may capture a significant fraction of terrestrial inputs from upstream agricultural sources. We have no direct quantification of the potential effects of the oligohaline wetlands lying between ET5.1 and ET5.2, and we suggest that this is a fruitful area for future research efforts in the Choptank and elsewhere.

In contrast to ET5.1, station ET5.2 lies near the chlorophyll maximum of the Choptank estuary (Berndt 1999). The chlorophyll maximum is a mesohaline region of estuaries where salinity stratification limits the depth of mixing and turbidity decreases, enabling phytoplankton in surface waters to access sufficient light to accumulate biomass using the abundant nutrient supplies from upstream or local sources (Fisher et al. 1988). While the location and size of the turbidity and chlorophyll maxima can shift up or downstream with river discharge (Fisher et al. 1988; Berndt 1999), ET5.2 appears to be a sensitive, well-placed station to indicate the collective effects of upstream, non-tidal, nutrient inputs, apparent modulation of those inputs by oligohaline marshes, and local, more downstream WWTP inputs.

An interesting question is why there is improving water quality at ET5.2 despite an overall increase in N and P inputs to the estuary as a whole. This response suggests that local actions matter; in this case greatly reducing local inputs from the largest WWTP improved adjacent estuarine water quality, even when the overall estuary was receiving more N and P. It is possible that other local water quality improvements occurred unrecorded near other waste water plants.

There are also other interpretations of the results reported here. The increasing chlorophyll a and decreasing Secchi depth at downstream station EE2.1 may have been due to unmeasured nutrient sources downstream of ET5.2 or tidal exchange with Chesapeake Bay. However, the major gradient of N and P concentrations is from freshwater to the Bay indicative of the importance of terrestrial nutrient sources (Table 4). Most agricultural inputs are upstream of EE2.1, along with most of the drainage basin, but two small wastewater treatment plants are close to EE2.1 (Trappe and Oxford; Table 2; Fig. 1). Although Trappe is a small plant, it only had secondary treatment during this period of record and is the third largest source of WWTP N. We do not have enough information to understand the differences between the interannual variations in water quality at stations EE2.1 and ET5.2, and this is clearly an important area for future research.

Test of the Hypothesis

The data presented provide partial support for our hypothesis that management efforts to reduce inputs from the three known nutrient sources in the Choptank watershed (agriculture, wastewater, atmospheric deposition) have improved water quality in the Choptank estuary. Both atmospheric deposition and WWTP inputs declined due to management actions, whereas overall inputs increased due to higher agricultural inputs, despite conservation efforts. An important question to explore is why efforts to reduce agricultural diffuse sources have not provided reductions in N and P inputs to the Choptank estuary. Many agricultural BMPs have been applied throughout the Choptank watershed (Staver and Brinsfield 2001), and many BMPs have been found within agricultural watersheds that we monitor (Fox et al. n.d.). There is a similar report by Kroon et al. (2016) on the lack of effectiveness of agricultural conservation actions in areas contributing to declining water quality near the Great Barrier Reef in Australia.

There are several potential explanations for the lack of an agricultural contribution to the improving water quality in the Choptank estuary. First, municipal WWTP are single entities funded by tax payers and subject to government regulations (NPDES) that are assessed through self-monitored data submitted to the US EPA (e.g., Fig. 5). In contrast, farmers are private business owners with a range of attitudes concerning voluntary conservation practices (Smith et al. 2007; van Grieken et al. 2019; Fox et al. n.d.). There is no routine monitoring of groundwater or surface discharges from farms with BMPs except for research projects such as this one or non-profit water quality organizations. This situation provides little incentive or feedback to farmers on voluntary BMP applications. Since at least some costs of the BMPs fall on the individual farmer, some hesitate to apply BMPs or apply additional BMPs because of the time and expense (Fox et al. n.d.). Cary and Roberts (2011) noted similar reservations by Australian farmers to adopt BMPs.

Second, there are time lags of years in the delivery of groundwater to streams (e.g., Meals et al. 2010; Sanford and Pope 2013). For example, median groundwater ages of 5–10 years have been reported in Choptank watersheds (Sutton et al. 2009). However, groundwater nitrate concentrations can be altered significantly in short time periods. Bunnell-Young et al. (2017) showed reductions in groundwater nitrate from ~ 700 μM (9.8 mg N L−1) to ~ 70 μM (0.98 mg N L−1) in the surface unconfined aquifer over 2–3 years after fields with grain production (corn-wheat-soy) were converted to unfertilized conservation lands in the Choptank Basin. However, this lower nitrate groundwater must be transported to a stream before we detect the impacts of positive agricultural management. It is possible that the small decreases in the baseflow TN of the agricultural watersheds reported in Table 1 are due to the slow emergence of groundwater with lower nitrate concentrations that are appearing as a result of conservation efforts that occurred 5–10 years ago. However, our widespread measurements of groundwater nitrate under agricultural fields with BMPs do not show systematic decreases in N concentrations (Table 1; Fox et al. 2014, n.d.).

Last, insufficient density, extent, and poor placement of voluntary conservation practices might be responsible for the lack of agricultural contribution to improving Choptank water quality. In active agricultural fields, current nitrate levels in groundwater are high (300–1200 μM, 4–17 mg NO3-N L−1, Fox et al. 2014, n.d.), similar to or higher than in streams draining mixed land uses (Figs. 3 and 4). In addition, agricultural soils typically have high P levels greatly exceeding crop needs (100–400 mg P kg−1; Sims et al. 2002; Fisher et al. 2018). These legacy nutrients will continue to deliver large quantities of N and P to downstream receiving waters such as the Choptank estuary until soil and groundwater N and P concentrations are reduced by more technical assistance and outreach to farmers on BMPs (Lemke et al. 2010). One possible solution to reduce the high concentrations in waters draining agricultural areas is a more systematic application of BMPs with follow-up monitoring of conservation effects to validate the expectations of BMP impacts (Lemke et al. 2011).

There is strong economic pressure on farmers to improve crop yields, regardless of low N and P use efficiency (Cassman et al. 2002). This economic pressure could be countered by equivalent economic pressure for good performance of conservation efforts (Winsten 2009; Winsten and Hunter 2011; van Grieken et al. 2019). Market or other financial rewards for validated conservation management could potentially result in more dense applications of BMPs and significant reductions in N and P losses from agricultural areas (Liu et al. 2018). This approach will require involvement of all stakeholders, including the broad range of food and feed producers, in order to balance the economic viability of agriculture and its impact on the environment.

Improvements in water quality have occurred in the middle of the Choptank estuary (Figs. 1 and 7). This appears to be due primarily to public tax investments and government regulations to reduce atmospheric deposition (Clean Air Act) and WWTP discharges (Chesapeake Bay Restoration Act). However, we have not yet detected evidence of systematic reductions in agricultural N and P inputs, although we did measure some reductions in N and P from some smaller, agricultural watersheds (Fig. 4; Table 1). If reductions in agricultural N and P inputs over broad areas do occur in the future, improvements in estuarine water quality larger than those reported here, and more consistently throughout the entire estuary, could be expected. For this reason, it is important to continue monitoring agricultural areas with enhanced management practices. Direct measurements are needed to evaluate whether reduced concentrations of N and P are occurring in groundwater and soils of areas with BMPs or whether inefficient or insufficient management practices have been implemented to make significant reductions in losses of N and P from agricultural fields.

Conclusions

The Choptank basin is 1756 km2 of land draining to 280 km2 of an estuarine tributary of the Chesapeake Bay. There are strong similarities between the Choptank basin and the Chesapeake as a whole, which enables the Choptank to be used as a model for progress in the Chesapeake. We updated an N and P input budget for the Choptank estuary previously constructed for reference year 1998 to evaluate progress towards water quality goals in reference year 2017. We found decreases in atmospheric N deposition and human wastewater inputs, but overall increases in agricultural N and P inputs, the largest contributor. Despite the lack of widespread reductions in agricultural N and P, we found local improvements at some small agricultural watersheds and at an estuarine station. The estuarine station is near a large WWTP outfall, and improvements in water quality occurred following reductions in wastewater and atmospheric N and P inputs from public investments in clean air and WWTP upgrades. While these local effects are encouraging, it is important to continue to improve agricultural practices to reduce losses of N and P from farms. Even greater improvements in estuarine water quality throughout the estuary could be expected following reductions from this largest source of N and P.