Comprehensive assessment of legacy organic contaminants and trends in lake trout from Cayuga Lake,New York: 2011–2017

A seven-year assessment of lake trout contaminants and statistical trend modeling was conducted on Cayuga Lake, New York as the reference monitoring site for the USEPA Great Lakes Fish Monitoring and Surveillance Program (GLFMSP). In all, over 200 individual compounds including polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), organochlorine pesticides (OCPs), polychlorinated napthalenes (PCNs), polychlorinated dibenzo-p-dioxins and polychlorinated dibenzofurans (PCDD/Fs) were monitored for the 2011–2017 study period. Legacy contaminants such as PCBs, OCPs, and PBDEs dominate the whole-fish lake trout contaminant profile in the 24–147 ng/g mean concentration range. Total dioxin-like contaminants monitored included PCDD/Fs, CP-PCBs, and PCNs with a mean concentration range of 2.0–516 pg/g. TEQ concentrations for the dioxin-like compounds ranged between 0.3 and 5.3 pg-TEQ/g, and total TEQ averaged 7.2 pg-TEQ/g for whole-fish lake trout. Contaminant trends were modeled using log-linear correlations for both whole-fish lake trout and lake trout eggs. All select compounds modeled for whole-fish lake trout showed significant decreases (average ?53.4%) over the study period. Tissue residue guidelines for the protection of wildlife were uniformly exceeded for both mammals and birds based on mean 2011–2017 total TEQ concentrations measured. All lake trout skin-on fillets from Cayuga Lake analyzed exceed the USEPA human health screening value of 0.15 pg-TEQ/g-ww. Whole-fish lake trout from Cayuga Lake have significantly lower mean total TEQs (6.5 pg-TEQ/g) as compared with the average total TEQ in the Great Lakes (range: 21.9–50.8 pg-TEQ/g) (U = 10.000, p < 0.001), although rates of contaminant decline, contaminant ratios, and half-life (t_1/2) are comparable.     

Relationship between lake trout spawning,embryonic survival,and currents: A case of bet hedging in the face of environmental stochasticity?

Lake trout, Salvelinus namaycush, spawning in the Great Lakes occurs primarily on cobble substrate at relatively shallow water depths that can experience strong water currents. Strong currents may limit embryonic survival by damaging or displacing eggs, but may also reduce the accumulation of fine material and limit foraging by potential egg predators. To better understand the importance of currents, we evaluated the role of currents in spawning habitat selection, egg density and survival, and egg predator density at a spawning reef in Lake Champlain (USA). Most spawning occurred one week after the largest storm event associated with the strongest currents and greatest upwelling. Highest spawning activity was associated with a relatively shallow part of the reef that had the highest current velocity and greatest potential for egg displacement. Within the interstices, the survival of naturally deposited eggs was unrelated to the concurrent loss of artificial eggs. We propose that the reproductive strategy of spawning on shallow areas of a reef that have the highest current velocity and high potential for egg loss represents a type of bet hedging to optimize survival of those embryos that remain within interstices. This strategy may have evolved in response to environmental stochasticity that resulted in higher egg survival.     

Results of the collaborative Lake Ontario bloater restoration stocking and assessment, 2012–2020

Bloater, Coregonus hoyi, are deepwater planktivores native to the Laurentian Great Lakes and Lake Nipigon. Interpretations of commercial fishery time series suggest they were common in Lake Ontario through the early 1900s but by the 1950s were no longer captured by commercial fishers. Annual bottom trawl surveys that began in 1978 and sampled extensively across putative bloater habitat only yielded one individual (1983), suggesting that the species had been locally extirpated. In 2012, a multiagency restoration program stocked bloater into Lake Ontario from gametes collected in Lake Michigan. From 2012 to 2020, 1,028,191 bloater were stocked into Lake Ontario. Bottom trawl surveys first detected stocked fish in 2015, and through 2020 ten bloater have been caught (total length mean = 129 mm, s.d. = 44 mm, range: 96–240 mm). Hatchery applied marks and genetic analyses confirmed the species identification and identified stocking location for some individuals. Trawl capture locations and acoustic telemetry suggested that stocked fish dispersed throughout the main lake within months or sooner, and the depth distribution of recaptured bloater was similar to historic distributions in Lake Ontario and other Great Lakes. Predicted bloater trawl catches, based on modeled population abundance and trawl survey efficiency, were similar to observed catches, suggesting that post-stocking survival is less than 20% and contemporary bottom trawl surveys can quantify bloater abundance at low densities and track restoration.     

