DOI: 10.5061/dryad.mpg4f4qzw
title: Time scales of ecosystem impacts and recovery under individual and serial invasions

description: Introductions of keystone or engineering species have complex and
 long-term impacts on multiple ecosystem features. Anticipating these
 consequences requires knowledge of the magnitude and pace of a species’
 impacts and whether they subside as resident species and communities
 adapt. Managing for the effects of invasive species is particularly
 complicated by a paucity of long-term ecosystem studies, the fact that
 invaders can represent novel ecological types, and as serial invasions by
 new species mediate the impacts of an original invader. We resolve the
 impacts of quagga (Dreissena rostriformis bugensis) and zebra mussels (D.
 polymorpha), two of the most widespread species that increasingly
 co-invade and re-engineer energy flows in freshwater ecosystems. Following
 many-decade time series and seven ecosystem features, we find remarkably
 similar ecosystem responses to mussel invasion across seven lakes in
 Europe and North America. Lakes invaded by zebra mussels only experienced
 the most severe changes within 5-10 years of invasion, followed by a
 partial recovery of several key ecosystem features. However, recovery
 disappeared, and initial impacts amplified in lakes where subsequent,
 serial invasions by quagga mussels competitively displaced zebra mussels.
 Our results suggest that the impact of quagga is stronger per unit biomass
 due to their greater effects on phytoplankton during spring and fall.
 Thus, we show that the effects of species introductions on many lake
 ecosystem features can manifest quickly and partially subside, except when
 multiple species invade a system.
