Search

For over 50 years, the conceptualisation of low-nutrient oligotrophic systems having longer food chains and thus lower energy transfer to fish than their high-nutrient eutrophic counterparts1 has achieved the status of an ecological paradigm. However, recent global assessments indicate global fish biomass could be much higher than previously thought2–4, suggesting that our traditional understanding of food webs may need to be revisited. Here, we challenge the classical paradigm by exploring the role of zooplankton in food webs across the world’s oceans. Using observed zooplankton size spectra, and output from a size-spectrum model that resolves nine zooplankton groups, we conclude that food chains in oligotrophic (low-nutrient) and eutrophic (high-nutrient) systems have similar lengths. We offer a compelling hypothesis to explain this emergent pattern: self-organisation of zooplankton groups across the global productivity gradient regulates food chain length. We find that in oligotrophic systems the increased carnivory and longer food chains are offset by relatively large gelatinous filter feeders eating the dominant small phytoplankton, resulting in shorter-than-expected food chains, but decreasing food quality for fish. Our findings highlight the pivotal role zooplankton play in regulating energy transfer. Better resolution of zooplankton groups, their feeding relationships and carbon content in models will increase our ability to estimate current global fish biomass 5, project future fish biomass under climate change6–8, and provide more-robust forecasts of nutrient9 and carbon cycling10.

Since the 1980s, the Common Fisheries Policy (CFP) has shaped European fisheries. It has often been criticised for being too prescriptive and, above all, for failing to protect either fishermen or ecosystems. The last reform dates back to the early 2010 s and has led to a slight but slow improvement in the state of ecosystems. Given that the CFP is in the process of evaluation, a group of French fishery scientists set up an initiative to add to the debate on what should be retained, reinforced or added to a possible new reform. This initiative came 10 years after a previous manifesto that presented their vision for fisheries in Europe. Four major issues emerged from the current initiative: (1) a need for transparency and simplification in fisheries management, (2) a need for more consultation and dialogue between stakeholders, (3) the urgency of the situation in the Mediterranean Sea, and (4) the necessity of putting into practice all research developments for an ecosystem approach to fisheries. Compared to 10 years ago, the response of scientists shows that the focus is no longer on achieving the maximum sustainable yield, but rather on the following steps to protect ecosystems and fisheries. An ecosystem approach to fisheries remains indispensable for both ecosystems and fishing activities. To this end, scientists put forward numerous proposals to improve the CFP, acknowledging that the final solutions should emerge from consultation with stakeholders. Climate change, an issue raised much more than in the manifest, reinforces the need to act.

Anchovy Engraulis encrasicolus distribution in European waters spans from the Mediterranean Sea to the North Sea, and is expected to expand further north with global warming. Observations from the eastern Mediterranean (North Aegean Sea), the Bay of Biscay and the North Sea reveal latitudinal differences in growth, maximum size, fecundity and timing of reproduction. We set up a mechanistic framework combining a bioenergetics model with regional physical-biogeochemical models providing temperature and zooplankton biomass to investigate the underlying mechanisms of variation in these traits. The bioenergetics model, based on the Dynamic Energy Budget theory and initially calibrated in the Bay of Biscay, was used to simulate growth and reproduction patterns. Environment partly explained the increased growth rate and larger body size towards the north. However, regional calibration of the maximum assimilation rate was necessary to obtain the best model fit. This suggests a genetic adaptation, with a pattern of cogradient variation with increasing resource towards the north, in addition to a countergradient thermal adaptation. Overall, the seasonal energy dynamics supports the pattern of body-size scaling with latitude, i.e. food-limited growth but low maintenance costs in the warm Aegean Sea, and larger size in the North Sea allowing sufficient storage capacity for overwintering. Further, the model suggests a synchronisation of reproductive timing with environmental seasonality as a trade-off between thresholds of temperature and reserves for spawning and overwintering, respectively. Finally, low temperature, short productive and spawning seasons, and insufficient reserves for overwintering appear to be current limitations for an expansion of anchovy to the Norwegian Sea.

Understanding discrepancies in particular years between acoustic and Daily Egg Production Method (DEPM) estimates is essential to better understand potential bias impacting both estimates, and develop a weighting per year assessing the coherence between survey estimates, to input in stock assessment models. Several plans for collaborative work have been proposed within WGACEGG to compare acoustic and DEPM indices (WGACEGG 2013, 2015, 2016), until the inclusion of such comparative approach as a formal group ToR in 2017 (ICES, 2017). Different approaches have been adopted so far, including the comparison between two egg production estimates (DEPM BIOMAN; CUFES PELGAS) and one acoustic estimate (PELGAS), analysis of the sensitivity of each estimate to its parameters, and the comparison between egg and adult maps. We here provide some rationale on how to investigate differences between acoustic and egg indices.

A prerequisite in fisheries management is to identify biological meaningful delineations of stocks, the fundamental exploited units. The European anchovy is one such species in need to better identify stock boundaries. Indeed, despite a spawning aggregation behaviour both in the southern Bay of Biscay and southern North Sea in spring and summer, it has a seemingly continuous distribution in autumn across a broad region ranging from the Bay of Biscay to the northern North Sea including the English Channel, thus with potential mixing of the spawning units. We therefore used genetic markers (single nucleotide polymorphisms) to quantify the degree of gene flow between the currently managed fish stock of the bay of Biscay and northern areas of European anchovy’s distribution. Our results confirm the clear distinction between the Bay of Biscay and northernmost populations, with assignment of all English Channel samples to the latter. We also found a clear overall pattern of isolation by distance that resulted primarily from an increasing differentiation with geographic distance in the Northern group magnified by the lack of gene flow with the anchovies of the Bay of Biscay. The small reduction in gene diversity towards northern latitudes may indicate partial isolation of the leading-edge component of the Northern population. Quotient plots relating anchovy’s distribution to environmental covariates showed that anchovies sampled in autumn in the English Channel originate from the summer spawning aggregation in the warm, low saline and plankton rich south-eastern North Sea. This change in the seasonal distribution of the northern population mimics, but in opposite direction, the one of the Bay of Biscay where anchovy spread towards the north from the spawning habitat in the south-eastern bay. The encounter, without mixing, of the two populations west of Brittany in autumn suggests strong spawning site fidelity. Finally, we identified for the first time anchovies belonging to the estuarine ecotype in the Loire estuary. Our results therefore support current management units, yet with some uncertainties for the catches in the transition zones between the Bay of Biscay and the English Channel, and within river plumes in the vicinity of estuaries hosting the estuarine ecotype.