Resuspension of Bottom Sediments in Lake Ontario During the Unstratified Period, 1992–1993

Time-series measurements of current velocity and water temperature made at a station in southwestern Lake Ontario between October 1992 and June 1993 are combined with analyses of material collected in a sequencing sediment trap moored at the same location. The results show that local resuspension of bottom material occurred several times during November and December. Analyses of total PCBs and Mirex concentrations were used to infer the origins of the material collected in the traps. Material collected during the fall has multiples sources; in addition to material already in suspension and bottom material resuspended from the immediate area, some material appears to be derived from the Niagara River delta. This is the first reported instance of bottom resuspension in the Great Lakes at depths below wave base during the unstratified period that has been confirmed by current velocity measurements.     

Spatial and seasonal variation in the phytoplankton community of Lake Victoria’s Mwanza Gulf,compared to northern parts of the lake

We investigated the phytoplankton species composition and abundance in two seasons in Mwanza Gulf, Lake Victoria (Tanzania). Phytoplankton was sampled and chlorophyll a content was measured in the dry and wet seasons of 2010–2011 at three stations, from the southern land-inward end of the Gulf towards the open lake. Cyanobacteria, mostly small colonial and filamentous species (e.g., Aphanocapsa spp., Planktolyngbya spp., Merismopedia spp.) dominated at each station (76–95 %), followed by Chlorophyta (5–21 %), whereas the contribution of Bacillariophyceae was small (0–6 %). Phytoplankton densities were generally higher in the rainy season and strongly increased going land-inward from the open lake. Low abundance of N-fixing phytoplankton species suggests that N-fixation was low. The chlorophyll a content in the mouth of the Gulf was low (mean values 4–6 µg/L) compared to values reported previously. Also, chlorophyll a values (means 11–14 µg/L) at land-inward stations of Mwanza Gulf were much lower than those in the northern gulfs (Napoleon Gulf, Murchison Bay and Nyanza Gulf). Between 2002 and 2009 the phytoplankton composition of Mwanza Gulf changed from a community mostly dominated by Bacillariophyceae into a community dominated by Cyanobacteria. In the open water of Lake Victoria, Bacillariophyceae and Cyanobacteria were both abundant. Cyanobacteria dominated both in the three northern gulfs and Mwanza Gulf, but all four showed substantial differences in species and genus compositions. Phytoplankton composition and abundance in Mwanza Gulf differs in many respects from the open water of Lake Victoria and its three northern gulfs.     

Reproductive Potential and Fecundity of Lake Trout Strains in Southern and Eastern Waters of Lake Ontario, 1977–1994

