Modelling the spatio-temporal distribution of age-1 Bay of Biscay anchovy (Engraulis encrasicolus) at spawning time


  • Leire Ibaibarriaga AZTI-Tecnalia, Marine Research Division
  • Andrés Uriarte AZTI-Tecnalia, Marine Research Division
  • Urtzi Laconcha AZTI-Tecnalia, Marine Research Division
  • Miguel Bernal General Fisheries Commission for the Mediterranean (GFCM), Food and Agriculture Organization of the United Nations (FAO)
  • María Santos AZTI-Tecnalia, Marine Research Division
  • Marina Chifflet AZTI-Tecnalia, Marine Research Division
  • Xabier Irigoien King Abdullah University of Science and Technology (KAUST)



spatial age composition, anchovy, Bay of Biscay, daily egg production method, generalized additive model, environmental effects


To efficiently manage and maintain fish stock abundance and age structure, it is necessary to understand variability in its spatial distribution. The main objective of this study was to describe the demographic structure of the Bay of Biscay anchovy during the main spawning period based on adult samples from 15 daily egg production method surveys. The proportion of age-1 individuals was modelled using generalized additive models with a binomial distribution and a logit link in relation to geographical and environmental covariates. The possibility of obtaining different models depending on the average age-1 proportion (low or high) was also explored. In general, age-1 individuals were found in shallow waters close to the coast, especially associated with the Gironde and Adour river plumes, whereas older individuals were prevalent on the shelf break and in oceanic waters. However, in years with a high age-1 proportion the younger individuals were also dominant in oceanic waters. These results could be used for management purposes, such as defining protected areas for particular age groups.


Download data is not yet available.


Albaina A., Irigoien X. 2007. Fine scale zooplankton distribution in the Bay of Biscay in spring 2004. J. Plankton Res. 29: 851-870.

Aldanondo N., Cotano U., Tiepolo M., Boyra G., Irigoien X. 2010. Growth and movement patterns of early juvenile European anchovy (Engraulis encrasicolus L.) in the Bay of Biscay based on otolith microstructure and chemistry. Fish. Oceanogr. 19: 196-208.

Allain G., Petitgas P., Lazure P., Grellier P. 2007. Biophysical modelling of larval drift, growth and survival for the prediction of anchovy (Engraulis encrasicolus) recruitment in the Bay of Biscay (NE Atlantic). Fish. Oceanogr. 16: 489-505.

Bailey M.C., Maravelias C. D., Simmonds E.J. 1998. Changes in the spatial distribution of autumn spawning herring (Clupea harengus L.) derived from annual acoustic surveys during the period 1984-1996. ICES J. Mar. Sci. 55: 545-555.

Barange M., Coetzee J.C., Twatwa N.M. 2005. Strategies of space occupation in anchovy and sardine in the southern Benguela: the role of stock size and intra-species competition. ICES J. Mar. Sci. 62: 645-654.

Begg G.A., Marteinsdottir G. 2002. Environmental and stock effects on spatial distribution and abundance of mature cod Gadus morhua. Mar. Ecol. Prog. Ser. 229: 245-262.

Bellier E., Planque B., Petitgas P. 2007. Historical fluctuations in spawning location of anchovy (Engraulis encrasicolus) and sardine (Sardina pilchardus) in the Bay of Biscay during 1967-73 and 2000-2004. Fish. Oceanogr. 16: 1-15.

Bernal M., Stratoudakis Y., Coombs S., Angelico M.M. Lago de Lanzós A., Porteiro C., Sagarminaga Y., Santos M., Uriarte A., Cunha E., Valdés L., Borchers D. 2007. Sardine spawning off the European Atlantic coast: Characterization of and spatiotemporal variability in spawning habitat. Prog. Oceanogr. 74: 210-227.

Boyra G., Martínez U., Cotano U., Santos M., Irigoien X., Uriarte A. (In press). Acoustic surveys for juvenile anchovy in the Bay of Biscay: abundance estimate as an indicator of the next year's recruitment and spatial distribution patterns. ICES J. Mar. Sci.

