Population dynamics of the tropical bivalve Cardita affinis from Málaga Bay, Colombian Pacific related to La Niña 1999–2000

Abstract

The population dynamics of the “cholga” Cardita affinis (Sowerby 1833) from Málaga Bay, Colombia, was studied from December 1999 to February 2002, which included the 1999–2000 La Niña event (LN) and the post-LN period 2001–2002. This climatic deviation caused precipitation anomalies in Málaga Bay. Salinity, precipitation, and sea surface temperature anomalies were highly correlated with the bivalve’s body mass cycle. Irregular spawning events were observed during LN by comparison with the regular period. Individual growth and mortality were found significantly higher during LN than during the post-LN period while longevity was almost twofold lower during LN. Increased mortality was probably related to environmental stress. Individual production and productivity were higher during LN, although the annual biomass was lower than during the post-LN period. These results may be related to higher food availability during LN, which agrees well with the results on growth performance. The observed changes provide a base line for future studies regarding effects of El Niño/LN events on population dynamics of tropical bivalves.

Introduction

Coastal marine ecosystems are largely governed by interactions of large-scale oceanographic conditions with climate and hydrology. This is particularly important in the high rainfall area of Pacific Colombia, where high river discharges characterize the environment (Restrepo and Kjerfve 2000; Kjerfve et al. 2001). The El Niño-Southern Oscillation (ENSO) however, causes a significant variability in regional rainfall, river discharge, and sediment load in this area (Restrepo and Kjerfve 2000). The warm episode El Niño (EN) causes drought conditions, while during the cold La Niña phase (LN) intense rainfall occurs in northern South America and in parts of Central America (Ropelewski and Halpert 1987). This is critical in estuarine environments, as ENSO significantly modifies local hydrology, in turn affecting ecosystem dynamics, and ultimately productivity (Peterson 2003). The hydrological anomalies associated with ENSO are therefore probably the most important forcing function for tropical marine ecosystems (Kjerfve et al. 2001; Urban 2001).

In South America, biological effects of ENSO have mainly been studied in economically important species from the highly productive Humboldt Current upwelling ecosystem (e.g. Wolff 1987; Arntz and Fahrbach 1996; Urban and Tarazona 1996). To various extents, marine resources are impacted negatively or positively by EN or LN, depending on the intensity of the event, the onset and the adaptation ability of the species (e.g. Arntz 1986; Wolff 1987; Yamashiro et al 2002; Arntz et al. 2006). However, information on how the population dynamics of estuarine species are affected is still scarce (Riascos 2006).

The bivalve Cardita affinis (Sowerby 1833) is an abundant species in tropical estuarine rocky shores (Libbey and Johnson 1997; Garcia 2002) with a wide distribution range (29°N to 4°S; Enquist et al. 1995). It is considered to be a suitable species for environmental studies, especially to assess effects of EN/LN for the following reasons: (1) this sessile bivalve lives attached in burrows drilled in a narrow zone of the lower estuarine intertidal rocky shore, therefore its population dynamics is likely to be affected by local hydrological and oceanographic changes (e.g. fresh water and nutrient input, sediment load, salinity, sea level, surface temperature), (2) C. affinis is a deposit feeding bivalve, its food availability is affected by hydrological and oceanographic changes, and (3) it is a broadcast spawner, therefore its settlement and recruitment largely depends on local circulation patterns. The aim of this study was to analyze and describe the population dynamics of the tropical bivalve C. affinis under changing environmental conditions, particularly the strong LN 1999–2000 event and the following post-LN period.

Materials and methods

Study area

