Trophic dynamics in the northern Humboldt Current system: insights from stable isotopes and stomach content analyses

The northern Humboldt Current system (NHCS) off Peru is one of the most productive world marine regions. It represents less than 0.1% of the world ocean surface but presently sustains about 10% of the world fish catch, with the Peruvian anchovy or anchoveta Engraulis ringens as emblematic fish resource. Compared with other eastern boundary upwelling systems, the higher fish productivity of the NHCS cannot be explained by a corresponding higher primary productivity. On another hand, the NHCS is the region where El Niño, and climate variability in general, is most notable. Also, surface oxygenated waters overlie an intense and extremely shallow Oxygen Minimum Zone (OMZ). In this context, the main objective of this study is to better understand the trophic flows in the NHCS using both stomach content and stable isotope analyses. The study focuses on a variety of organisms from low trophic levels such as zooplankton to top predators (seabirds and fur seals). The approach combines both long-term and specific studies on emblematic species such as anchoveta, and sardine Sardinops sagax and a more inclusive analysis considering the 'global' food web in the recent years (2008 – 2012) using stable isotope analysis. Revisiting anchovy and sardine we show that whereas phytoplankton largely dominated anchoveta and sardine diets in terms of numerical abundance, the carbon content of prey items indicated that zooplankton was by far the most important dietary component. Indeed for anchovy euphausiids contributed 67.5% of dietary carbon, followed by copepods (26.3%). Selecting the largest prey, the euphausiids, provide an energetic advantage for anchoveta in its ecosystem where oxygen depletion imposes strong metabolic constrain to pelagic fish. Sardine feed on smaller zooplankton than do anchoveta, with sardine diet consisting of smaller copepods and fewer euphausiids than anchoveta diet. Hence, trophic competition between sardine and anchovy in the northern Humboldt Current system is minimized by their partitioning of the zooplankton food resource based on prey size, as has been reported in other systems. These results suggest an ecological role for pelagic fish that challenges previous understanding of their position in the foodweb (zooplanktophagous instead of phytophagous), the functioning and the trophic models of the NHCS. Finally to obtain a more comprehensive vision of the relative trophic position of NHCS main components we used stable isotope analyses. For that purpose we analyzed the δ13C and δ15N stable isotope values of thirteen taxonomic categories collected off Peru from 2008 - 2011, i.e., zooplankton, fish, squids and air-breathing top predators. The δ15N isotope signature was strongly impacted by the species, the body length and the latitude. Along the Peruvian coast, the OMZ get more intense and shallow south of ~7.5ºS impacting the baseline nitrogen stable isotopes. Employing a linear mixed-effects modelling approach taking into account the latitudinal and body length effects, we provide a new vision of the relative trophic position of key ecosystem components. Also we confirm stomach content-based results on anchoveta Engraulis ringens and highlight the potential remarkable importance of an often neglected ecosystem component, the squat lobster Pleuroncodes monodon. Indeed, our results support the hypothesis according to which this species forage to some extent on fish eggs and larvae and can thus predate on the first life stages of exploited species. However, the δ13C values of these two species suggest that anchoveta and squat lobster do not exactly share the same habitat. This would potentially reduce some direct competition and/or predation.

End-to-end modelling for an ecosystem approach to fisheries in the Northern Humboldt Current Ecosystem

This work represents an original contribution to the methodology for ecosystem models' development as well as the rst attempt of an end-to-end (E2E) model of the Northern Humboldt Current Ecosystem (NHCE). The main purpose of the developed model is to build a tool for ecosystem-based management and decision making, reason why the credibility of the model is essential, and this can be assessed through confrontation to data. Additionally, the NHCE exhibits a high climatic and oceanographic variability at several scales, the major source of interannual variability being the interruption of the upwelling seasonality by the El Niño Southern Oscillation, which has direct e ects on larval survival and sh recruitment success. Fishing activity can also be highly variable, depending on the abundance and accessibility of the main shery resources. This context brings the two main methodological questions addressed in this thesis, through the development of an end-to-end model coupling the high trophic level model OSMOSE to the hydrodynamics and biogeochemical model ROMS-PISCES: i) how to calibrate ecosystem models using time series data and ii) how to incorporate the impact of the interannual variability of the environment and shing. First, this thesis highlights some issues related to the confrontation of complex ecosystem models to data and proposes a methodology for a sequential multi-phases calibration of ecosystem models. We propose two criteria to classify the parameters of a model: the model dependency and the time variability of the parameters. Then, these criteria along with the availability of approximate initial estimates are used as decision rules to determine which parameters need to be estimated, and their precedence order in the sequential calibration process. Additionally, a new Evolutionary Algorithm designed for the calibration of stochastic models (e.g Individual Based Model) and optimized for maximum likelihood estimation has been developed and applied to the calibration of the OSMOSE model to time series data. The environmental variability is explicit in the model: the ROMS-PISCES model forces the OSMOSE model and drives potential bottom-up e ects up the foodweb through plankton and sh trophic interactions, as well as through changes in the spatial distribution of sh. The latter e ect was taken into account using presence/ absence species distribution models which are traditionally assessed through a confusion matrix and the statistical metrics associated to it. However, when considering the prediction of the habitat against time, the variability in the spatial distribution of the habitat can be summarized and validated using the emerging patterns from the shape of the spatial distributions. We modeled the potential habitat of the main species of the Humboldt Current Ecosystem using several sources of information ( sheries, scienti c surveys and satellite monitoring of vessels) jointly with environmental data from remote sensing and in situ observations, from 1992 to 2008. The potential habitat was predicted over the study period with monthly resolution, and the model was validated using quantitative and qualitative information of the system using a pattern oriented approach. The nal ROMS-PISCES-OSMOSE E2E ecosystem model for the NHCE was calibrated using our evolutionary algorithm and a likelihood approach to t monthly time series data of landings, abundance indices and catch at length distributions from 1992 to 2008. To conclude, some potential applications of the model for shery management are presented and their limitations and perspectives discussed.

A behavioral ecology of shermen: hidden stories from trajectory data in the Northern Humboldt Current System

This work proposes an original contribution to the understanding of shermen spatial behavior, based on the behavioral ecology and movement ecology paradigms. Through the analysis of Vessel Monitoring System (VMS) data, we characterized the spatial behavior of Peruvian anchovy shermen at di erent scales: (1) the behavioral modes within shing trips (i.e., searching, shing and cruising); (2) the behavioral patterns among shing trips; (3) the behavioral patterns by shing season conditioned by ecosystem scenarios; and (4) the computation of maps of anchovy presence proxy from the spatial patterns of behavioral mode positions. At the rst scale considered, we compared several Markovian (hidden Markov and semi-Markov models) and discriminative models (random forests, support vector machines and arti cial neural networks) for inferring the behavioral modes associated with VMS tracks. The models were trained under a supervised setting and validated using tracks for which behavioral modes were known (from on-board observers records). Hidden semi-Markov models performed better, and were retained for inferring the behavioral modes on the entire VMS dataset. At the second scale considered, each shing trip was characterized by several features, including the time spent within each behavioral mode. Using a clustering analysis, shing trip patterns were classi ed into groups associated to management zones, eet segments and skippers' personalities. At the third scale considered, we analyzed how ecological conditions shaped shermen behavior. By means of co-inertia analyses, we found signi cant associations between shermen, anchovy and environmental spatial dynamics, and shermen behavioral responses were characterized according to contrasted environmental scenarios. At the fourth scale considered, we investigated whether the spatial behavior of shermen re ected to some extent the spatial distribution of anchovy. Finally, this work provides a wider view of shermen behavior: shermen are not only economic agents, but they are also foragers, constrained by ecosystem variability. To conclude, we discuss how these ndings may be of importance for sheries management, collective behavior analyses and end-to-end models.

Trophic ecology of jumbo squid and predatory fishes in the Northern Humboldt Current System

This work provides a contribution to a better understanding of the trophic ecology of important predators in the Northern Humboldt Current System, the jack mackerel (Trachurus murphyi), the chub mackerel (Scomber japonicus) and the jumbo squid (Dosidicus gigas) by the characterization of the highly variable feeding patterns of these species at different spatiotemporal scales. We provided new knowledge on the comparative trophic behaviour of these species, defined as opportunistic in previous investigations. For that purpose we applied a variety of statistical methods to an extensive dataset of 27,188 non-empty stomachs. We defined the spatial organization of the forage fauna of these predators and documented changes in prey composition according to predators’ size and spatiotemporal features of environment. Our results highligh the key role played by the dissolved oxygen. We also deciphered an important paradox on the jumbo squid diet: why do they hardly forage on the huge anchovy (Engraulis ringens) biomass distributed of coastal Peru? We showed that the shallow oxygen minimum zone present off coastal Peru could hamper the co-occurrence of jumbo squids and anchovies. In addition, we proposed a conceptual model on jumbo squid trophic ecology including the ontogenetic cycle, oxygen and prey availability. Moreover we showed that the trophic behaviour of jack mackerel and chub mackerel is adapted to forage on more accessible species such as for example the squat lobster Pleurocondes monodon and Zoea larvae. Besides, both predators present a trophic overlap. But jack mackerel was not as oracious as chub mackerel, contradictorily to what was observed by others authors. Fish diet presented a high spatiotemporal variability, and the shelf break appeared as a strong biogeographical frontier. Diet composition of our fish predators was not necessarily a consistent indicator of changes in prey biomass. El Niño events had a weak effect on the stomach fullness and diet composition of chub mackerel and jack mackerel. Moreover, decadal changes in diet diversity challenged the classic paradigm of positive correlation between species richness and temperature. Finally, the global patterns that we described in this work, illustrated the opportunistic foraging behaviour, life strategies and the high degree of plasticity of these species. Such behaviour allows adaptation to changes in the environment.