The research, led by Angelique D. Dahlberg of the Minnesota Aquatic Invasive Species Research Center at the University of Minnesota together with colleagues from the U.S. Geological Survey’s Upper Midwest Environmental Sciences Center and Earth Science Laboratories, focused on St. Albans Bay, a 66.3-hectare embayment of Lake Minnetonka, a heavily used recreational lake near Minneapolis-St. Paul. Zebra mussels were first confirmed in the lake in 2010 and have since become well established. In July 2019, the team applied EarthTec QZ, an acid-stabilized ionic copper product registered by the U.S. Environmental Protection Agency specifically for dreissenid mussel control, at a target concentration of 60 micrograms of copper per liter maintained over ten days, with an initial pulse of 100 micrograms per liter. Robinson Bay, an untreated embayment of similar size on the other side of the lake, served as the reference site for comparison.

The researchers then monitored the ecosystem for four consecutive field seasons, from 2019 through 2022, measuring chlorophyll-a as a proxy for phytoplankton abundance, sampling zooplankton and benthic invertebrate communities, and deploying cages of test organisms at five sites in each bay. The caged animals included young-of-year yellow perch, fathead minnows, bluegill, and largemouth bass, along with juvenile fatmucket mussels, a native unionid species, and adult zebra mussels collected from the lake. After the treatment period, tissues from these animals were analyzed for copper concentrations using inductively coupled plasma methods at Iowa State University’s Veterinary Diagnostics Laboratory, allowing the team to compare bioaccumulation across species and between treated and reference bays.

One of the study’s most counterintuitive findings involved phytoplankton. Chlorophyll-a concentrations rose in the treated bay one day after treatment began, accompanied by a drop in water clarity measured by Secchi disk transparency, a pattern indicating an algal bloom rather than an algal kill. The authors suggest this may have occurred because the treatment coincided with a bloom already underway, or because the copper dose was simply too low to broadly suppress algae, which previous research indicates typically require concentrations above 0.318 milligrams per liter to show measurable effects. Chlorophyll-a and clarity returned to near pretreatment levels within two weeks, and elevated chlorophyll in the treated bay in 2020 may have reflected the reduced filtering capacity of a diminished zebra mussel population.

Zooplankton, long known to be among the most copper-sensitive freshwater organisms, showed the clearest nontarget response. Abundance in the treated bay dropped significantly within one day of treatment, and diversity, measured with the Shannon index, fell sharply and remained low through two weeks posttreatment. The declines were not uniform across taxa: daphniids and diaptomids, groups with documented copper sensitivity, decreased most, while cyclopoid copepods tracked the reference bay almost exactly, with mean densities deviating by less than two individuals per liter. By August 2020, one year after treatment, zooplankton abundance in the treated bay had rebounded to levels exceeding the reference bay, and diversity had returned to reference-bay levels, indicating a substantial community recovery within a single year.

The statistical analysis also revealed how noisy lake ecosystems can be. Permutational analysis of variance showed that zooplankton community composition differed at nearly every sampling date in both bays, treated and untreated alike, underscoring the naturally chaotic dynamics of plankton communities documented in long-term experiments. This natural variability complicates attribution: some of the posttreatment shifts likely stemmed from the copper, but others reflected seasonal succession, nutrient dynamics, and predator-prey interactions operating independently of the treatment. The authors note that abundance may be a more responsive and practical metric than community structure for managers evaluating short-term treatment effects over a few years.

Benthic invertebrates, by contrast, showed no clear treatment effect. Abundance and diversity in the treated bay remained statistically similar across nearly all sampling occasions, and community structure was dominated by the same 13 families, including chironomids, hyalellid amphipods, and asellid isopods, throughout the study. The one notable signal was a decline in Hyalellidae immediately after treatment, which the authors attribute to copper exposure, given that the amphipod Hyalella azteca has a reported 48-hour lethal concentration 50 of 87 micrograms per liter in similar water chemistry. Because hyalellids are an important fish food source, such declines could ripple through the food web, although the family appeared to recover in subsequent years. Differences in substrate between the two bays also limited direct cross-bay comparisons of benthic communities.

