Overwintering aggregation patterns of European catfish Silurus glanis
Movement Ecology volume 11, Article number: 9 (2023)
Animal aggregation, particularly in large-bodied species, is both a fascinating and intriguing phenomenon. Here we analyzed the overwintering behavior of the European catfish, Silurus glanis Linnaeus, 1758, the largest freshwater fish in Europe. By tracking 47 subadults and adults in a shallow lake in southeastern France, we reported a consistent aggregative behavior across four successive winters. By implementing time series analysis and Cox proportional hazard models, we investigated the dynamics of these aggregations (formation, stability, dislocation), and the factors that govern it, whether external (temperature, time of the day) or specific to the fish (size, key individuals). These aggregations lasted 1.5–2 months and mainly took place in a single small 4 m-deep area whose environmental conditions (temperature, oxygen, substrate) did not differ from other parts of the lake. In some periods during winter, all tagged fish were aggregated, which suggests that a large proportion of the lake population gathered there. Low temperatures (below 9 °C) triggered the formation of aggregations. They became more stable with decreasing temperatures, while individuals more frequently left the aggregation, preferentially at dusk and at night, when temperatures increased. The largest individuals swam more frequently back and forth to the aggregation. Irrespective of their size, some individuals consistently arrived earlier in the aggregation in winter and left later. This predictable seasonal grouping of individuals and, more generally, the knowledge provided by such studies on how species use space have important operational value and are useful for species conservation as well as for species control.
Aggregations of individuals are widespread in the animal kingdom, spanning a large range of sizes and time durations . Though they have been observed for a long time, they still remain striking and intriguing for biologists . Aggregations are particularly common among fishes, both marine and freshwater, most of which form cohesive social groups at some stage of their life . They have been shown to bring numerous benefits, including protection from predators [26, 35], increased probability of encountering a mate [12, 16, 21], increased foraging efficiency [17, 25], reduction of energy expenditure [24, 38] and centralized information . These benefits are assumed to outweigh the costs incurred by aggregations, such as within-group competition for resources  and exposure to parasites [63, 67]. On an evolutionary time scale, aggregation behavior is thus expected to increase individual fitness . Aggregations in fish range from relatively small shoals, usually in freshwater streams and lakes (e.g. [15, 22, 31, 68], to structured schools of up to hundreds of thousands of individuals in marine systems (e.g. . In most fish species, aggregations are stable over time , but for others they can be transitory, for example when fish are attracted by aggregating devices (e.g. ) or when they spawn (e.g. [12, 16, 46]).
In addition to spawning, in some species, temporary aggregations have been described for overwintering, a phase in which fish individuals are exposed to numerous stressors, e.g. starvation, thermal stress, but this has rarely been studied . In freshwater systems, these wintering aggregations have long been described in common carp Cyprinus carpio [3, 29, 46, 48, 65] and more recently in another large and long-lived species, the lake sturgeon Acipenser fulvescens, the largest freshwater fish in North America . Winter aggregations also occur in centrarchid fishes during the light phase of the day . However, it is still an understudied aspect of their ecology. As regards to European catfish Silurus glanis, the largest European freshwater fish, only a single study has reported the existence of repeated aggregations involving 15–44 adults . These aggregations occurred throughout the year at the same place in a large river . Brevé et al.  also observed an aggregation of catfish adults in a river section underneath boats, but just once. Yet, in a small shallow eutrophic lake, Vanovac et al.  found that while common carp were aggregating during winter, a stocked population of European catfish was not. European catfish are reported to be solitary foragers , that preferentially feed at night. They expend more energy when in contact with conspecifics in preferred areas of habitat, probably to monopolize resources , at least in their native range. Aggregations in European catfish warrant further work .
While the ultimate causes of aggregations have been extensively studied and, in some cases, are well established, the proximate mechanisms that underlie the formation and dislocation of aggregations are less understood . The transitory nature of winter aggregations provides an opportunity to study these mechanisms. Moreover, most studies on animal aggregations have dealt with collective coordinated behaviors that emerge from interactions between individuals that are considered equivalent . Yet, there is a growing body of evidence that individual variability can play an important role in these aggregations  and recent advances in high resolution tracking of individuals now provide tools to investigate its role in aggregation dynamics .
In this paper, we present the results of an experiment in which 47 subadults and adults of the European catfish were tracked by acoustic telemetry for 4 years in a shallow eutrophic lake located in southeastern France. We focus on their overwintering behavior by analyzing their movements over four successive winters showing contrasted temperatures; we found that every winter, catfish exhibited an aggregative behavior. We investigated the dynamics of these aggregations (formation, stability, dislocation) as well as the factors that could govern it, whether related to external conditions (temperature, time of the day) or to the characteristics of the fish (size, key individuals).
Material and methods
“Etang des Aulnes” is a 104-ha shallow natural lake of 3.8 m mean depth and 6 m maximum depth, located in South-Eastern France in a protected natural area (Fig. 1). The fish assemblage in the lake, determined by fyke nets, fishing traps and electro fishing in October 2017, 2018 and 2019 was composed of 16 species. The most dominant species were common bream (Abramis brama, relative abundance 65%), European perch (Perca fluviatilis, 13%), pumpkinseed (Lepomis gibbosus, 8%), tench (Tinca tinca, 4%), pikeperch (Sander lucioperca, 4%), European catfish (Silurus glanis, 3%) and Northern pike (Esox lucius, 2%). In addition, two crayfishes were present: Procambarus clarkii and Faxonius limosus. Fishing is allowed but only during daytime from the eastern bank of the lake and no other activity is authorized.
Physical and chemical lake conditions
The “Etang des Aulnes” is eutrophic (see  for chemical element concentrations). Hourly temperature profiles were recorded with HOBO data loggers U22 (0.2 °C accuracy) at the deepest point in the lake (location 42 on Fig. 1) and at several other locations, among which points 5 and 38; they are presented later. During the winter period, temperatures are mainly homogenized in the water column; the characteristics of each winter are given in Table 1. At the deepest point in the lake, hourly profiles of dissolved oxygen concentration were also recorded with HOBO data loggers U26 (0.2 mg/L accuracy).
Specifically, a total of 47 catfish, subadults and adults, were caught by fyke nets, angling or electrofishing over four sampling campaigns: 10, 32, 2 and 3 individuals, respectively in October 2017, 2018, 2019 and 2020. Once caught, fish was anesthetized by immersion for 5–6 min in a benzocaine solution (80 mg/L). A 15–20 mm-long incision was made on the ventral side to insert the tag and was closed using 2–3 simple surgical sutures (3–0 Polydioxanon resorbable monofilament). The surgery itself took 5–6 min while the fish was in another solution of benzocaine (40 mg/L). Fish recovered in 10 min and were released after 3–6 h. The surgical procedure followed for fish tagging is detailed in Westrelin et al. . At the end of 2018 and 2019 winters, respectively 2 and 3 tags were stationary and, consequently, the corresponding individuals were discarded from the analyses as they were assumed to be dead or to have expelled their tag. Vemco V13-1L acoustic transmitters (length: 30.5 mm, weight: 9.2 g in the air, mean battery life: 1825 days, mean burst interval: 180 s—range 120–240 s—for the 12 used in 2017 and 320 s—range 260–80 s—in 2018) were used. The transmitter weight in the air did not exceed 2% of the fish body weight [55, 72]. The characteristics of the tracked fish are given in Table 1.
Fifty-two underwater omnidirectional Vemco acoustic receivers (20 VR2W 69 kHz and 32 VR2Tx 69 kHz) with their associated synchronization tag (additional V16-1L transmitter for VR2W and built-in V16-like transmitter for VR2Tx, 500–700 s), used to correct for receiver internal clock drift, were anchored to the bottom of the lake in October 2017 (Fig. 1). Seven reference tags (V13-1L, 840–960 s) were added to detect anomalies in the tracking system. On average, neighboring receivers were positioned 155 m from each other (range, 100–209 m), at a depth of 3.9 m (range, 1.5–6 m), 0.5 m above the bottom of the lake. Receivers were removed roughly every 6 months to download fish detections. From these detections, fish 2D positions were calculated with the Vemco Positioning System (VPS) . The horizontal position error, a dimensionless parameter calculated by the VPS for each position, gives information on the quality of the position estimate, and was used to filter the data set . Here, we retained only positions with horizontal position error below 100; this limit represented a good compromise between the mean position error (7.4 m, calculated on reference tags) and the percentage of positions kept (87%).
Space use metrics
Fish were continuously tracked from October 2017 onward. As our study focuses on the winter period, only data recorded from the 15th of October to the 15th of March of the next year were analyzed. To obtain synchronized individual tracks, individual raw positions were interpolated using the R package trajr  for each quarter hour between the first and the last position recorded. In order to visualize catfish space use in winter, all tracks were then plotted on the lake map to create videos of the catfish displacements (Additional file 1). The videos clearly show an aggregation localized in the western part of the lake.
To identify the aggregation zone over the 5-month periods, the home ranges of the pooled fish were estimated with an Epanechnikov kernel as the utilization distribution with probability levels 95%, 50% , and the level corresponding to the highest percentage delineating the aggregation zone identified on the videos. The home ranges were estimated for each phase of the daily cycle, i.e. dawn, day, dusk and night, defined at an hourly resolution. Dawn was defined as the period including the hour preceding the sunrise hour, the sunrise hour itself and the following hour. Dusk was defined as the period including the hour preceding the sunset hour, the sunset hour itself and the following hour. Dawn and dusk lasted 3 h each. Daytime was the period following dawn and preceding dusk and night was the period following dusk and preceding dawn. These spatial analyses were conducted using the R package adehabitatHR . To quantify the degree of aggregation over time, the mean distance between individuals was calculated each quarter hour using the R package spatstat  and the number of fish in the aggregation zone was counted.
To isolate the aggregation period, we applied an algorithm to detect possible breakpoints corresponding to structural changes in the 5-month time series of the mean distance between individuals (strucchange R package, [75, 76]. We considered that for each of the four winters, the periods with the lowest mean distances between individuals corresponded to the time at which aggregation occurred; these periods were confirmed by watching the videos. To investigate the formation and dislocation of the aggregation, these periods were extended to before, when no individual had joined the aggregation zone yet, and after, until all individuals had left it. These periods will hereafter be referred to as the extended aggregation periods.
To investigate whether certain individuals consistently joined the aggregation earlier than others at the beginning of the aggregation period, we performed a comparison of the ranks at which individuals first joined the aggregation over successive winters using a Friedman test. To highlight whether fish size influenced the timing of aggregation, we then compared the mean rank of joining the aggregation of the different size classes with a Kruskall-Wallis test. Fish size was defined from the total body length measured during tagging and was categorized into three classes: “Small”, “Large” and “Medium”, corresponding to total length < 850 mm, ≥ 1100 mm and in-between, respectively (Table 1). Fish were only measured during tagging as very few were recaptured; we hypothesized they did not switch to other size-classes during the study. The same analysis was performed at the end of the aggregation period, when the aggregation disbanded, by considering the rank at which individuals permanently left the aggregation. To quantify a possible link between the mean rank of arrival in and the mean rank of departure from the aggregation, a Spearman correlation coefficient was calculated. These analyses were only performed over the last three winters, when a common significant pool of individuals was present (38 individuals throughout these 3 winters).
To investigate the effect of temperature, time of day and fish size on the stability of winter aggregations, we used two multivariate mixed effects Cox proportional hazard models (coxme R package ). The first model was used to assess the effect of the covariates on the rate of temporarily leaving the aggregation for an excursion. While the second model was used to assess the effect of the covariates on the rate of returning to the aggregation after an excursion. The two models can be formulated as follows:
Survival(Start, Stop, Event) ~ TEMPERATURE + FISH SIZE + TIME OF DAY + (1|Fish identity)
where TEMPERATURE is the temperature at 3 m above the bottom of the lake at its deepest point, FISH SIZE is the fish class size (3 levels: Small, Medium, Large) and TIME OF DAY is the day period (4 levels: Dawn, Day, Dusk, Night).
In the first model, Event corresponds to the “leave the aggregation” behavior; for each Event, the time-to-event ranges from the time at which a given individual joined the aggregation (Start) to when it left it (Stop). In the second model, the Event corresponds to “join the aggregation”; for each Event, the time-to-event goes from the time at which a given individual left the aggregation (Start) to the time at which it joined it (Stop).
Fish identity was considered as a random effect to account for individual variability and repeated measurements made on the same individuals. The Cox models were run on the whole dataset which included the four pooled winters. Since two of the fixed effects (temperature and time of day) varied with time, the dataset was rearranged so that each quarter hour observation of a given individual fish was treated as a separate observation, i.e. containing a Start and Stop time and the corresponding Event type. The effects of significant categorical covariates were further analyzed by Tukey comparisons of pairwise estimated marginal means of their different levels (emmeans R package ).
The survival functions, which represent the probabilities of time-to-event over time, were estimated using the Kaplan–Meier method  and plotted for the different covariates (survival R package, . For the temperature effect, the following classes were used to plot the survival curves: [3; 5[, [5; 7[, [7; 9[ and [9; 14] °C; the range of the warmest class is larger to avoid a small sample size as temperatures above 11 °C rarely occurred in winter.
All statistical analyses and graphics were made using R version 3.6.3 .
The aggregation period lasted 41.9, 67.4, 55.2 and 59.7 days and the extended aggregation 49.8, 93.6, 113.0 and 76.2 days, respectively in 2017, 2018, 2019 and 2020 (Table 1) (Additional file 2 details the breakpoints analysis). The difference between these two periods corresponds to the time for the formation and dislocation of the aggregation, i.e., 7.8, 26.2, 57.8 and 16.5 days, respectively in 2017, 2018, 2019 and 2020.
Catfish exhibited strong fidelity to the same aggregation zone across the four winters, although they were more scattered in the milder 2019 winter (Fig. 2, Table 1). The aggregation area ranged between 2.1 and 3.2 ha, if we exclude 2019 when the aggregation was not as dense as in other winters (Additional file 1) and was split into two main zones yielding a total area of 4.4 ha. However, one of the two zones was the same as in the other winters and extended over 2.0 ha (Fig. 2). The aggregation zones corresponded to a probability of utilization of 35%, 50%, 30% and 40% over the 15 October–15 March period, respectively in 2017, 2018, 2019 and 2020, which means that fish spent 30–50% of their time in this zone over this 5-month period. If this spatial analysis is restricted to the extended aggregation period, the probability of utilization rises to 60%, 70%, 40% and 70%, respectively in 2017, 2018, 2019 and 2020. In the 2017, 2018 and 2020 winters, during the identified aggregation period, most individuals were indeed located within the aggregation zone (on average 68.0, 72.4 and 68.9% over the extended aggregation period in 2017, 2018 and 2020 respectively; Fig. 3a, b, d), but aggregation was weaker in 2019 (in average 38.7%, Fig. 3c). The identified aggregation periods corresponded to the time at which the coldest temperatures were recorded (Fig. 3 and Table 1). The aggregation zones were very stable across the different times of day (Fig. 4).
The main aggregation zone, which could be identified throughout the four winters, was on average 4.3 m deep (range [2.8; 4.8] m) and at a distance of 107.2 m from the bank (range [42.8; 196.6] m). The secondary aggregation zone, which appeared only in winter 2019, was on average 4.2 m deep (range [3.7; 4.6] m) and at a distance of 134.1 m from the bank (range [62.5; 215] m). The mean daily temperature differences between the main aggregation zone and the deepest point in the lake, which stands as a reference point, ranged between [− 0.2; 0.4] °C throughout the extended aggregation period and had a 0.0 °C mean over this period (Fig. 2, Additional file 3). The mean daily differences of temperature between a location close to the secondary aggregation zone and the deepest point in the lake ranged between [− 0.4; 0.5] °C throughout the extended aggregation period and had a − 0.1 °C mean over this period (Fig. 2, Additional file 4).
Over the extended aggregation period, at the deepest point in the lake, mean daily oxygen concentrations ranged between [9.8; 14.2], [6.1; 15.1], [9.4; 12.9] and [8.2; 13.7] mg/L, corresponding to saturation rates that ranged between [83.5; 113.9], [48.3; 134.8], [81.3; 112.8] and [70.2; 105.4] %, respectively in 2017, 2018, 2019 and 2020.
There was no significant consistency in the rank of first arrival of individuals in the aggregation between the three last studied winters (Friedman test: Chi2 = 51.56, df = 37, p = 0.056). However, the rank of first arrivals was comparatively consistent for the individual with the lowest mean rank (mean rank 4; range [1; 8] for individual 1030_18b) and for the two individuals with the highest mean ranks (mean rank 34.3 and 34.7; range [31; 37] and [34; 36], respectively for individuals 920_18 and 866_17) (Fig. 5a). Moreover, the rank of first arrival did not depend on fish size (Kruskall-Wallis test: Chi2 = 2.13, df = 2, p = 0.345), although the mean rank of arrival in the aggregation of 5 out of the 8 small fish was higher than that of 70% of all fish considered in the analysis (n = 38). Neither was there any consistency in the rank of last departure of the fish from the aggregation (Friedman test: Chi2 = 45.10, df = 37, p = 0.169), nor any size effect (Kruskall-Wallis test: Chi2 = 0.97, df = 2, p = 0.614). Nevertheless, there was much less variability in the rank of last departure of the three individuals with the lowest mean ranks (mean rank 3.33, 3.67 and 4.47; range [2; 5], [1; 9] and [3; 7], respectively for individuals 920_18, 866_17 and 839_17) and of the two individuals with the highest mean ranks (mean rank 30.7 and 33.7; range [28; 34] and [29; 36], respectively for individuals 873_18 and 1030_18b) (Fig. 5b). The correlation between the mean rank of arrival and mean rank of departure was significantly negative (Spearman rho = − 0.43, p = 0.007) meaning that, on average, the first individuals to arrive in the aggregation were the last to leave and vice versa.
There was a significant effect of temperature, fish size and time of day on the probability of leaving the aggregation for an excursion (Table 2a, b) and on the probability of coming back to the aggregation after an excursion (Table 3a, b). The probability of leaving the aggregation was significantly higher with increasing water temperature (Fig. 6a, Table 2c), meaning that, on average, fish stayed longer in the aggregation when the water temperature remained low. This probability was similar between small and medium fish but increased for large fish (Fig. 6b, Table 2c, d). Small and medium fish therefore made longer stays in the aggregation than large fish. Moreover, the probability of leaving the aggregation did not vary between dawn and day (Fig. 6c, Table 2d) but increased between day and dusk (Table 2c). It decreased between dusk and night and between night and dawn (Fig. 6c, Table 2c), which means that fish left the aggregation mostly at dusk and at night, but even more often at dusk.
The probability of returning to the aggregation after an excursion significantly decreased with increasing temperature (Table 3c), especially when water temperature exceeded 9 °C (Fig. 6d). Therefore fish made longer excursions out of the aggregation when the temperature was high. The probability of coming back to the aggregation did not differ between fish of small and medium size (Table 3d) but significantly decreased for large fish (Table 3c, Fig. 6e). Large fish thus made longer excursions out of the aggregation than small or medium fish. Finally, the probability for fish of coming back into the aggregation decreased between dawn and day and between dusk and night (Table 3c, Fig. 6f). It increased between day and dusk and between night and dawn (Table 3c, Fig. 6f). Ranked in order of importance, fish thus went back into the aggregation first at dawn, then at dusk and last at night.
Catfish aggregated across four winters for 1.5–2 months and showed a strong fidelity to the same zone representing 2–4% of the lake area. Even in other seasons, catfish have been shown to exhibit strong site fidelity [6, 10, 11, 52], but such an aggregation behavior has only been described in detail throughout a single year, and only once in a large river . Site fidelity in fish can generally be linked to environmental features that provide a benefit to individuals, such as warmer temperatures [6, 10] or a refuge area for catfish , deep and slow pools for lake sturgeon Acipenser fulvescens , proximity to beds of emergent vegetation or open water formed by turbulence from a lake aerator for common carp Cyprinus carpio  or current updrafts that reduce energy expenditure for grey reef sharks Carcharhinus amblyrhynchos . In our study however, the aggregation zone was not warmer than anywhere else in the lake and oxygen conditions were not limiting; moreover, a scan of the aggregation area with an acoustic camera (sonar 2D Oculus, 1.2 MHz ) on the 12th of February 2020 revealed no particular structure at the bottom of the lake. Finally, the main aggregation zone did not appear to be sheltered from the prevailing winds that blow from the north in this region and thus did not offer calmer water compared to the other parts of the lake.
The formation and stability of the catfish aggregation in our study were closely linked to temperature. Similar observations have been made in common carp that aggregate when water temperature drops below 8 °C . Fewer movements associated with the aggregation behavior is a way to save energy . Aggregations were fairly stable across different times of day. Most movements took place at dusk and at night when some individuals left the aggregation and at dawn and at dusk when they came back. This is in agreement with the preferential nocturnal activity reported for catfish in the literature [6, 11]. The largest individuals spent more time outside of the aggregation. As the whole-organism metabolic rate of an individual increases with its body size (allometric equation, see ), large fish may need to leave the aggregation for feeding more often.
Even if not significant, the rank of arrival of individuals in the aggregation was more or less consistent across winters (p = 0.056), which could suggest that some individuals had a leading position (e.g. ). However, this was not the case for the rank of departure. On one hand, there was a tendency for some individuals to spend as much time as possible within the aggregation by arriving earlier and leaving later, and, on the other hand, there was a tendency for some others to minimize this time. Thus, leaders in the aggregation formation are not likely to be leaders in its dislocation, quite the contrary. The possible leading position was not linked to fish size contrary to what was observed for leadership in roach shoals  or for dominance in the catfish itself .
Aggregation causes and consequences
Up to 100% of tagged fish gathered in the winter aggregation, which is much greater than the 23% reported for lake sturgeon  or the 70% reported for common carp . Based on capture-recapture data (unpublished data), the estimated catfish population (individuals greater than 600 mm) in “Etang des Aulnes” is around 770 individuals (95% confidence interval [184; 1356]). As tagged individuals can be considered as a representative sample of the subadult and adult catfish population and as all tagged individuals at one point or another took part in the aggregation, a great proportion of the population could potentially gather there. This would thus make a huge winter aggregation of several hundreds individuals in comparison to the 44 individuals observed by Boulêtreau et al. . In aggregations such as those that we observed, regrouping several hundreds individuals over 1.5–2 months, the level of competition between individuals could be considerably high . This could be particularly true for catfish that have been reported to actively defend their access to resources in their core area [14, 52] and to be solitary foragers  during active periods. However, this competition could be here considerably dampened due to the cold season when catfish usually feed very little .
The cause or function of the aggregative behavior we observed remains unknown. The temperature was far below suitable temperatures for spawning (between 20 and 25 °C ), and the aggregation period occurs far too early to allow fish to identify potential mates for the spawning season that occurs several months later (usually in May–June in this lake). Unlike size-assortative schooling (e.g. ), individuals of various sizes could be found in the aggregation (range [727; 2150] mm in our tagged individuals), among which the smallest had probably not reached a refuge size against the largest individuals yet, since the prey-to-predator length ratio for catfish can reach 0.57 . However, we do not know whether smaller individuals also aggregated. Some studies suggest that social interactions between conspecifics could play a role in fish aggregation, challenging the classical view of aggregation formation around floating devices . Moreover, the site fidelity for aggregating could favor social interactions  like non-random aggregations in sharks . As discussed above, regardless of their size, some individuals tended to stay longer in the aggregation while others seemed to shorten their stay. Irrespective of a possible hierarchy between individuals, longer stays in the aggregations could favor social interactions with conspecifics. These social interactions could regulate the stress between catfish, as was already shown in some species [1, 37]. However, Slavík and Horký  reported that catfish increase their energy consumption when in contact with conspecifics, assuming this was a stressful situation. This may not apply to our study though, since these authors made their observations on males only and in spring, when catfish return to normal activity and prefer to hunt solitarily . Catfish have also been reported to decrease their activity in the presence of familiar conspecifics . This raises the question of the role of animal personality and individual heterogeneity in collective behavior . For example, shy individuals could be more likely to cooperate whereas bold ones could act more independently [13, 34]. Based on the calculation of a proximity index, Vanovac et al.  concluded that catfish do not display within-species interactions whereas common carp do, especially in winter. We may however question the sensitivity of this index to the sample size of fish. In fact, in Fig. 4 of their paper, which represents the location of species across seasons with kernel densities, one can clearly see a tendency for catfish to cluster in winter, similar to our observations.
Due to adults being at least twice as large as native fish predators, catfish are usually considered as a ‘giant’ top predator  and are suspected to threaten the fish communities [18, 23, 64]. Therefore, in numerous ecosystems, managers aim to control their population. However, a fundamental constraint of control methods is their lack of selectivity in specifically removing non-native fish species . If our observations can be generalized to other locations and a considerable fraction of catfish populations also gather in the same zone during the coldest periods of the year, tagging only a few individuals could help locate and remove most individuals of a population (population control by the Judas technique) .
Long-lasting winter aggregations of catfish, composed of individuals of various sizes and that likely concern a large proportion of the population, were shown to occur consistently every year at the same place in a restricted area of the lake. The aggregation formation and stability was closely linked to low temperatures. The area where the fish gathered was moderately deep and not different from other parts of the lake. Some individuals seemed to spend a longer time in the aggregation consistently every winter. Further studies are needed to explain this result, but one hypothesis is that it could correspond to different behavioral types, with certain individuals seeking sociality while others are more independent. In the end, this predictable seasonal grouping of individuals could provide an opportunity for lake managers to efficiently control catfish population if needed. As for species conservation, if the overwintering habitat is a critical requirement, it could constitute a bottleneck habitat, crucial to secure and maintain. The knowledge provided by our study has thus both academic and operational implications.
Availability of data and materials
The datasets analyzed during the current study are available from the corresponding author on reasonable request.
Allen PJ, Barth CC, Peake SJ, Abrahams MV, Anderson WG. Cohesive social behaviour shortens the stress response: the effects of conspecifics on the stress response in lake sturgeon Acipenser fulvescens. J Fish Biol. 2009;74:90–104. https://doi.org/10.1111/j.1095-8649.2008.02112.x.
Baddeley A, Rubak E, Turner R. Spatial point patterns: methodology and applications with R. London: Chapman and Hall/CRC Press; 2015.
Bajer PG, Chizinski CJ, Sorensen PW. Using the Judas technique to locate and remove wintertime aggregations of invasive common carp. Fish Manag Ecol. 2011;18:497–505. https://doi.org/10.1111/j.1365-2400.2011.00805.x.
Boulêtreau S, Cucherousset J, Villéger S, Masson R, Santoul F. Colossal aggregations of giant alien freshwater fish as a potential biogeochemical hotspot. PLoS ONE. 2011. https://doi.org/10.1371/journal.pone.0025732.
Boulêtreau S, Santoul F. The end of the mythical giant catfish. Ecosphere. 2016. https://doi.org/10.1002/ecs2.1606.
Brevé NWP, Verspui R, de Laak GAJ, Bendall B, Breukelaar AW, Spierts ILY. Explicit site fidelity of European catfish (Silurus glanis, L., 1758) to man-made habitat in the River Meuse, Netherlands. J Appl Ichthyol. 2014;30:472–8. https://doi.org/10.1111/jai.12410.
Britton JR, Gozlan RE, Copp GH. Managing non-native fish in the environment. Fish Fish. 2011;12:256–74. https://doi.org/10.1111/j.1467-2979.2010.00390.x.
Brown JH, Gillooly JF, Allen AP, Savage VM, West GB. Toward a metabolic theory of ecology. Ecology. 2004;85:1771–89.
Calenge C. The package “adehabitat” for the R software: a tool for the analysis of space and habitat use by animals. Ecol Model. 2006;197:516–9. https://doi.org/10.1016/j.ecolmodel.2006.03.017.
Capra H, Pella H, Ovidio M. Individual movements, home ranges and habitat use by native rheophilic cyprinids and non-native catfish in a large regulated river. Fish Manag Ecol. 2018;25:136–49. https://doi.org/10.1111/fme.12272.
Carol J, Zamora L, García-Berthou E. Preliminary telemetry data on the movement patterns and habitat use of European catfish (Silurus glanis) in a reservoir of the River Ebro, Spain. Ecol Freshw Fish. 2007;16:450–6. https://doi.org/10.1111/j.1600-0633.2007.00225.x.
Clevenstine AJ, Lowe CG. Aggregation site fidelity and movement patterns of the protected marine predator giant sea bass (Stereolepis gigas). Environ Biol Fishes. 2021;104:401–17. https://doi.org/10.1007/s10641-021-01077-9.
Conrad JL, Weinersmith KL, Brodin T, Saltz JB, Sih A. Behavioural syndromes in fishes: a review with implications for ecology and fisheries management. J Fish Biol. 2011;78:395–435. https://doi.org/10.1111/j.1095-8649.2010.02874.x.
Cucherousset J, Horky P, Slavik O, Ovidio M, Arlinghaus R, Boulêtreau S, Britton R, Garcia-Berthou E, Santoul F. Ecology, behaviour and management of the European catfish. Rev Fish Biol Fish. 2018;28:177–90. https://doi.org/10.1007/s11160-017-9507-9.
Currie HAL, White PR, Leighton TG, Kemp PS. Collective behaviour of the European minnow (Phoxinus phoxinus) is influenced by signals of differing acoustic complexity. Behav Process. 2021;189:104416. https://doi.org/10.1016/j.beproc.2021.104416.
Daly R, Filmalter JD, Daly CAK, Bennett RH, Pereira MAM, Mann BQ, Dunlop SW, Cowley PD. Acoustic telemetry reveals multi-seasonal spatiotemporal dynamics of a giant trevally Caranx ignobilis aggregation. Mar Ecol Prog Ser. 2019;621:185–97. https://doi.org/10.3354/meps12975.
Day RL, MacDonald T, Brown C, Laland KN, Reader SM. Interactions between shoal size and conformity in guppy social foraging. Anim Behav. 2001;62:917–25. https://doi.org/10.1006/anbe.2001.1820.
De Santis V, Volta P. Spoiled for choice during cold season? habitat use and potential impacts of the invasive Silurus glanis L. in a deep, large, and oligotrophic lake (lake Maggiore, north Italy). Water. 2021;13:2549.
Domenici P, Herbert NA, Lefrançois C, Steffensen JF, McKenzie DJ. The effect of hypoxia on fish swimming performance and behaviour Swimming physiology of fish. Berlin: Springer; 2013. p. 129–59.
Espinoza M, Farrugia TJ, Webber DM, Smith F, Lowe CG. Testing a new acoustic telemetry technique to quantify long-term, fine-scale movements of aquatic animals. Fish Res. 2011;108:364–71. https://doi.org/10.1016/j.fishres.2011.01.011.
Fox RJ, Bellwood DR, Jennions MD. Why pair? Evidence of aggregative mating in a socially monogamous marine fish (Siganus doliatus, Siganidae). R Soc Open Sci. 2015;2:150252. https://doi.org/10.1098/rsos.150252.
Freeman MC, Grossman GD. Group foraging by a stream minnow: shoals or aggregations? Anim Behav. 1992;44:393–403. https://doi.org/10.1016/0003-3472(92)90050-J.
Guillerault N, Delmotte S, Boulêtreau S, Lauzeral C, Poulet N, Santoul F. Does the non-native European catfish Silurus glanis threaten French river fish populations? Freshw Biol. 2015. https://doi.org/10.1111/fwb.12545.
Hemelrijk C, Reid D, Hildenbrandt H, Padding J. The increased efficiency of fish swimming in a school. Fish Fish. 2015;16:511–21. https://doi.org/10.1111/faf.12072.
Herbert-Read JE, Romanczuk P, Krause S, Strömbom D, Couillaud P, Domenici P, Kurvers RHJM, Marras S, Steffensen JF, Wilson ADM, Krause J. Proto-cooperation: group hunting sailfish improve hunting success by alternating attacks on grouping prey. Proc R Soc B Biol Sci. 2016;283:20161671. https://doi.org/10.1098/rspb.2016.1671.
Herbert-Read JE, Rosén E, Szorkovszky A, Ioannou CC, Rogell B, Perna A, Ramnarine IW, Kotrschal A, Kolm N, Krause J, Sumpter DJT. How predation shapes the social interaction rules of shoaling fish. Proc R Soc B Biol Sci. 2017;284:20171126. https://doi.org/10.1098/rspb.2017.1126.
Heupel MR, Simpfendorfer CA. Quantitative analysis of aggregation behavior in juvenile blacktip sharks. Mar Biol. 2005;147:1239–49. https://doi.org/10.1007/s00227-005-0004-7.
Jacoby DMP, Papastamatiou YP, Freeman R. Inferring animal social networks and leadership: applications for passive monitoring arrays. J R Soc Interface. 2016. https://doi.org/10.1098/rsif.2016.0676.
Johnsen PB, Hasler AD. Winter aggregations of carp (Cyprinus carpio) as revealed by ultrasonic tracking. Trans Am Fish Soc. 1977;106:556–9. https://doi.org/10.1577/1548-8659(1977)106%3c556:WAOCCC%3e2.0.CO;2.
Jolles JW, King AJ, Killen SS. The role of individual heterogeneity in collective animal behaviour. Trends Ecol Evol. 2020;35:278–91. https://doi.org/10.1016/j.tree.2019.11.001.
Krause J, Hoare DJ, Croft D, Lawrence J, Ward A, Ruxton GD, Godin JGJ, James R. Fish shoal composition: mechanisms and constraints. Proc R Soc Lond B. 2000;267:2011–7. https://doi.org/10.1098/rspb.2000.1243.
Krause J, Reeves P, Hoare D. Positioning behaviour in roach shoals: the role of body length and nutritional state. Behaviour. 1998;135:1031–9. https://doi.org/10.1163/156853998792913519.
Krause J, Ruxton G. Living in groups. Oxford University Press; 2002.
Krause J, Ward AJ, James R, Croft DP. Group living and social networks. In: Magnhagen C, Braithwaite VA, Forsgren E, Kapoor BG, editors. Fish behaviour. Enfield: Science Publishers; 2008.
Larsson M. Possible functions of the octavolateralis system in fish schooling. Fish Fish. 2009;10:344–53. https://doi.org/10.1111/j.1467-2979.2009.00330.x.
Lenth RV. Least-squares means: the R package lsmeans. J Stat Softw. 2016;69:1–33. https://doi.org/10.18637/jss.v069.i01.
Magnhagen C. Personalities in a crowd: what shapes the behaviour of Eurasian perch and other shoaling fishes? Curr Zool. 2015;58:35–44. https://doi.org/10.1093/czoolo/58.1.35.
Marras S, Killen SS, Lindström J, McKenzie DJ, Steffensen JF, Domenici P. Fish swimming in schools save energy regardless of their spatial position. Behav Ecol Sociobiol. 2015;69:219–26. https://doi.org/10.1007/s00265-014-1834-4.
Martignac F, Daroux A, Bagliniere JL, Ombredane D, Guillard J. The use of acoustic cameras in shallow waters: new hydroacoustic tools for monitoring migratory fish population. A review of DIDSON technology. Fish Fish. 2015;16:486–510. https://doi.org/10.1111/faf.12071.
McLean DJ, Skowron Volponi MA. trajr: an R package for characterisation of animal trajectories. Ethology. 2018;124:440–8.
Moreno G, Dagorn L, Capello M, Lopez J, Filmalter J, Forget F, Sancristobal I, Holland K. Fish aggregating devices (FADs) as scientific platforms. Fish Res. 2016;178:122–9. https://doi.org/10.1016/j.fishres.2015.09.021.
Mourier J, Vercelloni J, Planes S. Evidence of social communities in a spatially structured network of a free-ranging shark species. Anim Behav. 2012;83:389–401. https://doi.org/10.1016/j.anbehav.2011.11.008.
Nathan R, Monk CT, Arlinghaus R, Adam T, Alós J, Assaf M, Baktoft H, Beardsworth CE, Bertram MG, Bijleveld AI, Brodin T, Brooks JL, Campos-Candela A, Cooke SJ, Gjelland KØ, Gupte PR, Harel R, Hellström G, Jeltsch F, Killen SS, Klefoth T, Langrock R, Lennox RJ, Lourie E, Madden JR, Orchan Y, Pauwels IS, Říha M, Roeleke M, Schlägel UE, Shohami D, Signer J, Toledo S, Vilk O, Westrelin S, Whiteside MA, Jarić I. Big-data approaches lead to an increased understanding of the ecology of animal movement. Science. 2022;375:eabg1780. https://doi.org/10.1126/science.abg1780.
Papastamatiou YP, Iosilevskii G, Di Santo V, Huveneers C, Hattab T, Planes S, Ballesta L, Mourier J. Sharks surf the slope: current updrafts reduce energy expenditure for aggregating marine predators. J Anim Ecol. 2021;90:2302–14. https://doi.org/10.1111/1365-2656.13536.
Parrish JK, Edelstein-Keshet L. Complexity, pattern, and evolutionary trade-offs in animal aggregation. Science. 1999;284:99–101. https://doi.org/10.1126/science.284.5411.99.
Penne CR, Pierce CL. Seasonal distribution, aggregation, and habitat selection of common carp in Clear Lake, Iowa. Trans Am Fish Soc. 2008;137:1050–62. https://doi.org/10.1577/T07-112.1.
Peuhkuri N, Ranta E, Seppä P. Size-assortative schooling in free-ranging sticklebacks. Ethology. 1997;103:318–24. https://doi.org/10.1111/j.1439-0310.1997.tb00021.x.
Pitcher TJ. Shoaling and shoaling behaviour in fishes. In Greenberg GMMH (ed) Comparative psychology: a handbook; 1998. p. 748–760.
R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2020.
Robert M, Dagorn L, Deneubourg JL. The aggregation of tuna around floating objects: what could be the underlying social mechanisms? J Theor Biol. 2014;359:161–70. https://doi.org/10.1016/j.jtbi.2014.06.010.
Schober P, Vetter T. Survival analysis and interpretation of time-to-event data: the tortoise and the hare. Anesth Analg. 2018;127:1. https://doi.org/10.1213/ANE.0000000000003653.
Slavík O, Horký P. When fish meet fish as determined by physiological sensors. Ecol Freshw Fish. 2009;18:501–6. https://doi.org/10.1111/j.1600-0633.2009.00370.x.
Slavík O, Horký P, Maciak M, Wackermannová M. Familiarity, prior residency, resource availability and body mass as predictors of the movement activity of the European catfish. J Ethol. 2016;34:23–30. https://doi.org/10.1007/s10164-015-0441-9.
Smith F. Understanding HPE in the VPS telemetry system. VEMCO Tutorials; 2013.
Snobl ZR, Koenigs RP, Bruch RM, Binkowski FP. Do tags exceeding 2% of total body weight impair lake sturgeon movement? N Am J Fish Manag. 2015;35:880–4. https://doi.org/10.1080/02755947.2015.1069425.
Souchon Y, Tissot L. Synthesis of thermal tolerances of the common freshwater fish species in large Western Europe rivers. Knowl Manag Aquat Ecosyst. 2012;405:03.
Sumpter DJ. Collective animal behavior. Princeton University Press; 2010.
Suski C, Ridgway M. Winter biology of centrarchid fishes. In: Cooke S, Philipp D, editors. Centrarchid fishes: diversity, biology, and conservation. Hoboken: Blackwell Publishing Ltd; 2009. p. 264–92.
Thayer D, Jonathan LW, Ruppert D, Watkinson TC, Poesch MS. Identifying temporal bottlenecks for the conservation of large-bodied fishes: lake sturgeon (Acipenser fulvescens) show highly restricted movement and habitat use over-winter. Glob Ecol Conserv. 2017;10:194–205. https://doi.org/10.1016/j.gecco.2017.03.008.
Therneau TM. Mixed effects Cox models [R package coxme version 2.2–17]; 2022a.
Therneau TM. A package for survival analysis in R. R package version 3.3–1; 2022b.
Tien JH, Levin SA, Rubenstein DI. Dynamics of fish shoals: identifying key decision rules. Evol Ecol Res. 2004;6:555–65.
Tobler M, Schlupp I. Influence of black spot disease on shoaling behaviour in female western mosquitofish, Gambusia affinis (Poeciliidae, Teleostei). Environ Biol Fishes. 2008;81:29–34. https://doi.org/10.1007/s10641-006-9153-x.
Vagnon C, Cattaneo F, Goulon C, Guillard J, Frossard V. The vulnerability of whitefish (Coregonus lavaretus) to the invasive European catfish (Silurus glanis) in a large peri-Alpine lake. Can J Fish Aquat Sci. 2022. https://doi.org/10.1139/cjfas-2021-0310.
Vanovac S, Howard D, Monk CT, Arlinghaus R, Giabbanelli PJ. Network analysis of intra-and interspecific freshwater fish interactions using year-around tracking. J R Soc Interface. 2021;18:20210445. https://doi.org/10.1098/rsif.2021.0445.
Vejřík L, Vejříková I, Blabolil P, Eloranta AP, Kocvara L, Peterka J, Sajdlova Z, Chung SHT, Smejkal M, Kiljunen M, Cech M. European catfish (Silurus glanis) as a freshwater apex predator drives ecosystem via its diet adaptability. Sci Rep. 2017. https://doi.org/10.1038/s41598-017-16169-9.
Ward AJW, Duff AJ, Krause J, Barber I. Shoaling behaviour of sticklebacks infected with the microsporidian parasite, Glugea anomala. Environ Biol Fishes. 2005;72:155–60. https://doi.org/10.1007/s10641-004-9078-1.
Ward AJW, Krause J. Body length assortative shoaling in the European minnow, Phoxinus phoxinus. Anim Behav. 2001;62:617–21. https://doi.org/10.1006/anbe.2001.1785.
Weber TC, Peña H, Jech JM. Consecutive acoustic observations of an Atlantic herring school in the Northwest Atlantic. ICES J Mar Sci. 2009;66:1270–7. https://doi.org/10.1093/icesjms/fsp090.
Webster MM, Hart PJB. Kleptoparasitic prey competition in shoaling fish: effects of familiarity and prey distribution. Behav Ecol. 2006;17:959–64. https://doi.org/10.1093/beheco/arl037.
Westrelin S, Boulêtreau S, Santoul F. European catfish Silurus glanis behaviour in response to a strong summer hypoxic event in a shallow lake. Aquat Ecol. 2022. https://doi.org/10.1007/s10452-022-09952-y.
Winter J. Advances in underwater biotelemetry. Fisheries techniques, 2nd edition American Fisheries Society, Bethesda, Maryland; 1996. p. 555–590.
Wolf JBW, Mawdsley D, Trillmich F, James R. Social structure in a colonial mammal: unravelling hidden structural layers and their foundations by network analysis. Anim Behav. 2007;74:1293–302. https://doi.org/10.1016/j.anbehav.2007.02.024.
Worton BJ. Kernel methods for estimating the utilization distribution in home-range studies. Ecology. 1989;70:164–8.
Zeileis A, Kleiber C, Krämer W, Hornik K. Testing and dating of structural changes in practice. Comput Stat Data Anal. 2003;44:109–23. https://doi.org/10.1016/S0167-9473(03)00030-6.
Zeileis A, Leisch F, Hornik K, Kleiber C. strucchange: an R package for testing for structural change in linear regression models. J Stat Softw. 2002;7:38. https://doi.org/10.18637/jss.v007.i02.
We are very grateful to Virginie Diouloufet, Julien Dublon, Dorian Milesi and Ange Molina for their investment in field work and to Yann Le Coarer and Nathalie Reynaud for dGPS treatments, including fieldwork support. We would like to thank Mélanie Debelgarric for her master internship which gave a very first overview of the topic and Léa Voisin, Maïlove Benoliel and numerous other people who occasionally helped in the field. We are very grateful to the two anonymous reviewers who made useful comments on this manuscript. We also thank Lionel Allègre, Stéphanie Bertrand and Dominique Ghione for their warm welcome and logistical support on the study site. And many thanks to Kurt Villsen for having checked English. This study has been conducted in accordance with the relevant animal ethics approvals. This project was in part funded by the “Département des Bouches-du-Rhône”.
This project was in part funded by the “Département des Bouches-du-Rhône”.
Ethics approval and consent to participate
The care and use of experimental animals complied with French animal welfare laws, guidelines and policies as approved by the French Ministry of Research through the authorisation number APAFIS#11294-2017091809143058 v2.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Synchronized tracks of European catfish over the four winters 2017 to 2020 in “Etang des Aulnes” https://doi.org/10.57745/U7UG5D.
Time series of the mean distance between individuals in winters 2017 to 2020. The dates of structural changes over the 5-month time series (15 October–15 March) and their 95% confidence interval are labelled on the x-axis and represented by vertical dotted lines and interval at their basis (very narrow intervals are not visible).
Temperature differences between the main aggregation zone and the deepest point in the lake. Time series are represented over winters 2017 to 2020 (a–d). The solid black line represents the mean daily temperature differences at 0.5 m above the bottom between locations 38 and 42 (see Fig. 1) over the 5-month time series (15 October–15 March). Location 38 is inside the aggregation zone that showed off each winter. Location 42 corresponds to the deepest point in the lake and stands as a reference point. The dates of aggregation given by the breakpoint detection algorithm are represented by vertical blue dotted lines. The period extended to the formation and dislocation of the aggregation is delimited by vertical blue solid lines.
Temperature differences between the secondary aggregation zone and the deepest point in the lake. Time series are represented over winter 2019. The solid black line represents the mean daily temperature differences at 0.5 m above the bottom between locations 5 and 42 (see Fig. 1) over the 5-month time series (15 October–15 March). Location 5 is very close to the secondary aggregation zone that showed off only in winter 2019. Location 42 corresponds to the deepest point in the lake and stands as a reference point. The dates of aggregation given by the breakpoint detection algorithm are represented by vertical blue dotted lines. The period extended to the formation and dislocation of the aggregation is delimited by vertical blue solid lines.
About this article
Cite this article
Westrelin, S., Moreau, M., Fourcassié, V. et al. Overwintering aggregation patterns of European catfish Silurus glanis. Mov Ecol 11, 9 (2023). https://doi.org/10.1186/s40462-023-00373-6
- Space use
- Fish behavior
- Silurus glanis