This paper presents a novel application of the geostatistical multivariate method known as min–max autocorrelation factors (MAFs) for analysing fisheries survey data in a space–time context. The method was used to map essential fish habitats and evaluate the variability in time of their occupancy. Research surveys at sea on marine fish stocks have been undertaken for several decades now. The data are time series of yearly maps of fish density, making it possible to analyse the space–time variability in fish spatial distributions. Space–time models are key to addressing conservation issues requiring the characterization of variability in habitat maps over time. Here, the variability in fisheries survey data series is decomposed in space and time to address these issues, using MAFs. MAFs were originally developed for noise removal in hyperspectral multivariate data and are obtained using a specific double principal components analysis. Here, MAFs were used to extract the most continuous spatial components that are consistent in time, together with the time series of their amplitudes. MAFs formed an empirical isofactorial model of the data, which served for kriging in each year using all available information across the data series. The approach was applied on the spawning distributions of sardine in the Bay of Biscay from 2000 to 2017. A multivariate approach for dealing with space–time data was adapted here, because the evolution in time was highly variable. Maps were classified using the amplitudes of the MAFs, and groups of typical distributions were identified, which showed different occurrence probabilities in different periods.

This paper is a practical contribution to two important debates raised by the implementation of marine ecosystem based management: (i) which ecosystem data can be provided by a fisheries survey optimised for ecosystem monitoring; and (ii) how to combine/select potential indicators to derive useful information on marine ecosystem status and dynamics? A suite of 143 potential indicators, including spatial indices, representing small pelagic fish and their biotic and abiotic environments are presented. Indicators were routinely derived from the PELGAS integrated ecosystem survey conducted in spring in the Bay of Biscay (BoB). The general patterns over time in this suite of 5–16 years, non-stationary time series are characterised using a methodology based on min–max autocorrelation factors (MAF), to select the most continuous indicators within, and across, several ecosystem components: hydrology, phytoplankton, mesozooplankton, small pelagic fish and megafauna. Potential interactions between selected indicators and external forcing variables, including climate and fishing, were assessed. The results confirm the importance of river discharges, bottom temperature, chlorophyll-a and mesozooplankton biomass in the dynamics of the BoB pelagic ecosystem. Small pelagic fish species appear to have followed distinct trajectories during the last 15 years. A marked decrease in anchovy and sardine mean weights at ages 1 and 2 over the last 15 years was highlighted; potentially caused by density-dependent competition. The quasi-absence of significant correlation between selected survey indicators and climate and fishing pressure proxies suggests a moderate exploitation rate of small pelagic fish resources, and confirms the so far limited effects of large-scale climate forcing on the BoB pelagic ecosystem. Perspectives for the assessment of marine ecosystem status on the basis of suites of indicators derived from integrated ecosystem surveys are discussed.

Similar or very contrasted puzzling population dynamics between anchovy and sardine occur worldwide. Underlying factors are not well understood, but insights towards different biological traits are suggested, in particular trophic specialisation, leading to different responses to environmental conditions. Based on most striking differences in biological and life history traits, i.e. size, spawning and feeding, we calibrated a bioenergetics model, based on the Dynamic Energy Budget theory, for Engraulis encrasicolus and Sardina pilchardus in the Bay of Biscay. Starting from the anchovy model, differences in traits were successively integrated to build the sardine model through a novel exploratory approach by scenarios. We used a robust method for parameter estimation, the Evolution Strategies, with a large dataset of length and mass at age, as well as energy density, which is the first time in such a model calibration. Energy density data proved to be particularly well suited to assess the quality of DEB model predictions and parameter set estimates. Insights in respective physiology were drawn from analysis of parameter values and predictions of the model. We showed that anchovy and sardine have distinct strategies with respect to energy acquisition and especially to allocation to spawning. Anchovy are characterised by higher metabolic rates and requirements. This species is more likely to benefit from periods of high food availability to carry out both growth, spawning and reserve storage. Sardine have less demanding food requirements and metabolic costs. Sardine take advantage of larger reserves storage capacity to decouple spawning and prey blooms and to lengthen spawning period, and thus display a more capital breeding spawning behaviour. Overall, our model outputs distinguish between anchovy that tend towards an almost “all or nothing” energetic strategy, and sardine that tend to carry out lower metabolic activities but on a more regular basis. This first modelling demonstration of a bioenergetics difference between these two species, and the explanation it brings in the understanding of their respective reproduction strategies, opens new perspectives in the interpretation of their differential responses at the population scale to environment variability.

The population of European anchovy of the Bay of Biscay collapsed at the beginning of the 21st century, causing the closure of its fishery between 2005 and 2010. In order to study both the human and environmental causes of the anchovy population dynamics, an approach coupling individual bioenergetics to an individual-based model was applied between 2000 and 2015. This modelling framework was forced with outputs from a physical–biogeochemical model. In addition to a base-case scenario with realistic forcing, alternative scenarios were run without inter-annual variability in either fishing mortality or environmental conditions. During the decrease in population biomass, a high fishing pressure coincided with a combination of environmental variables promoting the appearance of large individuals that could not survive severe winters because of their high energetic demands. The recovery of the population was favoured by a period of warm years with abundant food favouring the winter survival of age 1 individuals, in coincidence with the closure of the fishery. Our modelling approach also allows to test the consequences of a retrospective implementation of the current harvest control rule from 2000 which, according to our results, would have prevented the collapse of the population and avoided the fishery closure.

From 2000 to 2016, substantial changes in biomass and community structure of small pelagic fish and mesozooplankton have been reported in the Bay of Biscay. Since significant relationships have been found between phytoplankton chlorophyll a and mesozooplankton as well as between phytoplankton chlorophyll a and shifts in sardine body condition, it was hypothesized that phytoplankton communities may have also been affected during this period and may have played a role in these changes. However, the available data were insufficient to validate this hypothesis and the causes of these changes remained unexplained. The present study analyzed a spatio-temporal marine microphytoplankton dataset collected during the annual PELGAS (PELagiques GAScogne) surveys from 2003 to 2014. The thorough analysis of microphytoplankton taxonomic composition, with an approach integrating the relative role of environmental conditions as well as biotic interactions and applying the concept of ecological niche, confirmed that significant modifications in microphytoplankton community structure occurred during this period. Temporal changes were stronger than spatial differences at these sampling scales. Three main periods, 2003–2005, 2006 and 2007–2014, showing different community structure, diversity and dominant taxonomic units were highlighted. Twenty eight taxonomic units were involved in these community changes. Among them, five were identified as the protagonists (Pseudo-nitzschia spp., Gymnodinium spp. + Gyrodinium spp., Leptocylindrus danicus, Leptocylindrus minimus and Chaetoceros sp.). Variations in water temperature and equivalent freshwater depth constrained the realized ecological niches of these species and explained, at least in part, changes in community structure. This study stresses the need to improve our knowledge of phytoplankton species ecological niches and to take into account biotic interactions for a thorough understanding of the processes shaping plankton communities and the resulting diversity patterns.

