http://orcid.org/0000-0002-2491-6203
Figures
Abstract
Information on the buoyancy of eggs and larvae from deep-sea species is rare but necessary for explaining the position of non-swimming larvae in the water column. Due to embryonic morphology and ecology diversities, egg buoyancy has important variations within one species and among other ones. Nevertheless, it has hardly been explored if this buoyancy variability can be a strategy for deep-sea larvae to optimize their transport beyond their spawning areas. In the northwestern Mediterranean Sea, protozoea and mysis larvae of the commercial deep-sea shrimp Aristeus antennatus were recently found in upper layers, but to present, earlier stages like eggs and nauplii have not been collected. Using a Lagrangian transport model and larval characteristics, we evaluate the buoyancy and hydrodynamic effects on the transport of A. antennatus’ larvae in the northwestern Mediterranean Sea. The transport models suggested that 75% of buoyant eggs released between 500 and 800 m depth (i.e., known spawning area), reached the upper water layers (0–75 m depth). Then, according to the modeled larval drifts, three spawning regions were defined in the studied area: 1) the northern part, along a continental margin crossed by large submarine canyons; 2) the central part, with two circular circulation structures (i.e., eddies); and 3) the southern part, with currents flowing through a channel. The number of larvae in the most upper layer (0–5 m depth) was higher if the larval transport model accounted for the ascent of eggs and nauplii (81%) instead of eggs reaching the surface before hatching (50%). The larvae reaching the most water upper layer (0–5 m depth) had higher rates of dispersal than the ones transported below the surface layer (deeper than 5 m depth). The results of larval dispersal simulations have implications for the understanding of A. antennatus larval ecology and for management decisions related to the shrimp fisheries in the northwestern Mediterranean Sea.
Citation: Clavel-Henry M, Solé J, Kristiansen T, Bahamon N, Rotllant G, Company JB (2020) Modeled buoyancy of eggs and larvae of the deep-sea shrimp Aristeus antennatus (Crustacea: Decapoda) in the northwestern Mediterranean Sea. PLoS ONE 15(1): e0223396. https://doi.org/10.1371/journal.pone.0223396
Editor: Atsushi Fujimura, University of Guam, GUAM
Received: September 18, 2019; Accepted: December 11, 2019; Published: January 29, 2020
Copyright: © 2020 Clavel-Henry et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by the Ministerio de Economía, Industria y Competitividad, Gobierno de España (BES-2015-074126) to support her PhD within the CONECTA project (CTM2014-54648-C2-1-R) to MC-H. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Numerous species have a pelagic larval cycle which links the spawning places to the recruitment areas. Larval cycle is a relatively short time lapse compared to the life cycle of marine animal, but it is the phase when large dispersions occur [1]. For benthic species, the distribution of the species mostly relies on the transported larvae. Larvae have several mechanisms for positioning themselves in productive and favorable waters that optimize their growth and displacement [2]. According to those mechanisms, the larvae can be retained on the spawning places or connect to other areas that are of high interest for species with high commercial value. For that reason, many studies specifically addressed larval drifts in order to determine the efficiency of the fisheries management [3, 4].
Aristeus antennatus is a deep-sea shrimp with a high commercial value in the northwestern Mediterranean Sea. Since 1980, the reproductive cycle, biology, and the temporal and spatial dynamics of A. antennatus have been intensively studied based on data from commercial and scientific surveys [5–8]. The acquired knowledge has contributed to shape a local management plan which restricts fishing activity on A. antennatus since 2012 for the harbor with the highest landings (Palamós) [9, 10]. The current local management plan was implemented, partly assuming that the protected population of A. antennatus would increase the following years on the trawling grounds, and partly based on the idea that recruited juveniles are related to the spawners living on the trawling restricted area.
Nonetheless, like many deep-sea species, knowledge about the early-life stages of A. antennatus and their behaviors is scarce. In the case of A. antennatus, knowledge can be gathered from various larvae of Dendrobranchiata species, which is the taxonomic family of A. antennatus, e.g., larvae molt after the hatching of eggs with the following order: nauplii, protozoea, and mysis. Lecithotrophic stages (eggs and nauplii) of A. antennatus have not been observed yet, but its first planktotrophic stage (protozoea) is largely found in the surface layer [11]. Because ripe females spawn on the sea bottom, the embryonic and nauplius stages are assumed to have abilities (buoyancy, swimming ability) that take them to the surface. In [12], an Individual-Based Model (IBM) study showed that the near-bottom circulation has minor effect on passive individual dispersal (encompassing neutral eggs and inactive larvae of A. antennatus). The study indicated that near-bottom vertical currents do not advect individuals up to shallower layers, and suggested that the vertical distribution of eggs or larvae can only be explained by their capacity of moving across the water column. To our knowledge, the youngest stage of A. antennatus’ larvae found in superficial waters is the first substage of protozoea [11, 13]. This finding supports the hypothesis of positive buoyancy in the previous stages of A. antennatus’ larvae.
Buoyancy adjusts the vertical egg position in the water column by the difference between egg and water densities [14–16]. To date, because either some deep-sea species are gravid (i.e., eggs carried by the spawners) or egg information is unavailable, larval dispersal of deep-sea species rarely accounted for the simulation of the stages that interact with the water density. Within the Dendrobranchiata family or within other deep-sea Decapod species, diverse egg sizes and densities [17] were observed but were not comparable with the assumed characteristics of A. antennatus’ eggs. In advanced stages of Decapod larvae from various ecosystems (coastal, estuarine, and deep-sea), several patterns of vertical distribution and displacements have also been described [17]. For example, Penaeid nauplii can be found in superficial layers [18, 19] or near the bottom [20]. This biodiversity translates how challenging and delicate is the generalization of the larval behavior to each unknown larval stage of the deep-sea red shrimp.
Although little is known about the A. antennatus’ larvae, and as for other deep-sea species, it is likely that the vertical position of larvae is determinant for dispersal because the surface water has stronger velocity and higher temperature than deeper water layers [21, 22]. A number of hydrodynamic processes in the upper waters of the northwestern Mediterranean Sea can disperse or retain the larvae. The main circulation is driven by the Northern Current that flows southwestward between the surface and 250 m depth with maximum velocities around 0.30 m/s [23]. Several eddies circulate in a clockwise or anticlockwise direction and drive the distribution of larvae in the surface layer [24]. Those structures are unstably present during a few weeks and have different sizes and locations [25]. In summer, the Mediterranean Sea is also well-stratified by a thermocline around 15 m deep [26]. This thermocline forms a physical barrier for some small-scale ocean processes such as water mixing [27], and therefore can limit vertical displacement of Decapod larval stages [28, 29].
In this paper, we analyze the impact of buoyancy variability and three-dimensional water mass circulation on the drifts of eggs and larvae between spawning and recruitment areas for the benthic deep-sea shrimp A. antennatus species in the northwestern Mediterranean Sea. This work was an opportunity to test for the first time a protocol to define buoyancy values for undescribed larval stages of a deep-sea species and to provide new elements for the management plan makers.
Material and methods
An Individual-Based Model (IBM) for embryonic and larval stages of the deep-sea red shrimp (A. antennatus) was implemented with a Lagrangian particle-tracking framework. It simulated the early life-cycle behavior and dispersal patterns of shrimp larvae using 3D hydrodynamic model outputs of the northwestern (NW) Mediterranean Sea.
The hydrodynamic model
A climatological simulation of the NW Mediterranean Sea hydrodynamics was run using the Regional Ocean Modeling System (ROMS; [30]), a free-surface, terrain-following, primitive equation ocean model. A Mellor-Yamada 2.5 turbulence closure scheme is used for subgrid-scale mixing in the simulation [31]. Different models (IBM in [12], and a spatiotemporal food web model (Ecopath with Ecosim) in [32]) have successfully used the daily outputs of the ROMS model in the NW Mediterranean Sea and on the Valencian Gulf. The simulation domain ranged from 38° N to 43.69° N and from 0.65° W to 6.08° E (Fig 1), with a grid spacing of 2 km (with 256 x 384 grid points horizontally). The vertical domain is discretized using 40 vertical levels with a finer resolution near the surface (surface layer thickness between 0.49 m and 5.91 m). ROMS is built forced by a high resolution and accurate bathymetry of the western Mediterranean basin, which is fundamental for the drift study of A. antennatus’ eggs and larvae because adults are benthic and females spawn around 800 m depth [33].
The domain is divided by release zones (polygons from 1 to 12) and settlement zones (polygons from 1 to 24) along the continental margins of the management unit called Geographical Sub Area 6 (GSA 6) from the Food and Agriculture Organization. Eggs were spawned in the release zones at the seafloor between the 500 m and 800 m isobaths (black lines) where the mature female shrimp aggregates. All the zones (polygons from 1 to 24) are considered potential settlement zones of the deep-sea red shrimp. Arrows represent the average surface current (m/s) over July from the hydrodynamic model used in the study.
The ROMS outputs are validated in [12] and provided realistic products of the hydrodynamic and hydrographic characteristics (current velocities, salinity, and temperature) of the NW Mediterranean Sea. The upper water layers were characterized by the main southward current (i.e., the Northern Current) following the eastern coast of Spain and by several mesoscale circulations, like eddies. We estimated eddy size and center over monthly-averaged ROMS current velocities field with an algorithm described in [34]. Seawater density is a derived product of the salinity and temperature (equation 13 in [35]). On the horizontal dimension, average estimated temperature, salinity, and seawater density in summer (i.e., July to September) are 27.6°C, 37.5 PSU, and 1026.3 kg/m3 at the surface and 13.1°C, 38.2 PSU and 1035.35 kg/m3, at bottom, respectively. Those values are in the range of seawater measurements carried out in the NW Mediterranean Sea [23, 36].
The Individual-Based Model
The Lagrangian drift of individuals was calculated using the Python framework OpenDrift ([37]; available on https://github.com/OpenDrift). The equation for the advection of individuals in the three dimensions is resolved by a fourth-order Runge-Kutta scheme. Particles were defined as ‘stranded’ when they came into contact with the coastline. The modules for turbulence and buoyancy were activated by the modeler. The number of particles, the release coordinates and the duration of the drift setting the initial conditions of the Individual-Based Model are also described in the following sections.
Horizontal and vertical components.
The trajectories of individuals, representing virtual embryonic or larval stages, were based on the following equation:
(1)
where dX/dt is the 3D displacement of the individuals from their geographical and vertical position X on a time step dt; U is the advection component composed by the meridional, zonal and vertical velocities from ROMS; B is the biological component, and Dh and Dv are the small-scaled turbulent velocities in the horizontal and vertical dimensions, respectively.
The vertical velocity of individuals due to turbulent diffusivity (Dv) is adjusted by the random displacement scheme [38] with an internal time step of 20 s and a 1 m vertical resolution. Values of Kv were provided from the coefficient of salinity vertical diffusion from ROMS outputs. In our simulation, the vertical mixing was activated until individuals crossed the oceanic surface boundary layer (around 15 m deep [26]). Because the Mediterranean Sea is well-stratified in summer, above the oceanic surface boundary layer, the turbulence was accounted by simulating the horizontal mixing.
The horizontal velocity due to turbulent diffusivity was computed by a random walk scheme:
(2)
with rand representing a random value sampled in a Gaussian distribution G(0,1), dt the time step, and Kh the horizontal diffusivity coefficient. For Kh, we used the yearly constant and average value of 10 m2/s estimated by [39] in the Western Mediterranean Sea from a hydrodynamic model with similar horizontal resolution than ours.
The biological component B(X, t) represented the vertical terminal velocity due to the buoyancy force adjusted to the seawater density by the equilibrium of the Archimedes, gravitational, and friction forces under laminar flow (Reynolds number < 0.5) on a spherical object characterized by diameter and density. For the present study, this object size was randomly sampled in a Gaussian distribution with the average A. antennatus’ egg diameter of 3.3*10−4 m and its standard deviation at 0.45 *10−4 m [6]. The density of the sphere was selected after the analysis described in the sections below. The terminal velocity is computed by the equation provided in [40], or by Dallavalle equations [41] if the flow is transient (Reynolds number > 0.5).
Adaptation of the IBM to the shrimp larval cycle.
Following the temporal and spatial spawners distribution, virtual eggs [6, 12, 33] were set to be released in summer at midnight [42] on the bottom of the NW Mediterranean continental slope between 500 m and 800 m isobaths.
The IBM included activation and deactivation of the buoyancy force according to the individual stage during the drifts. The buoyancy module was initially activated at the release of the eggs. The eggs were positively buoyant by assuming that on the one hand, young larvae of A. antennatus are found in the surface layer [11, 13] and on the other hand, if the buoyancy is neutral or negative, the swimming abilities of the nauplius could not explain an average vertical rise on 600 m. When individuals reached a non-buoyant stage, the module was deactivated. For simplification, and due to lack of knowledge, we neglected the temporal change of egg size and density.
The drift duration was the sum of the pelagic egg and larval stage durations approximated by the temperature-dependent Pelagic Propagule Duration [43, 44]. Due to missing information regarding the embryonic and larval stages of A. antennatus, we assumed that its early-life development duration was similar to the near taxonomic Penaeid species. We reviewed 72 research articles (see S2 File) in which the larval stage duration of 42 Penaeid species was associated with the rearing water temperature. Then, we fitted a multiple linear model on those data to estimate A. antennatus larval stage durations (D, in day) according to the seawater temperature, such as:
(3)
where T is the rearing water temperature; Stage, a categorical variables characterizing the Penaeid larval stage and ε a random error. This model, whose initial assumptions were verified with a coefficient of determination R2 = 96%, had the following shape:
(4)
We estimated egg stage duration before the beginning of the drift simulations and larval stage duration was estimated each time an individual molted into the next stage (e.g., from nauplius to protozoea). The water temperature involved in Eq 4 was extracted at embryonic or larval position from ROMS model. Molting was allowed only if the simulation time was bigger than the cumulated duration of the previous and current stages.
Number of particles.
The methodology of [45] adjusted in [12] was used to determine the lowest number of individuals to be released and to guarantee that 95% of the dispersal variability is considered in our results. It is based on the average of the Fraction of Unexplained Variance [45], which is got by extracting randomly 100 subsamples of N drifts from 350,000 simulated trajectories and by computing the cross-correlation between the 100 subsamples. The subsample size N ranged between 1,000 and 300,000. The drifts of 350,000 individuals directly started at the surface with underneath sea bottom estimated between 500 and 800 m. Assuming to catch the full range of dispersal possibilities for larvae, the 350,000 drifts lasted the maximum PPD (38.8 days). The PPD was predicted from the relationship as described in the previous section, considering the coldest water near the bottom of the studied area (12.6°C) where eggs are spawned. Turbulent diffusivity was as well included for the drifts of 350,000 individuals. To avoid the underestimation of needed individual number, the 350,000 drifts also included a horizontal turbulent random walk with the maximum value of the coefficient Kh estimated at 100 m2/s in [39] instead of its average at 10 m2/s. Finally, the Fraction of Unexplained Variance was lower than 5% if 50,000 individuals are used (S1 Fig) to simulate the larval dispersal.
Larval dispersions
Modeled dispersion started with buoyant stages having predicted densities (in kg/m3) to rise toward upper water layers. When the stages were no more buoyant, the individuals were neutrally advected by the currents of the reached water mass layer.
Preliminary IBMs.
Beforehand, considering the few biological data available on taxonomically close species to A. antennatus and the important differences in larval ecology of deep-sea species (see S1 File), the density for the eggs reaching the surface needed to be estimated. This estimation was carried out relying on physical explanation (e.g., forces applied on the particles) to model the egg rise instead of unknown biological values.
To obtain the estimation (hence called optimal buoyancy), the depth reached by the individuals was analyzed when the buoyant stages ended. In the early summer (July 1), 50,000 particles were released and tracked up to the end of the nauplius stage (~ 7.5 days). The Lagrangian drift simulations of eggs were repeated for every 10 kg/m3 increment, between 800 and 1030 kg/m3. Then, for each particle reaching the surface (0–5 m depth layer), their densest value was kept. It represented the threshold where a denser value would not make the individual reach the surface and a lighter value would likely lead the individual to the surface. The average and standard deviation of the highest density values were kept to implement the drift simulation of the particles.
Other averages and standard deviations of the egg density were estimated by modifying one parameter in the configuration of the previous preliminary IBM. Independently, we computed those parameters i) for individuals reaching the upper boundary of the seasonal thermocline (about 15 m under the surface; [26]) instead of the surface layer, ii) when the buoyancy was applied during egg and nauplius stages, instead of egg stage only, iii) when the turbulence was activated, and iv) when the drift began at late summer (September 1). In case i), the method used to find the seasonal thermocline depth, we used the maximum slope of difference in the temperature profile during summer days [46]. In case ii), the particles were tracked up to the end of the nauplius stage (~7.5 days). In total, the preliminary experiment consisted of the analysis of simulations made with 24 values of egg density for six configurations of IBM (S1 Table).
Runs and sensitivity scenarios of IBMs.
We used six scenarios with different configurations of IBM. Sensitivity tests were performed introducing small or high variability due to the buoyant larvae or the hydrodynamics changes (Table 1). The scenario of reference IBM0 started in early summer at the beginning of the peak period of spawning and with buoyant eggs. We modified the average egg density to lead the individuals to different depths (surface or mixed layer depth), to use different buoyant stages (eggs or eggs + nauplii), to include the vertical and horizontal turbulent components, and to begin the drift simulation in either early or late summer (1st July or 1st September). In the last scenario IBMHot, the temperature in the upper 200 meters was incremented by 0.4°C, corresponding to the thermal increase expected in the upper layer of the western Mediterranean Sea over a decade [47].
Those scenarios were operated to consider the factors associated with larval ecology (IBM0, IBMPZ, and IBMMLD) and to the hydrodynamics (IBMLS, IBMTurb, IBMHot). Through IBMMLD, we explored the fact that the thermocline created a physical barrier for the larvae [28, 29]. In IBMPZ, we explored the possibility that buoyancy was still taken into account in the vertical movements of nauplius after hatching. The swimming abilities of the spheroid nauplius of Penaeid can hardly lead them to the surface [48, 49]. Yet, the nauplius stages of some coastal shrimps were found in the superficial layer [50], which means that an underlying process helps them to rise in the water column. Finally, the last scenarios allowed making a sensitivity analysis on the turbulent mesoscale impact on the main advection with the turbulent component (IBMTurb), and on the temporal changes in the physical environment (late summer hydrodynamic circulation) on the drifts (IBMLS). IBMhot was used to estimate the impact of a scenario of climate change on the larval drifts, through the rising of water temperature above average, which is expected to have a major impact on the larval duration and in seawater density.
For each scenario, the density (kg/m3) of the individuals was randomly sampled from a Gaussian distribution (Fig 2), using the average and standard deviation density values from the preliminary analysis. The simulation duration lasted the time that an egg developed into early juvenile according to the water temperature when larvae molted. The time step to advect individuals was one hour. The characteristics of the individuals (density, egg size, stages, duration of the stage) and their spatial position (latitude, longitude, and depth) were saved on a daily basis for further analyses.
Surfacing rate of buoyant eggs (black line) or buoyant eggs and nauplii (brown line) related to different tested densities (kg/m3). The average (vertical dashed lines) and standard deviation (colored ranges along the X-axis) of the water densities were used in different scenarios of Individual-Based Model (IBM). See details in Table 1.
Analysis of the drift simulations
The water characteristics during the drifts, the drift duration, the vertical rise of the individuals from spawning depths and the traveled distance were used to analyze the simulated drifts from the six scenarios. Distance d between two geographical positions given by the coordinates (X, Y) was computed by the haversine formula of great-circle distances;
(5)
with R as the radius of Earth. The drift distances corresponded to the aggregated distance during the simulation. The straight distance corresponded to the distance between the beginning and the end of the drifts.
In our study, we also analyzed the individual dispersals within and among areas with the connectivity matrices. Therefore, the NW Mediterranean Sea domain was divided into 24 zones (Fig 1) shaped by the main structures of topography like the Eivissa channel, the submarine canyons (Cap de Creus, Palamós, Blanes), and the gulfs (Valencian Gulf or Lion Gulf). In creating the zones we took account of the zoning criteria provided by the General Fisheries Commission for the Mediterranean (Geographical Sub-Areas 5, 6 and 7). Then, the connectivity rates between one release zone i and one settlement zone j were computed from the ratio of individuals in the settlement zone to the initial individuals from their release zone Nj|i/Ni. The release zones were the zones 1 to 12 while the settlement zones included all the 24 zones (Fig 1). In order to focus on the most relevant results, statistical tests and principal component analyses were implemented with averaged and scaled data with the basic packages of R.
Results
The simulated larval drifts significantly varied with the spawning places, the number of buoyant stages, the spawning period and the depth layer reached by buoyant phase.
Larval drifts in three spawning regions
The temporal hydrodynamic and individual buoyancy variations allowed setting relevant scenarios to use for connectivity analysis. Using the characteristics of the larval drifts ant the environmental influences (PPD, drift distances, water current, and water temperature) by scenarios, a Principal Component Analysis (PCA) differentiated three clusters of scenarios (Fig 3). One gathered all the IBM parameterized by small changes with respect to IBM0 (i.e., IBMTurb, IBMhot, and IBMMLD). The non-parametric Kruskal-Wallis test revealed that their PPD and drift distances were significantly different (p-values > 0.05). Nevertheless, their larval drifts expanded over a longer period (25.6 days) and over longer distances (118 km) even though larvae drifted at the thermocline depth (IBMMLD), or within small-turbulent water (IBMTurb), or in warmer water (IBMhot). The two other clusters were defined by the scenario IBMLS, in which larval drifts traveled the smallest distances (93 km), and by the scenario IBMPZ, in which the surfacing of larvae had higher rates (80%) and PPDs were the shortest (23.2 days). The present analysis focused on the drifts simulated in the scenarios IBM0, IBMPZ and IBMLS that were separated into distinct clusters and had important differences in their drift characteristics (S3 Fig).
The different scenarios of Individuals-Based Models used in this study are represented by blue symbols (circle, triangle or square) and visualized with the characteristics of larval drift and environmental influences (arrow). In the Principal Component Analysis (PCA), the correlation among the characteristics is indicated by the angle between arrows (i.e., an angle of 90° indicates no correlation and an angle of 180° indicates a negative correlation). Scenarios were correlated to the larval drift characteristics and environmental influences by the closeness of their Cartesian coordinates. The PCA grouped IBM0, IBMMLD, IBMHot, and IBMTurb (square), and separated IBMPZ (circle) and IBMLS (triangle). IBM scenarios are described in Table 1.
Generally, the larval drifts had similar trends for individuals released in specific latitudinal ranges of the NW Mediterranean Sea. For each scenario previously selected, a PCA (Fig 4) was carried out on the larval drift characteristics and environmental influences by release zones, gathering them into three regions. Besides, those three regions were driven by the same variables in each scenario. The region 1 gathered larvae released from zones 1 to 4 (approximately 41.39°N– 42.74°N) with the longest PPD and the straightest transports. The region 2 gathered individuals released between zones 5 and 10 (approximately 38.86°N– 41.39°N), which were abundant in the surface water but drifted the shortest distances. The region 3 gathered the individuals from the zones 11 and 12, and overlapped each side of the Eivissa Channel. Those individuals drifted the furthest and had low surfacing rates.
Release zones (colored symbols and numbers) linked to characteristics of the larval drift (arrows) in a Principal Component Analysis for A) the base scenario IBM0, B) the scenario with the buoyant stages up to Protozoea IBMPZ, and C) the scenario initialized at late summer IBMLS. The correlation among the characteristics of the larval drift is indicated by the angle between arrows (i.e., an angle of 90° indicates no correlation and an angle of 180° indicates a negative correlation). Release zones were correlated to the larval drift characteristics and environmental influences by the closeness of their Cartesian coordinates. The PCA grouped three regions: the release zones 1–4 (full green circle), the release zones 5–9 (full pink triangle) and the release zones 11 and 12 (full blue square). IBM scenarios are described in Table 1.
In each three regions, the individual drifts followed different circulation patterns. In region 1, Fig 5 shows an expanded and linear distribution of the individuals towards the Southwest. The dispersal of individuals from the region 2 shaped different eddies, which are almost-circular structures of the hydrodynamics. In all scenarios, those eddy-like structures were localized where the Gulf of Valencia widens (0.5° E; 39.5° N, and 2° E; 40.5° N). These eddies had only differences in their average radius and the position of their center changed between the early and late summer releases (IBM0 and IBMLS). For instance, these eddies get 6 to 8 km larger in IBMLS (S2 Table) and their center displaced 9.5 km to 14.7 km further (S2 Table) than eddies in IBM0. Last, individuals from region 3 were affected by the presence of an eddy with a center averagely positioning in the Eivissa channel (0.72° E; 38.70°N, S2 Table) in early summer (IBM0).
The early juveniles are positioned after the simulation with A) the base scenario IBM0, B) the scenario with the buoyant stages up to Protozoea IBMPZ, and C) the scenario initialized at late summer IBMLS. The regions are defined according to the results of a Principal Component Analysis (Fig 4). The black line represents the shelf break at 200 m. IBM scenarios are described in Table 1.
Advantage of two buoyant stages in the larval drifts
The highest surfacing rate regardless of the regions was when the ascent of larvae occurred with two buoyant stages. With an average density of 979 kg/m3 in IBMPZ, 93 kg/m3 denser than in IBM0, 81% of the eggs and nauplii rose to the shallower water layers, while in IBM0, 52% of the eggs reached the surface layer (Fig 6). Overall, the average time for the ascent of two buoyant stages (IBMPZ) lasted 7.5 days (i.e., 5.5 days longer than in IBM0). Nevertheless, at the end of the simulations in IBMPZ and for individuals at the surface, the drift durations increased by approximately two days and the drift distances were hardly 3.6 and 4.5 km longer (regions 1 and 2) or 12.6 km shorter (region 3) than in IBM0 (Table 2).
The depth position is recorded from the release day (Time = 0) to the end of the buoyant stages (grey arrow) in three Individual-Based Models. In IBM0 (black line) and IBMLS (grey dashed line), eggs were spawn at early and late summer, respectively, and buoyant eggs rose during 2 days. In IBMPZ (brown line), eggs and nauplii were buoyant and rose in the water column during 7.5 days.
The time lag between early and late arrival of individuals at the surface was small enough for not implying important divergences in the drift characteristics. Indeed, in the two scenarios IBM0 and IBMPZ, the surfaced larvae were not spatially exposed to different currents. When nauplii reached the upper 5 m layers in IBMPZ, they were 20.5 km away from the location where the eggs have surfaced in IBM0. However, this distance was mostly kept to the nearest kilometer between the individuals of IBM0 with the same age (7.5 days) than the surfaced individuals of IBMPZ. It illustrated that while the buoyant nauplii in IBMPZ were still rising in the water column, the surface current in which the individuals from IBM0 were advected, had not important amplitude or direction changes.
The overall consequences from tardive and abundant individuals at the surface were that the connectivity and retention strengthened between and within zones of the NW Mediterranean Sea. In both IBM0 and IBMPZ, the same settlement zones were connected but with different amplitudes (Fig 6). First, the dispersal rate between the release zones from the region 1 (zones 1 to 5) and from Blanes canyon and its neighboring zones (zones 5 to 7) rose from 42.6% in IBM0 to 56.7% in IBMPZ. Second, in the Gulf of Valencia’s zones (zones 7 to 9), corresponding to the south of region 2, the retention rates were from 33.4% to 43.6% higher in IBMPZ than IBM0. Last, each side of the Eivissa Channel (zone 11 and 12 in the region 3) connected with the other side if the surfacing was earlier like in IBM0. Indeed, in this scenario, the exchange of surfaced individuals between the two sides was bidirectional, with a tendency for individuals to cross northwardly the channel. A rate of 21.3% individuals from zone 11 connected to the northern zones of the Eivissa Channel (zone 9 and zone 12) and 11.5% individuals from the zone 12 crossed the channel and arrived on the southern zones of the Eivissa Channel (zones 10 and 11). While in IBMPZ, the exchange of late surfaced individuals across the channel was unidirectional and it was done by the northern zones of the channel. Under the influence of an eddy (Fig 4), the southern zones (zone 11 and 10) of the Eivissa Channel received 63.2% individuals from zone 12 and zone 11 kept 61.4% individuals (or 55.6% individuals more than in IBM0).
Larval drifts after early and late summer spawning
The present results distinguished the circulation in the water from the northernmost region (region 1) from the southern and warmer regions (regions 2 and 3). The drifts at the surface were exposed to higher temporal variability in the circulation fields if they started in the warmer regions of the NW Mediterranean Sea. In IBMLS, 49% of individuals (i.e., 3% fewer than in scenario IBM0) rose to the surface, and drifted 24.1 days (i.e., around one day less than in scenario IBM0). Yet, as shown in Fig 3, individuals from IBMLS traveled less. The temporal current changes were measured with the differences in drifted distance per day in scenario IBMLS. In region 2, surfaced larvae traveled 26 km less than in IBM0 (i.e., 3.1 km/day instead of 4.6 km/day) and in region 3, surfaced individuals traveled 96.4 km less (i.e., 3.7 km/day instead of 10.7 km/day). In contrast, the individuals in the colder region 1 had similar drifting velocity in IBMLS (i.e., 7.4 km/day) and in IBM0 (i.e., 7.6 km/day).
Besides the temporal changes in the intensity of the circulation fields, two circular structures in the NW Mediterranean Sea boosted the dispersal rates of the surface individuals. First, the diameter of an eddy-like structure between the latitudes 41°N and 42°N had modified the surface connectivity of individuals toward the Balearic Islands. In IBMLS, the north side of the Majorca Island (biggest island; zone 13 and 14) received 40.9% and 37.9% individuals from the Blanes canyon (zone 4) and its southern zone (zone 5 in Fig 7), while in IBM0, the reception rate of the individuals was lower than 3.5%. This new connection leaned on the changes in the current direction at the end of the summer and the size of the eddy-like structure. Individuals from the region 1 traveled a similar distance per day in IBM0 and IBMLS (7.5 km/day and 7.1 km/dayrespectively), but the circulation pattern had an angle closer to the south direction, which advected the individuals more southwardly in IBMLS (26° clockwise from geographical South) than in IBM0 (36° clockwise from geographical South). Fig 5 suggests that the eddy-like circulation over the north of the region 2 qualitatively had a bigger diameter. Second, the transport of surfaced individuals across the Eivissa Channel was limited without eddy structure. Indeed, the absence of an eddy in the Eivissa Channel (Fig 5C) paired with the fact that all surfaced individuals from zone 11 were transported unidirectionally across the Eivissa Channel. Additionally, most of the surfaced individuals from zone 12 mostly (25.2%) were advected in the eddy-like structure located at the south of the region 2.
Dispersal rates were calculated from release zones (X-axis) to potential settlement zones (Y-axis) from simulated larval drifts in A) the base scenario IBM0, B) the scenario with the buoyant stages up to Protozoea IBMPZ, and C) the scenario initialized at late summer IBMLS. For clarity purposes, we contoured the dispersal rates in regions 1 (green lines), 2 (pink lines) and 3 (blue lines). Identification of zones and regions as defined in Figs 1 and 4, respectively. Code color is related to the dispersal rate (%). The sum of the dispersal rates over a column gives the percentage of individuals in the surface layer.
Larval drifts in different depths of water layer
A part of the individuals did not reach the surface because of the variability in egg characteristics. For example, the decrease in the rate of surfaced larvae from region 3 between IBM0 and IBMLS was related to smaller egg diameters. Nonetheless, deeper drifts brought valuable knowledge on several features that are developed below.
First, the water temperature gradient in different depth layers was advantageous for longer transports. Due to the variability in egg size and egg density, the buoyant stages stabilized at different vertical positions from zero to 75 m for more than 75% of individuals. The drifts were 20 km longer because they lasted 11 supplementary days in the deeper layers regardless of the scenario involved. The highest gradient in the drift duration occurred in the upper 100 m layer where the drift lasted 12 supplementary days at 100 m than at the surface (Fig 8). Beyond 100 m, the drift duration was stabilized at 36–38 days. Nonetheless, the scale of the deeper drifts by region showed that the circulation was still stronger and weaker in the region 1 and 2, respectively, where the drifts were the longest (167.5 km) and the shortest (102.5 km), respectively and regardless the scenario. In a concrete case, the simulated drifts from region 3 in IBMLS indicated that the upper and deeper layer were decoupled because individuals drifted two times longer than the individuals in the upper layer (Table 2).
The average duration of propagules (from eggs to early juveniles) of Aristeus antennatus is in relation with the decreasing temperature (color gradient,°C) with increasing depths (Y-axis, m). In (A), PPD and water temperature from the IBM scenarios in early summer (IBM0 and IBMPZ) were averaged, and water temperature ranged between 13.1°C and 24.5°C. In (B), PPD and water temperature are from the IBM scenario in late summer (IBMLS), where water temperature ranged between 13.1°C and 23.5°C. IBM scenarios are described in Table 1.
Second, longer drifts and decoupled currents underneath the five meters depth favored the connectivity with the Balearic Islands (Fig 9). In other words, the velocity fields under the surface had slightly different directions. This was particularly seen by the drifts from the region 1 that, in two cases, provided higher arrival rates of individuals on the northwestern part of the Balearic Island grounds. One case is illustrated by the drifts in IBM0 showing 33.1% to 40% individuals from the zones 2 to 4 in the 5–330 m layer depth were transported toward the islands. The second case is showed in the IBMLS, with surfaced and deeper individuals reaching the zones 13 and 14, though the connection with zone 15 was established if only individuals were in the deeper layer (8 to 381 m depth). The same concept prevailed in the region of the Eivissa Channel. Within the first 100 m layer, longer drifts (160 to 400 km) from the Eivissa Channel connected to the south ground of the biggest Balearic Island (Majorca Island) by 39.5% individuals from zone 11 and 13.4% individuals from zone 12.
Dispersal rates were calculated from release zones (X-axis) to potential settlement zones (Y-axis) from simulated larval drifts in A) the base scenario IBM0, B) the scenario with the buoyant stages up to Protozoea IBMPZ, and C) the scenario initialized at late summer IBMLS. For clarity purposes, we contoured the dispersal rates in regions 1 (green lines), 2 (pink lines) and 3 (blue lines). Identification of zones and regions as defined in Figs 1 and 4, respectively. Code color is related to the dispersal rate (%). The sum of the dispersal rate over a column gives the percentage of individuals below the surface layer.
Finally yet importantly, the deeper and variable drifts of the individuals allowed broader connectivity. One release area could be connected up to eight other areas (Fig 9). Those connections were allowed by small numbers of individuals, which remained near the release areas, like for the case of individuals from region 1, which had a northwestward dispersal and thus, an opposite direction of dispersal than particles at the surface.
Discussion
Model results emphasized the importance of the buoyant phases and their variability in the larval cycle of deep-sea species. In the upper water layer of the northwestern Mediterranean Sea, the larval transport of A. antennatus relied on different circulation patterns and their spatiotemporal variability.
Importance of buoyant phases
The simulations of different buoyant phases in the larval drift of the red shrimp A. antennatus emphasized the importance of the lecithotrophic stage in the drifts. In the recommendation guide for IBMs by [51], buoyancy is a key parameter for explaining the changes in the larval drifts. However, prior to our study in the Mediterranean Sea, only [15] studied the impact of egg buoyancy on simulated drifts for the deep-sea shrimp Aristaeomorpha foliacea (taxonomically close to A. antennatus), using Atlantic anchovy egg density [52].
In our study, over half of the released individuals reached the layer where a high quantity of A. antennatus’ larvae is found (i.e., the surface layer [11]). We estimated an average egg density of 884 kg/m3, which was around 100 kg/m3 lighter than the ones described for other deep-sea shrimp eggs (1025–1082 kg/m3; [17]) and the seawater densities of the northwestern Mediterranean Sea (from 1025.3 kg/m3 in Cap de Creus to 1024.1 kg/m3 in the Eivissa channel; [23]). Nonetheless, egg density is partly species-dependent and therefore, egg density of deep-sea Decapod egg with bigger sizes than A. antennatus’ oocyte [17] may not be appropriate to use in our Lagrangian model. However, the displacement of A. antennatus’ eggs to the surface was comparable with some Penaeid eggs that were observed at the surface in laboratory water tanks and sea planktonic samples [50, 53], but this fact could not be quantitatively verified due to lack of Penaeid egg density data. To improve the next larval drift simulations with buoyant eggs of A. antennatus, the relationship between egg density and diameter (e.g., [17]) during the egg incubation should be considered. Indeed, when the eggs develop, several chemical reactions (i.e., the production of a perivitelline space by cortical activation of the eggs [54], the egg hydration [55], and embryonic development) likely modify the lipid rate, egg size, and egg density [53, 56]. Nevertheless, information about egg growth of deep-sea species and Penaeid eggs needs to be provided.
Extended to the nauplius stage, the rise of individuals illustrated the important role of the nauplius for increasing the probability to reach the surface layer. However, due to lack of knowledge, two parameters were not considered in the density values of IBMPZ: the loss of weight at hatching and the change of shape at hatching [57]. In the Penaeid family, the elliptical and hairy morphology of the nauplius leads to a change in the buoyancy effect [58], which may speed up or slow down the ascent of the larvae after hatching. Furthermore, at this stage, previous studies provided the implication of the swimming behavior in the control of the larval vertical migration. Indeed, although erratically [20], [59] and [60] observed a rise of the Penaeid nauplii followed by a resting time in which the nauplii slowly sank. The velocity of Penaeus larvae is poorly informed although [20] suggested that larvae rise quickly after hatching, guided by phototaxis. Nevertheless, all authors emitted doubts about the ability of nauplii to swim up to the surface. We underline the necessity to carry out studies on the locomotion of decapod shrimp larvae in the open ocean to define the potential balance between the buoyancy of the nauplii and their swimming ability.
Dispersal modulated by current variations
We noted that the main linear circulation along with eddy-like structures supported distinct dispersals of the red shrimp larvae according to their drift depths. Turbulence (in scenario IBMTurb) or individuals at the thermocline depth (in scenario IBMMLD) introduced small variations during the drift simulation that hardly modified the final dispersals of individuals. Dominant circulation patterns particularly controlled the larval dispersal rates according to three regions in our study (instead of two in [61]). In region 1, the spatial distribution of larvae near the surface was strongly influenced by the Northern Current and winds, as reported in other drift analyses of marine larval species with egg buoyancy on our study area [62, 63]. In regions 2 and 3, the larval drift was limited by local mesoscale eddies, which have high importance in retaining larvae [64]. Near Blanes Canyon, the beginning of Ebro shelf, and the Eivissa Channel, previous studies identified eddies potentially trapping larvae of various species such as Engraulis encrasicolus and Sardinella aurita close to the North side of the Ebro shelf [65–69]. Decapod larvae have never been reported trapped in eddies in our study area, but [13] and [70] showed that decapod larvae were also aggregated by eddies south of the Balearic Islands and on the Norway lobster (Nephrops norvegicus) in the Irish Sea.
New connection or reinforcement of connection between zones by larval transport and resulted from the temporal changes of the eddy characteristics (i.e., if present, the diameter and the spatial position). During summer, the displacement and intensity changes in eddies along the Iberian Coast [65, 68] modified the simulated drifts, and therefore, the connectivity between areas. The drifts were influenced by a cyclonic eddy in the western part of the Eivissa Channel, as described previously by [71]. Its presence changed the larval dispersal by preventing the larvae from crossing the channel. However, during our late summer simulation, the fishing ground of A. antennatus near the Cabrera Harbor of Majorca Island [72], was supplied by individuals from the Eivissa Channel. The present study modeled and described for the first time how the Palamós canyon or the southern Eivissa Channel fishing grounds could connect to the fishing grounds of Majorca Island. Nonetheless, exploring the temporal evolution of eddies and their influence on the larval drift is needed to consider the possibility of persistent connectivity.
Dispersals related to the location and position of larvae
Larval duration of surface A. antennatus’ individuals depended to the Modified Atlantic Water mass dynamics. In our study, the PPD of larvae drifting at the surface varied along with a latitudinal gradient of near-surface temperatures of the NW Mediterranean Sea, which were colder in the vicinity of Cap de Creus and warmer at the Eivissa Channel [73]. This gradient has already been outlined by the interface of the old and cold Modified Atlantic water with the hot and recent Modified Atlantic water [74]. Additionally, it conveyed that larvae spawn from one of those regions will have distinct drifts from the other and partially explained the delimitation of the three zones in our study. The late summer dispersal enlightened a change of seawater temperature, which, besides the change in atmospheric conditions, can be related to the position of those two water masses. Indeed, the old Modified Atlantic water goes southward driven by wind and the southern entrance of the recent Modified Atlantic water is modulated by the water circulation in the channels of the Balearic Sea [74].
The larval drift in deeper layers can be a positive strategy induced by the buoyant stages. First, it implied that when the hatching occurred in deeper water layers, the PPD was longer. Second, the larvae drifted in the core of the North Current, which is maximal between 10 and 100 m [75, 76]. Third, the underneath surface is advantageous to larval dispersal because it may provide a good growth environment with nutrients located near the deep chlorophyll maximum [77] and protection from predators [28]. Fourth, the dispersal rate was wider and revolved round close settlement zones as well as far away zones. Consequently, for the first time, our study presents a potential connectivity between the two distant fishing grounds of Palamós and Soller (in Majorca Island from the Balearic Islands) areas and the expected drift conditions to link those places, where catches of A. antennatus are relatively important [72, 77].
Conclusions
Our study analyzed the larval drift of A. antennatus under the hypothesis that either the nauplius stage or the protozoea stage is in the upper water masses. We considered concurrently the effects of two main factors on the connectivity that both contribute to the dispersal strategy of this benthic deep-sea species: buoyancy and water mass circulation. The buoyancy linked the eggs spawned at sea bottom to the larvae in the upper layers of the ocean. A consequence of the buoyancy was the split of the Lagrangian drifts into three latitudinal spawning areas in accordance with similar larval dispersal. Simultaneously, we found that connectivity patterns and retention rates were influenced by the presence of mesoscale structures such as meanders above the head of Blanes Canyon and eddies along the Valencian Gulf and in the Eivissa Channel. The drift of larvae in early and late summer followed the persistent main circulation pattern of the NW Mediterranean Sea, with temporally variable mesoscale structures adjusting the connectivity between the different zones. The application of buoyancy on different larval stages also highlighted that positions of the larvae in the water column influenced the dispersal direction and the intensity of connectivity between zones of the NW Mediterranean Sea. This study contributes to increasing the knowledge of A. antennatus’ egg and larval ecology, and illustrates the potential dispersal paths of the individuals during their pelagic life. The spatiotemporal influence of mesoscale circulations in the larval dispersal should be tackled in further studies to better understand the red shrimp population connectivity, and to upgrade advice towards fishing practices and fisheries management.
Supporting information
S1 Fig. Decreasing Pelagic Propagule Duration (PPD) with increasing water temperature.
Continuous lines are the fitted exponential curves on reviewed data (cross) for each stage: eggs (orange), nauplius (red), protozoa (brown), and mysis (black). Reviewed data of PPD and temperatures were extracted from published articles and referenced in S2 File.
https://doi.org/10.1371/journal.pone.0223396.s001
(TIF)
S2 Fig. Number of particles needed to decrease the Fraction of Unexplained Variance (FUV) below the statistical threshold of 0.05.
The FUV expressed the bias created by a low number of particles (FUV > 0.05) after dispersal simulation. Whiskers show the minimum and maximum of the 100 FUV calculated for each number of particles tested (1000, 2000, 5000, 10000, 50000, 100000, 200000, 300000). Horizontal dashed line is the statistical threshold of 0.05.
https://doi.org/10.1371/journal.pone.0223396.s002
(TIF)
S3 Fig. Kernel density estimates from 50000 simulated individuals in different scenarios.
The simulated transport characteristics were A) the drift distances, B) the Pelagic Propagule Duration, and C) the water temperature during the drifts. Black lines indicate the scenarios grouped together by the PCA (Fig 3) with IBM0, IBMMLD, IBMDiff and IBMHot in four different line types. Blue line indicates the scenario IBMPZ, and grey line indicates the scenario IBMLS. See Table 1 for a description of the IBM scenarios.
https://doi.org/10.1371/journal.pone.0223396.s003
(TIF)
S1 File. Information on eggs from deep-sea species and Penaeid species.
Deep-sea crustaceans (Decapoda) species and taxonomically close species to Aristeus antennatus with the depth ranges of adult distribution, egg diameter and egg density. Depth range extracted from Sealife Base (https://www.sealifebase.ca) or EOL (https://eol.org). * Studied species; n/a, non available data; in grey cells, species taxonomically close to Aristeus antennatus.
https://doi.org/10.1371/journal.pone.0223396.s004
(DOCX)
S2 File. Database on Penaeid embryonic and larval duration by stages.
T, water temperature (°C).
https://doi.org/10.1371/journal.pone.0223396.s005
(DOCX)
S1 Table. Frequence of surfacing individuals (in percentage) in the preliminary experiments.
Percentage of buoyant stages reaching the surface or the Mixed Layer Depth (MLD) for a given density (kg/m3) in Individual-Based Model modified by the release period (early summer and late summer) and by the activation of random turbulences.
https://doi.org/10.1371/journal.pone.0223396.s006
(DOCX)
S2 Table. Characteristics of eddies involved in the larval dispersals.
Three eddies were shaped by the larval dispersal simulation (see Fig 5). Coordinates of the eddy center were estimated with an algorithm from [34]. The edge of eddies was got from the streamlines. Computation of the eddy radius is based on the circle surface (Area = pi*Radius2 and approximated the estimated eddy size.
https://doi.org/10.1371/journal.pone.0223396.s007
(DOCX)
Acknowledgments
The main author thanks all the Norwegian researchers from NIVA (Oslo, NO) who contributed to this paper: Johannes Röhrs, Knut-Frode Dagestad, Trine Bekkby and her family. The main author also thanks the editing from Hanna Morrissette and the advice from Elizabeth North from HPL-UMCES in Maryland (US). This study was possible thanks to the Spanish project CONECTA (CTM2014-54648-C2-1-R) lead from the ICM-CSIC in Barcelona (ES).
References
- 1. Cowen RK, Paris CB, Srinivasan A. Scaling of connectivity in marine populations. Science (New York, NY). 2006;311(5760):522–27.
- 2.
Rothlisberg PC, Church JA. Processes Controlling the Larval Dispersal and Postlarval Recruitment of Penaeid Prawns. In: Sammerco PW, Heron ML, editors. The Bio-Physics of Marine Larval Dispersal. Washington DC: American Geophysical Union. 2013:235–52. https://doi.org/10.1029/CE045p0235
- 3. Criales MM, Cherubin LM, Browder JA. Modeling Larval Transport and Settlement of Pink Shrimp in South Florida: Dynamics of Behavior and Tides. Marine and Coastal Fisheries. 2015;7(1):148–76.
- 4. Siegel DA, Mitarai S, Costello CJ, Gaines SD, Kendall BE, Warner RR, et al. The stochastic nature of larval connectivity among nearshore marine populations. Proceedings of the National Academy of Sciences. 2008;105(26):8974–79. pmid:18577590
- 5. Company JB, Puig P, Sardà F, Palanques A, Latasa M, Scharek R. Climate Influence on Deep Sea Populations (Climate Influence on Deep Sea). PLoS ONE. 2008;3(1):e1431. pmid:18197243
- 6. Demestre M, Fortuno JM. Reproduction of the deep-water shrimp Aristeus antennatus (Decapoda: Dendrobranchiata). Marine ecology progress series. Oldendorf, 1992;84(1):41–51.
- 7. Maynou F. Environmental causes of the fluctuations of red shrimp (Aristeus antennatus) landings in the Catalan Sea. Journal of Marine Systems. 2008;71:294–302
- 8. Sardà F, Company JB, Bahamón N, Rotllant G, Flexas MM, Sánchez JD, et al. Relationship between environment and the occurrence of the deep-water rose shrimp Aristeus antennatus (Risso, 1816) in the Blanes submarine canyon (NW Mediterranean). Progress in Oceanography. 2009;82(4):227–38.
- 9.
Boletín Oficial del Estado (2013) Orden AAA/923/2013 de 16 de Mayo, nº 126, Sec. III, pp. 40016. http://www.boe.es/boe/dias/2013/05/27/pdfs/BOE-A-2013-5555.pdf.
- 10.
Boletín Oficial del Estado (2018) Orden APM/532/2018 de 26 de Mayo, nº 128, Sec. III, pp. 55045–51. https://www.boe.es/boe/dias/2018/05/26/pdfs/BOE-A-2018-7015.pdf
- 11. Carretón M, Company JB, Planella L, Heras S, García-Marín J-L, Agulló M, et al. Morphological identification and molecular confirmation of the deep-sea blue and red shrimp Aristeus antennatus larvae. PeerJ. 2010;7:e6063.
- 12. Clavel-Henry M, Solé J, Bahamon N, Briton F, Rotllant G, Company JB. Influence of the summer deep-sea circulations on passive drifts among the submarine canyons in the Northwestern Mediterranean Sea. Ocean Sci. 2019;15:1745–59.
- 13. Torres AP, Dos Santos A, Balbín R, Alemany F, Massutí E, Reglero P. Decapod crustacean larval communities in the Balearic Sea (western Mediterranean): Seasonal composition, horizontal and vertical distribution patterns. Journal of Marine Systems. 2014;138:112–26.
- 14. Lindley JA. Eggs and their incubation as factors in the ecology of planktonic crustacea. Journal of Crustacean Biology. 1997;17(4):569–76.
- 15. Palmas F, Olita A, Addis P, Sorgente R, Sabatini A. Modelling giant red shrimp larval dispersal in the Sardinian seas: Density and connectivity scenarios. Fisheries Oceanography. 2017;26(3):364–78.
- 16. Sundby S, Kristiansen T. The Principles of Buoyancy in Marine Fish Eggs and Their Vertical Distributions across the World Oceans. PLOS ONE. 2015;10(10):e0138821. pmid:26465149
- 17. Herring PJ. Observations on the embryonic development of some deep-living decapod crustaceans, with particular reference to species of Acanthephyra. Marine Biology. 1974;25(1):25–33.
- 18. Subrahmanyam CB. On the unusual occurrence of penaeid eggs in the inshore waters of madras. J. Mar. Biol Ass. India. 1965;7(1):83–8.
- 19. Temple RF, Fischer CC. Seasonal distribution and relative abundance of planktonic-stage shrimp (Penaeus spp.) in the northwestern Gulf of Mexico, 1961. Fisheries Bulletin. 1967;66:323–34.
- 20.
Gurney R. Larvae of decapod Crustacea. London: The Ray Society. 1942
- 21. Fiksen Ø, Jørgensen C, Kristiansen T, Vikebø F, Huse G. Linking behavioural ecology and oceanography: larval behaviour determines growth, mortality and dispersal. Marine Ecology Progress Series. 2007;347:195–206.
- 22. Kristiansen T, Jørgensen C, Lough R, Vikebø F, Fiksen Ø. Modeling rule-based behavior: habitat selection and the growth-survival trade-off in larval cod. Behavioral Ecology. 2009;20(3):490–500.
- 23. Salat J. The interaction between the Catalan and Balearic currents in the southern Catalan Sea. Oceanologica Acta. 1995;18(2):227–34.
- 24. Bakun A. Fronts and eddies as key structures in the habitat of marine fish larvae: opportunity, adaptive response and competitive advantage. Scientia Marina. 2006;70(S2):105–22.
- 25. Rubio A, Barnier B, Jordà G, Espino M, Marsaleix P. Origin and dynamics of mesoscale eddies in the Catalan Sea (NW Mediterranean): Insight from a numerical model study. Journal of Geophysical Research: Oceans. 2009;114:C06009.
- 26. Houpert L, Testor P, Durrieu de Madron X, Somot S, D’Ortenzio F, Estournel C, et al. Seasonal cycle of the mixed layer, the seasonal thermocline and the upper-ocean heat storage rate in the Mediterranean Sea derived from observations. Progress in Oceanography. 2015;132:333–52.
- 27.
Thorpe SA. An introduction to ocean turbulence. 1st ed. Cambridge University Press. 2007
- 28. Lloyd ÑJ, Metaxas A, deYoung B. Patterns in vertical distribution and their potential effects on transport of larval benthic invertebrates in a shallow embayment. Marine Ecology Progress Series. 2012;469:37–52.
- 29. Queiroga H, Blanton J. Interactions Between Behaviour and Physical Forcing in the Control of Horizontal Transport of Decapod Crustacean Larvae Advances in Marine Biology. 2004;47:107–214
- 30. Shchepetkin AF, McWilliams JC. The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model. Ocean Modelling. 2005;9(4):347–404.
- 31. Warner JC, Sherwood HG, Arango HG, Signell RP. Performance of four turbulence closure methods implemented using a generic length scale method. Ocean Modelling. 2005;8:81–113.
- 32. Coll M, Steenbeek J, Sole J, Palomera I, Christensen V. Modelling the cumulative spatial–temporal effects of environmental drivers and fishing in a NW Mediterranean marine ecosystem. Ecological Modelling. 2016;331:100–14.
- 33. Sardà F, Company JB, Castellón A. Intraspecific aggregation structure of a shoal of a western Mediterranean (Catalan Coast) deep-sea shrimp, Aristeus antennatus (Risso, 1816), during the reproductive period. Journal of Shellfish Research. 2003;22(2):569–79.
- 34. Nencioli F, Dong C, Dickey ., Washburn L, McWilliams JC. A Vector Geometry–Based Eddy Detection Algorithm and Its Application to a High-Resolution Numerical Model Product and High-Frequency Radar Surface Velocities in the Southern California Bight. Journal of Atmospheric and Oceanic Technology. 2010;27(3):564–79.
- 35. Fofonof NP, Millard RC Jr. Algorithms for the computation of fundamental properties of seawater. UNESCO Technical Papers in Marine Sciences. 1983;44:53 pp.
- 36. Manca B, Burca M, Giorgetti A, Coatanoan C, Garcia MJ, Iona A. Physical and biochemical averaged vertical profiles in the Mediterranean regions: An important tool to trace the climatology of water masses and to validate incoming data from operational oceanography. Journal of Marine Systems. 2004;48(1–4):83–116.
- 37. Dagestad K-F, Röhrs J. OpenDrift v1.0: a generic framework for trajectory modelling. Geosci. Model Dev. 2018;11:1405–20.
- 38. Visser A. Using random walk models to simulate the vertical distribution of particles in a turbulent water column. Marine Ecology Progress Series. 1997;158:275–81.
- 39. Sayol J-M, Orfila A, Simarro G, López C, Renault L, Galán A, et al. Sea surface transport in the Western Mediterranean Sea: A Lagrangian perspective. Journal of Geophysical Research: Oceans. 2013;118(12):6371–84.
- 40. Sundby S. A one-dimensional model for the vertical distribution of pelagic fish eggs in the mixed layer. Deep Sea Research Part A. Oceanographic Research Papers. 1983;30(6):645–61.
- 41.
Dallavalle JM. Micromeritics. The technology of fine particles. Pitman, New York. 1948:11–30.
- 42.
Tavares C, Martin JW. Suborder Dendrobranchiata Bate, 1988. In: Schram F, von Vaupel Klein C, editors, Treatise on Zoology—Anatomy, Taxonomy, Biology. The Crustacea, Volume 9 Part A. Leiden: Brill. 2010:99–164. https://doi.org/10.1163/ej.9789004164413.i-562
- 43. Hilário A, Metaxas A, Gaudron SM, Howell KL, Mercier A, Mestre NC, et al. Estimating dispersal distance in the deep sea: challenges and applications to marine reserves. Frontiers in Marine Science. 2015;2(6).
- 44. O'Connor MI, Bruno JF, Gaines SD, Halpern BS, Lester SE, Kinlan BP, et al. Temperature control of larval dispersal and the implications for marine ecology, evolution, and conservation. Proceedings of the National Academy of Sciences. 2007;104(4):1266–71. pmid:17213327
- 45. Simons RD, Siegel DA, Brown KS. Model sensitivity and robustness in the estimation of larval transport: A study of particle tracking parameters. Journal of Marine Systems. 2013;119–20:19–29.
- 46. Fiedler PC. Comparison of objective descriptions of the thermocline. Limnology and Oceanography: Methods. 2010;8(6):313–25.
- 47. Calvo EM, Coma R, Gili JM, Pascual J, Pelejero C, Ribes M, et al. Effects of climate change on Mediterranean marine ecosystems: the case of the Catalan Sea. Climate Research. 2011;50:1–29.
- 48. Chu KH, Sze CC, Wong CK. Swimming Behaviour during the Larval Development of the Shrimp Metapenaeus ensis (De Haan, 1844) (Decapoda, Penaeidae). Crustaceana. 1996;69(3):368–78.
- 49.
Martin JW, Criales MM, Dos Santos A. Dendrobranchiata. In: Martin JW, Høeg JT, Olesen J, editors. Atlas of Crustacean Larvae. Baltimore, MD: Johns Hopkins University Press. 2014
- 50. Price ARG. Distribution of Penaeid shrimp larvae along the Arabian Gulf coast of Saudi Arabia. Journal of Natural History. 1982;16(5):745–57.
- 51.
North EW, Gallego A, Petitgas P. Manual of recommended practices for modelling physical–biological interactions during fish early life. ICES Cooperative Research Report. 2009;295.
- 52. Parada C, van der Lingen C, Mullon C, Penven P. Modelling the effect of buoyancy on the transport of anchovy (Engraulis capensis) eggs from spawning to nursery grounds in the southern Benguela: An IBM approach. Fisheries Oceanography. 2003;12(3):170–84.
- 53. Dall W, Hill BJ, Rothlisberg PC, Staples DJ. The Biology of the Penaeidae. Advances in Marine Biology. 1990;27:489 pp.
- 54. Pongtippatee-Taweepreda P, Chavadej J, Plodpai P, Pratoomchart B, Sobhon P, Weerachatyanukul W, et al. Egg activation in the black tiger shrimp Penaeus monodon. Aquaculture. 2004;234(1–4):183–98.
- 55.
Pandian TJ. Reproduction and Development in Crustacea. In: CRC Press Taylor and Francis group, editors. Reproduction and Development in Aquatic Invertebrates. 2016:27–69.
- 56. Rosa R, Calado R, Narciso L, Nunes ML. Embryogenesis of decapod crustaceans with different life history traits, feeding ecologies and habitats: a fatty acid approach. Marine Biology. 2006;151(3):935–47.
- 57. Chu K, Ovsianico-Koulikowsky N. Ontogenetic changes in metabolic activity and biochemical composition in the shrimp, Metapenaeus ensis. Journal of Experimental Marine Biology and Ecology. 1994;183(1):11–26.
- 58. Coombs S, Boyra G, Rueda L, Uriarte A., Santos M, Conway D, Halliday N. Buoyancy measurements and vertical distribution of eggs of sardine (Sardina pilchardus) and anchovy (Engraulis encrasicolus). Marine Biology. 2004;145(5):959–70.
- 59. Heldt JH. La Reproduction chez les crustacés décapodes de la famille des Pénéides. Annales de l’Institut. Océanographique de Paris, 1938;28:1–206.
- 60.
Dobkin S. Larvae of the pink shrimp (Penaeus duorarum). Presented at Proceedings of the Gulf and Caribbean Fisheries Institute 12, Miami (United States). 1960
- 61. Rossi V, Ser-Giacomi E, López C, Hernández-García E. Hydrodynamic provinces and oceanic connectivity from a transport net- work help designing marine reserves, Geophys. Res. Lett. 2014;41:2883–91.
- 62. Mariani P, MacKenzie BR, Iudicone D, Bozec A. Modelling retention and dispersion mechanisms of bluefin tuna eggs and larvae in the northwest Mediterranean Sea. Progress in Oceanography. 2010;86(1–2):45–58.
- 63. Ospina-Álvarez A, Palomera I, Parada C. Changes in egg buoyancy during development and its effects on the vertical distribution of anchovy eggs. Fisheries Research. 2012;117–118:86–95.
- 64. Cowen RK, Lwiza KMM, Sponaugle S, Paris CB, Olson DB. Connectivity of Marine Populations: Open or Closed? Science. 2008;287(5454):857–9. pmid:10657300
- 65. Escudier R, Renault L, Pascual A, Brasseur P, Chelton D, Beuvier J. Eddy properties in the Western Mediterranean Sea from satellite altimetry and a numerical simulation. Journal of Geophysical Research: Oceans. 2016;121(6):3990–4006.
- 66.
Karimova S. Observation of the surface circulation of the Mediterranean Sea from space. Paper presented at the ESA Living Planet Symposium, Netherlands. 2016
- 67. Pinot JM, López-Jurado JL, Riera M. The CANALES experiment (1996–1998). Interannual, seasonal and mesoscale variability of the circulation in the Balearic Channels. Progress in Oceanography. 2002;55(3–4):335–70.
- 68. Rubio A, Arnau PA, Espino M, Flexas M, Jordà G, Salat J, et al. A field study of the behaviour of an anticyclonic eddy on the Catalan continental shelf (NW Mediterranean). Progress in Oceanography. 2005;66(2):142–56.
- 69. Sabatés A, Salat J, Raya V, Emelianov M. Role of mesoscale eddies in shaping the spatial distribution of the coexisting Engraulis encrasicolus and Sardinella aurita larvae in the northwestern Mediterranean. Journal of Marine Systems. 2013;111:108–19.
- 70. Hill AE, Brown J, Fernand L. The western Irish Sea gyre: a retention system for Norway lobster (Nephrops norvegicus)? Oceanologica Acta. 1996;19(2–3):357–68.
- 71. Lafuente J, Jurado J, Lucaya N, Yanez M, Garcia J. Circulation of water masses through the Ibiza Channel. Oceanologica Acta. 1995;18(2):245–54.
- 72. Carbonell A, Carbonell M, Demestre M, Grau A, Monserrat S. The red shrimp Aristeus antennatus (Risso, 1816) fishery and biology in the Balearic Islands, Western Mediterranean. Fisheries Research. 1999;44(1):1–13.
- 73.
La Violette PE. Overview of the major forcings and water masses of the western Mediterranean Sea. In La Violette PE, editors. Seasonal and Interannual Variability of the Western Mediterranean Sea. 1995:1–11. https://doi.org/10.1029/CE046p0001
- 74. López García MJ, Millot C, García-Ladona E. Surface circulation variability in the Balearic Basin. Journal of Geophysical Research. 1994;99:3285–96.
- 75. García-Ladona E, Castellón A, Font J, Tintore J. The Balearic current and volume transports in the Balearic basin. Oceanologica Acta. 1996;19(5):489–97.
- 76. Macías D, Stips A, Garcia-Gorriz E. The relevance of deep chlorophyll maximum in the open Mediterranean Sea evaluated through 3D hydrodynamic-biogeochemical coupled simulations. Ecological Modelling. 2014;281:26–37.
- 77. Gorelli G, B Company J, Sardà F. Management strategies for the fishery of the red shrimp Aristeus antennatus in Catalonia (NE Spain). Marine Stewardship Council Science Series. 2014;2:116–27.