Brochier T., Colas F., Lett C., Echevin V., Cubillos L.A., Tam, J., Chlaida M., Mullon C., Fréon P. 2009. Small pelagic fish reproductive strategies in upwelling systems: A natal homing evolutionary model to study environmental constraints. Prog. Oceanogr. 83: 261-269.

Brooks S., Tyler C.R., Sumpter J.P. 1997. Egg quality in fish: What makes a good egg? Rev. Fish Biol. Fish. 7: 387-416.

Cendrero O., Cort J.L., Cárdenas E. 1981. Revisión de algunos datos sobre la biología de la anchoa, Engraulis encrasicolus (L.) del Mar Cantábrico. Bol. Inst. Esp. Oceanogr. 311:118-123.

Cort J.L., Cendrero O., Iribar X. 1976. La anchoa, Engraulis encrasicolus L. del Cantábrico. Bol. Inst. Esp. Oceanogr. 220: 3-34.

Cort J.L., Cendrero O., Cárdenas, E. 1979. Nuevos datos sobre la anchoa del Cantábrico. Información Pesquera. Inst. Esp. Ocean. 9: 1-9.

Cotano U., Irigoien X., Etxebeste E., Álvarez P., Zarauz L., Mader J., Ferrer L. 2008. Distribution, growth and survival of anchovy larvae (Engraulis encrasicolus, L.) in relation to hydrodynamic and trophic environment in the Bay of Biscay. J. Plankton Res. 30: 467-481.

Craven P., Wahba G. 1979. Smoothing Noisy Data with Spline Functions. Numer. Math. 31: 377-403.

Dickey-Collas M., Nash R.D.M., Brunel T., van Damme C.J.G., Marshall C.T., Payne M.R., Corten A., Geffen A.J., Peck M.A., Hatfield E.M.C., Hintzen N.T., Enberg K., Kell L.T., Simmonds E.J. 2010. Lessons learned from stock collapse and recovery of North Sea herring: a review. ICES J. Mar. Sci. 67: 1875-1886.

Dormann C.F. 2007. Promising the future? Global change projections of species distributions. Basic Appl. Ecol. 8: 387-397.

Dunlop E.S., Baskett M.L., Heino M., Dieckmann U. 2009. Propensity of marine reserves to reduce the evolutionary effects of fishing in a migratory species. Evol. Appl. 2: 371-393. PMCid:PMC3352486

Ferrer L., Fontán A., Mader J., Chust G., González M., Valencia V., Uriarte Ad., Collins M. B. 2009. Low-salinity plumes in the oceanic region of the Basque Country. Cont. Shelf Res. 29: 970-984.

Gordoa A., Duarte C.M. 1991. Size-dependent spatial distribution of hake (Merluccius capensis and Merluccius paradoxus) in Namibian waters. Can. J. Fish. Aquat. Sci. 48: 2095-2099.

Hastie T.J., Tibshirani R.J. 1990. Generalized Additive Models. London: Chapman and Hall, 335 pp. PMCid:PMC332745

Ibaibarriaga L., Fernández C., Uriarte A., Roel B.A. 2008. A twostage biomass dynamic model for Bay of Biscay anchovy: a Bayesian approach. ICES J. Mar. Sci. 65: 191-205.

ICES. 2007. Report of the Workshop on Testing the Entrainment Hypothesis (WKTEST), 4-7 June 2007, Nantes, France. ICES C.M. 2007/LRC:10, 111 pp.

ICES. 2008a. Report of the Working Group on Anchovy (WGANC), 13 - 16 June 2008, ICES Headquarters, Copenhagen. ICES C.M. 2008 ACOM:04, 226 pp.

ICES. 2008b. Report of the Working Group on Acoustic and Egg Surveys for Sardine and Anchovy in ICES Areas VIII and IX (WGACEGG), 24–28 November 2008, Nantes, France. ICES C.M. 2008 LRC:17, 183 pp.

