© 2006 International Council for the Exploration of the Sea
Temporal and geographic variability of sardine maturity at length in the northeastern Atlantic and the western Mediterranean
a Instituto de Investigação das Pescas e do Mar (INIAP-IPIMAR), Av. Brasília s/n, 1449-006 Lisboa, Portugal
b Instituto Español de Oceanografía (IEO), Centro Oceanográfico de Vigo, Apdo. 1552, 36280 Vigo, Spain
c Museo do Mar de Galicia, Avenida Atlántida 160, 36280 Vigo, Spain
*Correspondence to A. Silva: tel: +351 213027000; fax: +351 213025948. e-mail: asilva{at}ipimar.pt.
We investigate spatial and temporal variations in sardine maturation patterns and discuss the implications for stock structure, using recent samples collected opportunistically across the species range and data collected regularly for the assessment of the Atlanto-Iberian stock. Maturity ogives were fitted to maturity-at-length data using Generalized Linear Models, and maturation patterns are summarized by length at first maturity (L50). Sardine attain sexual maturity at a total length of 10.916.8 cm. The percentage of mature small sardine increases from north to south across the Northeast Atlantic. L50 is low in the Mediterranean and Moroccan waters, but maturation is delayed towards the limits of the species' range. Sardine mature smaller in the Gulf of Cadiz than in other areas of the Atlanto-Iberian stock, although there is no indication of a different maturation pattern north of the stock boundary. These results are not sufficient to review the stock structure of sardine, but question the hypothesis of homogeneous biological properties of sardine populations within the stock area. Sardine maturation length reduced and the spawning period extended in western Iberia during the past two decades. Associated increases in fish condition and fat reserves during the summer feeding season suggest that variations in reproductive traits may have been caused by environmentally driven changes in food availability.
Keywords: condition, Iberian waters, maturity at length, sardine, stock structure
Received 15 August 2005; accepted 5 January 2006.
| Introduction |
|---|
|
|
|---|
Spatial and temporal trends in maturation have been documented for many fish species. On wide geographic scales, maturation tends to co-vary with growth in such a manner that the ratios between length/age at first maturity and maximum length at age remain approximately stable across the distribution range of species (Beverton, 1992). Populations inhabiting colder habitats generally grow more slowly, live longer, and delay maturation when compared with those from warmer areas of the distribution, the differences reflecting the ability of the species to adapt to large-scale patterns in environmental conditions (Beverton, 1992). As is the case for other life history traits, maturation expresses the interaction between the genetic background and environmental influences, and provides evidence of geographic and/or reproductive isolation of fish populations (Begg, 2005). Therefore, differences in maturation have been taken as evidence for population discreteness within holistic approaches for the identification of fish stocks (Begg and Waldman, 1999; Begg et al., 1999; Berg and Albert, 2003). Apart from their value in stock identification, spatial variations in maturation within the boundaries of a stock should be taken into account in stock assessment if they are associated with large spatial differences in population abundance (Hilborn and Walters, 1992).
Temporal trends in fish maturation have often been interpreted as compensatory density-dependent effects that regulate population growth (Rose et al., 2001). Density-dependent changes in maturation arise from food limitation attributable to increased intraspecific competition. The length or the age of maturation may be directly affected through changes in the quantity of energy reserves available for gonad development (Morgan, 2004), or indirectly via changes in growth which in turn influence the triggering of maturation (Engelhard and Heino, 2004). Density-dependent effects may explain the decrease in length/age at first maturity in many fish populations that have shown extensive declines in biomass (Armstrong et al., 1989; Nash et al., 2000; Engelhard and Heino, 2004). As most of these populations belong to heavily exploited fish stocks, fishing mortality has been suggested as the selective pressure inducing a rapid evolution of maturation characteristics (Barot et al., 2005). However, it has often proved difficult to disentangle the direct influence of population density from that of long-term environmental trends, and this is further confounded by the effects of environmental conditions on population abundance, which in turn can trigger density-dependent maturation (Rose et al., 2001).
The sardine (Sardina pilchardus) is a commercially important clupeoid distributed in the northeastern Atlantic from the North Sea to Senegal and throughout most of the Mediterranean Sea (Parrish et al., 1989). Its largest populations and the more productive sardine fisheries are situated off Morocco (stock biomass 15 million tonnes; catches around 600 000 t in recent years; FAO, 2003), whereas abundance and fishing intensity are low in the Mediterranean (http://www.faocopemed.org). Within Atlantic European waters, the sardine is considered to belong to a single stock for assessment purposes, the Atlanto-Iberian stock, with biomass around 500 000 t and catches around 100 000 t in recent years (ICES, 2005).
The Atlanto-Iberian sardine stock is delimited by the French/Spanish border in the north, and by the Strait of Gibraltar in the south (ICES Divisions VIIIc and IXa). It is exploited mainly by the purse-seine fleets of Spain and Portugal, and it has been assessed since the early 1980s assuming that it consists essentially of self-reproducing populations that share homogeneous life history properties. However, recent studies provide evidence of geographic differences in morphometric characters, age structure, and population dynamics within the stock area, which challenge the assumption of a single management unit (ICES, 2000; Carrera and Porteiro, 2003; Silva, 2003). The stock boundaries have also been questioned on the basis of evidence of a continuous distribution of sardine eggs and adults extending into French waters and of possible mixing among populations across the southern border of the stock (ICES, 2000; Stratoudakis et al., 2003). Uncertainty regarding the definition of the stock unit has affected the level of confidence in stock estimates and knowledge of long-term trends in abundance, and has highlighted the need for biological information that can both assist in defining management units and provide data to pursue stock assessment at finer spatial scales (ICES, 2005).
The main goal of this manuscript is to study the spatial and temporal variability of sardine maturation patterns and to discuss their contribution to the definition of stock structure. Data collected recently across most of the area of distribution of the species are used to describe large-scale geographic variability in life history traits such as longevity, maximum length, and maturity at length. Despite the narrow temporal scale and the limited sampling intensity, the samples provide an opportunity to compare sardine maturation over a wide geographic area, using data collected and analysed with similar methodology. Data on length, weight, sex, maturity, fat content, and age routinely collected by Spain and Portugal since the early 1980s for the assessment of the Atlanto-Iberian sardine stock are used to compare sardine maturation patterns among stock areas, and to assess the temporal stability of such patterns by following interannual changes in the length at first maturity (L50). Data on population abundance (from landings and estimated biomass) of the Atlanto-Iberian sardine stock are used to explore the existence of potential relationships with size at maturity (positive if both parameters represent a response to environmental conditions; negative if there is density-dependent competition). Finally, fish condition (weight and fat content at length) is used to explore the relationship between the accumulation of energy reserves during summer feeding and the size at first maturity.
| Material and methods |
|---|
|
|
|---|
Sardine samples collected within the remit of various international projects during the period 19992004, mainly for studying stock structure, are used to describe broad geographical patterns in sardine maturation. Overall, 51 samples of sardine (generally 100 fish) collected from the north (English Channel to southern France), south (northern Morocco to Mauritania), east (Azores and Madeira), and west (western Mediterranean) of the distribution of the Atlanto-Iberian stock were used (Figure 1). All samples were collected within the spawning season, either in the ports (five samples) or during research surveys, but the final distribution of samples among areas and years is very unbalanced (e.g. 28 samples in south France and only one sample in Mauritanian waters). Samples were transferred frozen to the laboratory of IPIMAR, where fish were thawed and processed following routine biological sampling (see below).
|
Biological samples from sardine commercial landings and acoustic surveys have been regularly collected from the area since the early 1980s, and are used to describe temporal trends and small-scale spatial variation in the maturation pattern of the stock. Samples from landings by purse-seiners were collected once or twice a month between 1981 and 2003 at the main ports of the Iberian Peninsula (Figure 1). Sampling started more recently and has been less regular in the Cantabrian Sea, North Galicia (1993), and Cadiz (2000), providing around 800 fish by port and year. Ports located in southern Galicia and in northern, southwestern, and southern Portugal account on average for 80% of the total landings of the stock. Sampling regularity in these ports was generally good, with an average of 1200 fish collected per port per year, except during periods of low sardine abundance or seasonal closures of the fishery. Each sample consists of six fish per half-centimetre class (Spanish ports) or a random collection of 100 fish (Portuguese ports). Samples from acoustic surveys have been collected irregularly since 1984 in spring (Spanish and Portuguese waters) and autumn (Portuguese waters). During such surveys, random collections of 40 fish (Spanish survey) or length-stratified samples of 1015 fish per half-centimetre length class (Portuguese surveys) were taken from each fishing haul.
Total length (cm), macroscopic maturity, and visceral fat stage were recorded for all fish. Total weight (g) was recorded for all market samples, but gutted weight (Wgut) and gonad weight (Wgon) are only available since the mid-1980s in sardine samples collected at Portuguese ports. Age data from survey samples collected within the stock area and from opportunistic samples collected outside this area are used to summarize the age structure of sardine populations across the study area. Macroscopic maturity stage was determined by visual examination of the gonads using a six-stage key: 1, immature; 2, developing; 3, pre-spawning; 4, spawning; 5, spawning/recovering; 6, post-spawning (Pinto, 1957; Pinto and Andreu, 1957). Fish at maturity stage 2 and above were considered to be mature (here used as a synonym for "adults", i.e. potentially part of the spawning population) and those in stages 35 were assumed to be showing spawning activity. The classification of the fat stage (1, no fat visible; 2, thin thread of fat surrounding part of the gut; 3, thicker layer of fat partially surrounding the gut; 4, thick layer of fat surrounding the gut) was based on the key proposed by Furnestin (1943).
Length-based maturity ogives were obtained from samples grouped into areas roughly according to ICES or FAO Divisions. For a broad geographic comparison of maturity, ogives were obtained from interannually pooled samples from the areas outside the stock and from pooled survey data (separately for the 1980s and the 1990s2000s) for areas within the stock region. To explore temporal trends in maturity within the stock area, maturity ogives were estimated from market samples pooled for each spawning season (OctoberJanuary, spanning the first part of the spawning season up to the peak). The sampling period used to estimate maturity ogives was based on an analysis of sardine spawning seasonality as a function of fish length in the 1980s and 1990s. The probability of spawning (sp) was modelled as an anisotropic bi-variate function (f(.)) of fish length class (l) 1218 cm and month (mo) (SeptemberAugust, coded with numbers 112) within each area and period using Generalized Additive Models with a binomial error distribution and a logit link (Wood, 2000). The general form of the models is
|
|
Generalized Linear Models (GLMs; McCullagh and Nelder, 1989) with a binomial error distribution and a logit link were fitted to the proportion of fish mature (m) by length class (l). The model has the general form
|
|
|
|
|
|
The standard error of L50 was estimated using the formula for the approximate variance of a ratio of two random variables (Mood et al., 1974):
|
|
|
|
Standard errors of the parameters were corrected to account for the inflation of variance (relative to the nominal binomial variance) attributable to the combination of data from different samples (see below). This correction was carried out by multiplying the parameter standard errors by the square-root of the replication estimate of overdispersion (McCullagh and Nelder, 1989; Millar et al., 2004). Asymptotic confidence intervals were calculated using the t-statistic. Geographical differences and temporal trends in L50 within western Iberian waters were tested for the period 19852003 by linear regression analysis. The effects of area (factor variable), year (continuous variable), and their interaction were tested by F-test using a forward procedure (McCullagh and Nelder, 1989). The inverse standard errors of L50 were used as weighting factors. The northern and southern Iberian areas were not included in this analysis owing to a lack of data. All calculations were carried out with R 1.9.1 (The R Development Core Team, 2004).
| Results |
|---|
|
|
|---|
Large-scale spatial variation in life history parameters
Sardine total length in all samples used in this study ranged between 7.0 and 29.5 cm, corresponding to ages of 013 years (Table 1). Length and age range decreased from north to south in the northeastern Atlantic and were generally smaller in the Mediterranean Sea. Cape Finisterre at the northwestern corner of the Iberian Peninsula seems to form a boundary in the length and age distribution of samples from the Atlantic, with higher median lengths and ages in areas to the north (
20 cm and 3 years) than to the south (
17 cm and 1 year old). The sample from Mauritania was an exception to this trend, showing a similar length distribution to northern areas but a reduced age range. Samples from the western Mediterranean and from the Atlantic islands of the Azores and Madeira had smaller length and age ranges than those from southern Iberia.
|
Table 2 summarizes information on sardine maturation by area. Point estimates of L50 varied between 10.9 and 16.8 cm (mean 14.3 cm, CV 10%), and slopes in L50 ranged from 0.13 to 0.53 cm1 (mean 0.30 cm1, CV 42%). There is some evidence of an overall increase in the percentage of small sardines mature from northern France to the Gulf of Cadiz but, unlike for length and age distributions, there is no clear boundary at Cape Finisterre. Instead, L50 is considerably higher in northern France (16.8 cm) but comparable among areas to the south of this region. At the southern extent of the Iberian Peninsula (Gulf of Cadiz), L50 was 2.4 cm lower than the overall mean, and transition to maturity is clearly steeper than in more northerly areas (slope in L50 0.52 cm1) in recent years. The area of transition from a higher L50 in the west to a lower L50 near Cadiz is not clear because maturity parameters from southern Portugal were closer to those from Cadiz in recent years but similar to the western Iberian areas during the 1980s. Therefore, the most evident discontinuity in the maturation pattern of sardine within the Atlanto-Iberian stock area occurred in the Gulf of Cadiz.
|
Maturity ogives could not be estimated for the western Channel, the Mediterranean, or North Moroccan areas, because immature fish were either absent or present in only small numbers within the samples. This would be expected for the Channel where the minimum length of fish sampled was 18 cm. Conversely, the small fish (1213.5 cm) in the Mediterranean and Moroccan samples were mostly adults, providing some indication that L50 values in those areas would be at least as low as those in northwestern Iberian areas. The samples from Mauritania (L50 16.2 cm) and Azores (L50 15.7 cm) suggest that sardine length at first maturity may increase towards the southern and the western limits of distribution of the species.
Spatio-temporal variation in maturity at length along the Iberian Peninsula
Annual maturity ogives are seldom defined for the northern stock areas (Cantabria and north Galicia) because samples are dominated by large fish (median length
21 cm, median age
4 years; Table 1). Pooled survey data provided estimates of L50 between 13.3 and 15.6 cm, within the range of variation of L50 in the more southerly areas (Table 2). Similarly, most of the sardine sampled in the Gulf of Cadiz were mature, in both survey (Table 2) and market samples. However, unlike in northern Spanish waters, small sardine were found in the area during all spawning seasons. Four maturity ogives were defined based on survey data from the Gulf of Cadiz, providing the lowest pooled estimate of L50 within the stock area (see above).
Off western and southern Iberia, maturity-at-length data covered both juvenile and adult fish, defining clear maturity ogives in most spawning seasons. Patterns of spawning seasonality suggested no major differences among size classes, areas, or periods that could hamper the comparison of maturity ogives (Figure 2). In general, spawning activity extended from late autumn to late spring, with a peak between December and March. Smaller fish have shorter spawning seasons, but maximum activity of the different size classes was synchronized. In addition, there was a pronounced increase in the duration of the season and in the intensity of peak spawning in all areas from the 1980s to the 1990s, particularly for smaller fish. Acoustic survey data (not shown) support the above changes in spawning seasonality between decades. Although the spawning season expanded between decades, neither the peak nor the overall pattern of activity showed marked temporal shifts. The above results confirm that samples used to construct maturity ogives generally cover the progress of the spawning season up to the peak of spawning (January) of fish of different size in each area and period. In two cases (south Galicia and south Portugal during the 1980s), the spawning season shifted slightly towards spring, possibly introducing a positive bias in L50.
|
Figure 3 shows annual estimates of L50 based on GLM parameters for the spawning seasons 19812003 off western and southern Iberia (from south Galicia to south Portugal). L50 was not significantly different among the three western areas (p = 0.18). Overall maturation patterns were comparable in the four areas, as indicated by the similarity among the lengths of 25% (
12.0 cm), 50% (
14 cm), and 75% (
15 cm) maturity, based on data pooled for the study period. Figure 3 also suggests that L50 decreased in the four areas; this is corroborated by decadal survey data (Table 2). In fact, when the period from 1985 is considered, there is a significant linear decline in L50 (p < 0.001), and differences in this trend are not significant among the three western areas (p = 0.85). Downward shifts of L50 were observed also for the northern areas of the stock (Cantabrian Sea and north Galicia), resulting in a decline in the average length at first maturity for the Atlanto-Iberian stock of sardine, from 14.5 cm in the 1980s to 13.3 cm in the 1990s.
|
Annual changes are similar in southern Galicia and northern Portugal, where L50 rose from the early to the mid-1980s and decreased towards the late 1990s. Off southern Galicia, there was an associated increase in the slope of the maturity ogives (not shown). Off southwestern Portugal, L50 declined after a period of high values during the mid-1980s. The series of L50 for southern Portugal has large temporal gaps owing to the lack of small immature fish in samples during most of the 1990s, precluding analysis of temporal trends. For the early 1980s, L50s are within the range of variation of values for the other stock areas. However, there is indication mainly from survey data that L50 was significantly lower in the late 1990s (Table 2). Data on gonad weight provide some support to maturity data, given the significant negative correlation between L50 and the gonadosomatic index of maturing fish (in the range L25L75) in western Portuguese areas (Figure 4).
|
Relationship between L50, condition, and abundance
The influence of fish body condition on maturation was explored using data from the western Iberian areas on the mean total weight and the percentage of fish with medium/high visceral fat (fat stage >2; % fat) in summer (July and August). This is when plankton production and sardine body condition peak. Therefore, fish weight and fat reserves in summer are a reliable index of the energy made available for reproduction in the current spawning season (Koops et al., 2004), and possibly also for the spawning season of the following year (Hunter and Leong, 1981). As the mean weight and the % fat showed similar temporal variations in the different length classes, 18-cm fish (the best sampled length class during the study period) were used, assuming to be representative of the average condition of the whole population. The two variables showed an increasing trend through the period in all areas: pooled data for 19821985 and 20002003 showed increases of 811% in mean weight and of 2658% in the percentage of fish with fat. L50 was significantly negatively related to the mean weight-at-length (Spearman rho = 0.35, p = 0.01) and to the percentage of fish with fat (Spearman rho = 0.47, p = 0.005) in the current summer season. However, the correlation improved when L50 was related to fish weight (Spearman rho = 0.51, p = 0.0003) and % fat (Spearman rho = 0.61, p = 0.0003) in the previous summer, indicating that condition has a prolonged effect on the length at first maturity (Figure 5).
|
The influence of population abundance on the length at first maturity was also studied, using data on sardine landings and biomass (from acoustic surveys) off western Iberia. The long-term trends in maturation do not show an association with local annual landings or survey abundance indices in Portuguese waters. Such a result was expected considering there were no major changes in sardine abundance within those areas during the study period. However, off southern Galicia L50 was significantly correlated with landings (r = 0.48, p = 0.02; Figure 6). Landings are considered to be representative of the large decrease in sardine abundance in this area (Table 1; Carrera and Porteiro, 2003), so the above relationship provides some evidence of a density-dependent effect on maturation. There is also some indication that population abundance may influence maturation on a shorter time scale, because recent shifts in sardine abundance have been followed by wide variations in L50. In fact, owing to the strong incoming 2000 year class, sardine abundance increased by a full order of magnitude in northern Portugal from 1999 to 2000 (ICES, 2002). Maturation length rose 2.0 cm in the area and increased also in the adjacent areas (south Galicia and southwest Portugal) the following spawning season, to where the cohort spread during its first year of life (Figure 3).
|
| Discussion |
|---|
|
|
|---|
All data analysed in this study were collected using the same macroscopic criteria for maturity, so the only concern with data quality relates to the accuracy of macroscopic staging. The reliability of macroscopic maturity data is supported by the significant correlation between L50 and the gonadosomatic index of maturing individuals observed since 1987 off western Portugal (Figure 4). GSI data further strengthen the hypothesis of a decline in the length at first maturity in the area. Additional information on the accuracy of maturity staging is limited to a few recent intercalibration exercises, which point to some uncertainties in maturity classification. A recent intercalibration of macroscopic staging of sardine ovaries with histological criteria showed that 24% of stage 1 gonads were misclassified (they were in fact stage 2), and 2736% of post-spawning fish were wrongly identified as pre-spawning (early spawning season) or spawning/recovering (late spawning season; ICES, 2003). These results indicate that both the proportion of mature fish and the proportion of active individuals may have been slightly underestimated in the past.
The samples analysed across the area of distribution of sardine suggest an increase in the percentage of small sardine that are mature from north to south within the Northeast Atlantic. There is also evidence of earlier maturation in the Gulf of Cadiz, and in Mediterranean and Moroccan waters and, possibly, delayed maturation towards the northern (north France), southern (Mauritania), and western (Azores) limits of distribution (Table 2). Earlier results on sardine maturation support the broad spatial trends in maturity described here, corroborating the upward shift in L50 from southern to northern France, and the low length at first maturity of sardine off Cadiz compared with northwestern Iberian areas (Table 3). Further, they support the increase in maturation length towards the southern (Saharan Morocco) and western (Canary Islands) extent of sardine range. Large-scale spatial variation of sardine length at maturity shows an overall positive relationship with longevity and maximum fish size (Tables 1 and 2) and bears some direct association with growth during the first year (Andreu and Fuster de Plaza, 1962; Morales-Nin and Pertierra, 1990; FAO, 2001). A similar geographic pattern of maturation has been reported for horse mackerel (Trachurus trachurus) (Abaunza et al., 2003): L50 varies from 15 to 23 cm in north African waters and 16 cm in northwestern Mediterranean to 1619 cm off southern Portugal, and increases gradually to the northern areas, attaining 2025 cm in the North Sea and English Channel. Length at maturity is also higher in sardine populations with greater maximum size. These broad spatial patterns suggest that maturation has an inverse relationship with the potential growth span of the species across its range, conforming to the general principles of life history strategy (Beverton, 1992; Rose et al., 2001).
|
In Atlantic waters, the most striking discontinuities in the length at first maturity of sardine are at the southern (south Morocco/Mauritania) and western (Atlantic Islands) limits of distribution, and also between the southern and the northern French waters. At the southern limit of the species, there is some evidence of differences in morphological characters, age structure, spawning areas, growth, and maturation between north/central Morocco and south Morocco/Mauritania, pointing to the existence of separate sardine populations (FAO, 2001). Additionally, ongoing studies on sardine genetics indicate that sardine from the Azores and Madeira (western) and Mauritania (southern) are distinct from those in the Atlantic between southern Iberia and the English Channel (Kasapidis et al., 2004). By contrast, the abrupt increase in L50 between southern and northern French waters is neither explained by genetic differentiation (Kasapidis et al., 2004) nor by major discontinuities in the distribution of sardine between the areas (Bellier et al., 2004; ICES, 2005). The predominance of large mature sardine in northern French waters and in the English Channel (Furnestin, 1943; Hickling, 1945) supports the hypothesis that nursery grounds of these sardine populations are located in the south, possibly in the Bay of Biscay (Furnestin, 1943). The French waters north of the Gironde Estuary (French coast, 45°30'N) are distinctive oceanographically from those of the Bay of Biscay, with considerably lower temperature attributable to a persistent cold-water mass (Puillat et al., 2004). Colder water conditions also prevail in the northern areas of the English Channel and North Sea, restricting sardine spawning activity to a few months during summer and autumn (Stratoudakis et al., 2004). Such conditions may influence the migration of juveniles to southern areas, and also constrain the reproductive development of the few small fish that remain in the areas.
Atlanto-Iberian sardine stock maturation is comparable among the western areas of distribution, the mean L25, L50, and L75 being around 12, 14, and 15 cm, respectively. There is limited information on maturation length in northern Spanish waters (north Galicia and Cantabria), which is dominated by large adults, possibly immigrants from adjacent recruitment areas off northwestern Iberia or the Bay of Biscay (Carrera and Porteiro, 2003). Nevertheless, decadal L50s for these areas are similar to those for western Iberia. The long-term average maturation pattern of sardine off southern Portugal is uncertain because of the scarcity of data, but recent information points to similarity with sardine maturation in the Gulf of Cadiz. In the latter area, sardine mature smaller than the general pattern, which in turn is similar to that observed in the southwestern Mediterranean. A similar spatial pattern was demonstrated for sardine morphometry, the similarity extending also to northern Morocco (Silva, 2003). However, earlier studies on meristic characters distinguished European Atlantic from Northwest African sardine populations (Parrish et al., 1989), whereas Mediterranean and Northwest African sardine have been considered a different subspecies from the European Atlantic sardine (Bauchot and Pras, 1980). The lack of genetic distinction among sardine populations around the Iberian Peninsula (Kasapidis et al., 2004) suggests that maturation differences are mainly driven by environmental conditions common to both Mediterranean and Atlantic waters of southern Iberia, which are themselves distinct from the prevailing conditions off western Iberia. Differences in maturation may be mediated by growth, because lower mean lengths at age have been observed in both the Gulf of Cadiz and Mediterranean waters than off Galicia and central Portugal (Rodríguez-Roda, 1970).
The long-term samples off western Iberia show an increasing trend in the percentage of small sardines maturing each spawning season since the early 1980s, along with associated increases in spawning activity (mainly of smaller fish) and in the duration of the spawning season between the 1980s and the 1990s. These changes suggest an improvement in the overall reproductive potential of sardine, and may be the result of enhanced fish condition during that period. In fact, for other clupeoid species, most of the energy spent in a specific spawning season comes from fat reserves stored during the current and the previous season (Hunter and Leong, 1981). This hypothesis is supported here by the significant correlation between mean sardine weight and % fat at length with the length at first maturity.
Temporal changes in sardine condition may have been caused by environmentally driven variation in food availability. There is evidence of environmental change off western Iberia during the last two decades (Santos et al., 2001, 2005; Borges et al., 2003; Guisande et al., 2004), and these seem to be part of wider temporal and spatial variations in the hydroclimatic conditions of the whole North Atlantic (Ottersen et al., 2001). Based on satellite-derived sea surface temperature data, Santos et al. (2005) document a shift in the upwelling regime off western Iberia, from a weak intensity in the 1980s to stronger in the 1990s, during both winter and summer. This increased upwelling intensity may have enhanced plankton productivity in the area and hence food availability for sardine and other clupeoid populations. Sardine condition and therefore L50 may also be partially regulated by density dependence, as observed for other pelagic species (Parrish and Mallicoate, 1995). The fact that the long-term decline in L50 was inversely related to abundance off southern Galicia, an area that experienced a great decrease in its sardine population during the last two decades (Carrera and Porteiro, 2003), strengthens this possibility. It is also likely that variations in both sardine condition and abundance were affected by changes in environmental conditions observed in the past two decades off western Iberia.
In spite of the above evidence of environmentally driven changes in the length at first maturity, the hypothesis of an evolutionary trend in sardine maturation is not excluded. In fact, the pronounced changes in the abundance and age structure of sardine populations in mainly the northern areas of the stock (Galician and Cantabrian waters; ICES, 2000) might have provided scope for the selection of early maturing phenotypes. It can be argued that selective changes would not likely take place over such a short time (<20 years), because most experimental studies suggest that 1020 generations are needed for selection to occur (Swain et al., 2005). However, only four generations were needed for the growth rate of Atlantic silverside (Menidia menidia), a short-lived pelagic species (Connover and Munch, 2002), to evolve. Sardine have moderate longevity (78 years), but they mature during the first two years of life, as do other short-lived pelagic species. Therefore, the possibility that the maturation of sardine evolved over just two decades warrants further investigation.
Maturation reaction norms (Heino et al., 2002; Barot et al., 2004) are a useful tool in exploring evolutionary trends in maturation, because they allow disentangling the influence of growth and mortality (assumed to reflect environmental variation) on maturation. Estimation of reaction norms for sardine maturation could complement these analyses, because variation in growth may explain the geographic and temporal changes in maturity at length described. Age data are available from commercial and survey length-stratified samples of sardine covering part of the area and period analysed here, and they can be used in future to pursue such an approach. Nevertheless, the estimation of reaction norms for sardine requires careful consideration in relation to sampling design, consistency of age determinations, and also to whether age and length distributions derived from commercial samples are representative of the whole populations.
Overall, there is clear evidence that sardine may attain sexual maturity at a wide range of lengths across its distribution area, varying from 11 cm to almost 17 cm (mainly age groups 0 and 1), in relation to broad patterns of growth and longevity. Within the Atlanto-Iberian stock, maturation patterns by length are comparable across the western areas (from southern Galicia to southwestern Portugal), and there is no evidence of pronounced differences towards the northern areas (northern Galicia and Cantabria). On the other hand, sardine mature smaller in the Gulf of Cadiz than in the remaining stock areas, and there is no indication of a different maturation pattern north of the stock boundary (southern French waters). As is the case with other life history traits, length at maturity does not provide direct evidence of stock structure, so these results are not sufficient to review the identity of the sardine stock in Atlantic European waters. However, they do assist in delineation of areas with distinct phenotypic groups of sardine, and this information can be utilized along with evidence from other approaches to stock identification for management purposes.
| Acknowledgements |
|---|
This study was developed within the framework of the EU Project SARDYN (Contract Q5RS-2002-000818). We thank Khallahi Ould Mohamed Fall (IMROP Mauritania), Jacques Masse (IFREMER Nantes), Jean-Louis Bigot (IFREMER Sète), Nick C. Halliday and the crew of the MBA/PML RVs "Sarsia" and "Plymouth Quest" (MBA, England), Dalila Carvalho (DRPM, Madeira, Portugal), Angela Canha (DOP, Azores, Portugal), and Malika Chlaida (INRH Morocco) for collaboration in the collection of samples (as part of SARDYN and the EU Study PELASSES 080/99), and E. Soares, A. Morais, and E. Peleteiro (IEO) for assistance in the laboratory. Recent port data from the Bay of Biscay were provided by the AZTI Foundation. The collection of regular biological samples was supported by National and European Biological Data Collection Programmes of IPIMAR, IEO, and AZTI. We thank R. Fryer for helpful assistance in the application of Generalized Linear Models, and R. H. Parrish, K. Ganias, and two anonymous referees for valuable suggestions on earlier drafts of the manuscript.
| References |
|---|
|
|
|---|
-
Abad R. and Giraldéz A. (1993) Reproducción, factor de condición y talla de primera madurez de la sardina, Sardina pilchardus, (Walb.), del litoral de Málaga, mar de Alborán (1989 a 1992). Boletín del Instituto Español de Oceanografía 8:145155.
Abaunza P., Gordo L., Karlou-Riga C., Murta A., Eltink A.T.G.W., García-Santamaría M.T., Zimmermann C., Hammer C., Lucio P., Iversen S.A., Molloy J., Gallo E. (2003) Growth and reproduction of horse mackerel, Trachurus trachurus (Carangidae). Reviews in Fish Biology and Fisheries 13:2761.[CrossRef][Web of Science]
Andreu B. (1955) The sexuality of sardines. Proceedings and Technical Paper of the C.G.P.M.FAO, Rome 3: pp. 4560.
Andreu B. and Fuster de Plaza M.L. (1962) Estudio de la edad y crecimiento de la sardine (Sardina pilchardus Walb.) del NW de España. Investigaciones Pesqueras 21:4995.
Armstrong M.J., Roel B.A., Prosch R.M. (1989) Long-term trends in patterns of maturity in the southern Benguela pilchard population: evidence for density-dependence? South African Journal of Marine Science 8:91101.
Barot S., Heino M., Morgan M.J., Dieckmann U. (2005) Maturation of Newfoundland America plaice (Hippoglossoides platessoides): long-term trends in maturation reaction norms despite low fishing mortality? ICES Journal of Marine Science 62:5664.
Barot S., Heino M., O'Brien L., Dieckmann U. (2004) Estimating reaction norms for age and size at maturation when age at first reproduction is unknown. Evolutionary Ecology Research 6:659678.[Web of Science]
Bauchot M.L. and Pras A. (1980) Guide des Poisson Marins d'Europe(Delachaux and Niestlé, Paris) 427 pp.
Begg G.A. (2005) Life history parameters. In Cadrin S.X., Friedland K.D., Waldman J.R. (Eds.). Stock Identification Methods, Applications in Fishery Science(Academic Press, San Diego) pp. 119150 719 pp.
Begg G.A., Hare J.A., Sheehan D.D. (1999) The role of life-history parameters as indicators of stock structure. Fisheries Research 43:141163.[CrossRef][Web of Science]
Begg G.A. and Waldman J.R. (1999) An holistic approach to fish stock identification. Fisheries Research 43:3544.[CrossRef][Web of Science]
Bellier E., Planque B., Petitgas P. (2004) Historical fluctuation of spawning area of anchovy (Engraulis encrasicolus) and sardine (Sardina pilchardus) in the Bay of Biscay from 1967 to 2004. ICES Document CM 2004/Q: 01.
Berg E. and Albert O.T. (2003) Cod in fjords and coastal waters of North Norway: distribution and variation in length and maturity at age. ICES Journal of Marine Science 60:787797.
Beverton R.J.H. (1992) Patterns of reproductive strategy parameters in some marine teleost fishes. Journal of Fish Biology 41:Suppl. B, 137160.[CrossRef][Web of Science]
Borges M.F., Santos A.M.P., Crato N., Mendes H., Mota B. (2003) Sardine regime shifts off Portugal: a time series analysis of catches and wind conditions. Scientia Marina 67:235244.[Web of Science]
Carrera P. and Porteiro C. (2003) Stock dynamic of the Iberian sardine (Sardina pilchardus, W.) and its implication on the fishery off Galicia (NW Spain). Scientia Marina 67:Suppl. 1, 245258.[Web of Science]
Connover D.O. and Munch S.B. (2002) Sustained fisheries yields over evolutionary time scales. Science 297:9496.
Engelhard G.H. and Heino M. (2004) Maturity changes in Norwegian spring-spawning herring before, during, and after a major population collapse. Fisheries Research 66:299310.[CrossRef][Web of Science]
FAO. (2001) Report of the FAO Working Group on the assessment of small pelagic fish off northwest Africa. FAO Fisheries Report 657:133.
FAO. (2003) Report of the FAO Working Group on the assessment of small-pelagic fish off northwest Africa. FAO Fisheries Report 723:54.
Furnestin J. (1943) Contribution à létude biologique de la sardine atlantique (Sardina pilchardus Walbaum). Revue des Travaux de lOffice Scientifique et Technique des Pêches Maritimes 13:211340.
Guisande C., Vergara A.R., Riveiro I., Cabanas J.M. (2004) Climate change and abundance of the Atlanto-Iberian sardine (Sardina pilchardus). Fisheries Oceanography 13:91101.[CrossRef][Web of Science]
Heino M., Dieckmann U., Godo O.R. (2002) Measuring probabilistic reaction norms for age and size at maturation. Evolution 56:669678.[CrossRef][Web of Science][Medline]
Hickling C.F. (1945) The seasonal cycle in the Cornish pilchard, Sardina pilchardus Walbaum. Journal of the Marine Biological Association of the UK 26:115138.
Hilborn R. and Walters C.J. (1992) Quantitative Fisheries Stock Assessment: Choice, Dynamics and Uncertainty(Chapman & Hall, New York) 570 pp.
Hjort J. (1910) Report on Herring investigations until January 1910. Publications de Circonstance du Conseil Internationale pour l' Exploration de la Mer, 53.
Fishery Investigations Hodgson W. C. and Richardson I. D. (1949) The experiments on the Cornish pilchard fishery in 19471948. 17: 22 pp.
Hunter J.R. and Leong R. (1981) The spawning energetics of female northern anchovy, Engraulis mordax. Fishery Bulletin US 79:215230.
ICES. (2000) Report of the Working Group on the assessment of mackerel, horse mackerel, sardine, and anchovy. ICES Document CM 2000/ACFM: 5. 546 pp.
ICES. (2002) Report of the Working Group on the assessment of mackerel, horse mackerel, sardine, and anchovy. ICES Document CM 2002/ACFM: 6. 470 pp.
ICES. (2003) Report of the study group on the estimation of spawning biomass of sardine and anchovy. ICES Document CM 2003/G: 17.
ICES. (2005) Report to the Working Group on the assessment of mackerel, horse mackerel, sardine, and anchovy. ICES Document CM 2005/ACFM: 8. 472 pp.
Kasapidis P., Planes S., Laurent V., Quinta R., Teia A., Lenfant P., Darivianakis S., Terzoglou T., Kotoulas G., Magoulas A. (2004) Stock discrimination and temporal and spatial genetic variation of sardine (Sardina pilchardus) in northeastern Atlantic, with a combined analysis of nuclear (microsatellites and allozymes) and mitochondrial DNA markers. ICES Document CM 2004/Q: 21. 22 pp.
Koops M.A., Hutchings J.A., McIntyre T.M. (2004) Testing hypotheses about fecundity, body size and maternal condition in fishes. Fish and Fisheries 5:120130.[CrossRef][Web of Science]
Krzeptowski M. (1983) Biological characteristics of the sardine (Sardina pilchardus Walb.) off West Sahara. Acta Ichthyologica et Piscatoria 13:1337.
Larrañeta M.G. (1959) Synopsis of biological data on Sardina pilchardus of the Mediterranean and adjacent seas. FAO Species Synopsis 4: 26 pp.
Larrañeta M. G. (1976) Size and age of first maturation, and relative fecundity in Sardina pilchardus off Castellon (Spanish Mediterranean coast). ICES Document CM 1976/J:4. 7 pp.
Le Gall J. (1930) Contribution à létude de la sardine des cotes Françaises de la Manche et de l'Atlantique. Review des Travaux des Pêches Maritimes 3:1947.
Mata M.M., Nespereira J.M.L., Pajuelo J.M.G., Aguilera R.S. (1997) Periodo reproductor y madurez sexual de la sardina Sardina pilchardus (Walbaum, 1792) en aguas de Gran Canaria (islas Canarias). Boletin de l'Instituto Español de Oceanografia 13:4755.
McCullagh P. and Nelder J.A. (1989) Generalized Linear Models(Chapman & Hall, London) 511 pp.
Millar R.B., Broadhurst M.K., MacBeth W.G. (2004) Modelling between-haul variability in the size selectivity of trawls. Fisheries Research 67:171181.[CrossRef][Web of Science]
Mood A.M., Graybill F.A., Boes D.C. (1974) Introduction to the Theory of Statistics 3rd edn. (McGraw-Hill, New York).
Morales-Nin B. and Pertierra J.P. (1990) Growth rates of the anchovy Engraulis encrasicolus and the sardine Sardina pilchardus in the northwestern Mediterranean Sea. Marine Biology 107:349356.[CrossRef]
Morgan M.J. (2004) The relationship between fish condition and the probability of being mature in American plaice (Hippoglossoides platessoides). ICES Journal of Marine Science 61:6470.[CrossRef][Web of Science]
Nash R.D.M., Withames P.R., Pawson M., Alesworth E. (2000) Regional variability in the dynamics of reproduction and growth of Irish Sea plaice, Pleuronectes platessa L. Journal of Sea Research 44:5564.[CrossRef]
Ottersen G., Planque B., Belgrano A., Post E., Reid P.C., Stenseth N.C. (2001) Ecological effects of the North Atlantic Oscillation. Oecologia 128:114.[CrossRef][Web of Science]
Parrish R.H. and Mallicoate D.L. (1995) Variation in the condition factors of California pelagic fishes and associated environmental factors. Fisheries Oceanography 4:171190.[Web of Science]
Parrish R.H., Serra R., Grant W.S. (1989) The monotypic sardines, Sardina and Sardinops: their taxonomy, distribution, stock structure, and zoogeography. Canadian Journal of Fisheries and Aquatic Sciences 46:20192036.
Pérez N., Porteiro C., Alvarez F. (1985) Contribución al conocimiento de la biología de la sardina de Galicia. Boletin de l'Instituto Español de Oceanografia 2:2737.
Pinto J.S. (1957) Essai d'une échelle de maturation sexuelle, pour les males de la sardine (Sardina pilchardus Walb.) basée sur l'histologie des testicules. Notas e Estudos do Instituto de Biologia Marítima 15: 13 pp.
Pinto J.S. and Andreu B. (1957) Echelle pour la caractérisation des phases evolutives de l'ovaire de sardine (Sardina pilchardus, Walb.) en rapport avec l'histophisiologie de la gonade. Proceedings and Technical Papers of the General Fisheries Council of the Mediterranean 46:393411.
Puillat I., Lazure P., Jégou A.M., Lampert L., Miller P.I. (2004) Hydrographical variability on the French continental shelf in the Bay of Biscay, during the 1990s. Continental Shelf Research 24:11431163.[CrossRef][Web of Science]
R Development Core Team. (2004) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-00-3, URL http://www.R-project.org.
Rodríguez-Roda J. (1970) La sardina, Sardina pilchardus (Walb.) del Golfo de Cádiz. Investigaciones Pesqueras 34:451476.
Rodríguez-Roda J. (1971) Estudio de la sardina, Sardina pilchardus (Walb.) capturada por la flota pesquera española en aguas de Larache (Marruecos). Investigaciones Pesqueras 35:655680.
Rose K.A., Cowan J.H., Winemiller K.O., Myers R.A., Hilborn R. (2001) Compensatory density dependence in fish populations: importance, controversy, understanding and prognosis. Fish and Fisheries 2:293327.[CrossRef]
Santos A.M.P., Borges M.F., Groom S. (2001) Sardine and horse mackerel recruitment and upwelling off Portugal. ICES Journal of Marine Science 58:589596.
Santos A.M.P., Kazmin A.S., Peliz A. (2005) Decadal changes in the Canary upwelling system as revealed by satellite observations: their impact on productivity. Journal of Marine Research 63:359379.[CrossRef][Web of Science]
Sedletskaya V.A. (1973) The reproduction of the Moroccan pilchard (Sardina pilchardus, Walb.). Journal of Ichthyology 13:848853.
Silva A. (2003) Morphometric variation among sardine (Sardina pilchardus) populations from the northeastern Atlantic and the western Mediterranean. ICES Journal of Marine Science 60:13521360.
Stratoudakis Y., Bernal M., Borchers D.L., Borges F. (2003) Changes in the distribution of sardine eggs and larvae off Portugal, 19852000. Fisheries Oceanography 12:4960.[CrossRef][Web of Science]
Stratoudakis Y., Coombs S., Halliday N., Conway D., Smyth T., Costas G., Franco C., Lago de Lanzós A., Bernal M., Silva A., Santos M. B., Alvarez P., Santos M. (2004) Sardine (Sardina pilchardus) spawning season in the North East Atlantic and relationships with sea surface temperature. ICES Document CM 2004/Q: 19. 19 pp.
Swain D. P., Hutchings J. A., Foote C. J. (2005) Environmental and genetic influences on stock identification characters. In Cadrin S. X., Friedland K. D., Waldman J. R. (Eds.). Stock Identification Methods, Applications in Fishery Science(Academic Press, London) pp. 4585 719 pp.
Wood S.N. (2000) Modelling and smoothing parameter estimation with multiple quadratic penalties. Journal of the Royal Statistical Society 62:413428.[CrossRef]
This article has been cited by other articles:
![]() |
A. Marcalo, T. A. Marques, J. Araujo, P. Pousao-Ferreira, K. Erzini, and Y. Stratoudakis Fishing simulation experiments for predicting the effects of purse-seine capture on sardine (Sardina pilchardus) ICES J. Mar. Sci., November 4, 2009; (2009) fsp244v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Silva, D. W. Skagen, A. Uriarte, J. Masse, M. B. Santos, V. Marques, P. Carrera, P. Beillois, G. Pestana, C. Porteiro, et al. Geographic variability of sardine dynamics in the Iberian Biscay region ICES J. Mar. Sci., April 1, 2009; 66(3): 495 - 508. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


), 15.5 (), and 18 cm (+) sardine in the 1980s and 1990s off (a) southern Galicia, (b) northern, (c) southwestern, and (d) southern Portugal.



