1. Introduction
Eutrophication of freshwater ecosystems is a global problem, resulting in both ecological and societal challenges, including potentially harmful algal blooms, hypoxia, loss of biodiversity, and economic impacts to home value and ecotourism [
1,
2,
3,
4]. These challenges are particularly problematic in heavily agricultural systems, where runoff can carry significant loads of sediment and nutrients, resulting in legacy issues in terrestrial soils and lake sediments [
5].
Another challenge is identifying which conservation practices (or best management practices: BMPs) are most appropriate and where are the most desirable locations to place them, in order to use limited resources most effectively to achieve restoration goals [
6]. In other words, as provocatively asked by Sowa et al. [
7], how much conservation is enough?
One of the most significantly impaired lakes in Michigan is Lake Macatawa, which is the terminus of a highly degraded watershed and has exhibited the symptoms of a hypereutrophic lake for more than 40 years [
8,
9]. Extremely high nutrient and chlorophyll concentrations, excessive turbidity, low dissolved oxygen, and a high rate of sediment deposition make it one of the most hypereutrophic lakes in Michigan [
8,
9]. Nonpoint source pollution from the watershed, particularly agricultural areas, is recognized as the primary source of the excess nutrients and sediment that fuel hypereutrophic conditions in Lake Macatawa [
8].
Because of this nutrient enrichment, Lake Macatawa and all of its tributaries are included on Michigan’s list of impaired water bodies, prompting the issuance of a phosphorus Total Maximum Daily Load (TMDL) for Lake Macatawa in 2000. The TMDL set an interim target total phosphorus (TP) concentration of 50 μg/L in Lake Macatawa [
10]. In recent years, monthly average TP concentrations were greater than 125 μg/L, and at times exceeded 200 μg/L [
8]. Thus, meeting the TMDL target represents a major challenge in the Macatawa watershed. The TMDL estimated that a 72% reduction in phosphorus loads from the watershed would be required to meet the TP concentration target [
10].
A large-scale, multidisciplinary, collaborative watershed remediation project aimed at improving water quality in Lake Macatawa was initiated in 2013 and called “Project Clarity”. This public–private partnership adopted a holistic approach that includes wetland restoration, in-stream remediation, BMPs, and community education, with the ultimate goal of improved water quality in Lake Macatawa.
In this paper, we analyze pre- and postrestoration water quality data of two wetland restoration projects, as well as lake status, to evaluate the performance of these BMPs, as well as overall watershed condition. We also speculate on why water quality has not significantly improved in the two years following construction of the restored wetlands.
2. Methods
2.1. Overall Site Description
The Macatawa watershed (464 km
2) is located in Ottawa and Allegan Counties (MI) and includes Lake Macatawa, the Macatawa River, and many tributaries. It is dominated by agricultural (46%) and urban (33%) land uses, which have contributed to the loss of 86% of the watershed’s natural wetlands [
8]. The watershed includes the cities of Holland and Zeeland and parts of 13 townships [
8]. Lake Macatawa is a 7.2 km
2 drowned river mouth lake. It is relatively shallow, with an average depth of 3.6 m and a maximum depth of 12 m in the western basin. The Macatawa River, the main tributary to the lake, flows into the lake’s shallow eastern basin. A navigation channel at the western end of the lake connects Lake Macatawa with Lake Michigan.
A watershed-wide sediment sampling study, conducted in 2011 and 2012 (Hope College, Holland, MI, USA, unpubl. data) identified the sub-basins within the Macatawa watershed with the highest amounts of sediment. Based on these results, the Peters Creek, Upper Macatawa River, and North Branch sub-basins (
Figure 1) were given the greatest priority for restoration action. As a consequence, two wetland restoration efforts were targeted for these regions: Middle Macatawa and Haworth projects. The specific sites were selected due to land owners’ willingness to donate or sell the parcels. Our monitoring initiative is focused on these two key wetland restoration areas in the Macatawa watershed and Lake Macatawa (
Figure 1). Details on these two efforts are provided below.
2.2. Wetland Restoration: Middle Macatawa & Haworth Properties
The Middle Macatawa and Haworth properties were purchased as part of Project Clarity and designated for wetland restoration. Restoration goals included slowing the flow of water in the Macatawa River and its tributaries, particularly during high flow events, thus trapping and retaining suspended sediments and nutrients. The Middle Macatawa project involved reconnecting ~16 ha of former pastureland to the adjacent river by placing a pipe from the river to an excavated area in the floodplain. The spoils from the excavated area create a berm that surrounds the excavated area and detains floodwater, allowing sediment and nutrients to settle out. The Haworth project restored a 17-ha wetland, composed of four basins that collect and store water from adjacent properties and the river during periods of high flow. Restoration construction at Middle Macatawa and Haworth was completed in late September and early October 2015, respectively.
Monitoring sites were established upstream and downstream of each restoration area. The Middle Macatawa study area (
Figure 1B) has two upstream sites (Macatawa River (Macatawa Up) and Peters Creek), which both flow into the Macatawa River and one downstream site (Macatawa River at the USGS gauging station (Macatawa Down)). The Haworth study area (
Figure 1C) consists of monitoring locations upstream and downstream of the restoration area on the North Branch of the Macatawa River. Our focus in this paper is on water quality and hydrologic monitoring from December 2016 through November 2017, although we also include prior monitoring data. Sampling occurred monthly during baseflow conditions and during three storm events in 2017 (
Table 1); storm sampling was triggered when the local USGS gage reached 300 ft
3/s (~8.5 m
3/s). During each monitoring event, general water quality parameters (dissolved oxygen (DO), temperature, pH, specific conductivity (SpCond), total dissolved solids (TDS), redox potential (ORP)) were measured using a YSI 6600 sonde. Grab samples were collected for analysis of phosphorus (soluble reactive phosphorus (SRP), TP) and nitrogen (ammonia (NH
3), nitrate (
), and total Kjeldahl nitrogen (TKN)) species. Water quality measurements and sample collection took place in the thalweg of the channel at permanently-established transects. Duplicate water quality samples and sonde measurements were taken every other month during baseflow conditions and all storm events. All samples were placed in a cooler on ice until returned to the lab, usually within four hours, where they were stored and processed appropriately (see below).
Water for SRP and
analyses was syringe-filtered through acid-washed 0.45-μm membrane filters into scintillation vials; SRP was refrigerated and
frozen until analysis. NH
3 and TKN were acidified with sulfuric acid and kept at 20 °C until analysis. SRP, TP, NH
3,
, and TKN were analyzed on a SEAL AQ2 discrete automated analyzer (SEAL Analytical, Inc.: Mequon, WI, USA) [
11]. Any values below detection were calculated as ½ the detection limit.
Stream hydrographs were generated at each monitoring location using water level loggers and staff gages that were installed at permanently established transects at four of the monitoring locations (the Macatawa Down site did not require one because we used the USGS gage).
Turbidity sensors were deployed at the upstream and downstream locations on the main branch of the Macatawa River before snowmelt in March 2017. The in situ turbidity sensors are YSI 600OMS V2 (Xylem Inc.: Rye Brook, NY, USA), which uses the same technology as the YSI 6600 through the same YSI 6136 optical turbidity sensor. The 6136 turbidity sensor has a range of 1 to 1000 NTU with a resolution 0.1 NTU, and accuracy of ± 2% of the reading or 0.3 NTU, whichever is greater. It features a mechanical self-wiping capability to help prevent biofilm covering the optics. We placed the turbidity sensors in a protective housing made from schedule 80 PVC with vents at the bottom to allow water to pass through. The measurements are recorded every 30 min, and stored internally until downloaded once a month. The YSI sensors are calibrated on a monthly basis. The turbidity sensors were removed in December 2017 to avoid possible ice damage and were returned to their former locations before the final snowmelt in spring of 2018.
2.3. Data Analysis
Our analysis focused on identifying (1) upstream vs. downstream differences during baseflow and storm flow conditions, and (2) pre- vs. postrestoration differences in nutrients and turbidity. Upstream–downstream differences between site pairs (e.g., North Up vs. North Down) within the 2017 sampling year at baseflow and at storm flow were statistically tested using either a two-tailed paired t-test (normally-distributed data) or Wilcoxon signed rank test (non-normally distributed data). Baseflow and storm flow conditions were evaluated separately for each site pair. A one-way analysis of variance test (ANOVA; normally distributed data) or Kruskal–Wallis test (one-way ANOVA on ranks; non-normally distributed data) was used to compare data from the three Middle Macatawa sites simultaneously. ANOVAs that detected significant differences were followed by post-hoc Tukey pairwise comparison tests.
Pre- and postrestoration differences were statistically tested separately for each site using two-tailed paired t-tests at baseflow and either two-tailed unpaired t-tests (normally distributed data) or Mann–Whitney rank sum tests (non-normally distributed data) at storm flow. In order to remove seasonality as a potentially biasing factor in analyses and because not all samples were taken at the same time from all sites, paired t-tests for baseflow incorporated an equal number of samples (n = 16) from identical months in pre- and postrestoration periods (Apr., Jun., Jul., Sep., Oct., Nov., Dec., Jan., Feb., Mar., Apr., May, Jun., Jul., Aug., Sep.). Storm flow analyses incorporated all possible sampled storm events (prerestoration: n = 4 (North Up) or n = 5 (North Down, all Middle Macatawa sites); postrestoration: n = 6 (all Haworth and Middle Macatawa sites)).
Normality was tested using the Shapiro–Wilk test and equal variance was tested using the Brown–Forsythe test. Data not meeting test assumptions of normality and equal variance were transformed prior to analysis. Statistical significance was indicated by p-values < 0.05. Trends of marginal significance were indicated by p-values < 0.10. All statistical tests were performed using SigmaPlot 13.0 (Systat Software, Inc.: San Jose, CA, USA).
2.4. Lake Macatawa
Water quality monitoring in the lake was conducted at five sites during spring, summer, and fall 2017 (
Figure 1D). At each sampling location, general water quality measurements (DO, temperature, pH, specific conductivity, TDS, ORP, and turbidity) were taken using a YSI 6600 sonde at the surface, middle, and near bottom of the water column. Water transparency was measured as Secchi disk depth. Water samples were collected from the surface and near-bottom of the water column using a Van Dorn Bottle and analyzed for SRP, TP, and chlorophyll
a (Chl
a). Additional Lake Macatawa water samples for
, NH
3, and TKN were collected for the first time in 2017. Samples also were taken for phytoplankton community composition and archived for possible future analysis.
Water for SRP and
analysis was syringe-filtered through acid washed 0.45-μm membrane filters into scintillation vials and refrigerated until analysis. SRP was refrigerated and
frozen until analysis. NH
3 and TKN were acidified with sulfuric acid and kept at 20 °C until analysis. SRP, TP,
, NH
3, and TKN were analyzed as previously described. Chl
a samples were filtered through GFF filters and frozen until analysis on a Shimadzu UV-1601 spectrophotometer (Shimadzu Corp.: Kyoto, Japan) [
12].
4. Discussion
Nonpoint source pollution, especially from agricultural areas, is a major factor causing eutrophication and associated algal blooms in regions throughout the world [
13,
14,
15]. Despite a wide variety of mitigation options [
15] and planning tools [
16], impaired water quality continues to plague water bodies on a global basis. The proliferation of cyanobacterial blooms, and their associated toxins, is anticipated to worsen when combined with climate change [
17].
Many lakes and rivers in the upper Midwest, including portions of the Laurentian Great Lakes, have experienced significant eutrophication over the past few decades. Excessive nutrient inputs, exacerbated by extreme precipitation events that mobilize recently applied, and readily bioavailable, fertilizer [
4,
18,
19], have contributed to potentially harmful algal blooms, hypoxic and anoxic conditions, and associated ecological and socioeconomic impairments. These challenges persist despite the investment of hundreds of millions of dollars in BMPs, begging several questions: Are we using ineffective BMPs? Are we locating BMPs in the wrong areas? Should we be more patient for the BMPs to become more effective? Does the intensity of agricultural land use overwhelm the assimilative capacity of the BMPs? Or perhaps, as implied by Kleinman et al. [
16], is there sufficient satisfaction with implementation of the management practice (an output) instead of its effectiveness (an outcome) that we do not push harder for better outcomes?
Lake Macatawa and its watershed is one of the most impaired systems in Michigan. To address excess sediment and P loading, a public–private partnership was formed. The project partners have implemented a number of BMPs including wetland creation, bank stabilization, two-stage ditches, tile drain control structures, and no-till rotations. Monitoring has been incorporated to assess the performance of BMPs, with the explicit recognition that the problems facing the watershed accumulated over 150 years, and they would not be resolved overnight. It is clear from our data that in the few years since BMPs have been implemented, there has been little improvement in water quality either downstream of the created wetlands or in the downstream receiving water body (Lake Macatawa). This is not surprising given the short time period since implementation; continued monitoring will reveal if water quality improves over time as BMPs both increase in number across the watershed and mature, resulting in greater effectiveness.
The data from tributary monitoring to date indicate that in terms of DO, stream conditions are passable, as DO concentrations <5 mg/L, indicative of stress to many fishes, were not observed in our sampling. However, conductivity levels above 600 µS/cm were common, and they are generally indicative of human-induced stress in aquatic ecosystems [
20]. In addition, phosphorus concentrations were extremely high, especially during storm events, and nitrate concentrations in excess of 10 mg/L, which can be toxic to warm-blooded animals under certain conditions, were measured in Peters Creek. These high nitrate levels are likely associated with the application of fertilizer on fields [
21,
22].
We identify several caveats in our results, which help place our findings in a broader context. First, because BMP construction occurred less than two years prior to our monitoring, it is possible that the short-term impacts of soil movement during restoration activities resulted in increased nutrients. If so, these concentrations should decline once the restored wetlands become mature and are fully functional. This point applies to other BMPs, such as two-stage ditches [
23] and bank stabilization efforts.
The TMDL for Lake Macatawa was established in 1999, before the impacts of climate change were widely recognized. Today, the concerns over episodic storm events, and their ability to move sediment and nutrients, are much more acute [
18,
24,
25]. Hence, the 50 µg/L TP target may be even harder to obtain than previously realized, and effectiveness of BMPs becomes ever more important. The mean surface TP in Lake Macatawa ranged from 78 to 221 µg/L (
Table 11), far above the target. Concentrations do decline at sites farther westward in the lake, presumably due to the settling out of particles and dilution from high quality Lake Michigan water advecting into the western end of Lake Macatawa. Despite occasional spikes in bottom water SRP and TP, overall there was no evidence of systemic internal P loading in Lake Macatawa, which if present, would be indicated by very high concentrations of SRP and/or TP (>400 µg/L) in bottom waters, as we have measured in other west Michigan drowned-river mouth lakes [
26,
27]. The influence of groundwater on the lake’s ecology has not been addressed and was not considered a major factor as part of the lake’s TMDL [
10], presumably because the septic systems in direct lake discharge sub-basin (
Figure 1) have been replaced by sanitary sewers.
Both surface and bottom Chl
a values frequently exceeded the 22 μg/L hypereutrophic threshold commonly used by MDEQ in its assessments of Lake Macatawa [
9]. Mean Secchi disk depths indicated low transparency throughout the year, less than 1 m, suggesting eutrophic to hypereutrophic conditions [
28]; high sediment loads may help keep algal blooms from becoming even more problematic due to reduced light transmission [
29]. The lack of improvement in lake condition is not surprising as it often takes years, if not decades, for lake conditions to improve once the stressors are removed, and in many cases, the stressors remain in place but at reduced levels, exacerbating lake impairment [
5,
30].
The lack of improvement in the watershed and lake in 2017 can be attributed to at least four reasons: (1) Restoration is still very recent, and until the restored sites are fully functional, which should take a number of years, it is unreasonable to expect a demonstrable change; (2) the two created wetland restoration sites have relatively small footprints and volume holding capacity compared to the entire watershed; the two sites have a combined area of 0.45 km2 compared to the watershed area of 464 km2. Given the volume of water moving through the Macatawa River, especially during storm events, the ability to detect a signal from the noise may be very difficult at any one particular site; (3) the natural environment is variable, so it will take a number of years to detect a robust trend at any site, regardless of direction; and (4) 2017 was a dry year (43% lower than long-term average), thereby resulting in fewer opportunities for the wetlands to serve as filtering and retention basins to remove transport of pollutants.
Future needs in this watershed include continued monitoring of water quality conditions and implementation of targeted agricultural BMPs. For example, tile drain effluent appears to be an additional source of P (and maybe N) [
31], but has not received adequate attention in the past and is now being recognized as a factor contributing to toxic algal blooms in the western basin of Lake Erie [
32,
33]. We are in the process of developing a SWAT model for the watershed (Iavorivska, in preparation) which will allow managers to test scenarios to assess potential efficacy of different BMPs placed in different locations.