Next Article in Journal
Hydrochemical Characteristics and Hydrogeochemical Simulation Research of Groundwater in the Guohe River Basin (Henan Section)
Next Article in Special Issue
Managing Water Level for Large Migratory Fish at the Poyang Lake Outlet: Implications Based on Habitat Suitability and Connectivity
Previous Article in Journal
Water Leakage and Crack Identification in Tunnels Based on Transfer-Learning and Convolutional Neural Networks
Previous Article in Special Issue
Effect of Subsurface Mediterranean Water Eddies on Sound Propagation Using ROMS Output and the Bellhop Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Short-Range Movement Pattern of Amphidromous Lagoon Fish Schools: Ecological Applications

1
Institut de Recherche pour le Développement, IRD, Univ Brest, CNRS, Ifremer, Lemar, Dakar BP 1386, Senegal
2
MARBEC, Univ Montpellier, CNRS, Ifremer, IRD, 34 203 Sète, France
3
Institut de Recherche pour le Développement, IRD, 34 203 Sète, France
4
Muséum National d’Histoire Naturelle, MNHN, IRD, Délégation Régionale Ouest, 29 280 Plouzané, France
5
Universidad Tecnológica de Chetumal, UT, Chetumal 77000, Quintana Roo, Mexico
6
Univ. Savoie Mont Blanc, INRAE, Alpine Center for Research on Trophic Networks and Limnic Ecosystems CARRTEL, 74200 Thonon-les-Bains, France
*
Author to whom correspondence should be addressed.
Water 2022, 14(9), 1463; https://doi.org/10.3390/w14091463
Submission received: 7 March 2022 / Revised: 26 April 2022 / Accepted: 28 April 2022 / Published: 3 May 2022
(This article belongs to the Special Issue Hydroacoustics in Marine, Transitional and Freshwaters)

Abstract

:
High-resolution multibeam sonar allows estimating movements of pelagic fish schools at short range. Taking advantage of this methodology, we calculated a Straightness Index (SI) to quantify the proportion of schools migrating actively from those residents in lagoon channels. This information enhances our knowledge of both fish school displacements and migration processes, which are essential to improve our understanding of ecosystem functioning. Most fish schools (65%) exhibited a SI value demonstrating oriented swimming behavior through the channels displayed by schools reaching the sea during fall migration. This trend appears as an intrinsic property of school movements, allowing monitoring of the school migration process in a channel to provide information for manager vs. fishing regulation measures or lagoon planning. The result strengthens the ‘multi-transit’ hypothesis, as 35% of schools show sinuous trajectories representative of schools staying in the channel or displaying high exploratory behaviors. Lastly, the fish school Exploration Swimming Speed (ESS) was tested as a fishery-independent sampling method to evaluate the proportion of different fish species monitored using hydroacoustics. This approach demonstrates the interest in using swimming behavioral characteristics of fish schools for ecological and management purposes.

1. Introduction

Fish movements are involved in a wide range of behavioral processes such as migration, space use, food searching, and reproduction [1,2]. Multibeam sonars in fisheries science [3] are usually applied to pelagic fish schools to scrutinize their tridimensional morphology [4,5,6]. Multibeam sonar observations in shallow waters allow the measurement of fish school swimming speeds and movements at different temporal scales [7,8,9], as well as calculation of the straightness index (SI) [1], which corresponds to the ratio between net and gross distances patrolled by a school. This index aims to characterize the path sinuosity of a school for given biological, ecological and social conditions. Moreover, the SI value could be an indicator of the efficiency for a fish school to orient its movement to reach a goal (e.g., [1,10,11]). Indeed, in a patchy food environment, foraging behaviors are often characterized by an increased turning rate and, therefore, by sinuous tracks, whereas migration behaviors should be more efficient when straight line movements are used [12,13,14]. In this work, we scrutinized the swimming patterns under the assumption that exploratory and migration behaviors show different trajectories, i.e., zigzag vs. straight line, respectively. Then, the path sinuosity of schools inferred from SI values is discussed with reference to the school migration behavior in lagoon channels and its consequences on school study methodologies using acoustic remote systems.
Otherwise, the fish swimming speed is often a key factor in ecological studies (e.g., [15,16,17]) and several approaches to measure this metric were experimentally developed in laboratories, water tanks and swim tunnels [18]. In situ, the swimming speed is usually estimated from electronic and ultrasonic tags on individual fish [19]. Studies of path orientation on individual fish were carried out using this technology [20]. For studying fish school velocity, one approach is to use multibeam acoustic systems [8,21,22] but very few studies are available in the literature, particularly at short range. Due to hydrodynamic constraints, the sizes of fish are homogeneous in a school [23]. Under this assumption, fish sizes by species obtained from both fishing data and the Exploratory Swimming Speed (ESS) [9] of each school estimated by the multibeam sonar tracking method were used to explore the relationship between school velocity and species. Thus, we proposed a new method to estimate the relative abundance of different groups of schooling species in a lagoon.

2. Materials and Methods

