Lévy trajectories of Peruvian purse-seiners as an indicator of the spatial distribution of anchovy (Engraulis ringens)

Transcription

1 ICES Journal of Marine Science, 62: 477e482 (2005) doi: /j.icesjms Lévy trajectories of Peruvian purse-seiners as an indicator of the spatial distribution of anchovy (Engraulis ringens) Sophie Bertrand, Julian M. Burgos, François Gerlotto, and Jaime Atiquipa Bertrand, S., Burgos, J. M., Gerlotto, F., and Atiquipa, J Lévy trajectories of Peruvian purse-seiners as an indicator of the spatial distribution of anchovy (Engraulis ringens). e ICES Journal of Marine Science, 62: 477e482. Analogous to other top predators foraging on patchy resources, the spatial behaviour of fishers provides information on the spatial organization of fish. Focusing on the Peruvian anchovy purse-seine fishery, vessel monitoring system data are used to describe fishing vessels trajectories, and acoustic survey data to characterize anchovy spatial distribution. Vessel trajectories were analysed in terms of move-length distribution, whereas fish distribution was characterized by spatial extent, concentration, and fractal dimension. Fishers perform Lévy flights that can be characterized by a single statistic (m); m is significantly correlated with the fractal dimension of fish distribution. It is argued that the Lévy statistic is a good candidate for an ecosystem indicator that might contribute to realtime monitoring of ecosystems. Ó 2004 International Council for the Exploration of the Sea. Published by Elsevier Ltd. All rights reserved. Keywords: anchovy spatial distribution, fishing vessel trajectories, fractal dimension, Lévy flights, predatoreprey relationships, vessel monitoring system. Received 1 April 2004; accepted 3 November S. Bertrand and F. Gerlotto: IRD (Research Institute for Development) Centre de Recherche Halieutique Méditerranéenne et Tropicale, Avenue Jean Monnet, BP 171, Sète Cedex, France. J. M. Burgos: University of Washington, School of Aquatic and Fishery Sciences, 1122 NE Boat Street (Box ), Seattle, WA , USA. J. Atiquipa: IMARPE (Peruvian Institute of the Sea), Esquina Gamarra y Gral. Valles s/n, Chucuito, Callao, Lima, Peru. Correspondence to S. Bertrand: tel: C ; fax: C ; Sophie.Bertrand@ird.fr. Introduction Fishers are predators, so we can hypothesize that part of their spatial behaviour reflects the spatial distribution of fish. However, predatoreprey relationships are composed of interactions across different scales (Steele, 1989; Levin, 1992). For instance, the spatial aggregation of small pelagic fish is characterized by a nested hierarchical organization (Kotliar and Wiens, 1990; Wu and Loucks, 1995; Fauchald, 1999): the school, the cluster of schools, the cluster of clusters, and ultimately the population (Fréon and Misund, 1999; Petitgas et al., 2001). This hierarchy leads to two major inferences: (i) from a functional perspective, fishers have to develop a spatial strategy for dealing with the multi-scale aggregation pattern (Russel et al., 1992); (ii) from a methodological perspective, studying the process of interaction between fish and fishers implies use of metrics that allow multiple scales of potential interactions to be taken into account. Several studies have focused on spatial predatoreprey relationships in a multi-scale, hierarchical patch dynamic framework, using various metrics (Rose and Leggett, 1990; Fauchald et al., 2000; Davoren et al., 2002; Fauchald and Erikstad, 2002; Bertrand et al., 2004a). Frontier (1987) conjectured that if the search for prey is the main determinant for the distribution of organisms, there should be a type of behavioural cascade in space. Choosing the fractal dimension as a natural and simple method for describing objects across a range of scales (Kotliar and Wiens, 1990; Sugihara and May, 1990; Halley et al., 2004), Russel et al. (1992) argued, and to some extent showed, that under the behavioural cascade hypothesis, the fractal dimensions of predator and prey distributions are correlated. During the past decade, vessel monitoring systems (VMS) have been increasingly implemented in fisheries. In Peru for instance, the entire industrial fleet is equipped with satellite transmitters that record vessel positions. This /$30.00 Ó 2004 International Council for the Exploration of the Sea. Published by Elsevier Ltd. All rights reserved.

