Estimating tag-shedding rates for skipjack tuna, Katsuwonus pelamis, off the Maldives

Estimating tag-shedding rates for skipjack tuna, Katsuwonus pelamis, off the Maldives – Statistical Data Included

M. Shiham Adam

One source of uncertainty in fishery assessments based on tag release and recapture data arises from tag shedding–the loss of tags from fish from the time of tagging until tag recovery. Independent estimation of tag-shedding rates from double tagging experiments is an integral part of well-designed tagging experiments. Failure to allow for tag shedding can result in biases in estimates of important parameters derived from tag-recapture data, such as fishery-induced and natural mortality, and migration rates. A wide variety of methods have been proposed for estimating shedding rates using data from double tagging experiments, the variety in part resulting from the nature of the data available from the experiments. The earlier literature has been well integrated and reviewed by Wetherall (1982). More recently, Xiao (1996) has developed a general model that unifies the estimation of tag-shedding rates from double tagging experiments with exact and pooled times at liberty. Some recent examples of application of these methods include Xiao et al. (1999) for school gummy shark and Fabrizio et al. (1999) for lake trout.

Tag shedding is of two types (Wetherall, 1982). Type-I shedding is a one-time event and occurs immediately after tagging, usually as a result of suboptimal placement of tags in the fish. Effectively, it reduces the number of tags initially put out to sea. Type-II shedding is the loss of a tag or both tags over a period of time after the fish has been tagged and released back into the sea. For long-lived species, it may not occur at a constant rate because some tags are likely to have been applied more effectively than others, and some may become firmly embedded (with growth of muscle tissues), such that they are very unlikely to be detached from the fish (Kirkwood, 1981).

The Ministry of Fisheries and Agriculture, Maldives, carried out two tagging experiments: one in 1990-91 (Yesaki and Waheed, 1992) and the other between 1993 and 1995 (Anderson et al., 1996). The latter tagging experiment included a double tagging experiment in which 504 skipjack tuna, Katsuwonus pelamis, were double tagged, by using the same type of tags and techniques of tagging as used in the single tagging experiment (Anderson et al., 1996).As of end 1996, a total of 53 of these fish had been recovered. These data are considered to be of immediate importance for designing large-scale experiments in the Indian Ocean, strategic tagging simulation studies (e.g. Bertignac, 1996), and for comparing the estimated rates with those obtained from similar tuna double tagging experiments conducted elsewhere. Our note reports estimates of tag-shedding rates carried out in the Maldives.

Materials and methods

Tagging methods

Tagging was conducted on board local vessels fishing for skipjack tuna with pole-and-line gear using livebait. Plastic dart tags (100 mm x 1.5 mm diameter), manufactured by Hallprint[TM], Australia, were used throughout the experiment. The captured fish were gently placed on deck, quickly slipped onto a wet wooden measuring board, and held in place by biologists wearing cotton gloves. The first tag was inserted dorsally on the left-hand side, at an acute angle adjacent to the second dorsal fin so that the barb was caught under the fin-ray extension or the neural spine. The second tag was inserted about 1-2 cm posterior to the first on the right-hand side in the same manner. Consecutively numbered pairs of tags were used; even numbers were inserted on the left-hand side and odd numbers on the right. Where possible, fish were returned to the water, facing the vessel’s bow, in a slightly headdown fashion. Tagging times (from hooking to release into the sea) ranged from 14 to 18 seconds. More details of the tagging program can be found in Yesaki and Waheed (1992) and Anderson et al. (1996).

Parameter estimation

The method of parameter estimation used here is the maximum likelihood approach introduced by Kirkwood and Walker (1984) and later extended by Hampton and Kirkwood (1990) and Hampton (1997). This method was developed for use with sets of data where recaptured fish are few, and with few fish that have shed a tag, provided that exact dates of recovery are known, as was the case with the data from the Maldives.

Following Kirkwood and Walker (1984), consider an originally single tagged fish. Assuming that [Alpha] is the type-I retention probability (1-type-I shedding probability) and [Lambda] is the type-II shedding rate, assumed to be constant for short-lived species such as the skipjack tuna, then for an originally single tagged fish, the probability of a tag being retained at time t is given by