The size and specific gravity of eggs of marine pelagic fish partly determine their dispersal and survival. Using an original dataset of anchovy and sardine eggs, sampled in spring over the last decade in the Bay of Biscay, we provide a parameterization of these properties on ambient water temperature and salinity. We used the density gradient column for measurement of egg specific gravity. The column was also filled with homogeneous water for sinking velocity experiments. For anchovy, these experiments confirm that the effect of egg permeability through the chorion could be neglected when modelling sinking, while it has to be considered for sardine, its perivitelline space representing 78.6% (±6.2%) of the total egg volume, as opposed to 5–10% for most teleosts species. We estimated a coefficient of permeability of the chorion of 0.0038 mm s−1. However, permeability should not affect the measurement of sardine egg specific gravity in a gradient column, provided a minimum duration before reading is respected for equilibrium to be reached. In relation to their environment, we found that the egg specific gravity is largely determined by sea surface salinity for both species, whereas egg size is weakly but significantly impacted by temperature, for sardine only. On average, the estimated difference in specific gravity between egg and surface water is −0.92 σT for anchovy and −1.06 σT for sardine. The detailed parameterization of the relationship between eggs and water properties should prove useful, in particular to modellers dealing with the dispersal of fish early life stages.

In the Bay of Biscay, mean body length and weight of sardines (Sardina pilchardus) have been decreasing since the early 2000s and could severely impact the fishing and seafood industry sector. These trends have no apparent link with fishing pressure, although the latter has been increasing since the late 2000s. As part of an effort to develop suitable assessment and management tools for this stock, we investigated the life-history traits of sardine and analyze its seasonal and inter-annual variations. Based on 14 years of morphometric data from both scientific surveys and professional samples, we analyzed the variability in sardine body condition and its responses to environmental changes. Generalized Additive Models revealed an age-sex specific decreasing trend in body length over the study period, with most of the variability explained by the age class. Linear Mixed Effect Models applied to the body condition evidenced its strong seasonality and an age class specific decreasing trend. Regardless of age class, maximal body condition is reached at the end summer, after the spawning and plankton productive periods. Overall, annual trends in body condition-at-age showed remarkable coherence, with a significant decrease since 2007 for all age classes, suggesting that factors influencing body condition operate at population level. The shift in sardine body condition towards lower values could be broken down into three periods, with a high dependence on surface Chlorophyll-a and sea surface temperature. However, this study highlights that the period supporting the main decrease in body condition is characterized by high Chlorophyll-a, the available proxy for food, which is counterintuitive. Such a result suggests more complex trophic responses involving secondary production, with potential shift in the timing of the production and/or the quality of the food. At the population level, those changes may have a long-term negative effect, with a decrease in body length and important changes in phenology (length at first maturity, reproductive phenology) and potential consequences on sardine population dynamics in the Bay of Biscay.

In this paper we present the ZooCAM, a system designed to digitize and analyse on board large volume samples of preserved and living metazooplankton (i.e. multicellular zooplankton) and fish eggs >300 µm ESD. The ZooCAM has been specifically designed to overcome the difficulties to analyse zooplankton and fish eggs in the framework of the PELGAS survey, and provide high frequency data. The ZooCAM fish eggs counts were comparable to those done with a dissecting microscope. The ZooCAM enabled the accurate prediction and fast on board validation of staged anchovy and sardine eggs in almost real time after collection. A comparison with the ZooScan, on a more complex zooplanktonic community, provided encouraging results on the agreement between the 2 instruments. ZooCAM and ZooScan enabled the identification of similar communities and produced similar total zooplankton abundances, size distributions, and size spectra slopes, when tested on the same samples. However these results need to be further refined due to the small number of samples used to compare the two instruments. The main ZooCAM drawback resides in a slight but sensible underestimation of abundances and sizes, and therefore individual and community biovolumes. The ZooCAM have been successfully deployed over 4 years, on numerous surveys without suffering any major failure. When used in line with the CUFES it provided high resolution maps of staged fish eggs and zooplanktonic functional groups. Hence the ZooCAM is an appropriate tool for the development of on board, high frequency, high spatial coverage zooplanktonic and ecosystemic studies.

Fishery annual surveys are shifting towards integrated ecosystem surveys. However, monitoring the whole ecosystem once a year entails special attention to the interpretation of observations, because they only represent a snapshot in the seasonal dynamics of the environment, of a species, or any dynamic process of interest. The dependency of this snapshot to the timing, but also to the duration, of the survey with respect to this seasonal dynamics needs to be considered. Fish stock assessment is only little impacted by the observed timing within a given season, whereas warming, stratification, plankton bloom occurrence, or fish fecundity, are processes rapidly changing especially in spring. Firstly, from independent satellite data, we described the seasonal and interannual Sea Surface Temperature (SST) variability in the Bay of Biscay in spring around the annual PELGAS surveys. Our results revealed that the in-situ surface temperature snapshot from a given survey may be strongly misleading on the existing spatial pattern, blurring the latitudinal temperature gradient, with potential impact in e.g. species habitat modelling. Secondly, based on these survey-independent SST data, we proposed a methodology to position annual surveys in the environment spring schedule. Our temporal indicator is based on a degree-day metric, computed as the cumulated temperature from the previous winter date with minimum annual temperature. The annual degree-day is then compared to the climatology of the degree-day metric to correct the actual date of survey to an environment date. The methodology was spatialised to take into account both latitudinal shift in phasing of seasonality, and rapidly changing conditions in the environment during the survey itself. Finally, we tested the methodology to refine the description of anchovy and sardine spawning phenology, as well as spawning spatial distribution for particular years such as hot year 2003. More generally, this is an example of how synoptic, independent data from operational oceanography, as a background history of the environment dynamics, can be used to improve the interpretation of biological observations in marine ecosystems

International audience; Mesozooplankton can be considered the most important secondary producers in marine food webs because they hold an intermediate position between the phytoplankton assemblage and the upper trophic levels. They also are a robust indicator of climatic and hydrological conditions. We conducted an analysis of the interannual variability of the spring mesozooplankton assemblage, as sampled by the PELGAS fisheries survey in the southern part of the Bay of Biscay (Northeast Atlantic Ocean) between 2003 and 2013. We examined hydrology and trophic drivers to explain the variability. Our results revealed that the subsurface temperature, the subsurface salinity, the biomasses of subsurface pico-, nano-, and microphytoplankton, and the copepod assemblage exhibited a recurrent spatial pattern that was driven mainly by freshwater and nutrient inputs from the main rivers. The mesozooplankton assemblage was dominated by copepods (82%), composed of coastal, neritic, and oceanic copepod genera that paralleled the various hydrological fronts converging in the southern Bay of Biscay. The copepod community displayed high temporal-variability; there were three periods of abundant adult copepods throughout the southern Bay of Biscay. The copepod community was structured primarily around the drive for resource control, especially by the microphytoplankton biomass (24.3% of the total variability), and to a lesser extent by hydrological features (13.7% of the total variability).

International audience; Coccolithophore blooms occur regularly from April to June in the Bay of Biscay where they have been observed for many years from ocean-colour imagery thanks to the ability of their calcite plates to scatter light. They are easily depicted on interpolated images of non-algal Suspended Particulate Matter (SPM) derived from satellite reflectance data in May, at the time of the PELGAS (PELagique GAScogne) annual surveys. Over the springs 2012 to 2015, the PELGAS surveys provided in-situ data on the turbidity throughout the water column and on the hydrological environment prevailing during the blooms in the vicinity of the continental shelf break. The satellite-derived SPM in the area of coccolithophore blooms was closely related to measured turbidity in the euphotic zone. The deep maxima of coccolithophore-derived biomass were well observed in the turbidity signal throughout the water column. Despite the recurrence of SPM patterns due to coccolithophores in the well-mixed waters over the shelf break in spring, the variability of the satellite-derived SPM and the turbidity profiles was particularly high during the four years studied. The year 2013 showed persistent blooms from mid-April to the end of May over large areas inside and outside the shelf waters, whereas no bloom was observed in 2014, neither on SPM images nor through turbidity profiles. The presence or absence of coccolithophore blooms during the years 2012-2015 seem related to the level of stratification of the water column, with more intense blooms occurring during mixed conditions. These results should improve classification of the PELGAS surveys from those occurring in a coccolith-rich environment, synonymous with mixed and cold waters, to those occurring in a coccolith-poor environment where surface waters are warmer and more stratified.

International audience; The Pélagiques Gascogne (PELGAS) integrated survey has been developed by a multidisciplinary team of Ifremer and La Rochelle University scientists since 2000, joined by commercial fishermen in 2007. Its initial focus was to assess the biomass and predict the recruitment success of anchovy in the Bay of Biscay in spring. Taking advantage of the space and versatility of R/V Thalassa II, sampling has been progressively extended to other ecosystem components. PELGAS therefore further developed the second objective of monitoring and studying the dynamic and diverse Biscay pelagic ecosystem in springtime. The PELGAS survey model has allowed for the establishment of a long-term time-series of spatially-explicit data of the Bay of Biscay pelagic ecosystem since the year 2000. Main sampled components of the targeted ecosystem are: hydrology, phytoplankton, mesozooplankton, fish and megafauna. The survey now provides two main ecosystem products: standard raster maps of ecosystem parameters, and a time series dataset of indicators of the Bay of Biscay pelagic ecosystem state. They are used to inform fish stock and ecosystem-based management, and support ecosystem research. The present paper introduces the PELGAS survey, as a practical example of an integrated, vessel-based, ecosystem survey. The evolution of the PELGAS scientific team and sampling protocols are presented and analysed, to outline factors crucial to the success of the survey. Data and results derived from PELGAS are reviewed, to exemplify scientific questions that can be tackled by integrated ecosystem survey data. Advantages and challenges of the survey are discussed and put into the context of marine ecosystem surveys in the European Marine Strategy Framework Directive and the Common Fisheries Policy.

Mapping habitats of species communities helps to inform on the ecological processes which drive their distribution. This information is critical to identify suitable areas for spatial management, aimed at preserving biodiversity, ecosystem functions or essential habitats. While demersal fish have been extensively studied at the community scale, small pelagic fish have mainly been characterised at the population scale. This paper presents a community-based approach on the biodiversity of small pelagic fish, with the aim to: (i) define small pelagic fish communities, (ii) characterise their spatial and interannual dynamics, and (iii) assess their habitats. We present a Multiple Factor Analysis (MFA)-based method that characterises the spatio-temporal variability in a series of multivariate maps. The main components of variability in these maps are identified to define “characteristic areas” displaying contrasting conditions in both MFA ordination and geographical spaces. Characteristic habitats of fish communities were defined by the correlation between the main components of variability from MFAs applied to maps containing fish information on one hand, and on hydrographic variables on the other hand. The analysis of the long term PELGAS survey series of fish biomass and hydrological indices maps resulted in the characterisation of mesoscale latitudinal gradients and coarse to mesoscale onshore-offshore gradients in both fish and hydrology datasets. A community with anchovy (E. encrasicolus) and chub mackerel (S. colias) as indicator species was consistently distributed in southeast Biscay. This area was associated with higher bottom temperatures, which likely affected the fish community through physiological processes. A second community with small clupeiforms as indicator species was found in coastal spawning habitats. These habitats were typically characterised by low salinity, a probable proxy for high productivity and good feeding grounds for those species. Mapping the habitats of small pelagic fish communities may inform marine spatial management, aimed at preserving biodiversity, ecosystem structure and function. In addition, it may help in achieving maximum sustainable yields of these commercially important species, and contribute to achieving and maintaining good environmental status of shelf seas ecosystems.

The present study describes the responses of summer phytoplankton biomass to changes in top-down forcing (expressed as zooplankton mortality) in three ecosystems (the North Sea, the Baltic Sea and the Nordic Seas) across different 3D ecosystem models. In each of the model set-ups, we applied the same changes in the magnitude of mortality (±20%) of the highest trophic zooplankton level (Z1). Model results showed overall dampened responses of phytoplankton relative to Z1 biomass. Phytoplankton responses varied depending on the food web structure and trophic coupling represented in the models. Hence, a priori model assumptions were found to influence cascades and pathways in model estimates and, thus, become highly relevant when examining ecosystem pressures such as fishing and climate change. Especially, the different roles and parameterizations of additional zooplankton groups grazed by Z1, and their importance for the outcome, emphasized the need for better calibration data. Spatial variability was high within each model indicating that physics (hydrodynamics and temperature) and nutrient dynamics also play vital roles for ecosystem responses to top-down effects. In conclusion, the model comparison indicated that changes in top-down forcing in combination with the modelled food-web structure affect summer phytoplankton biomass and, thereby, indirectly influence water quality of the systems.

Integrated marine survey series such as PELGAS provide comprehensive geo-referenced data over large sea areas in major ecosystem components with defined biological resolution and spatio-temporal sampling scale. These data enable to develop product at different levels of biological organization and spatial scale that are useful for ecosystem integrated assessments. Using the PELGAS integrated data series, we applied a generic procedure made of several steps to identify and map ecologically coherent ecosystem spatial units in the Bay of Biscay. First, the data were interpolated on a common spatial grid and organized as a time series of matrices containing at each time the variables as columns and the grid cell values as rows. The multi-table analysis method known as Multi-Factor Analysis (MFA) was then applied on the series of matrices, thus separating in the analysis the spatial from the temporal variability. Grouping the spatial grid cells based on their relative positions in the MFA space resulted in mapping ecosystem spatial entities based on common spatial patterns across ecosystem components. The result of the analysis is a map of ecosystem seascapes that are consistent over the years together with a map of their inter-annual variability. The ecosystem sub-units were in agreement with sub-regional production systems in the Bay of Biscay. This study thus highlights the possibility to characterize and monitor ecosystem spatial structure and develop indicators thereof for their use in ecosystem assessments. Also, it highlights the importance of identifying spatial limits of production systems for ecosystem description, assessment and management.