Two 24-h acoustic surveys were carried out consecutively in early autumn inside two Mediterranean lagoon channels in 1999 (Figure 1). The study area was located in the south of France (43°44′ N; 03°79′ E and 43°52′ N; 03°90′ E) in two coastal shallow water lagoons (Ingril, 549 ha and Prévost, 380 ha) connected to the sea by a channel [8]. These two lagoons are part of a series of shallow ponds with similar fish species communities [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]. The common species of these lagoons are Anguilla anguilla, Atherina boyeri, Dicentrarchus labrax, Sparus aurata, and various mullets (Liza ramada, Liza aurata, Mugil cephalus, Chelon labrosus) [8,25,26]. Only schooling species were considered in this work. During the acoustic surveys, cast net samplings were performed in the channel for fish identification [8]. During this sampling period, the dominant gregarious fish groups were represented by juveniles of mugilidae (fork length ‘FL1’~7.5 cm), of S. aurata and D. labrax (mean FL2~13 cm), and the last group was composed of sub-adults of S. aurata, D. labrax and of several mugilidae species (mean FL3~20 cm) [8], identified by local experts. Considering the ethical rules on animal welfare, no dissection or experiment was carried out on fish. These species encountered in lagoons are known to migrate to the sea in autumn for physiological (e.g., decrease in water lagoon temperature) and spawning reasons, whereas the spring migration to lagoons is considered to be linked to trophic motivations, lagoons being considered as a nursery for juveniles of several species [27].
A Reson Seabat 6012 multibeam sonar, operating at a frequency of 455 kHz with a pulse duration of 0.06 ms, was deployed at a fixed position on the channel bank. The sonar was used in horizontal beaming [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28] across both lagoon channels, width less than 25 m (20 m for the Ingril channel and 24 m for the Prévost channel) and 1.5 m depth (Appendix A). The efficient horizontal angle of echo reception is 90° (horizontal plan: 60 contiguous 1.5° beams; vertical angle of 15°), the sonar resolution was around 45 cm [29]. The sonar detection, i.e., echo traces, allows observing the kinematics of fish schools [8] over both 24-h surveys. Fish schools showing splitting/merging events were removed from the analysis to avoid double counting. The total time of continuous school detection (residence time) is variable due to the loss of school echo traces. Indeed, the swimming trajectory inside the sonar beams generates the variation in the observation time between fish schools.
The lowest elementary linear distances (l in meter) were estimated between two consecutive positions of the geometric center of fish school echo traces on a time interval of 2 s between consecutive position records. The sum ‘L’ of elementary ‘l’ distances for a given school corresponds to the gross displacement or path length of this school. Finally, the beeline distance ‘D’ (in meters) was measured [9] from the starting point (i.e., first detection) to the last position of the school in the sonar beams (net displacement). Fish school path lengths L and beeline distances D were estimated for each school with a minimum of three echo traces [9] reported from sonar detections.
To compare the gross to net path lengths, we calculated the straightness index (SI) [1], which corresponds to the ratio between the net and the gross distances patrolled by a fish school and computed as the ratio between D and L (SI = D/L). This index ranging between 0 and 1 is a reliable measure of the sinuosity of the path, i.e., a SI close to 1 illustrates directed movements linked to the efficiency of an orientation mechanism to reach a goal [1,30]. Moreover, the Exploratory Swimming Speed (ESS in m s−1) of each school was estimated in [9] by dividing the beeline distance D between the starting point to the last point of the school in the beam with the time interval between these two extreme records of the school in the sonar beams. Previous processing on the same data set have shown that the distance values between fish school echo traces were recorded for time intervals ranging from 1 to 8 s (average of 4 s) [9]. ESS mean values were 1.19 m s−1 (standard deviation: σ = 0.77) and 1.34 m s−1 (σ = 0.79) for the Ingril and Prévost channels, respectively; no difference was found between the Ingril and Prévost channels [9]. These Exploratory Swimming Speed (ESS) data were compared to the fish lengths by species obtained from fishing data [8] in order to explore the relationship between school velocity and species.
The difference in beeline distances (D), path length (L), Straightness Index (SI) and Exploratory Swimming Speed (ESS) of schools between both lagoons were tested using analysis of variance (ANOVA) tests. All the assumptions for the models were previously checked. If the assumptions were not met, a non-parametric Kruskal–Wallis test was used. Linear regression models were used to analyze the relationships between the straightness index and the total time of the school observation or the path length L in the sonar beams. All statistical tests were performed with R software (R Core Team, 2021; version 3.6.2, https://www.R-project.org/ (accessed on 16 July 2021)).

3. Results

3.1. Fish School Data Collection

A total of 164 fish schools was analyzed over both 24-h surveys [9], 41 originated from the Ingril channel and 123 from the Prévost one. The number of sampled school echo traces reached 621 (n = 174 in the Ingril channel and n = 447 in the Prévost channel). The residence time in the acoustic beam varies from 2 to 34 s (average of 10 s for both channels) [9]. No statistical difference (ANOVA, F(1,162) = 0.8469, p-value = 0.359) was found for the residence time between both channels.
The beeline value (D) of schools ranged between 0.34 and 22.1 m (average = 8.59 m; SD = 5.04) for the Ingril channel and between 1.57 to 37.1 m (average = 10.7 m; SD = 5.62) for the Prévost channel [9]. A statistical difference was found for the beeline distances (D) between channels (ANOVA, F(1, 162) = 4.37, p-value = 0.03815). However, no significant difference was found for the path length between channels (ANOVA, F (1, 162) = 2.811, p-value = 0.09553) (Figure 2b).

3.1.1. Straightness Index

For the Ingril and Prévost channels, SI values for schools ranged from 0.31 to 1 (average = 0.87, SD = 0.20) and from 0.12 to 1 (average = 0.89, SD = 0.16), respectively. No statistical difference for the SI was found between the Ingril and Prévost channels (Kruskal–Wallis, 0.50632, df = 1, p-value = 0.4767) (Figure 2c). Similar cumulated frequency distributions of SI values were found for both lagoons (Appendix C; Figure A4).
A linear significant relationship was found between the SI and the total time of the school observation in the sonar beams for the Prevost channel (p-value = 0.0127) but with a low adjusted R-squared (R2 = 0.04). No significant linear relationship was found for the Ingril channel (p-value = 0.162). Finally, a linear significant relationship was found between the SI and the path length (L) for the Prevost channel (R2 = 0.25, p-value = 3.36 × 10−9). No significant linear relationship was found for the Ingril channel (p-value = 0.273) (Figure 3).

3.1.2. Exploratory Swimming Speed (ESS) of Fish Schools

The scatterplot of ESS values versus the residence time showed the effect of a limiting factor (Figure 4). This limiting factor suggests maximal speed values of schools for an oriented path according to the residence time in the beams and regarding the fish community in concern. The theoretical limits of the observations given by the sampling protocol was a power trend curve (Figure 4).
According to the average sizes of the three groups of migratory gregarious fish sampled by fishing (FL1~7.5 cm, FL2~13 cm and FL3~20 cm), ESS values were divided into three groups (Figure 5) with an upper limit of ESS by group corresponding to a body length speed of about 19 Bl s−1. This limit value of 19 Bl s−1 was selected because, for a school with both a short residence time in the beams and an oriented path (i.e., major part of the observations in this study), ESS might be considered close to instantaneous speed values. Considering this limit value, individual fish swimming speeds for the three groups correspond to ESS values of 1.4 m s−1, 2.5 m s−1 and 3.8 m s−1, respectively.

4. Discussion

We used a high-resolution multibeam sonar in horizontal beaming to better understand the movements of pelagic fish schools during the migration period in shallow waters. Our approach demonstrates the interest in using swimming behavioral characteristics of fish schools for ecological and management purposes.

4.1. Straigtness Index and Fish Migration

If we assume that the oriented swimming behavior corresponds to a certain form of an active migration (continuous swimming activity), our results estimated that ~65% of fish schools (Appendix B, Figure A3) exhibited an active migration movement through the lagoon channels. Up to a threshold value of 0.9 for SI, the school was considered displaying an oriented swimming behavior. A significant linear relationship was observed between the path length (L) and the total time in the sonar beam for a given school. Rather than being dependent of one of the two previous variables considered, this movement pattern appeared as an intrinsic characteristic of the observed schools since such a general trend of school movements in both channels was observed regardless of the path length (Figure 3). Furthermore, as the time of the school observation in the beam depends essentially on the distance between the school and the beam-emitting source (the sonar transducer), this rules out the hypothesis of a likely relationship between SI values and the residence time of schools in the beam. The factor observed in Figure 4 represents the maximum distance covered by the sonar beam, knowing that it covers a maximum distance of 20 and 24 m perpendicular to the channels (the width of the lagoon channel at Ingrid and Prevost, respectively) and therefore (with a sonar opening angle of 90°) a maximum distance (parallel to the channel) between 40 and 48 m. Hence, the general pattern observed for the orientation of school displacements in lagoon channels would be independent of the sampling protocol.
This estimation is reliable under the hypothesis that schools exhibited this movement pattern throughout their transfer within the lagoon channels during the fall migration of mugilidae, sparidae and centrarchidae. Fish schools not having a well-defined migration behavior regarding their SI value are thus susceptible to be counted several times by the acoustic system (what has been called the multi-transit hypothesis [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30]) and thus overestimate the school passage in the channel. They reached ~35% of analyzed fish schools, which could represent resident schools in lagoon channels having an exploratory behavior rather than a migratory one.
No statistical difference was found for the SI values between the two channels. Hence, we cannot conclude about differences in swimming tracks of the fish species. Nevertheless, the beeline distance was significantly higher for the Prevost channel than for the Ingril channel. The path length was also higher for Prevost than for the Ingril channel although no statistical difference was found between the two channels, which is consistent with the fact that the SI were close to 1. Even if similar fish species communities were present in both the Prévost and Ingril channels [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24], this could be related to a change in abundance of fish species. Indeed, during the sampling period, the dominant gregarious fish species were represented by S. aurata and D. labrax and several mugilidae species. According to [8], we observed that among species that could be detected by the multibeam sonar, D. Labrax was in higher abundance in the catches during the fall survey (according to catches per unit of effort values, Figure A4). This result is consistent with the fact that the global school shapes were significantly different between the lagoons as shown in [8]. Furthermore, based on these data, 92.4% of the fish could display a migration behavior in fall for both channels, which is consistent with the SI values. However, further investigations would be required to establish a link between the fish species discriminations and the sonar detections.
We found according to the beeline distance that fish schools from both lagoons did not exhibit the same swimming behavior (Figure 2) in the lagoon channels. Such difference can be linked to the difference observed in fish species composition in the two lagoons (Appendix C, Figure A4). An alternative explanation for the difference in swimming behaviors, as seen for the beeline distance of fish schools from both lagoons could be the presence of resident fish in the lagoon channels. Indeed, we could assume that migratory fish schools (lagoon towards sea at the studied season) get a straighter displacement than the resident ones, which should lead to more exploration of the channels for e.g., trophic reasons. The absence of difference in the SI could be explained by the fact that the fall migration started in both channels. Nevertheless, additional information on the migration of each species is necessary to state such a hypothesis.

4.2. Do Exploratory Swimming Speed Values Can Help to Discriminate Fish Species?

Relative abundances of each size group in the lagoon channels were inferred from ESS distributions. On this basis, the proportion of each group (G1, G2 and G3) will be 60, 30 and 10% in the Prévost lagoon channel and 68, 20 and 12% in the Ingril one, respectively. Obviously, further surveys must be carried out to assess the relevance of this remote sensing approach, which is an interesting independent-fishery sampling method of gregarious fish communities in shallow-waters, i.e., an alternative method to usual lethal fish sampling like gillnets or fyke nets. Monitoring fish communities in sensitive zones like Marine Protected Areas or essential habitats using such methodology is obviously relevant. In combination with the e-DNA method [31], underwater video observations [32,33,34], such as complementary non-invasive approaches, might be used in the future for fish monitoring to gather diversity, quantification and behavioral knowledge.
To gather our data, the operating system carried out was time consuming (Appendix A, Figure A1). An automation of working sequences using dedicated software is suitable. As echo traces were easily identifiable, a discrimination algorithm of useful echoes might be technically easy to develop [35,36,37].

5. Conclusions

Combining fish school path sinuosity and exploratory swimming speed will bring interesting perspectives in ecological studies of pelagic fish schools. First, this information enhances our knowledge of both fish school displacements and migration processes, which are essential to improve our understanding of ecosystem functioning and therefore its modelling at the mesoscale [38,39]. The method is non-invasive and allows monitoring of the school migration process in a channel to provide information for manager vs. fishing regulation measures or lagoon planning. Second, the size structure of the fish community can be inferred from this approach forward as a fishery-independent sampling method. Thus, it will allow helping decision making for management measures at the seashore. Third, this kind of information remains rare in the literature but can be applied for other purposes, e.g., to estimate the fish school coefficient of diffusion in situ at a fine scale. Both straightness index (SI) and Exploratory Swimming Speed (ESS) should be proposed as indicators for various purposes: (i) to discriminate fish school species, (ii) characterize their behavioral motivation (feeding, spawning, and migration), (iii) identify the structure of the fish community into size groups, and (iv) even bring elementary information on fish activity rates for fish bioenergetic models.

Author Contributions

Conceptualization, P.B. (Patrice Brehmer); methodology, P.B. (Patrice Brehmer) and M.S.; software, P.B. (Patrice Brehmer); validation, P.B. (Pascal Bach) and J.G.; formal analysis, P.B. (Patrice Brehmer), V.D. and P.I.C.P.; investigation, P.B. (Patrice Brehmer); resources, P.B. (Patrice Brehmer); data curation, P.I.C.P.; writing—original draft preparation, P.B. (Patrice Brehmer) and M.S.; writing—review and editing, P.B. (Pascal Bach), P.B. (Patrice Brehmer), V.D., N.D. and J.G.; supervision, P.B. (Patrice Brehmer) and P.B. (Pascal Bach); project administration, P.B. (Patrice Brehmer); funding acquisition, P.B. (Patrice Brehmer) and M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by planning agreement between the national and regional governments CPER XI (Occitanie, France), and ended within the AWA project funded by French National Research Institute for Sustainable Development “IRD” and the German Federal Ministry of Education and Research BMBF grant number 01DG12073E and the European Maritime and Fisheries Fund “Acapela”.

Data Availability Statement

Not applicable.

Acknowledgments

We dedicate this work to Thang Do Chi † (Univ Montpellier). Thank to Simon Benhamou (CNRS) for early discussion, to David Mouillot (Univ Montpellier) for friendly work within this project, François Gerlotto (retired Orstom-IRD) for his early work on the use of multibeam sonar in shallow water to monitor fish school and Pascal Cotel (retired Orstom-IRD) as well as all volunteers for support during the field measurements and the French National Research Institute for Sustainable Development ‘IRD’ staff of IRD DR Ouest and DR Occitanie.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Fish school dynamic observations using high-resolution multibeam sonar (according to [8,9].
The RESON Seabat 6012 multi-beam sonar used for data acquisition emits on 60 contiguous beams of 1.5° each. The efficient horizontal angle of echo receptions is 90° with a vertical angle of 15°. The sonar frequency was 455 kHz with pulse duration of 0.06 ms; all the data were continuously stored on video recording supports. The sonar characteristics and the environmental parameters determine the threshold of the sonar resolution, in our case 45 cm [29]. A preliminary study of acoustic data, intended to quantify the migratory fish school flows collected from the school echo traces [40,41] count from acoustic imagery, was performed using the same acoustic equipment in horizontal beaming [8].
The video recordings were replayed at the laboratory to select the sonar sequences including fish school echo traces (Figure A1), which correspond to specific detections of homogeneous continuous responses clearly discriminated on the screen. For both sampling areas, we were able to observe both mobile and stationary echo traces [9]. Dynamic echo traces and characteristics of fish school detections were discriminated from fixed bottom echo traces (Figure A2). In this way, each selected series of sonar images corresponding to a detection of a school identified by an individual code was stored in a fish school library. Each separate fish school datum was extracted from the sonar images using the ‘Infobancs’ software [37]. For each fish school, we collected the number of consecutive echo traces ‘N’, the total time of observation of echo traces within the beam (in seconds) and the Euclidian position (x; y) of the center of the fish school, defined as the gravity center of the surface defining the detected biological structure. From this information and scale factor of observations on the screen [8] we estimated all distances traveled by the fish schools. The selected time interval was set at one second for the shortest observation, without restriction in the total time of observation above three seconds.
Figure A1. Scheme representing the sonar data collection, their treatment, which include several steps (selection of sonar sequence, digitalization, identification of echo traces on sonar images, data extraction and then exportation for final analysis on ad hoc software), and their analysis to obtain the swimming speed measurements.
Figure A1. Scheme representing the sonar data collection, their treatment, which include several steps (selection of sonar sequence, digitalization, identification of echo traces on sonar images, data extraction and then exportation for final analysis on ad hoc software), and their analysis to obtain the swimming speed measurements.
Water 14 01463 g0a1
Figure A2. Diagram of the experimental set up, showing the sonar beams covering the channel, bank to bank. On the right, an example of a sonar picture where the dynamic echoes of fish school (inside the circle, in t1, t2 and t3) were discriminated from the fixed echoes of the bottom (Reprinted with permission from [9]).
Figure A2. Diagram of the experimental set up, showing the sonar beams covering the channel, bank to bank. On the right, an example of a sonar picture where the dynamic echoes of fish school (inside the circle, in t1, t2 and t3) were discriminated from the fixed echoes of the bottom (Reprinted with permission from [9]).
Water 14 01463 g0a2

Appendix B

Figure A3. Cumulative frequency of the fish school net to gross displacement ratio ‘straightness index’; the schools having a SI above 0.9 (~65% of the schools) were considered to exhibit an oriented swimming behavior that corresponds to a certain form of active migration. In white, the Ingril lagoon schools, and in black the Prevost ones.
Figure A3. Cumulative frequency of the fish school net to gross displacement ratio ‘straightness index’; the schools having a SI above 0.9 (~65% of the schools) were considered to exhibit an oriented swimming behavior that corresponds to a certain form of active migration. In white, the Ingril lagoon schools, and in black the Prevost ones.
Water 14 01463 g0a3

Appendix C

Figure A4. Composition of 15 fish species commercially caught in the two Prévost and Ingril lagoons. Data extracted from [8].
Figure A4. Composition of 15 fish species commercially caught in the two Prévost and Ingril lagoons. Data extracted from [8].
Water 14 01463 g0a4

Appendix D

Table A1. Raw experimental data used for this study and extracted from [8,9], the Straightness index (SI) was the Beeline value divided (D) by the Path Length (L), calculated as SI = D/L. Below the theoretical function of limiting factor to the sampling protocol.
Table A1. Raw experimental data used for this study and extracted from [8,9], the Straightness index (SI) was the Beeline value divided (D) by the Path Length (L), calculated as SI = D/L. Below the theoretical function of limiting factor to the sampling protocol.
LagoonsResidence Time (s)ESS (m s−1)Beeline Value (D) (m)Path Length (L) (m)Straigthness Index (SI)
Ingril20.7771.551.570.988
Ingril21.0232.052.060.994
Ingril22.2754.554.570.996
Ingril32.6157.848.790.893
Ingril31.9685.906.050.976
Ingril31.2663.803.850.987
Ingril31.8975.695.750.989
Ingril32.7028.118.140.995
Ingril32.7188.158.161.000
Ingril41.6576.636.631.000
Ingril41.7206.887.160.961
Ingril50.1490.742.120.351
Ingril51.0225.117.430.688
Ingril51.1065.535.810.952
Ingril50.5922.962.961.000
Ingril52.00910.0410.250.980
Ingril50.4862.433.430.707
Ingril62.48414.9114.911.000
Ingril70.0480.341.080.313
Ingril70.8786.146.141.000
Ingril71.0647.4517.201.000
Ingril80.7856.287.250.867
Ingril82.58020.6422.090.934
Ingril82.21517.7217.721.000
Ingril90.8367.537.580.993
Ingril90.8327.497.491.000
Ingril91.57814.2017.510.811
Ingril111.07611.8313.240.894
Ingril110.92810.2110.211.000
Ingril120.7368.8427.920.316
Ingril130.5256.826.880.992
Ingril160.85413.6616.490.828
Ingril161.38322.1222.121.000
Ingril170.4988.4716.050.528
Ingril180.60510.9011.260.968
Ingril190.5249.9610.150.981
Ingril210.63913.4213.620.985
Ingril220.3637.9916.410.487
Ingril230.45010.3512.700.815
Ingril230.47310.8811.570.941
Ingril310.52316.2222.380.725
Prevost23.0686.146.440.953
Prevost22.1874.374.580.955
Prevost21.8863.773.820.988
Prevost32.2596.788.230.824
Prevost31.5154.544.960.915
Prevost31.8515.555.730.969
Prevost32.8548.568.561.000
Prevost31.0013.003.001.000
Prevost31.6194.864.861.000
Prevost33.0129.049.041.000
Prevost33.63710.9111.130.980
Prevost40.3911.562.760.568
Prevost40.7893.155.030.627
Prevost41.0934.376.780.645
Prevost40.4931.972.630.748
Prevost40.8193.283.950.829
Prevost43.77415.1017.860.845
Prevost43.69214.7715.900.928
Prevost40.5152.062.200.938
Prevost41.8017.207.670.939
Prevost41.0464.194.370.958
Prevost41.7006.807.050.964
Prevost41.9987.998.060.991
Prevost41.0434.174.200.993
Prevost43.08312.3312.410.993
Prevost41.7216.886.910.996
Prevost42.59610.3810.420.997
Prevost40.5802.322.530.916
Prevost50.5642.823.750.753
Prevost51.4267.137.510.949
Prevost50.9584.794.900.978
Prevost50.9064.534.630.980
Prevost51.9689.849.990.986
Prevost51.9049.529.640.987
Prevost51.2356.186.220.993
Prevost52.50612.5312.191.000
Prevost51.3876.939.040.767
Prevost62.26913.6114.200.959
Prevost61.1096.656.750.986
Prevost61.94211.6511.820.986
Prevost62.33714.0214.050.998
Prevost61.67610.0610.061.000
Prevost61.5579.349.341.000
Prevost63.25619.5419.541.000
Prevost60.8264.964.961.000
Prevost71.66311.6420.410.570
Prevost71.60211.2118.570.604
Prevost71.87013.0920.280.646
Prevost71.89813.2819.620.677
Prevost71.91413.4018.760.714
Prevost71.51210.5814.820.714
Prevost70.8265.786.810.850
Prevost70.6854.805.270.911
Prevost71.2448.719.070.959
Prevost71.83412.8415.490.829
Prevost71.73912.1713.440.905
Prevost72.26515.8617.280.918
Prevost71.2208.548.790.971
Prevost71.98713.9113.970.996
Prevost82.09616.7719.200.873
Prevost80.6104.884.920.992
Prevost80.4783.823.840.996
Prevost81.0648.518.540.997
Prevost81.0528.428.430.998
Prevost81.38511.0822.300.497
Prevost81.89515.1615.950.950
Prevost90.5184.665.330.874
Prevost90.3823.443.770.914
Prevost90.6936.246.460.966
Prevost91.37012.3312.560.982
Prevost92.61023.4927.380.858
Prevost91.56014.0414.800.949
Prevost92.06118.5519.000.976
Prevost101.69316.9317.200.985
Prevost100.4504.504.520.994
Prevost101.30013.0013.001.000
Prevost101.48914.8917.790.837
Prevost101.94219.4220.280.957
Prevost101.25312.5312.910.970
Prevost100.9419.419.420.999
Prevost101.33013.3013.301.000
Prevost110.5726.297.530.836
Prevost110.8579.439.660.976
Prevost121.07812.9314.700.879
Prevost120.8299.9510.050.990
Prevost120.7739.2710.630.872
Prevost121.72520.7020.701.000
Prevost121.06712.8112.811.000
Prevost130.6828.8710.090.879
Prevost131.09314.2014.370.988
Prevost130.7049.159.210.994
Prevost130.3494.546.120.742
Prevost130.88311.4813.940.823
Prevost130.81610.6112.200.870
Prevost131.36217.7127.670.640
Prevost141.05414.7514.741.000
Prevost141.40919.7319.780.997
Prevost150.70010.5010.800.972
Prevost150.96714.5115.600.930
Prevost151.11216.6818.490.902
Prevost160.6129.7813.110.747
Prevost160.4917.857.861.000
Prevost160.78412.5516.090.780
Prevost160.81513.0513.700.953
Prevost161.20819.3313.341.000
Prevost160.88114.1016.810.839
Prevost170.95116.1716.191.000
Prevost180.85315.3516.020.958
Prevost182.06337.1342.240.879
Prevost191.12421.3621.970.972
Prevost200.2454.917.260.676
Prevost210.2835.9450.830.117
Prevost210.2294.8137.740.127
Prevost220.82318.1124.290.746
Prevost250.64416.1016.111.000
Prevost270.74620.1320.101.000
Prevost270.67218.1420.850.870
Prevost290.54915.9216.090.990
Prevost330.41913.8315.860.872
Prevost330.48215.8917.280.920
Prevost340.30110.2427.840.368
Prevost340.2418.1910.470.782
Prevost340.49716.9018.110.933
Table A2. Assuming that fish can only be tracked over a maximum distance of 40 m (e.g., fish passing along the bank opposite to the sonar bank, i.e., where it was immersed) and setting an obvious limit constraint at ESS < 5m s−1 (above which swimming speed was not biologically realistic in our case study) and a Residence time ≥ 1 s, the theoretical function is as follows (maximum distance (m) of the sampling protocol divided by the residence time (s)).
Table A2. Assuming that fish can only be tracked over a maximum distance of 40 m (e.g., fish passing along the bank opposite to the sonar bank, i.e., where it was immersed) and setting an obvious limit constraint at ESS < 5m s−1 (above which swimming speed was not biologically realistic in our case study) and a Residence time ≥ 1 s, the theoretical function is as follows (maximum distance (m) of the sampling protocol divided by the residence time (s)).
Residence
time (s)
246810121416182022242628303234363840
ESS
(m s−1)
20.010.06.75.04.03.32.92.52.22.01.81.71.51.41.31.31.21.11.11.0

References

  1. Benhamou, S. How to reliably estimate the tortuosity of an animal’s path: Straightness, sinuosity, or fractal dimension? J. Theor. Biol. 2004, 229, 209–220. [Google Scholar] [CrossRef] [PubMed]
  2. Benhamou, S. Detecting an orientation component in animal paths when the preferred direction is individual-dependent. Ecology 2006, 87, 518–528. [Google Scholar] [CrossRef] [PubMed]
  3. Simmonds, E.J.; MacLennan, D.N. Fisheries Acoustics: Theory and Practice; Blackwell Publishing: London, UK, 2005. [Google Scholar]
  4. Brehmer, P.; Gerlotto, F.; Rouault, A. In Situ Inter-Standardization of Acoustics Data: An Integrated Database for Fish School Behaviour Studies. Acta Acust. United Acust. 2002, 88, 730–734. [Google Scholar]
  5. Guillard, J.; Fernandes, P.; Laloë, T.; Brehmer, P. Three-dimensional internal spatial structure of young-of-the-year pelagic freshwater fish provides evidence for the identification of fish school species. Limnol. Oceanogr. Methods 2011, 9, 322–328. [Google Scholar] [CrossRef] [Green Version]
  6. David, V.; Mouget, A.; Perrot, Y.; Le Goff, L.; Thiriet, P.; Diogoul, N.; Feunteun, E.; Acou, A.; Brehmer, P. Insights from a multibeam echosounder to survey pelagic fish shoals and their spatio-temporal distribution in ultra-shallow waters. Estuar. Coast. Shelf Sci. 2022, 264, 107705. [Google Scholar] [CrossRef]
  7. Brehmer, P. Fisheries Acoustics: Theory and Practice, 2nd edn. Fish Fish. 2006, 7, 227–228. [Google Scholar] [CrossRef]
  8. Brehmer, P.; Mouillot, D.; Chi, T.D. Amphidromus fish school diel flow in two Mediterranean lagoons by combining sonar and fishing data. J. Exp. Mar. Biol. Ecol. 2006, 334, 139–150. [Google Scholar] [CrossRef]
  9. Brehmer, P.; Guillard, J.; Pinzon, P.I.C.; Bach, P. Exploratory and Instantaneous Swimming Speeds of Amphidromous Fish School in Shallow-Water Coastal Lagoon Channels. Estuaries Coasts 2011, 34, 739–744. [Google Scholar] [CrossRef] [Green Version]
  10. Stasko, A.B.; Horrall, R.M.; Hasler, A. Coastal Movements of Adult Fraser River Sockeye Salmon (Oncorhynchus nerka) Observed by Ultrasonic Tracking. Trans. Am. Fish. Soc. 1976, 105, 64–71. [Google Scholar] [CrossRef]
  11. Hawkes, L.A.; Broderick, A.C.; Coyne, M.S.; Godfrey, M.H.; Lopez-Jurado, L.-P.; Lopez-Suarez, P.; Merino, S.E.; Varo-Cruz, N.; Godley, B.J. Phenotypically Linked Dichotomy in Sea Turtle Foraging Requires Multiple Conservation Approaches. Curr. Biol. 2006, 16, 990–995. [Google Scholar] [CrossRef] [Green Version]
  12. Weimerskirch, H.; Bonadonna, F.; Bailleul, F.; Mabille, G.; Dell’Omo, G.; Lipp, H.P. GPS tracking of foraging albatrosses. Science 2002, 295, 1259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Giacalone, V.; D’Anna, G.; Pipitone, C.; Badalamenti, F. Movements and residence time of spiny lobsters, Palinurus elephas released in a marine protected area: An investigation by ultrasonic telemetry. J. Mar. Biol. Assoc. U. K. 2006, 86, 1101–1106. [Google Scholar] [CrossRef]
  14. Klimley, A.P.; Beavers, S.C.; Curtis, T.C.; Jorgensen, S.J. Movements and Swimming Behavior of Three Species of Sharks in La Jolla Canyon, California. Environ. Biol. Fishes 2002, 63, 117–135. [Google Scholar] [CrossRef]
  15. Ware, D.M. Bioenergetics of Pelagic Fish: Theoretical Change in Swimming Speed and Ration with Body Size. J. l’Office Rech. Pêcheries Can. 1978, 35, 220–228. [Google Scholar] [CrossRef]
  16. Brochier, T.; Auger, P.-A.; Pecquerie, L.; Machu, E.; Capet, X.; Thiaw, M.; Mbaye, B.C.; Braham, C.-B.; Ettahiri, O.; Charouki, N.; et al. Complex small pelagic fish population patterns arising from individual behavioral responses to their environment. Prog. Oceanogr. 2018, 164, 12–27. [Google Scholar] [CrossRef]
  17. Videler, J.J. Fish Swimming; Chapman & Hall: London, UK, 1993. [Google Scholar]
  18. Marcinek, D.J.; Blackwell, S.B.; Dewar, H.; Freund, E.V.; Farwell, C.; Dau, D.; Seitz, A.C.; Block, B.A. Depth and muscle temperature of Pacific bluefin tuna examined with acoustic and pop-up satellite archival tags. Mar. Biol. 2001, 138, 869–885. [Google Scholar] [CrossRef]
  19. Zydlewski, G.B.; Haro, A.; Whalen, K.G.; McCormick, S.D. Performance of stationary and portable passive transponder detection systems for monitoring of fish movements. J. Fish Biol. 2001, 58, 1471–1475. [Google Scholar] [CrossRef]
  20. Dunlop, E.S.; Milne, S.W.; Ridgway, M.S.; Condiotty, J.; Higginbottom, I. In Situ Swimming Behavior of Lake Trout Observed Using Integrated Multibeam Acoustics and Biotelemetry. Trans. Am. Fish. Soc. 2010, 139, 420–432. [Google Scholar] [CrossRef]
  21. Trygonis, V.; Georgakarakos, S.; Dagorn, L.; Brehmer, P. Spatiotemporal distribution of fish schools around drifting fish aggregating devices. Fish. Res. 2016, 177, 39–49. [Google Scholar] [CrossRef]
  22. Fréon, P.; Misund, O.A. Dynamics of Pelagic Fish Distribution and Behaviour: Effects on Fisheries and Stock Assessment; Blackwell Science: London, UK, 1999. [Google Scholar]
  23. Mouillot, D.; Laune, J.; Tomasini, J.A.; Aliaume, C.; Brehmer, P.; Dutrieux, E.; Chi, T.D. Assessment of coastal lagoon quality with taxonomic diversity indices of fish, zoobenthos and macrophyte communities. Hydrobiologia 2005, 550, 121–130. [Google Scholar] [CrossRef]
  24. Bourquard, C. Structure et Mécanisme de Mise en Place, de Maintien et d’Evolution des Peuplements Ichthyologiques Lagunaires du Golfe du Lion. Ph.D. Thesis, Université Montpellier, Montpellier, France, 1985. [Google Scholar]
  25. Bach, P.; Legendre, P.; Amanieu, M.; Lasserre, G. Strategy of eel (Anguilla anguilla, L.) exploitation in the Thau lagoon. Estuar. Coast. Shelf Sci. 1992, 35, 55–73. [Google Scholar] [CrossRef]
  26. Quignard, J.P.; Mazoyer, C.; Vianet, R.; Wai, R.M.; Benharrat, K. Un exemple d’exploitation lagunaire en Languedoc: L’étang de l’Or (Mauguio)—Pêche et production halieutique. Sci. Pêche 1983, 336, 3–23. [Google Scholar]
  27. Brehmer, P.; Chi, T.D.; Laugier, T.; Galgani, F.; Laloë, F.; Darnaude, A.M.; Fiandrino, A.; Mouillot, D. Field investigations and multi-indicators for shallow water lagoon management: Perspective for societal benefit. Aquat. Conserv. 2011, 21, 728–742. [Google Scholar] [CrossRef] [Green Version]
  28. Brehmer, P.; Vercelli, C.; Gerlotto, F.; Sanguinède, F.; Pichot, Y.; Buestel, D.; Guénnegan, Y. Multibeam sonar three-dimensional monitoring of mussel culture grounds in open sea for management purpose. Aquaculture 2006, 252, 234–241. [Google Scholar] [CrossRef] [Green Version]
  29. Epstein, N. On tortuosity and the tortuosity factor in flow and diffusion through porous media. Chem. Eng. Sci. 1989, 44, 777–779. [Google Scholar] [CrossRef]
  30. Cronkite, G.; Mulligan, T.; Holmes, J.; Enzenhofer, H. Categorising salmon migration behaviour using characteristics of split-beam acoustic data. Aquat. Living Resour. 2007, 20, 205–212. [Google Scholar] [CrossRef]
  31. Thalinger, B.; Wolf, E.; Traugott, M.; Wanzenböck, J. Monitoring spawning migrations of potamodromous fish species via eDNA. Sci. Rep. 2019, 9, 15388. [Google Scholar] [CrossRef] [Green Version]
  32. Brehmer, P.; Soncho, G.; Trygornis, V.; Itano, D.; Dalen, J.; Fuchs, A.; Faraj, A.; Taquet, M. Toward an autonomous observatory of pelagic environment: Experiences from fish communities monitoring around drifting FADs. Thalass. Int. J. Mar. Sci. 2019, 35, 177–189. [Google Scholar] [CrossRef] [Green Version]
  33. Minart, C.; David, V.; Mouget, A.; Brehmer, P.; Acou, A.; Le Goff, L.; Feunteun, E.; Thiriet, P. An innovative sampling protocol for fish species identification methods in shallow waters: Towed diver, towed video and stereoscopic camera system. In Proceedings of the Oceans 2021, San Diego, CA, USA, 20–23 September 2021; pp. 1–10. [Google Scholar] [CrossRef]
  34. Gladstone, W.; Lindfield, S.; Coleman, M.; Kelaher, B. Optimisation of baited remote underwater video sampling designs for estuarine fish assemblages. J. Exp. Mar. Biol. Ecol. 2012, 429, 28–35. [Google Scholar] [CrossRef]
  35. Weill, A.; Scalabrin, C.; Diner, N. MOVIES-B: An acoustic detection description software. Application to shoal species’ classification. Aquat. Living Resour. 1993, 6, 255–267. [Google Scholar] [CrossRef]
  36. Perrot, Y.; Brehmer, P.; Habasque, J.; Roudaut, G.; Behagle, N.; Sarre, A.; Lebourges Dhaussy, A. Matecho: An Open-Source Tool for Processing Fisheries Acoustics Data. Acoust. Aust. 2018, 46, 241–248. [Google Scholar] [CrossRef]
  37. Brehmer, P.; Lafont, T.; Georgakarakos, S.; Josse, E.; Gerlotto, F.; Collet, C. Omnidirectional multibeam sonar monitoring: Applications in fisheries science. Fish Fish. 2006, 7, 165–179. [Google Scholar] [CrossRef]
  38. Cury, P.M.; Shin, Y.-J.; Planque, B.; Durant, J.M.; Fromentin, J.M.; Kramer-Schadt, S.; Stenseth, N.C.; Travers, M.; Grimm, V. Ecosystem oceanography for global change in fisheries. Trends Ecol. Evol. 2008, 23, 338–346. [Google Scholar] [CrossRef] [PubMed]
  39. Fulton, E.A.; Link, J.S.; Kaplan, I.C.; Savina-Rolland, M.; Johnson, P.; Ainsworth, C.; Horne, P.; Gorton, R.; Gamble, R.J.; Smith, A.D.M.; et al. Lessons in modelling and management of marine ecosystems: The Atlantis experience. Fish Fish. 2011, 12, 171–188. [Google Scholar] [CrossRef]
  40. Moreno, G.; Josse, E.; Brehmer, P.; Nøttestad, L. Echotrace classification and spatial distribution of pelagic fish aggregations around drifting fish aggregating devices (DFAD). Aquat. Living Resour. 2007, 20, 343–356. [Google Scholar] [CrossRef]
  41. Olsen, K. A note on estimating school size from echotraces. FAO Fish. Rep. 1969, 78, 37–48. [Google Scholar]
Figure 1. (A1) Map of the study area: the two French Mediterranean shallow water lagoons are located in the same watershed; they are part of a series of shallow ponds, and are both linked to the sea by a single channel. (A2,A3) On the right, an aerial view of the two lagoons (source: Google earth), with the channels (circle and zoom) where the sonar system was set up.
Figure 1. (A1) Map of the study area: the two French Mediterranean shallow water lagoons are located in the same watershed; they are part of a series of shallow ponds, and are both linked to the sea by a single channel. (A2,A3) On the right, an aerial view of the two lagoons (source: Google earth), with the channels (circle and zoom) where the sonar system was set up.
Water 14 01463 g001
Figure 2. (a) The beeline value (D) of fish schools for the Ingril and Prévost channels. A significant difference (*) was found between D values for channels. (b) However, no statistical difference was found for the path length (L) between the Ingril and Prévost channels. (c) Straightness index (SI) for the Ingril and Prévost channels, where no significant difference was found.
Figure 2. (a) The beeline value (D) of fish schools for the Ingril and Prévost channels. A significant difference (*) was found between D values for channels. (b) However, no statistical difference was found for the path length (L) between the Ingril and Prévost channels. (c) Straightness index (SI) for the Ingril and Prévost channels, where no significant difference was found.
Water 14 01463 g002
Figure 3. Relationship between the net to gross displacement ratio ‘straightness index’ and the fish school path length (Full triangles: the Ingril lagoon. Empty circles: the Prévost lagoon). The black line represents the single significant linear relationship found for the Prévost lagoon.
Figure 3. Relationship between the net to gross displacement ratio ‘straightness index’ and the fish school path length (Full triangles: the Ingril lagoon. Empty circles: the Prévost lagoon). The black line represents the single significant linear relationship found for the Prévost lagoon.
Water 14 01463 g003
Figure 4. Relationship between the fish school residence time (in seconds) and their exploratory swimming speed (ESS in m s−1). The scatterplot of ESS values versus the residence time showed the effect of a limiting factor, highlighted by the bold curve, which represents the theoretical limits (bold point) of the observations given by the sampling protocol, here a power trend curve of the theoretical function (see, Appendix D).
Figure 4. Relationship between the fish school residence time (in seconds) and their exploratory swimming speed (ESS in m s−1). The scatterplot of ESS values versus the residence time showed the effect of a limiting factor, highlighted by the bold curve, which represents the theoretical limits (bold point) of the observations given by the sampling protocol, here a power trend curve of the theoretical function (see, Appendix D).
Water 14 01463 g004
Figure 5. Histogram of average Exploratory Swimming Speed of fish schools observed in the Prevost (black) and Ingril (white) lagoons. Three groups (G1; G2 and G3) might be distinguished at a similar swimming speed interval (~1.25 m s−1), which were assumed to be related to main migratory gregarious species: mugilidae juveniles (G1: Fork length ‘FL’ < 7.5 cm), juveniles of Sparus aurata and Dicentrarchus labrax (G2: FL~13 cm), and G3 composed of sub-adults of several species (mullets, S. aurata and D. labrax) with an average FL size around 20 cm.
Figure 5. Histogram of average Exploratory Swimming Speed of fish schools observed in the Prevost (black) and Ingril (white) lagoons. Three groups (G1; G2 and G3) might be distinguished at a similar swimming speed interval (~1.25 m s−1), which were assumed to be related to main migratory gregarious species: mugilidae juveniles (G1: Fork length ‘FL’ < 7.5 cm), juveniles of Sparus aurata and Dicentrarchus labrax (G2: FL~13 cm), and G3 composed of sub-adults of several species (mullets, S. aurata and D. labrax) with an average FL size around 20 cm.
Water 14 01463 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brehmer, P.; Soria, M.; David, V.; Pinzon, P.I.C.; Bach, P.; Diogoul, N.; Guillard, J. Short-Range Movement Pattern of Amphidromous Lagoon Fish Schools: Ecological Applications. Water 2022, 14, 1463. https://doi.org/10.3390/w14091463

AMA Style

Brehmer P, Soria M, David V, Pinzon PIC, Bach P, Diogoul N, Guillard J. Short-Range Movement Pattern of Amphidromous Lagoon Fish Schools: Ecological Applications. Water. 2022; 14(9):1463. https://doi.org/10.3390/w14091463

Chicago/Turabian Style

Brehmer, Patrice, Marc Soria, Viviane David, Pablo Ivan Caballero Pinzon, Pascal Bach, Ndague Diogoul, and Jean Guillard. 2022. "Short-Range Movement Pattern of Amphidromous Lagoon Fish Schools: Ecological Applications" Water 14, no. 9: 1463. https://doi.org/10.3390/w14091463

APA Style

Brehmer, P., Soria, M., David, V., Pinzon, P. I. C., Bach, P., Diogoul, N., & Guillard, J. (2022). Short-Range Movement Pattern of Amphidromous Lagoon Fish Schools: Ecological Applications. Water, 14(9), 1463. https://doi.org/10.3390/w14091463

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop