INTRODUCTION
Garfish, Belone belone (Linnaeus, 1761), is widely distributed throughout the Mediterranean Basin and its adjacent seas, as well as in the eastern Atlantic. Three subspecies have been described along its area of distribution (Collette & Parin, Reference Collette, Parin, Whitehead, Bauchot, Hureau, Nielsen and Tortonese1986): B. belone belone (Linnaeus, 1761)—restricted to the north-eastern Atlantic; B. belone gracilis Lowe, 1839—distributed in the south of France to the Canary Islands and in the Mediterranean Sea; and B. belone euxini Günther, 1866, which is found in the Black Sea and the Sea of Azov. Garfish is an epipelagic species, which is found mainly in offshore areas. It migrates towards coastal waters during the spawning period, where it is commercially exploited (Zorica et al., Reference Zorica, Sinovčić and Čikeš Keč2011). In the Adriatic Sea, garfish is the main target of seine net (stretched mesh size 10–30 mm), but is also caught as a by-catch species in purse seine fisheries.
There are some published studies on its biology from the areas of its occurrence (Rosenthal & Fonds, Reference Rosenthal and Fonds1973; Dorman, Reference Dorman1989, Reference Dorman1991; Bedoui et al., Reference Bedoui, Garbi and El Abed2002; Fehri-Bedoui & Gharbi, Reference Fehri-Bedoui and Gharbi2004; Uçkun et al., Reference Uçkun, Akalin, Taşkavak and Toğulga2004; Korzelecka-Orkisz et al., Reference Korzelecka-Orkisz, Formicki, Winnicki, Bonisławska, Szulc, Biernaczyk, Tański and Wawrzyniak2005; Samsun et al., Reference Samsun, Samsun, Bilgin and Kalayci2006; Zorica et al., Reference Zorica, Sinovčić and Čikeš Keč2011). In the Adriatic Sea, data concerning garfish biology are restricted to its reproduction traits (Zorica et al., Reference Zorica, Sinovčić and Čikeš Keč2011), while knowledge about age and growth is lacking. Information concerning age and growth is essential for monitoring year-class strength, making stock assessments, documenting population recovery and constructing the demographic parameters fundamental to population models. Despite the fact that garfish is common in the Adriatic Sea, its total landings in Croatia are unknown (Croatian Bureau of Statistics), and its exploitation status is therefore undetermined. By improving understanding of the life history of the garfish population inhabiting the Adriatic Sea, fishery management will be provided with vital statistical and biological parameters that will allow definition of the present status and management strategy for this species.
MATERIALS AND METHODS
Garfish specimens were collected monthly in the eastern part of the middle Adriatic Sea from January 2003 to December 2008 (Figure 1) with the exception of a small number of months when, due either to logistical problems or weather conditions, sampling was not possible. No samples were obtained in June throughout the investigation, even though sampling and fishing effort did not cease.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-19460-mediumThumb-S002531541200149X_fig1g.jpg?pub-status=live)
Fig. 1. Study area with marked sampling station.
Samples of Belone belone were collected at night using seine net with stretched mesh of size 10 mm (the main net) to 34 mm (the net cod end). The seine net, which is approximately 10 to 30 m deep and 100 m long, is generally operated at 200 to 300 m from the coastline, depending on the sea floor configuration. Preferred sea depth for this type of fishing gear is deeper than 20 m. A total of 3393 individuals, 1166 males and 1195 females were analysed in the laboratory, immediately after landing. Total length (TL), standard length (SL) and fork length (FL) (±0.1 mm), as well as total body weight (W) (±0.01 g) of each specimen, were measured. Sex was macroscopically determined. The sagitta otoliths were removed by opening the otic bulla, washed in freshwater, dried, denoted and stored in plastic tubes prior to analysis.
Conversions between length measures were accomplished with linear regression models using the least squares method (Sokal & Rohlf, Reference Sokal and Rohlf1981), while possible significant differences in length–length relationships among sexes were determined by parametric or non-parametric tests (P < 0.05), depending on data normality. The length–weight relationship for both sexes and overall was calculated by applying the exponential regression W = aTLb , where W = total fish weight in g; TL = total length in cm; a = proportionality constant; and b = allometric growth parameter. The hypotheses of isometric growth, and abbreviation between male and female allometric growth parameter, were tested using a Student's t-test (P < 0.05).
Age was determined by reading the otoliths. The entire otolith was immersed in ethanol solution and examined against a black background with reflected light, using a binocular microscope with amplifications between 1.6 × 25.6 and 1.6 × 51. Hyaline and opaque rings were counted in the dorsoventral part of the concave side of the otoliths. The number of hyaline (translucent) rings was recorded over the rostrum area. Otolith reading was performed on two separate occasions, two months apart. Where there was an age discrepancy between two otolith readings, these otoliths were excluded from further analysis and determined as unreadable, due to false rings. The von Bertalanffy growth equation TL = L∞ [1 − e (−K)(t−t) 0] was fitted to length-at-age data, separately for each sex and overall, using non-linear least-squares parameter estimation, where TL = the length at time t, L∞ = the asymptotic length, K = the growth constant, t0 = the ‘age’ at which TL = 0 (Gulland, Reference Gulland1983). The growth performance index phi-prime (φ′; φ′ = lnK + 2lnL∞ ) was calculated using the equation described by Munro & Pauly (Reference Munro and Pauly1983) in order to compare growth parameters obtained in our study with those reported by other authors. The weights of the undamaged and cleaned otoliths (Wo ) were measured on a Mettler analytical balance (nearest 0.01 mg). The weights of left and right otolith (for 83 otolith pairs) were not significantly different (Mann–Whitney U-test, Z = 0.157; P = 0.876); 82% of otolith pairs had a weight within 0.0002 g of each other. Hence, only one otolith from subsequent pairs was weighed and used in further analysis. Independently of sex, the relationships of otolith weight to total length (Wo –TL) and otolith weight to age (W o –t) were analysed using the power function Wo = aTLb and W o = at b , respectively.
Natural mortality rate (M) was defined using the general regression equations below:
-
■ logM = −0.0066 − 0.279log L∞ + 0.6543logK + 0.4634logT (Pauly, Reference Pauly1980), where L∞ and K are the parameters of the von Bertalanffy growth equation and parameter T is the mean value of the water temperature (17.35°C) in the investigated area (Grbec et al., Reference Grbec, Vilibić, Bajić, Morović, Beg Paklar, Matić and Dadić2007);
-
■ M = 1.521/X0,72 − 0,155 (Rikhter & Efanov, Reference Rikhter and Efanov1976), where X refers to age 2+ at which most of the garfish specimens reach maturity, according to Zorica et al. (Reference Zorica, Sinovčić and Čikeš Keč2010);
-
■ M = 0.996/A0.95 (Taylor, Reference Taylor1959), where A 0.95 is the age at length of 95% L∞ . (18.85 years);
-
■ ln(M) = 0.55 − 1.6ln(L) + 1.44ln(L∞ ) + ln(K) (Gislason et al., Reference Gislason, Daan, Rice and Pope2010), where L is average length at age, L∞ is the asymptotic length and K is the growth constant.
An estimate of instantaneous total mortality rate (Z) was obtained using the linearized length-converted catch curve (Pauly, Reference Pauly1984) in the FISAT program (Gayanilo et al., Reference Gayanilo, Sparre and Pauly1994). The natural logarithm of the ratio between the number of fish in each length-class, and the time needed for the fish to grow through the length-class (ln Ni/Dti), was plotted against their corresponding relative age (t); total mortality was estimated from the descending slope b. Total mortality rates were calculated separately for males, females and overall garfish specimens. The fishing mortality was estimated by subtracting an average M from Z. The exploitation ratio (E) was calculated as E = F/Z (Beverton & Holt, Reference Beverton and Holt1957), while the survival rate (S) was estimated by the formula S = e −Z (Ricker, Reference Ricker1975).
RESULTS
Length–frequency distribution
Garfish from the eastern Adriatic ranged from 20.8 to 75.4 cm (mean ± standard deviation (SD): 38.3 ± 7.94 cm) in total length and 12.21 and 639.25 g (mean ± SD: 70.33 ± 68.53 g) in total body weight (Figure 2). Total length range of males was somewhat narrower (27.7 < TL < 62.6 cm; mean ± SD: 37.4 ± 5.37 cm) than that obtained for females (27.2 < TL < 75.4 cm; mean ± SD: 43.6 ± 9.12 cm). Accordingly, measurable range as well as mean total body weight of males (17.41 < W < 392.53 g; mean ± SD: 58.14 ± 39.74 g) was lower than for females (18.03 < W < 639.25 g; mean ± SD:109.77 ± 93.63 g). There was no statistically significant difference detected between length distribution of male and female garfish specimens collected during the whole investigation (Kolmogorov–Smirnov test, P < 0.05), although all specimens with TL > 62.5 cm were exclusively female.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-40807-mediumThumb-S002531541200149X_fig2g.jpg?pub-status=live)
Fig. 2. Size–frequency distribution of garfish specimens collected in the eastern Adriatic Sea, 2003–2008.
Length–length relationship
Conversions between total (TL) and fork (FL) length and total (TL) and standard (SL) length for each sex and overall are given in Table 1. All obtained length–length regressions were statistically significant (P < 0.05) with a high value of the coefficient of correlation (r > 0.99).
Table 1. The length–length relationships between total length (TL) and fork length (FL) and total (TL) and standard length (SL) separately given for males, females and overall garfish specimens, in the Adriatic Sea, January 2003–December 2008.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-61484-mediumThumb-S002531541200149X_tab1.jpg?pub-status=live)
N, number; r 2, coefficient of determination.
Length–weight relationship
Length–weight relationships for males, females and all individuals combined are shown on Figure 3. Positive allometrical growth was observed in all three cases: for each, a statistically significant difference from isometric value (b = 3) was obtained (P < 0.05). Although for both sexes positive allometry was established, statistically significant differences in slopes between males and females was noted (t = 137.128; df = 2359, P < 0.05).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-26703-mediumThumb-S002531541200149X_fig3g.jpg?pub-status=live)
Fig. 3. Length–weight relationships of garfish males (A), females (B) and overall (C).
Age and growth
Otoliths of 468 garfish specimens, with total length ranging from 23.3 cm to 75.4 cm, were used for aging. Out of 468 otoliths pairs 440, were successfully aged. According to the otolith reading, six age-classes were established for each sex (males: 2+ to 7+ ; females: 3+ to 8+), while in pooled material eight age-classes (1+ to 8+) were determined (Table 2). In the observed garfish population, specimens of 3+ to 4+ years old were the most numerous, while juvenile specimens (0+) were generally missing for the six years of our sampling period. Females dominated in the older age-classes (>7+). The von Bertalanffy growth curves (Figure 4) fitted well to the available length at age data. The estimated parameters of the equation were: L∞ = 85.2 cm (TL), K = 0.159, to = −0.3223, r 2 = 0.9317 males; L∞ = 89.5 cm (TL), K = 0.166, t o = −0.063, r 2 = 0.9476 females; L∞ = 90.3 cm (TL), K = 0.158, to = −0.109, r 2 = 0.9500 overall. In our consideration of obtained growth parameters, φ′ was calculated as 3.062, 3.124 and 3.110 for males, females and pooled material, respectively.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-95443-mediumThumb-S002531541200149X_fig4g.jpg?pub-status=live)
Fig. 4. The von Bertalanffy growth curves and its equations for male (A), female (B) and overall (C) garfish specimens, in the Adriatic Sea, 2003–2008.
Table 2. Age-class distribution of male, female and all examined garfish specimens with specified total length range and mean value for each age class, in the Adriatic Sea, January 2003–December 2008.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-21755-mediumThumb-S002531541200149X_tab2.jpg?pub-status=live)
N, number; x , mean TL; SD, standard deviation.
Weights of 208 undamaged garfish otoliths were noted. Garfish otolith weight ranged from 1.26 mg (TL = 28.2 cm; W = 18.9 g) to 18.12 mg (TL = 73.6 cm; W = 639.25 g), while the mean value of otolith weights was 4.66 ± 2.89 mg. The age–otolith weight key is given in Table 3. Within each age-class, considerable variation of otolith weight was noticed, which created some overlaps between the different age-classes. Significant exponential relationships between otolith weight and total length and mean otolith weight and age were established (Figure 5); larger garfish specimens have heavier otoliths than smaller ones of the same age.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-72285-mediumThumb-S002531541200149X_fig5g.jpg?pub-status=live)
Fig. 5. Garfish exponential relationship of: (A) otolith weights (Wo)–total body lengths (TL) and (B) mean otolith weight with its standard deviations (Wo ± SD)–age-classes (t) with associated equations and coefficients of determination.
Table 3. Age–otolith weight key for garfish individuals collected from the Adriatic Sea, January 2003–December 2008.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-71522-mediumThumb-S002531541200149X_tab3.jpg?pub-status=live)
N, number; SD, standard deviation.
Mortality
Natural mortality rates obtained by Rikhter & Efanov (Reference Rikhter and Efanov1976: M♂ = 0.590 yr−1; M♀ = 0.548 yr−1; M = 0.520 yr−1), Gislason et al. (Reference Gislason, Daan, Rice and Pope2010: M♂ = 0.417 yr−1; M♀ = 0.378 yr−1; M = 0.451 yr−1) and Pauly (Reference Pauly1980: M♂ = 0.321 yr−1; M♀ = 0.326 yr−1; M = 0.315 yr−1) seems reasonable for this pelagic fish species, whilst estimates of M derived from Taylor (Reference Taylor1959: M♂ = 0.054 yr−1; M♀ = 0.055 yr−1; M = 0.0053 yr−1) were too low for a species that links low and high trophic levels. So, for further analysis natural mortality rate was define as an average value of values provided by Rikhter & Efanov (Reference Rikhter and Efanov1976), Pauly (Reference Pauly1980) and Gislason et al. (Reference Gislason, Daan, Rice and Pope2010). All mortality rates (M, Z and F) for males, females and pooled material, as well as exploitation ratio (E) and survival rate (S), are given in Table 4. Total mortality rates (Z) were assessed by using the length converted catch method, but without garfish specimens that were not fully recruited (0–2+ years): this is shown in Figure 6.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-54931-mediumThumb-S002531541200149X_fig6g.jpg?pub-status=live)
Fig. 6. Length converted catch curve for males (A), females (B) and all (C) garfish specimens collected from the eastern Adriatic Sea, 2003–2008.
Table 4. Mortality and exploitation rates of garfish specimens given separately for each sex and overall, in the Adriatic Sea, January 2003–December 2008.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-53214-mediumThumb-S002531541200149X_tab4.jpg?pub-status=live)
DISCUSSION
Length–frequency distribution
Over the study period (2003–2008) samples from the garfish population inhabiting the eastern Adriatic were collected in order to study the life history of this species. Data concerning garfish reproduction patterns have been published (Zorica et al., Reference Zorica, Sinovčić and Čikeš Keč2011), while data describing its age and growth are presented here. During the sampling period, the length-range of analysed garfish specimens was considerably wider in comparison with earlier studies.
Length–length relationships
Length–length relationships were linear with r2 > 0.990 (Table 1). Estimated conversions among morphological measurements and proportional increases during the garfish growth were in accordance with those previously reported for the Adriatic Sea (Sinovčić et al., Reference Sinovčić, Franičević, Zorica and Čikeš Keč2004) and Mediterranean (Fehri-Bedoui & Gharbi, Reference Fehri-Bedoui and Gharbi2004).
Length–weight relationship
In general, positive allometry was established for garfish inhabiting the Adriatic Sea. Male garfish specimens gained in weight significantly more than females over their lifespan. This was also supported by the fact that average body weight of males was higher than those noted for females from the same length-classes (53.0 cm < TL < 62.0 cm; W♂ = 238.01 ± 58.721 g; W♀ = 228.86 ± 54.154 g). An overview of the literature for garfish allometrical coefficients (b) given in Table 5 generally corresponds with findings obtained in this paper. Slight allometric coefficient deviations between different geographical areas could be due to sampling strategy, degree of stomach fullness, gonad maturity, sex, size-range, health and general fish condition and preservation techniques (Tesch, Reference Tesch and Ricker1971).
Table 5. Total length (TL) or fork length (FL) range, allometric parameter (b) and von Bertalanffy growth parameters of Belone belone from various geographical regions.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170720053303-94824-mediumThumb-S002531541200149X_tab5.jpg?pub-status=live)
L ∞ , asymptotic length; K, the growth constant; t0 , the ‘age’ at which TL = 0.
Age and growth
Garfish estimation of age, using sagittal otoliths, was successful considering that only 6% of sampled otoliths were excluded from further analysis due to occurrence of false rings, which were difficult to interpret. In analysed material, age-class 3+ was dominant (45.9%). Specimens of age-class 0+ were absent, and a decreasing trend in numbers of older age-classes was noticed. This is probably correlated with selectivity of fishing gear and the fact that smaller and younger garfish specimens do not accompany their elders (Dorman, Reference Dorman1991; Zorica et al., Reference Zorica, Sinovčić and Čikeš Keč2010). The von Bertalanffy growth function provided an accurate description of somatic garfish growth at the ages encountered. The lack of small specimens may be a possible limitation of this study. Growth parameters (L∞ , K and t0 ) and growth performance index phi-prime (φ′) from the present study and previously published studies are given in Table 5. There were also some data concerning the age–frequency distribution of the garfish population inhabiting the Atlantic (Dorman, Reference Dorman1989, Reference Dorman1991) but as no length-range data or growth parameters were given, results were not comparable. Nevertheless, the values of growth constant (K) obtained in this study were very close to values estimated for the Mediterranean and Black Sea, which in general identified garfish as a slow-growing fish species. For our study, the highest value was recorded most probably due to wider length- and age-range. Garfish growth performance index (φ′) varied within a narrow limit over its area of distribution. Observed discrepancy of the performance index could be explained by the different sea temperatures between geographical areas, which affect the otolith growth and result in different formation of the annulus (Morales-Nin & Ralston, Reference Morales-Nin and Ralston1990). We also produced growth curves separately by sex. Namely, as females were more abundant in older age- classes they attained somewhat higher values of K and L∞ than males. This sex-related growth pattern was similar to one reported by Uçkun et al. (Reference Uçkun, Akalin, Taşkavak and Toğulga2004); males grow more slowly and reach a lower L∞ value than females.
Recent investigations of otolith weights confirmed them as potential age predictors, as it is known that their weight increases throughout the lifespan of an individual. Throughout this study garfish otolith weights were noted and statistically significant exponential relationships between otolith weight–total body length and otolith weight–age were established (Figure 5). Hence, the previously mentioned assumption—otolith weight as age predictor—was also confirmed also for this species. Only one other study dealing with garfish otolith weight has been published (Zorica et al., Reference Zorica, Sinovčić and Čikeš Keč2011), which reported linear correlation between otolith weight and total garfish length (Wo = 0.25TL–5.85, r2 = 0.904). Further differences might be observed if different sample size, in particular age/length distribution, is used.
Mortality
With a basic and reliable knowledge about age and growth characteristics, fundamental practice for good management can be established. Here, obtained estimates of natural mortality, mortality from fishing, and total mortality of garfish (Table 4) can only be compared with the values previously given by Samsun (Reference Samsun1996: M = 0.38; E = 0.66) and Samsun et al. (Reference Samsun, Samsun, Bilgin and Kalayci2006: M = 0.23; Z = 1.24; E = 0.81) on garfish populations inhabiting the Black Sea. Mortality values in these two mentioned studies pointed out that garfish in the Adriatic had a higher natural mortality value, but lower values for fishing mortality and total mortality than garfish in the Black Sea. Observed discrepancy between natural mortality rates is probably due to use of indirect methods that unfortunately generates some bias in natural mortality estimates as they employ different von Bertalanffy growth parameters and the age at which garfish are considered to reach maturity. In addition, different fishing gears are used in different geographical regions: fishing pressure is therefore unequal, resulting in dissimilar mortality rates. Fishing mortality was slightly higher than natural mortality (F > M) and the estimate of the exploitation rate was almost optimal (E = 0.513). According to Gulland (Reference Gulland1971), fish stock is optimally exploited when its exploitation ratio is lower than 0.5. Garfish populations in the Adriatic can therefore be considered endangered.
ACKNOWLEDGEMENTS
This study was supported by the Ministry of Science, Education and Sports of the Republic of Croatia, as a part of the research project ‘Biodiversity and Management of Pelagic and Demersal Resources of the Adriatic Sea (001-0013077-0532)'. The referees’ comments and suggestions are greatly appreciated.