Abstract
Many agricultural watersheds in the United States have impaired waterbodies due to nonpoint source pollution from agricultural activities and related processes. To understand the physical, chemical, and biological integrity of surface water in a coastal agricultural watershed, spatial and seasonal patterns of physicochemical and biological properties were investigated in Bayou Lacassine watershed (BLW) in Louisiana, United States. The relationship between the physicochemical and biological properties were also investigated. Sampling sites were located in the Bayou Chene and Lacassine Bayou subwatersheds within the BLW. Dissolved oxygen (DO), turbidity, conductivity, temperature, pH, total suspended solids (TSS), total dissolved solids (TDS), total solids (TS), five-day biological oxygen demand (BOD5), nitrate and nitrite-nitrogen (NO3/NO2-N), total Kjeldahl nitrogen (TKN), soluble reactive phosphorus (SRP), total phosphorus (TP), chloride (Cl−), fluoride (F−), and sulfate (SO4) were determined weekly from samples collected during 2012 to 2015. Fish and benthic invertebrate diversity and abundance in the two subwatersheds were determined in early summer and in fall of 2012 and 2013 at nine sites. Water quality was generally better at the most downstream site than at the most upstream site where agricultural intensity was highest, with significant differences in turbidity, TSS, TDS, TS, NO3/NO2-N, TKN, TP, and BOD5. There was also seasonal variation for the water quality parameters due to variability in agricultural activities and climatic conditions within the watershed. Results of the relationship between physicochemical properties and fish community variables showed that species richness, diversity, and abundance were negatively affected by elevated TS, NO3/NO2-N, and conductivity. For the benthic invertebrates, diversity was negatively related to BOD5. This study demonstrated unexpected longitudinal and seasonal patterns in physicochemical and biological properties of surface waters in a coastal agricultural watershed. This information is valuable in developing nonpoint source pollution control strategies for these subwatersheds.
- Key words
- agriculture
- fish and macroinvertebrates
- nonpoint source pollution
- physicochemical and biological properties
- water quality
- watershed
Introduction
Many agricultural watersheds in the United States have impaired waterbodies due to nonpoint source pollution from agricultural activities and related processes, such as soil erosion, fertilizer and pesticide application, agricultural drainage, and surface runoff (Lombardo et al. 2000; Poudel et al. 2010; Perez and Walker 2014). Water quality impairment can also be due to inputs from household septic systems (Poudel 2016) and natural sources such as wildlife, mineral deposits, and algal blooms. The impairment of waterbodies in agricultural watersheds can exhibit strong seasonality, corresponding with cropping seasons and related agricultural activities (Poudel et al. 2013). Similarly, agricultural watersheds can manifest spatial variation in water quality impairment due to differences in land use types, soils, nutrient inputs, and agricultural activities across the landscape (Demcheck et al. 2004; Mueller-Warrant et al. 2012; Poudel et al. 2013). The impairment of physical and chemical properties of surface waters can negatively impact biological communities due to, for example, hypoxia and harmful algal blooms (Zhou et al. 2008; Broussard and Turner 2009; Riseng et al. 2011; Budria 2017; Breitburg et al. 2018), increased levels of fecal bacteria (Brendel and Soupir 2017), elevated levels of suspended sediment (Basnyat et al. 1999; Riseng et al. 2011), or the presence of pesticides (Echeverría-Sáenz et al. 2012; Anderson et al. 2014, 2018).
The implementation of best management practices (BMPs) in agricultural watersheds has lowered nutrient loadings and improved biotic index scores (Maret et al. 2008), increased water clarity (Kronvang et al. 2005), and increased seagrass abundance (Greening et al. 2014). While these examples include cases where effects on biota were assessed, the effectiveness of BMPs is typically assessed solely by evaluating the physicochemical properties of surface waters (Yeung et al. 2017). However, measuring biological properties of surface waters is likely to provide more conclusive information on the effectiveness of BMPs in restoring the integrity of waterbodies (Karr 1993). As an example, invertebrate community indices were responsive to the adoption of stream health BMPs at deer farms in New Zealand (Rhodes et al. 2007). These BMPs included fencing to exclude deer from the streams, water troughs, culverts, buffer strips, and grazing management. Other studies have assessed the impacts of habitat, land use, and water quality on fish and/or benthic invertebrates in coastal watersheds (Helson and Williams 2013; Sawyer et al. 2004). While both fish and benthic invertebrates respond to environmental changes (Pilière et al. 2014), it appears that these two taxonomic groups are often affected by different water quality and habitat characteristics (Backus-Freer and Pyron 2015; Johnson and Ringler 2014). This makes monitoring of fish and benthic invertebrates complementary, such that assessments and monitoring for effectiveness of BMPs are best done by considering both of these taxonomic groups.
Water quality and nonpoint source pollution control issues are different in coastal watersheds in the southern United States than they are in most of the previously studied watersheds, due to the unique combination of specific agricultural activities and environmental conditions in coastal areas. Coastal agricultural watersheds in southwestern Louisiana (Huner et al. 2002; Poudel and Jeong 2009; Poudel et al. 2013; Poudel 2016) are good examples of this uniqueness. The main agricultural products in these watersheds include rice (Oryza sativa L.), soybeans (Glycine max [L.] Merr.), and sugarcane (Saccharum officinarum L.). Crawfish production and cattle farming (pasture) are other agricultural activities. Rice planting occurs from late February to early May, and the rice growing season continues, with ratoon harvest, until September, while crawfish production occurs from October to May (Huner et al. 2002). Rice fields are drained after crop harvesting through mid-spring depending on management practices such as winter-holding of water for weed control, water leveling of rice fields, water conservation, and crawfish production in rice fields. Crawfish ponds are drained from late spring to early summer. These watersheds are characterized by the lack of a substantial elevation gradient, small and slow-moving streams referred to as “bayous” in this part of the United States, tidal influence at the coastal part of the watersheds, the presence of swamps, and frequent flooding.
The specific objectives of this study were (1) to assess the presence of spatial variation (including a longitudinal gradient) in physicochemical water quality variables and in fish and benthic invertebrate abundance and diversity, (2) to assess the presence of temporal variation (including seasonal differences) in these variables, and (3) to assess the relationships between physicochemical properties of surface waters and fish and benthic invertebrate abundance and diversity. Information generated from this study will be beneficial for watershed managers and planners evaluating water quality and/or developing watershed management plans in coastal agricultural watersheds.
Materials and Methods
Study Area. This study was conducted in a coastal agricultural watershed in southwestern Louisiana, United States—the Bayou Lacassine watershed (BLW) in the Mermentau River Basin. The BLW consists of two subwatersheds: Lacassine Bayou and Bayou Chene subwatersheds (figure 1). Their key characteristics are given in table 1. Both subwatersheds were first included on the 303(d) list of impaired water bodies in 1999. The suspected causes of impairment were excess nutrients, sediments, turbidity, lead, organic enrichment, and low dissolved oxygen (DO) (table 1) (Louisiana Department of Environmental Quality, 1999 court ordered 303[d] list) (LDEQ 2000). Both waterbodies were not meeting fish and wildlife propagation uses as required by the Clean Water Act.
According to the Bayou Lacassine Watershed Plan (LDEQ 2009), BMPs covering a total of 21,030 ha were implemented in this watershed between 2004 and 2009. Major BMPs implemented included conservation crop rotation, seasonal residue management, nutrient and pest management, and irrigation land leveling. Similarly, BMPs covering a total of 58,274 ha were implemented in the Bayou Chene subwatershed between 2005 and 2015 (LDEQ 2016). These included conservation crop rotation, irrigation land leveling, prescribed grazing, nutrient management, residue and tillage management, no-till, irrigation water management, integrated pest management, dry seeding, and pesticide management. Agencies and programs involved in the implementation of BMPs in BLW included the Office of Soil and Water Conservation, Section 319 funds, USDA Natural Resources Conservation Service, Louisiana Department of Agriculture and Forestry, local soil and water conservation districts, resource conservation and development districts, and Louisiana Master Farmer Program. It is very likely that more than one BMP was implemented on the same piece of land.
Sampling, Field Measurements, and Laboratory Analyses. Weekly water quality monitoring was done at 10 sampling locations (sites 1C to 10C) in the Bayou Chene subwatershed and five sampling sites (sites 1L to 5L) in the Lacassine Bayou subwatershed (figure 1). Sites were selected on the basis of factors including land use type, implementation of BMPs, and accessibility. Sampling duration and the number of sampling events are presented in table 2. To achieve more intensive sampling of areas with major agricultural activities, six additional sites that represented the tributaries of Bayou Chene were added in August of 2014. While water quality sampling for sites 1L, 2L, 3L, 4L, 5L, and 7C ended on March 26, 2015, sampling at the remaining sites was continued until May 21, 2015, to capture conditions during the late spring season. Site 7C was dropped due to sampling safety risks. Consequently, the variability in the number of sampling events among sites reflects differences in the sampling period rather than differences in sampling frequency. Surface water temperature, turbidity, pH, conductivity, and DO were determined in the field using a YSI Sonde attached to a handheld data logger (model 6820 with 650MDS, YSI Incorporated, Yellow Springs, Ohio). Other water quality parameters were determined from water samples collected weekly with the use of a Van Dorn sampler. Laboratory analyses, using Standard Methods (Clesceri et al. 1998), quantified total suspended solids (TSS), total dissolved solids (TDS), total solids (TS), five-day biological oxygen demand (BOD5), nitrate/nitrite-nitrogen (NO3/NO2-N), total Kjeldahl nitrogen (TKN), soluble reactive phosphorus (SRP), total phosphorus (TP), chloride (Cl−), fluoride (F−), and sulfate (SO4). During each sampling event, a single in-stream measurement and a single grab sample were collected at each site at the center of the main flow of the stream and at one-half the total water depth when the stream was less than 2 m deep, or at 1 m depth when the stream was over 2 m deep. The average water depth of sampling locations ranged from 0.27 m at site 9C, to 4.33 m at site 2L (table 2). While sites 9C, 7C, and 4C had less than 1 m average water depth, sites 3L, 8C, 3C, 5L, 1L, 10C, and 4L had average water depth between 1 to 2 m, and sites 1C, 2C, 6C, 5C, and 2L had average water depth greater than 2 m. Except for flood events, water depths did not vary much during the year. Prior to the collection of water samples, sampling equipment was rinsed with ambient water. Based on the Jennings weather station, which is located within the watershed, average annual rainfall during this study (2012 to 2015) was 1,742 mm, which exceeded the longer-term (1980 to 2011) annual rainfall average of 1,524 mm.
Biological sampling was conducted for benthic invertebrates and for fish. Both groups of organisms were collected at nine sites (1L to 5L, 1C to 4C) during four sampling periods. These sampling periods covered two seasons (summer and fall) during each of two years (2012 and 2013). Summer samples were collected during May and/or June, and fall samples during September and/or October. In order to get better representation of the fish community, fish sampling was done using both cast and dip nets. A 3.66 m radius and 9.53 mm mesh size net was used to collect fish samples. Fish from 20 net casts (5 casts in 4 different areas at a site) were combined into one sample. For each dip net sample, we sampled for 20 minutes (5 minutes at each of 4 different spots at a site). Fish were identified to the species level (except for one individual for which only genus-level identification was possible). A total of 9,390 fish were collected. These fish belonged to 26 species and 12 families; western mosquitofish (Gambusia affinis) were very abundant throughout, gizzard shad (Dorosoma cepedianum) were common at a single site, and the other species were generally encountered in low numbers (table 3). Fish that could not be identified in the field or that were kept as voucher specimens were euthanized; other fish were released upon identification and quantification.
To collect benthic invertebrates, sediment samples were taken with a Petite Ponar grab (four replicates per site during each sampling period) and sieved in the field (500 μm screen). The >500 μm material was fixed in formalin (10%) with a Rose Bengal stain. Samples were later preserved (70% ethanol with 5% glutaraldehyde), sorted (animals separated from the detritus), and identified. A total of 9,731 individual organisms were collected and identified. All individuals were identified to at least the taxonomic level of order. The benthic invertebrates collected belonged to 19 orders; the most abundant were aquatic worms, leaches, clams, flies (mainly chironomid larvae), mussels, and snails (table 4).
Data Analyses. Means and ranges were determined for all of the physicochemical variables at each of the monitoring sites. The presence of differences among sites and among seasons was determined for each of the physicochemical variables with analyses of variance (one-way ANOVA), followed (where overall effects were significant) by pairwise comparisons using Tukey-Cramer honestly significant difference (HSD) tests. For the fish and macroinvertebrate data, spatial differences were analyzed by one-way ANOVA, while differences between seasons and years were analyzed by two-way ANOVA (with the interaction factor dropped from the model when this factor was not significant). The relationships between the physicochemical and biological properties (abundance and diversity data for fish and benthic invertebrates) and the presence of temporal changes for the fish and benthic invertebrate data were determined by regression analysis. All statistical analyses were conducted in JMP (SAS Institute Inc. 2009).
Results and Discussion
Spatial Patterns. In Bayou Chene, average values for TSS, TDS, TS, NO3/NO2-N, TKN, TP, F−, SO4, and conductivity were lower at the most downstream site (1C) than at the most upstream site (4C), by 18.7% to 60.3% (table 5). The same pattern was observed for Lacassine Bayou, where average values for turbidity, TSS, TDS, TS, NO3/NO2-N, TKN, TP, BOD5, F−, and conductivity at downstream site 2L were lower by 23.0% to 50.7% compared to upstream site 4L (table 6). A similar spatial pattern, with water quality improving with distance away from the more-developed upstream areas, has been observed for enteric bacteria (Mallin et al. 2000) and BOD (Yoon et al. 2015). The present study's data are indicative of a better water quality downstream. Potential reasons for the poorer water quality upstream include the upstream presence of concentrated agricultural activities, the downstream presence of riparian buffers and wetlands that improve surface water quality by enhancing nutrient uptake, sediment retention, litter decomposition (Whigham et al. 1988; Johnston 1991), and the dilution effect due to increased volume of water downstream.
In Bayou Chene, the downstream reduction in nutrients and suspended sediment was not accompanied by an increase in DO (e.g., the DO level at the furthest downstream site was 53.4% lower than it was at our furthest upstream site). We believe that site-specific factors were responsible for the latter. Water depth was very low (table 2), and the water flow appeared to be relatively fast at upstream site 4C, providing more opportunity for air/water exchange of oxygen (O2). In contrast, DO values were low at the downstream site 1C, a site with typically slow water flow and a buildup of decomposing vegetation.
The analysis of fish community and benthic invertebrate densities also showed spatial variation and a longitudinal gradient. Statistically significant differences among sites were detected for several of the fish and benthic invertebrate density and/or diversity measures (table 7). Though trends were variable, densities and diversities were, for example, typically higher at sites 3C and 1C than they were at site 4C and higher at sites 3L and 1L than at site 4L (table 7). These results show that the biological integrity for water bodies in this watershed was generally higher downstream than it was at upstream sites—opposite from the more typical pattern of water quality being highest at upstream sites (Wynes and Wissing 1981) though in line with the physicochemical data.
Temporal Patterns. Physicochemical properties of surface water in BLW showed very strong seasonal patterns (figures 2a to 2j). Dissolved O2 levels were highest (averaging 7.7 mg L−1) in January and lowest in August (averaging 2.1 mg L−1). Monthly averages declined sharply from February to May, and then remained low until levels started to increase in October, with averages of around 5.0 mg L−1 measured for November and December. Thus, April to October is the period that DO conditions require special attention in this watershed, as DO levels decline to levels associated with hypoxia and impact invertebrate and fish communities (Baden et al. 1990; Briggs et al. 2017). Although lower DO levels coincided with higher water temperatures, a closer look at conditions during the February to May period reveals that lower DO values also coincided with high turbidity, TDS, and nutrient levels. Consequently, the lower DO levels during these months are expected to be due to a combination of variables—work on the hypoxia situation in the northern Gulf of Mexico has demonstrated that underlying mechanisms can be quite complex (Bianchi et al. 2010). While reducing water temperatures (e.g., by increasing riparian vegetation) would aid in improving DO levels in these waterbodies, minimizing inputs of sediments, nutrients, and other dissolved substances could also be very beneficial.
Various factors likely contributed to the high turbidity and high levels of solids and nutrients during the March to May period. While agricultural activities during the preceding winter season are typically limited, rainfall is quite high during the months of January and February. Consequently, this may have resulted in runoff from dry and barren fields. Winter-killed grasses in pasturelands and the livestock there may also have contributed (Poudel 2016), along with the large flocks of overwintering geese in the subwatersheds. Early spring is also a time for the drainage of crawfish ponds and rice fields. Wet planting of rice requires drainage from rice fields during these months, which contributes to high sediment loads (LSU AgCenter 2014). Another likely contributing factor is runoff from fertilizer applied to the rice fields being prepared for spring planting.
With regard to the seasonality of fish community, very high numbers of mosquitofish were caught at some sites during summer fish sampling, which appeared to have resulted in a lower overall fish community diversity during summer than fall (table 8). This diversity difference disappeared when mosquitofish were excluded from the analysis. It is well known that mosquitofish can handle low-O2 conditions (Stoffels et al. 2017), in part due to their ability to gulp air from the air/water interface. No other obvious seasonal differences existed for the fish community metrics. While it was anticipated that especially the fish abundance would have shown a seasonal effect in response to the low summer DO levels in these bayous, the fall sampling conducted during September/October may have been too early in the fall to allow the fish community to recover from the hypoxia period. It is also possible that the resident biota are adapted to low-O2 conditions. Research with both natural and laboratory populations has demonstrated that the evolution of resistance to environmental stressors can occur rapidly (Oziolor et al. 2014; Xie and Klerks 2003). The presence of an overall temporal trend in the fish community was assessed by comparing 2012 with 2013 data. No significant differences were detected between the two years (table 8). There were also no significant trends detected by regression analyses using the four time points (data not shown).
For the benthic invertebrates, both abundance and diversity generally increased during the two-year period that this sampling was conducted (figure 3). Benthic invertebrate density increased with time for both bayous combined (R2 = 0.929; p = 0.0361), but these trends were not statistically significant for either Bayou Chene (R2 = 0.210, p = 0.542) or Lacassine Bayou (R2 = 0.759, p = 0.129) separately. Biodiversity parameters showed a more consistent increase with time. For Bayou Chene, the regression of the Shannon diversity index, H’, was not quite statistically significant (R2 = 0.815, p = 0.097), while it was significant for the number of taxonomic orders (R2 = 0.916, p = 0.043). In Lacassine Bayou, the biodiversity increased over time for both the diversity index parameter (R2 = 0.904, p = 0.049) and the number of taxonomic orders (R2 = 0.991, p = 0.004). While the four sampling points are insufficient for providing solid evidence for a temporal trend, this provides some indication that water quality is improving as a consequence of introduced management practices. Evidence provided below (“Relationship between Physicochemical Characteristics and Fish and Benthic Invertebrates”) supports the notion that the organismal data agree with results for established physicochemical water quality measures. There may be a lag time to such an agreement. Research evaluating the effectiveness of forestry BMPs indicates that results for water quality variables may be more informative than results from biomonitoring when monitoring is done shortly after implementation of BMPs, since the ecological changes may have a slower response time (Yeung et al. 2017).
Relationship between Physicochemical Characteristics and Fish and Benthic Invertebrates. The biomonitoring results were correlated to those of the physicochemical monitoring. For the fish community variables, regressions were statistically significant for three of the physicochemical variables (TS, NO3/NO2-N, and conductivity) and three of the fish community metrics (number of species, H’, and abundance on basis of cast net sampling), with all the relationships indicating lower fish abundance and diversity at higher levels of TS and nutrients (table 9). For example, the biodiversity index H’ was clearly lower for sites with TS levels exceeding 300 mg L−1 (figure 4a). It is well documented that fish community diversity is typically negatively affected by eutrophication (Heiskary and Bouchard Jr. 2015), though some tolerant species may replace some of the sensitive ones. While eutrophication is more likely to result in increased biomass, certain effects on fish abundance appears to depend on specific habitat conditions (Artigas et al. 2013). High levels of suspended sediment or high turbidity are known to negatively affect fish diversity and abundance (Richardson and Jowett 2002).
For the benthic community data, relationships with the physicochemical data were statistically significant only for the biodiversity (as expressed by the number of taxonomic orders represented in the samples) and the BOD5 (table 10). It is not surprising that there was a significant effect for this biodiversity variable; taxon richness measures are particularly responsive to environmental stressors (Carlisle and Clements 1999). The diversity was inversely related to the BOD5 (figure 4b). A high BOD5 value is likely to result in low DO levels in the water, and anoxic sediments typically have a benthic community dominated by the few groups (notably oli-gochaetes and tolerant chironomid taxa) that can handle these conditions. It is noticeable that the relationships between biomonitoring and physicochemical properties were more pronounced for the fish community than the benthic community. This may be a function of the physicochemical monitoring's focus on the water column rather than the sediment, though it may also be a consequence of the lower detail in benthic invertebrate identification (identification to order level rather than the species-level identification done for fish). Benthic invertebrates generally respond more directly to water quality measures than is the case for fish communities (Johnson and Ringler 2014).
Summary and Conclusions
In the present study's coastal watershed, water quality was better at the downstream sites than at the upstream sites. This unusual longitudinal pattern was observed for both the physicochemical variables measured and the fish and benthic invertebrates. Water quality appeared inversely related to land use intensity and the prevalence of agricultural activities in the watershed and positively related to the presence of wetlands and marshes. Temporal differences were also observed in this watershed. Strong seasonal patterns were observed for the physicochemical variables, with water quality being lowest in spring and summer, due to a combination of anthropogenic factors (chiefly agricultural activities) and climatic factors such as temperature and rainfall. Abundance and biodiversity data for the benthic invertebrates provided some indication that water quality is improving. Observed relationships between physicochemical variables and biological properties confirmed that the physicochemical variables were reflective of ecological integrity of the watershed. This study confirms the need for considering both spatial and temporal variation in water quality when assessing the integrity of surface water bodies in agricultural watersheds.
Acknowledgements
We acknowledge Louisiana Department of Environmental Quality (LDEQ) and US Environmental Protection Agency Region 6 for funding this project. We are grateful to University of Louisiana at Lafayette staff and students who assisted in this project. We also thank USDA Natural Resources Conservation Service, local soil and water conservation districts, resource conservation and development councils, Louisiana Department of Agriculture and Forestry, and LDEQ staff and other stakeholders for helping us in this project. We would also like to thank the anonymous reviewers for their very helpful comments on the manuscript.
- © 2020 by the Soil and Water Conservation Society