The study was carried out in La Despensa (3°55′57′′N, 77°20′55′′W) in Málaga Bay (Fig. 1). This small bay of 126 km² with a mean depth of 14 m is located in the Pacific central zone of Colombia. The zone is strongly influenced by the Inter-Tropical Convergence Zone (ITCZ): as the low pressure belt of the ITCZ moves north and south between 10°N and 3°S, it passes over Málaga Bay and causes two rain periods, one from September to November and another one from April to June, interrupted by periods with lower precipitation. The study area is one of the regions with the highest annual precipitation in the western hemisphere (7,000–11,000 mm; Restrepo et al. 2002). The bay is characterized by a semidiurnal tide regime with an average range of 3.5 m. Sea surface temperature (SST) remains nearly constant year-round, ranging between 26 and 29°C while salinity ranges between 18 and 27 (Rubio 1984).

Fig. 1

a Investigation area and monthly mean precipitation (mm month⁻¹; 0.5° resolution) averaged between 1961 and 1990 (Precipitation data from DWD/GPCC) b Location of the studied Cardita affinis belt at La Despensa

Sampling

Quantitative monthly samples of C. affinis were taken from December 1999 to February 2002. Four stations (three replicates, using a 0.5 × 0.5 m size quadrate) were sampled at 100 m intervals of an along-shore transect covering the entire zone, where C. affinis colonizes the rocky cliff (lower intertidal). Sampling was carried out during spring low tide. Maximum individual anterior–posterior shell length (SL) was measured to the lower 0.1 mm with a vernier calliper, thus obtaining monthly length–frequency distributions (LFD). From the quantitative samples, two sub-samples were selected (A: 30–35 individuals covering the available size range of C. affinis, 5–40 mm; B: 30–35 individuals SL > 25 mm). Sub-sample A was used to study body condition index (body mass cycle), while sub-sample B was used to study the reproductive cycle.

Salinity and SST were recorded in situ three times each sampling date by using a YSI sensor. Additionally, monthly SST data from Niño region 1 + 2 (0–10°S; 90–80°W) provided by the Climate Prediction Center-NOAA was included in the analysis to account for associations between large scale and local environmental variability. Monthly mean precipitation (MMP) was recorded by a meteorological station of the “Corporación Autónoma del Valle del Cauca” near the study area (4°N, 77°02′W).

Body mass and reproductive cycle

To describe the body mass cycle, all soft parts were removed and dried to constant mass at 70°C in order to determine shell-free dry mass (SFDM). Ash-free dry mass (AFDM) was obtained by ignition at 550°C for 5 h. A linear regression analysis was carried out to estimate the constants “a” and “b” of the allometric equation

$$ \log ({\text{AFDM}}) = \log (a) + b \times \log ({\text{SL}}). $$

(1)

This equation was used to estimate monthly variability of the shell size–mass relationship for a standard individual of 30 mm SL and to identify spawning events. Decreasing mass between two successive months was used as indication of a spawning event, while an increase over a longer period was interpreted as the developing phase before the spawning season (Urban and Tarazona 1996). Additionally, the monthly gonado-somatic index (GSI, Eq. 2), defined as the relative volume of gonadal tissue (V _gon) from the total body volume (V _body), was estimated using a method based on linear measurements of the gonad region, which forms a sheath around the digestive gland (for details see Urban and Riascos 2002):

$$ {\text{GSI}} = V_{{{\text{gon}}}} /V_{{{\text{body}}}} \times 100. $$

(2)

Body mass may also be affected by short-term changes in food availability. To avoid this potential bias, we construed the body mass cycle in connection with gonad studies. For this we determined the gonad stage from microscopic observations on fresh gonadal material (smear samples). We used a semi quantitative scale proposed by Guillou et al. (1990), which allows to classify both males and females into four gonad stages: indifferent, ripe I, ripe II and spent. A Spearman correlation analysis was performed to evaluate the relationship between environmental parameters and variables describing the body mass and reproductive cycle.

Growth

