Impacts of zebra mussel (Dreissena polymorpha) invasion in an inland lake with suboptimal water chemistry
Data files
Apr 06, 2026 version files 189.88 KB
-
calcium.csv
528 B
-
Comartin_and_Ricciardi_Script.R
25.07 KB
-
mussels.csv
153.04 KB
-
README.md
4.66 KB
-
u_density_mortality.csv
5.84 KB
-
zm_size.csv
730 B
Abstract
Zebra mussel (Dreissena polymorpha) biofouling has caused rapid declines of native unionid mussel populations in North American inland waters with calcium [Ca2+] concentrations at or above 25 mg/L. We surveyed zebra mussel fouling and unionid mortality across sites and time in a large invaded lake (Lake Memphrémagog, Québec) and its connected waterbodies. The lake has a mean calcium concentration of 20.3 (range 18.6–23.2) mg/L, considered suboptimal for supporting a zebra mussel population large enough to exert lethal impacts on native mussels. Nevertheless, fouling intensities in the lake are of the same magnitude as those associated with rapid unionid declines elsewhere. Substantial fouling and unionid mortality in the lake began ~7 years after invasion, longer than the typical 3–5 year pattern recorded in calcium-rich waterbodies. We hypothesize that suboptimal water chemistry delayed the onset of lethal fouling but is insufficient to provide a refuge for unionid populations within the lake. Our results demonstrate that risk assessment must give more nuanced consideration to how the biodiversity impacts of zebra mussel invasion vary with water chemistry.
Dataset DOI: 10.5061/dryad.866t1g22v
The files provided are:
1) calcium.csv = Mean calcium concentration for every site in the different bodies of water sampled. Includes standard deviation, maximum and minimum calcium concentrations across sites.
2) mussels.csv = Sampling data for unionid mussels across year, season, body of water, and site. Includes data on quadrat number, unionid number, species ID, status, unionid mass, unionid length, zebra mussel fouling intensity, total unionid + zebra mussel mass, zebra mussel mass, and the zebra mussel/unionid mass ratio.U_density_mortality = Data on unionid density and mortality across year, season, body of water, and site. Includes data on the quadrat number.
3) u_density_mortality.csv = Unionid density and mortality data for every year, season, body of water, and site. Includes data on the quadrat number and the number of unionids that were found alive and recently dead.
4) zm_size.csv = Data on the size of zebra mussels that were found fouling on unionids in 2024 across sites and body of water. Includes the zebra mussel life stage and the zebra mussel fouling intensity.
5) Comartin_and_Ricciardi_Script.R = R Script used for data summary, analyses, and visualization using the .csv files above.
Descriptions
.csv files:
Year: Sampling year; 2022, 2023 or 2024.
Season: Season sampling occurred; Spring, Summer, or Fall.
BOW: Body of water where sampling took place;
ME = Lake Memphrémagog
MA = Lake Magog
MR = Magog River
Site: Site where sampling took place;
IL = Île Lord
SI = Sherbrooke Water Intake
BM = Baie de Magog
RM = Rue des Manoirs
RC = Rue Colby
BA = Baie de l'Anse
BC = Baie Cummins
PC = Pointe Cummins
IA = Île à l'Aigle
QK = Quai Knowlton
PF = Penfield
FB = Fitch Bay
CV = Cedarville
IR = Île Ronde
DB = Descente Bournival
PD = Plage Deauville
CO = Rue des Colombes
CN = Club Nautique
TG = Halte Turgeon-Gaudrault
Mean_Ca: Mean calcium concentration calculated from all calcium concentrations measured at that site (in mg/L).
SD_Ca: Standard deviation for the mean calcium concentration calculation.
Max_SD: Maximum calcium concentration that was measured at that site (in mg/L).
Min_SD: Minimum calcium concentration that was measured at that site (in mg/L)
Q_no: Number of the quadrat deployed (ranges between 1 to 10).
U_no: Number assigned to a unionid collected from a given site.
Species ID: Identification code for the species of unionid collected;
ELCO = Elliptio complanata
LARA = Lampsilis radiata
PYCA = Pyganodon cataracta
PYGA = Pyganodon grandis
Status: Status of the unionid;
A = Alive
RD = Recently dead
OD = Old dead
U_mass: Unionid mass (in g).
U_length: Unionid length (in cm).
ZM: Number of zebra mussels found fouled onto a given unionid (fouling intensity).
U_ZM_mass: Mass of the unionid + mass of attached zebra mussels (in g).
ZM_mass: Zebra mussel mass (in g).
Mass ratio: Mass of attached zebra mussels (g) / mass of unionid host (g).
U_density: Unionid density (individuals/m2).
U_alive: Number of alive unionids.
U_RD: Number of recently dead unionids.
Life_stage: Zebra mussel life stage; Adult or Juvenile
Missing values:
In mussels.csv, quadrats where unionids were absent during sampling are marked NA across the row.
In u_density_mortality.csv where unionids were absent (and unionid density is 0), NA is marked for U_alive and U_RD.
Code/Software
Comartin_and_Ricciardi_Script.R:
R is required to run Comartin_and_Ricciardi_Script.R The script was created using version 4.4.2 and uses the following packages: dyplr (v. 1.2.0), tidyr (v. 1.3.1), ggplot2 (v. 3.5.2) lme4 (v. 1.1.37), lmerTest (v. 3.1.3) glmmTMB (v. 1.1.13), DHARMa (0.4.7), rstatix (0.7.2), scales (v. 1.4.0), and ggbreak (v. 0.1.6). Annotations are provided throughout the script for data cleaning and summary, statistical tests and data visualization.
Note: JR_data.csv and Ricciardi_data.csv are used in the R Script but not provided here.
Data in the JR_data.csv file can be obtained by contacting the corresponding author: doi:10.1111/j.1365-2427.2008.02009.x Data in the Ricciardi_data.csv file can be obtained by contacting the corresponding author: doi:10.1046/j.1365-2427.2003.01071.x.