Among the caged animals, the picture was mixed and, in some respects, sobering. Adult zebra mussels in the treated bay suffered 32 percent mortality and accumulated the highest tissue copper concentrations of any species tested, averaging 40.85 micrograms per gram, consistent with laboratory toxicokinetic studies. Native fatmuckets accumulated substantial copper as well, averaging 26.41 micrograms per gram, yet remarkably, no fatmucket mortality occurred in the treated bay; the only deaths were four individuals in the untreated reference bay. The authors caution, however, that laboratory studies report juvenile fatmucket toxicity thresholds well below the treatment concentration, so delayed or sublethal effects on native mussels remain a concern worth further investigation. Among fish, fathead minnows were the only species showing treatment-related mortality, with survival falling from 84 percent in the reference bay to 38 percent in the treated bay, and their tissue copper concentrations averaged an order of magnitude above reference levels. The authors note that caged fish cannot escape the plume, whereas free-ranging fish typically avoid contaminated water, so the caged results likely represent a worst-case exposure scenario.

The findings carry practical weight for the growing number of managers contemplating open-water copper treatments. The study suggests that a low-dose copper approach can suppress zebra mussels while causing limited and largely reversible collateral damage: zooplankton recovered within a year, benthic communities within two, and phytoplankton effects were transient. But the work also highlights trade-offs that managers must weigh deliberately. Timing treatments to avoid periods when larval fish depend heavily on zooplankton, allowing recovery intervals of one to three years between applications, and clarifying whether the management goal is eradication or suppression below a threshold of ecological impact are all decisions the data can now inform. Complete eradication of an established population may not be achievable at low doses, but reducing mussels below damaging levels while preserving native community function appears within reach. As zebra and quagga mussels continue their spread across North America, this four-year experiment provides a rare, ecosystem-scale evidence base for choosing when, where, and how to deploy copper against one of the continent’s most consequential aquatic invaders.

Subject of Research: Field evaluation of low-dose copper molluscicide effects on nontarget freshwater organisms during zebra mussel suppression in Lake Minnetonka.

Article Title: Assessing Low-Dose Copper Treatment for Dreissenid Mussels: Effects on Nontarget Organisms

Article References: Assessing Low-Dose Copper Treatment for Dreissenid Mussels: Effects on Nontarget Organisms. (n.d.). https://doi.org/10.1007/s00267-025-02354-4

Image Credits: AI Generated

DOI: 10.1007/s00267-025-02354-4

Keywords: zebra mussels, copper molluscicide, invasive species, Lake Minnetonka, nontarget effects, zooplankton, benthic invertebrates, bioaccumulation, ecotoxicology, aquatic invasive species, dreissenid control, lake management

The research, published in Environmental Earth Sciences, was led by Yong Li and Zhongyao Liang of Xiamen University together with colleagues at the Tianjin Municipal Engineering Design and Research Institute and the Nanjing Institute of Geography and Limnology. The team built a multi-scale interpretable machine-learning framework and applied it to nearly four years of daily water-quality monitoring at Shahu Lake, a shallow, 45-square-kilometre water body lying at 1,093 to 1,102 metres above sea level where the foothills of the Helan Mountains meet the alluvial plain of the Yellow River. Like many lakes in arid and semi-arid China, Shahu Lake experiences pronounced swings in its eutrophication status, driven by artificial water replenishment and a warming, variable climate.

At the heart of the study is a model built on the XGBoost algorithm, a gradient-boosted tree method prized for its efficiency with noisy, multivariate environmental data. The researchers fed the model daily measurements of water temperature, chemical oxygen demand measured by the manganese method, ammonia nitrogen, total phosphorus, total nitrogen, turbidity, and the nitrogen-to-phosphorus ratio, all drawn from the China National Environmental Monitoring Centre between November 2021 and August 2025. After screening out eighteen anomalous readings, just over one percent of the record, the team retained a complete set of 1,370 daily observations, splitting them eighty-twenty into training and independent testing sets. Hyperparameters were tuned with Bayesian optimization combined with five-fold cross-validation.

The resulting model performed remarkably well. It achieved a coefficient of determination of 0.889 on the training data and 0.810 on the unseen testing data, with root-mean-square errors of 3.78 and 5.13 micrograms per litre respectively and Kling-Gupta efficiencies above 0.81 in both cases. Those figures indicate strong predictive accuracy and, crucially, good generalisation, meaning the model had not simply memorised the training record but had learned relationships that hold for data it had never seen. Yet prediction alone was never the goal. The real innovation lies in how the researchers interrogated the model, using the TreeSHAP recursive algorithm, an exact method for attributing each prediction to individual input variables, to peel apart the drivers of chlorophyll-a, the standard proxy for phytoplankton biomass.

What makes the study distinctive is that the interpretation is deliberately multi-scale. Most previous applications of SHAP in aquatic science have averaged feature importance across an entire dataset, a global view that can mask drivers that matter only during particular seasons or events. Li and colleagues computed global importance across all 1,370 days, regional importance within the summer bloom window of July to September, and local importance for the ten most extreme chlorophyll-a events on record. At the regional scale they normalised SHAP values within each feature to remove scale differences, a methodological refinement designed to make season-specific comparisons meaningful.