We assessed the reproductive potential of various genetic strains of hatchery lake trout (Salvelinus namaycush) in southern and eastern Lake Ontario from indices of fecundity and indices of male abundance. Indices were constructed from catches of mature lake trout in gill nets during September 1980 to 1994 after correcting for mortality from sea lampreys (Petromyzon marinus) which occurred between September sampling and late fall spawning. Strain and age were assigned to individual lake trout based on clipped fins and maxillary bones or coded wire tags. Fecundity-length relationships for fish of the same age, determined from mature females collected in 1977 to 1981 and 1994, were not different (P > 0.05) among genetic strains. For all strains combined, fecundity-length relationships in 1977 to 1981 were not different among fish of various ages but in 1994, age-5 and -6 fish had fewer eggs (P < 0.003) than age-7 fish, and age-7 fish had fewer eggs (P < 0.003) than fish of age 8, 9, or 10. Annual indices of fecundity varied 19 fold and indices of mature males varied 11 fold; both indices were low in the early 1980s, increased sharply in the mid 1980s, and peaked in 1993. The strain which dominated fecundity and mature male indices shifted during the study from Seneca Lake strain to Lake Superior strain and then back to Seneca Lake strain. However, changes in either reproductive potential or genotypes do not appear responsible for the abrupt appearance of naturally-produced yearling lake trout throughout southern and eastern Lake Ontario in 1994–1995, the first widespread occurrence of juveniles produced by hatchery lake trout in Lake Ontario.     

Assessment of lake–groundwater interactions and anthropogenic stresses,using numerical groundwater flow model,for a Rift lake catchment in central Ethiopia

A steady-state groundwater flow model (MODFLOW) was used to study lake and groundwater interactions in a complex rift volcanic catchment. It also was used to assess the effects of water pumping from wells, and of variable recharge rates associated with climate and lake level changes, on the dynamics of the volcanic aquifers surrounding Lake Awassa. The model simulations were made after first developing a reasonable conceptual model, on the basis of conventional hydrogeological mapping, pumping test and hydrometeorological data analyses, and from ancillary information obtained from hydrochemical and isotope techniques. The model results indicated that the lakes and Rift aquifers are fed by large groundwater inputs that originate in the highlands. The lakes and rivers have important roles in recharging the aquifers in some locations. Lake Awassa receives a major groundwater inflow from its southern and eastern shorelines, while substantial water leakage from the lake occurs along the northern shoreline. The annual groundwater outflow from the catchment is estimated to 52.5 × 10⁶ m³. Scenario analyses revealed that increasing the current pumping rate from wells by fourfold will substantially reduce the groundwater level substantially, although the regional flow pattern would remain the same. There appears to be no immediate danger to the Rift aquatic environment from the current water pumping rate. Drying the small Lake Shalo and associated swamps, however, will cause a large change in the water balance of the larger Lake Awassa. Slight changes in groundwater recharge can cause large differences in groundwater levels for most of the Rift caldera floor far from the lake shores. This study provides a reasonable foundation for developing detailed transient predictive models, which can then readily be used as a decision support tool for development and implementation of sustainable water resources practices.     

Spatial and temporal distributions of lake whitefish spawning stocks in Northern lakes Michigan and Huron, 2003–2008

Adult lake whitefish were tagged and released from the Big Bay de Noc (BBN) and Naubinway (NAB) stocks in northern Lake Michigan, and the Detour (DET) and Cheboygan (CHB) stocks in northern Lake Huron during 2003–2006 to describe their spatial and temporal distributions. The contemporary spatial distributions were compared with past distributions of the BBN and NAB stocks. Sixty-two percent of BBN tag recoveries occurred in Wisconsin waters during winter, spring and summer, but 83% of fall tag recoveries were made near the tagging site. Eighty-eight percent of the NAB tag recoveries were made in the management unit of tagging and 7% occurred into northern Lake Huron. Over 90% of the DET stock remained in the vicinity of the tagging sites regardless of the season, while 75% of the CHB tag recoveries were made in northwestern Lake Huron and 17% were made in Ontario. Based on regression tree analysis, there were strong stock, season, and year effects on movement distances, with weaker effects due to sex and length at tagging. Spatial distribution of the BBN stock changed from 1978–1982 to 2003–2008, but spatial distribution of the NAB stock did not. Substantial differences in movement and distribution existed among the four stocks, large seasonal differences in spatial distribution were found within some stocks, and lake whitefish exhibited strong spawning site fidelity. Present management unit boundaries are inappropriate for managing three of our four stocks, and agencies should consider developing single harvest limits for both northern Lake Huron and western Lake Michigan.     