Irigoien X., Fiksen Ø., Cotano U., Uriarte A., Alvarez P., Arrizabalaga H., Boyra G., Santos M., Sagarminaga Y., Otheguy P., Etxebeste E., Zarauz L., Artetxe I., Motos L. 2007. Could Biscay Bay Anchovy recruit through a spatial loophole? Prog. Oceanogr. 74: 132-148.

Irigoien X, Cotano U., Boyra G., Santos M. Alvarez P., Otheguy P., Etxebeste E., Uriarte A., Ferrer L., Ibaibarriaga L. 2008. From egg to juvenile in the Bay of Biscay: spatial patterns of anchovy (Engraulis encrasicolus) recruitment in a non-upwelling region. Fish. Oceanogr. 17: 446-462.

Irigoien X., Chust G., Fernandes J.A., Albaina A., Zarauz L. 2011. Factors determining the distribution and betadiversity of mesozooplankton species in shelf and coastal waters of the Bay of Biscay. J. Plankton Res. 33: 1182-1192.

Jørgensen C., Dunlop E.S., Opdal A.F., Fiksen Ø. 2008. The evolution of spawning migrations: state dependence and fishinginduced changes. Ecology 89: 3436-3448. PMid:19137949

Kennedy J., Geffen A, Nash, R. 2007. Maternal influences on egg and larval characteristics of plaice (Pleuronectes platessa L.). J. Sea Res. 58, 65-77.

Lasker R. (ed.). 1985 An egg production method for estimating spawning biomass of pelagic fish: Application to the northern anchovy, Engraulis mordax. US Department of Commerce, NOAA Technical Report NMFS 36, 99 pp.

Lehuta S., Mahévas S., Petitgas P., Pelletier D. 2010. Combining sensitivity and uncertainty analysis to evaluate the impact of management measures with ISIS–Fish: marine protected areas for the Bay of Biscay anchovy (Engraulis encrasicolus) fishery. ICES J. Mar. Sci. 67: 1063-1075.

Loots C., Vaz S., Planque B., Koubbi P. 2011. Understanding what controls the spawning distribution of North Sea whiting (Merlangius merlangus) using a multi-model approach. Fish. Oceanogr. 20: 18-31.

Marteinsdottir G., Steinarsson A. 1998. Maternal influence on the size and viability of Icelandic cod (Gadus morhua L.) eggs and larvae. J. Fish Biol. 52: 1241-1258.

Marteinsdottir G., Thorarinsson K. 1998. Improving the stockrecruitment relationship in Icelandic cod (Gadus morhua) by including age diversity of spawners. Can. J. Fish. Aquat. Sci. 55: 1372-1377.

Morgan M.J., Murua H., Kraus G., Lambert Y., Marteinsdóttir G., Marshall C.T., O'Brien L., Tomkiewicz, J. 2009. The evaluation of reference points and stock productivity in the context of alternative indices of stock reproductive potential. Can. J. Fish. Aquat. Sci. 66: 404-414.

Motos L., Uriarte A., Valencia V. 1996. The spawning environment of the Bay of Biscay (Engraulis encrasicolus L.). Sci. Mar. 60(Suppl. 2): 117-140.

Mullon C., Cury. P., Penven P. 2002. Evolutionary individual-based model for the recruitment of anchovy (Engraulis capensis) in the southern Benguela. Can. J. Fish. Aquat. Sci. 59: 910-922.

Parker K. 1980. A direct method for estimating northern anchovy, Engraulis mordax, spawning biomass. Fish. B-NOAA 78: 541–544.

Petitgas P., Massé J., Grellier P., Beillois P. 2003. Variation in the spatial distribution of fish length: a multi annual geostatistics approach on anchovy in Biscay 1983-2002. ICES C.M. 2003/Q:15, 14 pp.

Petitgas P., Secor D. H., McQuinn I., Huse G., Lo N. 2010. Stock collapses and their recovery: mechanisms that establish and maintain life-cycle closure in space and time. ICES J. Mar. Sci. 69: 1841-1848.

Planque B., Bellier E., Lazure P. 2007. Modelling potential spawning habitat of sardine (Sardina pilchardus) and anchovy (Engraulis encrasicolus) in the Bay of Biscay. Fish. Oceanogr. 16: 16-30.