Each month, 25–30 individuals were marked in situ with a carborundum disc of 1 mm thickness mounted in the mandrel of a portable electric grinder. In order not to disturb the bivalve too much, only a shallow groove was applied at the exposed anterior end of the shell, still leaving the bivalve attached in its burrow. The length from the groove to the anterior margin was measured with a vernier calliper to the nearest 0.1 mm. A number tag for individual identification was fixed close to the borrow of each marked specimen. After different time intervals ranging from 1 to 6 months individuals were extracted from their cavities and SL at recapture time (SL_t2) recorded. SL at marking time (SL_t1) was estimated by subtracting the length increment (distance between the groove and the anterior shell margin) from SL_t2. Two sets of tagging-recapture data were analysed independently. The first set consisted of 84 individuals recaptured during the period November 1999 to December 2000 (LN), and the second set of 97 individuals recovered from January 2001 until February 2002 (post-LN).

Growth was described by the von Bertalanffy growth function (VBGF; von Bertalanffy 1938):

$$ L_{t} = L_{\infty } {\left( {1 - e^{{ - K{\left( {t - t_{0} } \right)}}} } \right)} $$

(3)

where L _t is the length at age t (year), t ₀ the age at zero length, L _∞ is the asymptotic length (mm) and K, the growth constant (year⁻¹). Growth parameters K and L _∞ were estimated using Fabens’ method (1965) by fitting a rearranged function of Eq. (3) to tagging-recapture data using the iterative non-linear least-square method (Simplex algorithm, Press et al. 1986):

$$ L_{2} = L_{1} + {\left( {L_{\infty } - L_{1} } \right)} {\left( {1 - e^{{ - K{\left( {t_{2} - t_{1} } \right)}}} } \right)} $$

(4)

where L ₁ is the length at marking time (t ₁) and L ₂ the length at recapture time (t ₂). As our data lack size-increment measurements referring to larger individuals L _∞ was not determined iteratively, but set to the maximum length observed.

To compare the VBGF obtained for each period (LN and post-LN) an analysis of the residual sum of squares (RSS) was performed (Chen et al. 1992). The index of growth performance phi prime (Φ′; for a review see Brey 2001) was used as a measure of growth performance:

$$ \Phi^{\prime} = \log {\text{ }}K + 2 \times {\text{log }}L_{\infty } $$

(5)

Mortality and longevity

Total mortality Z (year⁻¹) was calculated independently for the two periods (LN and post-LN) with the single negative exponential model on pooled length–frequency data:

$$ N_{t} = N_{0} e^{{ - Zt}} $$

(6)

(where t is the time and N ₀ is the number of individuals at t = 0) and the length converted catch curve on the pooled length–frequency data (Pauly 1983). Thereby, using the parameters of the VBGF, length of pooled length–frequency samples was converted into age as:

$$ N_{i} /\Updelta t_{i} = {N^{\prime}_{0}} e^{{ - Zt_{i} }} $$

(7)

where N _i is the number of individuals in length class i, Δt _i is the time required to grow through this size class, and t _i is the age of the middle size class i. Mortality Z is calculated by linear regression analysis:

$$ {\text{Ln}}\;{\left( {N_{i} /\Updelta t_{i} } \right)} = a + bt_{i} ;\quad Z = - b $$

(8)

Longevity (Tm) was calculated as:

$$ {\text{Ln }}(Z) = a + b \times {\text{ln}}({Tm}) $$

(9)

where a = 1.23 and b = −0.832 are constants, estimated from different mollusc stocks (see Hoenig 1983).

Production

Total annual somatic production P _S (in g AFDM m⁻² year⁻¹) was calculated for the intertidal C. affinis belt by the mass specific growth rate method (Crisp 1984) from the mean of quantitative samples, pooled length–frequency data, the VBGF parameters and the size–mass relationship (Eq. 1) for each period (LN and post-LN):

$$ P_{{\text{s}}} = {\sum {N_{i} M_{i} G_{i} } } $$

(10)

where N _i is the mean abundance of individuals (N m⁻²) estimated from quantitative samples, M _i the mean body mass (g AFDM) in length class i and G _i , the mass-specific growth rate (year⁻¹):

$$ G_{i} = b{\text{ }}K{\text{ }}[(L_{\infty } /L_{i} ) - 1]{\text{ }} $$

(11)