The global picture largely confirms classical eutrophication theory. Water temperature emerged as the single most influential driver, with the model showing chlorophyll-a responses turning sharply positive as waters warm from about 15 to 24 degrees Celsius before plateauing at higher temperatures. Chemical oxygen demand, a proxy for organic matter and its mineralisation, ranked second, pointing to the role of internal nutrient recycling in shallow lakes, where decaying organic matter on the lakebed releases nitrogen and phosphorus back into the water column. Among nutrients, phosphorus consistently outweighed nitrogen in long-term importance, echoing the foundational whole-lake experiments of David Schindler and supporting the long-standing paradigm that phosphorus control is the backbone of lasting eutrophication management.

But the regional analysis revealed something the global view had concealed. During the July-to-September bloom period, the hierarchy of drivers shifted markedly. Nitrogen-related variables, particularly ammonia nitrogen and the nitrogen-to-phosphorus ratio, rose in importance, while the system effectively transitioned from chronic phosphorus limitation to transient nitrogen limitation under elevated thermal conditions. The response to water temperature peaked narrowly between 23 and 24 degrees Celsius, matching the optimal thermal niche for rapid cyanobacterial proliferation. The authors interpret this seasonal nitrogen stress not as an independent cause of blooms but as a consequence of them: as phosphorus fuels rapid growth, phytoplankton draw down available nitrogen until it becomes the immediate bottleneck, constraining bloom persistence and even influencing how much Microcystis biomass sinks out of the water column.

The local-scale analysis added a further layer of nuance. When the researchers examined the ten highest chlorophyll-a events individually, they found that each followed a distinct decision pathway through the model, even though the samples shared broadly similar environmental backgrounds. Water temperature contributed a strong, consistent 15.6 to 17.5 micrograms per litre in every case, functioning as the essential prerequisite for extreme blooms, with chemical oxygen demand adding a stable 3.65 to 5.08 micrograms per litre. Nutrient contributions, however, diverged sharply from event to event. Under persistently high phosphorus, the influence of total nitrogen versus the nitrogen-to-phosphorus ratio shifted depending on immediate nitrogen availability, and in a few instances nitrogen variables actually registered negative contributions, which the authors attribute to Liebig-style co-limitation by light, temperature, or phosphorus rather than any evidence that nitrogen depletion suppresses blooms under natural conditions.

These findings carry direct management implications. The framework argues for a dual-tiered strategy: long-term recovery still demands sustained reduction of both external and internal phosphorus loads, but during the summer bloom window managers may need adaptive controls on short-term nitrogen fluxes as well. Ignoring the transient nitrogen limitation, the authors warn, could foster nitrogen-depleted conditions that favour cyanobacterial dominance, a concern consistent with recent dual-nutrient reduction studies on Lake Erie and other systems. More broadly, the study positions interpretable machine learning as a state-dependent decision-support tool, capable of telling lake managers not just what matters on average but what matters now, in this season, under these conditions.

The authors are careful to note the limits of their approach. Because the analysis rests on observational data from a single monitoring station, the identified relationships are statistical associations rather than strict causal mechanisms, and the specific feature attributions may not extrapolate to lakes with different data distributions. Nonetheless, the framework’s central lesson is likely to resonate far beyond Shahu Lake: eutrophication is not governed by a single rule but by rules that change with time and state, and only analytical tools that respect those shifts, combining high-frequency monitoring, gradient-boosted models, and exact multi-scale attribution, can hope to keep pace with blooms that are becoming more frequent and more intense in a warming world.

Subject of Research: Multi-scale interpretable machine learning analysis of chlorophyll-a dynamics and eutrophication drivers in Shahu Lake

Article Title: Unraveling multi-scale drivers of Chl-a dynamics via interpretable machine learning: a case study of Shahu Lake

Article References: Li, Y., Huang, Y., Lin, L., Kong, X., & Liang, Z. (2026). Unraveling multi-scale drivers of Chl-a dynamics via interpretable machine learning: a case study of Shahu Lake. Environmental Earth Sciences, 85(15), Article 398. https://doi.org/10.1007/s12665-026-13133-7

Image Credits: AI Generated

DOI: 10.1007/s12665-026-13133-7

Keywords: chlorophyll-a, algal blooms, eutrophication, Shahu Lake, interpretable machine learning, XGBoost, SHAP, phosphorus limitation, nitrogen limitation, water quality monitoring, cyanobacteria, lake management