(1) Q(t) = [Alpha] exp (-[Lambda]t).

Suppose that fish are double tagged with identical tags and released at time t = 0, and let [p.sub.i](t) be the probability that the fish is alive and at liberty retaining i (i=0, 1, 2) tags at time t. Then, under the assumption that both tags are retained after immediate shedding and have independent and identical probabilities,

(2) [p.sub.0](t): [1 – Q[(t)].sup.2]

[p.sub.1](t) = 2Q(t)[1 – Q(t)]

[p.sub.2](t) = Q[(t).sup.2].

In practice, identifiable recaptures will consist only of fish retaining either one tag or two tags. If [P’.sub.i](t) is the probability that an originally double-tagged fish, recaptured at time [t.sub.i], is reported to have retained i tags (i=1, 2), then conditional on the retention of at least one tag, the probability of capturing a fish retaining two tags at time t is

(3) [P’.sub.2](t) = [p.sub.2](t)/1 – [p.sub.0](t)

and the probability of capturing a fish retaining only one tag at time t is

(4) [P’.sub.1](t) = [p.sub.1](t)/1 – [p.sub.0](t).

Suppose n fish were recaptured in the experiment and reported to have retained at least one tag on recapture, and that the ith fish was recaptured at time ti. Define indicator variables [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] if two tags were reported upon recapture, and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] if only one tag was reported upon recapture.

It follows from Equations 3 and 4 that if all fish suffer the same risks of mortality, the log likelihood ([Psi]) of the data conditional on recapture times {[t.sub.i], i:1, 2, … , n} is given by