where b is the constant of the allometric Eq. 1, K and L _∞ are VBGF parameters and L _i is the mean length in size class i. Individual somatic production P _i (in g AFDM year⁻¹) was calculated as:

$$ P_{i} = {\sum {M_{i} G_{i} } } $$

(12)

and mean annual biomass $\overline{B}$ (in g AFDM m⁻² year⁻¹) of the population as:

$$ \overline{B} = {\sum {N_{i} M_{i} } } $$

(13)

The ${\text{P}}_{{\text{s}}} /\overline{B}$ ratio (productivity) of the C. affinis population was calculated from annual somatic production P _s and annual mean biomass $\overline{B}.$

Results

Reproduction

Body mass cycle in relation to GSI and the distribution of reproductive stages is shown in Fig. 2. Animals in the spent stage were found throughout the entire study period. However, a main spawning period, indicated by a clear body mass decrease, a declined GSI and increasing microscopic observations of spent stages between January and March 2001 is evident (Fig. 2). Two minor spawning events occurred from April to June 2000 and 2001. An additional spawning event seems to occur between October and November 2000 as indicated by the strong increase of spent stages, and a decrease of ripe animals. However, no reduction of body mass was observed. In general, higher proportions of ripe animals were observed during LN compared with post-LN period, whereas GSI values were consistently higher during the post-LN period.

Fig. 2

Body mass (g ash-free dry mass) cycle for a 30 mm standard individual (calculated from size-mass relationship, Eq. 1) and reproductive cycle of Cardita affinis based on the gonado-somatic index (GSI) and gonadal stages from microscopic observations of smear samples. Spawning events are indicated by black bars on top: Bar size corresponds to the extent of the spawning events indicated by body mass decreases; GSI decreases and spent stage increases between successive months. La Niña 1998–2000 chronology based on data from NOAA/National Weather Center-Climate Prediction Center

The Spearman correlation analysis showed significant correlations between reproductive and environmental variables: a positive correlation exists between the body mass cycle, salinity and SST anomalies of the Niño 1 + 2 region; whereas a negative relationship was observed between body mass cycle and monthly mean precipitation (Table 1; Fig. 3). The latter was higher during LN; in fact, during 2000 it reached 13,939 mm in comparison to an annual historical (1984–2004) average of 11,435 mm (Corporación Autónoma del Valle del Cauca 2000). The body mass cycle did not show a clear annual pattern.

Table 1 Spearman correlation matrix of reproductive factors and abiotic variables

Fig. 3

Relationship between body mass cycle of Cardita affinis and environmental variability during the study period. Salinity, sea surface temperature anomalies of Niño region 1 + 2 (SST anomalies), monthly mean precipitation in the study area (MMP) from a government agency in Málaga Bay were plotted. La Niña 1998–2000 chronology based on data from NOAA/National Weather Center-Climate Prediction Center

Growth, mortality and longevity

In total 815 individuals were marked, of which 181 (22.2%) were recovered. Estimated VBGF parameters (L _∞ and K) and the corresponding Φ′ index are given in Table 2. RSS analysis showed that growth during LN was significantly higher (F _{3, 175} = 9.24; P < 0.05) than during the post-LN period. The results of the size-converted catch curves (Fig. 4) indicate a higher mortality rate of C. affinis during LN (Z = 0.53) compared to the post-LN period (Z = 0.32). Longevity was lower during LN (9.21 year) than during post-LN (17.25 year).

Table 2 Cardita affinis. Von Bertalanffy growth function parameters (± standard error) estimated with Fabens’ method (1965) from two marking-recapture experiments (LN and post-LN)

Fig. 4

Length-converted catch curves of Cardita affinis for two periods (LN 1999–2000 and post-LN 2001–2002), based on pooled length–frequency data. Von Bertalanffy growth parameters required for this method are given in Table 2. Filled data points: used for regression; open data points: excluded from regression. Regression equation, r ² and Z values are given

Production

