Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Resilience of Key Biological Parameters of the Senegalese Flat Sardinella to Overfishing and Climate Change

  • Kamarel Ba,

    Affiliations Institut Sénégalais de Recherches Agricoles (ISRA)/Centre de Recherches Océanographiques de Dakar-Thiaroye (CRODT), Dakar, Sénégal, Université Cheikh-Anta-Diop (UCAD)/Institut Universitaire de Pêche et d’Aquaculture (IUPA), Dakar, Sénégal

  • Modou Thiaw,

    Affiliation Institut Sénégalais de Recherches Agricoles (ISRA)/Centre de Recherches Océanographiques de Dakar-Thiaroye (CRODT), Dakar, Sénégal

  • Najih Lazar,

    Affiliation Department of Fisheries, University of Rhode Island East Farm, Kingston, Rhode Island, United States of America

  • Alassane Sarr,

    Affiliation Université Cheikh-Anta-Diop (UCAD)/Institut Universitaire de Pêche et d’Aquaculture (IUPA), Dakar, Sénégal

  • Timothée Brochier,

    Affiliations Institut Sénégalais de Recherches Agricoles (ISRA)/Centre de Recherches Océanographiques de Dakar-Thiaroye (CRODT), Dakar, Sénégal, Institut de Recherche pour le Développement (IRD), UMR 195 Lemar, Dakar, Sénégal

  • Ismaïla Ndiaye,

    Affiliation Université Cheikh-Anta-Diop (UCAD)/Institut Universitaire de Pêche et d’Aquaculture (IUPA), Dakar, Sénégal

  • Alioune Faye,

    Affiliation Université Cheikh-Anta-Diop (UCAD)/Institut Universitaire de Pêche et d’Aquaculture (IUPA), Dakar, Sénégal

  • Oumar Sadio,

    Affiliation Institut de Recherche pour le Développement (IRD), UMR 195 Lemar, Dakar, Sénégal

  • Jacques Panfili,

    Affiliation Institut de Recherche pour le Développement (IRD), UMR Marbec, Université Montpellier 2, Montpellier Cedex 5, France

  • Omar Thiom Thiaw ,

    thiawo@hotmail.com

    Affiliation Université Cheikh-Anta-Diop (UCAD)/Institut Universitaire de Pêche et d’Aquaculture (IUPA), Dakar, Sénégal

  • Patrice Brehmer

    Affiliations Institut Sénégalais de Recherches Agricoles (ISRA)/Centre de Recherches Océanographiques de Dakar-Thiaroye (CRODT), Dakar, Sénégal, Institut de Recherche pour le Développement (IRD), UMR 195 Lemar, Dakar, Sénégal

Abstract

The stock of the Senegalese flat sardinella, Sardinella maderensis, is highly exploited in Senegal, West Africa. Its growth and reproduction parameters are key biological indicators for improving fisheries management. This study reviewed these parameters using landing data from small-scale fisheries in Senegal and literature information dated back more than 25 years. Age was estimated using length-frequency data to calculate growth parameters and assess the growth performance index. With global climate change there has been an increase in the average sea surface temperature along the Senegalese coast but the length-weight parameters, sex ratio, size at first sexual maturity, period of reproduction and condition factor of S. maderensis have not changed significantly. The above parameters of S. maderensis have hardly changed, despite high exploitation and fluctuations in environmental conditions that affect the early development phases of small pelagic fish in West Africa. This lack of plasticity of the species regarding of the biological parameters studied should be considered when planning relevant fishery management plans.

1. Introduction

In Senegal (West Africa), coastal pelagic fishes are the most important marine resource, accounting for about 64% of annual catches [1]. Small pelagic fishes, such as sardinella (Sardinella aurita and Sardinella maderensis), have long been known to play a key role in the food supply for local communities. The catches of marine fisheries reached 448,000 t in 2012, the highest on record in Senegal [2]. More than 288,000 t of the catch was made up of small pelagic fish, of which 86% were sardinella (119,000 t for S. aurita and 128,700 t for S. maderensis) [2]. Sardinella are short-lived, zooplankton feeders with similar shape and size (maximum fork-length ~30 cm). Sardinella stocks have been considered to be overexploited since 2006 [3] by small-scale fishing canoes, a few local purse seiners, large foreign fishing vessels operating within the framework of bilateral agreements and often illegal, unreported, unregulated (IUU) fishing fleets [4, 5].

S. aurita undertakes extensive north/south seasonal migrations, whereas S. maderensis is less migratory [6, 7, 8, 9], suggesting that this species is able to adapt more readily to environmental variations. The trophic and ontogenetic migrations of sardinella are related to the strong seasonal variations in the southeastern edge of the Canary Upwelling System, an eastern boundary upwelling current [7, 8]. Sardinella are found along the north and west African coast from the Mediterranean Sea to Cape Frio in Angola, mainly in coastal upwelling areas off northwest Africa (10–26°N), in the Gulf of Guinea off Ghana and the Ivory Coast and in the south from Gabon to Angola (0–18°S) [6]. This paper focuses on the northwest population.

Past studies described the S. maderensis condition index, reproduction and growth based on data collected more than 25 years ago [69]. This study updated these parameters using more recent data from small-scale fisheries and published datasets to detect any changes that may have occurred as a result of climate change and/or over-fishing. The amount of food available from the environment may well change the amount of food reserves fish can build up or the energy they can invest in growth or reproduction [10]. Adaptive responses to overfishing have already been shown for various species, for example by reproducing earlier (e.g. smaller length at first maturity) [11].

The results of this study were compared and discussed in relation to relevant past studies in various north tropical Atlantic areas.

2. Materials and Methods

2.1. Study area

The Senegalese Exclusive Economic Zone (EEZ) (West Africa) extends from 18° 00’ N and 20° 00’ W to 12° 15’ N and 16° 30’ W. The Cape Verde Peninsula divides the EEZ into two areas that have distinct topographic characteristics (Fig 1). In the northern part of the zone, the continental shelf is fairly narrow, with the edge running north–north-east. In the southern part of the zone, the continental shelf widens and the edge runs north to south. Senegal has a cold season (November to April) and a hot season (May to October) [9]. The monthly AVHRR Pathfinder sea surface temperature (SST) records (http://coastwatch.pfeg.noaa.gov) with a resolution of 5 km show that the mean SST along the northern and southern coast increased by 0.03°C per year over the 1985–2009 period [12]. Furthermore, the mean SST was higher along the southern coast (24.55°C ± 3.06°C) than the northern coast (23.41°C ± 2.85°C).

thumbnail
Fig 1. Map of the study area, fishing areas (dark ellipses) covered by sampling sites (Kayar, Mbour, and Joal) and overlay of spatial distribution and migration pattern of Sardinella maderensis (adapted from [13]).

“Grande Côte” is the northern part (Dakar to Saint-Louis including Kayar) with a narrow continental shelf and “Petite Côte” is the southern part (Dakar to Sine-Saloum) with a wider continental shelf. Casamance is further south between Gambia and Guinea-Bissau. The nursery areas are delimited by a long dashed line (starting at the coast) and the spawning areas are delimited by a continuous bold line with arrows and numbers (months).

https://doi.org/10.1371/journal.pone.0156143.g001

Seasonal upwelling occurs along the Senegalese coast between November and May, maintaining a relatively low sea surface temperature depending on the strength of the wind (~15–18°C) [13]. The upwelling becomes permanent further to the north (in Mauritania). This upwelling is closely related to atmospheric dynamics and is determined by the seasonal changes in the direction of the trade winds associated with the migration of the Intertropical Convergence Zone (ITCZ) [14]. The trade winds stop in Senegal from June to October and a moderate northward current brings warm (> 25°C) water to the coast, which can be seen as an extension of the North Equatorial Countercurrent.

2.2. Biological data

Flat sardinella catches landed by small-scale fisheries were sampled monthly from March 2012 to February 2013. Purse seine and gill nets were the main fishing gear used. This study sampled only specimens caught using purse seine nets to have the same fishing gear selectivity for monthly comparisons of the length distribution. When available, up to 500 individuals were sampled from the three main landing sites: Kayar, Mbour and Joal (Fig 1). Interviews with fisherman (36 interviews; with one per month per site) performed during the sampling campaign, indicated that canoes fishing with purse seine nets at Kayar covered the area 15° 50’ N to 14° 80’ N. Mbour and Joal had overlapping fishing zones ranging from 14° 30’ N (in Senegal) to 13°48’ N (in Gambian waters).

The fork length (FL) of the fish sampled was measured to the nearest 1 mm and the body weight (W) was determined to the nearest 0.01 g. One hundred fish were dissected and sexed each month at each site. The maturity stages were determined macroscopically using the Fontana [15] maturity scale (Table 1). Eviscerated fish and gonads were also weighed.

thumbnail
Table 1. A brief description of the maturity stages according to [15], used in this study.

https://doi.org/10.1371/journal.pone.0156143.t001

2.3. Allometric relationships and condition factor

Length-weight relationships of the form: (Eq 1) were used [16], where W is the weight and FL is the fork length as defined in section 2.2, and a and b are estimated parameters.

The condition factor was defined using the Le Cren index [16] from: (Eq 2) where

Wobs is the observed weight and Wth is the theoretical weight (the weight predicted by Eq 1).

2.4. Estimating the von Bertalanffy growth parameters

The available length frequency data were grouped into 2 cm class intervals arranged by month from March 2012 to February 2013 (Table 2). The data were then analyzed using the R ELEFAN package, a recently updated Electronic LEngth Frequency ANalysis method (Daniel Pauly, unpublished). The von Bertalanffy Growth Function (VBGF) modified to allow for seasonal growth oscillations [17] was used to describe the growth of S. maderensis: (Eq 3) where

Lt is the predicted length at age t, L is the asymptotic length, the mean length that fish in a given population would reach if they were to grow continuously, K is the curvature parameter of the VBGF [18] and t0 is the “age” of the fish at zero length if they had always grown as predicted by the Equation, and where: (Eq 3a) (Eq 3b) (Eq 3c) where

ts is the start of every seasonal cycle with respect to t = 0, C is the amplitude of seasonal growth oscillations [18] and WP is the winter point (the period of the year when growth is most reduced). C was fixed at 0.1 and WP at 0.15, in this study.

thumbnail
Table 2. Size distribution of Sardinella maderensis specimens (n = 7708) collected in Senegalese waters from March 2012 to February 2013.

https://doi.org/10.1371/journal.pone.0156143.t002

ELEFAN determines the growth curve passing through a maximum number of peaks by adjusting L and K to maximize the goodness of fit Rn: (Eq 4) where ASP (available sum of peaks) is the sum of all values of available peaks, and ESP (explained sum of peaks) is the sum of all peaks and troughs through which the growth curve passes [1920]. The growth performance index φ’ [21] was used to evaluate growth parameter estimates, as its values are similar for the same species [22]. φ’ is defined as: (Eq 5)

Because the value of t0 cannot be estimated from length-frequency data, a very approximate value of t0 was estimated by substituting L (in cm) and K (year-1) in the following Equation [23]: (Eq 6)

2.5. Reproductive biology: sex ratio, sexual maturity, and spawning periods

Size-dependent sex ratios were computed from monthly sample data and a chi-square test was used to assess seasonal variations in the sex ratio of S. maderensis.

The mean length at first maturity (L50) was estimated using only samples collected during the reproductive periods. L50 is the length at which 50% of the individuals are mature (≥ stage III, maturing, on the maturity scale), as estimated using the logistic function: (Eq 7) where “%M” is the percentage of mature fish,”FL” is the mid-point of each length class and”a” and”L50” are estimated parameters. The model was fitted using R software and least squares (quasi-Newton method).

The gonadosomatic index (GSI) is useful for analyzing the state of sexual maturation of female and male fish gonads. GSI was calculated for each sampled individual based on [24], where: (Eq 8) where wg is gonad weight and WEv is eviscerated fish weight.

2.6. Statistical analysis

All statistical analyses were performed using the “stats”, “pgirmess” and “agricolae” R packages, with a significance level of α < 0.05. In general, the statistical analyses were taken from [2527].

The Shapiro test was used the assumption of normality and Bartlett’s test was used for the assumption of homogeneity of variance as Student’s t-test assumes that the observations within each group are normally distributed and that the variances are equal in both groups. Student’s t-test was used to assess the effect of season (hot versus cold) on individual sizes, when the Shapiro and Bartlett’s tests confirmed normality and homogeneity of variance. Otherwise, non-parametric Wilcoxon tests were used.

In order to compare the three sites, the data were either analyzed using one-way analysis of variance (ANOVA) followed by Tukey’s honest significant difference (HSD) post-hoc tests or using Kruskal-Wallis tests depending on whether the data was normally distributed. The chi-square test (χ²) was used to compare the sex ratios and sexual maturity stages.

3. Results

3.1. Length-weight relationship and condition factor

The analysis of the relationship between the total weight and fork length (Eq 1) using a total of 5,660 specimens showed that it was isometric: a = 0.02 ± 0.04, b = 3.00 ± 0.01 and r2 = 0.93 (Fig 2). For the hot season only (from June to October), S. maderensis exhibited negative allometry (a = 0.02 ± 0.07; b = 2.95 ± 0.02; R2 = 0.89 and n = 2043) and for the cold season (November to May), growth showed also slightly negative allometry (a = 0.02 ± 0.04; b = 2.98 ± 0.01; R2 = 0.94 and n = 3617).

thumbnail
Fig 2. Length-weight relationships for Sardinella maderensis in Senegalese waters.

Blue triangles and red dots represent the length-weight relationships during the cold season and hot season, respectively. The black line is the model fitted to all measurements regardless of season.

https://doi.org/10.1371/journal.pone.0156143.g002

The monthly variation in the condition factor (CF) is shown in Fig 3. The average CF calculated from a total of 5,660 individuals was equal to 1.01 ± 0.11. Comparisons of CF between months using the Kruskal-Wallis test suggested significant differences (χ² = 1814.64, df = 11, p-value < 0.05), which were confirmed for several monthly pairs by pairwise comparisons using the Kruskal-Wallis Multiple Comparison test (Fig 3). During April and July, the CF was higher CF but no real seasonal trend was identified.

thumbnail
Fig 3. Condition Factor (CF) for Sardinella maderensis from March 2012 to February 2013.

Means are represented by green triangles inside the notched box plots. Red figures under the error bars are the sample sizes. Months with the same blue letters (A-E) are not significantly different according to the Kruskal-Wallis test.

https://doi.org/10.1371/journal.pone.0156143.g003

3.2. Size spectrum

The length-frequency distributions from Kayar, Mbour and Joal had two modes at each site (Fig 4). The first minor peak length occurred at 19 cm at Kayar and at 16 cm for both Mbour and Joal. The second major peak length occurred at 24 cm at Kayar and at 22 cm for both Mbour and Joal. Individuals caught at “Petite Côte” were smaller than those caught at “Grande Côte” (Fig 1).

thumbnail
Fig 4. Sardinella maderensis size spectra for Kayar (dashed horizontal), Mbour (black), and Joal (grey).

Sizes are represented by the fork length.

https://doi.org/10.1371/journal.pone.0156143.g004

3.3. Growth parameters

L was estimated at 33.40 cm using the Wetherall method [28] as implemented in ELEFAN in R (S1 Fig). This estimate was used with the fixed values of C and WP giving the best fit value of K = 0.35 year-1 (S2 Fig). Fig 5 shows the length-frequency plots in Table 2, which were re-expressed as “peaks” and “troughs”, with the growth curve superimposed. The growth performance index (φ'), corresponding to the estimated growth parameters, was 2.59, which is close to the mean φ' (Δ = 0.09) for this species (Table 3). The estimated values of L and K were plotted on an auximetric plot (Fig 6). The estimated t0 value obtained by substituting L and K into Eq 6 [23] was -0.46 year. S. maderensis grows fast in the first year: it measures 13 cm at 6 months and reaches 17 cm after one year (Fig 5). The growth rate decreases gradually with age and reaches 25 cm at 2 years and 29.5 cm at 3 years.

thumbnail
Fig 5.

Seasonalized von Bertalanffy growth curve (with L = 33.40 cm, K = 0.34 year-1, C = 1 year-1, WP = 0.15 and Rn = 0.26) of Sardinella maderensis superimposed over (a) the original length-frequency histograms and (b) the restructured length-frequency histograms. The black and grey bars are positive and negative deviations from the “weighted” moving average of two size classes and represent pseudo-cohorts. The red dot is the starting point through which the curve passes to fit the model by maximizing Rn (index of goodness of fit which is analogous, but not equivalent to the correlation coefficient r in linear regression (see Eq 4)).

https://doi.org/10.1371/journal.pone.0156143.g005

thumbnail
Fig 6. Auximetric plot for Sardinella maderensis in Senegalese waters using the data from Table 2.

Dots represent the growth performance index (φ') fitted against Log10 (K) on the y-axis and Log10 (L) on the x-axis. Light grey dots represent fish species in the same family (Clupeidae) while dark grey dots are other species. Red dots represent values for Sardinella maderensis from the literature and the black dot represents the value from this study.

https://doi.org/10.1371/journal.pone.0156143.g006

thumbnail
Table 3. Growth parameters of Sardinella maderensis, based on the results of various publications in West Africa.

N/A means not available; L-F means length-frequencies. The length is expressed as fork length.

https://doi.org/10.1371/journal.pone.0156143.t003

3.4. Sex ratio and length at first maturity

The sex ratio (SR) for all samples of S. maderensis was biased in favor of females (55.81%; n = 1068). The proportion of females was significantly higher in the hot season (SR = 58.97%, χ ² = 13.82, df = 1, p-value < 0.05) compared to the cold season (SR = 53.68%, χ ² = 3.45, df = 1, p-value = 0.06). The sex ratio broken down by size class was male-biased between 15 to 22 cm while females dominated under 15 cm and over 22 cm (Fig 7). Broken down by sexual maturity stages, the differences were only statistically significant in stages II and VI (Table 4). In stage VI, all individuals were females (n = 4). In stages I and V, there were more males while there were more females in stages II (p-value < 0.5, n = 166), III and IV.

thumbnail
Fig 7. Female (light grey) and male (dark grey) Sardinella maderensis sex ratio in 2 cm size classes.

Sizes are represented by the fork length.

https://doi.org/10.1371/journal.pone.0156143.g007

thumbnail
Table 4. The proportions of male and female Sardinella maderensis individuals at the different sexual maturity stages.

https://doi.org/10.1371/journal.pone.0156143.t004

The body length of sexually mature females ranged from 12 cm to 26 cm, while that of mature males ranged from 13 cm to 23 cm. The cumulative frequency curve of individuals that had reached or passed maturity stage 3 indicated that nearly 50% (L50) of females were mature at 17.50 cm and 50% of males were mature at 16.60 cm (Fig 8).

thumbnail
Fig 8. Length at sexual maturity of Sardinella maderensis.

Measured male (white dots), measured female (black dots). Expected distribution of mature individuals is represented by a solid line for females and a dashed line for males.

https://doi.org/10.1371/journal.pone.0156143.g008

3.5. Spawning seasons

This study confirmed that S. maderensis spawns throughout the year in Senegal. However, two main reproductive periods can be defined from the analysis of sexual maturity stages and GSI. In this study, the reproductive period was considered to be the whole cycle of gonad maturation up to spawning. The first reproductive period occurred from April to October. Within this period, gonad maturation occurred in April, June and September while spawning occurred mainly from June to August and in October. The analysis of sexual maturity stages showed that individuals in stages I, II, and VI were present in November and December, after the first reproductive period. A second reproductive period, with continuous spawning (more intense than during the first period), occurred from January to the end of February with the gonad maturation period in December.

The monthly record of sexual maturity stages showed that females in stage V were more abundant in January-February, June-August, and, to a lesser extent, October. The proportion of females at stages ≥ IV peaked in February, June, August, and October (Fig 9A). Stage III females were mainly present in May and September. Stage V males were present from March to July and from November to February with higher incidence in February. Stage III and IV males were found throughout the year mostly from May to November (Fig 9B).

thumbnail
Fig 9.

Female (a) and male (b) sexual maturity stages against month of year.

https://doi.org/10.1371/journal.pone.0156143.g009

Monthly changes in GSI (Fig 10) followed sexual maturity stage V (Fig 9). The GSI was significantly different between months both for females (χ² = 115.94, df = 11, p-value < 0.05) and for males (χ² = 203.99, df = 11, p-value < 0.05). GSI was highly variable between August and September and between January and February. There were two spawning peaks in June and February for males and females with lower levels of spawning from June to August. In October, there was a higher presence of stage V and higher GSI simultaneously compared to September and November both for males and females and, therefore, there was probably a short reproductive period. There was a resting period in September and from November to December.

thumbnail
Fig 10.

(a) male and (b) female gonadosomatic index (GSI) of Sardinella maderensis between March 2012 and February 2013. The means are represented by green triangles inside the notched box plots. Red figures under the error bars are the sample sizes. Months with the same blue letters (A-E) are not significantly different according to the Kruskal-Wallis test.

https://doi.org/10.1371/journal.pone.0156143.g010

4. Discussion

The current study was based on more than 6,000 samples collected along the northern and southern Senegalese coast in 2012–2013 aiming to update the state of knowledge of S. maderensis and compare findings with previous observations made in Senegal [79, 2931], Mauritania [32, 33], Cabo Verde [34], Gambia [35], Sierra Leone [36], Ghana and Ivory Coast [37], Cameroon [3840], Congo [41, 4244] and Nigeria [41].

The other locations where comparable studies have been accomplished had different characteristics: (a) coastal ecosystems without upwelling (Cabo Verde, Sierra Leone, Cameroon and Nigeria), (b) areas where the intensity and duration of the upwelling is different from Senegal (Mauritania, Ghana and Ivory Coast, and Congo) and (c) estuarine ecosystems with either brackish water (the Gambia river [35]) or hyperhaline water (Sine-Saloum estuary, Senegal [31]).

The spatial distribution and exploitation of S. maderensis, as well as its growth and reproductive patterns, are closely related to the dynamics of the upwelling systems [6]. Adaptation of the biological parameters of the species to environments with different upwelling intensity should be expected. On the northwest African coast, the upwelling is stronger than in the other areas, permanent in the north of Mauritania and from November to May further south in Senegalese waters. The results were also compared with earlier measurements for Senegal to determine whether the biological parameters had changed with an increase in SST of 0.03°C per year over the past 30 years [14] and with the dramatic increase in fishing pressure, with annual landings increasing ten-fold from ~10 000 tons in 1980 to 128 700 tons in 2012 [2].

The growth performance index of S. maderensis estimated in this study (φ' = 2.59) fell in the range (2.40 to 2.97) previously found in Senegal and in the other areas, with the exception of one very low value found in Congo [42] (Table 3). However, it would seem that the low value in the Congo was due to an error because this was the only value outside the range of all other studies including studies in the same area. Our results also confirm other observations in Senegal in 1986 [8] that S. maderensis reached 17 cm to 18 cm body length at the age of 12 months in Senegal based on scale aging, with the growth rate decreasing until the fish approached their asymptotic length at around 3 years [8]. The results of different methods of aging fish (length-frequency analysis (this study and [8, 2930]), scales or other calcified structures [8, 42]) converge closely to the same value of the growth performance index. This suggests that this index was not significantly affected by climate variability or increased fishing pressure.

Large individuals (> 26 cm) were common in catches from Kayar (northern coast) but rare in catches from Mbour and Joal (southern coast) (Fig 1). This confirms results obtained more than 25 years ago [89]. This is probably due to the migration pattern with larger fish, with greater swimming efficiency, being concentrated along the northern coast where the southward and offshore currents are stronger. Such size segregation in the migration process has been suggested by biophysical modeling in West Africa [4546] and observed in other places [4748].

The condition factor is an indication of fish well-being and reproductive activity. This factor was evaluated over the whole year and showed just one trough in March and two peaks, one in April (during the upwelling period) and one in July (period of recovering from spawning). Cury and Fontana [6] suggested that the condition factors are higher in the cold season when upwelling induces plankton blooms. The condition factors were lower during the spawning period and during periods without upwelling.

The allometric scaling exponent reported in this study was similar to that recorded in other studies (Cabo Verde in 1987, Gambia in 2005, Mauritania in 1979 and Senegal in 1982 and 1986; see Table 5). The results also confirmed the measurements obtained in Senegal in 1980 [7] and 1988 [8].

thumbnail
Table 5. Length-weight parameters in West Africa for Sardinella maderensis.

The length is expressed as fork length.

https://doi.org/10.1371/journal.pone.0156143.t005

The overall sex ratio of S. maderensis obtained in this study was slightly female biased (55.81%), confirming several previous studies undertaken about 30 years ago in Senegal (Table 6). That was also the case in Cameroun and Mauritania (Table 6), even if locally in the northern coast of Senegal it was reported a sex ratio slightly male biased (~51%). Two different mechanisms have been proposed to explain this. On the one hand, for several fish species, warmer waters increase the female bias [4950]. This is in line with our observations in Senegal and with an even higher female bias being found in Cameroon (warm waters) but it does not explain the high female bias in Mauritania, where the temperature is usually lower. On the other hand, the “Petite Côte” is a well-known S. maderensis spawning area [7] and favorable area for fish larval development due to high phytoplankton biomass production and favorable environmental conditions, e.g. high coastal retention for eggs and larvae, in particular during the hot season [51, 52]. Thus in this area and this season, the female bias may be due to an evolutionary strategy giving a better reproduction efficiency for short-lived fish species. This mechanism may also explain the high female bias in Mauritania and Cameroon because the samples were collected in the local reproduction area [6, 7].

thumbnail
Table 6. Sardinella maderensis sex ratios based on the published literature for West Africa.

https://doi.org/10.1371/journal.pone.0156143.t006

The L50 values (16.60 cm for male and 17.50 cm for female) were in the same range as previous studies in West African aquatic ecosystems (Table 7). The L50 does not, therefore, appear to be influenced by geographical location [29, 40]. For Senegal, previous estimates of L50 were also similar, with the exception of [53] who reported a higher value (expressed in Total Length TL). There did not appear to be any reduction in L50 for S. maderensis for Senegal. However, a reduction in the L50 has been observed for S. aurita off the coast of Ghana [54].

thumbnail
Table 7. L50 for Sardinella maderensis based on the literature for West African areas.

L50 represents the size at which 50% of individuals in a given population attain maturity. (N/A means not available).

https://doi.org/10.1371/journal.pone.0156143.t007

This study found two spawning periods for S. maderensis, confirming previous reports [6]. The first spawning period was from April to October, with two spawning peaks, one in June–August and the other in October. The second spawning period occurred from early January to the end of February, with a single intense spawning peak. Several studies [79, 38, 5556] reported that S. maderensis almost spawned all year round in Senegal, with two spawning peaks, [one short (~ one month) and intense in cold season (February) and another longer (~four months) in hot season (June to October)], nevertheless we have to notice that it was reported a single spawning period in Mauritania by Pascual-Alayón et al. [57] from May to July.

In regard to our results S. maderensis has only limited plasticity to environmental variations. S. maderensis is less sensitive to climatic fluctuations and tolerates larger environmental fluctuations than S. aurita. S. maderensis, therefore, exhibits a smaller variability of biological and demographic parameters than S. aurita [6]. S. maderensis has a less flexible adaptive strategy resulting in a smaller plasticity of the biological parameters studied. This does not allow for better utilization of transitory resources but, on the other hand, the population can persist when conditions are unfavorable.

Previous studies have used the size at first sexual maturity to set the size at first capture for fisheries management [58, 59]. Based on our findings, we recommend that this should be set at 18 cm. The sexual maturity rest period of S. maderensis should be established in relation to the spawning periods, particularly during spawning from June to October. To draw up a fishery management plan, biological and ecological studies should be expanded to include Mauritania, Gambia and Guinea Bissau. In addition, socio-economic and behavioral aspects of fishing communities must be understood to draw up a management plan that would be acceptable.

The biological parameters of sardinella were expected to change with the effect of climate change and overfishing [13]. It was expected that either the effect of fishing pressure (which has increased considerably in Senegal over the last decades [3]) or the effect of the environment differences would have been detectable. Nevertheless the results of this study are consistent with those previously reported in the literature studied and suggest that S. maderensis has not been affected by overfishing or global climate change despite its habitat covering very different marine ecosystems including the Canary Current and the Guinea Current Large Marine Ecosystem. This narrow plasticity of the biological parameters studied should be taken into account in fishery management plans.

Supporting Information

S1 Fig. Wetherall [28] plot of Sardinella maderensis, leading to an estimate of L = 33.40 cm.

https://doi.org/10.1371/journal.pone.0156143.s001

(DOCX)

S2 Fig. K-scan routine for the size-frequency data of Sardinella maderensis in ELEFAN (Daniel Pauly, unpublish).

https://doi.org/10.1371/journal.pone.0156143.s002

(DOCX)

S3 Fig. Seasonalized von Bertalanffy growth curve (with L = 33.9 cm, K = 0.4 year-1, C = 1 year-1, WP = 0.15 and Rn = 0.21) of Sardinella maderensis off the northern coast superimposed over the restructured length-frequency histograms.

The black and grey bars are positive and negative deviations from the “weighted” moving average of two size classes and represent pseudo-cohorts. The red dot is the starting point through which the curve passes to fit the model by maximizing Rn (index of goodness of fit which is analogous, but not equivalent to r in linear regression (see Eq 4)).

https://doi.org/10.1371/journal.pone.0156143.s003

(DOCX)

S4 Fig. Seasonalized von Bertalanffy growth curve (with L = 32.1 cm, K = 0.38 year-1, C = 1 year-1, WP = 0.15 and Rn = 0.21) of Sardinella maderensis off the southern coast superimposed over the restructured length-frequency histograms.

The black and grey bars are positive and negative deviations from the “weighted” moving average of two size classes and represent pseudo-cohorts. The red dot is the starting point through which the curve passes to fit the model by maximizing Rn (index of goodness of fit which is analogous, but not equivalent to r in linear regression (see Eq 4)).

https://doi.org/10.1371/journal.pone.0156143.s004

(DOCX)

Author Contributions

Conceived and designed the experiments: KB MT PB TB OT. Performed the experiments: KB AF IN MT. Analyzed the data: KB OS. Contributed reagents/materials/analysis tools: OT AS MT. Wrote the paper: KB MT PB TB AF AS IN NL OS JP OT.

References

  1. 1. Sow F, Samb B, Cheng H. Diagnoses of the state of exploitation of small pelagic stocks in Senegal by the use of size frequencies: case Sardinella. Resources Sciences. 2010; 32: 620–626.
  2. 2. Anonymous. Résultats généraux des Pêches maritimes du Sénégal. Ministère de la Pêche et de l'Économie Maritime, Direction des Pêches Maritimes. 2010.
  3. 3. Food and Agriculture Organization of the United Nations. Report of the FAO Working Group on the Assessment of Small Pelagic Fish off Northwest Africa. Casablanca. Morocco, 24–28 May 2011. Rome: FAO Fisheries and Aquaculture. 2013; Available: http://www.fao.org/docrep/017/i3135b/i3135b.pdf
  4. 4. Dobo A. Illegal Chinese fishing in West African waters: A study on Chinese IUU activities and its consequences to socio-ecological systems. Master’s thesis, Stockholm University. 2009. Available: http://su.diva-portal.org/smash/record.jsf?pid=diva2%3A329241&dswid=-7159
  5. 5. Belhabiba D, Koutob V, Sall A, Lam VWY, Pauly D. Fisheries catch misreporting and its implications: The case of Senegal. Fisheries Research. 2014; 151: 1–11.
  6. 6. Cury P, Fontana A. Compétition et stratégies démographiques comparées de deux espèces de sardinelles (Sardinella aurita et Sardinella maderensis) des côtes ouest-africaines. Aquatic Living Resources. 1988; 1: 165–180.
  7. 7. Boëly T. Biologie de deux espèces de sardinelles (Sardinella aurita (Vallencienes, 1847) et Sardinella maderensis (Lowe, 1841) des côtes sénégalaises. Thèse de Doctorat d’Etat, Université de Paris VI. 1980. Available: http://horizon.documentation.ird.fr/exl-doc/pleins_textes/pleins_textes_5/b_fdi_02-03/03735.pdf
  8. 8. Camarena-Luhrs T. La croissance de Sardinella maderensis (Lowe, 1841) au Sénégal. Océanographie tropicale. 1986; 21: 143–151. ISSN 0245-9418
  9. 9. Fréon P. Réponses et adaptations des stocks de clupéidés d'Afrique de l'ouest à la variabilité du milieu et de l'exploitation: Analyse et réflexion à partir de l'exemple du Sénégal. PhD, Université d'Aix-Marseille. 1986. Available: http://horizon.documentation.ird.fr/exl-doc/pleins_textes/pleins_textes_2/etudes_theses/25601.pdf
  10. 10. Nisbet RM, Muller EB, Lika K, Kooijman SALM. From molecules to ecosystems through dynamic energy budget models. Journal of Animal Ecology. 2000; 69: 913–926.
  11. 11. Ernande B, Dieckmann U, Heino M. Adaptive changes in harvested populations: plasticity and evolution of age and size at maturation. Proceedings of the Royal Society of London B: Biological Sciences. 2004; 271: 415–423.
  12. 12. Thiaw M. Dynamique des ressources halieutiques à durée de vie courte: cas des stocks de poulpe et de crevettes exploités au Sénégal. Thèse de Doctorat mention Halieutique, Agrocampus Ouest Rennes, France. 2010. Available: www.halieutique.agrocampus-ouest.fr/pdf/943.pdf
  13. 13. Roy C, Cury P, Fontana A, Belvèze H. Stratégies spatio-temporelles de la reproduction des clupéidés des zones d'upwelling d'Afrique de l'Ouest. Aquatic Living Resources. 1989; 2: 21–29. http://dx.doi.org/10.1051/alr:1989003
  14. 14. Faure V. Environnement et variabilité des populations de poulpes Octopus vulgaris en Afrique de l'Ouest. In A. Caverivière A, Thiam M and Jouffre D. editors. Le poulpe Octopus vulgaris: Sénégal et côtes nord-ouest africaines. IRD Éditions, Paris. 2002; pp. 129–142.
  15. 15. Fontana A. Étude de la maturité sexuelle des sardinelles Sardinella eba (Val) et Sardinella aurita (c. et v.) de la région de Pointe-Noire. Cahier ORSTOM, Série Océanographie. 1969; 7: 101–113.
  16. 16. Le Cren ED. The length weight relationship and seasonal cycle in gonad weight and condition in perch (Perca fluvialitis). Journal of Animals and Ecology. 1951; 20: 201–219.
  17. 17. Somers IF. On a seasonally-oscillating growth function. Fishbyte (ICLARM). 1988; 6: 8–11.
  18. 18. Pauly D, Moreau J, Abad N. Comparison of age-structured and length-converted catch curves of brown trout Salmo trutta in two French rivers. Fisheries Research. 1995; 22: 197–204.
  19. 19. Etim L, Sankare Y. Growth and mortality, recruitment and yield of the fresh-water shrimp, Macrobrachium völlenhovenii, Herklots 1851 (Crustacea, Palaemonidae) in the Fahe reservoir, Côte d'Ivoire, West Africa. Fisheries Research. 1988; 38: 211–223.
  20. 20. Rueda M, Santos-Martínez A. Population dynamics of the striped mojarra Eugerres plumieri from the Ciénaga Grande de Santa Marta, Colombia. Fisheries Research. 1999. 42: 155–166.
  21. 21. Pauly D, Munro JL. Once more on the comparison of growth in fish and invertebrates. Fishbyte. 1984; 2: 1–21.
  22. 22. Piñeiro C, Saínza M. Age estimation, growth and maturity of the European hake (Merluccius merluccius (Linnaeus, 1758)) from Iberian Atlantic waters. ICES Journal of Marine Science. 2003; 60: 1086–1102.
  23. 23. Pauly D. Some simple methods for the assessment of tropical fish stocks. FAO Fisheries Technical Paper 1983; 234: 52.
  24. 24. Le Bec C. Cycle sexuel et fécondité de la sole Solea vulgaris (quensel, 1806) du golfe de Gascogne. Revues des Travaux de l’Institut des Pêches maritimes. 1985; 47: 181–187.
  25. 25. McDonald JH. Handbook of Biological Statistics. 2nd ed. Sparky House Publishing: Baltimore, Maryland; 2009.
  26. 26. Sokal RR, Rohlf FJ. Biometry: the principles and practice of statistics in biological research. 3rd ed. New York, USA: W. H. Freeman; 1995.
  27. 27. Zar JH. Biostatistical analysis. 4th ed. Upper Saddle River: Prentice Hall, New Jersey; 1999.
  28. 28. Wetherall JA. A new method for estimating growth and mortality parameters from length-frequency data. Fishbyte (ICLARM). 1986; 4: 12–14.
  29. 29. Samb B. Seasonal growth, mortality and recruitment pattern of Sardinella maderensis off Senegal. FAO Fisheries Report. 1988; 389: 257–271.
  30. 30. FAO. Contributions to tropical fisheries biology. FAO Fisheries Report. 1988. Available: http://www.fao.org/docrep/014/s8200e/s8200e.pdf
  31. 31. Diouf PS. Les peuplements de poisson des milieux estuariens de l’Afrique de l’Ouest: l’exemple de l’estuaire hyper salin du Sine-Saloum. Thèse de Doctorat, Université de Montpellier II. 1996 Available: http://horizon.documentation.ird.fr/exl-doc/pleins_textes/pleins_textes_7/TDM_7/010008130.pdf
  32. 32. FAO. Rapport du séminaire pour l’unification des techniques de lecture d’âge et d’interprétation de Sardina pilchardus (Walb.). COPACE/TECH 79/18: 32. 1979.
  33. 33. Lawal HS, Mylnikov N. Contribution à l'étude de la relation taille-poids, de la reproduction et du sex-ratio des principales espèces pélagiques en Mauritanie. Groupe de Travail, 88 p. 1988.
  34. 34. FAO. Rapport du troisième groupe de travail ad hoc sur l'évaluation des stocks de céphalopodes dans la région nord du COPACE. COPACE/PACE Séries 86/41. 1987;101. Available: http://www.fao.org/docrep/003/x6836f/X6836F00.htm#TOC
  35. 35. Écoutin JM, Albaret JJ, Trape S. Length-weight relationships for fish populations of a relatively undistributed tropical estuary: the Gambia. Fisheries Research. 2005; 72: 347–351.
  36. 36. Showers PAT. Comparative growth performance for species of the family Clupeidae of Sierra Leone. Naga ICLARM Q. 1996; 19: 42–44.
  37. 37. Marchal É. Biologie et écologie des poissons pélagiques côtiers du littoral ivoirien. In Environnement et ressources aquatiques de Côte-d’Ivoire. Tome I—Le milieu marin. Le Loeuff P. Marcha1 E, Kothias J-BA editors. Editions de l’ORSTOM, Paris. 1993; pp. 237–269.
  38. 38. Djama T, Gabche C, Njifonju O. Growth of Sardinella maderensis in the Lobe Estuary, Cameroon. Fishbyte.1989; 7: 8–10.
  39. 39. Gabche C, Hockey H-UP. Growth, mortality and reproduction of Sardinella maderensis (Lowe, 1841) in the artisanal fisheries off Kribi, Cameroon. Fisheries Research. 1995; 24: 331–344.
  40. 40. Djama T, Gabche C, Youmbi-Tienctheu J. Comparisons of the growth of West African Stock of Sardinella maderensis with Emphasis on Cameroon. ICLARM Fishbyte. 1989; 7: 13–14.
  41. 41. Marcus O. Breeding, age and growth in Sardinella maderensis (Lowe 1839) Pisces: Clupeidae from coastal waters around Lagos, Nigeria. Journal of Science. 1989; 23: 1–5.
  42. 42. Gheno Y, Le Guen JC. Détermination de l'âge et croissance de Sardinella eba (Val.) dans la région de Pointe-Noire. Cahier ORSTOM, série Océanographie. 1968; 6: 69–82.
  43. 43. Postel E. Résumé des connaissances acquises sur les clupéidés de l'ouest africain. Rapport du Conseil d’Exploration de La Mer. 1955; 137: 7.
  44. 44. Rossignol N. Première observations sur la biologie des Sardinelles dans la région de Pointe Noire (Sardinella eba, Val., S. aurita Val). Rapport du Conseil d’Exploration de La Mer. 1955; 137: 17–21.
  45. 45. Youmbi JT, Djama T, Gabche C. Reproductive patterns of Sardinella maderensis (Lowe, 1841) off Cameroon. Journal of Applied Ichthyology. 1991; 7: 60–63.
  46. 46. Brochier T, Auger PA, Thiaw M, Pecquerie L, Machu E, Mbaye BC et al. Investigating Sardinella aurita seasonal migratory pattern off North-West Africa with a biophysical model. Book of abstract IV Congress of Marine Sciences (ISBN—10: 84-697-0471-0), pp 458–458, June 11th to 13th 2014, Las Palmas de Gran Canaria, Spain. 2014.
  47. 47. Nøttestad L, Giske , Holst JC, Huse G. A length-based hypothesis for feeding migrations in pelagic fish. Canadian Journal of Fisheries Aquatics Science. 1999; 56: 26–34.
  48. 48. Lingen CD van der, Coetzee JC, Hutchings L. Overview of the KwaZulu-Natal sardine run. African Journal of Marine Science. 2010; 32: 271–277.
  49. 49. Ospina-Álvarez N, Piferrer F. Temperature-Dependent Sex Determination in Fish Revisited: Prevalence, a Single Sex Ratio Response Pattern, and Possible Effects of Climate Change. PLoS ONE. 2008; 3: 28–37.
  50. 50. Valenzuela N, Adams DC, Janzen FJ. Pattern Does Not Equal Process: Exactly When Is Sex Environmentally Determined? The American Naturalist. 2003; 161: 676–683. pmid:12776892
  51. 51. Demarcq H, Faure V. Coastal upwelling and associated retention indices derived from satellite SST. Application to Octopus vulgaris recruitment. Oceanologica Acta. 2000; 23: 391–408.
  52. 52. Mbaye BC, Brochier T, Echevin V, Lazar A, Lévy M, Mason E, et al. Do Sardinella aurita spawning seasons match local retention patterns in the Senegalese–Mauritanian upwelling region? Fisheries Oceanography 2015; 24: 69–89.
  53. 53. FAO. Science and Management of Small Pelagics. Symposium on Science and the Challenge of Managing Small Pelagic Fisheries on Shared Stocks in Northwest Africa, 11–14 March 2008, Casablanca, Morocco. Garcia S; Tandstad M and Caramelo AM editors. FAO Fisheries and Aquaculture Proceedings No. 18, 2012. Available: http://www.fao.org/docrep/017/i2523b/i2523b.pdf
  54. 54. Quaatey SNK, Maravelias CD. Maturity and spawning pattern of Sardinella aurita in relation to water temperature and zooplankton abundance off Ghana, West Africa. Journal of Applied Ichthyology. 1999. 15: 63–69.
  55. 55. Conand F, Fagetti E. Description et distribution saisonnière des larves de sardinelles des côtes du Sénégal et de la Gambie en 1968 et 1969, Cahier ORSTOM, Série Océanographie. 1971; 10(3): 293–318.
  56. 56. Kromer JL. Rio Grande de Buba: Bio-écologie et paramètres environnementaux. UICN/Ministère des pêches de Guinée-Bissau. 1994; 119p.
  57. 57. Pascual-Alayón P, Sancho A, Hernández E, Santamaría MTG, Duque V, Balguerías E et al. Reproductive aspects of sardine (Sardina pilchardus Walbaum, 1792), round sardinella (Sardinella aurita Valenciennes, 1847), flat sardinella (Sardinella maderensis Lowe, 1838) and mackerel (Scomber colias Gmelin, 1789) off Mauritanian waters, Northwest Africa. In S. Garcia, M. Tandstad and A.M. Caramelo Editors. Science and Management of Small Pelagics. Symposium on Science and the Challenge of Managing Small Pelagic Fisheries on Shared Stocks in Northwest Africa. 2012 Available:ftp://ftp.fao.org/fi/Cdrom/P18_CDROM/root%5Cnew%5CSYMPOSIUM_PART2THEME1-POSTERSPAPERS/P1.3_PASCUAL-ALAYON.pdf
  58. 58. Froese R. Keep it simple: three indicators to deal with overfishing. Fish and Fisheries. 2004; 5: 86–91.
  59. 59. Froese R, Binohlan C. Empirical relationships to estimate asymptotic length, length at first maturity, and length at maximum yield per recruit in fishes, with a simple method to evaluate length frequency data. Journal of Fish Biology. 2000; 56: 758–773.