Modeling the spatial and temporal population dynamics of the copepod Centropages typicus in the northwestern Mediterranean Sea during the year 2001 using a 3D ecosystem model
Introduction
The distribution of marine zooplankton results from the combined effects of abiotic and biotic environmental factors (e.g., hydrodynamics, temperature, trophic conditions) and the biological processes occurring in each species throughout its life history (e.g., physiology and growth, buoyancy, behavior). Modeling the spatial distributions of zooplankton taking into account features of their ontogenetic development was initiated in the 1980s, such as in the models of Calanus marshallae in the Oregon upwelling area (Wroblewski, 1980, Wroblewski, 1982), and increased in the framework of Global Ocean Ecosystem Dynamics (GLOBEC) program (Carlotti and Poggiale, 2010, Carlotti et al., 2000, Gentleman, 2002).
Zooplankton population models take into account different types of structuring parameters, such as age (Wroblewski, 1980); age within stage, as for Pseudocalanus sp. (Davis, 1984) and Euterpina acutifrons (Carlotti and Sciandra, 1989); weight, as for Calanus finmarchicus (Bryant et al., 1997, Slagstad, 1981); and stage and weight, as for Pseudocalanus sp. (Eisenhauer et al., 2009, Fennel, 2001, Stegert et al., 2007). The level of complexity required depends on the objectives of the study. On one hand, a reduction in model complexity simultaneously allows computation of population models in a 3D numerical environment and avoidance of numerical diffusion within the demographic structure (Gentleman et al., 2008). On the other hand, a mechanistic representation of biological processes related to environmental factors including biological properties such as physiological (Caparroy and Carlotti, 1996, Carlotti and Sciandra, 1989, Carlotti and Wolf, 1998, Fennel and Neumann, 2003, Neumann and Fennel, 2006) and behavioral aspects (Slagstad and Tande, 2007) is also required to correctly represent population dynamics. Consequently, environmental factors differentially modify the biological processes occurring in each developmental stage in a species. For instance, temperature modifies stage-specific physiological rates and, consequently, the rates of growth and development. It also triggers a shift in biological phases from full activity to dormancy and vice versa (Halsband-Lenk et al., 2004). Temperature-dependent rates must be taken into account when copepod population dynamics are simulated at a seasonal scale and in areas characterized by contrasting temperature gradients. Trophic interactions and nutrition utilization represent complex processes that are stage dependent because they depend strongly on the predator–prey size relationship (Hansen et al., 1994). The best survival conditions for different developmental stages may occur in different habitats, and zooplanktonic population maintenance therefore has to be understood in the context of the regional dynamics of physical and ecosystems dynamics, including the role of various mesoscale features.
Zooplankton populations in the northwestern Mediterranean Sea (NWMS) are also quite dependent on hydro- and ecosystem dynamics in different ecoregions (Alcaraz et al., 2007, Boucher et al., 1987, Champalbert, 1996, Gaudy and Champalbert, 1998, Molinero et al., 2008, Sabatès et al., 1989). The local wind regimes and topography lead to cyclonic circulation of Modified Atlantic Waters from the Ligurian to the Balearic Sea (Millot, 1990). The Northern Current (NC) flows westward from the Ligurian Sea to the Balearic Sea, inducing permanent seasonally fluctuating Ligurian and North Balearic fronts (Birol et al., 2010), acting either as a temporary dynamic barrier between the shelf and the open sea or allowing cross-shelf exchanges between the Gulf of Lions and the offshore domain (Qiu et al., 2010). Moreover, the region features coastal upwelling (Albérola and Millot, 2003, Millot, 1979, Pinazo et al., 1996) and anticyclonic eddies in the western part of the Gulf of Lions, depending on the occurrence of Mistral and Tramontane winds (Campbell et al., 2013, Hu et al., 2009). In addition, the Rhone River represents an important source of fresh water, locally modifying the seawater density in the plume extension (Estournel et al., 1997), and it is also an important source of nutrients and organic matter, enhancing biological productivity within the shelf region (e.g., Blanc et al., 1969, Diaz et al., 2008, Tusseau et al., 1998).
To our knowledge, previous 3D ecosystem models applied to the NWMS have represented mesozooplankton based only on a single generic variable (for instance, Crise et al., 1998, Fontana et al., 2009, Lévy et al., 1998). Studying trophic interactions between the copepod population and its prey and predator species at relevant spatial and time scales is possible using a copepod population model coupled to a 3D-resolved ecosystem model. Fennel and Neumann (2003) and Neumann and Fennel (2006) have successfully replaced a bulk zooplankton variable in a biogeochemical model of the Baltic sea with a stage-resolved copepod model.
To model zooplankton distributions in the NWMS, the pelagic plankton ecosystem model (Auger et al., 2011) implemented in the Eco3M platform (Baklouti et al., 2006a, Baklouti et al., 2006b) has been used, coupled with the eddy-resolving regional circulation model Symphonie (Marsaleix et al., 2008), which reproduces the mesoscale features of an area. The copepod population module (Eisenhauer et al., 2009) has been applied to the calanoid Centropages typicus (noted C. t. below), a dominant neritic species in the Mediterranean Sea (Mazzocchi et al., 2007, Siokou-Frangou et al., 2010). Its biology and ecology are rather well documented (see special issue by Carlotti and Harris, 2007) compared to that of other copepods. The general seasonal distribution is briefly summarized here. It usually shows a wide peak of abundance in the spring–summer period at coastal stations monitored in Corsican and Balearic waters, the Ligurian Sea, Gulf of Lions and Catalan Sea (see Mazzocchi et al., 2007). A few spatial studies have been carried out at different times of the year. In the Ligurian Sea, reports are available for surveys conducted in December 1990 (Licandro and Icardi, 2009), without any occurrence of C. t., in April–May 1986 (Pinca and Dallot, 1995), when copepodite and adult stages were distributed over a large area between the French Côte d'Azur and Corsica, and in April 2011 (Molinero et al., 2008), during a short survey across the Ligurian front from Nice. Although it is difficult to compare data from different years and collected using different sampling strategies and methods, these studies show that adult and copepodite stages present spatially heterogeneous distributions, with incomplete overlap between the stages. Molinero et al. (2008) measured patches with dimensions of 4 km for copepodites and of 10 km for adults. Similarly, along a transect across the NC front from Barcelona in May 1989, Saiz et al. (1992) observed the highest densities of C. t. adults concentrated inshore from the density front. In five surveys performed by Kouwenberg (1994) between September 1986 and May 1988 and two surveys conducted in March–April 1998 and in January 1999 by Youssara (2002) and Gaudy et al. (2003), respectively, all in the Gulf of Lions, it was shown that (i) the winter population was maintained in coastal areas of the Gulf of Lions at low densities, except in the Rhone mouth area, and (ii) the spring populations were higher on the shelf than in the near open sea. A few other surveys covering the Gulf of Lions and the Algero-Provençal Basin (Furnestin, 1960, Plounevez and Champalbert, 2000) showed that C. t. has the highest densities in summer on the shelf and in the open sea and the lowest densities in the Northern Current. In late summer, the highest densities were observed in the open sea. On the shelf, C. t. is usually replaced by Temora stylifera in late summer and fall. Concerning its biological characteristics, C. t. exhibits a generation time of two weeks under optimal food and temperature conditions, which can extend for more than one month at 13 °C. At the end of fall and winter, when the population shows its lowest densities, it is maintained with a few individuals in the late copepodite and adult stages surviving with a flexible diet (Calbet et al., 2007, Carlotti et al., 2007). C. t. neither displays diapausing eggs, nor diapausing copepodite stages.
The objectives of the present work are to provide a 3D perspective of C. t. population dynamics in response to hydrographic events and the associated trophic supplies from spring to fall and to identify the potential habitats of overwintering late copepodite stages and adults. Our main goal is to understand the mechanisms guiding the seasonal distribution patterns of the developmental stages using realistic patterns of circulation and trophic conditions. The annual dynamics are presented for two main oceanographic stations that are representative of shelf and open sea conditions, while the spatial distributions over the whole area are presented for a few key dates during the year.
Section snippets
Setup of the coupled model
The setup of the biological model comprises the stage-structured copepod population (SSP) model (Eisenhauer et al., 2009) and the pelagic plankton ecosystem (PPE) model (Auger et al., 2011). The method of coupling between the two models is one-way coupling from the PPE to the SSP model. The prey and predator fields for the SSP model are state variables from the PPE model (see details below), but the SSP model outputs do not influence the plankton fields. The state variables of the SSP and PPE
Results
Understanding the simulated C. t. population dynamics requires a multidimensional approach, including the spatial and temporal dimensions of biomasses and abundance patterns as well as stage development and body growth dimensions. To present our results in a step-by-step manner, we first detail the annual dynamics in two contrasting typical environments, one in a deep sea area, the other in a shelf area. Then, we analyze spatial patterns for three key periods.
Model parameterization
Compared to the parameter values of the SSP model applied to C. t. by Eisenhauer et al. (2009), a few parameter values were updated in the present work. The stage-dependent half-saturation constant (K) of the functional response was modified in accordance with the prey field constituted by low trophic levels in the biogeochemical model. Bacterial and microzooplankton biomasses contribute to the available food concentration for C. t., particularly during summer, when the vertical stratification
Conclusion and perspectives
The results of the model developed in the present work show the complex coupling between hydrodynamics, prey and predator distribution patterns and the population dynamics of a target species. The SSP model considering five stage classes is an acceptable trade-off between a tractable representation and taking into account the main ontogenic characteristics of the target species coupled with hydrodynamics and the environmental ecosystem. The setup of the one-way coupled ecosystem-SPP model
Acknowledgments
This research was supported by the SESAME project (EC Contract No. GOCE-036949, funded by the European Commission's Sixth Framework Program under the priority ‘Sustainable Development, Global Change and Ecosystems’) and by the CNRS-INSU-MERMEX program (in particular, the WP2 actions: SPECIMED and Ecosystem Modeling). The authors thank the referees for their constructive comments and recommendations, which have improved the readability and quality of the paper.
References (87)
- et al.
Circulation in the French Mediterranean coastal zone near Marseilles: the influence of wind and the Northern Current
Cont. Shelf Res.
(2003) - et al.
Physical control of zooplankton communities in the Catalan Sea
Prog. Oceanogr.
(2007) - et al.
The penetration of the Northern Current over the Gulf of Lions (Mediterranean) as a downscaling problem
Oceanol. Acta
(2001) - et al.
Investigation of mechanistic formulations depicting phytoplankton dynamics for models of marine pelagic ecosystems and description of a new model
Prog. Oceanogr.
(2006) - et al.
Investigation and sensitivity analysis of a mechanistic phytoplankton model implemented in a new modular numerical tool (Eco3M) dedicated to biogeochemical modelling
Prog. Oceanogr.
(2006) - et al.
Aspects of the seasonal variability of the Northern Current (NW Mediterranean Sea) observed by altimetry
J. Mar. Syst.
(2010) - et al.
The seasonal dynamics of Calanus finmarchicus: development of a three-dimensional structured population model and application to the northern North Sea
J. Sea Res.
(1997) - et al.
The feeding ecology of the copepod Centropages typicus (Kröyer)
Prog. Oceanogr.
(2007) - et al.
Nutrients and plankton spatial distributions induced by a coastal eddy in the Gulf of Lion. Insights from a numerical model
Prog. Oceanogr.
(2013) - et al.
The biology and ecology of Centropages typicus: an introduction
Prog. Oceanogr.
(2007)
Towards methodological approaches to implement the zooplankton component in “end-to-end” food-web models
Prog. Oceanogr.
Modeling zooplankton dynamics
ICES Zooplankton Methodology Manual
Development and growth rates of Centropages typicus
Prog. Oceanogr.
A seasonal three-dimensional study of the nitrogen cycle in the Mediterranean Sea: Part I. Model implementation and numerical results
J. Mar. Syst.
Biogeochemical and ecological functioning of the low-salinity water lenses in the region of the Rhone River freshwater influence, NW Mediterranean Sea
Cont. Shelf Res.
Zooplankton population model coupled to a biogeochemical model of the North Western Mediterranean Sea ecosystem
Ecol. Model.
The plume of the Rhone: numerical simulation and remote sensing
Cont. Shelf Res.
The Rhone River Plume in unsteady conditions: numerical and experimental results
Estuar. Coast. Shelf Sci.
Assimilation of SeaWiFS chlorophyll data into a 3D-coupled physical–biogeochemical model applied to a freshwater-influenced coastal zone
Cont. Shelf Res.
Space and time variations in zooplankton distribution south of Marseilles
Oceanol. Acta
Biomass, metabolism and nutrition of zooplankton in the Gulf of Lions (NW Mediterranean)
Oceanol. Acta
Numerical simulations of eddies in the Gulf of Lion
Ocean Model.
Reproduction, hatching success, and early naupliar survival in Centropages typicus
Prog. Oceanogr.
Copepod distribution in relation to seasonal hydrographics and spatial structures in the northwestern Mediterranean (Gulf of Lions)
Estuar. Coast. Shelf Sci.
The onset of a bloom after deep winter convection in the northwestern Mediterranean Sea: mesoscale process study with a primitive equation model
J. Mar. Syst.
River discharges of water and nutrients to the Mediterranean and Black Sea: major drives for ecosystem changes during past and future decades
Prog. Oceanogr.
Energy conservation issues in sigma-coordinate free-surface ocean models
Ocean Model.
Seasonal and interannual variations in phytoplankton production at DYFAMED time-series station, northwestern Mediterranean Sea
Deep-Sea Res. II
Temporal variability of Centropages typicus in the Mediterranean Sea over seasonal-to-decadal scales
Prog. Oceanogr.
The Gulf of Lions' hydrodynamics
Cont. Shelf Res.
Spatial distribution of the copepod Centropages typicus in Ligurian Sea (NW Mediterranean). Role of surface currents estimated by Topex–Poseidon altimetry
C. R. Biol.
Surface patterns of zooplankton spatial variability detected by high frequency sampling in the NW Mediterranean. Role of density fronts
J. Mar. Syst.
A method to represent vertical migration of zooplankton in 3D-Eulerian Models
Ocean Model.
Spatial and temporal variability of phytoplankton biomass in upwelling areas of the northwestern Mediterranean: a coupled physical and biogeochemical modelling approach
J. Mar. Syst.
Diet, feeding behavior and trophic activity of the anchovy (Engraulis encrasicolus L.) in the Gulf of Lion (Mediterranean Sea)
Oceanol. Acta
The influence of hydrodynamic processes on zooplankton transport and distributions in the Northwestern Mediterranean: Estimates from a Lagrangian model
Ecol. Model.
Structure and resilience of overwintering habitats of Calanus finmarchicus in the Eastern Norwegian Sea
Deep-Sea Res. II Top. Stud. Oceanogr.
21st century climate change scenario for the Mediterranean using a coupled Atmosphere-Ocean regional Climate Model
Glob. Planet. Chang.
Parameterisation of a zooplankton population model for Pseudocalanus elongatus using stage durations from laboratory experiments
Ecol. Model.
Interaction of currents and vertical migration in maintaining Calanus marshallae in the Oregon upwelling zone—a simulation
Deep Sea Res. A. Oceanogr. Res. Pap.
Zooplankton community during the transition from spring bloom to oligotrophy in the open NW Mediterranean and effects of wind events. 1. Abundance and specific composition
J. Plankton Res.
Functioning of the planktonic ecosystem on the Gulf of Lions shelf (NW Mediterranean) during spring and its impact on the carbon deposition: a field data and 3-D modelling combined approach
Biogeosciences
Cited by (4)
Implementation of an end-to-end model of the Gulf of Lions ecosystem (NW Mediterranean Sea). I. Parameterization, calibration and evaluation
2019, Ecological ModellingCitation Excerpt :The E2E approach developed in this study is based on a fully dynamic coupling (i.e. two-ways coupling) of two pre-existing sub-models representing the dynamics of LTL organisms driven by hydrodynamics and climate processes on the one hand, and the dynamics of HTL organisms impacted by fishing activities on the other hand. The first model component is the Eco3M-S/Symphonie model that has been successfully used in the North-Western Mediterranean Sea to advance our understanding of the influence of hydrodynamics and atmospheric drivers on the distribution of plankton at different spatial and temporal scales (e.g. Auger et al., 2011, 2014; Campbell et al., 2013; Carlotti et al., 2014). The second component model is the individual- and size-based model OSMOSE (Shin and Cury, 2004).
Assessing and predicting marine ecosystem changes in the Southern European Seas - SESAME
2014, Journal of Marine SystemsZooplankton community response to the winter 2013 deep convection process in the NW Mediterranean Sea
2017, Journal of Geophysical Research: Oceans