(5) [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Maximum likelihood estimates of the parameters ([Alpha], [Lambda]) and their asymptotic standard errors were found by maximizing [Psi] with respect to them, by using the nonlinear minimization routines in AD Model Builder (Otter Research Ltd., 1996).

Results

The numbers of recoveries of originally double tagged fish reported as retaining one (indicating left- or the right-side tag) or two tags on recapture are shown in Table 1, with their times at liberty.

Table 1 Numbers of recoveries of originally double tagged fish retaining one (indicating which tag, left or right) or two tags on recapture, and times at liberty.

Days at

liberty 1 tag 2 tags

6 1

7 1

8 1

9 1

10 2

11 1 (L) 1

13 1

14 3

15 1

16 2

17 1

18 1 (L)

19 1

20 2

21 1

23 1

36 1

50 1

52 1

59 2

60 1 (R)

61 1

63 1

66 1

71 1

73 1

78 1

90 1

91 2

98 1

142 1 (R)

146 1

154 1

191 1

210 1

230 1

240 1

258 1 (L)

279 1 (R)

335 1

341 1

433 1 (L)

The maximum likelihood estimate of [Lambda] was 0.22/yr (SE=0.13) and that of a was 0.97 (SE=0.03). If [Alpha] 1, the maximum likelihood estimate of [Lambda] was 0.30/yr (SE:0.065/ yr). A likelihood ratio test (Cox and Hinkley, 1974) shows that the full model does not provide a significantly better fit to the data (at 5% level) than its special case ([Alpha]=1) (P:0.146).

Discussion

One of the most important assumptions in double-tagging experiments and in the model used in our analysis is that shedding rates of the first and second tags is the same (Hearn et al., 1991). A simple way to investigate possible differences in the shedding rates of the first and second tags is to examine the number of returns of fish that have retained either the first or the second tag and their temporal distribution. There were seven recoveries of fish reported to have retained a single tag, of which four fish retained the left tag and three retained the right tag (Table 1). Although few, the similar numbers of recaptured fish with tags on the left side as those on the right side and the temporal distribution of these fish (Table 1) show clearly that there is no evidence of differences in their shedding rates.

Another important assumption, but one that might easily be violated, and that is difficult to test, relates to the way in which double-tagged recoveries were reported. The commonassumption(asmadehere)isthatrecoveriesoffishretaining double tags were always reported as a pair, and never as a “single” tag recovery. It is also assumed that the probability of reporting fish recovered with single or double tags is the same. Hampton (1997) showed that under this assumption, the reporting probability has no influence on the maximum likelihood estimates of [Alpha] and [Lambda]. Given the publicity and incentives to return all the recaptured tags (Anderson et al., 1996), and the procedures adopted in tagging, it is highly unlikely that these reporting-rate assumptions were violated in the Maldives tagging program.

In their original paper describing the method of analysis and application to small data sets, Kirkwood and Walker (1984) noted that potential bias arises because the resulting estimates are conditional on the times of recapture; a different time sequence of recapture times would result in different parameter estimates. Clearly, this potential bias arose in our experiment. Ideally, for robust estimation of the parameters, one would need a data set with large numbers of fish that have shed a tag. This point needs to be borne in mind when designing future tagging experiments for skipjack tuna in the Indian Ocean.

The type-II shedding rate estimate obtained from our experiment is similar to those obtained by Hampton and Kirkwood (1990) for southern bluefin tuna tagged and released during the 1960s and 1970s. However, recent estimates of type-II shedding rates obtained for skipjack tuna from the tropical western Pacific were much lower than those estimated in our study, as were estimates for southern bluefin tuna tagged in 1984 (Table 2). These differences in shedding estimates may be due to differences in the tagging method used and the use of measuring boards instead of cradles. Estimates of [Lambda] obtained for fish tagged on measuring boards ranged from 0.17 to 0.22/yr. For fish tagged in cradles, the corresponding range of estimates of [Lambda] was 0.05-0.09/yr (Table 2).

[TABULAR DATA 2 NOT REPRODUCIBLE IN ASCII]

In the South Pacific Commission (SPC) double tagging experiments of 1989-92 (Hampton, 1997) and in the CSIRO experiment of 1984 (Hampton and Kirkwood, 1990), specially designed vinyl cradles were used to position the fish dorsal side up, making it easier to insert the tags. It is also likely that such cradles may have helped to keep fish calm, perhaps because the eyes of the fish were covered. In the Maldives experiment, the CSIRO experiments of 1963-77 and the Western Australian Department of Fisheries’ experiment of 1970-78 (Hampton and Kirkwood, 1990), measuring boards were used. In many instances in the Maldivian experiment, the fish, particularly large individuals, struggled quite a lot while being tagged. This may have caused additional trauma to the fish and at the same time may have increased the likelihood of the tag being placed suboptimally, thereby increasing the apparent tag-shedding rate. Some fish handlers of the Maldivian tagging program observed that fish struggled less if they were gently laid on the board and the eye was quickly covered with one hand. This technique, however, was not applied universally.

Hampton and Kirkwood (1990) suggested that differences in shedding rates observed in the southern bluefin tuna experiments might also be due to differences in the quality of the tags used. Tags used in the earlier experiments (1960s and 1970s) were inferior, because the streamers were prone to detach. Tags used in later experiments had their heads heat-fused to the shafts and thus were almost impossible to detach under normal conditions (Hampton, 1997). The tags used in the Maldivian experiment were recently manufactured and therefore this source of bias was unlikely. The higher tag shedding rate observed in the Maldivian data appears most likely to be due to the tagging method employed.

Ideally, experiment-specific double tagging experiments should be carried out to estimate independent experiment-specific tag-shedding rates. The above comparisons suggest that, in cases where independent estimates are not available, care must be exercised when applying estimates of tag-shedding rates obtained from other tagging experiments, especially where different tagging techniques have been used.

Acknowledgments

The authors would like to thank Charles Anderson, Graham Pilling, and Robert Wakeford for their comments on the manuscript. In addition to providing the data, the staff of Marine Research Center, particularly Zaha Waheed and Ahmed Haifz, helped on numerous queries regarding the data. We thank the Hon. Hassan Sobir, former Minster of Fisheries and Agriculture, Maizan Hassan Maniku, Director General Marine Research Center, for their support in this study. The manuscript was improved by comments from two anonymous reviewers. The study was part of the Ph.D. research of the senior author at Imperial College, London, primarily funded by the Islamic Development Bank’s Merit Scholarship Program.

Literature cited

Anderson, R. C., M. S. Adam, and A. Waheed. 1996. Tuna tagging activities in the Maldives. In Proceedings of the expert consultations on Indian Ocean tunas, 6th session, Colombo, Sri Lanka, 25-29 September 1995 (A. A. Anganuzzi, K. A. Stobberup, and N. J. Webb, eds.), p. 333-347. [Reprinted with minor revisions in Mald. Mar. Res. Bull. 2:83-116.]

Bertignac, M. 1996. A simulation model of tagging experiments for yellowfin tuna in the western Indian Ocean. In Status of interactions of Pacific tuna fisheries in 1995 (R. S. Shomura, J. Majkowski, and R. F. Harman, eds.), p. 162-177. FAO Fish. Tech. Paper 365.

Cox, D. R., and D. V. Hinkley. 1974. Theoretical statistics. Chapman and Hall, London, 511 p.

Fabrizio, M.C., J. D. Nicholas, J. E. Hines, B. L. Swanson, and S. T. Schram. 1999. Modelling data from double-tagging experiments to estimate heterogeneous rates of tag shedding in lake trout (Salvelinus mamaycush). Can. J. Aquat. Fish. Sci. 56(8):1409-1419.

Hampton, J. 1997. Estimation of tag reporting and tag-shedding rates in a large-scale tuna tagging experiment in the western tropical Pacific Ocean. Fish. Bull. 95:68-79.

Hampton, J., and G. P. Kirkwood. 1990. Tag shedding by southern bluefin tuna Thunnus maccoyii. Fish. Bull. 88:313-321.

Hearn, W. S., G. M. Leigh, and R. J. Beverton. 1991. An examination of a tag-shedding assumption, with application to southern bluefin tuna. ICES J. Mar. Sci. 48:41-51.

Kirkwood, G. P. 1981. Generalised models for the estimation of rates of tag shedding by southern bluefin tuna Thunnus maccoyii. J. Cons. Int. Explor. Mer 39:256-260.

Kirkwood, G. P., and M. H. Walker. 1984. A new method for estimating tag shedding rates, with application to data for Australian salmon, Arripis trutta esper Whitely. Aust. J. Mar. Freshwater Res. 35:601-606.

Otter Research Ltd. 1996. An introduction to AD Model Builder for use in non-linear modeling and statistics. Otter Research Ltd., Nanaimo, B.C., Canada, 68 p.

Wetherall, J. A. 1982. Analysis of double tagging experiments. Fish. Bull. 80(4):687-701.

Xiao, Y. 1996. A general model for estimating tag-specific shedding rates and tag interactions from exact or pooled times at liberty for a double tagging experiment. Can. J. Fish. Aquat. Sci. 53(8):1852-1861.

Xiao, Y., L. P. Brown, T. I. Walker, and A. E. Punt. 1999. Estimation of instantaneous rates of tag shedding for school shark, Galeorhinus galeus and gummy shark, Mustelus antarcticus, by conditional likelihood. Fish. Bull. 97(1):170-184.

Yesaki, M., and A. Waheed. 1992. Results of the tagging program conducted in the Maldives during 1990. IPTP (Indo Pacific Tuna Development and Management Program), 23 p.

Manuscript accepted 21 July 2000.

M. Shiham Adam Renewable Resources Assessment Group T. H. Huxley School of Environment, Earth Sciences and Engineering Imperial College of Science Technology and Medicine, RSM Building Prince Consort Road London SW7 2BP, United Kingdom Present address: Pelagic Fisheries Research Program Joint Institute of Marine and Atmospheric Research University of Hawaii at Manoa Honolulu, Hawaii 96822 E-mail address: msadam@soest.hawaii.edu

Geoffrey P. Kirkwood Renewable Resources Assessment Group T. H. Huxley School of Environment, Earth Sciences and Engineering Imperial College of Science, Technology and Medicine RSM Building, Prince Consort Road London SW7 2BP, United Kingdom

COPYRIGHT 2001 National Marine Fisheries Service

COPYRIGHT 2004 Gale Group