INTRODUCTION
⌅The fishes of the family Triglidae, commonly known as gurnards or sea robins, are medium-sized demersal fishes that are found on soft bottoms in coastal areas up to 500 m depth in all oceans, in both tropical and temperate regions ( Hureau 1986 Hureau J-C. 1986. Triglidae. In: Whitehead P.J.P., Bauchot M.L., et al. (eds), Fishes of the North-Eastern Atlantic and the Mediterranean. Unesco, pp. 1230-1238. https://doi.org/10.2307/1444931 , Richards and Jones 2002 Richards W., Jones D. 2002. Preliminary classification of the gurnards (Triglidae: Scorpaeniformes). Mar. Fresh. Res. 53: 274-282. https://doi.org/10.1071/MF01128 ). Gurnards are commercially important in European waters. In France, they constitute the most important fishery in the western English Channel (ICES Division VIIe), where the market is well established ( ICES 2006 ICES. 2006. Report of the working group on the assessment of new MOU species (WGNEW). ICES (CM Papers and Reports), CM2006/ACFM:11, 240 pp.); in the North Sea, the grey gurnard, Eutrigla gurnardus (Linnaeus, 1758) is used for human consumption ( ICES 2014 ICES 2014. Report of the Working Group on Assessment of New MoU Species (WGNEW). ICES (CM Papers and Reports), CM2014/ACOM:21, 162 pp.); in UK coastal waters especially in the English Channel, gurnards are targeted commercially ( Marriott et al. 2010 Marriott A.L., Latchford J.W., McCarthy I.D. 2010. Population biology of the red gurnard (Aspitrigla cuculus L.; Triglidae) in the inshore waters of Eastern Anglesey and Northwest Wales. J. Appl. Ichth. 26: 504-512. https://doi.org/10.1111/j.1439-0426.2010.01455.x ); and in the Mediterranean, tub gurnard, Chelidonichthys lucerna (Linnaeus, 1758), streaked gurnard, Chelidonichthys lastoviza (Bonnaterre, 1788) and longfin gurnard, Chelidonichthys obscurus (Walbaum, 1792) are a valuable by-catch of mixed demersal fisheries ( Colloca et al. 2019 Colloca F., Milisenda G., Capezzuto F., et al. 2019. Spatial and temporal trend in the abundance and distribution of gurnards (Pisces: Triglidae) in the northern Mediterranean Sea. Sci. Mar. 83: 101-116. https://doi.org/10.3989/scimar.04856.30A ). When landed, they are mainly reported as “mixed gurnards” making it difficult to interpret official landings, but an average of 3905 t of gurnards were landed yearly between 2010 and 2017 in European Union-28 ( Eurostat 2021 Eurostat. 2021. Catches - major fishing areas (from 2000 onwards). https://ec.europa.eu/eurostat/databrowser/view/FISH_CA_MAIN__custom_560273/bookmark/table?lang=en&bookmarkId=40fc2884-8708-401b-9700-c9977a5bfca6/ Accessed online 11 February 2021.). In 2010, legislation in the EU changed, requiring gurnards to be landed under their exact species names ( ICES 2014 ICES 2014. Report of the Working Group on Assessment of New MoU Species (WGNEW). ICES (CM Papers and Reports), CM2014/ACOM:21, 162 pp.). Except for red gurnard, Chelidonichthys cuculus (Linnaeus, 1758), for which the precautionary approach was applied and landings were recommended not to exceed 2894 t in 2020 and 2021 in ICES subareas 3-8 ( ICES 2019 ICES. 2019. ICES Fisheries Overviews: Bay of Biscay and Iberian Coast ecoregion - Fisheries overview, including mixed-fisheries considerations. ICES (CM Papers and Reports), ICES Advice 2019:46, 35 pp.), there is currently no technical measure specific to gurnard species in EU waters, and their exploitation is only subject to the general regulation in the areas where they are caught ( ICES 2014 ICES 2014. Report of the Working Group on Assessment of New MoU Species (WGNEW). ICES (CM Papers and Reports), CM2014/ACOM:21, 162 pp.). Nevertheless, the increasing interest in gurnards as potential commercial species has promoted the development of monitoring programmes for stock assessment, including the one recommended by ICES for C. cuculus within the new Memorandum of Understanding between the European Commission and ICES, aimed at obtaining information on biological parameters ( ICES 2006 ICES. 2006. Report of the working group on the assessment of new MOU species (WGNEW). ICES (CM Papers and Reports), CM2006/ACFM:11, 240 pp.). However, the impact of fishing on gurnard populations has rarely been examined, and studies on key biological parameters and population dynamics used for fisheries management have focused on targeted species, ignoring by-catch or discarded ones, resulting in a lack of biological data and incomplete information about fishing exploitation ( Ordines et al. 2014 Ordines F., Farriols M.T., Lleonart J., et al. 2014. Biology and population dynamics of by-catch fish species of the bottom trawl fishery in the western Mediterranean. Mar. Sci. 15: 613-325. https://doi.org/10.12681/mms.812 ).
This is specifically the case for gurnards off the Portuguese coast, where eight species are caught by the artisanal (mainly polyvalent) and trawler (mainly bottom) fleets: C. cuculus, C. lastoviza, C. lucerna, C. obscurus, E. gurnardus, large-scaled gurnard, Lepidotrigla cavillone (Lacepède, 1801), spiny gurnard, Lepidotrigla dieuzeidei Blanc and Hureau 1973, and piper gurnard, Trigla lyra Linnaeus, 1758. Gurnards are within the top 25 fish species most landed in Portuguese waters, with an average of 371 t landed per year between 2010 and 2019 ( INE and DRGM 2011-2020 INE, DGRM. 2011-2020. Estatísticas da Pesca 2011 to 2020. https://www.ine.pt/xportal/xmain?xpid=INE&xpgid=ine_publicacoes&PUBLICACOEStipo=ea&PUBLICACOEScoleccao=107656&selTab=tab0&xlang=pt/ Accessed online 11 February 2022.). Among the four most landed species are C. cuculus and T. lyra ( Rocha et al. 2018 Rocha A., Feijó D., Gonçalves P. 2018. Gurnards: species landings’ composition in ICES Division 27.9a. ICES WGWIDE Report 2018. https://doi.org/10.13140/RG.2.2.24277.47841 / Accessed online 11 February 2022.). For C. cuculus, information available on life-history parameters (age, growth and reproduction), abundance and distribution, length and weight relationships, fisheries data, stock identification and management is scattered in European waters ( ICES 2006 ICES. 2006. Report of the working group on the assessment of new MOU species (WGNEW). ICES (CM Papers and Reports), CM2006/ACFM:11, 240 pp., Marriott et al. 2010 Marriott A.L., Latchford J.W., McCarthy I.D. 2010. Population biology of the red gurnard (Aspitrigla cuculus L.; Triglidae) in the inshore waters of Eastern Anglesey and Northwest Wales. J. Appl. Ichth. 26: 504-512. https://doi.org/10.1111/j.1439-0426.2010.01455.x , Ordines et al. 2014 Ordines F., Farriols M.T., Lleonart J., et al. 2014. Biology and population dynamics of by-catch fish species of the bottom trawl fishery in the western Mediterranean. Mar. Sci. 15: 613-325. https://doi.org/10.12681/mms.812 ). For T. lyra, studies on age and growth, length-age composition, reproduction and mortality were published for the Saronikos Gulf ( Papaconstantinou 1981 Papaconstantinou C. 1981. Age and growth of piper, Trigla lyra, in Saronikos Gulf (Greece). Cybium 5: 73-87.) and the Aegean Sea ( Papaconstantinou et al. 1992 Papaconstantinou C., Petrakis G., Caragitsou E., Mytilineou C. 1992. Preliminary study on the biology of piper (Trigla lyra, L. 1758) in the Aegean Sea. FAO Fisheries Report 477: 127-137., Ismen et al. 2016 Ismen A., Ihsanoglu M.A., Ismen P., Yigin C.C. 2016. Age, growth and reproduction of piper gurnard, Trigla lyra in Saros Bay, the northern Aegean Sea. Rapp. Comm. int. Mer Médit. 41: 340.).
Key biological parameters routinely used in fisheries to estimate exploitation level ( ICES 2018 ICES. 2018. Working Group on Biological Parameters (WGBIOP). ICES (CM Papers and Reports), CM2018/EOSG:07, 186 pp.) and still absent from gurnards were analysed in the present study. Therefore, the main goals were (1) to estimate the sex ratio and the weight-length relations, (2) to validate the assignment of age, (3) to model growth, and (4) to estimate mortality and obtain an initial approximation of the exploitation level for C. cuculus and T. lyra in the coastal waters of western Portugal.
MATERIALS AND METHODS
⌅Sampling
⌅A total of 558 T. lyra specimens and 425 C. cuculus specimens were collected between July 2018 and December 2019 from commercial fishing vessels operating off mainland Portugal (Peniche, west coast, 39°03′20.4′′N, 9°39′54.2′′W to 38°00′32.7′′N, 8°55′30.7′′W; Fig. 1 ) using trawls, trammel nets, traps and pots. 21 T. lyra individuals under 22 cm total length (11.2 cm to 15.3 cm) that are usually rejected on board were directly obtained from fishermen. No C. cuculus individuals were collected in June. Total length (TL, to the nearest 0.1 cm) and eviscerated weight (EW, to the nearest 0.01 g) were recorded, and the sagittal otoliths were removed, rinsed with water, air-dried and stored in tagged vials. Sex was recorded by macroscopic observation of the gonads. For mortality estimates, the length composition of landings was recorded monthly for each species at Peniche during the study period.
Sex ratio and length-weight relationship
⌅Significant departures from the 1:1 sex ratio were investigated using chi-square tests ( Legendre and Legendre 1998 Legendre P., Legendre L. 1998. Numerical ecology. Elsevier, Amsterdam, 989 pp.). The Mann-Whitney test (U) was applied to analyse variation in TL and EW between sexes. The relationship between TL (cm) and EW (g) was calculated using a power function:
where a is the intercept and b is the allometric exponent. Growth allometry between sexes was investigated using the Student t-test ( Zar 1999 Zar J.H. 1999. Biostatistical analysis. Prentice Hall, NY, 663 pp.).
Ageing methodology and validation
⌅Using a digital camera (DFC 290; Leica, Wetzlar, Germany) coupled to a stereomicroscope (SMZ 745T; Nikon, Kanagawa, Japan), digital images were acquired from whole right sagittal otoliths previously immersed in water with the medial (convex) side down under reflected light and a dark background at 18× magnification. Otolith increment widths were measured along the growth axis from the nucleus to the posterior edge using the “Marking growth rings in fish otoliths” project ( https://sils.fnwi.uva.nl/bcb/objectj/examples/otoliths/MD/otoliths.html ), for ObjectJ plugin ( https://sils.fnwi.uva.nl/bcb/objectj/ ) of the ImageJ freeware image processing software (National Institute of Mental Health, Bethesda, MD, USA). The posterior edge was chosen because this region has a longer radius and therefore offers better differentiation and individualization of the growth increments ( Fig. 2 ).
The annual periodicity of annulus formation (annual growth increments) was validated semi-directly ( Panfili and Morales-Nin 2002 Panfili J., Morales-Nin B. 2002. Growth and growth analysis. In: Panfili J., Pontual H., Troadec H., Wright P. (eds), Manual of Fish Sclerochronology. Ifremer-IRD coedition, pp. 146-166.) through the analysis of the type of edge and the marginal increment ratio (MIR) ( Samamé 1977 Samamé M. 1977. Determinación de la edad y crecimiento de la sardina Sardinops sagax sagax (J). Boletín-Instituto del Mar del Perú 3: 95-112.). The edge of the otoliths was classified as opaque (highly calcified) or translucent (less calcified). The edge type and the mean MIR and standard deviation (SD) were plotted per month to assess annual trends in the formation of growth annuli.
Precision and bias in otolith reading
⌅To evaluate the precision of age assignments and to ensure consistency, representative subsamples of sagitta otolith (75 T. lyra and 148 C. cuculus) covering all length classes and both sexes were read by two authors. Whenever the age estimates did not coincide, a third reading was performed. The average per cent error (APE) ( Beamish and Fournier 1981 Beamish R.J., Fournier D.A. 1981. A method for comparing the precision of a set of age determinations. Can. J. Fish. Aquat. Sci. 38: 982-983. https://doi.org/10.1139/f81-132 ), the coefficient of variation (CV) and the index of precision (D) ( Chang 1982 Chang W.Y.B. 1982. A statistical method for evaluating the reproducibility of age determination. Can. J. Fish. Aquat. Sci. 39: 1208-1210. https://doi.org/10.1139/f82-158 ) were used to compare age assignments between readers. Bias was evaluated on the basis of age-bias plots ( Campana et al. 1995 Campana S.E., Annand M.C., McMillan J.I. 1995. Graphical and statistical methods for determining the consistency of age determinations. T. Am. Fish. Soc. 124:131-138. https://doi.org/10.1577/1548-8659(1995)124<0131:GASMFD>2.3.CO;2 ), which help identify departures from the 1:1 equivalence line, and systematic differences in assigned ages between readers were assessed with the Bowker-type test for symmetry ( Hoenig et al. 1995) Hoenig J.M., Morgan M.J., Brown C.A. 1995. Analysing differences between two age determination methods by tests of symmetry. Can. J. Fish. Aquat. Sci. 52: 364-68. https://doi.org/10.1139/f95-038 using the FSA package of R (cran.r-project.org) ( Ogle 2020 Ogle D.H. 2020. FSA: Simple Fisheries Stock Assessment Methods. R package version 0.8.32.).
As there was a good agreement between readers (CV<7.6% and APE<5.5%; Campana 2001 Campana S.E. 2001. Accuracy, precision and quality control in age determination, including a review of the use and abuse of age validation methods. J. Fish Biol. 59: 197-242. https://doi.org/10.1111/j.1095-8649.2001.tb00127.x ), the remaining otoliths were read twice by the first author and readings were used to assign a modal age assuming the 1 January as the birthdate.
Growth model
⌅Because there was a low representation of small individuals (<22 cm) in the sample of T. lyra, the von Bertalanffy growth parameters ( von Bertalanffy 1938 von Bertalanffy L. 1938. A quantitative theory of organic growth (inquiries of growth laws II). Hum. Biol. 10: 181-213.) for this species were estimated using (1) whole-otolith reading, (2) back-calculation ( Francis 1990 Francis R.I.C.C. 1990. Back-calculation of fish length: a critical review. J. Fish Biol. 36: 883-902. https://doi.org/10.1111/j.1095-8649.1990.tb05636.x , Jones 2000 Jones C.M. 2000. Fitting growth curves to retrospective size-at-age data. Fish. Res. 46: 123-129. https://doi.org/10.1016/S0165-7836(00)00139-9 , Wilson et al. 2009 Wilson J.A., Vigliola L., Meekan M.G. 2009. The back-calculation of size and growth from otoliths: validation and comparison of models at an individual level. Exp. Mar. Biol. Ecol. 368: 9-21. https://doi.org/10.1016/j.jembe.2008.09.005 ), and (3) a combination of the previous methods ( Gordo et al. 2016 Gordo L.S., Neves A., Vieira A.R., Paiva R.B., Sequeira V. 2016. Age, growth and mortality of the comber Serranus cabrilla (Linnaeus, 1758) in the Eastern Atlantic. Mar. Biol. Res. 12: 656-662. https://doi.org/10.1080/17451000.2016.1169295 ).
In (1) age was directly estimated for individuals between 10.1 and 44.4 cm. In (2) the mean length-at-age was estimated using a power function between otolith radius (OR) and TL:TL=a OR b , where a is the intercept and b is the slope ( Folkvord and Monsegaard 2002 Folkvord A., Mosegaard H. 2002. Growth and growth analysis. In: Panfili J., Pontual H., Troadec H., Wright P. (eds), Manual of Fish Sclerochronology. Ifremer-IRD. pp. 146-166.).
In (3) the two previous approaches were used, with back-calculation providing information on length-at-ages 1 and 2 years and whole-otolith readings data for ages 3 to 11 years. The three models were judged using the Akaike information criteria (AIC) ( Akaike 1973 Akaike H. 1973. Information theory and an extension of the maximum likelihood principle. In: Petrov B.N., Csaki F. (eds), Second International Symposium on Information Theory. Akademiai Kiado, pp. 267-281., Burnham and Anderson 2002 Burnham K.P. and Anderson D.R. 2002. Model Selection and Inference: A Practical Information-Theoretic Approach. Springer-Verlag, New York, 355 pp.). The relative fit of each model was assessed via ΔAIC, with models for which ΔAIC<2 having substantial support and models for which ΔAIC>10 having essentially no support. The model showing the lowest AIC value was chosen as the best approach to estimate the von Bertalanffy growth parameters.
To estimate C. cuculus growth, sex-specific length-at-age data based on whole-otolith readings were fitted to the von Bertalanffy growth model ( von Bertalanffy 1938 von Bertalanffy L. 1938. A quantitative theory of organic growth (inquiries of growth laws II). Hum. Biol. 10: 181-213.). The likelihood ratio tests ( Kimura 1980 Kimura D.K. 1980. Likelihood methods for the von Bertalanffy growth curve. Fish. Bull. 77: 765-76.) were used to evaluate the significance of differences in growth parameters between sexes in both species. The growth performance index (ɸ’=log10 K+2 log10 L ∞; Munro and Pauly 1983 Munro J.L., Pauly D. 1983. A simple method for comparing growth of fishes and invertebrates. ICLARM Fishbyte 1: 5-6.) was used to compare the growth of fish between this and other studies previously published, and the Student t-test was used to compare mean ɸ’ for the Atlantic and Mediterranean C. cuculus populations. Insufficient data precluded this comparison for T. lyra.
Analyses were performed in R software ( RStudio Team 2020 RStudio Team 2020. RStudio: Integrated development environment for R. Boston, MA: PBC.) using the FSA ( Ogle 2020 Ogle D.H. 2020. FSA: Simple Fisheries Stock Assessment Methods. R package version 0.8.32.) and Rcpp ( Eddelbuettel et al. 2020 Eddelbuettel D., Francois R., Allaire J.J., et al. 2020. Rcpp: Seamless R and C++ Integration. R package version 1.0.4.6.) packages for growth models and the fishmethods package ( Nelson 2019 Nelson G.A. 2019. fishmethods: Fishery Science Methods and Models. R package version 1.11-1.) for the likelihood ratio tests.
Mortality
⌅Estimates of mortality for both species were obtained considering only fish caught by trammel nets, the main fishing gear used to capture T. lyra and C. cuculus in Portuguese waters. The instantaneous total mortality rate (Z, year-1) was estimated based on the age-converted catch curve using the Chapman-Robson mortality estimator ( Chapman and Robson 1960 Chapman D.G., Robson D.S. 1960. The analysis of a catch curve. Biometrics 16: 354-368. https://doi.org/10.2307/2527687 ) and the first age group as 1 year older than the age of peak ( Smith et al. 2012 Smith M.W., Then A.Y., Wor C., et al. 2012. Recommendations for Catch-Curve Analysis. N. Am. J. Fish. Manange. 32: 956-967. https://doi.org/10.1080/02755947.2012.711270 ). The instantaneous natural mortality rate (M, year-1) was estimated using Paulynls-T estimator ( Pauly 1980 Pauly D. 1980. On the interrelationship between natural mortality, growth parameters and mean environmental temperature in 175 fish stocks. ICES J. Mar. Sci. 39: 195-212. https://doi.org/10.1093/icesjms/39.2.175 ), M=4.118 K 0.73 L ∞ -0.33 and one-parameter t max, M est=5.109/t max ( Then et al. 2015 Then A.Y., Hoenig J.M., Hall N.G., Hewitt D.A. 2015. Evaluating the predictive performance of empirical estimators of natural mortality rate using information on over 200 fish species. ICES J. Mar. Sci. 72: 82-92. https://doi.org/10.1093/icesjms/fsu136 ). Fishing mortality (F, year-1) was estimated as F=Z-M and the exploitation rate as E=F/Z. Mortalities were estimated using the FSA package of R ( Ogle 2020 Ogle D.H. 2020. FSA: Simple Fisheries Stock Assessment Methods. R package version 0.8.32.).
All results were considered statistically significant at the 5% level.
RESULTS
⌅Length-weight relationships
⌅A total of 420 females and 138 males of T. lyra and 243 females and 182 males of C. cuculus were sampled ( Table 1 ). The overall sex ratio (M:F) differed significantly from 1:1 for both species (T. lyra, 0.33:1; χ 2=142.52, df=1, P<0.001; C. cuculus, 0.75:1; χ 2=8.76, df=1, P<0.05), with females prevailing. Significant differences for TL and EW between sexes and for both species were found (Mann-Whitney test, T. lyra TL, U=39905, P<0.001; T. lyra EW, U=40528, P<0.001; C. cuculus TL, U=30941, P<0.001; C. cuculus EW, U=30316, P<0.001), with females being larger and heavier.
| T. lyra | C. cuculus | |||
|---|---|---|---|---|
| Females | Males | Females | Males | |
| N | 420 | 138 | 243 | 182 |
| TL range (cm) | 11.2-44.4 | 10.9-41.0 | 10.2-42.1 | 13.8-33.4 |
| EW range (g) | 15.09-756.92 | 14.63-679.22 | 7.83-679.22 | 19.32-326.09 |
Both species showed a strong relationship between TL and EW expressed by the following equations: T. lyra females EW=0.006 TL3.0756 (r=0.98, P<0.056, n=420) and males EW=0.006 TL3.0648 (r=0.98, P<0.05, n=138); C. cuculus females EW=0.006 TL3.1252 (r=0.95, P<0.05, n=243) and EW=0.008 TL3.0107 (r=0.95, P<0.05, n=182). For both species, there were no significant differences in growth between females and males (paired t-test, T. lyra t-test=-0.207, P>0.05; C. cuculus t-test=-0.026, P>0.05), and b was not significantly different from 3 in both sexes (paired t-test, T. lyra females t-test=2.645, P>0.05; T. lyra males t-test=1.190, P>0.05; C. cuculus females t-test=0.04, P>0.05; C. cuculus males t-test=0.004, P>0.05), indicating that growth was isometric.
Ageing method validation
⌅A regular pattern of deposition of annuli was visible in the otoliths of both species, with alternate opaque and translucent concentric increments depositing around an opaque nucleus. Increments were wider close to the nucleus and narrower towards the edge of the otoliths ( Fig. 2 ). False increments were more frequent in T. lyra, hindering interpretation.
The deposition patterns of the growth increments showed that the first and second increments were well separated ( Fig. 3 ), particularly in T. lyra. From age 3 onwards, separation between increments decreased with age, being more evident for T. lyra, which attained greater ages than C. cuculus. The distances (mean±SD) from the nucleus to the first ten (T. lyra) and six (C. cuculus) increments were recorded (supplementary Table S1 ).
For the analysis of the proportion of opaque and translucent edges and MIRs, 385 otoliths of T. lyra (n=385) and 422 of C. cuculus (n=422) were used, suggesting similar annual patterns of growth increment formation ( Fig. 4 ): opaque zones occurred more frequently in spring and summer, with MIR tending to increase from May to September, especially for T. lyra; translucent edges dominated in autumn and winter, with the lowest values of MIR between October and April, indicating the formation of new increments during this period. Therefore, it is possible to assume that an annual growth increment corresponds to consecutive opaque and translucent growth zones, and age was chosen to be assigned by counting the translucent zones.
Precision and bias in age readings
⌅The indices of precision APE, CV, and D in age assignments between readers for both species are presented in Table 2 , and the age-bias plots in Figure S1. There was good precision in age readings, though better results were achieved for C. cuculus than for T. lyra, with a total agreement of 87% in the former and of 60% in the latter. No evidence of systematic disagreement between readers was obtained for either species (T. lyra, χ2 R1vsR2=19, df=5, P>0.05; C. cuculus, χ2 R1vsR2=13.10, df=7, P>0.05) indicating simple random error (supplementary Fig. S1 ).
| Indices of precision | Trigla lyra | Chelidonichthys cuculus |
|---|---|---|
| APE (%) | 3.943 | 1.819 |
| CV (%) | 5.576 | 2.573 |
| D (%) | 0.200 | 0.064 |
| Total agreement (%) | 60 | 87 |
| Agreement (0±1) (%) | 100 | 100 |
Growth models
⌅In total 558 T. lyra and 425 C. cuculus otoliths were used for ageing. In T. lyra age ranged from 1 to 11 years in females and 1 to 10 years in males, and in C. cuculus from 0 to 9 years in females and 0 to 5 years in males. The age-length keys obtained from direct otolith reading for both species are shown in supplementary Tables S2 and S3 , respectively.
Results for T. lyra from the three methods used to estimate the von Bertalanffy growth parameters are shown in Table 3 . No significant differences in growth parameters between females and males were found using whole-otolith readings (likelihood ratio test: χ2=1.09, df=3, P>0.05), and sexes were pooled for subsequent analysis. There was a strong relationship between TL and OR (TL=101.63 OR1.1457, r=0.855, P<0.05, n=455). Figure 5A presents the von Bertalanffy growth curves obtained from the three methodologies. According to the AIC, the combined (AIC=40.9) and back-calculation (AIC=43.3) models had substantial support as ΔAIC<2, but as the combined method showed the lowest value it was chosen as the best approach for estimating the growth of this species.
| Species | Area | Reference | Reading method | Ageing method | Sex | N | Size [cm] | Age [yr] | L∞ (cm) | K (yr-1) | t0 (yr) | ɸ’ |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| T. lyra | Portuguese coast (A) | This study | WO | DE | All | 558 | 10.9-44.4 | 1-11 | 53.7 | 0.13 | -1.918 | 2.56 |
| T. lyra | Portuguese coast (A) | This study | WO | BC | All | - | - | - | 45.9 | 0.20 | -1.134 | 2.63 |
| T. lyra | Portuguese coast (A) | This study | WO | CM | All | - | - | - | 44.7 | 0.16 | -2.781 | 2.50 |
| T. lyra | Aegean Sea (M) | Papaconstantinou et al., 1992 Papaconstantinou C., Petrakis G., Caragitsou E., Mytilineou C. 1992. Preliminary study on the biology of piper (Trigla lyra, L. 1758) in the Aegean Sea. FAO Fisheries Report 477: 127-137. | BW | BC | All | 2180 | 4.3-47.2 | 0-8 | 57.4 | 0.02 | -1.245 | 1.75 |
| T. lyra | Saronikos Gulf (M) | Papaconstantinou, 1981 Papaconstantinou C. 1981. Age and growth of piper, Trigla lyra, in Saronikos Gulf (Greece). Cybium 5: 73-87. | BW | BC | All | 279 | 4-46 | 0-7 | 74.0 | 0.11 | -0.811 | 2.78 |
| C. cuculus | Portuguese coast (A) | This study | WO | DE | All | 425 | 10.2-42.1 | 0-9 | 40.2 | 0.28 | -1.074 | 2.66 |
| C. cuculus | Western Mediterranean (M) | Ordines et al., 2014 Ordines F., Farriols M.T., Lleonart J., et al. 2014. Biology and population dynamics of by-catch fish species of the bottom trawl fishery in the western Mediterranean. Mar. Sci. 15: 613-325. https://doi.org/10.12681/mms.812 | WO | DE | All | 873 | 8-28 | 0-5 | 37.6 | 0.18 | -1.825 | 2.40 |
| C. cuculus | Northwest Wales and eastern Anglesey, UK (A) | Marriott et al., 2010 Marriott A.L., Latchford J.W., McCarthy I.D. 2010. Population biology of the red gurnard (Aspitrigla cuculus L.; Triglidae) in the inshore waters of Eastern Anglesey and Northwest Wales. J. Appl. Ichth. 26: 504-512. https://doi.org/10.1111/j.1439-0426.2010.01455.x | WO | DE | M | 188 | 15.4-35.0 | 1-6 | 41.3 | 0.21 | -1.250 | 2.55 |
| C. cuculus | Northwest Wales and eastern Anglesey, UK (A) | Marriott et al., 2010 Marriott A.L., Latchford J.W., McCarthy I.D. 2010. Population biology of the red gurnard (Aspitrigla cuculus L.; Triglidae) in the inshore waters of Eastern Anglesey and Northwest Wales. J. Appl. Ichth. 26: 504-512. https://doi.org/10.1111/j.1439-0426.2010.01455.x | WO | DE | F | 552 | 10-5-43.1 | 1-7 | 40.9 | 0.24 | -1.140 | 2.60 |
| C. cuculus | Northwest Wales and eastern Anglesey, UK (A) | Marriott et al., 2010 Marriott A.L., Latchford J.W., McCarthy I.D. 2010. Population biology of the red gurnard (Aspitrigla cuculus L.; Triglidae) in the inshore waters of Eastern Anglesey and Northwest Wales. J. Appl. Ichth. 26: 504-512. https://doi.org/10.1111/j.1439-0426.2010.01455.x | WO | DE | All | 740 | 10.5-43.1 | 1-7 | 42.4 | 0.21 | -1.210 | 2.58 |
| C. cuculus | English Channel (A) | Dorel (1986) Dorel D. 1986. Poissons de L’ Atlantique nord est relations taillepoids. Institut Français de Recherche pour L’exploitation de la Mer, France, 183 pp., ICES (2006) ICES. 2006. Report of the working group on the assessment of new MOU species (WGNEW). ICES (CM Papers and Reports), CM2006/ACFM:11, 240 pp. | - | DE | M | - | - | 1-4 | 35.6 | 0.23 | -3.370 | 2.46 |
| C. cuculus | English Channel (A) | Dorel (1986) Dorel D. 1986. Poissons de L’ Atlantique nord est relations taillepoids. Institut Français de Recherche pour L’exploitation de la Mer, France, 183 pp. ICES (2006) ICES. 2006. Report of the working group on the assessment of new MOU species (WGNEW). ICES (CM Papers and Reports), CM2006/ACFM:11, 240 pp. | - | DE | F | - | - | 1-4 | 41.1 | 0.24 | -2.570 | 2.61 |
| C. cuculus | English Channel (A) | Dorel (1986) Dorel D. 1986. Poissons de L’ Atlantique nord est relations taillepoids. Institut Français de Recherche pour L’exploitation de la Mer, France, 183 pp. ICES (2006) ICES. 2006. Report of the working group on the assessment of new MOU species (WGNEW). ICES (CM Papers and Reports), CM2006/ACFM:11, 240 pp. | - | DE | All | - | 9-43 | 1-4 | 36.3 | 0.24 | -0.170 | 2.50 |
| C. cuculus | Tyrrhenian Sea (M) | Colloca et al., 2003 Colloca F., Cardinale M., Ardizzone G.D. 2003. Tracing the life history of red gurnard (Aspitrigla cuculus) using validated otolith annual rings. J. Appl. Ichth. 19: 1-9. https://doi.org/10.1046/j.1439-0426.2003.00342.x | WO | BC | M | 122 | 7-27 (ALL) | 1-3 | 23.4 | 0.59 | -0.380 | 2.51 |
| C. cuculus | Tyrrhenian Sea (M) | Colloca et al., 2003 Colloca F., Cardinale M., Ardizzone G.D. 2003. Tracing the life history of red gurnard (Aspitrigla cuculus) using validated otolith annual rings. J. Appl. Ichth. 19: 1-9. https://doi.org/10.1046/j.1439-0426.2003.00342.x | WO | BC | F | 127 | 1-3 | 24.2 | 0.74 | -0.070 | 2.64 | |
| C. cuculus | Southwest Adriatic Sea (M) | Maisan et al. (1998) Maisan R., Ungaro N., Marzano M.C., Martino M. 1998. Growth of Aspitrigla cuculus (Osteichthyes, Triglidae) in the south-western Adriatic area: preliminary results. Biol. Mar. Medit. 5: 694-696. | - | DE | M | - | 9.5-34.5 | 1-4 | 29.3 | 0.41 | -0.490 | 2.55 |
| C. cuculus | Southwest Adriatic Sea (M) | Maisan et al. (1998) Maisan R., Ungaro N., Marzano M.C., Martino M. 1998. Growth of Aspitrigla cuculus (Osteichthyes, Triglidae) in the south-western Adriatic area: preliminary results. Biol. Mar. Medit. 5: 694-696. | - | DE | F | - | 9.5-34.5 | 1-5 | 29.2 | 0.40 | -0.680 | 2.53 |
| C. cuculus | Gulf of Lion (M) | Campillo (1992) Campillo A. 1992. Les pêcheries françaises de Méditeranée: synthèse des connaissances. Institut Français de Recherche pour L’exploitation de la Mer, France, 231 pp. | - | DE | All | - | 10-30 | 1-5 | 34.9 | 0.28 | -0.660 | 2.53 |
| C. cuculus | Baie de Douarnenez (A) | Baron (1985) Baron J. 1985. Les Triglides (Teleosteens, Scorpaeniformes) de la Baie de Douarnenez. 1. La Croissance de Eutrigla gurnardus, Trigla lucerna, Trigla lastoviza et Aspitrigla cuculus. Cybium 9: 127-144. | CBU | DE | M | 118 | max 50 | 1-13 | 37.1 | 0.51 | -0.080 | 2.85 |
| C. cuculus | Baie de Douarnenez (A) | Baron (1985) Baron J. 1985. Les Triglides (Teleosteens, Scorpaeniformes) de la Baie de Douarnenez. 1. La Croissance de Eutrigla gurnardus, Trigla lucerna, Trigla lastoviza et Aspitrigla cuculus. Cybium 9: 127-144. | CBU | DE | F | 232 | max 50 | 1-21 | 41.7 | 0.46 | -0.050 | 2.90 |
| C. cuculus | Aegean Sea (M) | Papaconstantinou (1983) Papaconstantinou K.A. 1983. Aspects of the biology of Aspitrigla cuculus (L., 1758) (Pisces, Scorpaeniformes) in the Gulf of Saronikos. Thalassographica 6: 49-75. | - | DE | M | 507 (all) | 4-18 | 1-5 | 20.4 | 0.51 | 0.008 | 2.33 |
| C. cuculus | Aegean Sea (M) | Papaconstantinou (1983) Papaconstantinou K.A. 1983. Aspects of the biology of Aspitrigla cuculus (L., 1758) (Pisces, Scorpaeniformes) in the Gulf of Saronikos. Thalassographica 6: 49-75. | - | DE | F | - | - | 1-7 | 27.6 | 0.22 | 0.008 | 2.22 |
The von Bertalanffy growth parameters estimated for C. cuculus are also presented in Table 3 . No significant differences between males and females were found (likelihood ratio test: χ2=1.36, df=3, P>0.05) and sexes were also pooled. The von Bertalanffy growth curve obtained is presented in Figure 5B . Table 3 also summarizes the von Bertalanffy growth parameters and the growth performance indices obtained for other geographical areas. Mean ɸ’ for the Atlantic and Mediterranean C. cuculus populations was significantly different (Atlantic ɸ’=2.677±0.162; Mediterranean ɸ’=2.463±0.136; t-test=2.690; df=9.754; P<0.05).
Mortality
⌅Estimates of Z were 0.61 and 1.67 year-1 for T. lyra and C. cuculus, respectively. Estimates of M varied with the method used: for T. lyra, between 0.31 (Paulynls-T estimator) and 0.46 year-1 (one-parameter tmax); for C. cuculus between 0.48 (Paulynls-T estimator) and 0.57 year-1 (one-parameter tmax). F varied between 0.14 and 0.30 year-1 for the former, and between 1.11 and 1.19 year -1 for the latter. The exploitation rate (E) was higher for C. cuculus, ranging between 0.66 and 0.71 year-1, while in T. lyra it ranged between 0.23 and 0.49 year-1 ( Table 4 ).
| Trigla lyra | Chelidonichthys cuculus | |||||||
|---|---|---|---|---|---|---|---|---|
| Method | Z (yr-1) | M (yr-1) | F (yr-1) | E (yr-1) | Z (yr-1) | M (yr-1) | F (yr-1) | E (yr-1) |
| 0.61 | 1.67 | |||||||
| Paulynls-T | 0.31 | 0.30 | 0.49 | 0.48 | 1.19 | 0.71 | ||
| One-parameter tmax | 0.46 | 0.14 | 0.23 | 0.57 | 1.11 | 0.66 | ||
DISCUSSION
⌅This study presents the first age, growth and mortality characterization for T. lyra and C. cuculus from the northeast Atlantic (Portuguese) coast, two by-catch species caught by artisanal and trawler fleets. According to ICES (2014) ICES 2014. Report of the Working Group on Assessment of New MoU Species (WGNEW). ICES (CM Papers and Reports), CM2014/ACOM:21, 162 pp., the quality of landings data is poor since the species have usually not been well separated when landed, and updated biological parameters are lacking, particularly for the most important commercial gurnards such as C. cuculus ( ICES 2006 ICES. 2006. Report of the working group on the assessment of new MOU species (WGNEW). ICES (CM Papers and Reports), CM2006/ACFM:11, 240 pp.). For a correct stock assessment and sustainable management of fisheries, detailed biological information from all the different species is needed.
Individuals of both species were obtained from commercial fishing vessels using the same fishing gear, but the length range varied between species. Samples of C. cuculus included individuals of all length classes between 10 and 43 cm, though most were larger than 23 cm, as expected from landings. Nevertheless, the length distribution gathered from the present study falls within those observed for other European areas (see Table 3 ). T. lyra individuals under 22 cm were underrepresented, and no individuals between 16 and 22 cm were collected from landings; the few individuals gathered were obtained directly from fishermen. This could be related to species-specific ecological features, as they have differeent bottom type preferences, different feeding patterns, a specific depth distribution ( Papaconstantinou 1983 Papaconstantinou K.A. 1983. Aspects of the biology of Aspitrigla cuculus (L., 1758) (Pisces, Scorpaeniformes) in the Gulf of Saronikos. Thalassographica 6: 49-75., Lopez-Lopez et al. 2011 Lopez-Lopez L., Preciado I., Velasco F., et al. 2011. Resource partitioning amongst five coexisting species of gurnards (Scorpaeniforme: Triglidae): Role of trophic and habitat segregation. J. Sea Res. 66: 58-68. https://doi.org/10.1016/j.seares.2011.04.012 ) and contrasting reproductive strategies ( Neves et al. 2021 Neves A., Sousa I., Sequeira V., et al. 2021. Enhancing knowledge on low-value fishing species: the distinct reproductive strategy of two gurnard species. J. Fish Biol. 1-12. https://doi.org/10.1111/jfb.14849 ), which may influence their selectivity in the fishery. According to Lopez-Lopez et al. (2011) Lopez-Lopez L., Preciado I., Velasco F., et al. 2011. Resource partitioning amongst five coexisting species of gurnards (Scorpaeniforme: Triglidae): Role of trophic and habitat segregation. J. Sea Res. 66: 58-68. https://doi.org/10.1016/j.seares.2011.04.012 , on the Galician and Cantabrian Sea continental shelves, mature T. lyra over 26 cm are restricted to shallower depths (around 150 m), which is suggested to be related to an ontogenetic shift associated with a change in diet. The fact that individuals under 26 cm are underrepresented in the present study can be related to their higher depth distribution, the size less commercially valued, and the selectivity of the gears regarding the size and the depth where they are used on the Portuguese coast (artisanal gears are used in shallow waters).
Females proved to be more frequent in landings than males for both species in the sample, contrary to findings in the Saronikos Gulf for T. lyra ( Papaconstantinou 1981 Papaconstantinou C. 1981. Age and growth of piper, Trigla lyra, in Saronikos Gulf (Greece). Cybium 5: 73-87.) and in the Aegean Sea for C. cuculus ( Papaconstantinou 1983 Papaconstantinou K.A. 1983. Aspects of the biology of Aspitrigla cuculus (L., 1758) (Pisces, Scorpaeniformes) in the Gulf of Saronikos. Thalassographica 6: 49-75.), where males dominated. Differences in length and weight between sexes were also found, with females being larger and heavier. This is consistent with records for C. cuculus in northwest Wales and eastern Anglesey ( Marriott et al. 2010 Marriott A.L., Latchford J.W., McCarthy I.D. 2010. Population biology of the red gurnard (Aspitrigla cuculus L.; Triglidae) in the inshore waters of Eastern Anglesey and Northwest Wales. J. Appl. Ichth. 26: 504-512. https://doi.org/10.1111/j.1439-0426.2010.01455.x ), although in the Saronikos Gulf no differences in length between sexes were observed for T. lyra ( Papaconstantinou 1981 Papaconstantinou C. 1981. Age and growth of piper, Trigla lyra, in Saronikos Gulf (Greece). Cybium 5: 73-87.).
There were no significant differences in length-weight relationships between males and females, and growth proved to be isometric in both studied species, contrarily to what was found on the south coast of Portugal, where different length-weight relations for both sexes were found, and species showed positive allometric growth ( Olim and Borges 2006 Olim S., Borges T.C. 2006. Weight-length relationships for eight species of the family Triglidae discarded on the south coast of Portugal. J. Appl. Ichthyol. 22: 257-259. https://doi.org/10.1111/j.1439-0426.2006.00644.x ). This can be related to the lower length range sampled in that study (C. cuculus<15.8 cm and T. lyra<24 cm) covering mainly juveniles ( Ismen et al. 2016 Ismen A., Ihsanoglu M.A., Ismen P., Yigin C.C. 2016. Age, growth and reproduction of piper gurnard, Trigla lyra in Saros Bay, the northern Aegean Sea. Rapp. Comm. int. Mer Médit. 41: 340., Neves et al. 2021 Neves A., Sousa I., Sequeira V., et al. 2021. Enhancing knowledge on low-value fishing species: the distinct reproductive strategy of two gurnard species. J. Fish Biol. 1-12. https://doi.org/10.1111/jfb.14849 ). In other geographic areas, C. cuculus showed positive allometric growth with similar length-weight relations between sexes ( Olim and Borges 2006 Olim S., Borges T.C. 2006. Weight-length relationships for eight species of the family Triglidae discarded on the south coast of Portugal. J. Appl. Ichthyol. 22: 257-259. https://doi.org/10.1111/j.1439-0426.2006.00644.x , Marriott et al. 2010 Marriott A.L., Latchford J.W., McCarthy I.D. 2010. Population biology of the red gurnard (Aspitrigla cuculus L.; Triglidae) in the inshore waters of Eastern Anglesey and Northwest Wales. J. Appl. Ichth. 26: 504-512. https://doi.org/10.1111/j.1439-0426.2010.01455.x , Vallisneri et al. 2014 Vallisneri M., Tommasini S., Stagioni M., et al. 2014. Distribution and some biological parameters of the red gurnard, Chelidonichthys cuculus (Actinopterygii, Scorpaeniformes, Triglidae) in the north-central Adriatic Sea. Acta Ichth. Pisc. 44: 173-180. https://doi.org/10.3750/AIP2014.44.3.01 ). These differences are not unusual and may vary within species depending on the season, population and environmental conditions, which may affect the weight of adults through the maturation of gonads or stomach fullness ( Froese 2006 Froese R. 2006. Cube law, condition factor and weight-length relationships: history, meta-analysis and recommendations. J. Appl. Ichthyol. 22: 241-253. https://doi.org/10.1111/j.1439-0426.2006.00805.x ). This variation was minimized in the present study by using eviscerated instead of total weight, but it was not possible to account for differences among studies in the length range of individuals.
Edge analysis and MIR validated the annual periodicity of annulus formation in T. lyra and C. cuculus, which is crucial for age and growth studies. The two age-validation methodologies proved that deposition in both species corresponds to one annulus per year with opaque increments formed in spring/summer. A similar pattern was found for C. cuculus in the Aegean Sea ( Papaconstantinou 1983 Papaconstantinou K.A. 1983. Aspects of the biology of Aspitrigla cuculus (L., 1758) (Pisces, Scorpaeniformes) in the Gulf of Saronikos. Thalassographica 6: 49-75.) and the Central Tyrrhenian Sea ( Colloca et al. 2003 Colloca F., Cardinale M., Ardizzone G.D. 2003. Tracing the life history of red gurnard (Aspitrigla cuculus) using validated otolith annual rings. J. Appl. Ichth. 19: 1-9. https://doi.org/10.1046/j.1439-0426.2003.00342.x ) and is likely associated with higher food availability and therefore more energy available for allocation to growth.
In both species, there was a good separation of the first two increments followed by an increasing overlap with distance to the nucleus. This may be associated with energy investment in both growth and reproduction after maturation, leading to narrower annual increments ( Wootton 1992 Wootton R.J. 1992. Constraints in the evolution of fish life histories. Neth. J. Zool. 42: 291-303. https://doi.org/10.1163/156854291X00342 ). Indeed, length at first maturity estimated for T. lyra off the northern Aegean Sea and C. cuculus off the Portuguese coast varied between 21.9 and 22.7 cm for males ( Ismen et al. 2016 Ismen A., Ihsanoglu M.A., Ismen P., Yigin C.C. 2016. Age, growth and reproduction of piper gurnard, Trigla lyra in Saros Bay, the northern Aegean Sea. Rapp. Comm. int. Mer Médit. 41: 340.) and between 19.9 and 22.1 cm for females ( Neves et al. 2021 Neves A., Sousa I., Sequeira V., et al. 2021. Enhancing knowledge on low-value fishing species: the distinct reproductive strategy of two gurnard species. J. Fish Biol. 1-12. https://doi.org/10.1111/jfb.14849 ), respectively, which according to the present study should correspond to a maturation age of around two years.
Analysis of precision and bias between readers indicated good consistency for both species, with CV and APE values below the limits suggested by Campana (2001) Campana S.E. 2001. Accuracy, precision and quality control in age determination, including a review of the use and abuse of age validation methods. J. Fish Biol. 59: 197-242. https://doi.org/10.1111/j.1095-8649.2001.tb00127.x . Nevertheless, T. lyra otoliths were more difficult to read than those of C. cuculus, with higher values for APE, CV and D and a lower percentage of total agreement between readers. Otoliths of T. lyra displayed some discontinuities and false rings, making it more difficult to identify the annuli, especially in larger individuals, which attained older ages than C. cuculus, and thus made the ageing process more difficult. But whole otoliths proved to be good structures for evaluating age in these species, without the need for additional treatment, such as burning in T. lyra ( Papaconstantinou 1981 Papaconstantinou C. 1981. Age and growth of piper, Trigla lyra, in Saronikos Gulf (Greece). Cybium 5: 73-87., Papaconstantinou et al. 1992 Papaconstantinou C., Petrakis G., Caragitsou E., Mytilineou C. 1992. Preliminary study on the biology of piper (Trigla lyra, L. 1758) in the Aegean Sea. FAO Fisheries Report 477: 127-137.) and crack and burning in C. cuculus ( Baron 1985 Baron J. 1985. Les Triglides (Teleosteens, Scorpaeniformes) de la Baie de Douarnenez. 1. La Croissance de Eutrigla gurnardus, Trigla lucerna, Trigla lastoviza et Aspitrigla cuculus. Cybium 9: 127-144.). The maximum age estimated was 11 years (TL=44.4 cm) for T. lyra and 9 years (TL=42.1 cm) for C. cuculus, values that were slightly above those obtained in other studies ( Table 3 ).
To estimate von Bertalanffy growth parameters, different approaches were needed for the two species. The scarcity of young individuals of T. lyra precluded the adjustment of the von Bertalanffy equation to the age-length key and justified the use of two additional approaches: the back-calculation method and a combined method ( Paiva et al. 2013 Paiva R.B., Neves A., Sequeira V., et al. 2013. Age, growth and mortality of Pontinus kuhlii (Bowdich, 1825) (Scorpaeniformes: Scorpaenidae) in the Gorringe, Ampère, Unicorn and Lion seamounts. Sci. Mar. 77: 95-104. https://doi.org/10.3989/scimar.03632.28B , Gordo et al. 2016 Gordo L.S., Neves A., Vieira A.R., Paiva R.B., Sequeira V. 2016. Age, growth and mortality of the comber Serranus cabrilla (Linnaeus, 1758) in the Eastern Atlantic. Mar. Biol. Res. 12: 656-662. https://doi.org/10.1080/17451000.2016.1169295 ). The combined methodology showed the best adjustment to the data and the most realistic growth parameters.
T. lyra attains a lower TL, and K is naturally higher in this species. Differences in estimated parameters for both species from other European waters were obvious and reflected in different values of ɸ’ ( Table 3 ). T. lyra from the Mediterranean showed lower K values, suggesting slower growth rates than in the present study. Results for C. cuculus are more difficult to interpret given the data variability among different studies, but the significant difference in mean ɸ’ between the Atlantic and Mediterranean populations suggests that they may have different growth strategies, as already observed by Marriott et al. (2010) Marriott A.L., Latchford J.W., McCarthy I.D. 2010. Population biology of the red gurnard (Aspitrigla cuculus L.; Triglidae) in the inshore waters of Eastern Anglesey and Northwest Wales. J. Appl. Ichth. 26: 504-512. https://doi.org/10.1111/j.1439-0426.2010.01455.x . Several factors can contribute to these differences, such as the use of different ageing methodologies, sample heterogeneity, different environmental conditions, latitudes and even fishing pressures.
Mortality is an important parameter to understand population dynamics and one of the most influential parameters in stock assessment. Although evaluation of mortality is of great concern among important exploited target species, this does not always happen for by-catch and low-value species such as gurnards. Given the difficulty of estimating M directly and reliably, two estimators based on one and two parameters were used. As expected, some variation was found depending on the method ( Then et al. 2015 Then A.Y., Hoenig J.M., Hall N.G., Hewitt D.A. 2015. Evaluating the predictive performance of empirical estimators of natural mortality rate using information on over 200 fish species. ICES J. Mar. Sci. 72: 82-92. https://doi.org/10.1093/icesjms/fsu136 ). For the Portuguese coast, T. lyra fishing mortality is lower than for C. cuculus, also facing a lower exploitation rate. The exploitation rate estimated for northwest Wales and eastern Anglesey (E=0.64 yr-1, Marriott et al. 2010 Marriott A.L., Latchford J.W., McCarthy I.D. 2010. Population biology of the red gurnard (Aspitrigla cuculus L.; Triglidae) in the inshore waters of Eastern Anglesey and Northwest Wales. J. Appl. Ichth. 26: 504-512. https://doi.org/10.1111/j.1439-0426.2010.01455.x ) is similar to those obtained in the present study which, according to Gulland (1971) Gulland J.A. 1971. The fish resources of the ocean. Fishing News (Books) Ltd, West Byfleet (UK), 255 pp. are considered high because they are over 0.50 yr-1. Species-specific time series with information on length composition in commercial landings could allow us to evaluate the evolution in exploitation rates over time and infer its influence on age composition, which was not possible given the absence of information. On the other hand, these results are in line with those of other studies because the official landings statistics indicate that C. cuculus is the second most landed gurnard species on the Portuguese coast ( Rocha et al. 2018 Rocha A., Feijó D., Gonçalves P. 2018. Gurnards: species landings’ composition in ICES Division 27.9a. ICES WGWIDE Report 2018. https://doi.org/10.13140/RG.2.2.24277.47841 / Accessed online 11 February 2022.), justifying the higher fishing mortality and exploitation rate estimated in the present study.
Data on age and growth are essential for the understanding of vital traits of species and populations ( Pontual et al. 2002 Pontual H. (de), Panfili J., Wright P.J., Troadec H. 2002. General introduction. In: Panfili J., Pontual H(de), et al. (eds), Manual of fish sclerochronology. Ifremer-lRD coedition, pp. 19-28.), they are the most frequently used life-history parameters to identify putative stocks ( Begg 2005 Begg G.A. 2005. Life history parameters. In: Cadrin S.X., Friedland K.D., Waldman J.R. (eds) Stock Identification Methods: Applications in Fishery Science. Elsevier Academic Press, pp. 119-150. https://doi.org/10.1016/B978-012154351-8/50007-1 ) and, together with mortality, they are important for understanding population dynamics ( Ralston and Williams 1988 Ralston S.V., Williams H.A. 1988. Depth distributions, growth, and mortality of deep slope fishes from the Mariana Archipelago. NOAA-Technical Memorandum, NMFS-SWFC 113: 1-47.). Although commercially important in European waters, gurnards are mainly by-catch, often discarded, and reported under “mixed gurnards”, with incomplete biological and fishing exploitation information by species. The results of this study will contribute to knowledge of the biological parameters of T. lyra and C. cuculus, helping to define stock characteristics and future management strategies for sustainable exploitation, as recommended by the ICES ( ICES 2006 ICES. 2006. Report of the working group on the assessment of new MOU species (WGNEW). ICES (CM Papers and Reports), CM2006/ACFM:11, 240 pp.). However, additional information on the reproductive strategy of these and other gurnard species and their mixed landings should be obtained for better decision-making in fisheries management.