Individual somatic production during LN increased to a maximum of 0.031 g AFDM year⁻¹ at 26.5 mm SL and decreased thereafter, whereas it was lower during the post-LN period reaching 0.028 g AFDM year⁻¹ at 28.5 mm SL (Fig. 5).

Fig. 5

Individual somatic production and annual population production of Cardita affinis during 1999–2000 (LN) and 2001–2002 (post-LN); a individual somatic production for different size classes; b annual population production and mean abundance for different size classes

The LFD of the C. affinis population showed a unimodal size distribution, with medium-sized individuals prevailing. However, during LN larger individuals dominated the population, although the mean abundance was almost the same for the two periods (Fig. 5; Table 3). Taking into account that medium to large sized individuals supported the highest individual somatic production, the shift to smaller individuals during 2001–2002 (post-LN) leads to a reduction of total annual production. The estimates for annual production, annual biomass and $P_{s} /\overline{B}$ ratio are summarised in Table 3.

Table 3 Summary of production estimation for Cardita affinis during La Niña and the following post-La Niña period

Discussion

Cardita affinis showed continuous gametogenic activity (Fig. 2) as it is commonly reported for bivalves from tropical regions, where interannual temperature variability is low (Urban 2001). However, periods of increased activity are still identifiable and often correlated to climatic variation (Urban 2000, 2001; Riascos 2006).

During the study period the reproductive cycle of C. affinis was characterized by: (1) the absence of an evident seasonal spawning pattern (Fig. 2), (2) a lack of correlation between GSI/ body mass versus the percentage of spent animals, and (3) the fact that body mass was the only parameter being significantly correlated with environmental factors (Table 1).

The absence of a seasonal spawning pattern may indicate an abnormal environmental condition the species has to cope with and a low ability to recover the normal rhythm. Some authors described changes in spawning patterns in marine bivalves related to environmental anomalies during ENSO (Birket and Cook 1987; Arntz et al. 1987; Urban and Tarazona 1996). In the study area interannual variability of water discharge and sediment load are associated with ENSO (Restrepo and Kjerfve 2000). During LN Southeast trade winds are well developed and the ITCZ remains north of its typical position (Kjerfve et al. 2001). As a result, rainfall intensifies, thus increasing MMP. This effect is clearly visible in Fig. 3, the MMP showing significant differences (P < 0.05) between the LN and the post-LN period.

The lacking correlation between both GSI/body mass and the percentage of spent animals may be an indication of the species mechanism to face atypical conditions: gamete re-absorption has been observed as a mechanism to reallocate energy to metabolic maintaining under environmental stress (Urban and Tarazona 1996). This mechanism may explain the irregular spawning event observed in September–October 2000, in which the percentage of spent stages increased, but did not correspond with a decrease of GSI and body mass (i.e. gametes were not released). Alternatively, spawning could have taken place, but may have been masked by the increasing trend of body condition and GSI. This explanation is supported by the quick recovery of ripe individuals in December 2000 and January 2001 (Fig. 2). In the last case, we would expect a recruitment event during LN (as recorded during the post LN period) but it did not occur. However, competent larvae may have been recruiting elsewhere or may have died under those anomalous environmental conditions.

The high correlation between body mass changes and the environmental parameters (Table 1; Fig. 3) indicates that environmental variability may be related to food availability. Terrestrial runoff, river discharge and resuspension of bottom sediments by tidal currents have been identified as the principal source of nutrients in tropical estuarine ecosystems (von Prahl et al. 1990; Cantera and Blanco 2001; Medeiros et al. 2001), which have a fast turnover rate (Alongi 1994). Thus, it can be proposed that during the rainy season with high river discharge, primary productivity (and thus food availability) raises, as has been observed in tropical estuaries in Costa Rica (Gocke et al. 2001) and coastal lagoons of the Caribbean coast of Colombia (Hernández and Gocke 1990). Therefore, the body mass increase observed during November 2000 and January 2001 and the main spawning event in January/March 2001 may be related to the abnormally high MMP between August and November 2000 (LN) and the subsequent increase in primary production.