description: <strong>Study area and history of
 invasion</strong> Studied lakes differed in the duration
 of <em>Dreissena</em> spp. invasion and the amount of
 information available. Lake Lukomskoe has the longest dataset of zebra
 mussel invasion (1972-2008) but has the least amount of available data,
 especially for the pre-invasion period. Due to the lack of primary data
 for Lake Lukomskoe we used published data, mostly summer or growing season
 averages (see below). Much more information (including primary data) is
 available for Narochanskie lakes and Oneida Lake for both the pre- and
 post-invasion periods. All information on lakes Eem and Veluwe was
 retrieved from Noordhuis et al. (2016). Lake
 Lukomskoe is the
 5<sup>th</sup> largest lake in Belarus with a wide littoral
 zone covered predominantly with sandy or silty sand sediments and
 occasionally with silt, and a profundal zone (&gt; 6 m) with silt
 (Karatayev 1983). In 1969, Lake Lukomskoe became a cooling reservoir of
 the Lukoml Thermal Power Station. Zebra mussels were first found in this
 lake in 1972; no <em>Dreissena</em> were found during benthic
 studies in 1968 and 1969 (Lyakhnovich et al. 1982). Pre-invasion data for
 Lake Lukomskoe are limited. Before zebra mussel invasion, the lake was
 eutrophic with summer Secchi depth of 1.8 m in 1955 and 2.0 m in 1965
 (reviewed in Lyakhnovich et al. 1988). Phytoplankton summer biomass for
 the pre-invasion period was estimated as 10.9 gm<sup>-3</sup>,
 macrophyte coverage as 5% of the lake surface area, zooplankton biomass at
 1.7 gm<sup>-3</sup> and zoobenthos biomass (excluding
 molluscs) as 0.6 gm<sup>-2</sup> (review in Lyakhnovich et al.
 1988). Following the invasion of <em>Dreissena</em> in late
 1970s - early 1980s, the lake became meso-oligotrophic (Karatayev 1983,
 Lyakhnovich et al. 1988). However, after the establishment of fish
 hatchery in cages targeted at 1000 tons of fish production per year, 3000
 tons of organic matter with large amount of phosphorous (about 1% of fish
 food) started entering the lake, causing water quality deterioration,
 eutrophication, and oxygen decline (Mitrakhovich et al.
 2008). To
 evaluate <em>Dreissena</em> population dynamics and their
 ecological impacts on Lake Lukomskoe, we analyzed data available for the
 following years: Secchi depth: 1972, 1973-1981, 1983-1985, 1987,
 1989-1992, 1994, 2003-2005, 2008, phytoplankton biomass: 1972, 1974-1979,
 1983, 1984, 1987, 1989-1992, 1994, 2004-2005, macrophyte coverage: 1973,
 1975, 1976, 1978, 1989, 1992, 2006, zooplankton biomass: 1972, 1974-1984,
 1989-1992, 1994, 2003-2005, biomass of benthic invertebrates excluding
 molluscs: 1972, 1974-1978, 1981-1983, 1985, 1987 1989-1992, 2003-2005,
 and <em>Dreissena</em> lake-wide biomass: 1972, 1974-1978,
 1989, and 2005 (Karatayev 1983; Karatayev et al. 1997; Lyakhnovich et al.
 1982, 1988, Mitrakhovich et al. 2008). Oligo-mesotrophic Lake
 Naroch and mesotrophic Lake Myastro were
 colonized by zebra mussels in mid-1980s (Burlakova et al. 2006). These
 lakes are among the best-studied lakes in the former Soviet Union. Data on
 water chemistry and major aquatic communities (phytoplankton, zooplankton,
 and benthos) are available since 1978 (Winberg 1985, Ostapenya et al.
 1993, 1994, 2012). Zebra mussels were found in 1984 in Lake Myastro
 (Burlakova 1998) and in 1989 in Lake Naroch (Ostapenya et al.
 1993). Oneida
 Lake at
 207 km<sup>2</sup> is the largest lake within the borders of
 New York State. It is a shallow, polymictic, meso to eutrophic lake with
 some bays and shoal areas covered with macrophytes down to 6 m, although
 coverage varies with fetch and substrate (Fitzgerald et al. 2016). The
 bottom substrate is variable and consists of rocky, sandy and silty areas,
 with most of the bottom below 9 m in silt. Few mussels were found in this
 deeper area until the expansion of quagga mussels in 2008 – 2009
 (Hetherington et al. 2019). The lake is an important fishing lake, with
 large populations and fisheries for walleye (<em>Sander
 vitreus</em>) and black bass (<em>Micropterus sp</em>.).
 The lake has been sampled for a variety of variables, for some as far back
 as 1956 (Rudstam et al. 2016). All data used in this paper are available
 online through the Knowledge Network for Biocomplexity (Rudstam 2020a, b,
 c, d, e) and these references also includes details of sampling and
 analysis methods. Analysis of the ecosystem effects of zebra mussels goes
 back to the early years of the invasion (Mellina et al. 1995, Horgan and
 Mills 1999, Idrisi et al. 2001, Mayer et al. 2002) and mussel effects are
 a central theme of several chapters of a 2016 book on Oneida Lake (Rudstam
 et al. 2016). Lake
 Balaton with its 596
 km<sup>2</sup> surface, and 3.25 m average water depth (based
 on bathymetric data set, Zlinszky et al. 2010), is the largest shallow
 lake in Central-Eastern Europe. Functionally it is a tourist lake, serving
 recreation for over millions of people in the summer. The lake elongates
 in the west-east direction with the maximum length of 78 km. The dominant
 sediment type is soft mud with a low (2-4%) organic content (Máté 1987)
 and a high content (50-60%) of Mg-bearing calcite (Müller and Wagner 1978,
 Pósfai et al. 2019). Hard substrates occur in small areas with the
 exception of the rocky shoreline protection structures occupying about
 half of the total shoreline lengths of 236 km. Stemming from the lake
 morphometry, a strong west-east trophic gradient can be measured in summer
 along the longitudinal axis resulting in eutrophic conditions
 (Chlorophyll-a &gt; 20 µg/L) in the western basin, whereas
 oligotrophic conditions (Chlorophyll-a &lt; 5 µg/L) prevail in the
 eastern basins (Sebestyén et al. 2017). Lake Balaton was undergoing a
 rapid eutrophication starting in the 1970s, due to the enlarged load of
 phosphorus from agricultural fields and communal sources of the tributary
 (Herodek 1986). Control of the nutrient input by building reservoirs and a
 drainage system along the lake residential areas have likely reversed the
 eutrophication after 1990s (Padisák and Reynolds 1998, Istvánovics et al.
 2007). Zebra mussels were introduced to Lake Balaton much earlier (around
 the early 1930s, Sebestyén, 1937) than quagga mussels (first observation
 was in 2008, Balogh and Purgel 2012,  Balogh et al. 2018), and by 2018,
 zebra mussels almost disappeared from the eastern basin, while a longer
 period of co-occurrence was observed in the western basin (Balogh et al.
 2019). While data are available for both the eastern and western basins of
 Lake Balaton, the western basin experiences unusually high levels of
 organic matter input. Therefore, we limit our analysis to the eastern
 basin data throughout. <strong><em>Dreissena</em></strong><strong> sampling protocol</strong> The most detailed <em>D. polymorpha</em> study in Lake Lukomskoe was conducted in the summer of 1978, when samples were collected from 14 transects regularly distributed throughout the lake running perpendicular to the shoreline (Karatayev 1983). At each transect samples were collected from 0.5, 1, 1.5, 2, 3, 4, 5, 6, and 8 m depth. From 0.5 to 6 m depth samples were collected by SCUBA divers using quadrat (0.25 m<sup>2</sup>), and an Ekman grab (0.025 m<sup>2</sup> sampling area) was used at 8 m depth. In 1989 and 2005, <em>Dreissena</em> samples were collected from 7 of the 14 transects sampled in 1978 using quadrats and bottom grabs (Peterson grab on hard substrates and Ekman grab on soft unconsolidated sediments, both 0.025 m<sup>2</sup>) (Karatayev 1983, Mitrakhovich et al. 2008). At each transect 8 – 10 stations were sampled at depths similar to those sampled in 1978. All pre-1978 <em>Dreissena</em> surveys were conducted using bottom grabs. We do not have the information on the amount of sampling stations for most of these years, however most likely they ranged between 20 and 60 stations per sampling event.  Detailed sampling protocol for <em>Dreissena</em> density and biomass assessment in the Narochanskie Lakes can be found in Burlakova et al. (2006). In brief, samples were collected in Lake Naroch from eight permanent transects in 1990, 1993, 1994, 1995, and 1997 (Fig. 1) and latter in 2005 (Mastitsky et al. 2006). In Lake Myastro samples were collected from 5 transects per lake in 1993 and 1995. For all years and all lakes, transects were initiated on the shore and ran perpendicular to the shore toward the center of the lake. For each transect, up to 10 replicate samples were collected at 0.5, 1, 1.5, 2, and 3 m depth and then at an interval of 1 or 2 m down to the maximum depth where zebra mussels were found. In 1990-1997 samples down to 2 m were collected using the quadrat method (0.25 m<sup>2</sup>). In Lake Naroch in 1997 quadrat samples were collected down to 7 m depths with the aid of a surface supplied air diving system (Pioneer 230 X, Brownie’s Third Lung). A similar sampling design was implemented in Narochanskie lakes in 2016-2018 when divers collected samples from 0.5 to 8 m depth using 0.25 m<sup>2</sup> quadrat. Deeper samples in all lakes were collected with an Ekman grab on soft sediments or a Peterson grab on hard sediments (both 0.025 m<sup>2</sup>). All quadrat and grab samples collected from all Belarusian lakes were washed through a 500 µm mesh, and within 48 h of sampling all zebra mussels larger than 1 mm maximum dimension were counted, opened with a scalpel to remove water from the mantle cavity, and the total sample was weighed to the nearest 0.01 g after being blotted dry on absorbent paper (wet weight, soft tissue plus shell). Oneida Lake dreissenid sampling developed over time (details in Hetherington et al. 2019). Mussels were collected at 8 - 10 sites (1992-2002) and at 15-26 sites (2003-2013) by SCUBA divers. Triplicates were taken at each site (1992-2002) and at most sites (2003-2013). Since 2013, Ekman grabs (1-3 replicates per site) were used on soft and sandy substrates and the number of sites was increased. SCUBA divers sampled rocky substrates as in past years. Mussels were returned to the lab the same day and frozen until processing time. Mussels were then thawed, identified to species, counted, and measured. Shell-on wet weights were calculated based on the mid-point of the length group using the regression shell-on wet weight = exp(2.97*ln(L, mm)-8.99) derived for zebra mussels in 1993 (Rudstam 2020a). Whole lake abundance estimates are based on the proportion of the lake bottom in different substrate and depth regions (shallower than 9m – rock, sand and soft substrate, deeper than 9 m). In Lake Balaton, a detailed dreissenid sampling procedure was described earlier (Balogh et al., 2008). Briefly and referring to the differences, macroinvertebrate samples were collected in May and July in Keszthely, Szigliget/Fonyód, Tihany and Balatonalmádi from 2003-2005 and 2009-2014 and in Kesztely and Tihany 2008 and 2015-2018. Three stones were collected randomly by hand from the upper (near the surface), about 10–20 cm below water level and lower (near the bottom of the lake) portions from the rip-rap. The deepest station is located at Tihany, where samples were collected by a scuba diver. Mussels were removed by knife, sieved (mesh size: 300 lm) and preserved in 70% ethanol. Dreissenids were identified by species and size. To calculate the surface area of sampled stones, the entire surface of the stones were traced onto wrapping paper, which were cut, weighed, and an algorithm for paperweight vs. surface area was used. The density of dreissenids was represented as ind/m<sup>2</sup> stone surface. For biomass calculation, we used the density and the body size-body mass relationship applying to the lengths of mussels (Muskó and Bakó 2005, Balogh et al. 2019, Báldi et al. 2019). Mussel biomass is provided as g of total wet weight of mussels with shells per m<sup>2</sup> of lake bottom surface. <strong><em>Environmental variables</em></strong> Environmental variables analyzed in this paper as response variables (potentially affected by <em>Dreissena</em>) to evaluate the impacts of <em>Dreissena</em> spp. invasion included transparency (as estimated by Secchi depth, m), total phosphorus concentration (mg/L, not studied in Lake Lukomskoe), macrophyte coverage (% of the lake bottom), chlorophyll <em>a</em> (g/m<sup>3</sup>, not studied in Lake Lukomskoe), and wet biomass of phytoplankton (g/m<sup>3</sup>), zooplankton (g/m<sup>3</sup>), and benthos excluding molluscs (g/m<sup>2</sup>). In Oneida Lake benthos biomass was not determined and for the model density was used instead. Dry mass of zooplankton calculated for Oneida Lake was multiplied by 5 to obtain wet biomass. Lake Lukomskoe No information was available on the amount of sampling stations and sampling events per year for Lake Lukomskoe, especially for pre-1977 period. Starting in 1977, phytoplankton and zooplankton in the lake was sampled at two pelagic stations at 8.0 m depth from surface, 2, 4, 6, and 7.5 m depths, usually 3 – 10 times per growing season. Phytoplankton was sampled using 1L Ruthner bathometer. Immediately after collection 0.5 L samples were preserved with acid Lugol’s solution and kept in amber glass bottles in dark for about 10 days. Then, the top layer of water was siphoned, and samples were concentrated to 100 mL. Zooplankton was collected with 10 L Vovk sampler and filtered through 64 mm mesh. Samples were concentrated to ca. 150 mL and fixed with formalin. Zooplankon groups (cladocerans, copepods, and rotifers) were counted under a dissection scope and identified to the lowest practical level using compound microscope. Biomass was estimated by multiplying the density of each species by their individual standard wet weight (same weight was used for the whole time period). Zoobenthos was sampled from 15 – 20 stations (occasionally up to 46 stations) 3 to 10 times per year and washed through a 500 mm mesh. All macroinvertebrates collected were fixed with 10% buffered formalin, identified to the lowest practical taxonomic level (usually species, genus or family), counted, blotted dry on absorbent paper and weighed to the nearest 0.1 mg (total wet mass). Narochanskie Lakes To study chlorophyll, phytoplankton, and zooplankton in Narochanskie lakes, samples were collected monthly since 1978 during the vegetation period (May–October) from standard monitoring stations: two pelagic stations in Lake Naroch, and one station in Lake Myastro. Samples were collected using two-liter Ruttner sampler from six depth layers (0.5, 3, 6, 8, 12, and 16 m) in Lake Naroch and from four depths (0.5, 4, 7 and 9 m) in Lake Myastro. Water samples from all the depths for each lake were mixed proportionally to the fraction of the depth layer in the total volume of the lake to obtain an integral sample representing the mean composition of lake water. Secchi disk transparency was measured with a 30 cm diameter white disk. To determine chlorophyll concentration (without correction for the presence of pheopigments) lake water was filtered using “Nucleopore” nuclear membrane with a pore diameter 1.5 μm. Chlorophyll analysis was performed spectrophotometrically after extracting pigments in 90% acetone (SCOR-UNESCO, 1996; Kovalevskaya et al. 2020). Phytoplankton samples (0.5 L) were fixed with Utermöhl’s solution (Mikheyeva 1989) and analyzed using Zeiss Axiolab microscope. A Fuchs-Rosenthal chamber with a volume of 3.2 mm<sup>3</sup> was used to count small phytoplankton; larger species were counted in a 1 mL chamber, while large colonial organisms (such as <em>Gloeotrichia</em> <em>echinulata</em>, <em>Volvox</em>) were counted using Bogorov chamber. The phytoplankton density was expressed in cell number (number of one-celled species, number of cells in filaments and colonies) per liter, and phytoplankton biomass was estimated using the biovolume of individual algal species assuming a specific density of 1 g/cm<sup>3</sup>. Biovolumes were obtained by measuring cell sizes and comparing cell shapes with geometric figures (Hillebrand et al. 1999, Mikheyeva 1999). For zooplankton samples, 10 L of water was filtered using 64 µm mesh size and fixed in 4% formalin. Individual zooplankton were identified to species or genera, counted and measured using Zeiss Axiolab and Zeiss Stemi 2000 microscopes. The body mass of Cladocera and Copepoda was determined using length/volume relationship assuming a specific density of 1 g/cm<sup>3</sup>. Rotifer species body volume was calculated using geometric figures most comparable to each species’ body shape using equations in Balushkina and Winberg (1979). The specific density of zooplankton was set to 1 g/cm<sup>3</sup>. Zoobenthos samples were taken usually three times per year along a nearshore to offshore transect. For Lake Naroch sample depths along the transect ranged from 1 to 16 m and for Lake Myastro from 1 to 10 m. In lakes Naroch and Myastro, samples were collected at every two meters (a total of 9 stations were sampled in Lake Naroch, 6 stations in Lake Myastro). Samples were collected with Petersen grab on hard and Ekman-Burge sampler with the Borutsky modification on soft substrates (sampling areas of both grabs 0.025 m<sup>2</sup>), washed through a 265 mm mesh, and fixed with 10% neutral buffered formalin. Macroinvertebrates were identified to the lowest possible taxonomic level (usually species, genus or family) and counted. The biomass of individuals or groups of zoobenthos was determined by weighing on a torsion balance after blotted dry on absorbent paper. Oneida Lake Oneida Lake limnology data has been collected since 1975 from four or five sites at weekly intervals from soon after ice-out through late fall. To make data comparable with other lakes (e.g., Naroch, Myastro), we used the average of weekly measures collected in May through October, while for comparison of zebra vs. quagga mussel impact in Oneida Lake for Secchi depth, chlorophyll and total phosphorous we used data collected in March through November. Integrated water samples were collected with a 1.9 cm inner diameter tygon tube (Nalgene) lowered from the surface to 0.5 m above the bottom, then closed at the top, retrieved, and emptied into 4L plastic bottles or into a plastic carbuoy. Within 2 to 4 h after collection, samples for chlorophyll <em>a</em> were filtered onto a Whatman 934-AH filter and frozen. Filters were extracted with 90% acetone and measured spectrophotometrically (Strickland and Parsons 1972). Secchi disk transparency was measured with a 20 cm diameter black and white disk. Subsamples of integrated lake water were preserved in Lugol's solution (1975 to 1995) or glutaraldhyde (1996 to present). From 1975 to 1995, the phytoplankton samples were settled in Utermöhl settling tubes at 10 ml aliquots for bloom phase samples and 25 ml aliquots for clear water phase samples. The samples were left to settle over a 24 hr period, after which the settling slide was placed in a Wild inverted microscope for identification and enumeration. Identification was made at two different magnifications (400x and 100x) to the species level when possible; otherwise, phytoplankton cells were identified to the genus. Phytoplankton were identified to species or higher taxonomic level and abundances were reported in #/mL. Biovolumes for this time period was estimated from standard species specific values. From 1996 onwards, phytoplankton samples from different stations were pooled prior to analyses and 8 to 14 weekly samples were counted by Phycotech Inc. No samples were processed from 1998 and 1999. Phycotech counted and measured 400 natural units including all algal cells that were alive at the time of sampling (cells with content) using magnifications of 250x, 500x and 1250x. The 500x magnification was the primary one used. The algae were mounted on slides and the amount of water mounted varied with the density of algal cells (0.1 to 100 mL/slide). Algal taxa were identified to species when possible and measured to provide sample specific biovolumes for the sample, by cell and by natural units (colonies, paired or other multiples of cells). Zooplankton were collected from the same sites using a 153-um mesh nylon net (0.5 m diameter) towed vertically from approximately 0.5 m off the sediment surface to the water surface. The efficiency of the net was measured with flow meters from 1999 to present. If flow meter readings indicate a malfunction or human error (efficiencies below 50% and above 125%) and when no flow meter was used, we assumed an efficiency of 87.4% (average of the 1999 to 2010 sampling period 87.4% SD 9.5%, N=1655). Flow meters were calibrated each year. Samples were preserved in 8 % buffered sugar-formalin solution (1975-1996) or 70 % ethyl alcohol (1997-present). Crustacean zooplankton were identified to species, counted and measured. A minimum of 100 animals were counted and measured from each sample. Biomass for individual species was calculated using length-weight regressions based mainly on Bottrell et al. (1976) and summarized in Watkins et al. (2011). Rotifers were not counted in Oneida Lake. Benthos samples were collected from three sites using a 152x152 mm square (0.023 m<sup>2</sup>) Ekman grab sampler. Samples were passed through a 253 µm mesh screen, preserved in ethanol, and stained with rose bengal. Organisms were removed from samples without a hand lens or microscope, and categorized to various taxonomic levels, mostly family or higher level. Lake Balaton Water quality parameters were recorded using a Horiba U-10 water multi-parameter measuring instrument and WTW Profiline Multi 3320 parameter device (Xylem Analytics, Germany). We used Secchi disc to evaluate water transparency. Depth, temperature, conductivity and pH were measured during the sampling. Chlorophyll-a concentration was determined spectrophotometrically after hot methanol extraction using the absorption coefficients determined by Wellburn (1994) following concentration of water samples on glass fiber filters (Macherey-Nagel; GF-5; nominal pore size is 0.4 μm). Collected by authors monthly chlorophyll-a data covered the vegetation period (May-October) from standard basins/monitoring area. Collected hourly to monthly Chl-a data were available in the National Water Quality Database, supplemented with Chl-a data from the ecological observatories (Istvánovics and Honti, 2021). <strong>References</strong> Báldi, K., C. Balogh, O. Sztanó, K. Buczkó, I. B. Muskó, L. G-Tóth, and Z. Serfőző. 2019. Sediment contributing invasive dreissenid species in a calcareous shallow lake – Possible implications for shortening life span of lakes by filling. <em>Elementa</em> <strong>10</strong>: 1-22. https://doi.org.10.1525/elementa.380 Balogh, C., and S. Purgel. 2012. A kvagga kagyló (<em>Dreissena bugensis</em>) térhódítása a Balatonban. [Rapid spread of quagga mussel (<em>Dreissena bugensis</em>) in Lake Balaton.] <em>Hidrologiai Közlöny</em> <strong>92(5–6)</strong>: 6–7 (in Hungarian) Balogh, C., I. B. Muskó, L. G.-Tóth, and L. Nagy. 2008. Quantitative trends of zebra mussels in Lake Balaton (Hungary) in 2003–2005 at different water levels. <em>Hydrobiologia</em> <strong>613</strong>: 57–69. https://doi.org/10.1007/s10750-008-9472-3 Balogh, C., A. Vláčilová, L. G.‐Tóth, and Z. Serfőző. 2018. <em>Dreissenid</em> colonization during the initial invasion of the quagga mussel in the largest Central European shallow lake, Lake Balaton, Hungary. <em>Journal of Great Lakes Research </em><strong>44</strong>: 114–125. https://doi.org/10.1016/j.jglr.2017.11.007 Balogh. C., Z. Serfőző, A. Bij de Vaate, R. Noordhuis, and J. Kobak. 2019. Differences in biometry, shell resistance and attachment strength of dreissenid mussels during the invasion of the quagga mussel (<em>Dreissena rostriformis bugensis</em>) in large European lakes. <em>Journal of Great Lakes Research </em><strong>45</strong>: 777-787. https://doi.org/10.1016/j.jglr.2019.05.011 Balushkina, E. V., and G. G. Winberg. 1979. The relationship between zooplankter weight and body length. <em>In:</em> General foundations for the study of aquatic ecosystems (pp. 169-172). Nauka, Leningrad, USSA (in Russian). Bottrell, H. H., A. Dunkan, Z. M. Gliwicz, E. Grygierek, A. Herzig, A. Hillbricht-Ilkowksa, H. Kurasawa, P. Larsson, and T. Weglenska. 1976. A review of some problems in zooplankton production studies. <em>Norwegian Journal of Zoology</em> <strong>24</strong>: 419–456. Burlakova, L. E. 1998. Ecology of <em>Dreissena polymorpha</em> (Pallas) and its role in the structure and function of aquatic ecosystems. Ph.D. Dissertation, Zoology Institute of the Academy of Science, Republic Belarus, Minsk, Belarus (in Russian). Burlakova. L. E., A. Y. Karatayev, and D. K. Padilla. 2006. Changes in the distribution and abundance of <em>Dreissena polymorpha </em>within lakes through time. <em>Hydrobiologia</em> <strong>517</strong>: 133–146. Fitzgerald, D. G., B. Zhu, L. G. Rudstam, S. B. Hoskins, D. E. Haddad, N. R. Burtch, J. T. Coleman, D. L. Crabtree, and E. L. Mills. 2016. Dynamics of aquatic vegetation in Oneida Lake, 1915-2005: A response to ecosystem change. <em>In:</em> L. G. Rudstam, E. L. Mills, J. R. Jackson, D. J. Stewart editors. Oneida Lake: Long-term dynamics of a managed ecosystem and its fisheries (pp. 181–199). American Fisheries Society, Bethesda, Maryland, USA. Herodek, S. 1986. Phytoplankton changes during eutrophication and P and N metabolism. pp. 183-204. In: Modelling and Managing Shallow Lake Eutrophication. (Eds: L. Somlyódy, G. van Straten). Springer, Berlin. <strong>            </strong>Hetherington, A., L. G. Rudstam, R. L. Schneider, K. T. Holeck, C. W. Hotaling, J. E. Cooper, and J. R. Jackson. 2019. Invader invaded: population dynamics of zebra mussels <em>(Dreissena polymorpha</em>) and quagga mussels (<em>Dreissena rostriformis bugensis</em>) in polymictic Oneida Lake, NY, USA (1992–2013). <em>Biological Invasions</em> <strong>21</strong>: 1529–1544. Hillebrand, H., C. Dürselen, D. Kirschtel, U. Pollingher, and T. Zohary. 1999. Biovolume calculation for pelagic and benthic microalgae. <em>Journal of Phycology</em>. <strong>35(2)</strong>: 403–424. DOI: 10.1046/j.1529-8817.1999.3520403.x. Horgan, M. J., and E. L. Mills. 1999. Zebra mussel filter feeding and food-limited production of <em>Daphnia</em>: recent changes in lower trophic level dynamics of Oneida Lake, New York, USA. <em>Hydrobiologia</em> <strong>411</strong>: 79-88. Idrisi, N., E. L. Mills, L. G. Rudstam, and D. J. Stewart. 2001. Impact of zebra mussels, <em>Dreissena polymorpha</em>, on the pelagic lower trophic levels of Oneida Lake, New York<em>. Canadian Journal of Fisheries and Aquatic Sciences</em> <strong>58</strong>:1430-1441. Istvánovics, V., A. Clement, L. Somlyódy, A. Specziár, L. G.-Tóth, and J. Padisák. 2007. Updating water quality targets for shallow Lake Balaton (Hungary), recovering from eutrophication. <em>Hydrobiologia</em> <strong>581</strong>: 305–318. https://doi.org/10.1007/s10750-006-0509-1 Istvánovics, V., and M. Honti. 2021. Stochastic simulation of phytoplankton biomass using eighteen years of daily data-Predictability of phytoplankton growth in a large, shallow lake. Sci. <em>Total Environment</em> <strong>764</strong>: 143636. https://doi.org/10.1016/j.scitotenv.2020.143636 Karatayev, A. Y. 1983. Ecology of <em>Dreissena polymorpha</em> Pallas and its effects on macrozoobenthos of the thermal power plant's cooling reservoir.  Candidate Dissertation, Zoology Institute of Academy of Science Belorussian SSR, Minsk, Belarus (in Russian). Karatayev, A. Y., L. E. Burlakova, and D. K. Padilla. 1997. The effects of <em>Dreissena polymorpha</em> (Pallas) invasion on aquatic communities in eastern Europe. <em>Journal of Shellfish Research</em> <strong>16</strong>: 187–203. Kovalevskaya R. Z., H. A. Zhukava, and B. V. Adamovich. 2020. Modification of the method of spectrophotometric determination of chlorophyll a in the suspended matter of water bodies. <em>Journal of Applied Spectroscopy</em> <strong>87(1)</strong>: 72–78. DOI 10.1007/s10812-020-00965-9 Lyakhnovich, V. P., S. I. Gavrilov, A. Y. Karatayev, I. V. Karatayeva, and T. I. Nekhaeva. 1982. Long-term changes in macrozoobenthos of Lukomskoe lake.  <em>Vestsi Akademii navuk Belaruskai SSR. Seryia bialahichnykh navuk </em><strong>1</strong>: 91–93 (in Belorussian). Lyakhnovich, V. P., A. Y. Karatayev, P. A. Mitrakhovich, L.V. Guryanova, and G. G. Vezhnovets. 1988. Productivity and prospects for utilizing the ecosystem of Lake Lukoml, thermoelectric station cooling reservoir.  <em>Soviet Journal of Ecology</em> <strong>18</strong>: 255–259. Mastitsky S. E., Y. K. Veres, O. A. Nayarovich, and S. Y. Kondobarov. 2006. Role of zebra mussel (<em>Dreissena polymorpha</em>) in structure of malacological complex of Lake Naroch. <em>In: </em>A. E. Kundas, A. E. Okeanova, and C. C. Poznyak editors. Materials of the 6th International scientific conference “Sakharovskie Readings 2006: ecological problems of the 21st century”. Part 1, (pp. 322–324). Minsk, Belarus (in Russian). Mate, F. 1987. A Balaton-meder recens iiledekeinek terkepezese (Mapping of modern Lake Balaton bottom sediments). <em>In:</em> Magyar Allami Foldtani Intezet Evi Jelentese az 1985. Evrol (Annual Report of the Hungarian Geological Institute of 1985). Hungarian Geological Institute, Budapest, pp. 367-379. Mayer, C. M., R. A. Keats, L. G. Rudstam, and E. L. Mills. 2002. Scale-dependent effects of zebra mussels on benthic invertebrates in a large eutrophic lake. <em>Journal of the North American Benthological Society </em><strong>21</strong>:616-633. Mellina, E., J. B. Rasmussen, and E. L. Mills. 1995. Impact of the zebra mussel (<em>Dreissena polymorpha</em>) on phosphorus cycling and chlorophyll in lakes. <em>Canadian Journal of Fisheries and Aquatic Sciences</em> <strong>52</strong>: 2553-2573. Mikheyeva, Т. М. 1989. Methods of quantitative accounting of nanophytoplankton (review). <em>Hydrobiological Journal. [Gidrobiologicheskiy zhurnal]</em> <strong>25(4)</strong>: 3–21 (in Russian). Mikheyeva, Т. М. 1999. Algal flora of Belarus. Taxonomic Catalogue. BSU Publishing, Minsk, Belarus (in Russian). Mitrachovich, P. A., V. M. Samoilenko, Z. K. Kartashevich, A. A. Svirid, E. A. Kozlov, G. N. Koolev, and N. A. Papko. 2008. Ecosystem of Lukoml thermoelectric station cooling reservoir. Pravo and Economica Press, Minsk, Belarus. Müller. G., and F. Wagner. 1978. Holocene carbonate evolution in Lake Balaton (Hungary): a response to climate and impact of man. – In: Matter, A., Tucker, M. E. (ed.): Modern and ancient lake sediments. Special publications of the International Association of Sedimentologists. Blackwell Scientific Publications, 57–81. https://doi.org/10.1002/9781444303698.ch4 Muskó, I. B., and B. Bakó. 2005. The density and biomass of <em>Dreissena polymorpha</em> living on submerged macrophytes in Lake Balaton (Hungary). <em>Archiv für Hydrobiologie</em> <strong>162(2)</strong>: 229-251. https://doi.org/10.1127/0003-9136/2005/0162-0229 Noordhuis, R., B. G. van Zuidam, E. T. H. M. Peeters, and G. J. van Geestet. 2016. Further improvements in water quality of the Dutch Border lakes: two types of clear states at different nutrient Levels. <em>Aquatic Ecology</em> <strong>50</strong>: 521–539. Ostapenya, A. P, A. A. Kovalev, T. V. Zhukova, T. M. Mikheeva, N. M. Kryuchkova, V. A. Babitsky, R. Z. Kovalevskaya, S. B. Kostyukevich, G. A. Inkina, T. A. Makarevich, E. P. Zhukov, V. F. Ikonnikov, A. M. Samusenko, A. F. Orlovsky, A. N. Rachevsky, and O. F. Jakushko, 1993. Ecological Passport of Naroch Lake. Ekomir Press, Minsk, Belarus (in Russian). Ostapenya, A. P, A. A. Kovalev, T. M. Mikheeva, V. A. Babitsky, T. V. Zhukova, N. M. Kryuchkova, R. Z. Kovalevskaya, S. B. Kostyukevich, G. A. Inkina, T. A. Makarevich, E. P. Zhukov, V. F. Ikonnikov, A. M. Samusenko, A. F. Orlovsky, A. N. Rachevsky, and O. F. Jakushko, 1994. Ecological Passport of Myastro Lake. Ekomir Press, Minsk, Belarus (in Russian). Ostapenya, A. P, T. V. Zhukova, T. M. Mikheyeva, R. Z. Kovalevskaya, T. A. Makarevich, H. A. Zhukava, E. V. Lukyanova, L. V. Nikitina, O. A. Makarevich, N. V. Dubko, V. S. Karabanovich, I. V. Savich, and Yu. K. Veras. 2012. Bentification of the lake ecosystem: causes, mechanisms, possible consequences, research prospects. <em>Proceedings</em> <em>of the</em> <em>Belarusian State University</em>  <strong>7(1):</strong> 135–148 (in Russian). Padisák. J., and C. S. Reynolds. 1998. Selection of phytoplankton associations in Lake Balaton, Hungary, in response to eutrophication and restoration measures, with special reference to the cyanoprokaryotes. <em>Hydrobiologia</em> <strong>384</strong>: 41-53. https://doi.org/10.1023/A:1003255529403 Pósfai, M., Z. Molnár, P. Pekker, I. Dódony, S. Frisia, and P. Meister. 2019. Microstructure of magnesium-bearing carbonates precipitating from shallow freshwater lakes. – GSA Annual Meeting in Phoenix, Arizona, 335892. https://doi.org/10.1130/abs/ 2019am-335892 Rudstam, L. G., E. L. Mills, J. R. Jackson, and D. J. Stewart. 2016. An introduction to the Oneida Lake research program and data sets. <em>In:</em> L. G. Rudstam, E. L. Mills, J. R. Jackson, D. J. Stewart editors. Oneida Lake: Long-term dynamics of a managed ecosystem and its fisheries (pp. 3–12). American Fisheries Society, Bethesda, Maryland, USA. Rudstam, L. G. 2020a. Lakewide zebra and quagga mussel summary, Oneida Lake, New York, 1992 to present. Web Data on Knowledge Network for Biocomplexity. Available: https://knb.ecoinformatics.org/view/kgordon.23.73, Original data set published in 2008.  Updated through 2019. Rudstam, L.G., 2020b. Limnological data and depth profile from Oneida Lake, New York, 1975-present. Web Data on Knowledge Network for Biocomplexity. Available http://knb.ecoinformatics.org/#view/kgordon.35.106, Original data set published in 2008.  Updated through 2019. Rudstam, L. G. 2020c. Zooplankton survey of Oneida Lake, New York, 1964 to present. Web Data on Knowledge Network for Biocomplexity. Available:<strong> </strong>https://knb.ecoinformatics.org/view/kgordon.17.89: Original data set published in 2008.  Updated through 2019. Rudstam, L. G. 2020d. Benthic invertebrates in Oneida Lake, New York, 1956-present. Web Data on Knowledge Network for Biocomplexity Available https://knb.ecoinformatics.org/#view/kgordon.4.75: Original data published in 2008. Updated through 2019. Rudstam, L. G. 2020e. Phytoplankton survey in Oneida Lake, New York, 1975 -present. Web Data on Knowledge Network for Biocomplexity. Available: https://knb.ecoinformatics.org/view/kgordon.31.87:Original data set published in 2008.  Updated with year 2018. SCOR-UNESCO Working Group no. 17. Determination of Photosynthetic Pigments in Sea-Water, Monographs on Oceanologic Methodology (1996). UNESSCO. Paris, 9–18. Sebestyen, O.  1937.  Colonization of two new fauna-elements of Pontus-origin (Dreissena polymorpha Pall. and Corophium curvispinum G. O. Sars forma devium Wundsch) in Lake Balaton.  Verhandlungen der Internationalen Vereinigung für Theoretische und Angewandte Limnologie 8: 169-182. Sebestyén, V., J. Németh, T. Juzsakova, E. Domokos, Z. Kovács, and  Á. Rédey. 2017. Aquatic environmental assessment of Lake Balaton in the light of physical-chemical water parameters. <em>Environmental Science and Pollution Resierch</em> <strong>24</strong>: 25355–25371. https://doi.org/10.1007/s11356-017-0163-3 Strickland, J. D. H., and T. R. Parsons. 1972. A practical handbook of seawater analysis. <em>Fisheries Research Board of Canada, Bulletin</em> <strong>167</strong>, 2nd edition. Watkins, J. M., L. G. Rudstam, and K. T. Holeck. 2011. Length-weight regressions for zooplankton biomass calculations – A review and a suggestion for standard equations. eCommons Cornell http://hdl.handle.net/1813/24566. Wellburn, A. R. 1994. The spectral determination of chlorophylls a and b, as well as total carotenoids, using various solvents with spectrophotometers of different resolution. <em>Journal of Plant Physiology</em> <strong>144</strong>: 307–313. https://doi.org/10.1016/S0176-1617(11)81192-2 Winberg, G. G. 1985. Ecological System of the Narochanskie Lakes.  Universitetskoe Press, Minsk, Belarus (in Russian). Zlinszky, A., G. Molnár, and B. Székely. 2010. Mapping bathymetry and thickness of lacustrine deposits of Lake Balaton (Hungary), using lake seismic profiles. <em>Földtani Közlöny</em> <strong>140</strong>: 429–438.
description: Please see details in the ReadMe file.