Planque B., Loots C., Petitgas P., Lindstrøm U., Vaz S. 2011a. Understanding what controls the spatial distribution of fish populations using a multi-model approach. Fish. Oceanogr. 20: 1-17.

Planque B., Bellier E., Loots C. 2011b. Uncertainties in projecting spatial distributions of marine populations. ICES J. Mar. Sci. 68: 1045-1050.

R Development Core Team. 2011. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0.

Santiago J., Sanz A. 1992. Egg production estimates of the Bay of Biscay anchovy, Engraulis encrasicholus (L) spawning stock in 1987 and 1988. Bol. Inst. Esp. Oceanogr. 8: 225-230.

Santos M., Uriarte A., Ibaibarriaga L. 2011. Spawning stock biomass estimates of the Bay of Biscay anchovy (Engraulis encrasicolus, L.) in 2010 applying the Daily Egg Production Method. RIM 18: 76-91.

Shchepetkin A.F., McWilliams J.C. 2005. The Regional Ocean Modeling System: A split-explicit, free-surface, topography following coordinates ocean model, Ocean Model. 9: 347-404.

Simmonds J.E., MacLennan D.N. 2005. Fisheries acoustics: Theory and practice, 2nd edition. Oxford: Wiley-Blackwell, 456 pp.

Smith W.H.F., Sandwell D.T. 1997. Global seafloor topography from satellite altimetry and ship depth soundings. Science 277: 1957-1962.

Somarakis S., Palomera I., Garcia A., Quintanilla L., Koutsikopoulos C., Uriarte A., Motos L. 2004. Daily egg production of anchovy in European waters. ICES J. Mar. Sci. 61: 944-958.

Stratoudakis Y., Bernal M., Borchers D., Borges, M.F. 2003. Changes in the distribution of sardine eggs and larvae off Portugal, 1985-2000. Fish. Oceanogr. 12: 49-60.

Stratoudakis Y., Bernal M., Ganias K., Uriarte A. 2006. The daily egg production method: recent advances, current applications and future challenges. Fish Fish. 7: 35-57.

Thorrold S.R., Latkoczy C., Swart P.K., Jones C.M. 2001. Natal homing in a marine fish metapopulation. Science 291: 297-299. PMid:11209078

Uriarte A. 2001. Anchovy population at age estimates and variances from the application of the daily egg production method. ICES C.M. 2001/P:25, 6 pp.

Uriarte A., Prouzet P., Villamor B. 1996. Bay of Biscay and Ibero Atlantic anchovy population and their fisheries. Sci. Mar. 60(Suppl. 2): 237-255.

Vaz S., Petitgas P. 2002. Study of the Bay of Biscay anchovy population dynamics using spatialised age-specific matrix models. ICES C.M. 2002/O:07, 20 pp.

Werner E.E., Gilliam J.F. 1984. The ontogenetic niche and species interactions in size-structured populations. Ann. Rev. Ecol. Syst. 15: 393-425.

Wood S.N. 2000. Modelling and smoothing parameter estimation with multiple quadratic penalties. J. R. Stat. Soc. Ser. B Stat. Methodol. 62: 413-428.

Wood S.N. 2004. Stable and efficient multiple smoothing parameter estimation for generalized additive models. J. Amer. Statist. Assoc. 99: 673-686.

Wood S.N. 2006. Generalized Additive Models: An introduction with R. Chapman and Hall/CRC, Boca Ratón, 392 pp.

Wood S.N., Augustin N.H. 2002. GAMs with integrated model selection using penalized regression splines and applications to environmental modelling. Ecol. Model. 157: 157-177.




How to Cite

Ibaibarriaga L, Uriarte A, Laconcha U, Bernal M, Santos M, Chifflet M, Irigoien X. Modelling the spatio-temporal distribution of age-1 Bay of Biscay anchovy (Engraulis encrasicolus) at spawning time. Sci. mar. [Internet]. 2013Sep.30 [cited 2024Apr.21];77(3):461-72. Available from:




Most read articles by the same author(s)