Advances in eutrophication science within Lake of the Woods,a complex transboundary lake – Forward to the special section

Outside of the five Great Lakes, Lake of the Woods (LoW) is the largest transboundary lake shared by Canada and the United States. Over the last two decades, an international consortium of researchers has advanced our understanding of eutrophication science in this complex watershed. This introductory paper outlines eleven research articles included in this special section, within three main themes pertinent to LoW: external phosphorus sources to LoW and internal phosphorus loading, lake and watershed modeling, and new insights into cyanobacterial and harmful algal blooms (cHABs). We conclude with a roadmap to guide future transboundary water quality management in LoW, including remaining research gaps and future monitoring needs.     

Biomarkers of algal populations in phytoplankton,filamentous algae,and sediments from the eastern basin of Lake Erie 2003–2005

Algal carotenoids (n:16) and chlorophylls (n:7) were determined in phytoplankton (n:60), filamentous algae (n:14) and sediments (n:44), collected from two nearshore sites in eastern Lake Erie (Van Buren Point (VBP), 6.5 m; Point Gratiot (PG), 17 m), during summer and fall of 2003, 2004 and 2005. The most prominent biomarkers for diatoms (fucoxanthin), cryptophytes (alloxanthin), chlorophytes (chlorophyll-b) and cyanobacteria (zeaxanthin) revealed temporal variations in phytoplankton community composition, which were correlated to water temperatures: i) seasonal succession, from diatoms with some cryptophytes in June (cool), to an increased percentage of chlorophytes and cyanobacteria in August and September (warm); ii) differences between 2004 (cool; cryptophytes more abundant) and 2005 (warm; chlorophytes more abundant). Filamentous algae (chlorophytes, epiphytic diatoms; some cyanobacteria) varied in condition, according to levels of chl-a, pheopigments, and class biomarkers: high, at VBP (growing; decaying; bont/E suspect); low, at PG (mostly dead). Relative to phytoplankton, sediments were depleted in several biomarkers (chl-a; diadino-, neo- and violaxanthin) but enriched with others, particularly at PG (pheopigments; diatoxanthin; canthaxanthin, echinenone; alloxanthin). Sediment composition was characterized by strong differences between sites (chl-a, chl-b and fucoxanthin nearly 10-fold greater at VBP than PG) and increasing accumulation of biomarkers from year-to-year. A linear log–log function, relating total carotenoids to total chlorophylls (intercept −0.516, slope 1.054, r² 0.96), implied increased biodegradation among specimen types: levels of pigments (pmol/g ww) decreased three orders of magnitude, from phytoplankton, through filamentous algae, to sediments. Deviation from a 1:1 relationship indicated 1.7-fold depletion of carotenoids relative to chlorophylls at mid-range.     

Status of the amphipod Diporeia spp. in Lake Superior, 2006–2016

The amphipod Diporeia spp. has historically been an important component of the benthic food web of the Laurentian Great Lakes. The Great Lakes Water Quality Agreement included its population density as an indicator of ecological condition for Lake Superior, with target values of 220–320 m^?2 in nearshore areas (≤100 m depth) and 30–160 m^?2 in offshore areas (>100 m). To assess the status of Diporeia in Lake Superior, we used a probability-based lake-wide survey design to obtain estimates of Diporeia density and biomass in 2006, 2011 and 2016. A PONAR grab sampler was used to collect Diporeia at 50–53 sites each year, with approximately half in the nearshore (<100 m depth) region of the lake and half in the offshore. The mean area-weighted lake-wide density was 395 ± 56 (SE) m^?2 in 2006, 756 ± 129 m^?2 in 2011, and 502 ± 60 m^?2 in 2016. For all years, both density and biomass were greater in the nearshore than in the offshore stratum. The densities for 2006–2016 were 3–5 times higher than those reported from a lake-wide survey conducted in 1973 by the Canada Centre for Inland Waters. The severe declines in Diporeia populations observed in the other Great Lakes during recent decades have apparently not occurred in Lake Superior. Further research is needed to understand spatial and temporal variability of Diporeia populations in Lake Superior to enhance the utility of Diporeia density as an indicator of benthic condition.     

