Agricultural Best Management Practice Abundance and Location does not Influence Stream Ecosystem Function or Water Quality in the Summer Season

Best management practices (BMPs) are tools commonly used to mitigate negative impacts of agriculture on water quality; however, the relationship between BMPs and aquatic ecological function is unknown. Our research goal was to determine the association between both stream ecosystem metabolism and water quality, and the abundance and location of four different BMPs in agricultural catchments. Dissolved oxygen was measured over a two-week period in mid-June and used to estimate ecosystem metabolism of 13 headwater streams representing a gradient of BMP implementation in Southern Ontario, Canada. Stepwise regression analyses were used to associate stream metabolism and water quality with metrics describing the abundance and location of BMPs within each catchment. Studied streams exhibited rates of metabolism comparable to catchments from other agricultural regions. However, metrics of BMP implementation were not associated with either stream metabolism or water quality. Our results suggest that BMPs in the studied agricultural catchments are not improving water quality or mitigating water quality impacts on stream metabolism during the summer season. We propose that seasonality of catchment hydrology and time lag effects associated with past agricultural land use may be masking the mitigation benefits of BMPs on stream ecosystem conditions during the summer season.


Introduction
Streams draining catchments dominated by agricultural land uses are often exposed to a variety of non-point source contaminants known to degrade water quality and alter key stream ecosystem processes [1].Agricultural practices can release nutrients and sediments to streams in agricultural areas through use of fertilizers, soil tillage, livestock husbandry and grazing [2][3][4].Nutrients from fertilizers and animal waste (i.e., nitrogen and phosphorus) are essential limiting resources for the growth of aquatic vegetation and elevated inputs from agricultural sources have been reported to increase primary productivity in streams [5].Excess primary productivity may cause periodic hypoxic conditions at night due to plant respiration, as well as from heterotrophic respiration following the die off of plant biomass, stressing ecological communities [6].However, agricultural practices also have the potential to limit primary production through increased sediment loading and associated reductions in water clarity limiting light penetration to benthic habitats [7].Growing evidence of the impacts of agriculture on the health of receiving streams and associated downstream ecosystems has made the need to implement effective management strategies to mitigate agricultural effects a major issue in agricultural regions worldwide [8].
In response to the growing recognition of the economic and environmental consequences of nutrient and sediment release from agricultural lands, managerial and structural conservation strategies, or Best Management Practices (BMPs), have been developed to help mitigate the effects of these common stressors.Structural BMPs are designed to reduce loss of pollutants to aquatic ecosystems by enhancing filtration, cementation, and uptake of pollutants in runoff on, or at the edge, of agricultural lands [9].Common structural BMPs include planting and maintenance of riparian buffer strips, livestock exclusion fences, manure storage facilities, and erosion control structures.Additionally, structural BMPs can be strategically implemented in locations identified as significant contributors of surface water pollutants to help maximize mitigation effects [10][11][12][13].
Many past studies have examined the effectiveness of BMPs in mitigating the impacts of agriculture at small scales and noted improved quality of agricultural runoff [14][15][16].However, past studies exploring the relationship between BMPs and stream ecological condition have shown more equivocal results, despite also frequently reporting reduced concentrations of nutrients and sediments associated with BMP use [14,15,[17][18][19][20][21].To date, past studies have largely focused on taxonomic based measures of community condition, particularly aquatic invertebrates.However, potential for rehabilitation of ecological diversity in agricultural streams may be limited by historical depletion of the regional species pool frequently associated with widespread human activity [22][23][24].As such, ecological benefits associated with BMP use may be most apparent in rehabilitation of ecological processes, such as nutrient dynamics and organic matter input, that are more directly linked to the nutrient and sediment input rates that BMPs are designed to manage [25].There has, however, been limited examination of the response of river ecosystem processes, such as stream metabolism, to the implementation of BMPs.
Stream metabolism is the combination of primary productivity and ecosystem respiration and is thus a measure of the balance between production and consumption of organic matter in a stream.As such, stream metabolism is a critical ecological process providing insight into ecosystem biomass and trophic dynamics [26,27].Combined, gross primary productivity (GPP) and ecosystem respiration (ER) can be used to determine the Net Ecosystem Metabolism (NEM; NEM = GPP ´ER) and the production to respiration ratio (P/R; P/R = GPP/ER) of the stream.These parameters indicate the trophic status (i.e., net heterotrophic or net autotrophic) of a stream providing an estimate of the food base and the flow of energy at a catchment scale [27].Stream metabolism has been widely recommended as a biological indicator to assess the functional status of riverine environments [27][28][29].Furthermore, the response of stream metabolism to key stressors commonly associated with agricultural activities (i.e., temperature, sediment, nutrient concentration, and organic matter) has been extensively analyzed and reviewed [27].Stream metabolism is thus very well suited to serve as an indicator of the response of stream processes to mitigation effects of agricultural BMPs.
The goal of this study was to examine the mitigation effects of agricultural BMPs on water quality and stream metabolism in the summer season.This goal was achieved by assessing variation in stream water quality and ecosystem metabolism among agricultural headwater streams and associating these ecosystem metrics with the number and location of four common structural BMPs (i.e., manure storage facilities, livestock access restriction fences, erosion control structures, and riparian buffers) present in each catchment.It was predicted that headwater catchments containing a larger number of BMPs and BMPs located within areas of greater connection to hydrologic flow routes would show improved ecosystem condition, as indicated by decreasing concentrations of nutrients and sediments, as well as rates of stream metabolism.

Study Area
Our study was situated within the Grand River Watershed (GRW), which is located in Southwestern Ontario, Canada (Figure 1a).The watershed covers an area of approximately 7000 km 2  Water 2015, 7, 6861-6876 containing the major urban centers of Kitchener, Waterloo, Guelph, Cambridge, and Brantford.Settlement in the GRW occurred in the early 1800s along with the onset of deforestation and agricultural practices.Agricultural practices intensified from the 1960s onward with the introduction of synthetic fertilizers and pesticides.Contemporary land use patterns in the GRW continue to be dominated by agriculture, covering approximately 80% of the landscape, with fragmented, natural vegetation covering about 19% of the total area [30].Major agriculture types in the watershed include a combination of cash crops (i.e., soybeans and corn) and livestock operations (i.e., dairy cattle, beef, hog, and poultry) [31].Climate in the GRW is temperate, with the central area (43.45 ˝N, 80.38 ˝W) having an average annual temperature, precipitation, and relative humidity of 7.0 ˝C, 916.5 mm, and 87.8%, respectively [32].
Water 2015, 7page-page 3 dominated by agriculture, covering approximately 80% of the landscape, with fragmented, natural vegetation covering about 19% of the total area [30].Major agriculture types in the watershed include a combination of cash crops (i.e., soybeans and corn) and livestock operations (i.e., dairy cattle, beef, hog, and poultry) [31].Climate in the GRW is temperate, with the central area (43.45°N, 80.38° W) having an average annual temperature, precipitation, and relative humidity of 7.0 °C, 916.5 mm, and 87.8%, respectively [32].

Site Selection and BMP Description
For our study, 13 sampling sites were established on small agricultural streams within the headwaters of the Nith and Conestoga subwatersheds (Figure 1b,c).Both subwatersheds are characterized by an agricultural landscape (~83% agriculture), till soils (~76% of the soils), and a flat to rolling topography [33].Study sites were selected to incorporate regional variation in BMP employment and comparable physiography and land cover among the catchments.Selected catchments were limited to those with areas between 5 and 12 km 2 , till based soils (>65%), low natural (<25%) and urban (<5%) land covers, high agricultural (>85%) land use, locations accessible for sampling, and those with streams unimpeded by dams or other structures.The streams selected for sampling have not been subject to routine chemical or ecological monitoring in the past and thus historical or pre-BMP data was not available.
This study focused on four structural BMPs: (1) manure storage systems (MS); (2) livestock access restrictions (LAR); (3) erosion control structures (EC); and, (4) riparian buffer areas (RBA) (Table 1).Structural BMPs have the advantage of being able to be strategically located on the landscape in hydrologically connected areas to maximize mitigation benefits.Furthermore, these BMPs are readily identifiable on the landscape using orthoimages.BMPs in the study area were predominantly implemented through collaboration between landowners and the Rural Water Quality Program administered by the Grand River Conservation Authority (GRCA).BMPs in this database were located in the study catchments using a Geographic Information System (GIS).Additional BMPs implemented by individual farmers were identified using high-resolution (0.5 m) orthoimages [34].Sampling sites encompassed the full range of BMP abundance and location of the four different BMPs present in the region.A set of BMP metrics were calculated for each of the 13 selected catchments based on the type and location of BMPs present in the catchment (Table 2).Abundance descriptors included # of BMPs per farm, % RBA, length of LAR, % river with LAR, and % river protected.Location was characterized by hydrologic distance to stream, hydrologic distance to sampling point, flow accumulation, and % flow accumulation.The calculated BMP metrics were placed into one of two groups, BMP summary and BMP type.The BMP summary group consisted of five metrics that were calculated as a summation of all BMPs (i.e., MS, LAR, and EC) present in the catchment, excluding riparian area.In the BMP type group there was a total of 28 metrics describing the abundances and locations of the individual BMP types (i.e., MS, LAR and EC) within the catchments.The extent of RBAs was also included in the BMP Water 2015, 7, 6861-6876 type group as a metric describing the vegetated proportion of the total riparian area.Location based metrics were not calculated for RBAs due to difficulties in differentiating the stream channel from the riparian zone given the small stream channel widths (<2 m) and level of resolution available for land cover and DEM layers.EC metrics were not individually examined because of low abundance (n = 2); however, EC BMPs were included in the BMP summary group metrics.

Field Sampling
Daily rates of stream metabolism were estimated at our 13 study sites from instream dissolved oxygen (DO) measurements collected over a two-week period (10 June to 26 June) in summer of 2014.An early summer sampling period was chosen in order to encompass a time period where plants in both aquatic and terrestrial systems were expected to be actively growing [35].Concentrations of DO were recorded at 15-minute intervals over the two-week sampling period using a D-Opto logger (Zebra-Tech Ltd).Stream metabolism was estimated for each study reach using the open-system, single station method [36].In this method a single DO logger is used to measure DO changes over time at a specific point along the study reach.The DO change is assumed to be uniform throughout the upstream reach.Therefore, potential drivers occurring along this reach are presumed to influence metabolism measurements.Probe deployment followed the recommendations of Grace and Imberger [36].Logger malfunction at one site resulted in this site being excluded from analyses leaving 12 sites included in the study.
To determine whole stream metabolism of each study site, the nighttime regression method was used to calculate the reaeration coefficient (k) as recommended by Young et al. [37].This regression was performed using the River Metabolism Estimator v1.2 spreadsheet macro in Microsoft Excel v14.4.6 [38] developed by Young and Collier [39].Regressions were considered reliable if the calculated k coefficient was positive, and the value of R 2 was greater than 0.4 [sensu 39].Abiding to these constraints, the calculation of the k coefficient for one site was modified due to the occurrence of weak R 2 values in the nighttime regressions throughout the sampling period.As the recorded data appeared reliable, the k coefficient of this site was determined using the O'Connor-Dobbins empirical equation as recommended by Young et al. [37] for streams with relatively low reaeration coefficients (i.e., <50 day ´1).Calculated k coefficients were used to determine daily values of gross primary production (GPP) and ecosystem respiration (ER) based on the following general equation: where ∆DO is the change in DO concentration (g O 2 m ´3) over the 15 min interval, GPP is volumetric , k is the oxygen reaeration coefficient, and D is the DO deficit (g O 2 m ´3) over the measurement interval.Rates of GPP and ER were then corrected and standardized into areal daily rates using equations 45 through 49 from Grace and Imberger [36].Two other metabolism parameters were used to describe the relationship between GPP and ER, the NEM and the P/R ratio.NEM was calculated as the difference between the values of GPP and ER and the P/R ratio was calculated as the quotient of GPP over ER.Daily values of the calculated metabolism parameters (i.e., GPP, ER, NEM, and P/R) from all reliable sampling days were then averaged by site to give an estimation of metabolism rates for the two-week sampling period.Physicochemical characteristics of each stream were assessed and recorded at the time and location of logger deployment.Measurements collected were nutrients, pH, conductivity, flow velocity, channel depth, channel width, discharge, temperature, and canopy cover.Instream nutrients were measured twice at each of the 13 study reaches by collecting grab water samples at the time of logger deployment and retrieval.Samples were collected according to the Canadian Aquatic Biomonitoring Network (CABIN) protocol [40] and analyzed for nitrogen and phosphorous species (i.e., total nitrogen (TN), ammonium (NH 4 ), nitrite and nitrate (NO 2 -NO 3 ), total phosphorus (TP), total dissolved phosphorus (TDP), and soluble reactive phosphorus (SRP)) using a flow injection analysis [41].Water samples were also analyzed for turbidity and the amount of total suspended solids (TSS) following the nephelometric method and a standard dry weighting method, respectively [41].pH and conductivity were measured in a well-mixed area of each reach using a Hoskin Scientific Professional Plus Probe (Model: Pro 10102030).Flow velocity was measured at a series of evenly spaced points (0.1 m or 0.2 m increments channel width pending) across a bank to bank transect perpendicular to the thalweg with a current velocity meter (Swoffer Instruments Inc. Model: 2100).Measurements were averaged to provide an overall estimate of velocity for each sampled reach.The channel cross-sectional depth was measured using a similar approach.Stream depths were measured from 10 evenly spaced points across five transects spaced 10 m apart upstream of the logger.The depths from all measured points along each reach were then averaged to provide an estimate of the average stream depth along the study reach.Stream water temperature at each site was recorded every 15-min alongside the DO measurements by the D-Opto logger.Temperature measurements were averaged over the entire sampling period as descriptors of each site.The immediate riparian canopy coverage of each site at the location of the D-opto logger was measured using a densiometer (Forestry Suppliers Inc. No. 43888).Canopy measurements were taken in the middle of the stream channel at waist height at the most downstream, middle and most upstream transect locations.Measurements of the proportion of canopy cover were taken facing the upstream, downstream, right bank, and left bank directions.All 12 canopy measures were averaged for the site as an overall measure of reach canopy cover.

Data Analysis
Stream metabolism may be affected by a variety of external environmental variables that can mask the potential influence of BMPs.In small, low flow streams proximate variation in light Water 2015, 7, 6861-6876 availability is known to be an important factor that can alter the local metabolism rates [42].Thus, in order to isolate the potential impact of BMPs on stream metabolism rates for among site differences in channel shading of the sampling reach, the effect of light availability was removed by calculating the residuals of the relationship between the canopy cover metric, and both GPP and ER stream metabolism metrics using a least squared regression analysis.Residuals from these regressions were outputted as the new values for each metabolism parameter for use in further statistical analyses.These regressions were completed using SYSTAT 13 v.13.00.05 [43].
Multiple least squared regression analysis was used to determine the strength and direction of the association between the non-collinear BMP metrics and the stream metabolism residuals, as well as water quality parameters, for the sampled catchments.Separate regression analyses were performed on each residual metabolism metric and water quality species with each of the BMP metric groupings (i.e., BMP summary and BMP type) serving as independent predictive variables in their respective analyses.For the analysis, a backwards-stepwise regression method was used (confidence interval = 0.95; probability = 0.10).Regression models were considered significant if p < 0.10.This method was also used to determine the association between the stream metabolism residuals and the non-collinear water quality parameters.To reduce the likelihood of TYPE I error associated with stepwise regressions we followed recommendations by Blanchet et al. [44] and first tested the statistical significance of a global model and ensured that the resultant R 2 from the stepwise process was not larger than that of the global model.Non-significant global models were not analyzed further using stepwise analyses.Regression analyses were performed in SYSTAT 13 v.13.00.05 [43].

Results
The majority of measured physical parameters exhibited minimal variation among the study catchments (i.e., CV < 0.5; Table 3).Sampled headwater streams were similar in regards to average channel depth (mean = 0.14 ˘0.06 m), average channel width (mean = 1.65 ˘0.56 m), water temperature (mean = 17.96 ˘1.42 ˝C), pH (mean = 7.87 ˘0.27) and conductivity (mean = 575.84˘87.74 µS/cm).However, flow velocity, discharge and canopy cover revealed greater variation among sites (i.e., large = CV > 0.5).Flow velocity ranged from <0.01 m/s to 0.13 m/s, with a mean velocity of 0.05 m/s.Average discharge was 0.02 m 3 /s and ranged from <0.01 m 3 /s to 0.05 m 3 /s.Canopy cover ranged from a minimum of 6.92% to a maximum of 85.75% cover.Mean canopy coverage was 29.27%.Water quality parameters (i.e., TN, NH 4 , NO 2 -NO 3 , TP, SRP, TDP, turbidity, and TSS) sampled during logger deployment and retrieval generally exhibited large variation among the study catchments.The exception was TN concentrations collected at logger retrieval (CV = 0.33).Furthermore, concentrations of all examined nutrient species and TSS were at least 1.3 fold different between sampling events with samples taken during logger retrieval having a consistently larger mean, median, maximum, and standard deviation among sites.
A total of 60 BMP structures were identified in our 13 study catchments.This total consisted of 40 MS structures, 18 LAR fences, and 2 EC structures (Table 4).The total number of BMPs per study catchment ranged from 1 to 16 with an average BMP implementation of 4.62 structures per catchment.On average each catchment contained 3.08 MS structures, 1.38 LAR fences, and 0.15 EC structures.Average BMP use was less than one BMP per farm (mean = 0.41).BMP implementation was on average 0.29, 0.09, and 0.02 BMPs per farm for MS, LAR, and EC, respectively.Abundance metrics of LAR and EC structures exhibited larger variation among the study sites, as these BMPs were less common and not found in all catchments.Of the 13 study catchments selected, LAR fences were present in six and EC structures in two.BMP location showed substantial variation in regards to the flow accumulation area and hydrologic distance from the BMP to the river in both the BMP summary and BMP type groups, with the exception of LAR flow distance to the river (CV = 0.46).In contrast, location metrics involving the hydrologic distance from the BMP to the instream sampling point exhibited minimal variation among the catchments for all BMP metrics except EC (CV = 0.76).The length of LAR adjacent to the stream was variable among catchments and ranged from 63 m to 2846 m (CV = 1.03).Likewise, the percentage of river protected by LAR structures exhibited a large variation among catchments (CV = 1.65).The percentage of RBA surrounding the study streams ranged from 0% to 42% among the catchments with an average of 20%.Combining the proportion of river protected through LAR and RBA, the minimum percentage of protected channel length was 5% and the maximum was 61%.

VIF Analysis
Of the five metrics in the BMP summary group, none were found to be collinear (i.e., BMPs/farm, BMP % flow accumulation, BMP mean flow accumulation, BMP flow distance to sampling point, and BMP flow distance to river).In the BMP type grouping 19 of 28 metrics were found to be collinear, leaving only nine remaining metrics for further analysis, which included MS/farm, MS median flow accumulation, MS % flow accumulation, MS flow distance to sampling point, MS flow distance to river, LAR/farm, LAR median flow accumulation, LAR % flow accumulation, and % RBA.Among the water quality parameters, all phosphorus species were collinear leaving six of eight independent variables for further analysis including TN, NH 4 , NO 2 -NO 3 , TP, turbidity, and TSS.

Stream Ecosystem Metabolism
Estimated k coefficients were moderately variable among catchments (CV = 0.53) and had a mean value of 21.99 day ´1 (Table 5).The calculated magnitudes of mean daily GPP and ER exhibited more variation than did the k coefficients and both varied approximately 40-fold across the study sites.GPP rates ranged from 0.88 to 35.93 g O 2 /m 2 /day with a mean of 11.75 g O 2 /m 2 /day for all selected sites.ER occupied a comparable range from 0.73 to 26.78 g O 2 /m 2 /day and had an average of 11.11 g O 2 /m 2 /day.Sites were typically production dominated (GPP > ER) with an average NEM of 0.64 and P/R ratio of 1.48.However, negative NEM values and P/R values of less than one were observed at 4 of the 12 study catchments.
GPP and ER were negatively associated with the amount of riparian canopy cover at each sampling site.GPP exhibited the strongest association with canopy cover having an R 2 value of 0.342 (n = 12, p = 0.027), whereas ER had an R 2 value of 0.169 (n = 12, p = 0.102).Based on the significance of the relationships between canopy cover, and GPP, residuals were used in subsequent regressions concerning GPP.Despite the lack of statistical significance between canopy cover and ER, we considered there to be sufficient variance in ER explained by light availability to warrant using residuals in further regressions regarding ER.
Table 5. Summary of the reaeration coefficient (k), gross primary productivity (GPP), ecosystem respiration (ER), net ecosystem metabolism (NEM), and production to respiration ratio (P/R) of 12 headwater streams in the Grand River Watershed.Global regression analyses between stream metabolism residuals (i.e., GPP and ER) and BMP summary metrics showed that BMP abundance and location did not explain a significant amount of variation in stream metabolism among the study catchments.Furthermore, global models between stream metabolism residuals and BMP type metrics also showed no significant associations.Similarly, global regression analyses between water quality parameters (i.e., TN, NH 4 , NO 2 -NO 3 , TP, SRP, TDP, turbidity, and TSS) and both BMP summary and type metrics revealed non-significant results.Lastly, the results of global regression models between stream metabolism residuals and non-collinear water quality parameters also revealed no significant associations.

Discussion
Ecosystem metabolism of streams draining agricultural land is known to be elevated due to the impact of cultivation related stressors, such as increased light and nutrient inputs [45].The average metabolic rates of the 12 streams in our study exhibited a range of 0.88-35.93g O 2 /m 2 /day for daily GPP and 0.73-26.78g O 2 /m 2 /day for daily ER.Metabolism values from the studied agricultural catchments (i.e., 78%-97% cover) are comparable to values reported by previous studies in agricultural streams in Central USA (e.g., Wiley et al. [45]: GPP = 0.1-44.2g O 2 /m 2 /day, ER = 6.2-41.6 g O 2 /m 2 /day; Bernot et al., 2010 [46]: GPP = 0.1-16.2g O 2 /m 2 /day, ER = 0.9-15.7 g O 2 /m 2 /day; Frankforter et al. [47]: GPP = 0-8.4g O 2 /m 2 /day, ER = 0.03-36 g O 2 /m 2 /day).Although the maximum rates of GPP and ER from our study sites are almost two-fold higher than rates from Bernot et al. [46] and Frankforter et al. [47], this difference may be attributed to the shallow (mean = 0.14 m) and relatively slow moving (mean = 0.05 m/s) streams assessed in our study.Shallow depths are associated with increased light availability and associated primary production at the streambed in all areas of the stream [48].Likewise, low flow rates can allow for the accumulation of nutrients and organic matter in the sediments and water column, facilitating uptake by biota and increasing instream production and respiration [49].Our study also examined metabolism over a two-week period in summer.Temporal differences associated with seasonality have been seen to alter the magnitude of stream metabolism, with greatest rates generally occurring from spring through autumn [50].However, the greater light intensity and duration, and resultant increase in water temperature in the summer months can promote a greater activity of both photosynthetic and heterotrophic organisms facilitating higher metabolic rates.
Concentrations of major nutrients sampled in our study were variable among catchments, and frequently exceeded both national and provincial water quality guidelines.The average TP concentration of all samples was approximately 9-fold greater than the provincial threshold for the protection of aquatic life (0.03 mg/L) [51] and reflective of a hyper-eutrophic status (>0.1 mg/L) [52].Likewise, the average NH 4 concentration was approximately 10-fold greater than the proposed national guideline (0.019 mg/L) [53].However, the average nutrient concentrations of all collected samples (TN = 2.69 mg/L, NH 4 = 0.185 mg/L, NO 2 -NO 3 = 1.66 mg/L, TP = 0.27 mg/L, SRP = 0.18 mg/L, TDP = 0.19 mg/L) were comparable to results of studies on other small agricultural streams in the North American Great Lakes region.For example, Gabel et al. [21] reported concentrations of 0.087, 0.179, and 1.101 mg/L for TDP, NH 4 , and NO 3 respectively, among seven agricultural streams with abundant BMPs.Additionally, Chambers et al. [54] found an average TN concentration of 2.68 mg/L and TP concentration of 0.08 mg/L along a gradient of agricultural catchments in Southern Ontario (6.7%-100% cover).Likewise, sediment descriptors (Turbidity = 7.66 NTU, TSS = 8.84 mg/L) in our study were also similar to measures reported in streams from other agricultural catchments containing a variety of BMPs (e.g., Yates et al. [18]: Turbidity = 8.10 NTU, TSS = 21.3 mg/L).Although differences are noted between phosphorus and TSS parameters, this may be attributed to the type of agricultural practice, and BMP characteristics within the catchments.
Our study found BMP abundance and location metrics to be unassociated with rates of stream metabolism.The lack of association is not surprising in light of the fact that BMPs were not associated with summer concentrations of major nutrient forms and that nutrient concentrations were not related to stream metabolism.Our finding that both stream metabolism residuals and water quality parameters are not associated with BMP metrics may be attributable to one, or a combination, of three potential hypotheses: (1) BMP implementation has not exceeded the threshold required to initiate an ecosystem response; (2) limited connection between BMPs and the hydrologic system during the summer season; and (3) legacy effects of past agricultural activities and practices.
Stream ecosystems often display a threshold response to a variety of environmental factors [1,55], including management practices [18].Although our study catchments comprised a substantial range of BMP implementation, a significant change in stream ecological condition was not associated with the increases in abundance or improved location of BMPs attained in the Grand River Watershed.As such it may be that the combined effects of BMPs in our study area are insufficient to reduce nutrient losses to streams and initiate a measurable change in water quality and stream metabolism.This explanation is supported by predicted nutrient loss reductions from models testing BMP implementation strategies in Great Lakes river basins [56].In this past study, a scenario of BMP implementation levels considered feasible through conservation efforts similar to the Rural Water Quality Program were insufficient to generate substantive reductions in nutrient losses.However, a scenario of more extensive and more targeted BMP implementation was found to generate significant reductions in nutrient yield.Based on these findings the Rural Water Quality Program may need to adopt strategies that encourage more intensive implementation and improved targeting of BMPs toward hydrologically connected areas in order to enhance the condition of streams in the Grand River watershed.
BMPs assessed in past studies have been shown to reduce nutrient levels in agricultural streams (e.g., [21,[57][58][59][60][61][62]).However, results from our study did not suggest that BMPs were enhancing water quality in the studied streams.The absence of associations between BMPs and water nutrients may be attributed to seasonal climate variation [63,64].Greater evapotranspiration rates in summer months can remove soil water and dry out the land resulting in precipitation infiltrating the soil in all but the largest of events.Under these warm, dry conditions the land is less likely to generate significant surface runoff and BMPs designed to capture nutrients and sediments in runoff have less influence on surface water quality.Furthermore, the availability of nutrients on the landscape is abated during the summer season from actively growing crops and reduced agricultural inputs.In contrast, during spring and autumn agricultural inputs intensify and soils are more often saturated due to lower evapotranspiration rates and large amounts of precipitation or melt water.Saturated conditions will result in greater overland connection between the land and hydrologic system where structural BMPs may more effectively reduce transport of contaminants.Indeed, Holmes [65], found significant negative associations between instream nutrient concentrations and BMP abundance and location metrics based on samples collected in November from similar sized catchments in the GRW.
As our study did not find any of the calculated BMP metrics to be significant predictors of nutrient concentrations, it is possible that instream nutrient concentrations are being driven by processes other than surface runoff (e.g., groundwater flow) during the drier, summer season.Additional studies are needed to assess the seasonal effectiveness of BMPs in mitigating agricultural stressors known to impact ecosystem function.
A third potential hypothesis for the lack of association between current BMP use and a response in ecosystem condition is the potential for agricultural stressors, in this case nutrients, to persist through legacy effects.Many studies have outlined the importance of legacy effects in ecological assessments of stream ecosystems [66][67][68][69].Alterations to the ecological condition of stream ecosystems as a result of agriculture have been observed to influence the current physical, chemical, and biological characteristics, including the whole-stream GPP, even after 55-years of agriculture abandonment and succession [68].In regards to BMPs implemented in active agricultural catchments such as the ones in our study, Meals et al. [70] proposed that the lack of response in water quality and affiliated indicator parameters is likely a result of lag time.Lag time realizes the persistence and legacy of agricultural stressors and encompasses the time required for the BMP to produce a beneficial effect, for that effect to be delivered to the aquatic system, and for the aquatic system to respond to the effect [70].Depending on the aim of a BMP, this lag time can vary from years to several decades differing with the BMP treatment and mitigation target, spatial scale, hydrologic flow pattern, flushing rate, and intensity and residence time of the stressor [70].For example, long-term effects of nitrogen fertilizer use in small agricultural watersheds has been proposed to last greater than 30 years if leeched into soil and ground water pools, which can be a particularly important non-point source input into streams during dry seasons [71].Likewise, the presence of phosphorus in agricultural soils can be evident for approximately 14 years [72,73].As the BMPs implemented through the Rural Water Quality Program were instituted from 3 to 15 years prior to sampling it is possible that not enough time has passed for the beneficial effects of the BMPs to impact stream ecosystem function.Although BMPs are advocated as effective conservation strategies [74], the abundance and persistence of agricultural stressors may be overwhelming and delaying the effectiveness of BMPs in improving the function of agricultural streams in the GRW.As BMPs in our study are still relatively young, it is recommended that continued BMP implementation and long term monitoring be employed as a strategy to achieve and realize ecosystem improvements.

Figure 1 .
Figure 1.Map showing the location of the Grand River Watershed within the Great Lakes region of North America (a); and locations of the 13 headwater study catchments (grey) within the Nith and Conestoga subwatersheds of the Grand River Watershed (b).Map showing examples of the distributions of analyzed BMPs (Yellow: manure storage structure, Purple: livestock access restriction fence, and Red: erosion control structure) within study catchments (c).Orange surface areas represent agricultural land use.

Figure 1 .
Figure 1.Map showing the location of the Grand River Watershed within the Great Lakes region of North America (a); and locations of the 13 headwater study catchments (grey) within the Nith and Conestoga subwatersheds of the Grand River Watershed (b).Map showing examples of the distributions of analyzed BMPs (Yellow: manure storage structure, Purple: livestock access restriction fence, and Red: erosion control structure) within study catchments (c).Orange surface areas represent agricultural land use.

Table 1 .
Descriptions and expected benefits of four structural best management practices included in this study.

Table 2 .
BMP metrics and associated descriptions used to quantify abundance and location of agricultural BMPs in small catchments in southern Ontario.

Table 3 .
Descriptive statistics for physico-chemical parameters collected from 13 study streams.Water quality parameters were measured during logger deployment (shaded rows) and after the sampling period during logger retrieval (unshaded rows).Values of <0.01 indicate concentrations or velocities below detection limits.

Table 4 .
Descriptive statistics of BMP metrics for the 13 sampled catchments.BMP descriptions are separated based on metric grouping (i.e., BMP summary or BMP type).