Peu de modèles d'évaluation de stock intègrent une description fine des processus biologiques régissant les populations exploitées par la pêche au sein des écosystèmes. Ce manque est particulièrement regrettable dans le cas des petits pélagiques dont la dynamique est fortement influencée par les conditions environnementales. En effet, la pêche et les fluctuations environnementales et climatiques impactent fortement la dynamique de population de la sardine du golfe de Gascogne. Pour comprendre la dynamique de cette population au sein de l'écosystème et les effets des pêcheries, il est donc nécessaire d'identifier et de caractériser cet ensemble de pressions. Afin d'améliorer la gestion et le suivi de la dynamique de ce stock, le développement d'un modèle d'évaluation écosystémique a été initié. Sa formalisation nécessite notamment la prise en compte de facteurs endogènes et exogènes structurant la croissance et la condition de cette espèce. Les indices de condition, notamment l'indice de condition relative de Le Cren, permettent notamment de comparer l'état physiologique individuel des poissons et donc par extrapolation de caractériser l'état de la population. Les études précédentes chez les sardines (mais également chez d'autres espèces de poissons) montrent qu'il existe d'importantes variations spatiotemporelles dans la condition, la croissance et la structure en taille et en âge. En nous basant sur les données récoltés lors de la campagne PELGAS, nous mesurons dans un premier temps les variations de ces paramètres clés dans le temps et dans l'espace (gradient latitudinal) afin d'identifier les facteurs environnementaux régissant l'expression de ces caractères morphologiques qui sont liés à la croissance et la reproduction des individus. Dans un second temps, pour aborder ces mécanismes d'un point de vue quantitatif, une approche bioénergétique individuelle à l'aide du DEB (Dynamic Energy Budget) sera développée pour analyser les processus bottom-up. Son utilisation en IBM (Individual-based model) permettra d'explorer également les processus de densitédépendance.

An acoustic survey (PELGAS) is carried out every year in the Bay of Biscay in spring onboard the French research vessel Thalassa. The objective of PELGAS survey is to study the abundance and distribution of pelagic fish in the Bay of Biscay. The main target species are anchovy and sardine but they are considered in a multi-specific context and within an ecosystemic approach as they are located in the centre of pelagic ecosystem. The Pelgas17 acoustic survey has been carried out with medium weather conditions (regular wind, low atmospheric temperatures) for the whole area, from the South of the bay of Biscay to the west of Brittany. The help of commercial vessels (two pairs of pelagic trawlers and a single one) during 18 days provided about 110 valid identification hauls instead of about 60 before 2007 when Thalassa was alone to identify echotraces. Their participation increased the precision of identification of echoes and some double hauls permitted to confirm that results provided by the two types of vessels (R/V and Fishing boats) were comparable and usable for biomass estimate purposes. These commercial vessels participated to the PELGAS survey in a very good spirit of collaboration. Vessels (and the scientific observer onboard) are founded by EMFF (European Maritime and Ficheries Found) for the period 2017- 2019, with the financial help of "France Filière Pêche" which is a groupment of French fishing organisations. The PELGAS17 survey observed a relatively high level of anchovy biomass (134 500 tons), which seems to be higher to previous year, comparable to 2011 and far away from the 2015 biomass The biomass estimate of sardine observed is 465 022 tons, which constitutes a strong increase from the last survey.

Hydrodynamic Ocean Circulation Models and Lagrangian particle tracking models are valuable tools e.g. in coastal ecology to identify the connectivity between offshore spawning and coastal nursery areas of commercially important fish, for risk assessment and more for defining or evaluating marine protected areas. Most studies are based on only one model and do not provide levels of uncertainty. Here this uncertainty was addressed by applying a suite of 11 North Sea models to test what variability can be expected concerning connectivity. Different notional test cases were calculated related to three important and well-studied North Sea fish species: herring (Clupea harengus), and the flatfishes sole (Solea solea) and plaice (Pleuronectes platessa). For sole and plaice we determined which fraction of particles released in the respective spawning areas would reach a coastal marine protected area. For herring we determined the fraction located in a wind park after a predefined time span. As temperature is more and more a focus especially in biological and global change studies, furthermore inter-model variability in temperatures experienced by the virtual particles was determined. The main focus was on the transport variability originating from the physical models and thus biological behavior was not included. Depending on the scenario, median experienced temperatures differed by 3 °C between years. The range between the different models in one year was comparable to this temperature range observed between modelled years. Connectivity between flatfish spawning areas and the coastal protected area was highly dependent on the release location and spawning time. No particles released in the English Channel in the sole scenario reached the protected area while up to 20% of the particles released in the plaice scenario did. Interannual trends in transport directions and connectivity rates were comparable between models but absolute values displayed high variations. Most models showed systematic biases during all years in comparison to the ensemble median, indicating that in general interannual variation was represented but absolute values varied. In conclusion: variability between models is generally high and management decisions or scientific analysis using absolute values from only one single model might be biased and results or conclusions drawn from such studies need to be treated with caution. We further concluded that more true validation data for particle modelling are required.

Assessing the validity of a model is essential for its credibility especially when the model is used as decision making tool. Complex dynamic fishery models are recommended to investigate the functioning of fisheries and to assess the impact of management strategies, particularly spatial fishing regulations. However, their use is limited due to the difficulty and computational cost of parameterizing and gaining confidence, particularly for parameter rich models. These difficulties are compounded by uncertainty regarding parameter values, many of which are often taken from literature or estimated indirectly. Here we propose a methodology to improve confidence and understanding in the model, easily transferable to any complex model. The approach combines sensitivity analysis, scalability of parameters, optimization procedures, and model skill assessment in order to parameterize, validate and achieve the most plausible formulation of a model given the available knowledge while reducing the computational load. The methodology relies on five steps: (1) sensitivity analysis, (2) classification of parameters into a hierarchy according to their sensitivity and the nature of their uncertainty, (3) building of alternative formulations, (4) calibration and (5) skill evaluation. The approach is illustrated here by reviewing the parameterization of the ISIS-Fish model of the anchovy fishery in the Bay of Biscay. By using this approach, it is possible to make a thorough assessment of lacking information (e.g. accessibility to fishing and adult mortality) and to identify the strengths and weaknesses of the model in the context of different hypotheses. When applied to the ISIS-Fish model, the results suggest higher egg and adult mortality than formerly estimated, as well as new estimates for the migration towards spawning areas. They show the reliability of the model in terms of correlations with observations and the need for further efforts to model purse seiner catches. The methodology proved to be a cost-efficient tool for objectively assessing applied model validity in cases where parameter values are a mix of literature, expert opinion and calibration.

