Ontogenetic shifts in growth and habitat use of freshwater Coastal Cutthroat Trout revealed through otolith analysis of multiple elements
Data files
Apr 13, 2026 version files 6.09 MB
-
adult_primordia_measurements.csv
7.54 KB
-
Comox_cutthroat_growth_data_annuli_Zinc_2025.csv
32.82 KB
-
Comox_lake_unadjustedformolarmass.csv
4.34 MB
-
Comox_OT_2025.R
50.39 KB
-
Comox_Sources.csv
3.51 KB
-
Comox_stream_unadjustedformolarmass.csv
1.64 MB
-
juvenile_primordia_measurements.csv
8.23 KB
-
README.md
10.14 KB
Abstract
This dataset contains otolith microchemistry transects and derived products for adfluvial Coastal Cutthroat Trout (Oncorhynchus clarkii clarkii) sampled in the Comox Lake–watershed system (Vancouver Island, BC). Raw transects for stream-caught juveniles and lake-caught adults were acquired as elemental concentrations (ppm) and converted in this workflow to molar ratios (element, mmol/mol) for Mg, Mn, Zn, Sr, Ba (Li included where available). Primordia and check positions (e.g., first winter) were compiled from microscope-based measurements for both adults and juveniles, enabling extraction of life-stage-specific windows (e.g., natal/rearing sections) along each transect.
We use multiple analytical tracks.
(1) Natal origin discrimination: After cleaning and unit conversion, we compute OTF-level means over a standardized rearing window and log-transform elements for MANOVA/ANOVA with Tukey HSD across capture sources; letters are overlaid on boxplots to summarize pairwise differences. We then fit linear discriminant functions (LDFA) on juvenile rearing means (Sr, Ba, Mn, Mg; Zn available in the data and plots) and apply the trained model to adult lake fish to estimate natal contributions, reporting jackknife confusion matrices, posterior probabilities, and assignment summaries.
(2) Movement/habitat signal detection: We assess Mn spikes as putative indicators of habitat shifts (e.g., stream→lake entrance) using a hybrid rule (rolling mean + individual mean fallback), summarizing spike count, magnitude, start/end locations, and before/after contrasts within fish.
(3) Zn oscillations as annulus proxies by linking Zn inflection distances to visually measured annulus starts (1–5AL) and fitting per-annulus simple linear models, reporting coefficients, r², model p-values, and slope-≈1 tests.
Together, the dataset and scripts support reproducible conversion, QC, visualization of transects, hypothesis testing across sources, discriminant assignment of adults, and mechanistic exploration of Mn and Zn features as indicators of habitat transitions and age structure.
We have submitted our raw elemental transect data for stream-caught juveniles (Comox_stream_unadjustedformolarmass.csv), raw elemental transect data for lake-caught adults (Comox_lake_unadjustedformolarmass.csv), adult otolith primordia and check position data (adult_primordia_measurements.csv), juvenile otolith primordia and check position data (juvenile_primordia_measurements.csv), capture metadata linking fish to source tributaries (Comox_Sources.csv), annuli growth measurements and zinc oscillation data (Comox_cutthroat_growth_data_annuli_Zinc_2025.xlsx), and R script (Comox_OT_2025.R).
Descriptions
Comox_stream_unadjustedformolarmass
Raw elemental transect data (ppm) for stream-caught juvenile Coastal Cutthroat Trout (Oncorhynchus clarkii clarkii) collected from natal tributaries.
- OT: Unique scan identifier
- Seconds: Ablation time position along transect
- Mg: Magnesium concentration (ppm)
- Mn: Manganese concentration (ppm)
- Zn: Zinc concentration (ppm)
- Sr: Strontium concentration (ppm)
- Ba: Barium concentration (ppm)
Comox_lake_unadjustedformolarmass
Raw elemental transect data (ppm) for lake-caught Coastal Cutthroat Trout collected from Comox Lake via gillnet. Values recorded as "UD" indicate concentrations that were Under Detection limit (i.e., below the instrument's detection threshold for that element at that measurement point).
- OT: Unique scan identifier
- Seconds: Ablation time position along transect
- Li: Lithium concentration (ppm); UD indicates the value was below the detection limit
- Mg: Magnesium concentration (ppm)
- Mn: Manganese concentration (ppm); UD indicates the value was below the detection limit
- Zn: Zinc concentration (ppm); UD indicates the value was below the detection limit
- Sr: Strontium concentration (ppm)
- Ba: Barium concentration (ppm)
adult_primordia_measurements
Otolith primordia and growth check position data (microns) for lake-caught adult fish. Source codes refer to gill net set identifiers used during lake sampling.
- Source: Gill net set identifier (see Source Codes below)
- OTF: Unique fish identifier
- Primordia: Micron position of the otolith core (primordium), determined visually
- Img_reversed: Whether the image orientation was reversed during measurement
- Exogenous_Feeding: Micron position of the exogenous feeding check
- First_Winter: Micron position of the first winter annulus check
juvenile_primordia_measurements
Otolith primordia and hatch check position data (microns) for stream-caught juvenile fish. Source codes follow a year.site.species.number convention (e.g., 2018.UC.CT.049) or site-number convention (e.g., UP-1).
- Source: Sample identifier combining site code and fish number (see Source Codes below)
- OTF: Unique fish identifier
- Primordia: Micron position of the otolith core (primordium), determined visually
- Img_reversed: Whether the image orientation was reversed during measurement
- Hatch_Check: Micron position of the hatch check
Comox_Sources
Capture metadata linking each fish (OTF) to its source tributary or lake capture location.
- Source: Tributary or capture location code (see Source Codes below)
- OTF: Unique fish identifier
Comox_cutthroat_growth_data_annuli_Zinc_2025
Annuli growth measurements and zinc oscillation data for lake-caught fish (xlsx format, sheet: ABLDATA). Contains otolith dimensions, annulus positions, zinc inflection points, manganese and barium spike locations, back-calculated fork lengths, posterior assignment probabilities from discriminant function analysis, and associated capture metadata.
- Otolith #: Unique fish identifier (corresponds to OTF)
- Gill Net ID: Gill net set identifier (e.g., GN1, GN2, GNUN)
- FL (mm): Fork length at capture (mm)
- Source: Assigned natal tributary group from DFA (e.g., Lake, UP, UCR, RECO)
- Stream: Assigned natal tributary group
- Posterior: Posterior probability of natal assignment from discriminant function analysis
- Center of Primordia (µm): Micron position of the otolith core
- OTW: Otolith total width (µm)
- OR: Otolith radius (µm)
- 1AL–5AL: Annulus locations (µm from primordium) for ages 1–5
- 1AW–5AW: Annulus widths (µm) for ages 1–5
- 1–5: Zinc inflection points corresponding to annuli (µm)
- 1D–5D: Distance between consecutive zinc inflection points (µm)
- Mangenese Spike Location: Micron position of the manganese spike indicating lake entry
- Barium Spike Location: Micron position of the barium spike
- OTW@LEM: Otolith width at the lake entry marker (µm)
- OTR@LEM: Otolith radius at the lake entry marker (µm)
- Regression (FL mm): Back-calculated fork length at lake entry using linear regression
- Fraser Lee: Back-calculated fork length at lake entry using the Fraser-Lee method
Comox_OT_2025
Full data-processing workflow including:
- Data import and cleanup (ppm → mmol/mol conversion relative to Ca)
- Reformatting and merging of lake + stream datasets
- MANOVA and ANOVA (Sr, Ba, Mg, Mn, Zn)
- Tukey HSD group comparisons and boxplots
- Linear Discriminant Function Analysis (LDFA) for natal tributary assignment
- Mn-spike detection workflow with rolling mean thresholds for lake entry inference
- Zn oscillation analysis for age-at-lake-entry estimation
- Plot generation and CSV exports
Source Codes
Natal Tributary Codes (Comox_Sources and juvenile files)
- BE: Beech Creek
- CO: Comox Creek
- CR: Upper Cruickshank Creek
- CRCO: Cruickshank River mainstem (near Comox Creek confluence)
- CRRE: Cruickshank River mainstem (near Reese Creek confluence)
- ER: Eric Creek
- PE: Pearce Creek
- RE: Reese Creek
- UP: Upper Puntledge River
- WL: Willemar Lake (excluded from final analysis)
Gill Net Identifiers (adult primordia and xlsx files)
- GN1–GN11: Gill net set identifiers used during the 2019 summer lake sampling campaign; numbers correspond to individual net positions in Comox Lake
- GNUN: Gill net unknown — fish captured as bycatch during fall 2018 in a separate study; associated capture metadata (e.g., fork length, weight, net location) were unavailable
DFA Group Codes (xlsx file)
- UCR: Upper Cruickshank + Eric Creek (combined group)
- RECO: Reese Creek + Comox Creek (combined group)
- UP: Upper Puntledge River (retained as a standalone group)
Juvenile Primordia Source Format
Juvenile primordia source codes follow two conventions:
- Year.Site.Species.Number (e.g., 2018.UC.CT.049), where site codes are: UC = Upper Cruickshank, EC = Eric Creek, RC = Reese Creek, CK = Comox Creek, UP = Upper Puntledge, BC = Beech Creek, PC = Pearce Creek
- Site-Number (e.g., UP-1, CO-3, CR(943)-1, UPCR-5), where UPCR = Upper Puntledge / Cruickshank River junction
Data Collection Details
- Collection Period: 2017–2019 (natal tributary juveniles collected fall 2017 and fall 2018; lake-based adults collected fall 2018 and summer 2019)
- Location: Comox Lake and associated natal tributaries, Vancouver Island, British Columbia, Canada
- Species: Coastal Cutthroat Trout (Oncorhynchus clarkii clarkii)
- Life Stages: 0-age juveniles (stream-caught from natal tributaries) and mixed-age adults (lake-caught via gillnet)
- Instrument: Laser Ablation Inductively Coupled Plasma Mass Spectrometer (LA-ICP-MS), School of Earth and Ocean Sciences, University of Victoria
- Units: All raw concentrations in ppm; converted to mmol/mol relative to Ca using molecular mass ratios (assuming 40% Ca by weight at 400,000 ppm). Values recorded as "UD" were below the instrument's detection limit.
Key Information Sources
- Wells, B.K. et al. (2003). Stock identification of Chinook salmon using otolith microchemistry.
- Olley, R., et al. (2011). Otolith microchemistry applications in biogeochemical tracer studies.
- Limburg, K.E. & Elfman, M. (2010). Patterns and magnitudes of Zn:Ca in otoliths.
- Lough, J.M. et al. (2014). Cutthroat Trout Stock Assessment Report for Comox Lake, BC. BC Ministry of Environment.
- Bottlenecks to Survival Project: British Columbia Conservation Foundation (BCCF), 2018–2024.
Code/Software
R is required to run Comox_OT_2025.R; the script was created using R version 2024.04.1, "Chocolate Cosmos" (R Development Core Team 2024). Annotations are provided throughout the script through 1) library loading, 2) data import and cleanup, 3) reformatting and merging, 4) analyses (MANOVA, ANOVA, Tukey HSD, LDFA), 5) Mn-spike detection, 6) Zn oscillation and annuli analysis, and 7) plot generation and CSV exports. Required R packages: ggplot2, reshape2, grid, gridExtra, lattice, plyr, dplyr, broom, MASS, zoo, multcomp, Hmisc, ellipse.
Accessibility and Versioning
- Authors: Jamieson B. Atkinson (ORCID: 0000-0002-6272-7761), Brendan J. Anderson, James P. Losee
- Affiliations: Coastal Ecosystems Research Collective, Nanaimo, BC, Canada; BC Ministry of Water, Land and Resource Stewardship, Nanaimo, BC, Canada; Wild Salmon Center, Portland, OR, USA
- Contact: Jamieson Atkinson — jatkinson@coastal-research.com
- Citation: Atkinson, J.B., Anderson, B.J., Losee, J.P. (2025). Comox Lake–Comox River Cutthroat Trout Otolith Microchemistry Dataset (2017–2019). Dryad Digital Repository. https://10.5061/dryad.w9ghx3g42
- Version: 1.0 (November 2025)
- Repository: Dryad Digital Repository (https://datadryad.org)
- Persistent Identifier (DOI): https://10.5061/dryad.w9ghx3g42
- License: CC0 1.0 Universal (CC0) public domain dedication, as required by Dryad
- Acknowledgements: Funding was provided by the Habitat Conservation Trust Foundation, the Provincial Fisheries program of BC, and the British Columbia Conservation Foundation.
