INTRODUCTION
⌅Distribution shifts of marine fishes towards polar areas due to global warming have been reported or predicted in the northeast (Perry et al. 2005Perry A.L., Low P.J., Ellis J.R., Reynolds J.D. 2005. Climate Change and Distribution Shifts in Marine Fishes. Science 308: 1912-1915. https://doi.org/10.1126/science.1111322 , Rjinsdorp et al. 2009Rijnsdorp A.D., Peck M.A., Engelhard G.H., et al. 2009. Resolving the effect of climate change on fish populations. ICES J. Mar. Sci. 66: 1570-1583. https://doi.org/10.1093/icesjms/fsp056 ) and northwest Atlantic (Murawski 1993Murawski S.A. 1993. Climate Change and Marine Fish Distributions: Forecasting from Historical Analogy. Trans. Am. Fish. Soc. 122: 647-658. https://doi.org/10.1577/1548-8659(1993)122<0647:CCAMFD>2.3.CO;2 ), in Australian waters (Cheung et al. 2012Cheung W.W.L., Meeuwig J.J., Feng M., et al. 2012. Climate-change induced tropicalisation of marine communities in Western Australia. Mar. Freshw. Res. 63: 415-427. https://doi.org/10.1071/MF11205 and references therein) and in the Pacific Ocean (Alabia et al. 2018Alabia I.D., García Molinos J., Saitoh S.I., et al. 2018. Distribution shifts of marine taxa in the Pacific Arctic under contemporary climate changes. Divers Distrib. 24: 1583-1597. https://doi.org/10.1111/ddi.12788 , Cheung et al. 2015Cheung W.W.L., Brodeur R.D., Okey T.A., Pauly D. 2015. Projecting future changes in distributions of pelagic fish species of Northeast Pacific shelf seas. Prog. Oceanogr. 130: 19-31. https://doi.org/10.1016/j.pocean.2014.09.003 ). Moreover, changes in marine biodiversity have been foreseen as a result of climate change (Cheung et al. 2009Cheung W.W.L., Lam V.W.Y., Sarmiento J.L., et al. 2009. Projecting global marine biodiversity impacts under climate change scenarios. Fish. Fish. 10: 235-251. https://doi.org/10.1111/j.1467-2979.2008.00315.x ). Variations in temperature, pH and oxygen dissolved in the ocean affect the physiology and dynamics of all marine organisms (Pörtner et al. 2005Pörtner H.O., Storch D., Heilmayer O. 2005. Constraints and trade-offs in climate-dependent adaptation: energy budgets and growth in a latitudinal cline. Sci. Mar. 69: 271-285. https://doi.org/10.3989/scimar.2005.69s2271 , Rijnsdorp et al. 2009Rijnsdorp A.D., Peck M.A., Engelhard G.H., et al. 2009. Resolving the effect of climate change on fish populations. ICES J. Mar. Sci. 66: 1570-1583. https://doi.org/10.1093/icesjms/fsp056 , Ottersen et al. 2010Ottersen G., Kim S., Huse G., et al. 2010. Major pathways by which climate may force marine fish populations. J. Mar. Syst. 79: 343-360. https://doi.org/10.1016/j.jmarsys.2008.12.013 ) and, together with distributional shifts, alter trophic structures, the composition of marine communities and thus ecosystem balances (Vergés et al. 2014Vergés A., Steinberg P.D., Hay M.E., et al. 2014. The tropicalization of temperate marine ecosystems: climate-mediated changes in herbivory and community phase shifts. Proc. R. Soc. B Biol. Sci. 281(1789): 20140846. https://doi.org/10.1098/rspb.2014.0846 ). While any alteration in marine ecosystems will have a major impact on the services they provide, it is especially evident that changes affecting species exploited by fisheries will have a direct impact on the economy and human health because of the importance of fisheries in the economic development of coastal communities and in the provision of healthy food (Hollowed et al. 2013Hollowed A.B., Barange M., Beamish R.J., et al. 2013. Projected impacts of climate change on marine fish and fisheries. ICES J. Mar. Sci. 70: 1023-1037. https://doi.org/10.1093/icesjms/fst081 ).
Atlantic chub mackerel (Scomber colias, Gmelin, 1789) is a medium-sized migratory coastal pelagic fish distributed in warm and temperate waters on both sides of the Atlantic Ocean (Castro and Santana 2000Castro J.J., Santana A.T. 2000. Synopsis of Biological Data on The Chub Marckerel (Scomber japonicus Houttuyn, 1782). FAO Fish Synopsis. 39: 1-77.), including the Macaronesian archipelagos and reaching the Mediterranean and southern Black seas (Collette and Nauen 1983Collette B.B., Nauen C.E. 1983. FAO Species Catalogue Vol. 2 Scombrids of the world an annotated and illustrated catalogue of Tunas, Mackerels, Bonitos and related species know to date, Whitehead et al. 1984Whitehead P.J.P., Bauchot M.L., Hureau J.C., et al. 1984. Scombridae. In: Fishes of the North-eastern Atlantic and the Mediterranean. UNESCO, Paris. 981-997. ISBN: 92-3-002308-6 (v.2), Collette et al. 2011Collette B., Amorim A.F., Boustany A., et al. 2011. Scomber colias. The IUCN Red List of Threatened Species 2011: e.T170357A6767497. https://doi.org/10.2305/IUCN.UK.2011-2.RLTS.T170357A6767497.en ). Although for many years the Atlantic chub mackerel was considered the same species as the Indo-Pacific chub mackerel (Scomber japonicus, Houttuyn 1789), morphologic and genetic studies carried out in the first decade of the 20th century demonstrated that they are two different species (Infante et al. 2007Infante C., Blanco E., Zuasti E., et al. 2007. Phylogenetic differentiation between Atlantic Scomber colias and Pacific Scomber japonicus based on nuclear DNA sequences. Genetica. 130: 1-8. https://doi.org/10.1007/s10709-006-0014-5 , Catanese et al. 2010Catanese G., Manchado M., Infante C. 2010. Evolutionary relatedness of mackerels of the genus Scomber based on complete mitochondrial genomes: Strong support to the recognition of Atlantic Scomber colias and Pacific Scomber japonicus as distinct species. Gene. 452: 35-43. https://doi.org/10.1016/j.gene.2009.12.004 , Cheng et al. 2011Cheng J., Gao T., Miao Z., Yanagimoto T. 2011. Molecular phylogeny and evolution of Scomber (Teleostei: Scombridae) based on mitochondrial and nuclear DNA sequences. Chinese J. Oceanol. Limnol. 29: 297-310. https://doi.org/10.1007/s00343-011-0033-7 ).
In the Canary Islands (FAO area 34.1.2), small pelagic fish are targeted by the artisanal purse-seine fleet, and the Atlantic chub mackerel is the species most commonly caught all year around (FAO 2020aFAO. 2020a. Report of the FAO Working Group on the Assessment of Small Pelagic Fish off Northwest Africa. Casablanca, Morocco, 8-13 July 2019. Casablanca, Morocco.). In the Gulf of Cadiz (FAO area 27.9.a.s), Atlantic chub mackerel is jointly with anchovy and sardine one of three target species fished by the purse seine fleet operating in this area. Additionally, a specifically dedicated seasonal purse-seine fishery targets the species during the second and third quarters of the year, supplying the regional canning industry and providing food for caged bluefin tuna. On the western Portuguese coast (FAO area 27.9.a.c), the northwestern Spanish coast (FAO area 27.9.a.n) and the Cantabrian Sea (FAO area 27.8.c), Atlantic chub mackerel has gained economic importance in recent decades, becoming an important resource for the fishing fleet, especially purse seiners, which are responsible for the bulk of the total landings of this species (ICES 2020ICES. 2020. Workshop on Atlantic chub mackerel (Scomber colias) (WKCOLIAS). https://doi.org/10.17895/ICES.PUB.5970 , 2021aICES. 2021a. Second Workshop on Atlantic Chub Mackerel (Scomber colias) (WKCOLIAS2). https://doi.org/10.17895/ICES.PUB.8142 ). In these last three areas, increasing landings of Atlantic chub mackerel have coincided with the decrease in European sardine (Sardina pilchardus) landings and restrictions on Atlantic mackerel (Scomber scombrus) catches.
Despite its wide distribution, 90% of the Atlantic chub mackerel catches in the northeast Atlantic proceed from northwest African waters (FAO 2020bFAO. 2020b. Fisheries and aquaculture software. FishStatJ - software for fishery statistical time series. In: FAO Fish. Aquac. Dep. https://www.fao.org/fishery/statistics/software/fishstatj/en ). More northerly, in the Atlantic Iberian waters, landings have increased significantly in the most recent years (ICES 2020ICES. 2020. Workshop on Atlantic chub mackerel (Scomber colias) (WKCOLIAS). https://doi.org/10.17895/ICES.PUB.5970 ). This increase has been associated with a higher availability of the fish, as it has expanded its distribrution northwards (Martins et al. 2013Martins M.M., Skagen D., Marques V., et al. 2013. Changes in the abundance and spatial distribution of the Atlantic chub mackerel (Scomber colias) in the pelagic ecosystem and fsheries off Portugal. Sci. Mar. 77: 551-563. https://doi.org/10.3989/scimar.03861.07B , Punzón et al. 2016Punzón A., Serrano A., Sánchez F., et al. 2016. Response of a temperate demersal fish community to global warming. J. Mar. Syst. 161: 1-10. https://doi.org/10.1016/j.jmarsys.2016.05.001 , ICES 2020ICES. 2020. Workshop on Atlantic chub mackerel (Scomber colias) (WKCOLIAS). https://doi.org/10.17895/ICES.PUB.5970 ). This phenomenon has been suggested to be probably driven by the increase in the sea temperature (Costoya et al. 2015Costoya X., de Castro M., Gómez-Gesteira M., Santos F. 2015. Changes in sea surface temperature seasonality in the Bay of Biscay over the last decades (1982-2014). J. Mar. Syst. 150: 91-101. https://doi.org/10.1016/j.jmarsys.2015.06.002 , Jurado-Ruzafa et al. 2019Jurado-Ruzafa A., González-Lorenzo G., Jiménez S., et al. 2019. Seasonal evolution of small pelagic fish landings index in relation to oceanographic variables in the Canary Islands (Spain). Deep Res. Part II Top Stud. Oceanogr. 159: 84-91. https://doi.org/10.1016/j.dsr2.2018.07.002 ) as a result of climate change (Tasker 2008Tasker M.L. 2008. The effect of climate change on the distribution and abundance of marine species in the OSPAR Maritime Area. ICES Coop. Res. Rep., 293. 45 pp., Reid and Valdés 2011Reid P.C., Valdés L. 2011. ICES status report on climate change in the North Atlantic. ICES Cooperative Research Report No. 310. 262. Copenhagen). This distribution expansion makes Atlantic chub mackerel an interesting case study for the analysis of the impact of climate change on exploited species. In fact, Atlantic chub mackerel could compete in the new habitats with other important commercial pelagic species such as the congener Atlantic mackerel (S. scombrus), horse mackerel (Trachurus trachurus), sardine (Sardina pilchardus) and anchovy (Engraulis encrasicolus) (Bachiller and Irigoien 2015Bachiller E., Irigoien X. 2015. Trophodynamics and diet overlap of small pelagic fish species in the bay of biscay. Mar. Ecol. Prog. Ser. 534: 179-198. https://doi.org/10.3354/meps11375 , Garrido et al. 2015aGarrido S., Ben-Hamadou R., Santos A.M.P., et al. 2015a. Born small, die young: Intrinsic, size-selective mortality in marine larval fish. Sci. Rep. 5: 17065. https://doi.org/10.1038/srep17065 , bGarrido S., Silva A., Pastor J., et al. 2015b. Trophic ecology of pelagic fish species off the Iberian coast: Diet overlap, cannibalism and intraguild predation. Mar. Ecol. Prog. Ser. 539: 271-286. https://doi.org/10.3354/meps11506 , Veiga-Malta et al. 2019Veiga-Malta T., Szalaj D., Angélico M.M., et al. 2019. First representation of the tropic structure and functioning o the Portuguese continental shelf ecosystem: insights into the role of sardine. Mar. Ecol. Prog. Ser. 617-618: 323-340 https://doi.org/10.3354/meps12724 ), especially in areas where its presence has increased recently. This competition could affect the fishing yields as well as the pelagic ecosystem balance.
Genetic, morphologic and parasite studies have shown significant differences between Atlantic chub mackerel from the western Atlantic and the eastern Atlantic (Scoles et al. 1998Scoles D.R., Collette B.B., Graves J.E. 1998. Global phylogeography of mackerels of the genus Scomber. Fish. Bull. 96: 823-842., Roldán et al. 2000Roldán M.I., Perrotta R.G., Cortey M., Pla C. 2000. Molecular and morphologic approaches to discrimination of variability patterns in chub mackerel, Scomber japonicus. J. Exp. Mar. Bio. Ecol. 253: 63-74. https://doi.org/10.1016/S0022-0981(00)00244-6 , Costa et al. 2011Costa G., Cavallero S., D’Amelio S., et al. 2011. Helminth parasites of the Atlantic chub mackerel, Scomber colias Gmelin, 1789 from Canary Islands, Central North Atlantic, with comments on their relations with other Atlantic regions. Acta Parasitol. 56: 98-104. https://doi.org/10.2478/s11686-011-0006-1 ), supporting the existence of two different populations at either side of the Atlantic Ocean. However, whether this species constitute a single stock or not in the eastern Atlantic is still unknown, as the studies performed in this area provide conflicting results (Scoles et al. 1998Scoles D.R., Collette B.B., Graves J.E. 1998. Global phylogeography of mackerels of the genus Scomber. Fish. Bull. 96: 823-842., Mele et al. 2014Mele S., Pennino M.G., Piras M.C., et al. 2014. Parasites of the head of Scomber colias (Osteichthyes: Scombridae) from the western Mediterranean Sea. Acta Parasitol. 59: 173-183. https://doi.org/10.2478/s11686-014-0207-5 , Muniz et al. 2020Muniz A.A., Moura A., Triay-Portella R., et al. 2020. Population structure of the chub mackerel (Scomber colias) in the North-east Atlantic inferred from otolith shape and body morphometrics. Mar. Freshw. Res. 72: 341-352. https://doi.org/10.1071/MF19389 ). The information available on the northeast Atlantic suggests that in both African and European waters the species migrates from the wintering areas (mainly located in Mauritanian waters, southern Portugal and the inner part of the Bay of Biscay) towards northern waters in summertime and, in the case of the Bay of Biscay, towards the western Iberian Peninsula (ICES 2021aICES. 2021a. Second Workshop on Atlantic Chub Mackerel (Scomber colias) (WKCOLIAS2). https://doi.org/10.17895/ICES.PUB.8142 ). However, knowledge regarding population dynamics, life-history traits, migration patterns and connectivity of Atlantic chub mackerel in the northeast Atlantic and environmental drivers of their variation is scarce and fragmented. Therefore, large-scale studies such as the preseent one are essential to understand not only the population structure of Atlantic chub mackerelin in the eastern Atlantic but also the influence of climate change on the dynamics of this expanding commercial species and, subsequently, on the ecosystems it colonizes and the fisheries targeting it.
The aim of this study was to analyse the spatial variation of the biological parameters of Atlantic chub mackerel between the traditional areas of distribution (Canary Islands and Gulf of Cadiz) and those of recent expansion further north on the Atlantic coast of the Iberian Peninsula (W Portuguese Coast, NW Spanish waters and Cantabrian Sea), as well as the role of sea surface temperature (SST) as a climate change-linked driver of life-history traits. The results of the present work will be useful for the future assessment of the species in European Atlantic waters, which has been promoted by ICES (2020)ICES. 2020. Workshop on Atlantic chub mackerel (Scomber colias) (WKCOLIAS). https://doi.org/10.17895/ICES.PUB.5970 .
MATERIAL AND METHODS
⌅The study area covered most of the distributional range of Atlantic chub mackerel in the northeast Atlantic (Fig.1), from the Canary Islands (FAO area 34.1.2) to the Cantabrian Sea (FAO area 27.8.c). Table 1 contains information about the characteristics of the sampling procedures in each study area.
The Spanish Institute of Oceanography (IEO, CSIC, Spain) and the Instituto Português do Mar e da Atmosfera (IPMA, Portugal) collect biological data of Atlantic chub mackerel (length, weight and maturity) in the Atlantic Iberian waters from regular (monthly or quarterly) sampling of the commercial fleet (including landings, observers on board and biological sampling in the laboratory), as well as from the acoustic and demersal research surveys carried out within the EU Data Collection Framework. In the Canary Islands, the IEO monitors fleet activity and collects biological information on Atlantic chub mackerel monthly from commercial landings exclusively. For the present study, biological information on Atlantic chub mackerel collected between 2011 and 2019 from the monitoring programmes mentioned above was analysed.
Study area | Subdivision | Period | N | Length range | Age range | Origin |
---|---|---|---|---|---|---|
Cantabrian Sea | 27.8.c. | 2011-2019 | 10173 (4888) | 11-51 | 0-14 | C-S |
NW Spanish Coast | 27.9.a.n | 2011-2019 | 3979 (2018) | 12-42 | 0-6 | C-S |
W Portuguese waters | 27.9.a.c | 2011-2019 | 15075 (7458) | 13-54 | 0-10 | C-S |
Gulf of Cadiz | 27.9.a.s | 2011-2019 | 12318 (5263) | 11-40 | 0-7 | C-S |
Canary Islands | 34.1.2 | 2013-2019 | 6778 (3426) | 13-43 | n.a. | C |
For all collected specimens, total length (TL), total weight (TW) and sex were recorded. TL was recorded to the nearest 0.1 cm and TW to the nearest 0.1 g in the Canary Islands, 0.01 g for laboratory sampling and 0.5 g for scientific surveys in the remaining areas. Maturity stages were determined through macroscopic examination of the gonads in all areas. In the Canary Islands, maturity classification was based on the Holden and Raitt (1974)Holden M.J., Raitt D.F.S. 1974. Manual of fisheries science. Part 2: Methods of resource investigation and their application. FAO Fish. Tech. Rep. 115: maturity key (I, immature; II, maturing virgin or resting; III, maturing; IV, spawning; V, spent), while in the remaining areas it was determined according to the Walsh et al. (1990)Walsh M., Hopkins P., Witthames P.R., et al. 1990. Estimation of total potential fecundity and atresia in the western mackerel stock in 1989. ICES Document CM. 1990/H: 31. 22 pp. maturity key (I, immature; II, maturing; III, mature; IV, spawning; V, partial post-spawning; VI, final post-spawning).
More detailed information about the sampling schedules is provided in the supplementary material, including quarterly length distribution of sampled individuals in each study area (Table S1, Fig. S1).
In all study areas otoliths were collected from both males and females. Whole otoliths were mounted on black plastic slides, concave side up, covered with transparent resin and observed with a binocular microscope under reflected light for age estimation. Age readings from otoliths sampled from fish caught in the Canary Islands did not provide consistent results, so Canary data were not considered in the analyses involving age.
Since no sexual dimorphism has been found in this species (Techetach et al. 2019Techetach M., Ajana R., Saoud Y. 2019. Reproductive parameters of Atlantic chub mackerel Scomber colias in M’diq Bay, Morocco. J. Mar. Biol. Assoc. United Kingdom. 99: 957-962. https://doi.org/10.1017/S0025315418000930 ), all analyses were performed for both sexes combined, except for the study of the onset and duration of the spawning season (see details below).
Length-weight relationship
⌅The length-weight (TL-TW) model used was
where TW i and TL i are the total weight and length of the individual i, respectively, N the total number of individuals, and a and b the parameters to be estimated. The model was linearized taking logarithms (with base 10) on both sides, log10(TWi) = log10(a) + blog10(TLi), for fitting in R software (R Core Team 2021Team R.C. 2021a. R: A language and environment for statistical computing.Team Rs. 2021b. RStudio: Integrated Development for R. RStudio.) using the common function lm for linear models.
In addition to TL, other covariables can be included in the model predictor. In particular, the TL-TW parameters a and b can be specific for each of the K categories of a factor (Note this K has not any relationship with Fulton’s condition factor). Then the model including the interaction between length and the factor is formulated as follows.
where and are the a and b estimates, respectively, for category j (j=1,…,K) of the factor covariable, and the variables I j , j=1,…,K are binary dummy variables which take the value of 1 if individual i corresponds to category j. Then, the model leads to a and b specific estimates for each category of the factor covariable. In the current analysis, the model has been implemented for estimating a specific TL-TW relationship in each study area.
Condition
⌅Condition was approached using Fulton’s K condition index (Fulton 1902Fulton T. W. 1902. The rate of growth of fishes. 20th Annual Report of the Fishery Board of Scotland 1902 (3): 326-446.) to facilitate comparison with previous studies that used this proxy of condition.
Spawning season
⌅The spawning period was determined for each area based on maturity stage data of females, as males sampled for the present study were in the spawning capable stage almost all year around. Furthermore, the existence of spawning capable females is the limiting factor for effective breeding. For the analysis, the proportion of spawning capable over the total number of mature females was the selected index. In the case of the Canary Islands, only females in stage IV (Holden and Raitt’s scale) were considered as spawning capable females, whereas in the rest of areas, they corresponded to females in stages IV and V (Walsh’s maturity scale). Spawning season was defined as the period of consecutive months when the proportion of spawning capable females was higher than 25%.
Length at maturity
⌅The percentage of mature females at length class was fitted to the logistic equation:
where is the predicted mature proportion, a and b the estimated coefficients of the logistic equation and TL the total length.
The length at which 50% of specimens are mature (length at first maturity, L 50) was estimated as the minus ratio of the coefficients (-a/b). All maturity ogive parameters were estimated using the sizeMat (Torrejon-Magallanes 2020Torrejon-Magallanes J. 2020. sizeMat: Estimate Size at Sexual Maturity. https://cran.r-project.org/web/packages/sizeMat/sizeMat.pdf ) package of the R software (R Core Team, 2021Team R.C. 2021a. R: A language and environment for statistical computing.Team Rs. 2021b. RStudio: Integrated Development for R. RStudio.), applying a frequentist generalized linear model with 500 iterations in the boostrap resampling.
Modelling
⌅Geographical variability
⌅Comparisons of length, age, weight and condition between study areas were delivered by the Kruskall-Wallis rank sum test (Kruskal and Wallis 1952Kruskal W.H., Wallis W.A. 1952. Use of Ranks in One-Criterion Variance Analysis. J. Am. Stat. Assoc. 47: 583-621 https://doi.org/10.1080/01621459.1952.10483441 ; kruskal.test function of the R software), because variances were not homogeneus between areas.
For post hoc comparisons between group levels, the Wilcoxon pairwise Rank sum Test was used (pairwise.wilcox.test function of the R software) with the Benjamini and Hochberg (1995)Benjamini Y., Hochberg Y. 1995. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B. 57: 289-300. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x method for multiple testing corrections. These functions are from the stats R package (R Core Team, 2021Team R.C. 2021a. R: A language and environment for statistical computing.Team Rs. 2021b. RStudio: Integrated Development for R. RStudio.).
Temperature effects
⌅To explain the potential geographical differences in biological parameters between study areas, SST was chosen as an explanatory environmental variable because it is directly linked to global warming. SST data for the study areas were downloaded from the SST Aqua MODIS, NPP, 4km, Daytime (11 microns), 2003-2019 (Monthly Composite) dataset at NOAA West Coast Regional Node ERDDAP data server (https://coastwatch.pfeg.noaa.gov/erddap). After data download, a 0.25×0.25-degree grid was used to average the SST, and subsequently the monthly averaged SST was assigned to the centroid of each grid cell. The R script for data download and averaging replication is available at https://git.csic.es/jtornero/erddap-sst-download.
A correlation analysis showed that TL and TW were highly correlated and that the condition factor, K, was also correlated with both length and weight (Fig. S2). We used generalized linear models (GLMs; Wood 2017Wood S.N. 2017. Generalized additive models: An introduction with R. 2nd edition. CRC Press, New York. 496 pp. https://doi.org/10.1201/9781315370279 ) to model the effect of SST and geographical distribution (study area) on TL exclusively, assuming that their variability would be reflected directly in TW and partially in K, because these variables were correlated.
GLMs are an extension of linear models for which the distribution of the response variable can be other than Gaussian. For this reason, a link function, g, is required between the linear predictor and the conditional expectation of the response variable Y, µ(X)=E(Y|X), and the GLM is formulated as
where X=(X 1,…, X p) are the covariables, µ(X)=E(Y|X) is the conditional expectation of the response variable Y, g is a link function between µ(X) and the linear predictor, and ( , ,…, ) are the unknown model parameters.
A GLM assumes that the response variable follows a distribution belonging to the exponential distribution family; in our analysis, it is assumed that our response variable (i.e. TL) follows a Gamma distribution with a natural logarithm link because it is a strictly positive continuous variable.
In GLM, the effects of categorical variables are considered fork-1 of the J factor levels, with the remaining one being considered the base level. Hence the estimated coefficient of each factor level will indicate the deviation from the value of the base level.
The presence of influence points, outliers that greatly affect the regression estimates, was checked using the influence function of the mgcv package (Wood 2011Wood S.N. 2011. Mgcv: GAMs with GCV/AIC/REML smoothness estimation and GAMMs by REML/PQL), which was also used to perform the GLMs.
All statistical analyses were performed with R v.4.0.5 (R Core Team 2021Team R.C. 2021a. R: A language and environment for statistical computing.Team Rs. 2021b. RStudio: Integrated Development for R. RStudio.) and RStudio v.1.4.1717 (RStudio Team 2021Team R.C. 2021a. R: A language and environment for statistical computing.Team Rs. 2021b. RStudio: Integrated Development for R. RStudio.). In all analyses, results were considered statistically significant when the p-value was <0.01. All figures were plotted using the ggplot2 R package (v.3.3.5. Wickham 2016Wickham H. 2016. ggplot2: Elegant graphics for data analysis. 2nd edition. Springer International Publishing, Cham. 260 pp. http://dx.doi.org/10.1007/978-3-319-24277-4 ).
RESULTS
⌅The TL of specimens collected for this study ranged between 11 and 54 cm (between 6 and 1498 g TW), corresponding to estimated ages between 0 and 14 years (Table 1). The largest and oldest specimens were collected on the western Portuguese coast and the Cantabrian Sea.
Results of the Kruskal-Wallis test demonstrated that the medians of all biological parameters (age, TL, TW and Fulton’s K) were significantly different between the study areas (p<0.01. Fig. 2). In fact, the Wilcoxon test revealed that all variables differed significantly between all areas (p<0.01; Fig. S3). In general, age, TL and TW progressively increased northwards from the Canary Islands to the Cantabrian Sea, with the exception of the Gulf of Cadiz, where the values were lower than expected according to the geographical pattern mentioned above. Mean age increased by 86.6% from the Gulf of Cadiz (1.5 years) to the Cantabrian Sea (2.8 years); mean TL decreased by 21.8% from the Canary Islands (25.7 cm) to the Gulf of Cadiz (21.8 cm) and then increased by 43.6% from Cadiz to the Cantabrian Sea (31.3 cm; by 21.8% compared with the Canary Islands); finally, mean TW decreased by 70.2% from Canary Islands (150.1 g) to the Gulf of Cadiz (88.2 g) and increased by 335% from here to the Cantabrian Sea (2955 g; by 96.9% compared with the Canary Islands). Table 2 shows the basic descriptive statistics of the biological variables analysed in this study.
Area: | 27.8.c | ||||||
---|---|---|---|---|---|---|---|
N | mean | sd | median | min | max | se | |
SST | 10055 | 15.85 | 3.02 | 14.4 | 10.88 | 21.95 | 0.03 |
totL | 10173 | 31.31 | 5.99 | 32.1 | 11.2 | 50.8 | 0.06 |
totW | 10163 | 295.46 | 166.07 | 275 | 8 | 1237 | 1.65 |
K | 10163 | 0.84 | 0.12 | 0.84 | 0.15 | 3.25 | 0 |
gonW | 2828 | 4.46 | 9.92 | 1.51 | 0.01 | 128.2 | 0.19 |
Age | 8084 | 2.75 | 1.37 | 3 | 0 | 14 | 0.02 |
Area: | 27.9.a.n | ||||||
N | mean | sd | median | min | max | se | |
SST | 3979 | 16.38 | 2.15 | 16.05 | 12.84 | 20.7 | 0.03 |
totL | 3979 | 28.4 | 4 | 28.9 | 12.5 | 41.5 | 0.06 |
totW | 3977 | 218.45 | 95.89 | 214.34 | 35.6 | 757.71 | 1.52 |
K | 3977 | 0.89 | 0.11 | 0.88 | 0.54 | 4.04 | 0 |
gonW | 1454 | 1.39 | 2.76 | 0.78 | 0.01 | 66.48 | 0.07 |
Age | 2356 | 2.28 | 1.02 | 2 | 0 | 6 | 0.02 |
Area: | 27.9.a.c | ||||||
N | mean | sd | median | min | max | se | |
SST | 14899 | 18.04 | 2.32 | 18.36 | 13.65 | 22.36 | 0.02 |
totL | 15074 | 27.94 | 5.46 | 28.5 | 13.9 | 53.8 | 0.04 |
totW | 15072 | 207.56 | 127.54 | 192.55 | 18.14 | 1497.92 | 1.04 |
K | 15072 | 0.82 | 0.11 | 0.82 | 0.22 | 2.74 | 0 |
gonW | 10662 | 2.49 | 5.01 | 1.02 | 0.01 | 124.45 | 0.05 |
Age | 6072 | 2.71 | 1.76 | 2 | 0 | 10 | 0.02 |
Area: | 27.9.a.s | ||||||
N | mean | sd | median | min | max | se | |
SST | 12318 | 21.05 | 1.7 | 21.74 | 16.08 | 22.84 | 0.02 |
totL | 12318 | 21.75 | 3.87 | 21.4 | 11.4 | 39.9 | 0.03 |
totW | 12318 | 88.22 | 57.12 | 74 | 6.5 | 554 | 0.51 |
K | 12318 | 0.75 | 0.09 | 0.75 | 0.26 | 1.45 | 0 |
gonW | 1418 | 1.46 | 2.77 | 1 | 0.04 | 32.9 | 0.07 |
Age | 3008 | 1.51 | 0.96 | 1 | 0 | 7 | 0.02 |
Area: | 34.1.2 | ||||||
N | mean | sd | median | min | max | se | |
SST | 6721 | 22.1 | 1.81 | 22.1 | 18.73 | 25.66 | 0.02 |
totL | 6778 | 25.65 | 3.81 | 25.2 | 13.2 | 42.5 | 0.05 |
totW | 6778 | 150.06 | 80.44 | 128.9 | 15.4 | 674.1 | 0.98 |
K | 6778 | 0.82 | 0.09 | 0.8 | 0.52 | 1.29 | 0 |
gonW | 6778 | 1.7 | 4.2 | 0.5 | 0.1 | 89.5 | 0.05 |
Age | 0 | NA | NA | NA | NA | NA | NA |
The condition factor K ranged between 0.15 and 4.04. Individuals in the best condition were observed in the northern areas (NW Spanish waters and the Cantabrian Sea), but the geographical trend in K was not as clear as that observed in age, TL and TW. Similar mean K values were obtained in Atlantic chub mackerel from the Cantabrian Sea (27.8.c), the western Portuguese coast (27.9.a.c) and the Canary Islands (34.1.2) ( between 0.82 and 0.84). The highest ( =0.89) and lowest ( =0.75) mean values were found on the northwestern Spanish coast (27.9.a.n) and the Gulf of Cadiz (27.9.a.s), respectively (Table 2). As expected, TW increased with TL and the values fitted to a power relationship, but this relationship differed significantly between study areas (p<0.01, Fig. 3), even when only the overlapping size range was considered. Though the intercept and slope was significantly different between all study areas (p<0.01; Fig. S4), the length-weight relationship was more similar among the groups from the Cantabrian Sea, the western Portuguese coast and the Canary Islands than for the groups from the northwestern Iberian Peninsula and the Gulf of Cadiz, which were more similar to each other, as observed in the condition factor. The final model predicts that the largest specimens from the northwestern Spanish waters and the Gulf of Cadiz are heavier than those from the other three areas.
Spawning activity was detected in all study areas based on monthly variability in the proportion of actively spawning females, although on the northwestern Spanish coast it was negligible (<25%) compared with the remaining areas. In the Canary Islands the peak of spawning activity was recorded in January, when 58% of females were active. In the Gulf of Cadiz a peak of 70% was recorded in February, on the western Portuguese coast a peak of 61% was recorded in March, on the northwestern Spanish coast a peak of 20% was recorded in April-May, and in the Cantabrian Sea a peak of 27% was recorded in May. Assuming that spawning period corresponds to the consecutive months when the proportion of actively spawning females is greater than 25%, a certain prolongation of the spawning season is also observed towards the north. While in the Canary Islands spawning is concentrated in a single month (January), in the Gulf of Cadiz it lasts three months (February-April) and on the western Portuguese coast it lasts six months (January-June). Spawning activity on the northwestern Spanish coast did not reach the threshold established in this study and in the Cantabrian Sea. Though spawning females were detected from January to September, their proportion was less than 25% only in May and July (Fig. 4). In summary, the peak of spawning was delayed and the spawning season extended northwards.
Regarding size at first maturity, significant differences in the maturity ogive were detected between areas in both the intercept and the slope (p<0.01, Fig. 5). Maturity ogives were estimated based on macroscopic maturity classification of ovaries, tso ahus certain degree of innacuracy is expected owing to the difficulty in distinguishing macroscopically between immature and recovering ovaries. The best fit of maturity ogives was obtained in the Cantabrian Sea (R2=0.73), while in the remaining areas the goodness of fit was lower (R2<0.42), especially in the Gulf of Cadiz (Table 3). However, the proportion of mature specimens was significantly related to length in all areas and fitted to the logistic model (p<0.01), although the adjustment can only be considered really good in the Cantabrian Sea.
Study area | a | b | L50 | R2 |
---|---|---|---|---|
27.8.c | -17.84 | 0.71 | 25.23 | 0.73 |
27.9.a.n | -10.21 | 0.47 | 21.70 | 0.41 |
27.9.a.c | -6.91 | 0.40 | 17.49 | 0.36 |
27.9.a.s | -5.71 | 0.24 | 23.98 | 0.19 |
34.1.2 | -13.98 | 0.90 | 15.49 | 0.38 |
All these limitations imply that the results should be interpreted with caution. However, according to our analysis, a progressive increase in the L50 towards the north was observed, from 15.5 cm in the Canary Islands to 25.2 cm in the Cantabrian Sea, with the exception, once again, of the Gulf of Cadiz (L50=24 cm).
The SST in the study areas fluctuated between 10.9°C and 25.7°C (Table 2), with a significant and progressive decrease in mean SST from the Canary Islands (22.1°C) to the Cantabrian Sea (15.9°C) (p<0.01, Fig. S5). Spawning season was detected to take place between 14°C and 21°C in all study areas. More especifically, mean SST during the spawning season decreased progressively from the Canary Islands (20.5°C) to northwestern Spanish waters (14.6ºC), increasing again in the Cantabrian Sea (18.1°C). Our results demonstrated that SST influences total fish length. In fact, our model, which includes the study area, mean annual SST and interaction between them as covariates, explains 33.3% of TL deviance (Table 4). In the Atlantic Iberian waters, the SST-TL relationship showed a positive slope in the Cantabrian Sea (mean SST =15.8°C), which decreased and became negative southwards, coinciding with the increase in SST: 16.4°C on the northwestern Spanish coast, 18°C on the western Portuguese coast and 21°C in the Gulf of Cadiz. Whereas in the Cantabrian Sea and on the northwestern Spanish coast the slope was positive, on the westsern Portuguese coast and in the Gulf of Cadiz it was negative. On the other hand, in the Canary Islands, with the highest mean annual SST (22.1°C), the SST-TL relationship showed the steepest slope (Fig. 6).
Estimate | Std. Error | t value | p-value | |
---|---|---|---|---|
Intercept | 3.34 | 0.009 | 348.8 | < 0.01 |
27.9.a.n | -0.07 | 0.024 | -2.93 | < 0.01 |
27.9.a.c | 0.01 | 0.015 | 0.721 | 0.47 |
27.9.a.s | -0.03 | 0.022 | -1.174 | 0.24 |
34.1.2 | -0.42 | 0.029 | -14.834 | < 0.01 |
SST | 0.01 | 0.001 | 10.655 | < 0.01 |
27.9.a.n:SST | -0.001 | 0.001 | -1.264 | 0.21 |
27.9.a.c:SST | -0.01 | 0.001 | -8.795 | < 0.01 |
27.9.a.s:SST | -0.02 | 0.001 | -15.648 | < 0.01 |
34.1.2:SST | 0.01 | 0.001 | 6.184 | < 0.01 |
R2 | 0.329 | Dev.Expl. | 33.3% |
DISCUSSION
⌅Our results show a clear latitudinal gradient of biological parameters of Atlantic chub mackerel from the Canary Islands (division 34.1.2) to the Cantabrian Sea (27.8.c), with the exception of the Gulf of Cadiz (27.9.a.S), where the estimated values in general break this trend. Both the age range and mean age increase northwards, leading to an increase not only in length and weight, but also in length at first maturity and condition factor, though in the latter case the progression is not as evident as for the other parameters. The smallest, lightest and youngest specimens were collected in the Gulf of Cadiz. Growth was not estimated in this work because age data were not available for all study areas. However, the differences in length suggest that the age structure in Canary waters is likely to be biased towards younger ages than in the Iberian Peninsula groups (with the exception of the Gulf of Cadiz). Although the sampling in the Gulf of Cadiz was biased towards the second half of the year (Fig. S1), the lengths of the specimens caught in this area were smaller than in the other areas in all quarters of the year, so these differences were not due to sampling bias but rather to the biological attributes of this component of the northeast Atlantic population of Atlantic chub mackerel, as we discuss below.
The increase in size with latitude coincides with Bergman’s rule (Bergmann 1847Bergmann C. 1847. Über die Verhältnisse der Wärme ökonomie der Thiere zu ihrer Grösse. Göttinger Stud. 3: 595-708.). Several works have attempted to verify Bergman’s rule on ectothermic species, with contradictory results for fishes (Belk and Houston 2002Belk M.C., Houston D.D. 2002. Bergmann’s Rule in Ectotherms: A Test Using Freshwater Fishes. Am. Nat. 160: 803-808 https://doi.org/10.1086/343880 and references therein). Rypel (2014)Rypel A.L. 2014. The cold-water connection: Bergmann’s rule in North American freshwater fishes. Am. Nat. 183: 147-156. https://doi.org/10.1086/674094 suggested that the expresion of Bergman’s rule in fishes is linked to the species’ thermal niches and is only applicable to coldwater species. However, Atlantic chub mackerel is a temperate species (Castro and Santana 2000Castro J.J., Santana A.T. 2000. Synopsis of Biological Data on The Chub Marckerel (Scomber japonicus Houttuyn, 1782). FAO Fish Synopsis. 39: 1-77.), so our results suggest that Rypel’s hypothesis could also be extended to this type of species or, at least, to those with a wide distribution range, such as Atlantic chub mackerel (Whitehead et al. 1984Whitehead P.J.P., Bauchot M.L., Hureau J.C., et al. 1984. Scombridae. In: Fishes of the North-eastern Atlantic and the Mediterranean. UNESCO, Paris. 981-997. ISBN: 92-3-002308-6 (v.2)). Recent studies show that pelagic fish, including Atlantic mackerel, also belonging to the genus Scomber, follow Bergmann’s rule, confirming that temperature is one of the most important determinants of body size (Hattab et al. 2021Hattab T., Gucu A., Ventero A., et al. 2021. Temperature strongly correlates with regional patterns of body size variation in Mediterranean small pelagic fish species. Mediterr. Mar. Sci. 22(4): 800-811. https://doi.org/10.12681/mms.26525 ). Similarly, latitudinal differences in length at maturity have been observed in the horse mackerel (Trachurus trachurus) in the Mediterranean Sea (Ferreri et al. 2019Ferreri R., McBride R.S., Barra M., et al. 2019. Variation in size at maturity by horse mackerel (Trachurus trachurus) within the central Mediterranean Sea: Implications for investigating drivers of local productivity and applications for resource assessments. Fish. Res. 211, pp.291-299. https://doi.org/10.1016/j.fishres.2018.11.026 ), which some authors have associated with primary production and temperature (Bonanno et al. 2016Bonanno A., Barra M., Basilone G., et al. 2016. Environmental processes driving anchovy and sardine distribution in a highly variable environment: the role of the coastal structure and riverine input. Fish. Oceanogr. 25: 471-490. https://doi.org/10.1111/fog.12166 ; Basilone et al. 2017Basilone G., Mangano S., Pulizzi M., et al. 2017. European anchovy (Engraulis encrasicolus) age structure and growth rate in two contrasted areas of the Mediterranean Sea: the paradox of faster growth in oligotrophic seas. Mediterr. Mar. Sci. 18: 504-516. https://doi.org/10.12681/mms.2059 ). Pörtner et al. (2005)Pörtner H.O., Storch D., Heilmayer O. 2005. Constraints and trade-offs in climate-dependent adaptation: energy budgets and growth in a latitudinal cline. Sci. Mar. 69: 271-285. https://doi.org/10.3989/scimar.2005.69s2271 concluded that eurytherms are able to dynamically adjust the range of their tolerance windows according to temperature fluctuations, so adapted eurytherms show better fitness than native stenotherm species. These authors further suggested that natural selection should favour the adapted eurytherm species, which are energy-efficient, grow fast and reproduce successfully. These results are especially relevant when one considers global warming and its combined impact with fisheries exploitation, as the most commonly exploited pelagic species have shown a drastic decline in growth, condition and size (Van Beveren et al. 2014Van Beveren E., Bonhommeau S., Fromentin J.M., et al. 2014. Rapid changes in growth, condition, size and age of small pelagic fish in the Mediterranean, Mar. Biol. 161: 1809-1822. https://doi.org/10.1007/s00227-014-2463-1 ; Brosset et al. 2017Brosset P., Fromentin J.M., Van Beveren E., et al. 2017. Spatio-temporal patterns and environmental controls of small pelagic fish body condition from contrasted Mediterranean areas, Prog. Oceanogr. 151: 149-162. https://doi.org/10.1016/j.pocean.2016.12.002 ). Given the wide distribution of Atlantic chub mackerel, it can be considered a highly adaptive eurytherm, so its ability to colonize new habitats is likely to be high, and it may continue to expand northwards as long as water becomes warmer within its temperature tolerance or preference, so it is able to adapt to global warming even better than native species.
This type of geographical gradients of life-history traits has been reported in other species from both freshwater and marine habitats (Vila-Gispert et al. 2002Vila-Gispert A., Moreno-Amich R., García-Berthou E. 2002. Gradients of life-history variation: An intercontinental comparison of fishes. Rev. Fish. Biol. Fish. 12: 417-427. https://doi.org/10.1023/A:1025352026974 , Blanck and Lamouroux 2006Blanck A., Lamouroux N. 2007. Large-scale intraspecific variation in life-history traits of European freshwater fish. J. Biogeogr. 34: 862-875. https://doi.org/10.1111/j.1365-2699.2006.01654.x ), as well as different life stages, from early ones (Castro et al. 2002Castro L.R., Llanos A., Blanco J., et al. 2002. Influence of Latitude Variations in Spawning Habitat Characteristics on the Early Life History Traits of the anchoveta, Engraulis ringens, off northern and central Chile. In: Van der Lingen C.D., Roy C., Fréon P., Barange M., Castro L., Gutierrez M., Nykjaer L. and Shillington F. (eds). Report of a GLOBEC-SPACC/IDYLE/ENVIFISH workshop on spatial Approaches to the Dynamics of Coastal Pelagic Resources and their Environment in Upwelling Areas. GLOBEC Report 16: 42-45., Takahashi et al. 2012Takahashi M., McCormick M.I., Munday P.L., Jones G.P. 2012. Influence of seasonal and latitudinal temperature variation on early life-history traits of a coral reef fish. Mar. Freshw. Res. 63: 856-864. https://doi.org/10.1071/MF11278 ) to adults (Huret et al. 2019Huret M., Tsiaras K., Daewel U., et al. 2019. Variation in life-history traits of European anchovy along a latitudinal gradient: a bioenergetics modelling approach. Mar. Ecol. Prog. Ser. 617-618: 95-112. https://doi.org/10.3354/meps12574 ). Hughes et al. (2017)Hughes J.M., Stewart J., Lyle J.M., et al. 2017. Influence of latitudinal variation in environmental gradients and population structure on the demography of a widespread pelagic fish, Arripis trutta (Forster, 1801). Environ. Biol. Fishes. 100: 121-135. https://doi.org/10.1007/s10641-016-0565-y relate the latitudinal gradient of life-history traits to temperature, while other authors relate it to other environmental processes such as upwelling events (Gertseva et al. 2017Gertseva V., Matson S.E., Cope J. 2017. Spatial growth variability in marine fish: example from Northeast Pacific groundfish. ICES J. Mar. Sci. 74: 1602-1613. https://doi.org/10.1093/icesjms/fsx016 ), meso- and large-scale oceanographic processes (Castro et al. 2009Castro L., Claramunt G., Krautz M.., et al. 2009. Egg trait variation in anchoveta Engraulis ringens: a maternal response to changing environmental conditions in contrasting spawning habitats. Mar. Ecol. Prog. Ser. 381: 237-248. https://doi.org/10.3354/meps07922 , Stocks et al. 2014Stocks J.R., Gray C.A., Taylor M.D. 2014. Synchrony and variation across latitudinal gradients: The role of climate and oceanographic processes in the growth of a herbivorous fish. J. Sea Res. 90: 23-32. https://doi.org/10.1016/j.seares.2014.03.002 ), fishing pressure (Gertseva et al. 2010Gertseva V.V., Cope J.M., Matson S.E. 2010. Growth variability in the splitnose rockfish Sebastes diploproa of the northeast Pacific Ocean: pattern revisited. Mar. Ecol. Prog. Ser. 413: 125-136. https://doi.org/10.3354/meps08719 ), food quality and availability (Perrotta et al. 2005Perrotta R.G., Carvalho N., Isidro E. 2005. Comparative study on growth of Chub Mackerel (Scomber japonicus Houttuyn, 1782) from three different regions: NW Mediterranean, NE and SW Atlantic*. Rev. Investig. y Desarr. Pesq. 17: 67-79) or a combination of several of them (Huret et al. 2019Huret M., Tsiaras K., Daewel U., et al. 2019. Variation in life-history traits of European anchovy along a latitudinal gradient: a bioenergetics modelling approach. Mar. Ecol. Prog. Ser. 617-618: 95-112. https://doi.org/10.3354/meps12574 ). The reality is that energy trade-offs in fish are complex and depend on the balance between energy inputs and the metabolic costs of maintenance, growth or maturity. Additionally, all these metabolic processes are strongly influenced by the physical and biological environment in which individuals develop their life stages (Kooijman 2009Kooijman S.A.L.M. 2009. Dynamic energy budget theory for metabolic organisation, third edition. Cambridge University Press. https://doi.org/10.1017/CBO9780511805400 ) and by their life strategies. Because of this, linking geographical clines of biological parameters to specific environmental drivers is difficult, especially in widely distributed species such as the Atlantic chub mackerel.
In the case of pelagic fish, latitudinal changes of biological parameters have been documented in several life-history traits. For example, in Arripis trutta from Australia the number of large and old fish was reported to increase northwards, from the cold waters of Tasmania to the warm waters of southern Queensland, although the initial growth rate was higher in the southern areas (Hughes et al. 2017Hughes J.M., Stewart J., Lyle J.M., et al. 2017. Influence of latitudinal variation in environmental gradients and population structure on the demography of a widespread pelagic fish, Arripis trutta (Forster, 1801). Environ. Biol. Fishes. 100: 121-135. https://doi.org/10.1007/s10641-016-0565-y ). Furthermore, Huret et al. (2019)Huret M., Tsiaras K., Daewel U., et al. 2019. Variation in life-history traits of European anchovy along a latitudinal gradient: a bioenergetics modelling approach. Mar. Ecol. Prog. Ser. 617-618: 95-112. https://doi.org/10.3354/meps12574 observed an increase in growth rate, body size and length at maturity with latitude in European anchovy, and suggested a combination of factors (food availability, maintenance costs, environmental seasonality and temperature) as drivers of these latitudinal gradients.
The present study demonstrates that geographical differences in age and length structure, condition, spawning and maturity of NE Atlantic chub mackerel are significant, showing a clear latitudinal gradient with small and early maturing individuals in the south (the Canary Islands). As the maturity classification was based on macroscopic staging, the results for this parameter must be interpreted with caution; however, considering that we found significant differences between areas in the other biological parameters, the differences in length at first maturity detected in this study can be expected to be real. Latitudinal gradient is also detected in reproductive phenology, with earlier, shorter and more intense spawning activity in the southern study areas than in the northern ones. Latitudinal variability of reproductive season onset and duration has been reported in other fish species too. Huret et al. (2019)Huret M., Tsiaras K., Daewel U., et al. 2019. Variation in life-history traits of European anchovy along a latitudinal gradient: a bioenergetics modelling approach. Mar. Ecol. Prog. Ser. 617-618: 95-112. https://doi.org/10.3354/meps12574 reported a delay and shortening of the reproductive season of European anchovy in more northern latitudes, linked to seasonal variability of temperature and food availability and the ability of individuals to accumulate energy reserves. Barbee et al. (2011)Barbee N.C., Hale R., Morrongiello J., et al. 2011. Large-scale variation in life history traits of the widespread diadromous fish, Galaxias maculatus, reflects geographic differences in local environmental conditions. Mar. Freshw. Res. 62: 790-800. https://doi.org/10.1071/MF10284 found that, in the diadromous fish Galaxias maculatus, spawning season onset occurred earlier at high latitude (New Zealand), where spawning fish were older and larger than on the south coast of Australia. These authors linked this geographical variability of life-history traits to key environmental parameters that vary seasonally, such as temperature, day-light duration and productivity. In our case, latitudinal differences in reproductive phenology could be related not only to temperature and food availability, but also to population structure. However, more analyses are required to corroborate this hypothesis.
The maximum length observed in the present study was between 40 and 54 cm (mean between 21.8 and 31.3 cm), while the maximum age recorded was between 6 and 14 years (mean between 1.5 and 2.8 years). These results are in accordance with the observations of previous studies (Velasco et al. 2011Velasco E.M., del Arbol J., Baro J., Sobrino I. 2011. Edad y crecimiento del estornino Scomber colias del sur de España: una comparación entre muestras procedentes del Atlántico NE y del SW Mediterráneo. Rev. Biol. Mar. Oceanogr. 46: 27-34. https://doi.org/10.4067/S0718-19572011000100004 , Jurado-Ruzafa et al. 2021Jurado-Ruzafa A., Sotillo B., Hernández E., et al. 2021. The Atlantic chub mackerel (Scomber colias) in the Canary Islands (Spain), Fishery and Biological data Update. Second Workshop on Atlantic chub mackerel (Scomber colias) Second Workshop on Atlantic Chub Mackerel (Scomber colias) (WKCOLIAS2). ICES Working Documents WD3. ICES Sci. Rep. 2(20): 148-164., Navarro et al. 2021bNavarro M.R., Landa J., Villamor B., et al. 2021b. First approach to the growth and age corroboration of Northeast Atlantic chub mackerel (Scomber colias) in Northern Iberian waters. Estuar. Coast. Shelf. Sci. 107433. https://doi.org/10.1016/j.ecss.2021.107433 ).
Mean Fulton’s K ranged between 0.82 and 0.89 considering all study areas, except in the Gulf of Cadiz, where K was slightly lower (mean = 0.75). In Atlantic chub mackerel, condition factor is significantly related to TL (Fig. S6). Therefore, the low values observed in the Gulf of Cadiz may be due to the small size of specimens from this area. Fulton’s K values between 0.7 and 1.1 have been reported for the Atlantic chub mackerel from the Atlantic Portuguese coast (Alves 2016Alves M.F. 2016. Survey of parasites of Atlantic Chub Mackerel (Scomber colias) with economic and public health impact: MSc thesis. 1-37., Santos et al. 2017Santos M.J., Castro R., Cavaleiro F., et al. 2017. Comparison of anisakid infection levels between two species of Atlantic mackerel (Scomber colias and S. scombrus) off the Atlantic Portuguese coast. Sci. Mar. 81: 179-185. https://doi.org/10.3989/scimar.04552.26A , Barboza et al. 2020Barboza L.G.A., Lopes C., Oliveira P., et al. 2020. Microplastics in wild fish from North East Atlantic Ocean and its potential for causing neurotoxic effects, lipid oxidative damage, and human health risks associated with ingestion exposure. Sci. Total. Environ. 717: 134625. https://doi.org/10.1016/j.scitotenv.2019.134625 ), suggesting that the mean condition indices obtained in the study period and areas are in the lower part of this range. Fish body condition is closely related to feeding, metabolic rates and the capacity to store energy reserves, i.e. to energetic trade-offs (Saborido-Rey and Kjesbu 2005Saborido-Rey F., Kjesbu O.S. 2005. Growth and maturation dynamics. 26 pp. http://hdl.handle.net/10261/47150 ), all of which are modulated by environmental factors such as temperature or prey availability and diversity, i.e. marine community biodiversity. Changes in the pelagic ecosystem induced by global warming, fishing pressure or any other environmental driver may lead to a decrease in the health status of native pelagic fishes (Shephard et al. 2014Shephard S., Rindorf A., Dickey-Collas M., et al. 2014. Assessing the state of pelagic fish communities within an ecosystem approach and the European Marine Strategy Framework Directive. ICES J. Mar. Sci. 71: 1572-1585. https://doi.org/10.1093/icesjms/fsu005 , Muhling et al. 2017Muhling B., Lindegren M., Clausen L.W., et al. 2017. Impacts of Climate Change on Pelagic Fish and Fisheries. Clim. Chang. Impacts Fish. Aquat. 771-814. https://doi.org/10.1002/9781119154051.ch23 ).
The biological characteristics of Atlantic chub mackerel from the Gulf of Cadiz (young, small and low-condition specimens) should be also addressed from the environmental-induced point of view. In this area, there has been an increase in landings of Atlantic chub mackerel since 2007. While this increase coincides with the recent expansion of the species throughout Ibero-Atlantic waters, it also coincides with the decline of the Ibero-Atlantic sardine stock, which led to severe restrictions on fishing opportunities for this species. These limitations may have led to a change in the behaviour of the fleet, which may have redirected its effort towards Atlantic chub mackerel to compensate for the decline in sardine catches (ICES 2021aICES. 2021a. Second Workshop on Atlantic Chub Mackerel (Scomber colias) (WKCOLIAS2). https://doi.org/10.17895/ICES.PUB.8142 ). This increase in fishing effort could be behind the juvenescence of Atlantic chub mackerel in this area. However, as explained below, the Gulf of Cadiz could also be either a mixing area between the Atlantic and Mediterranean populations or a hinge area in the latitudinal and longitudinal biological gradient observed for Atlantic chub mackerel in the NE Atlantic and Mediterranean.
Regarding the onset and duration of the spawning season, most of the existing studies of Atlantic chub mackerel are based on the analyses of the monthly prevalence of active females and the gonado-somatic index, showing some variability among regions. The spawning season in North African waters has been reported from December to March in Morocco (Techetach et al. 2010Techetach M., Hernando-Casal J.A., Saoud Y., Benajiba M.H. 2010. Reproductive biology of chub mackerel Scomber japonicus in Larache area, Moroccan North Atlantic coast. Cybium 34: 159-165., Wahbi et al. 2011Wahbi F., Errhif A., Ettahiri O. 2011. Cycle de reproduction et variabilité du régime alimentaire du maquereau Scomber japonicus (Houttuyn, 1782) débarqué au port de Casablanca. In: Garcia S., Tandstad M., Caramelo A.M. (eds) Science and Management of Small Pelagic, Fisheries and Aquaculture proceedings. FAO, pp 127-138) and from January to March in Mauritania and Senegal (ICES 2020ICES. 2020. Workshop on Atlantic chub mackerel (Scomber colias) (WKCOLIAS). https://doi.org/10.17895/ICES.PUB.5970 ). In the Atlantic islands the spawning season of Atlantic chub mackerel has been reported from November/December to March in the Canary Islands (Lorenzo and Pajuelo 1996Lorenzo J.M., Pajuelo J.G. 1996. Growth and reproductive biology of chub mackerel Scomber japonicus off the Canary Islands. Southafrican J. Mar. Sci. 17: 275-280. https://doi.org/10.2989/025776196784158635 , Jurado-Ruzafa et al. 2021Jurado-Ruzafa A., Sotillo B., Hernández E., et al. 2021. The Atlantic chub mackerel (Scomber colias) in the Canary Islands (Spain), Fishery and Biological data Update. Second Workshop on Atlantic chub mackerel (Scomber colias) Second Workshop on Atlantic Chub Mackerel (Scomber colias) (WKCOLIAS2). ICES Working Documents WD3. ICES Sci. Rep. 2(20): 148-164.), from January to April in Madeira (Vasconcelos et al. 2012Vasconcelos J., Afonso-Dias M., Faria G. 2012. Atlantic chub mackerel (Scomber colias) spawning season, size and age at first maturity in Madeira waters. Life Mar. Sci. 29: 43-51) and from March to July/August in the Azores (Carvalho et al. 2002Carvalho N., Perrotta R.G., Isidro E. 2002. Age, growth and maturity in the chub mackerel (Scomber japonicus Houttuyn, 1782) from the Azores. Life Mar. Sci. 19A: 93-99.). Finally, in the Atlantic Iberian waters spawning season has been reported from December to March in waters of southern Portugal waters (ICES 2020ICES. 2020. Workshop on Atlantic chub mackerel (Scomber colias) (WKCOLIAS). https://doi.org/10.17895/ICES.PUB.5970 ), from December to May in western Portuguese waters (Nunes et al. 2019Nunes C., Silva A.V., Feijó D., et al. 2019. Atlantic chub mackerel (Scomber colias) growth and reproduction off the Portuguese coast in relation to the population dynamics. Front. Mar. Sci. Conference Abstract: XX Iberian Symposium on Marine Biology Studies (SIEBM XX). https://doi.org/10.3389/conf.fmars.2019.08.00022 ) and, lastly, from March to July on the northwestern Spanish coast and the Cantabrian Sea (Villamor et al. 2017Villamor B., Carrera P., Castro J., et al. 2017. The Chub Mackerel (Scomber colias) in the Atlantic Spanish Waters (ICES divisions 8.c and 9.a): Biological, fishery and survey data. Report of the Working Group on Widely Distributed Stocks (WGWIDE). ICES Working Document WD XVI. ICES CM 2017/ACOM: 23. 873-934. Navarro et al. 2021aNavarro M.R., Domínguez-Petit R., Landa J., et al. 2021a. Preliminary observations on sexual maturity of chub mackerel (Scomber colias) in the Northern Iberian Atlantic waters (ICES Divisions 27.8.c and 27.9.aN). Second Workshop on Atlantic chub mackerel (Scomber colias) (WKCOLIAS2). 11 pp.). The spawning period of Atlantic chub mackerel in the Mediterranean Sea has been reported from April to August (Rizkalla 1998Rizkalla S. 1998. Some biological characters of chub mackerel (Scomber japonicus, Houttuyan, 1782) from the Mediterranean waters of Egypt. Egypt J. Aquat. Biol. Fish. 2: 101-116. https://doi.org/10.21608/ejabf.1998.1627 , Cengiz 2012Cengiz Ö. 2012. Age, growth, mortality and reproduction of the chub mackerel (Scomber japonicus Houttuyn, 1782) from Saros Bay (Northern Aegean Sea, Turkey). Turkish J. Fish. Aquat. Sci. 12: 799-809. https://doi.org/10.4194/1303-2712-v12_4_08 , Allaya et al. 2013Allaya H., Hattour A., Hajjej G., Trabelsi M. 2013. Biologic characteristics of Scomber japonicus (Houttuyn, 1782) in Tunisian waters (Central Mediterranean Sea). Afr. J. Biotechnol. 12: 3040-3048. https://doi.org/10.5897/AJB12.2723 ). Our results show a clear latitudinal gradient of spawning season onset and duration that matches with previous studies. Reviewing all this information together, there seems to be a latitudinal trend from the Strait of Gibraltar towards the poles, but also a longitudinal one from the Strait of Gibraltar towards the Mediterranean Sea and the central Atlantic. Changes in reproductive phenology are usually associated with environmental drivers (Rogers and Dougherty 2019Rogers L.A., Dougherty A.B. 2019. Effects of climate and demography on reproductive phenology of a harvested marine fish population. Glob. Chang. Biol. 25: 708-720. https://doi.org/10.1111/gcb.14483 , Slesinger et al. 2021Slesinger E., Jensen O.P., Saba G. 2021. Spawning phenology of a rapidly shifting marine fish species throughout its range. ICES J. Mar. Sci. 78: 1010-1022. https://doi.org/10.1093/icesjms/fsaa252 ), although they could be related to the age structure of populations, as spawning season onset and duration may change depending on the age of reproductive specimens (Lambert 1987Lambert T.C. 1987. Duration and intensity of spawning in herring Clupea harengus as related to the age structure of the mature population. Mar. Ecol. Prog. Ser. 39: 209-220. https://doi.org/10.3354/meps039209 , Wright and Trippel 2009Wright P.J., Trippel E.A. 2009. Fishery-induced demographic changes in the timing of spawning: consequences for reproductive success. Fish. Fish. 10: 283-304. https://doi.org/10.1111/j.1467-2979.2008.00322.x ). However, with the available data we cannot distinguish whether the phenological differences observed in this study are due to environmental drivers or to the differences in the population structure of each study area. This is especially important in the context of global warming, which is leading to distribution shifts of fish species because reproductive phenology affects population resilience (Lowerre-Barbieri et al. 2017Lowerre-Barbieri S., De Celles G., Pepin P., et al. 2017. Reproductive resilience: a paradigm shift in understanding spawner-recruit systems in exploited marine fish. Fish. Fish. 18: 285-312. https://doi.org/10.1111/faf.12180 ), and more attention should be paid to this issue in the near future.
Regarding the size at first maturity (L50) of Atlantic chub mackerel, our results based on macroscopic observations showed a latitudinal trend, with the lowest value in the Canary Islands (15.5 cm) and increasing progressively towards the north (17.5 cm on the western Portuguese coast and 25.2 cm in the Cantabrian Sea). Once again, the L50 values for the Gulf of Cadiz (24 cm) break with this trend, despite the fact that in general the individuals in this area are the smallest. Previous studies reported higher values in the Bay of Biscay (L50=29 cm; Lucio 1997Lucio P. 1997. Biological aspects of Spanish chub mackerel (Scomber japonicus, Houttuyn, 1782) in the Bay of Biscay from the Basque Country catches. ICES CM l997/BB: 10. 31 pp. https://www.ICES.dk/sites/pub/CM%20Doccuments/1997/BB/1997_BB10.pdf ) and in Portuguese waters (L50=27 cm; Martins 1996Martins M. 1996. New biological data on growth and maturity of Spanish Mackerel (Scomber japonicus) off the Portuguese coast (ICES Division IXa). ICES CM 1996/H: 23, 17 pp.), although more recent investigations showed L50 values in the Bay of Biscay similar to those reported here (25 cm; Villamor et al. 2017Villamor B., Carrera P., Castro J., et al. 2017. The Chub Mackerel (Scomber colias) in the Atlantic Spanish Waters (ICES divisions 8.c and 9.a): Biological, fishery and survey data. Report of the Working Group on Widely Distributed Stocks (WGWIDE). ICES Working Document WD XVI. ICES CM 2017/ACOM: 23. 873-934., Navarro et al. 2021aNavarro M.R., Domínguez-Petit R., Landa J., et al. 2021a. Preliminary observations on sexual maturity of chub mackerel (Scomber colias) in the Northern Iberian Atlantic waters (ICES Divisions 27.8.c and 27.9.aN). Second Workshop on Atlantic chub mackerel (Scomber colias) (WKCOLIAS2). 11 pp.). In western Portuguese waters, Nunes et al. (2019)Nunes C., Silva A.V., Feijó D., et al. 2019. Atlantic chub mackerel (Scomber colias) growth and reproduction off the Portuguese coast in relation to the population dynamics. Front. Mar. Sci. Conference Abstract: XX Iberian Symposium on Marine Biology Studies (SIEBM XX). https://doi.org/10.3389/conf.fmars.2019.08.00022 reported an L50 equal to 22.6 cm and 19.2 cm for females and males, respectively; while Gonçalves et al. (2016)Gonçalves J.M.S., Blanc N., Brandão C., et al. 2016. Valorização de recursos pesqueiros: Cavala Algarvia. Relatório final. Universidade do Algarve, CCMAR, Faro. 44pp.+Anexos. http://hdl.handle.net/10400.1/8699 in southern Portugal reported an L50 for both sexes combined of 18.6 cm. These Portuguese values are slightly higher than our estimations (17.5 cm). In the case of the Canary Islands, our estimate of L50 is the smallest one and lower than previous estimations reported in the same area (18-19 cm; Lorenzo and Pajuelo 1996Lorenzo J.M., Pajuelo J.G. 1996. Growth and reproductive biology of chub mackerel Scomber japonicus off the Canary Islands. Southafrican J. Mar. Sci. 17: 275-280. https://doi.org/10.2989/025776196784158635 , Jurado-Ruzafa et al. 2021Jurado-Ruzafa A., Sotillo B., Hernández E., et al. 2021. The Atlantic chub mackerel (Scomber colias) in the Canary Islands (Spain), Fishery and Biological data Update. Second Workshop on Atlantic chub mackerel (Scomber colias) Second Workshop on Atlantic Chub Mackerel (Scomber colias) (WKCOLIAS2). ICES Working Documents WD3. ICES Sci. Rep. 2(20): 148-164.). This difference is likely due to differences in the data used to calculate the L50, which in previous studies included post-spawning and spawning specimens collected exclusively during the spawning period. In Atlantic Moroccan waters, L50 values increased from 23 cm in the north (Techetach et al. 2010Techetach M., Hernando-Casal J.A., Saoud Y., Benajiba M.H. 2010. Reproductive biology of chub mackerel Scomber japonicus in Larache area, Moroccan North Atlantic coast. Cybium 34: 159-165.) to 25 to 27.5 cm in the south (Wahbi et al. 2017Wahbi F. 2017. Bio-écologie de Trachurus trachurus (Linnaeus, 1758) et de Scomber colias Gmelin, 1789 dans l’écosystème pélagique de la zone atlantique sud marocaine (21oN - 26o30’N). PhD thesis, Univ. Hassan II, Faculté des Sciences Ain Chock, Casablanca. 230 pp.). High values of L50 were reported in the Azores (27.8 cm; Carvalho et al. 2002Carvalho N., Perrotta R.G., Isidro E. 2002. Age, growth and maturity in the chub mackerel (Scomber japonicus Houttuyn, 1782) from the Azores. Life Mar. Sci. 19A: 93-99.), followed by Madeira (22 cm; Vasconcelos et al. 2012Vasconcelos J., Afonso-Dias M., Faria G. 2012. Atlantic chub mackerel (Scomber colias) spawning season, size and age at first maturity in Madeira waters. Life Mar. Sci. 29: 43-51) and the Mediterranean Sea (16.8 and 19 cm; Cengiz 2012Cengiz Ö. 2012. Age, growth, mortality and reproduction of the chub mackerel (Scomber japonicus Houttuyn, 1782) from Saros Bay (Northern Aegean Sea, Turkey). Turkish J. Fish. Aquat. Sci. 12: 799-809. https://doi.org/10.4194/1303-2712-v12_4_08 , Cikeš and Zorica 2012Cikeš Keč V., Zorica B. 2012. The reproductive traits of Scomber japonicus (Houttuyn, 1782) in the Eastern Adriatic Sea. J. Appl. Ichthyol. 28: 15-21. https://doi.org/10.1111/j.1439-0426.2011.01893.x , Techetach et al. 2019Techetach M., Ajana R., Saoud Y. 2019. Reproductive parameters of Atlantic chub mackerel Scomber colias in M’diq Bay, Morocco. J. Mar. Biol. Assoc. United Kingdom. 99: 957-962. https://doi.org/10.1017/S0025315418000930 ). Spatial patterns of size at first maturity do not seem to be as clear as spawning onset and duration, although the L50 values presented here should be interpreted with caution, as stated above. Nevertheless, certain latitudinal and longitudinal trends of maturity patterns can be glimpsed and should be corroborated by histological studies. In any case, the spatial differences observed in size and age structure and condition suggest the existence of different life-history patterns, so spatial differences in maturation can also be expected. This aspect would be especially relevant for stock assessment purposes.
The spatial population structure of the Atlantic chub mackerel in the eastern Atlantic, the migration processes through the distribution range and the connectivity, including the Mediterranean Sea, remain unknown (ICES 2021aICES. 2021a. Second Workshop on Atlantic Chub Mackerel (Scomber colias) (WKCOLIAS2). https://doi.org/10.17895/ICES.PUB.8142 ). The geographical differences in life-history traits estimated in this and previous studies suggest that there may be at least two population components: one in the Canary Islands and one in the Iberian Peninsula. Furthermore, the Iberian unit could include two subcomponents in the Cantabrian Sea (Bay of Biscay) and the Portuguese coast, whose degree of connection could vary depending on the abundance and migration intensity of Atlantic chub mackerel in each area (ICES 2020ICES. 2020. Workshop on Atlantic chub mackerel (Scomber colias) (WKCOLIAS). https://doi.org/10.17895/ICES.PUB.5970 and references therein). Additionally, the Gulf of Cadiz could be considered a hinge area, or even a mixing area between the northeast and central east Atlantic populations of Atlantic chub mackerel as well as between the notheast Atlantic and Mediterranean stocks, as suggested by ICES (2021a)ICES. 2021a. Second Workshop on Atlantic Chub Mackerel (Scomber colias) (WKCOLIAS2). https://doi.org/10.17895/ICES.PUB.8142 . The results of the growth analysis carried out by Velasco et al. (2011)Velasco E.M., del Arbol J., Baro J., Sobrino I. 2011. Edad y crecimiento del estornino Scomber colias del sur de España: una comparación entre muestras procedentes del Atlántico NE y del SW Mediterráneo. Rev. Biol. Mar. Oceanogr. 46: 27-34. https://doi.org/10.4067/S0718-19572011000100004 support the hypothesis of a mixing area, as they found no differences in growth between the Gulf of Cadiz and the Alboran Sea, suggesting that the Strait of Gibraltar is not a geographical barrier for this species. However, in general, results of studies targeting the stock structure of Atlantic chub mackerel in the northeast Atlantic and the Mediterranean Sea show differences depending on the methodological approach. Some genetic studies showed no significant differences between the Atlantic chub mackerel from the Mediterranean and eastern Atlantic areas (Scoles et al. 1998Scoles D.R., Collette B.B., Graves J.E. 1998. Global phylogeography of mackerels of the genus Scomber. Fish. Bull. 96: 823-842., Zardoya et al. 2004Zardoya R., Castilho R., Grande C., et al. 2004. Differential population structuring of two closely related fish species, the mackerel (Scomber scombrus) and the chub mackerel (Scomber japonicus), in the Mediterranean Sea. Mol. Ecol. 13: 1785-1798. https://doi.org/10.1111/j.1365-294X.2004.02198.x ). However, significant regional differences have been found between smaller areas in the east Atlantic Ocean, including the Mediterranean Sea and Macarronesian islands, based on the analysis of morphology, meristic characteristics (Allaya et al. 2016Allaya H., Zrelli S, Hajjej G. 2016. Identification of Atlantic Chub mackerel Scomber colias population through the analysis of body shape in Tunisian waters. Cah. Biol. Mar. 57: 195-207, Bouzzammit and El Ouizgami 2019Bouzzammit N., El Ouizgani H. 2019. Morphometric and meristic variation in the Atlantic chub mackerel Scomber colias Gmelin, 1789 from the Moroccan coast. Indian J. Fish. 66: 8-15. https://doi.org/10.21077/ijf.2019.66.2.78488-02 , Muniz et al. 2020Muniz A.A., Moura A., Triay-Portella R., et al. 2020. Population structure of the chub mackerel (Scomber colias) in the North-east Atlantic inferred from otolith shape and body morphometrics. Mar. Freshw. Res. 72: 341-352. https://doi.org/10.1071/MF19389 ) and associated parasites (Mele et al. 2014Mele S., Pennino M.G., Piras M.C., et al. 2014. Parasites of the head of Scomber colias (Osteichthyes: Scombridae) from the western Mediterranean Sea. Acta Parasitol. 59: 173-183. https://doi.org/10.2478/s11686-014-0207-5 ). Muniz et al. (2020)Muniz A.A., Moura A., Triay-Portella R., et al. 2020. Population structure of the chub mackerel (Scomber colias) in the North-east Atlantic inferred from otolith shape and body morphometrics. Mar. Freshw. Res. 72: 341-352. https://doi.org/10.1071/MF19389 , based on otolith morphology analysis, suggested the existence of one group in the northeast Atlantic islands (Azores, Madeira and the Canary Islands) and another one on the Iberian Portuguese coast. Correia et al (2021)Correia A.T., Moura A., Triay-Portella R., et al. 2021. Population structure of the chub mackerel (Scomber colias) in the NE Atlantic inferred from otolith elemental and isotopic signatures. Fish. Res. 234: 105785. https://doi.org/10.1016/j.fishres.2020.105785 , according to microchemical analysis of the same otoliths, suggested the existence of at least four stock components in the northeast Atlantic (Azores, Madeira, Canary Islands and Iberian Portuguese coast). There is no information on the migratory behaviour of this species in the study area, although by analogy with the previous study on the species in northwest African waters (García 1982García S. 1982. Distribution, migration and spawning of the main fish resources in the Northern CECAF Area. Rome), we assume that the Atlantic chub mackerel performs annual reproductive migrations all along its distributional range. If so, connectivity between the different population components may exist, but the extent of this migration and therefore the degree of potential connectivity remain unknown. Many of these studies have been approached on a local or regional scale, and thus lack the global vision required to investigate the population structure of a species as widely distributed as the Atlantic chub mackerel. Jansen et al. (2013)Jansen T., Gislason H., Goldstien S.J. 2013. Population Structure of Atlantic Mackerel (Scomber scombrus). PloS One. 8(5): e64744. https://doi.org/10.1371/journal.pone.0064744 proposed for the northeast Atlantic mackerel (S. scombrus) that the population structure is more a dynamic cline rather than connected contingents, and this could also be the case of the Atlantic chub mackerel. However, more studies targeting migration behaviour, genetics, morphometry and meristic characterization from a global perspective should be carried out to confirm meaningful biological management units for assessment purposes.
Considering that the increase in abundance of Atlantic chub mackerel in the Atlantic Iberian waters is assumed to be a relatively recent event (Villamor et al. 2017Villamor B., Carrera P., Castro J., et al. 2017. The Chub Mackerel (Scomber colias) in the Atlantic Spanish Waters (ICES divisions 8.c and 9.a): Biological, fishery and survey data. Report of the Working Group on Widely Distributed Stocks (WGWIDE). ICES Working Document WD XVI. ICES CM 2017/ACOM: 23. 873-934.), it is striking that individuals in the northernmost (i.e. most recently occupied) study areas show the greatest diversity of age and size and the best condition. There are two possible explanations for this: a) the presence of older and larger individuals of Atlantic chub mackerel could be due to the fact that they have a greater capacity for movement and therefore a greater capacity for expansion to new suitable habitats in a global warming scenario; and/or b) fishing pressure or predation could be less intense in the north than in northwest African waters, reducing the presence of larger and older individuals in this area. With the information available, we cannot determine the causes of the age and size structure observed in each study area. According to Saunders and Tarling (2018)Saunders R.A., Tarling G.A. 2018. Southern ocean mesopelagic fish comply with Bergmann’s Rule. Am. Nat. 191: 343-351. https://doi.org/10.1086/695767 , greater body size in mesopelagic fishes is a necessary attribute to reach colder regions, which supports our first hypothesis that the largest individuals are those with the greatest migratory capacity.
In most studied fish species, great size, age and condition are associated with high reproductive potential owing to the existence of maternal effects (Green 2008Green B.S. 2008. Chapter 1: Maternal Effects in Fish Populations. Adv. Mar. Biol. 54: 1-105. https://doi.org/10.1016/S0065-2881(08)00001-1 ). In addition, a prolonged reproductive season increases the probability of offspring finding optimal environmental windows for their development, thus increasing their survival (Mertz and Myers 1994Mertz G., Myers R.A. 1994. Match/mismatch predictions of spawning duration versus recruitment variability. Fish. Oceanogr. 3: 236-245. https://doi.org/10.1111/j.1365-2419.1994.tb00101.x ). Given the biological characteristics of the groups analysed in this study and the duration of the reproductive season in each area, the geographical variability of the reproductive potential of this species and its impact on stock resilience should be analysed in the future. This variability may be a mechanism of adaptation to the environmental conditions (temperature, productivity and turbulence) that could favour or jeopardize offspring survival and hence recruitment.
Our model indicates that SST does indeed play a relevant role in the size of Atlantic chub mackerel and, presumably, on its growth. While on the western Portuguese coast and the Gulf of Cadiz this effect is negative, in the remaining areas it is positive. The rates of physiological processes usually increase with temperature up to a certain limit, above which the temperature effect becomes negative. This temperature effect is also observed in the growth of fish, so at low temperatures growth increases as temperature increases until a tipping point (optimum temperature) when growth decreases with increasing temperature (Jobling 1995Jobling M. 1995. The influence of environmental temperature on growth and conversion efficiency in fish. 1995: 1-26.). The temperature range at which fish populations exhibit this behaviour changes geographically owing to acclimatation mechanisms (Sunday et al. 2011Sunday J.M., Bates A.E., Dulvy N.K. 2011. Global analysis of thermal tolerance and latitude in ectotherms. Proc. R. Soc. B Biol. Sci. 278: 1823-1830. https://doi.org/10.1098/rspb.2010.1295 ). Therefore, growth is expected to increase with temperature in areas that are below their optimum temperature and to decrease in areas that have exceeded this threshold. Assuming the total length as a proxy of the growth rate of Atlantic chub mackerel for each study area, and based on the slope of the SST-TL relationship, we could conclude that S. colias in the Cantabrian Sea and on the northwestern Spanish coast is farther from its optimum temperature than it is on the western Portuguese coast, where it is slightly above it, and in the Gulf of Cadiz, where it is clearly above it. In the Canary Islands, where mean annual SST is considerably higher than in the rest of the study areas, the effect of SST on the TL of Atlantic chub mackerel was the most markedly positive. We do not know the reasons behind this relationship but it could be related to the comparatively small specimens of Atlantic chub mackerel inhabiting the Canary waters. Further knowledge of this relationship between SST and growth would help to better understand the dynamics and population structure of Atlantic chub mackerel in the northeast Atlantic in order to properly calibrate the assessment and manage this resource adequately under future global warming scenarios.
Many fish species have been reported to be moving to higher latitudes as a consequence of climate change (Perry et al. 2005Perry A.L., Low P.J., Ellis J.R., Reynolds J.D. 2005. Climate Change and Distribution Shifts in Marine Fishes. Science 308: 1912-1915. https://doi.org/10.1126/science.1111322 , Rijnsdorp et al. 2009Rijnsdorp A.D., Peck M.A., Engelhard G.H., et al. 2009. Resolving the effect of climate change on fish populations. ICES J. Mar. Sci. 66: 1570-1583. https://doi.org/10.1093/icesjms/fsp056 ). One of the most emblematic cases in recent decades is that of the Atlantic mackerel (S. scombrus), whose distribution has expanded northwards and westwards in the NE Atlantic Ocean, reaching the waters of Svalvard and Iceland (Berge et al. 2015Berge J., Heggland K., Lønne O.J., et al. 2015. First Records of Atlantic Mackerel (Scomber scombrus) from the Svalbard Archipelago, Norway, with Possible Explanations for the Extensions of Its Distribution, Arctic 68: 54-61. https://doi.org/10.14430/arctic4455 ) driven by stock size and temperature (Astthorsson et al. 2012Astthorsson O.S., Valdimarsson H., Gudmundsdottir A., Óskarsson G.J. 2012. Climate-related variations in the occurrence and distribution of mackerel (Scomber scombrus) in Icelandic waters. ICES J. Mar. Sci. 69: 1289-1297. https://doi.org/10.1093/icesjms/fss084 , Olafsdottir et al. 2019Olafsdottir A.H., Utne K.R., Jacobsen J.A., et al. 2019. Geographical expansion of Northeast Atlantic mackerel (Scomber scombrus) in the Nordic Seas from 2007 to 2016 was primarily driven by stock size and constrained by low temperatures. Deep Sea Res. Part II Top Stud. Oceanogr. 159: 152-168. https://doi.org/10.1016/j.dsr2.2018.05.023 ). Similarly, landings of Atlantic chub mackerel in the Atlantic Iberian waters have increased during the last few decades (Martins et al. 2013Martins M.M., Skagen D., Marques V., et al. 2013. Changes in the abundance and spatial distribution of the Atlantic chub mackerel (Scomber colias) in the pelagic ecosystem and fsheries off Portugal. Sci. Mar. 77: 551-563. https://doi.org/10.3989/scimar.03861.07B , Villamor et al. 2017Villamor B., Carrera P., Castro J., et al. 2017. The Chub Mackerel (Scomber colias) in the Atlantic Spanish Waters (ICES divisions 8.c and 9.a): Biological, fishery and survey data. Report of the Working Group on Widely Distributed Stocks (WGWIDE). ICES Working Document WD XVI. ICES CM 2017/ACOM: 23. 873-934., ICES 2020ICES. 2020. Workshop on Atlantic chub mackerel (Scomber colias) (WKCOLIAS). https://doi.org/10.17895/ICES.PUB.5970 ), likely associated with an increase in stock size and its northward expansion linked to global warming. The arrival of new species in ecosystems can have a major impact on trophic balances and the dynamics of other species (Cheung et al. 2009Cheung W.W.L., Lam V.W.Y., Sarmiento J.L., et al. 2009. Projecting global marine biodiversity impacts under climate change scenarios. Fish. Fish. 10: 235-251. https://doi.org/10.1111/j.1467-2979.2008.00315.x , Hollowed et al. 2013Hollowed A.B., Barange M., Beamish R.J., et al. 2013. Projected impacts of climate change on marine fish and fisheries. ICES J. Mar. Sci. 70: 1023-1037. https://doi.org/10.1093/icesjms/fst081 , Muhling et al. 2017Muhling B., Lindegren M., Clausen L.W., et al. 2017. Impacts of Climate Change on Pelagic Fish and Fisheries. Clim. Chang. Impacts Fish. Aquat. 771-814. https://doi.org/10.1002/9781119154051.ch23 ). The Atlantic chub mackerel is an opportunistic pelagic species that feeds mainly on euphausiids and decapod crustaceans, but also significantly on eggs, larvae and juveniles of other pelagic and demersal fish species (Torres et al. 2013Torres M.A., Coll M., Heymans J.J., Christensen V., Sobrino I. 2013. Food-web structure of and fishing impacts on the Gulf of Cadiz ecosystem (South-western Spain). Ecol. Modell. 265: 26-44. https://doi.org/10.1016/j.ecolmodel.2013.05.019 ). In fact, Garrido et al. (2015a)Garrido S., Ben-Hamadou R., Santos A.M.P., et al. 2015a. Born small, die young: Intrinsic, size-selective mortality in marine larval fish. Sci. Rep. 5: 17065. https://doi.org/10.1038/srep17065 identified sardine eggs as an important component of the diet of Atlantic chub mackerel and reported a considerable overlap in the prey spectrum between European sardine and juvenile chub mackerel. Taking this into account, it is probable that the expansion of Atlantic chub mackerel in Atlantic Iberian waters played an important role in the decline of Iberian sardine populations observed in the same period. To corroborate this hypothesis, the trends of both species in the Atlantic Iberian waters in the coming years should be monitored, especially when the sardine population has started to recover (ICES 2021bICES. 2021b. Working Group on Southern Horse Mackerel Anchovy and Sardine (WGHANSA).). Another question that would be worth investigating further is whether the change in the distribution of S. scombrus could have reduced the occupation of certain pelagic niches in Atlantic Iberian waters, which could be being occupied by warmer-water species such as the Atlantic chub mackerel. Furthermore, its expansion in Atlantic Iberian waters may also impact on predators, as it has been reported as an important prey for other commercial species such as European hake, sharks, bluefin tuna and marine mammals (Torres et al. 2013Torres M.A., Coll M., Heymans J.J., Christensen V., Sobrino I. 2013. Food-web structure of and fishing impacts on the Gulf of Cadiz ecosystem (South-western Spain). Ecol. Modell. 265: 26-44. https://doi.org/10.1016/j.ecolmodel.2013.05.019 , Varela et al. 2013Varela J.L., Rodríguez-Marín E., Medina A. 2013. Estimating diets of pre-spawning Atlantic bluefin tuna from stomach content and stable isotope analyses. J. Sea Res. 76: 187-192. https://doi.org/10.1016/j.seares.2012.09.002 , Giménez et al. 2017Giménez J., Marçalo A., Ramírez F., et al. 2017. Diet of bottlenose dolphins (Tursiops truncatus) from the Gulf of Cadiz: Insights from stomach content and stable isotope analyses. PLoS One 12: e0184673. https://doi.org/10.1371/journal.pone.0184673 ).
This paper outlines the population structure of the Atlantic chub mackerel in the northeast Atlantic from a broad geographical perspective and provides information on the regional variability of key biological parameters commonly required for assessment, such as length-weight relationships, the maturity ogive and size composition. However, there are still many gaps to be filled. In relation to population structure, it would be desirable to extend the study to North African and western Mediterranean waters in order to carry out not only analyses of life-history parameters such as growth and maturation, but also other stock identification studies such as genetic, morphometric and parasite analyses with a wider perspective. Similarly, it would be necessary to analyse how the biological differences detected between the study areas affect the productivity and resilience of the species. This would require time series of fecundity, egg and larval quality and recruitment. Additionally, analysing biodiversity changes in the pelagic ecosystem of the study areas would allow us to understand the impact of new species arrivals on fisheries yields. Ultimately, from an ecological point of view, this species is a promising case study for monitoring the impact of climate change on pelagic ecosystems.