International audience; Multiple year oceanographic simulations (hindcast) are identified as a priority oceanography product for fisheries and environment studies since they provide a unique continuous long-term dataset allowing integrated assessment of the ocean state and evolution. We performed a 37 year (1972-2008) hindcast run with a coupled physical-biogeochemical model in the Bay of Biscay. The coupled model and the hindcast configuration are described. A model skill assessment is performed with a large set of in-situ data. Average seasonal currents show major circulation patterns over the shelf. Among tracers, temperature and salinity have the best agreement, ahead of nitrates and silicates, chlorophyll, and finally phosphates and ammonium. For chlorophyll, improved pattern statistics are found when compared to monthly composites of satellite-derived chlorophyll data. From the hindcast, we derived indices related to mesoscale activity (eddies, plumes, fronts, stratification) and production (chlorophyll and primary production). They help characterise the evolution of the environment in a functional way, on both the seasonal and multi-decadal scales. From these indices, first, a multivariate analysis reveals an increasing number of years that deviate from the mean seasonal pattern. Second, we propose interpretations of the simulated increasing trends detected in several of them (temperature, thermocline depth and primary production). We also recommend further developments to confirm these simulated evolutions, from addition of open boundary forcing with a global circulation model, to the improvement of the dynamics of nutrient regeneration and of the seasonal variability of secondary production. As a perspective, we review the different applications made from our hindcast in relation to anchovy life cycle, a species of major interest in the Bay of Biscay.

International audience; Simulating fish larval drift helps assess the sensitivity of recruitment variability to early life history. An individual based model (IBM) coupled to a hydrodynamic model was used to simulate common sole larval supply from spawning areas to coastal and estuarine nursery grounds at the meta-population scale (4 assessed stocks), from the southern North Sea to the Bay of Biscay (Western Europe) on a 26-yr time series, from 1982 to 2007. The IBM allowed each particle released to be transported by currents, to grow depending on temperature, to migrate vertically depending on development stage, to die along pelagic stages or to settle on a nursery, representing the life history from spawning to metamorphosis. The model outputs were analysed to explore interannual patterns in the amounts of settled sole larvae at the population scale; they suggested; (i) a low connectivity between populations at the larval stage, (ii) a moderate influence of interannual variation in the spawning biomass, (iii) dramatic consequences of life history on the abundance of settling larvae and (iv) the effects of climate variability on the interannual variability of the larvae settlement success. (C) 2015 Elsevier B.V. All rights reserved.

To anticipate the response of fish populations to climate change, we developed a framework that integrates requirements in all life stages to assess impacts across the entire life cycle. The framework was applied on plaice (Pleuronectes platessa) and Atlantic herring (Clupea harengus) in the North Sea, Atlantic cod (Gadus morhua) in the Norwegian/Barents Seas and European anchovy (Engraulis encrasicolus) in the Bay of Biscay. In each case study, we reviewed habitats required by each life stage, habitat availability, and connectivity between habitats. We then explored how these could be altered by climate change. We documented environmental processes impacting habitat availability and connectivity, providing an integrated view at the population level and in a spatial context of potential climate impacts. A key result was that climate-driven changes in larval dispersion seem to be the major unknown. Our summary suggested that species with specific habitat requirements for spawning (herring) or nursery grounds (plaice) display bottlenecks in their life cycle. Among the species examined, anchovy could cope best with environmental variability. Plaice was considered to be least resilient to climate-driven changes due to its strict connectivity between spawning and nursery grounds. For plaice in the North Sea, habitat availability was expected to reduce with climate change. For North Sea herring, Norwegian cod and Biscay anchovy, climate-driven changes were expected to have contrasting impacts depending on the life stage. Our review highlights the need to integrate physiological and behavioural processes across the life cycle to project the response of specific populations to climate change.

For many marine fish species, recruitment is strongly related to larval survival and dispersal to nursery areas. Simulating larval drift should help assessing the sensitivity of recruitment variability to early life history. An individual-based model (IBM) coupled to a hydrodynamic model was used to simulate common sole larval supply from spawning areas to coastal and estuarine nursery grounds at the population scale in the eastern Channel on a 14-yr time series, from 1991 to 2004. The IBM allowed each particle released to be transported by currents from the hydrodynamic model, to grow with temperature, to migrate vertically giving stage development, and possibly to die according to drift duration, representing the life history from spawning to metamorphosis. Despite sensitivity to the larval mortality rate, the model provided realistic simulations of cohort decline and spatio-temporal variability of larval supply. The model outputs were analysed to explore the effects of hydrodynamics and life history on the interannual variability of settled sole larvae in coastal nurseries. Different hypotheses of the spawning spatial distribution were also tested, comparing homogeneous egg distribution to observation and potential larval survival (PLS) maps. The sensitivity analyses demonstrated that larval supply is more sensitive to the life history along larval drift than to the phenology and volume of spawning, providing explanations for the lack of significant stock–recruitment relationship. Nevertheless, larval supply is sensitive to spawning distribution. Results also suggested a very low connectivity between supposed different sub-populations in the eastern Channel.

The abundance and spatial occupation of European anchovy Engraulis encrasicolus have in- creased in the North Sea since the mid-1990s. We use a cross-disciplinary approach combining genetics, transport modelling, survey time series analyses and physical oceanographic modelling to investigate 3 hypotheses on the reasons for this change. Evidence from connectivity studies suggests that the population of North Sea anchovy is separate from that in the Bay of Biscay. The recruitment pulses observed in survey data fit a life cycle which includes spawning in early summer and larval development in late summer. This also supports the concept of population expansion originating from local remnant population(s). In terms of growth physiology, suitable thermal windows have expanded, making conditions more favourable for life cycle closure and population persistence/productivity. In addition to the increased frequency of warm sum- mers, which favour larvae and juvenile growth, the decrease in the number of severe winters is also likely to improve overwinter survival. Overall, the evidence supports the hypothesis that the increase in anchovy abundance originated from the improved productivity of existing populations. This increase was associated with an expansion in thermal habitats and is probably not due to a northward shift in the distribution of southern conspecifics.

