Alternative migratory strategies related to life history differences in the Walleye (Sander vitreus)
Movement Ecology volume 10, Article number: 10 (2022)
While Pace of Life Syndrome predicts behavioural differences between individuals with differential growth and survival, testing these predictions in nature is challenging due to difficulties with measuring individual behaviour in the field. However, recent advances in acoustic telemetry technology have facilitated measurements of individual behaviour at scales not previously possible in aquatic ecosystems.
Using a Walleye (Sander vitreus) population inhabiting Black Bay, Lake Superior, we examine whether life history characteristics differ between more and less mobile individuals as predicted by Pace of Life Syndrome. We tracked the movement of 192 individuals from 2016 to 2019 using an acoustic telemetry study, relating patterns in annual migratory behaviour to individual growth, and seasonal changes in optimal thermal-optical habitat.
We observed two consistent movement patterns in our study population—migratory individuals left Black Bay during late summer to early fall before returning to the bay, whereas residents remained within the bay year-round. The average maximum length of migrant Walleye was 5.5 cm longer than residents, and the sex ratios of Walleye caught during fall surveys was increasingly female-biased towards the mouth of Black Bay, suggesting that a majority of migrants were females. Further, Walleye occupancy outside of Black Bay was positively associated with increasing thermal-optical habitat.
Walleye in Black Bay appear to conform to Pace of Life Syndrome, with migrant individuals gaining increased fitness through increased maximum size, which, given size-dependent fecundity in this species, likely results in greater reproductive success (via greater egg deposition vs. non-migrants). Further, apparent environmental (thermal) controls on migration suggest that migratory Walleye (more so than residents) may be more sensitive to changing environmental conditions (e.g., warming climate) than residents.
The theoretical concept of the Pace of Life Syndrome (POLS) provides a framework for understanding how animal behaviour relates to individual life history variation . Life history variability results in differences in age-structured productivity and mortality, with bold and highly mobile individuals expected to have greater food intake, growth rates, and reproductive output than more cautious individuals . These advantages of bold behaviour, however, come with trade-offs of higher energetic costs and predation risk [36, 37]. While ecosystems with scarce or abundant resources often drive evolution to highly mobile or sedentary strategies respectively, unstable or intermediately favourable conditions are well suited to behavioural variability in movement patterns . Under such conditions, the POLS predicts that individuals will vary in a suite of behavioural and physiological characteristics associated with life history [36, 37].
Continually shifting environmental conditions can lead to challenges in identifying behavioural patterns and individual personalities, highlighting the importance of long term study [31, 36]. Disparate movement behaviours have been observed in spatial ecology studies with links to other personality features [12, 17], and in some cases have identified links between movement patterns with life history . Documenting movement patterns in natural environments has historically proved challenging without intensive time and monetary investments, but advancements in technology are making such studies increasingly more feasible, particularly for fishes [11, 28].
The emerging field of passive acoustic telemetry, which makes use of hydrophone-receiver arrays to record locations of fish fitted with acoustic transmitters, provides detailed insight into movement patterns found within populations across spatial scales [23, 25]. Traditional space use studies using mark recapture techniques that physically re-capture animals are often highly resource intensive and cannot gather the same resolution of data as telemetry [22, 41]. This ability to passively and accurately track the movement of large numbers of fish allows researchers to compare individual movement patterns with life history outcomes such as growth, survival, and reproduction (e.g., ).
The Great Lakes Acoustic Telemetry Observation System (GLATOS) is a large scale acoustic telemetry network in the world’s largest freshwater lakes, the Laurentian Great Lakes, which allows researchers to address questions related to fish spatial ecology . Walleye (Sander vitreus), a piscivorous top predator whose productivity is tightly tied to water temperature and clarity (preferring turbid environments between 18 and 22 °C; [9, 27], are frequently studied using the GLATOS network due to their economic importance within the Great Lakes , and throughout its range .
Within the Great Lakes, Walleye have been observed undertaking long distance migrations , and in some instances females travel greater distances than males [14, 35]. Access to different habitats may underlie these different movement patterns, with males preferring warmer, shallower habitats than females . Further, energetics modeling indicates that lower feeding activity in male percids results in female-biased sexual size dimorphism . These observations suggest that larger female Walleye may undergo migration in order to seek out higher energy prey as a strategy to achieve larger size, and therefore, fecundity .
To determine whether differential behavioural strategies are associated with life history differences in Walleye, our study makes use of a 4-year acoustic telemetry program in Black Bay, Lake Superior (Fig. 1). Black Bay is shallower, warmer, more turbid, and supports a different fish community than the deep, clear, and cold surrounding regions of Lake Superior . Building on evidence elsewhere indicating differential movement patterns within Walleye populations [14, 18, 35], we aimed to first verify that this pattern of differential movement exists in Black Bay. Second, we sought to use this system of differential movement to test the POLS prediction that migratory Walleye (that leave Black Bay) have both greater maximum size (and therefore fecundity), and greater mortality than resident Walleye. Given that Walleye display sexual size dimorphism, as well as sex-based variation in mobility observed elsewhere [14, 35], we used netting data collected from Black Bay to test the hypothesis that migratory and resident Walleye life history strategies were related to sex (i.e., evidence of a greater proportion of females migrating than males). Finally, given previous research documenting the role of the environment in shaping Walleye habitat occupancy and production [9, 27], we made use of a thermal-optical habitat model to test the hypothesis that Walleye migration is driven by the amount of available optimal thermal-optical habitat within different regions of their range.
Fish capture and tagging
Two datasets were used in this study; Walleye movement, growth, and thermal-optical habitat use were assessed using individuals that were captured and acoustically tagged, while the spatial distributions of male and female Walleye as both juveniles and adults were assessed using data gathered in the Ontario Ministry of Natural Resources and Forestry (MNRF) Fall Walleye Index Netting (FWIN) program (see ).
Adult Walleye to be acoustically tagged were captured via trapnets, electrofishing, short set gill nets, and angling within and immediately outside of Black Bay, and along the lower Black Sturgeon River from 2016 to 2018 (Additional file 1: Table S1; Figure S1). Fish were measured (total length, mm), and the second dorsal spine was taken for age and growth assessment. Fish were intra-coelomically fitted with Vemco V16 (2016–2018: n = 180, random transmission intervals between 60 and 180 s) and V13 (2017: n = 8, random transmission intervals between 120 and 240 s; 2018: n = 4, random transmission intervals between 60 and 180 s) acoustic tags and released at their respective capture sites. Two experienced biologists tagged 192 Walleye in 2016 (n = 94), 2017 (n = 61), and 2018 (n = 37) following Canadian Council on Animal Care approved procedures (Lakehead University file # 1465777; Additional file 1: Supplemental Information).
The movement of acoustically tagged Walleye was tracked within and immediately outside of Black Bay using acoustic receivers. The shallower, relatively protected waters of Black Bay remain ice covered longer than the main portion of Lake Superior, with ice up typically beginning in late November and persisting until early May. The ice-free open water season of late May- early November was used to deploy, retrieve, and maintain acoustic receivers in the area. During this open water season in 2016, acoustic receivers (n = 32; Vemco VR2W, VR2AR; 69 kHz) with omnidirectional hydrophones were deployed within and outside of Black Bay to detect acoustically tagged Walleye. Receivers were suspended 0.5–1.5 m off bottom with an anchor and float system and arranged in 5 single row gates at significant ecological boundaries, with 800 m spacing between receivers (Fig. 1). These gates mark the boundary between the north and south end of Black Bay at Bent Island (BEI; n = 9), the mouth of Black Bay at George Point (GEP; n = 5), Edward Island outside of the mouth of Black Bay (EDI; n = 7), and the tips of the Sibley (SIP; n = 2) and Black Bay Peninsulas (BLP; n = 3). At each of six positions throughout the bay, individual receivers independent of gates were also deployed (hereafter referred to as sentinel receivers) to provide assessment of Walleye use at specific points of interest; two sentinel receivers were placed at the mouths of the Black Sturgeon and Wolf Rivers, where Walleye are locally known to congregate each spring, and the other sentinel receivers were placed at other points of suspected Walleye congregation throughout the bay.
During winter (November–May), receivers in water shallower than 5 m were removed to prevent damage, leaving 24 receivers. All gate receivers were re-deployed for the open water season of 2017 (May–November), as well as a grid-work of receivers in the north end of Black Bay (n = 13, 5 km spacing), and sentinel positions (n = 4; Fig. 1). All gates were redeployed again in 2018 with lower receiver coverage (BEI n = 5, GEP n = 3, EDI n = 4, SIP n = 2, BLP n = 2). In 2019, only the peninsular gates (SIP, BLP) were redeployed, each with one receiver (used for evaluating long-range movement). Data was downloaded each spring and fall, and detection efficiency was assessed at the BEI, GEP, and EDI receiver gates by simulating receiver line performance based on range testing carried out on the receivers in the study area (Additional file 1 Supplemental Information).
Movement and survival
Walleye movement and survival was assessed from May 2016 to October 2017 using a multistate mark-resight model [6, 21] (Additional file 1: Supplemental Information) analyzed in Program MARK . These years were chosen from the larger acoustic telemetry study as they preceded changes in gate deployment for other studies focusing on the inner bay not described here (Fig. 1). States that could be occupied by Walleye were designated as the areas between receiver gates (Fig. 1), and during each sighting occasion, every fish was assigned an occupied state or indicated as not sighted. The occupied state assigned to each fish for a given period was determined using a weighted average of detections of each fish at all gates and sentinel receivers at the north end of Black Bay. A transition to a new state occurred when the weighted average of detections across the gates and northern sentinel receivers fell within a new state during the next sighting occasion. This produced a condensed encounter history for each acoustically tagged Walleye describing movement throughout the study area.
To determine the appropriate time scale for evaluating patterns, movement and survival were assessed with sighting occasions at monthly and bi-weekly time intervals. Candidate models were constructed and run in Program MARK, and top explanatory models for both time intervals were selected on the basis of AICc, where models within ΔAICc of 2 of the top model were considered to be equivalent . The top models were used to determine time- and state-dependence on transition, survival, and sighting probabilities. Parameter estimates from the top candidate models were used to assess the probability of transitioning between states and the proportion of surviving fish in each state during each occasion.
Repeatability in individual Walleye movement patterns was assessed using the acoustic telemetry deployment from 2016 to 2018 (Fig. 1). Walleye tagged in 2016 and 2017 were assigned as either ‘migrant’ or ‘resident’ for each year that they were detected during the period when receiver gates were present (2016–2018) based on the maximum outbound extent from the northern end of Black Bay; Walleye detected beyond the George Point (GEP) receiver line (mouth of Black Bay) were considered migratory for a given year, while Walleye whose maximum outbound detection was at, or within the George Point receiver line were considered resident. In order to evaluate individual movement pattern repeatability across all 3 years of observation, generalized linear mixed models (GLMMs) framework fit with a logit link was applied. The full model we evaluated considered annual observations of Walleye migration status (stay, leave) as a function of fish length at capture (as a fixed effect) and fish ID as a random effect. Using log-likelihood ratio tests, we evaluated the significance of both the fixed and random effect comparing the top model to reduced models excluding each of these terms individually. The repeatability of individual movement patterns across years was then assessed by calculating the intra-class correlation coefficient (R;  on the model including significant variables only (see results).
Differences in age and total length between Walleye captured in 2016 with clearly defined ageing structures (n = 53) and identified as either migratory or resident were evaluated with Welch two sample t tests.
Growth patterns of Walleye defined as migratory or resident in both 2016 and 2017 were determined from dorsal spine annuli by back calculating length at age using the Fraser-Lee method . Only spines with clearly defined annuli were used in this analysis (n = 53). Dorsal spine annuli measurements were taken from the focus to the edge of each annulus along the horizontal elongated transect using Image J . Von-Bertalanffy growth curves were fitted to both migratory and resident Walleye back calculations of length (Additional file 1: Supplemental Information). For both migratory and resident groups of Walleye, error around parameter estimates were generated using bootstrapping, drawing each size at age estimate within a group independently with replacement. Parameter estimates and associated confidence intervals were compared to determine if significant differences existed between the growth patterns of resident and migratory Black Bay Walleye.
Optimal Walleye habitat area was defined as the available optimal Thermal-Optical Habitat Area (TOHA), or the benthic area where temperature and light conditions were optimal for Walleye productivity . For three states within the study system (north of Bent Island, between Bent Island and George Point, and the region beyond Edward Island) monthly Secchi disk readings were taken during the 2017 open water season. Surface illuminance for 2017 was collected from the IISD Experimental Lakes Area (IISD-ELA) near Kenora, Ontario using a Kipp and Zonen SP Lite Sensor. The IISD-ELA site shares similar latitude and climate patterns to Black Bay , and no closer illuminance data were available. The hourly maximum and minimum depths at which optimal light conditions existed for Walleye in Black Bay were calculated as in  (Additional file 1: Supplemental Information). Hourly maximum and minimum optical depth preferences were averaged for each monthly period corresponding to the multistate mark-resight model to determine the depth range providing Walleye with preferred light intensities. Vertical temperature profiles were also created for the northernmost state (max depth = 14 m), the state between Bent Island and George Point (max depth = 36 m), and the southernmost state (max depth = 70 m) using temperature loggers deployed in these regions at 1 m intervals from 2 to 20 m depth, 5 m intervals from 20 to 40 m depth, and a logger at 50 m depth. Monthly temperature averages at each depth were used to determine the depth range at which conditions were optimal for Walleye (18–22 °C) [9, 27]. Depths of optimal thermal and optical habitat were then used to calculate TOHA for the three states during each monthly sighting occasion by calculating the bottom area that fell within the optimal depth ranges using ArcGIS 10.5 (Environmental Systems Research Institute, Redlands, California) using a 30 m by 30 m cell raster digital bathymetry model. The state between George Point and Edward Island was excluded from this analysis due to lost equipment, but the three states used here cover the range of habitats available to Black Bay Walleye.
Walleye habitat occupancy for each monthly sighting occasion was determined by the number of tagged fish assigned to each state in the multi-state mark-resight model. To evaluate the role of TOHA in driving Walleye occupancy, linear regressions were applied to monthly occupancy estimates within individual states, as well as a linear regression across combined data for all three states.
FWIN program data were used to assess if spatial differences existed in the relative abundance of male and female Black Bay Walleye. The relative number of male and female Walleye captured in each state of the study area were compared using Pearson’s chi-squared tests for both juvenile (0–3 years old) and adult (4 + years old) Walleye captured in FWIN nets set during September and October of 2002 to 2017. Because only one juvenile Walleye was captured in the southernmost state in any of the FWIN netting programs, this state was removed from analysis in juvenile fish to avoid violating the chi-square assumption that expected frequency of at least 80% of cells be at least five . Where chi-square tests showed statistical significance, standardized residuals were calculated to assess the contribution of each cell to the significance of the model . Following significant chi-squared results, the contribution of each cell (corresponding to each state for each sex) to this significance was determined by calculating cell-standardized residuals. Residuals greater than 2 indicate that the cell value is significantly higher than expected by random distribution, while residuals less than -2 indicate the cell value is significantly less than expected .
Walleye tagged across all years (2016–2018) were similar in age and length (Additional file 1: Table S2). Of 192 Walleye fitted with acoustic transmitters between 2016 and 2018, 184 were detected at least once on the receiver array between 2016 and 2019. The annual proportion of fish detected at each receiver gate remained consistent throughout the study period (Table 1). 179 of the tagged Walleye were detected north of the Bent Island gate at least once. Walleye were detected within the north end of Black Bay throughout the tracking period, but most of the fish congregated in this area over winter with most final fall and initial spring detections (between which north end receivers were removed) occurring north of the Bent Island each year. Walleye detections on the George Point and Edward Island receiver gates began in June, with occasional June detections as far as the peninsular gates. Walleye detections on the peninsular receiver gates (SIP, BLP) increased in July, and detections on all receivers outside of Black Bay remained highest between August and October of each year (Additional file 1: Figure S2). Between 2016 and 2019, 86 Walleye were detected beyond the George Point receiver gate. Individual movement patterns were highly variable with some fish exhibiting frequent movement back and forth between gates, while others traveled in a more linear fashion (Additional file 1: Figure S3).
The top candidate models from the multistate mark-resight modelling for both monthly and bi-weekly sighting occasions were identical, indicating that monthly time bins captured sufficient resolution of movement in Black Bay (Additional file 1: Table S3). Models from both time periods indicated that Walleye movement was not time dependent (Additional file 1: Table S3). Transition probabilities from one state to another did, however, depend on the previous state (Fig. 2). The probability of acoustically tagged Walleye being detected was not state dependent, but varied with time, with detection probabilities remaining above 80% except for the months of October–February (Additional file 1: Figure S4). Finally, the monthly survival probability of acoustically tagged Walleye was 97.9% and was consistent across states and time (78% annual survival).
While almost all Walleye were detected in the north end of Black Bay during the winter months, two distinct sub-groups—migrators and residents—were identified within the population based on their outbound extent of migration. Log-likelihood tests indicated that both the random effect of individual ID and the fixed effect of fish length explained significantly more variation than the models from which they were absent (individual ID, p < 0.0001; fish length, p < 0.0001). Using this full model, we found significant repeatability of individual movement patterns between years, accounting for variation in fish size (R = 0.89, p < 0.0001). Migratory and resident Walleye did not differ in mean age (t = -0.03, df = 50.66, p = 0.97) however, total length of migratory Walleye at the time of tagging (mean = 629 ± 10.67 mm) was greater than that of residents (mean = 588 ± 9.21 mm; t = -2.91, df = 50.21, p = 0.005).
Asymptotic lengths (L∞) differed significantly between resident and migratory Walleye based on bootstrapped 95% confidence intervals of back-calculated size at age (Migratory: L∞ = 655 mm 95% CI 642–667 mm; Resident: L∞ = 600 mm 95% CI 591–610 mm; Fig. 3). However, the curvature parameter (Brody’s k) of residents and migrators overlapped significantly (Migratory: k = 0.37 95% CI 0.34–0.39; Resident: k = 0.40 95% CI 0.38–0.43).
Walleye increased their use of habitat outside of Black Bay when the amount of available TOHA outside of the bay increased (linear regression, F1,3 = 99.26, R2 = 0.96, p = 0.002; Fig. 4). Although turbidity increased from the southernmost to the northernmost states, Secchi readings did not vary much over time (apart from deep measurements indicating clear water at the end of summer, Additional file 1: Figure S5), indicating that changing water temperature was the main driver of Walleye optimal habitat availability. While Walleye occupancy did not significantly vary with available TOHA in the north end of the bay (F1,3 = 0.03, p = 0.87), or the region between Bent Island and George Point (F1,3 = 5.10, p = 0.11) separately, the occupancy values fall near to those predicted by the TOHA model from outside of Black Bay (Fig. 4). A linear regression of combined occupancy and TOHA data from all states with temperature data was significant (F1,13 = 98.54, R2 = 0.87, p < 0.0001), but violated the assumption of homoscedasticity that log-transformation did not correct.
Adult FWIN captured Walleye showed non-random spatial distribution based on sex (χ23 = 72.43, p < 0.0001), with standardized residuals indicating Walleye captured outside of Black Bay were mostly female and Walleye captured in the north end of Black Bay were mostly male (Table 2). Ratios of male:female juvenile Walleye were more even across states (χ22 = 1.65, p = 0.44; Table 2). Furthermore, a larger majority of juvenile Walleye were captured in the northernmost state compared to adults.
Using an acoustic telemetry network established in North America’s largest lake, we identified two distinct and repeatable Walleye movement strategies corresponding to different life history modalities. These movement strategies were highly conserved between years and split with almost 50% of the tagged population adopting either migratory or resident strategies. Consistent with the POLS hypothesis, migrant fishes achieved larger body sizes than residents. Migratory Black Bay Walleye not only achieve a larger maximum size (L∞) than residents but share a similar rate of approach to this larger maximum size (Brody’s k), indicating that migratory Walleye also have a faster growth rate than residents (i.e., necessary for a similar rate of approach to result in a larger maximum size). Larger body size in Walleye typically confers greater fecundity , with the largest or oldest females in the population having disproportionately high reproductive success . Thus, fitness advantages in this migratory Walleye population may be achieved by reaching larger body sizes more quickly, conferring greater fecundity.
However, our study did not support all life history patterns predicted by PLOS. Whereas POLS predicts lower survival in migrants, multistate mark-resight modelling indicated no differences in survivorship between migrant and resident Black Bay Walleye. Because Black Bay is a fish sanctuary north of Bent Island (where all fish overwinter) and is closed throughout the bay to commercial Walleye fishing , this removes a major source of adult mortality. Similar survivorship between groups should lead to an increased proportion of migratory Walleye in the Black Bay population due to the fitness advantage of increased size. This, perhaps, contributed to the Walleye population growth in the region observed in the FWIN netting program . While unstable environments create ideal conditions for behavioral variability predicted by POLS, long term changes to anthropogenic pressures may lead to a shift in fitness advantages of different behaviours, and eventual change in the proportion of these behaviours within a population. Humans have a long history of altering animal behaviour, including by modifying the real and perceived risks faced by animals in the face of the human “super predator” .
Changing climate conditions are also likely to have lasting impacts on the variability in behavioural strategies used by animals. Our findings indicate that thermal habitat conditions act as an environmental control on the timing of long-range Walleye migrations. Habitat use of migratory and resident Black Bay Walleye differed noticeably from August to October, when the water outside of Black Bay reached thermal and optical optima and Walleye occupancy outside the bay increased. Outside of late summer/early fall, water temperatures outside of Black Bay were below optimal levels for Walleye, and both migratory and resident fish shared the inner basin habitat. The apparent strong dependence of Walleye migration on thermal conditions suggests that major changes in weather or climate may greatly impact their degree of migration, such that in warmer years with longer summers (predicted to occur more frequently with climate change), migration opportunities into historically colder waters could increase. Behavioural variability expected by POLS may allow populations to adjust to changing environmental conditions however, extreme shifts associated with some climate change scenarios may result in the loss of a behavioural characteristic and lead to reduced population stability .
But what is the driver of annual migratory patterns in Black Bay Walleye? Given that migratory fish tend to grow larger and faster than residents, we speculate that migration is driven by access to larger-bodied, energetically-dense coregonids such as Cisco (Coregonus artedi), as opposed to the need to evade poor thermal-optical habitat conditions in Black Bay. Indeed, large regions inside Black Bay always fell within the optimal thermal-optical regime for Walleye between May and October, and we found no significant relationship between TOHA and Walleye occupancy within Black Bay. Access to larger, more energy dense prey has been shown elsewhere to provide increased growth efficiency and maximum size in Walleye . Optimized forage intake accessed through improved food quality and larger prey size can lead to a greater energetic surplus despite increased metabolic costs to migration [38, 42]. Interestingly, the period (Aug-Oct) where migratory Walleye leave Black Bay and are able to access these different food resources coincides with the period when female Walleye increase ovary development .
Male and female Black Bay Walleye showed systematic sex-based space use differences in detailed Walleye netting surveys. We found that adult female Walleye were more likely to be captured outside of Black Bay than males, whereas adult male Walleye were more likely to be captured inside Black Bay than females. The timing of the netting survey coincides with the period when migratory Walleye leave Black Bay, supporting the hypothesis that females are more common among migratory fish than males. Indeed, the difference in growth exhibited between migratory and resident Walleye closely resembles the sexual size dimorphism displayed in this species [19, 39].
Notably, immature Walleye (which lack sexual growth differences) did not show the same differential spatial use pattern as sexually mature Walleye, suggesting a greater exploratory behaviour in mature (female) fish only. In previous studies examining Walleye movement, females travelled greater distances than males to access cooler, deeper water [30, 35]. Behaviours resulting in trophic niche separation have been documented in a variety of populations and communities [8, 16, 40] and one plausible explanation for female migration is that the sex benefiting most from resource acquisition (females in the case of Walleye) is the one that migrates to find those resources. Our findings indicate that hypotheses surrounding sex-biased differences in life history strategy should be considered in future work on POLS.
The tendency of mature female Walleye to be captured outside Black Bay (with no similar pattern observed in juveniles), combined with patterns in growth differences strongly suggest that migrating Walleye in our study are predominantly female. Because immature male and female Walleye grow at similar rates before maturation , survival should only differ for mature fish. However, despite apparent sex-based differences adult Black Bay Walleye movement patterns, survival rates appear to be similar. This may be due to the lack of significant commercial harvest of Black Bay Walleye—which was historically a major source of mortality for this population until 2003 [2, 15]—providing an equalizing effect on death rates between male/female or resident/migrant individuals.
Distinct movement patterns of migration and residency corresponded to life history differences that were largely predicted by POLS. Spatial state transitions depended on the state to which fish were previously classified, but there was no time dependence of movement patterns based on multistate mark-resight modeling. This was despite the predictable observed patterns of out-migration observed between August and October, and the congregation in the north end of the bay each winter and spring. Space use by migratory Black Bay Walleye did not differ from that of residents for much of the year, with all Walleye congregating in the north end of the bay during the winter months and spring spawning. Within both migratory and resident Black Bay Walleye groups we found variation in movement patterns; some fish made direct movements between positions, while others made frequent back and forth trips. Given our telemetry array design and study objectives, we evaluated the maximum outbound extent of movement, but not the total movement (which in highly mobile resident individuals could potentially be greater than migrants). In winter, receivers were removed from shallow regions of Black Bay to avoid ice damage, contributing to the marked decline in resight probability during this time. Detection simulations along receiver lines (which were in place over winter) indicated that a fish crossing them would be detected 99% of the time based on 2016 and 2017 deployments, suggesting Walleye movement throughout the Bay was minimal during this period and that few Walleye left the north end of Black Bay during winter.
Future work using telemetry to address questions related to animal life history will benefit by exploring multiple spatial and temporal scales. As observed in this study, variation in individual movement may occur in annual migration patterns, or in the frequency of smaller movements on the scale of hours to days. Using telemetry arrays with greater spatial resolution may reveal more exploratory vs more cautious behavioral patterns within resident populations. Greater spatial resolution in spawning areas will also allow researchers to compare behavioural differences in spawning behavior or site selection to life history outcomes such as spawning success.
By combining acoustic telemetry and multistate mark-resight models with a back-calculation method of estimating fish growth, our study supports the POLS hypothesis that migratory individuals are likely to achieve larger body sizes at a given age and therefore achieve greater fitness than their resident counterparts. Further, we suspect that energetic rewards of females accessing more energy-dense prey drives this pattern of migrant fishes achieving larger maximum body size. While differences in survival were not observed, we believe that this is related to the recent removal of a major source of mortality in the form of commercial fishing. Finally, because migratory opportunities appear to be temperature-dependent, we speculate that conditions associated with climate warming may result in larger differentials in migrant vs. resident Walleye, and if these patterns are sex-dependent, increased differentials of sexual size dimorphism of Lake Superior Walleye.
Availability of data and materials
The acoustic telemetry dataset analysed during the current study are available upon request in the GLATOS (Great Lakes Acoustic Telemetry Observation System) repository, https://glatos.glos.us/home/project/BBWAT (Date range: May 2016 to November 2019). The FWIN netting data analyzed in this study is available upon request from the Upper Great Lakes Management Unit, Ontario Ministry of Natural Resources and Forestry, Thunder Bay, Ontario, P7E6S7. Light data analyzed in this study was provided by Ken Sanilands (International Institute for Sustainable Development Experimental Lakes Area, Winnipeg, Manitoba, Canada, R3B0Y4). Other datasets and code generated and/or analyzed in the current study have been included in the Supplementary Information (Additional files 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12).
Akaike Information Criterion
Bent Island (Receiver Gate)
Black Bay Peninsula (Receiver Gate)
Edward Island (Receiver Gate)
Fall Walleye Index Netting
George Point (Receiver Gate)
Great Lakes Acoustic Telemetry Observation System
International Institute for Sustainable Development-Experimental Lakes Area
Ontario Ministry of Natural Resources and Forestry
Pace of Life Syndrome Hypothesis
Sibley Peninsula (Receiver Gate)
Thermal-Optical Habitat Area
Baccante DA, Reid DM. Fecundity changes in two exploited walleye populations. N Am J Fish Manag. 1988;8(2):199–209. https://doi.org/10.1577/1548-8675(1988)008.
Berglund E. Assessment and monitoring of Black Bay, lake superior walleye using fall walleye index netting (FWIN) 2002–2014. Upper Great Lakes Management Unit, Ontario Ministry of Natural Resources and Forestry, Thunder Bay, Ontario, 37 pp; 2014.
Bewick V, Cheek L, Ball J. Statistics review 8: qualitative data—tests of association. Crit Care. 2004;8(1):46–53. https://doi.org/10.1186/cc2428.
Biro PA, Stamps JA. Are animal personality traits linked to life-history productivity? Trends Ecol Evol. 2008;23(7):361–8. https://doi.org/10.1016/j.tree.2008.04.003.
Borkholder BD, Edwards AJ. Comparing the use of dorsal fin spines with scales to back-calculate length-at-age estimates in walleyes. N Am J Fish Manag. 2001;21(4):935–42. https://doi.org/10.1577/1548-8675(2001)021%3c0935:CTUODF%3e2.0.CO;2.
Brownie C, Hines JE, Nichols JD, Pollock KH, Hestbeck JB. Capture-recapture studies for multiple strata including non-Markovian transitions. Biometrics. 1993;49(4):1173–87. https://doi.org/10.1007/S00442-01.
Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach (2nd ed). In: Ecological modelling, vol 172; 2002. https://doi.org/10.1016/j.ecolmodel.2003.11.004.
Cabrol J, Trombetta T, Amaudrut S, Aulanier F, Sage R, Tremblay R, Nosais C, Starr M, Plourde S. Trophic niche partitioning of dminant North-Atlantic krill species, Meganyctiphanes norvegica, Thysanoesssa inermis, and T. raschii. Limnol Oceanogr. 2018;64(1):165–81. https://doi.org/10.1002/lno.11027.
Chu C, Minns CK, Moore JE, Millard ES. Impact of oligotrophication, temperature, and water levels on walleye habitat in the Bay of Quinte, Lake Ontario. Trans Am Fish Soc. 2004;133(4):868–79. https://doi.org/10.1577/T03-095.1.
Columbo SJ, McKenney DW, Lawrence KM, Gray PA (2007) Climate change projections for Ontario: policymakers and planners. Applied Research and Development Branch, Ontario Ministry of Natural Resources and Forestry, Sault Ste. Marie, Ontario, 38pp.
Cooke SJ, Martins EG, Struthers DP, Gutowsky LFG, Power M, Doka SE, Dettmers JM, Crook DA, Lucas MC, Holbrook CM, Krueger CC. A moving target—incorporating knowledge of the spatial ecology of fish into the assessment and management of freshwater fish populations. Environ Monit Assess. 2016. https://doi.org/10.1007/s10661-016-5228-0.
Cote J, Clobert J, Brodin T, Fogarty S, Sih A. Personality-dependent dispersal: characterization, ontogeny and consequences for spatially structured populations. Philos Trans R Soc B Biol Sci. 2010;365(1560):4065–76. https://doi.org/10.1098/rstb.2010.0176.
Dhellemmes F, Finger J, Smukall MJ, Gruber SH, Guttridge TL, Laskowski KL, Krause J. Personality-driven life history trade-offs differ in two subpopulations of free-ranging predators. J Anim Ecol. 2021;90:260–72.
Fielder DG. Mortality, exploitation, movement, and stock size of saginaw bay walleyes, 1981–2011; 31 years of tag return analysis. Alpena Fisheries Research Station, Michigan Department of Natural Resources, Alpena, Michigan, p. 58; 2016.
Furlong P, Foster RF, Colby PJ, Friday M. Black Sturgeon River dam: a barrier to the rehabilitation of Black Bay walleye. Upper Great Lakes Management Unit, Ontario Ministry of Natural Resources and Forestry; Northern Bioscience Consulting, 27pp; 2006. http://s3.documentcloud.org/documents/552802/blacksturgeonreport.pdf.
Gavrilchuk K, Lesage V, Ramp C, Sears R, Bérubé M, Bearhop S, Beauplet G. Trophic niche partitioning among sypatric baleen whale species following the collapse of groundfish stocks in the Northwest Atlantic. Mar Ecol Prog Ser. 2014;497:285–301.
Harrison PM, Gutowsky LFG, Martins EG, Patterson DA, Cooke SJ, Power M. Personality-dependent spatial ecology occurs independently from dispersal in wild burbot (Lota lota). Behav Ecol. 2015;26(2):483–92. https://doi.org/10.1093/beheco/aru216.
Hayden TA, Holbrook CM, Fielder DG, Vandergoot CS, Bergstedt RA, Dettmers JM, Krueger CC, Cooke SJ. Acoustic telemetry reveals large-scale migration patterns of walleye in Lake Huron. PLoS ONE. 2014;9(12):1–20. https://doi.org/10.1371/journal.pone.0114833.
Henderson BA, Collins N, Morgan GE, Vaillancourt A. Sexual size dimorphism of walleye (Stizostedion vitreum vitreum). Can J Fish Aquat Sci. 2003;60(11):1345–52. https://doi.org/10.1139/f03-115.
Henderson BA, Wong JL, Nepszy SJ. Reproduction of walleye in Lake Erie: allocation of energy. Can J Fish Aquat Sci. 1996;53(1):127–33. https://doi.org/10.1139/f95-162.
Hestbeck JB, Nichols JD, Malecki RA. Estimates of movement and site fidelity using mark-resight data of wintering Canada geese. Ecology. 1991;72(2):523–33.
Hussey NE, Kessel ST, Aarestrup K, Cooke SJ, Cowley PD, Fisk AT, Harcourt RG, Holland KN, Iverson SJ, Kocik JF, Mills Flemming JE, Whoriskey FG. Aquatic animal telemetry: a panoramic window into the underwater world. Science. 2015;348(6240):1255642. https://doi.org/10.1126/science.1255642.
Huveneers C, Simpfendorfer CA, Kim S, Semmens JM, Hobday AJ, Pederson H, Stieglitz T, Vallee R, Webber D, Heupel MR, Peddemors V, Harcourt RG, Reynolds J. The influence of environmental parameters on the performance and detection range of acoustic receivers. Methods Ecol Evol. 2016;7(7):825–35. https://doi.org/10.1111/2041-210X.12520.
Kaufman SD, Morgan GE, Gunn JM. The role of ciscoes as prey in the trophy growth potential of walleyes. N Am J Fish Manag. 2009;29(2):468–77. https://doi.org/10.1577/m07-117.1.
Kessel ST, Cooke SJ, Heupel MR, Hussey NE, Simpfendorfer CA, Vagle S, Fisk AT. A review of detection range testing in aquatic passive acoustic telemetry studies. Rev Fish Biol Fish. 2013;24(1):199–218. https://doi.org/10.1007/s11160-013-9328-4.
Krueger CC, Holbrook CM, Binder TR, Vandergroot C, Hayden TA, Hondorp DW, Nate NA, Paige K, Riley S, Fisk AT, Cooke SJ. Acoustic telemetry observation systems: challenges encountered and overcome in the Laurentian Great Lakes. Can J Fish Aquat Sci. 2018;75(10):1755–63.
Lester NP, Dextrase AJ, Kushneriuk RS, Rawson MR, Ryan PA. Light and temperature: key factors affecting walleye abundance and production. Trans Am Fish Soc. 2004;133(3):588–605. https://doi.org/10.1577/T02-111.1.
Lucas MC, Baras E. Methods for studying spatial behaviour of freshwater fishes in the natural environment. Fish Fish. 2000;1(4):283–316. https://doi.org/10.1046/j.1467-2979.2000.00028.x.
Luttbeg B, Sih A. Risk, resources and state-dependent adaptive behavioural syndromes. Philos Trans R Soc B Biol Sci. 2010;365(1560):3977–90. https://doi.org/10.1098/rstb.2010.0207.
Matley JK, Faust MD, Raby GD, Zhao Y, Robinson J, MacDougall T, Hayden TA, Fisk AT, Vandergoot CS, Krueger CC. Seasonal habitat-use differences among Lake Erie’s walleye stocks. J Great Lakes Res. 2020;46(3):609–21. https://doi.org/10.1016/j.jglr.2020.03.014.
Montiglio PO, Dammhahn M, Dubuc Messier G, Réale D. The pace-of-life syndrome revisited: the role of ecological conditions and natural history on the slow-fast continuum. Behav Ecol Sociobiol. 2018. https://doi.org/10.1007/s00265-018-2526-2.
Nakagawa S, Schielzeth H. Repeatability for Gaussian and non-Gaussian data: a practical guide for biologists. Biol Rev. 2010;85:935–56.
Nakayama S, Rapp T, Arlinghaus R. Fast-slow life history is correlated with individual differences in movements and prey selection in an aquatic predator in the wild. J Anim Ecol. 2017. https://doi.org/10.1111/1365-2656.12603.
Poe GL, Lauber TB, Connelly NA, Creamer S, Ready RC, Stedman RC. Net benefits of recreational fishing in the Great Lakes Basin: a review of the literature (Issue HDRU Series No 13–10); 2013. http://www2.dnr.cornell.edu/hdru/pubs/HDRUReport13-10.pdf.
Raby GD, Vandergoot CS, Hayden TA, Faust MD, Kraus RT, Dettmers JM, Cooke SJ, Zhao Y, Fisk AT, Krueger CC. Does behavioural thermoregulation underlie seasonal movements in Lake Erie walleye? Can J Fish Aquat Sci. 2018;75:488–96. https://doi.org/10.1139/cjfas-2017-0145.
Réale D, Dingemanse NJ, Kazem AJN, Wright J. Evolutionary and ecological approaches to the study of personality. Philos Trans R Soc B Biol Sci. 2010;365(1560):3937–46. https://doi.org/10.1098/rstb.2010.0222.
Réale D, Garant D, Humphries MM, Bergeron P, Careau V, Montiglio P-O. Personality and the emergence of the pace-of-life syndrome concept at the population level. Philos Trans R Soc London Ser B Biol Sci. 2010;365(1560):4051–63. https://doi.org/10.1098/rstb.2010.0208.
Rennie MD, Ebener MP, Wagner T. Can migration mitigate the effects of ecosystem change? Patterns of dispersal, energy acquisition and allocation in Great Lakes lake whitefish (Coregonus clupeaformis). Adv Limnol. 2012;63:455–76. https://doi.org/10.1127/advlim/63/2012/455.
Rennie MD, Purchase CF, Lester N, Collins NC, Shuter BJ, Abrams PA. Lazy males? Bioenergetic differences in energy acquisition and metabolism help to explain sexual size dimorphism in percids. J Anim Ecol. 2008;77(5):916–26. https://doi.org/10.1111/j.1365-2656.2008.01412.x.
Robillard MM, Casselman JM, McLaughlin RL, Mackereth RW. Alternative growth histories in populations of Lake Supereior brook trout: critical support for partial migration. Biol Cons. 2011;144(7):1931–9. https://doi.org/10.1016/j.biocon.2011.03.022.
Robson DS, Regier HA. Sample size in Petersen mark-recapture experiments. Trans Am Fish Soc. 1964;93(3):215–26. https://doi.org/10.1577/1548-8659(1964)93.
Roff DA. The evolution of migration and some life history parameters in marine fishes. Environ Biol Fishes. 1988;22(2):133–46. https://doi.org/10.1007/BF00001543.
Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5. https://doi.org/10.1038/nmeth.2089.
Sharpe D. Your chi-square test is statistically significant: now what? Pract Assess Res Eval. 2015;20(8):1–10.
Shaw AK. Drivers of animal migration and implications in changing environments. Evol Ecol. 2016;30(6):991–1007. https://doi.org/10.1007/s10682-016-9860-5.
Venturelli P, Murphy C, Shuter B, Johnston T, van Coeverden de Groot P, Boag P, Casselman J, Montgomerie R, Wiegand M, Leggett W. Maternal influences on population dynamics: evidence from an exploited freshwater fish. Ecology. 2010;91(7):2003–12. https://doi.org/10.1890/09-1218.
Wilson MW, Ridlon AD, Gaynor KM, Gaines SD, Stier AC, Halpern BS. Ecological impacts of human induced animal behaviour change. Ecol Lett. 2020;23(10):1522–36. https://doi.org/10.1111/ele.13571.
White GC, Burnham KP. Program MARK: survival estimation from populations of marked animals. Bird Study. 1999;46(sup001):S120–39. https://doi.org/10.1080/00063659909477239.
We thank the staff at the Upper Great Lakes Management Unit of the Ontario Ministry of Natural Resources and Forestry (OMNRF) in Thunder Bay. Particularly we thank Captain Dave Montgomery for assistance with field work, Paul Drombolis for assistance with field work and training on fish ageing, Tony Chiodo for assistance with field work and data management, and Marilee Chase for assistance with fish surgeries. Thanks to Ken Sandilands (IISD-ELA) for providing light sensor data. We also thank Brian Shuter for his insight during this project, and input and encouragement from two anonymous reviewers and the Associate Editor to pursue further analysis which improved the manuscript.
Project funding was provided by the Canada Research Chairs program and NSERC Discovery awarded to Michael Rennie; the OMNRF; and the Canadian-Ontario Agreement for Water Quality and Ecosystem Health.
Ethics approval and consent to participate
Procedures for tagging Walleye followed Canadian Council on Animal Care approved procedures (Lakehead University file # 1465777; Supplemental Information).
Consent for publication
No personal information is included as data in this manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Supplemental text and figures referenced in main text.
. R-script for determining the number of fish at each receiver gate, the number of fish detected each year, the annual migratory status of fish, and the secchi measurements of the Black Bay study system.
. R-script for back calculating the lengths of acoustically tagged Walleye.
. R-script for t-tests to assess differences in age and length of migratory and resident Walleye.
. R-script for assessing the detection efficiency of receiver gates.
. R-script for assessing the distribution of male and female Walleye captured throughout the Black Bay study system in FWIN nets.
. R-script for converting GLATOS detections into the encounter history format for Program MARK.
. R-script for assessing repeatability in Walleye migratory status using GLMMs.
. R-script for determining the TOHA available to Black Bay Walleye through the open water season.
. Migratory status (1=migratory, 0=resident) of Black Bay Walleye in 2016 (M16), 2017 (M17), and 2018 (M18) with fish length to assess the influence of length and individual on migratory status.
. Maximum outbound extent of 2016 tagged Black Bay Walleye in 2016 (G16) and 2017 (G17). 1=BEI, 2=GEP, 3=EDI, 4=BLP and SIP.
. Secchi depth measurements of the four states of the Black Bay study system through the open water season of 2017.
About this article
Cite this article
McKee, G., Hornsby, R.L., Fischer, F. et al. Alternative migratory strategies related to life history differences in the Walleye (Sander vitreus). Mov Ecol 10, 10 (2022). https://doi.org/10.1186/s40462-022-00308-7