Growth was significantly higher during LN than during the post-LN period, though the asymptotic length was slightly higher during the latter (Table 2). Accordingly, growth performance was higher during LN. Increased growth performances of suspension feeding bivalves could be mainly related to food availability (Peterson 1982; Sastre 1984; Nakaoka 1992; Laudien et al. 2003; Heilmayer et al. 2004). As ENSO-induced increases in freshwater flow seem to boost estuarine phytoplankton biomass and production in tropical estuarine ecosystems (e.g. Abreu et al. 1994; Ciotti et al. 1995), the observed differences in growth performance could be related to increased food availability during LN.

Results on growth and reproduction support the hypothesis of favourable feeding conditions during LN. However, the results on mortality seem to contradict this theory. Mortality was higher during LN, precisely when lowest mortality is expected due to presumably favourable conditions. In a study on the influence of run-off on growth and survival of the Sydney rock oyster, Saccostrea glomerata, Paterson et al. (2003) observed a similar pattern: elevated seston loads associated with low salinities were reflected in an initial increased oyster growth. However, higher mortality rates were reported when run-off was too high. While storms, which temporarily increase turbidity and sedimentation in estuaries, do not necessarily affect bivalves, continued sediment inputs (e.g. river discharge) significantly affect their physiology (Norkko et al. 2006). Anomalous high run-off lasted until the end of 2000, which may initially facilitate growth of C. affinis, but was detrimental when associated stress surpassed the tolerance range of this species, thus leading to both high growth and a higher mortality rate in 2000. Extended periods of wave and salinity-stress may also promote detrimental effects on sedentary fauna (Ritter et al. 2005; Zardi et al. 2006). Nevertheless, physiological experiments are needed to determine the tolerance to different levels of these stressors and disentangle their relative importance on growth and mortality. Results on mortality are reflected in the observed longevity, which was almost twofold during 2001–2002 in comparison with the LN period. While a life span of 17 years seems too long, it is not unusual for tropical bivalves (see Stern-Pilot and Wolff 2006).

The higher growth rate observed during LN (K = 0.23) compared to the post-LN period (K = 0.13) led to a higher individual production, except for the largest individuals (SL ≥ 34 mm) (Fig. 5a). The LFD shows that only medium-sized individuals dominated the population and juveniles were scarce (Fig. 5b). This might be the result of strong recruitment variability associated to EN 1997–1998 and LN 1999–2000. The strong spawning event observed in January–March 2001 and the return to normal environmental conditions can explain the observed increase of small individuals (7–8 mm) during 2001–2002 (Fig. 5b).

During LN the C. affinis population consisted of slightly larger individuals, even when the mean abundance was almost the same for the two periods (Table 3). Taking into account that medium to large sized individuals supported the highest individual somatic production, the shift to larger individuals during 2001–2002 (post-LN) lead to a reduction of total annual production. Finally, although the annual biomass was lower during LN, the population productivity $(P_{s} /\overline{B} {\text{ ratio}})$ was higher, thus giving further support to the assumption of favourable food conditions. The observed differences in productivity of C. affinis are consistent with changes in productivity of the surf clam Donax dentifer inhabiting sandy beaches at Málaga Bay (Riascos 2006): productivity during LN 1999–2000 (2.34) was higher by almost a factor of two in relation to that of the post-LN period (1.21). This suggests that LN brings about strong changes in coastal energy flow.

Population dynamics during ENSO have been mostly studied at higher latitudes, where changes are mainly related to temperature anomalies. Our study shows that in estuarine ecosystems LN modifies freshwater input, which in turn may affect biogeochemistry and primary production in coastal waters (Ciotti et al. 1995), thus promoting changes on growth, mortality, reproduction and somatic production of C. affinis. Unfortunately, similar studies on comparable ecosystems and related species are still lacking, thus comparisons are precluded. Owing to this and the fact that interannual variability in population dynamics is expected under normal periods, our study provides a base line for future studies on ENSO-related effects in tropical estuarine ecosystems.