Hinrichsen, H-H., Dickey-Collas, M., Huret, M., Peck, M. A., and Vikebø, F. B. 2011. Evaluating the suitability of coupled biophysical models for fishery management. – ICES Journal of Marine Science, 68: 1478–1487. The potential role of coupled biophysical models in enhancing the conservation, management, and recovery of fish stocks is assessed, with emphasis on anchovy, cod, herring, and sprat in European waters. The assessment indicates that coupled biophysical models are currently capable of simulating transport patterns, along with temperature and prey fields within marine ecosystems; they therefore provide insight into the variability of early-life-stage dynamics and connectivity within stocks. Moreover, the influence of environmental variability on potential recruitment success may be discerned from model hindcasts. Based on case studies, biophysical modelling results are shown to be capable of shedding light on whether stock management frameworks need re-evaluation. Hence, key modelling products were identified that will contribute to the development of viable stock recovery plans and management strategies. The study also suggests that approaches combining observation, process knowledge, and numerical modelling could be a promising way forward in understanding and simulating the dynamics of marine fish populations.

Monitoring the environment of fish is a key component of the ecosystem approach to fisheries management. Here, we propose a methodology to statistically monitor time series of maps of environmental indices. These maps were derived from a 37-year hindcast of a coupled physical–biogeochemical model. The space–time variability in the maps was decomposed using empirical orthogonal functions into time-invariant spatial patterns and time-varying amplitudes for these patterns. A statistical process control chart was then applied to the time series of the amplitudes. In that way, changes were detected with known statistical performance. The monitoring system also specified how long the changes lasted and which years, seasons and zones were affected. Results for all indices were assembled in an integrative dashboard of the detected deviations. For illustration, the procedure was applied to the Bay of Biscay pelagic zone. The selected environmental indices characterized the evolution of hydrological structures such as fronts and river plumes, as well as changes in temperature, water column stratification, horizontal current flow and primary production. A major result was that, in the last decade, sea surface temperature showed repeated significant shifts towards warming, which were largest in the northern half of the Bay, while the spatial extension of river plumes over the shelf alternated between wet and dry years. From 2005, several other indices showed repeated significant deviations: increase in sea bottom temperature, increase in the depth of the pycnocline and changes in coastal currents. The procedure provided an integrated view of ecosystem variability and change for all its components and their spatial organization.

Monitoring and predicting the biogeochemical state of the ocean and marine ecosystems is an important application of operational oceanography that needs to be expanded. The accurate depiction of the ocean’s physical environment enabled by Global Ocean Data Assimilation Experiment (GODAE) systems, in both real-time and reanalysis modes, is already valuable for various applications, such as the fishing industry and fisheries management. However, most of these applications require accurate estimates of both physical and biogeochemical ocean conditions over a wide range of spatial and temporal scales. In this paper, we discuss recent developments that enable the coupling new biogeochemical models and assimilation components with the existing GODAE systems, and we examine the potential of such systems in several areas of interest: phytoplankton biomass monitoring in open oceans, ocean carbon cycle monitoring and assessment, marine ecosystem management at seasonal and longer time scales, and downscaling in coastal areas. A number of key requirements and research priorities are then identified for the future. The GODAE systems will need to improve their representation of physical variables that currently are not yet considered essential, such as upper-ocean vertical fluxes that are critically important to biological activity. Further, the observing systems will need to be expanded in terms of in situ platforms (with intensified deployments of sensors for O2 and chlorophyll, and inclusion of new sensors for nutrients, zooplankton, micronekton biomass, and others), satellite missions (e.g., hyperspectral instruments for ocean color, light detection and ranging (LIDAR) systems for mixed-layer depths, wide-swath altimeters for coastal sea levels), and improved methods to assimilate these new measurements.

Processing SeaWiFS (Sea-viewing Wide Field-of-view Sensor) data provides useful information for the observation and modelling of the phytoplankton production of the Bay of Biscay. Empirical algorithms allow the retrieval of chlorophyll a and non-living Suspended Particulate Matter (SPM) concentrations. These data are used to constrain a coupled 3D physical–biogeochemical model of the Bay of Biscay continental shelf. Two issues are investigated, depending on the variable used, to constrain the winter to spring phytoplankton production for the year 2001. First, SPM data is used as forcing data to correct the corresponding state variable of our model. This allows the realistic simulation of the light limited bloom at the end of February 2001, as observed with SeaWiFS chlorophyll a images and from the NUTRIGAS field cruise. Second, chlorophyll a data is used for parameter estimation of the biogeochemical model. The ability of assimilating these data is tested to improve the simulation of strong blooms observed in late May 2001 in the Loire and Gironde plumes. A global optimization method (Evolutive Strategies) is adapted to the complete 3-D coupled model, in order to find the best set of parameters. The hydrological conditions during the bloom can be validated with data from the PEL01 field cruise. After selection of the most sensitive parameters, the method is tested with twin experiments. Then, the use of real SeaWiFS data reduces the model/data misfit by a factor of two, improving the simulation of bloom intensities and extensions. The sets of parameters retrieved in each plume are discussed.

Adult anchovies in the Bay of Biscay perform north to south migration from late winter to early summer for spawning. However, what triggers and drives the geographic shift of the population remains unclear and poorly understood. An individual-based fish model has been implemented to explore the potential mechanisms that control anchovy's movement routes toward its spawning habitats. To achieve this goal, two fish movement behaviors – gradient detection through restricted area search and kinesis – simulated fish response to its dynamic environment. A bioenergetics model was used to represent individual growth and reproduction along the fish trajectory. The environmental forcing (food, temperature) of the model was provided by a coupled physical–biogeochemical model. We followed a hypothesis-testing strategy to actualize a series of simulations using different cues and computational assumptions. The gradient detection behavior was found as the most suitable mechanism to recreate the observed shift of anchovy distribution under the combined effect of sea-surface temperature and zooplankton. In addition, our results suggested that southward movement occurred more actively from early April to middle May following favorably the spatio-temporal evolution of zooplankton and temperature. In terms of fish bioenergetics, individuals who ended up in the southern part of the bay presented better condition based on energy content, proposing the resulting energy gain as an ecological explanation for this migration. The kinesis approach resulted in a moderate performance, producing distribution pattern with the highest spread. Finally, model performance was not significantly affected by changes on the starting date, initial fish distribution and number of particles used in the simulations, whereas it was drastically influenced by the adopted cues.

Proceedings of the 15th French-Japanese Oceanography Symposium; International audience; Considerable research has been conducted in the past 15 years around the Eastern English Channel ecosystem. Data collected since the 1970s on the biotic and abiotic compartments have been collated and mapped out in the mid-2000. This spatially explicit information formed a sound basis to improve knowledge on, and model, the functioning and dynamics of key ecosystem compartments, with a focus on flatfish species and fisheries and their interactions with other sectors of activity (aggregate extractions, maritime traffic). The more recent ongoing works are dedicated to the integration of those processes into several complementary end-to-end ecosystem models.