2 478 S. Bertrand et al. system produces real-time data with high spatial resolution, which has largely been used for law enforcement purposes. However, the information has great potential for monitoring exploited marine resources, and may be applied in an ecosystem context. In parallel, the Peruvian marine institute (IMARPE) usually performs between one and four acoustic surveys per year to assess biomass and spatial distribution of pelagic resources. Based on these spatially explicit data and within the conceptual framework of interacting hierarchical multiscale processes, our purpose here is to investigate whether fishing vessel trajectories can be used as an indicator of pelagic fish distribution. To test this hypothesis, we (i) characterize fishing vessel trajectories, (ii) characterize anchovy spatial distribution, and (iii) investigate if there is a relation between the spatial behaviour of fish and fishers. Material and methods Vessel monitoring system data Since 1999, a VMS (SISESAT: SIstema de SEguimiento SATelital) has been operational for the entire Peruvian industrial fishing fleet (about 800 vessels). Each vessel is equipped with an Argos satellite transmitter (CLS Argos, Toulouse, France) furnished with its own Global Positioning System (GPS) receiver. In principle, each transmitter sends its time-referenced position with an accuracy of 100 m on an hourly basis. When a GPS position (marked as class G) is not available, Argos provides a position with a lower level of accuracy (labeled class 1, 2, or 3). The data considered cover the period November 1999eMarch SISESAT transmits continuously, providing vessel positions without any distinction between positions at sea or in port. A first algorithm for preprocessing data aimed to identify positions at sea by excluding consecutive emissions located within a 2-nautical-mile radius of any port. Trips were defined as series of emissions originating at sea, bounded by emissions at port. We chose to work only with accurate data, retaining only trips for which all positions belonged to the G class and whose duration between consecutive emissions was %2 h. A second algorithm aimed to identify those trips that were most likely intended for fishing. Fishing trips usually last between 24 and 48 h, and cannot last much longer because the purse-seiners do not have freezing systems aboard. Therefore, we arbitrarily set maximum trip duration to 120 h, and rejected longer trips. Moreover, because fishing sets impose speeds!3 knots, we discarded trips for which speed was always greater. A third algorithm split the trajectory of each remaining fishing trip into elementary moves, moves being defined (Turchin, 1998) as the sum of consecutive positions of a vessel in the same direction. Instantaneous direction at sea is liable to be highly variable, so we arbitrarily defined a change of direction as when the angular difference between two consecutive pairs of positions exceeded 10(. Move lengths were calculated as orthodromic distances. A common method used to characterize animal trajectories (Turchin, 1998) is to describe the statistical distribution of move lengths between successive changes of direction. The probability density functions of these distributions can be formulated in a generalized form as an inverse power-law distribution (Viswanathan et al., 1996): PðlÞwl ÿm ð1þ where P(l) is the probability for a move length equal to l. Values of m are related to different types of redistribution kernel (Mandelbrot, 1977; Viswanathan et al., 1999): m R 3 corresponds to a Gaussian distribution or Brownian motion; m % 1 corresponds to probability distributions that cannot be normalized, but as m tends towards 1, the fraction of short moves becomes negligible and the forager moves along straight lines (i.e. close to a ballistic motion); 1! m! 3 defines the domain of long-tailed distributions named Lévy flights. This last type of motion has a high probability of large moves, and the variance (and even the mean) of the move-length distribution diverges with increasing number of moves, indicating that there is no characteristic scale for the observed distribution (scale invariance). The Lévy flight type of spatial behaviour has been documented for a variety of foraging natural predators: Drosophila (Cole, 1995), albatross (Viswanathan et al., 1996), reindeer (Mårell et al., 2002), jackal (Atkinson et al., 2002), and seals (Austin et al., 2004). To estimate m from VMS observations, first we built the histogram of the move-length frequency distribution of the fishing trips for each vessel and month (Figure 1a). Histogram shape, and any fitting procedure, is sensitive to the bin size. The most commonly applied rule for selecting the number of classes in constructing a histogram is that proposed by Sturges (1926). However, two problems arise with this method: (i) data are supposed to have a normal distribution, and (ii) the method tends to oversmooth the histogram. To overcome these limitations, histograms were constructed with Scott s (1979) rule for estimating bin width (h): hz3:5sn ÿ1=3 ð2þ where s is the standard deviation of the sample and n is total sample size. Finally, a least-squares linear regression was fitted through the tail (Viswanathan et al., 1999) of the loglog plot of the move-length frequency distribution (Figure 1b), the slope being an estimate of the inverse power-law exponent (m). Acoustic survey data IMARPE research vessels conducted 13 pelagic acoustic surveys between December 1999 and March 2003.

3 Lévy trajectories of Peruvian purse-seiners targeting anchovy 479 (a) Latitude (S) (b) Frequency(l) Longitude (W) Log 10 (Move length frequency) Log 10 (Move length) Move lengths (l) in km 120 Figure 1. (a) A fishing trip trajectory (black dots correspond to VMS positions), and (b) a move-length frequency distribution of the same vessel in April 2000 (histogram for 9 trips corresponding to 113 moves), with inset a log-log scale plot plus fitted regression line (m Z 2.01, r 2 Z 0.96, p! 0.001). Acoustic data were collected with a scientific echosounder system (38 and 120 khz split-beam SIMRAD EK500, Kongsberg SIMRAD AS, Norway), calibrated using standard procedures (Foote et al., 1987). The water column was sampled to depths of 250 and 500 m, for 120 and 38 khz, respectively. Standard survey design consisted of parallel transects running from 3.7 km from the coast to about 185 km offshore, with inter-transect distance varying between 26 and 30 km. Acoustic backscattered energy by surface unit (s A ) was recorded in each elementary sampling distance unit (ESDU) of 1852 m. Species composition was identified through trawling, echo-trace characteristics, and abiotic conditions. Data corresponding to anchovy were selected, and s A was considered to represent an index of anchovy local abundance. For each survey, anchovy spatial distribution was characterized by three metrics: (i) IMARPE routinely produces an estimate of the total surface occupied by anchovy (S in km 2 ); (ii) an index of spatial concentration (Ss; Petitgas, 1998; Bertrand et al., 2004b) was calculated as the difference between the experimental concentration profile of s A and a theoretical homogeneous spatial distribution; (iii) the fractal dimension of anchovy distribution was estimated using the method proposed by Bruno and Raspa (1989) and Halley et al. (2004). These authors showed that fractal dimension (D) is linearly related to the slope at the origin (a) of the log-transformed semivariogram (Matheron, 1965) for a pattern on a twodimensional landscape: DZ2 ÿ a=2 ð3þ For each survey, we computed the Napierian log-log semivariogram on presenceeabsence data for anchovy (Figure 2), and a was estimated by least-squares regression. The number of points included in the linear regression was chosen so as to maximize r 2. Results Raw data consisted of 41 months, 809 vessels, and emissions, corresponding to trips. Preprocessing of the data identified fishing trips ( moves) performed by 756 vessels during the 41 months studied. Mean m values among vessels by month ranged between 1.13 and 2.09 (Figure 3). Consequently, fishing trip trajectories are in the domain of Lévy flights motion. Months exhibiting particularly high variability of m correspond to the austral winter, when relatively few fishing trips were identified (Table 1). At that time of year, fishing closures are frequent, and the aberrant values may correspond to exploratory trips, or other atypical behaviour. Over the period studied, the total distribution area of anchovy (S) ranged between and km 2, spatial concentration (Ss) varied by a ratio of 13%, and fractal dimension (D) by a ratio of 4.7% (Table 1). Among the 13 surveys carried out, the five potentially coinciding with periods of fishing closure (threshold set at 100 fishing trips) Log (Variance (Distance)) Log (Distance) Figure 2. Napierian log-log transformation of the semi-variogram of the anchovy distribution for the survey October/November 2000 (x-axis corresponds to the distance between pairs of points for which variance is computed; slope-derived D Z 1.639, r 2 O 0.999).

4 480 S. Bertrand et al Dec99 Feb00 were excluded. Therefore, only eight surveys could be used to compare vessels trajectories and anchovy distribution metrics (Figure 4). Only the correlation between D and m was significant (Table 2). Discussion Apr00 Jun00 Aug00 Oct00 Dec00 Feb01 Jun01 Aug01 Oct01 Dec01 Feb02 Apr02 Jun02 Aug02 Oct02 Dec02 Figure 3. Time-series of estimated m by month (black squares: means; vertical bars: 0.95 confidence intervals; grey arrows indicate acoustic survey dates, dashed when corresponding to a period of fishing closure). The three metrics used to describe the spatial distribution of the fish were chosen as satisfying a trade-off between Table 1. Metrics describing anchovy distribution and vessel information by time period (S, stock range in km 2 ; Ss, spatial concentration, without dimension; D, fractal dimension based on presenceeabsence data, without dimension; m,lévy flights statistic of vessel trajectories, without dimension; N, number of vessels for which m was computed; n, number of representative fishing trips; information in parenthesis was not included in the correlation analyses because data refer to periods of potential fishing closures). Period Anchovy distribution Feb03 Vessel information S Ss D m N n Nov.eDec Jan.eFeb Jun.eJul (Aug.eSep. 2000) (23 911) (0.466) (1.650) (1.26) (12) (50) Oct.eNov Mar.eApr (Jul.eAug. 2001) (80 955) (0.434) (1.615) (1.94) (3) (18) (Aug.eSep. 2001) (42 244) (0.446) (1.652) (1.13) (16) (91) Oct.eNov Feb.eMar (Aug. 2002) (35 644) (0.456) (1.613) (1.89) (9) (45) Oct.eNov (Feb.eMar. 2002) (43 613) (0.452) (1.546) (1.71) (14) (91) 74.1 S Ss D Dec99 Feb00 Jun00 Oct00 available data and the need to describe the system over a continuum of scales. While spatial extent (S) has a direct interpretation, spatial concentration (Ss) and fractal dimension (D) have no intuitive interpretation. Ss can be taken as a measure of the way by which a given fish biomass fills the available space (Bertrand et al., 2004b), but it provides no information on the relative location of aggregations. Geometrically speaking, D can be seen as the degree of occupation of the physical space by a contorted or fragmented surface (Mandelbrot, 1977; Frontier, 1987), or as the degree of spatial irregularity of the analysed spatial function. Ecologically speaking, D is usually interpreted in terms of degree of patchiness or clustering of aggregations (Frontier, 1987). Results revealed that trajectories of fishers can be described as Lévy flights, but what does this mean? As m decreases, long moves are made longer, and short moves are made shorter (Mandelbrot, 1977). This has two related consequences: from a Lagrangian perspective, as m decreases, sinuosity of the trajectories decreases, whereas from an Eulerian perspective, as m decreases, diffusion of the individuals increases, so the area explored increases (Berkolaiko et al., 1996). Lévy flights are also characterized by their scale invariance (Mandelbrot, 1977). Following this idea, the spatial distribution of the points visited by a Lévy flyer (corresponding to the points of changing Apr01 Oct01 Mar02 Oct Figure 4. Time-series comparisons between the sinuosity of vessel trajectories (m, dashed line) and (a) anchovy spatial extent (S), (b) spatial concentration (Ss), and (c) fractal dimension (D). Table 2. Values of r 2 and of the F-statistic, and the associated probability (p) for the regression analysis between vessel trajectory metric (m) and fish distribution metrics (S, Ss, and D). Correlation r 2 F [1,6] p m w S n.s. m w Ss n.s. m w D !0.05

5 Lévy trajectories of Peruvian purse-seiners targeting anchovy 481 direction) has been shown to form a fractal with dimension D Z m ÿ 1 (Viswanathan et al., 1996; Atkinson et al., 2002). Mandelbrot (1977) indicated that as this fractal dimension (or m) decreases, clusters of visited points by the Lévy flyer become increasingly tightly packed, and thus increasingly separate from neighbouring clusters, and the number of apparent hierarchical levels seems to increase. In other words, m indicates the degree of clustering of the visited points. The significant correlation between D and m confirms the hypothesis that fishers spatial behaviour reflects to some degree the spatial distribution of fish. The absence of correlation between m and either S or Ss for the study period suggests that fishers spatial searching behaviour is mainly determined by the degree of clustering of fish. This, it is thought, is due to the process of information acquisition by fishers, and to their perception of the environment. Although significant, the correlation between D and m may appear weak. This may be explained by three factors: (i) the range of variation of D for the study period is about 4.7%, so although the linear regressions for estimating D had very high r 2 (O0.99), the estimated value is sensitive to the regression used; (ii) variability in behaviour among vessels may weaken the relationship; (iii) the pelagic ecosystem is characterized by strong interannual variability owing to El Niño events, which affect patchiness and connectivity of highly productive zones (Sugihara and May, 1990). Because there were no strong El Niño events during the period investigated, we do not think that we have seen the entire range of variability of the system. This would argue in favour of analysing longer time-series with this approach. We conclude that the m statistic of fishing vessel trajectories constitutes a good candidate as an ecosystem indicator because: (i) Lévy flights are intuitively and statistically appropriate for describing the spatial behaviour of fishers; (ii) they can easily be related to fractal geometry, and are appropriate for describing multi-scale processes; (iii) m is easy to compute and its interpretation is straightforward; (iv) the method can be applied in any ecosystem where VMS data are being collected; (v) owing to the ongoing development of VMS, its application should be cost-effective. If analysis of a longer time-series confirms the observed trend between m and D, m may be a promising indicator of the degree of clustering of fish aggregations that could contribute to real-time monitoring of ecosystems. Acknowledgements We thank IMARPE for providing the VMS and acoustic survey data, Arnaud Bertrand for his active participation during the study, and Pierrick Penven and François Colas for fruitful conversations about Lévy flights. This work is a contribution of IRD s Research Unit ACTIVE UR061 and the Interdepartemental Thematic Action Humboldt Current System (ATI SCH ) funded by IRD. References Atkinson, R. P. D., Rhodes, C. J., Macdonald, D. W., and Anderson, R. M Scale-free dynamics in the movement patterns of jackals. Oikos, 98: 134e140. Austin, D., Bowen, W. D., and McMillan, J. I Intraspecific variation in movement patterns: modelling individual behaviour in a large marine predator. Oikos, 105: 15e30. Berkolaiko, G., Havlin, S., Larralde, H., and Weiss, G. H The expected number of distinct sites visited by N Lévy flights on a one-dimensional lattice. Physical Review E, 53: 5774e5778. Bertrand, A., Barbieri, M. A., Córdova, J., Hernández, C., Gómez, F., and Leiva, F. 2004a. Diel vertical behaviour, predatoreprey relationships and occupation of space by jack mackerel (Trachurus murphyi) off Chile. ICES Journal of Marine Science, 61: 1105e1112. Bertrand, S., Díaz, E., and Ñiquen, M. 2004b. Interactions between fish and fisher s spatial distribution and behaviour: an empirical study of the anchovy (Engraulis ringens) fishery of Peru. ICES Journal of Marine Science, 61: 1127e1136. Bruno, R., and Raspa, G Geostatistical characterization of fractal models of surfaces. Geostatistics, 1: 77e89. Cole, B. J Fractal time in animal behaviour: the movement activity of Drosophila. Animal Behaviour, 50: 1317e1324. Davoren, G. K., Montevecchi, W. A., and Anderson, J. T Scale-dependent associations of predators and prey: constraints imposed by flightlessness of common murres. Marine Ecology Progress Series, 245: 259e272. Fauchald, P Foraging in a hierarchical patch system. The American Naturalist, 153: 603e613. Fauchald, P., and Erikstad, K. E Scale-dependent predatore prey interactions: the aggregative response of seabirds to prey under variable prey abundance and patchiness. Marine Ecology Progress Series, 231: 279e291. Fauchald, P., Erikstad, K. E., and Skarsfjord, H Scaledependent predatoreprey interactions: the hierarchical spatial distribution of seabirds and preys. Ecology, 81: 773e783. Foote, K. G., Knudsen, H. P., Vestnes, D. N., MacLennan, D. N., and Simmonds, E. J Calibration of acoustic instruments for fish density estimation: a practical guide. ICES Cooperative Research Report, pp. Fréon, P., and Misund, O. A Dynamics of Pelagic Fish Distribution and Behaviour: Effects on Fisheries and Stock Assessment. Fishing News Books, London. 348 pp. Frontier, S Applications of fractal theory to ecology. In Developments in Numerical Ecology, pp. 335e378. Ed. by P. Legendre, and L. Legendre. Springer, Berlin. NATO ASISeries, 14. Halley, J. M., Hartley, S., Kallimanis, A. S., Kunin, W. E., Lennon, J. J., and Sgardelis, S. P Uses and abuses of fractal methodology in ecology. Ecology Letters, 7: 254e271. Kotliar, N. B., and Wiens, J. A Multiple scales of patchiness and patch structure: a hierarchical framework for the study of heterogeneity. Oikos, 59: 253e260. Levin, S. A The problem of pattern and scale in ecology. Ecology, 73: 1943e1967. Mandelbrot, B. B Fractals, Form, Chance and Dimension. W. H. Freeman and Co., San Francisco. 365 pp. Mårell, A., Ball, J. P., and Hofgaard, A Foraging and movement paths of female reindeer: insights from fractal analysis, correlated random walks, and Lévy flights. Canadian Journal of Zoology, 80: 854e865. Matheron, G La Théorie des Variables Régionalisées et Ses Applications. Masson et Cie, Paris. 306 pp. Petitgas, P Biomass-dependent dynamics of fish spatial distributions characterized by geostatistical aggregation curves. ICES Journal of Marine Science, 55: 443e453. Petitgas, P., Reid, D., Carrera, P., Iglesias, M., Georgakarakos, S., Liorzou, B., and Massé, J On the relation between

6 482 S. Bertrand et al. schools, clusters of schools, and abundance in pelagic fish stocks. ICES Journal of Marine Science, 58: 1150e1160. Rose, G. A., and Leggett, W. C The importance of scale to predatoreprey spatial correlations: an example of Atlantic fishes. Ecology, 71: 33e43. Russel, R. W., Hunt, G. L., O Coyle, K., and Cooney, R. T Foraging in a fractal environment: spatial patterns in a marine predatoreprey system. Landscape Ecology, 7: 195e209. Scott, D. W On optimal and data-based histograms. Biometrika, 66: 605e610. Steele, J. H The ocean landscape. Landscape Ecology, 3: 185e192. Sturges, H. A The choice of a class interval. Journal of the American Statistical Association, 21: 65e66. Sugihara, G., and May, R. M Applications of fractals in ecology. Trends in Ecology and Evolution, 5: 79e86. Turchin, P Quantitative Analysis of Movement: Measuring and Modelling Population Redistribution in Plants and Animals. Sinauer Associates, Sunderland, MA. 336 pp. Viswanathan, G. M., Afanasyev, V., Buldyrev, S. V., Murphy, E. J., Prince, P. A., and Stanley, H. E Lévy flight search patterns of wandering albatrosses. Nature, 381: 413e415. Viswanathan, G. M., Buldyrev, S. V., Havlin, S., da Luz, M. G. E., Raposo, E. P., and Stanley, H. E Optimizing the success of random searches. Nature, 401: 911e914. Wu, J., and Loucks, O. L From balance of nature to hierarchical patch dynamics: a paradigm shift in ecology. The Quarterly Review of Biology, 70: 439e466.