Rethinking the Concepts of Virtual Water and Water Footprint in Relation to the Production–Consumption Binomial and the Water–Energy Nexus

In the field of Ecological Economics, the need of using physical indicators to analyse economic processes, at the same time they serve as tools in decision making, has been lately highlighted. Virtual Water (VW) and Water Footprint (WF) are two useful indicators in achieving this objective, the first one from the perspective of production, the second one from that of consumption. This difference between them is interesting inasmuch as it allows to identify the subjects who are responsible for water consumption, whether producers or consumers, and proves both indicators’ potential when designing water management policies. In this work, we consider a hypothesis according to which there is a clear difference between the two concepts—Virtual Water and Water Footprint—and this difference, although evident in their respective conceptualizations, is not reflected in their estimations and applications. This is true to the point that the two concepts are often used as synonyms, thus wasting the enormous potential associated to their difference. Starting from this hypothesis, our objective is, first of all, to highlight this evident but ignored difference between VW and WF through a deep and thorough literature review of the conceptual definitions and contributions, the methodologies developed and the applications made regarding the two concepts. Second, we intend to make a conceptual and methodological proposition aimed at underlining the differences already mentioned and to identify responsibilities in water consumption. We do it by broadening the context of analysis and by integrating the production–consumption binomial and water–energy nexus.     

Concentrations and loads of nutrients and major ions in the Niagara River, 1975–2018

Environment and Climate Change Canada has monitored Niagara River water quality in support of the Great Lakes Water Quality Agreement since establishing a fixed site at Niagara-on-the-Lake in 1975. Using over 40 years of data from this site along with the Fort Erie location added in 1983, we examine the status and trends of concentrations and loadings of nutrients and major ions and assess evidence of sources between the two stations. Trends were observed for the majority of measured parameters and there is strong agreement between trends in concentrations and loadings which are generally higher at the downstream site; however, upstream/downstream differences indicate relatively little loading occurs along the length of the river itself. For total phosphorus (TP), inputs from Lake Erie via the Niagara River account for the majority of loading to Lake Ontario and, in some years, exceed the 7000 MTA Lake Ontario target. Between 2014 and 2018, we calculate the mean Niagara River TP loading to be 5275 MTA. We highlight the major changes in water quality constituents over time, including TP, and reveal increased seasonal consumption of TP and SiO₂, reflecting potential increases in the biological productivity in Lake Erie. The long and rich Niagara River dataset, which comprises year round sampling (including rare winter data), provides detailed tracking of changing Great Lakes water quality and could be further utilized to assess the impacts of climate change, improve understanding of diatom and harmful algal bloom dynamics, and enhance knowledge of in-lake major ion and nutrient dynamics.     

Reporting on the status,trends, and drivers of algal blooms on Lake of the Woods using satellite-derived bloom indices (2002–2021)

Despite significant declines in external phosphorus loads, Lake of the Woods continues to experience severe recurring cyanobacterial harmful algal blooms (cHABs) covering as much as 80% of the lake surface area. Satellite-derived bloom indices were used to assess the status, trends, and drivers of cHAB conditions for the period 2002 to 2021 in support of developing ecosystem objectives and response indicators for the lake. Areas of greatest potential concern, with the most prolonged bloom occurrences, were in the southeast of the lake. Significant decreases in bloom indices suggest the lake may now be responding to historical nutrient reductions. The greatest rates of decrease were within the main water flow paths, with little change in the more isolated embayments, suggesting flushing plays a key role in regulating regional bloom severity. Significant inter-annual variability in bloom phenology was observed, with blooms peaking later in recent years, which may be in response to climate-induced changes in the lake and watershed. The absence of a direct relationship between external phosphorus loads and annual bloom severity reflects the complexity of the lake’s response to eutrophication and the potential roles of other drivers including climate and a strong legacy effect of sedimentary nutrients. A case study of the 2017 bloom season captures the compounding interaction of meteorological variability and seasonal nutrient delivery in regulating the bloom response. Results highlight the need for greater understanding of seasonal and regional variability of bloom drivers to aid in forecasting the lake’s recovery under both nutrient management and climate change scenarios.     

Study of Drawdown–Drain Discharge Relationship and its Application in Design of Cost Effective Subsurface Drainage System in Mugogo Swamp, Busogo, Rwanda

Mostly the swamps in Rwanda are surrounded by volcanic hills with small streams flowing to discharge runoff and seepage water. Mugogo swamp is located in Busogo sector, Musanze district, North province. Total area of the swamp is approximately 50 ha. The swamp is surrounded by hills and elevated volcanic rocky terrains. Potato is the main crop cultivated in the swamp. The average production rate of potato is 7 MT/ha which is very low compared to 12 MT/ha in well drained areas. During rainy season seepage water and runoff water from the surrounding hills cause the waterlogged condition of the swamp and affecting the potato cultivation and land productivity. The remedial measure for this swamp is to divert separately the runoff and seepage water from surrounding catchment area and then remove the recharge water by pumping through a system of subsurface drains. Hydraulic head–drain discharge relationship can be fitted with quadratic equation. Equivalent drainable porosity and equivalent hydraulic conductivity are determined as 0.105 m/day and 0.34% respectively for drain depth of 40 cm from soil surface. Effective hydraulic conductivity in the soil profile shows that its average value in the top 15 cm of soil layer is 0.17 m/day and that in the remaining depth up to impermeable layer is 0.015 m/day. Third degree polynomial expressions are made for Head–hydraulic conductivity and head–drainable porosity relationships. The nonlinear relation of hydraulic conductivity and drainable porosity with drawdown shows that the proximity of Kinoni stream does not affect drainage parameters of the area because of less seepage from the stream. The study also reveals that adoption of 7 m drain spacing is very less if crop parameter is not considered and will result higher drain cost. Drainage coefficient of 5 mm/day is arrived considering the rainfall distribution, infiltration rate of soil, allowable water logging tolerance of potato crop. Required drain spacings are calculated for different drainage coefficients of 1, 2, 3, 4 and 5 mm/day under different drawdown conditions to plot subsurface drainage characteristic curves of the swamp. These curves are useful to directly read the drain spacing and drain depth for the required drainage coefficient without going for tedious calculations. Cost analysis shows that the ratio of drain spacing to drain depth can be a decisive factor to select best combination of drain depth and drain spacing. For drainage coefficient of 5 mm/day, optimum drain spacing-depth ratio is found as 7.2 with a cost of 0.689 million Frw/ha. For different drainage coefficients in the swamp, the drain depth of 1.5 m is crucial and optimum cost occurs at this depth. It is also found that any increase in drawdown beyond the drawdown at critical drain depth will not reduce the cost significantly.     

Biotic and Abiotic Factors Related to Rainbow Smelt Recruitment in the Wisconsin Waters of Lake Superior, 1978–1997

Lake Superior rainbow smelt (Osmerus mordax) recruitment to 12–13 months of age in the Wisconsin waters of Lake Superior varied by a factor of 9.3 during 1978–1997. Management agencies have sought models that accurately predict recruitment, but no satisfactory models had previously been developed. In this study, modeling was conducted to determine which factors best explained recruitment variability. The Ricker stock-recruitment model derived from only the paired stock and recruit data accounted for 63% of the variability in recruitment data. The functional relationship that accounted for the greatest amount of recruitment variation (81%) included rainbow smelt stock size, May rainfall, and bloater (Coregonus hoyi) biomass. Model results were interpreted to mean that recruitment was affected negatively by increased river flows from increased rainfall, and affected positively by the biomass of bloater, and those results were interpreted to mean that bloater mediated the effects of lake trout predation on rainbow smelt recruits. Model results were also interpreted to mean that stock size caused compensatory, density-dependent mortality on rainbow smelt recruits. Correlations observed here may be of value to managers seeking approaches to either enhance or control populations of this species, which is not indigenous to the Great Lakes.     

Is reduced benthic flux related to the Diporeia decline? Analysis of spring blooms and whiting events in Lake Ontario

Benthic monitoring by USGS off the southern shore of Lake Ontario from October 1993 to October 1995 provides a detailed view of the early stages of the decline of the native amphipod Diporeia. A loss of the 1994 and 1995 year classes of Diporeia preceded the disappearance of the native amphipod at sites near Oswego and Rochester at depths from 55 to 130 m. In succeeding years, Diporeia populations continued to decline in Lake Ontario and were nearly extirpated by 2008. Explanations for Diporeia's decline in the Great Lakes include several hypotheses often linked to the introduction and expansion of exotic zebra and quagga mussels (Dreissena sp.). We compare the timeline of the Diporeia decline in Lake Ontario with trends in two sources of organic matter to the sediments — spring diatom blooms and late summer whiting events. The 1994–95 decline of Diporeia coincided with localized dreissenid effects on phytoplankton in the nearshore and a year (April 1994 to May 1995) of decreased flux of organic carbon recorded by sediment traps moored offshore of Oswego. Later declines of profundal (> 90 m) Diporeia populations in 2003 were poorly associated with trends in spring algal blooms and late summer whiting events.     

Nexus meets crisis: a review of conflict,natural resources and the humanitarian response in Darfur with reference to the water–energy–food nexus

Darfur has been widely used as a case study by both those arguing for causality between environmental scarcity and war and those disputing it. This article challenges that approach by drawing on debates taking place within Darfur, reflecting on both the conflict and the humanitarian response. It argues that reviewing Darfur on its own terms makes a stronger basis to identify transferable lessons for interventions elsewhere. It considers water, food and energy, and finds that supporting governance is an essential theme for promoting economic recovery and laying a foundation for a well-managed water–energy–food nexus.     

Long-term changes in hypolimnetic dissolved oxygen in a large lake: Effects of invasive mussels,eutrophication and climate change on Lake Simcoe, 1980–2012

Lake Simcoe has been influenced by multiple environmental drivers over the past decades, especially by reductions in phosphorus (P) loading, climate change, and invasive species such as dreissenid mussels (DM) which became firmly established in 1996. We examined the cumulative impact of these drivers on the volume-weighted hypolimnetic dissolved oxygen concentration (VWHDO) below 18?m at station K42 in Kempenfelt Bay during ice-free seasons from 1980 to 2012. Hypolimnetic DO depletion began in early spring when thermal stratification was observable but weak and continued throughout the ice-free season until cooling sufficiently lowered water column stability. In comparison to the pre-DM invasion period (1980–1995), mean annual VWHDO_init was 2.4?mg?L^?1 higher in the post-DM period (1996–2012), VWHDO_min was 1.54?mg?L^?1 higher and the mean duration of the depletion period (L) was 16?days longer. Mean DO depletion rate (DR) and temperature adjusted DO depletion rate (DR_adj) were slightly lower (7% and 5%, respectively) after 1996. P controls and DM had a positive effect on VWHDO, presumably by lowering productivity and diverting organic matter away from the hypolimnon. However, longer L apparently offset improvements in VWHDO_min. If lengthening of L associated with regional warming continues, then additional efforts to reduce P loads will be necessary to achieve the goal of maintaining VWHDO_min above the target of 7?mg?O₂?L^?1 throughout the summer and fall.