A coupled three-dimensional physical–biogeochemical model was developed in order to simulate the ecological functioning of the Rio de la Plata estuary and plume. The biogeochemical model reproduces the nitrogen cycle between five compartments: dissolved inorganic nitrogen, phytoplankton, zooplankton, detritus and dissolved organic nitrogen. The coupling is tested in seasonal climatological configurations and for the particular year 1999. The circulation is forced with Parana and Uruguay rivers discharges, NCEP wind and tide. The biogeochemical model includes loads of inorganic and organic nitrogen from both rivers. The model reproduces the correct tidal amplitudes in the estuary, as well as the most outstanding features of the observed horizontal and vertical structures of the salinity plume. Simulated surface chlorophyll a concentrations exhibit maximum values all year long seaward of the turbidity front, between the 0.5 and 15 isohalines, in agreement with SeaWiFS images of the area. The model simulates well the low primary production in the light-limited highly turbid tidal river ( 20 gC / m 2 / yr ), the high production area in the frontal zone where it can reach 500 gC / m 2 / yr , and the nutrient-limited production in the outer estuary and inner shelf ( 300 gC / m 2 / yr ), with realistic values in each case. According to the 1999 model simulation, the tidal river is the location of organic nitrogen remineralization with a consequent increase of the inorganic pool. At the entrance of the frontal zone, inorganic nitrogen represents about 75% of the whole nitrogen pool, it is reduced to 50% at its sea end-member. The outer estuary has the same sink role for inorganic nitrogen, suggesting that organic nitrogen is the major form exported to the shelf.

Characterizing the space–time variability in spatial distributions as well as understanding its drivers is basic to designing robust spatial management plans. As a prerequisite, we analyse here how this variability relates to population dynamics in conjunction with environmental conditions. For that, spatio-temporal statistical approaches are needed but seldom used in fisheries science. To fill this gap, we showcase the usefulness of the method of empirical orthogonal functions (EOFs). Guidelines are given to apply the method on a series of gridded maps as derived from fisheries survey dataseries that now span over decades. The method is applied to the series, 2000–2012, of the spatial distributions of European anchovy in the Bay of Biscay at spawning time. Across the series, the EOF decomposition allowed to identify three main types of spatial distributions. One type corresponded to an extended distribution, another to a restricted distribution in core areas, and the third to a very coastal distribution. The coastal spawning distribution corresponded to a low population growth rate as it was never followed by a large recruitment in the subsequent year. We did not attempt to explain the spatial patterns per se but the drivers of change from one type of distribution to another. Stock size and fish size as well as bottom temperature and water column stratification were the covariates that controlled the variability in the spatial distributions over time. Further, the spatial distribution at spawning time related to recruitment in the following year, meaning that variability in the spatial distribution of spawning affected population dynamics. The typology of maps based on EOF decomposition summarized this spatial variability into spatial spawning configurations, which may serve spatial planning.

Acoustic surveys are carried out every year in the Bay of Biscay in spring onboard the French research vessel Thalassa. The objective of PELGAS surveys is to study the abundance and distribution of pelagic fish in the Bay of Biscay. The main target species are anchovy and sardine but they are considered in a multi-specific context and within an ecosystemic approach as they are located in the centre of pelagic ecosystem., la campagne PELGAS se déroule tous les ans entre la fin avril et le début juin pour évaluer les stocks de poisson pélagique (sardine, anchois, chinchard, sprat). Cette campagne est écosystémique et observe tous les compartiments de écosystème pélagique, depuis le plancton jusqu'au prédateurs supérieurs en passant par les paramètres physiques de l'eau (température, salinité, densité, oxygéne dissous...). Ce rapport détaille la partie résultats de l’évaluation des stock ainsi que les principales caractéristiques de l'environnement tel que rencontrés cette année.

Daewel, U., Hjøllo, S. S., Huret, M., Ji, R., Maar, M., Niiranen, S., Travers-Trolet, M., Peck, M. A., van de Wolfshaar, K. E. 2014. Predation control of zooplankton dynamics: a review of observations and models. – ICES Journal of Marine Science, 71: 254–271. We performed a literature review to examine to what degree the zooplankton dynamics in different regional marine ecosystems across the Atlantic Ocean is driven by predation mortality and how the latter is addressed in available modelling approaches. In general, we found that predation on zooplankton plays an important role in all the six considered ecosystems, but the impacts are differently strong and occur at different spatial and temporal scales. In ecosystems with extreme environmental conditions (e.g. low temperature, ice cover, large seasonal amplitudes) and low species diversity, the overall impact of top-down processes on zooplankton dynamics is stronger than for ecosystems having moderate environmental conditions and high species diversity. In those ecosystems, predation mortality was found to structure the zooplankton mainly on local spatial and seasonal time scales. Modelling methods used to parameterize zooplankton mortality range from simplified approaches with fixed mortality rates to complex coupled multispecies models. The applicability of a specific method depends on both the observed state of the ecosystem and the spatial and temporal scales considered. Modelling constraints such as parameter uncertainties and computational costs need to be balanced with the ecosystem-specific demand for a consistent, spatial-temporal dynamic implementation of predation mortality on the zooplankton compartment.

Linking lower and higher trophic levels requires special focus on the essential role played by mid-trophic levels, i.e., the zooplankton. One of the most relevant pieces of information regarding zooplankton in terms of flux of energy lies in its size structure. In this study, an extensive data set of size measurements is presented, covering parts of the western European continental shelf and slope, from the Galician coast to the Ushant front, during the springs from 2005 to 2012. Zooplankton size spectra were estimated using measurements carried out in situ with the Laser Optical Plankton Counter (LOPC) and with an image analysis of WP2 net samples (200 μm mesh size) performed following the ZooScan methodology. The LOPC counts and sizes particles within 100–2000 μm of spherical equivalent diameter (ESD), whereas the WP2/ZooScan allows for counting, sizing and identification of zooplankton from ~ 400 μm ESD. The difference between the LOPC (all particles) and the WP2/ZooScan (zooplankton only) was assumed to provide the size distribution of non-living particles, whose descriptors were related to a set of explanatory variables (including physical, biological and geographic descriptors). A statistical correction based on these explanatory variables was further applied to the LOPC size distribution in order to remove the non-living particles part, and therefore estimate the size distribution of zooplankton. This extensive data set provides relevant information about the zooplankton size distribution variability, productivity and trophic transfer efficiency in the pelagic ecosystem of the Bay of Biscay at a regional and interannual scale.

Modelling of fish movement behaviour within a heterogeneous marine environment is a challenging but also key issue for understanding the effect of environmental factors and climatic change on fish processes (growth, distribution, mortality, reproduction). Fish movement models have the capability to encompass the combined effect of environment and empirical knowledge of fish individuals as energy requirements, known preys and predators, swimming capacities into a unified framework. Following an Individual Based Model (IBM) approach, a fish movement model has been developed to simulate the active movement of adult anchovy in the Bay of Biscay (BoB) in response to the spatio-temporal variations of both biotic and abiotic factors, as well as its internal conditions based on a bioenergetics model.