Macrozoobenthic community structure in a poikilohaline Mediterranean lagoon (Laki Lagoon, northern Aegean)

Abstract

The community structure of the macrobenthic fauna was studied in the overall area of Laki Lagoon in September 1997 (salinity 32–35 psu) and monthly from February 1998 to February 1999 in the innermost part of the lagoon (salinity 0.1–6.8 psu). Community structure was analyzed by means of uni- and multivariate methods. In September 1997, the macrofauna in the outer part of the lagoon was characterized by a higher diversity and the occurrence of both lagoonal and marine species, and in the innermost part by a higher total abundance and the occurrence of a few lagoonal species. The combination of distance from the sea, depth, salinity and sediment organic matter correlated best with the spatial distribution pattern of the macrobenthic fauna. Community structure in the innermost part of the lagoon showed a seasonal periodicity. Species composition during spring 1998, at 0.1–2.0 psu, was similar to that in September 1997. During summer the macrobenthic fauna became impoverished, but recovered from late summer onwards. The salinity increase during summer (up to 5–7 psu) was followed by the appearance of marine species in the innermost part of the lagoon. Total abundance displayed a peak in late spring and a lower one in mid-autumn. The seasonal dynamics of the faunal assemblage was mainly governed by water temperature. Predation pressure by Atherina boyeri may have contributed to quantitative community changes during autumn.

Introduction

Coastal lagoons, which are of considerable naturalistic and economic interest, are usually characterized by marked daily and seasonal fluctuations in environmental conditions (e.g. in salinity) depending on their geomorphology and hydrological regime. Additionally, in many cases, these coastal environments suffer severe dystrophic crises during summer, being followed by mortality events. Although coastal lagoons have been the subject of extensive studies, our knowledge of the factors determining the spatial and seasonal variation in their macrozoobenthic communities, especially in the microtidal Mediterranean Sea, is restricted and mainly related to abiotic rather than biotic factors (e.g. predation pressure). Useful information on the distribution and seasonal dynamics of the macrobenthic fauna in Mediterranean coastal brackish habitats has been given in some recent publications (e.g. Reizopoulou et al. 1996; Gouvis et al. 1997; Kevrekidis 1997; Lardicci et al. 1997; Tagliapietra et al. 1998; Arvanitidis et al. 1999; Bachelet et al. 2000; Koutsoubas et al. 2000; Mistri et al. 2001; Mistri 2002). However, very little information exists on the seasonal dynamics of macrozoobenthic communities in low-salinity parts of Mediterranean lagoons.

Laki Lagoon, one of the three lagoons of the Evros Delta (northern Aegean), which is protected as a wetland of international value according to the Ramsar Convention, is of notable natural history and recreational interest. Nevertheless, information on the macrobenthic fauna of Laki Lagoon is mainly limited to polychaetes, amphipods and molluscs (Kevrekidis and Koukouras 1988; Gouvis and Koukouras 1993; Kevrekidis et al. 1996; Gouvis et al. 1998).

The main goal of this study is a better understanding of the structure and function of coastal brackish ecosystems through (1) the description of the distribution of the macrobenthic fauna in Laki Lagoon and of the seasonal variation in macrozoobenthic community structure and dynamics at very low salinities, and (2) the investigation of the key variables affecting this spatial and seasonal variation.

Methods

Study area

Evros Delta is located at the NE end of the Aegean Sea (Fig. 1). Fresh water flows into the delta area primarily from the eastern branch of the Evros river; additional fresh water enters the delta area through the western branch of the Evros and the streams Mikri Maritsa and Loutron, usually from autumn to early summer. Three islets and some lagoons have been formed in the delta area. Laki Lagoon, occupying an area of about 100 ha, communicates with the sea via two openings, both situated on the western side of the lagoon (Fig. 1).

Fig. 1

A Geographical location of the study site. Maps of B the Evros Delta and C Laki Lagoon showing the sampling stations

Sampling and laboratory techniques

Samples of the benthic fauna were collected in September 1997 at sampling sites LA₁ and LA₂ located on a transect across the innermost, northern part of Laki; at sampling sites LC₁, LC₂ and LC₃ located on a transect across the southernmost part of the lagoon; and at sampling site LB located at the central part of the lagoon (Fig. 1). At each site, two random replicate units were taken using a van Veen-type special sampler with handles (Larimore 1970); the sampler covered a surface of 400 cm² (20×20 cm²) and penetrated to a depth of 20 cm. The samples were sieved on a 0.5-mm mesh sieve and the animals were kept in 5% formalin solution. Depth and water salinity near the bottom were also measured. Sediment samples were taken using a small corer for particle size analysis and estimation of the amount of organic matter.

Monthly samples were collected from February 1998 to February 1999 at station LA₁ (Fig. 1). Four random replicate units were taken using the modified van Veen grab. The samples were sieved through a 0.5-mm screen. Sediment samples were taken using a small corer. Depth, salinity, temperature and dissolved O₂ of the water near the bottom were also measured. Additionally, two samples of epibenthic decapods and small fish were taken using a special net (with a 40×40 cm opening) which was pulled on the sediment for a distance of 10 m.

Sediment analysis and estimation of organic matter were made according to the methods described by Buchanan (1984). All animals were identified to the lowest possible taxon and counted.

Data analysis

Macrobenthic community structure was analyzed by the total number of individuals (n), number of species (S), and Margalef’s species richness (d), Shannon–Wiener species diversity (H′) (natural logarithm) and evenness (J) indices.

The spatial and seasonal macrofaunal patterns were tested using the Analysis of Similarities (ANOSIM) test. Non-metric multidimensional scaling (MDS) was used to investigate faunal similarities among samples (data transformed to fourth root). Species responsible for similarities and for differences between sample groups were investigated using the similarity percentages routine (SIMPER) (Clarke 1993). Environmental variables which correlated best with the multivariate pattern of the macrobenthic community were identified by means of harmonic Spearman coefficient, ρ_w (BIO–ENV analyses) (Clarke and Ainsworth 1993). Spearman’s rank correlation coefficient was applied to identify the highly correlated environmental variables. The PRIMER package developed at Plymouth Marine Laboratory was used.

Results

The environment

Abiotic variables

The distance of the sampling stations from the sea and the values of abiotic variables of water and sediment at these stations in September 1997 are shown in Table 1. Station LB and transect LC were located closer to the sea (Fig. 1, Table 1). Transect LA was characterized by shallow waters, while along transect LC sediments were coarser and sediment organic matter content was lower (Table 1). Salinity near the bottom varied between 32.2 psu and 35.1 psu throughout the lagoon (Table 1).

Table 1 Sampling stations in Laki: distance from the sea and physicochemical parameters of water and sediment in September 1997

An increased freshwater inflow during winter and spring 1998, following a period of intense rainfall on the Evros catchment area, resulted in a sharp decline in salinity in Laki. Salinity near the bottom was close to 0 psu in February/March 1998; it increased during spring, peaked in late summer (6.8 psu) and then gradually decreased until February 1999 (0.2 psu) (Fig. 2). The values of abiotic factors of water and sediment at station LA₁ from February 1998 to February 1999 are shown in Table 2. Depth varied between 5 cm and 50 cm. Water temperature near the bottom was lowest in February/March 1998 (~10.5°C) and from November 1998 to January 1999 (6.7–9.6°C), and highest in July (25.9°C) (Table 2; Fig. 2). Dissolved oxygen showed relatively high values in early spring, autumn and winter and was minimal in June (5.3 mg/l). The sediment consisted of very fine sand almost throughout the sampling period; sediment organic matter showed comparatively high values from July to October (1.17–1.66%) and in February 1999 (1.48%) (Table 2). These values were much lower than those reported from several other Mediterranean lagoons (e.g. Reizopoulou et al. 1996; Lardicci et al. 1997; Tagliapietra et al. 1998).

Fig. 2

Monthly variation in salinity and water temperature at station LA₁ from February 1998 to February 1999

Table 2 Monthly variation of physicochemical parameters of water and sediment at station LA₁ during the sampling period

Biotic variables

Macroalgae were present in the innermost part of Laki Lagoon (station LA₁) (Fig. 1) throughout the study period (February 1998–February 1999): their biomass (as total algal dry weight) varied seasonally, showing main peaks in August (135.43 g m⁻²) and in November (101.20 g m⁻²) (unpublished data; see also Table 3).

Four species of decapods and five species of small fish were found in the monthly samples taken with the special net at station LA₁ from February 1998 to February 1999. The monthly variation in the number of individuals/10 m² is given in Table 3. Most of them were mainly found in summer and/or autumn; Atherina boyeri was collected in September and in November.

Table 3 Monthly variation in total biomass of macroalgae (g DW m⁻²) (unpublished data) and in number of individuals/10 m² of epibenthic decapods and small fish collected with a special net at station LA₁ during the sampling period

Faunal composition

A total of 25 macrobenthic invertebrates was found along transects LA and LC and at station LB in September 1997. Sixteen of these species displayed a cumulative mean dominance of 99% (Table 4). The polychaete Streblospio shrubsolii and Tubificidae were found at all stations showing a cumulative mean dominance of 48% (Table 4). The amphipods Corophium orientale and Gammarus aequicauda, the polychaete Capitomastus minimus and the gastropod Ventrosia maritima also showed a remarkable mean dominance although they were less frequently collected (Table 4).

Table 4 Taxa of the assemblages at the sampling stations in Laki lagoon in September 1997 showing the highest values of presence (P), mean dominance (D_m), and cumulative dominance (D_c)

Nineteen taxa were collected at station LA₁ throughout the annual cycle (February 1998–February 1999). Thirteen of these taxa showed a cumulative mean monthly dominance of 99.6%; four species were constantly found throughout the study period (Table 5). The amphipod C. orientale dominated the assemblage, showing a mean monthly dominance of 54.6%. The polychaetes S. shrubsolii and Hediste diversicolor and the bivalve Abra ovata displayed a more or less remarkable mean monthly dominance. Four taxa (Tubificidae, V. maritima, G. aequicauda and Cerastoderma glaucum) were frequently found (in 10 or 9 months); however, only the first two taxa showed a notable mean monthly dominance.

Table 5 Taxa of the assemblage at station LA₁ showing the highest values of presence (P) and mean monthly dominance (D_m) from February 1998 to February 1999, and cumulative mean monthly dominance (D_c)

Structural analyses

Univariate analysis

The major biological parameters (total density n, number of species S, species richness d, species diversity H′ and evenness J) at each station in September 1997 are given in Table 6. The n showed its highest values on transect LA, especially at station LA₁, and was minimal at station LB; S and d had comparatively higher values at stations LC₁ and LC₂ and the lowest values at station LB. H′ had its highest values on transect LC and was minimal at station LB, as well; J showed its highest value at station LC₃.

Table 6 Number of individuals (n, individuals m⁻²), species number (S), species richness (d), species diversity (H′) and evenness (J) of the macrobenthic fauna at the sampling stations in September 1997

Total macrobenthic density at station LA₁ increased sharply from February/March 1998 (mean values 8,544 and 5,406 individuals m⁻², respectively) to May (22,225 individuals m⁻²); then it decreased sharply to September (2,906 individuals m⁻²) (Fig. 3). Density temporarily increased in October (9,463 individuals m⁻²) and then again gradually from November to January/February 1999 (7,231–7,719 individuals m⁻²). Mean value of S ranged between 6.75 and 7.75 from February 1998 to July (with the minimum in July), between 7.5 and 9 from August to November, and between 10.25 and 11.75 from December 1998 to February 1999 (Fig. 4). Mean value of d varied between 0.605 and 0.748 from February 1998 to July (with the minimum in July), between 0.799 and 1.007 from August to November, and between 1.041 and 1.249 from December 1998 to February 1999 (Fig. 4). Mean values of H′ and J fluctuated between 0.787 and 1.7 and between 0.363 and 0.804, respectively, being lowest from May to July and in October (Fig. 4). A temporary sharp increase in abundance of C. orientale contributed to the low values of H′ and J in October.

Fig. 3

Monthly variation in the total number of individuals of the macrobenthic fauna from February 1998 to February 1999 at station LA₁ (means and SE)

Fig. 4

Monthly variation in the number of species, species richness, species diversity and evenness of the macrobenthic fauna from February 1998 to February 1999 at station LA₁ (means and SE)

Multivariate analysis

Two significantly distinct groups were formed by the samples collected at the sampling stations in September 1997 (ANOSIM, global R=0.89, P<0.001). The first group was formed by the samples collected in the outer, southernmost part of the lagoon (stations LC₁–LC₃) along with a sample collected at station LB, and the second group by the samples collected in the shallower, innermost northern part of the lagoon (stations LA₁, LA₂). The MDS ordination plot of the samples collected at the sampling stations in September 1997 is shown in Fig. 5.

Fig. 5

Multidimensional scaling (MDS) ordination plot of the samples collected at Laki in September 1997

Four taxa (S. shrubsolii, Tubificidae, Aricidea cerrutii and Glycera tridactyla) contributed >65% to the similarity in the first sample group, and five species (C. orientale, V. maritima, A. ovata, S. shrubsolii and C. glaucum) in the second group (SIMPER) (Table 7). The differences between these groups were mainly due to (1) the exclusive occurrence of three species in the first group (A. cerrutii, G. tridactyla, C. minimus) and of two species in the second one (V. maritima, H. diversicolor), and (2) the higher average density of four taxa (C. orientale, A. ovata, Tubificidae, C. glaucum) in the second group (SIMPER, cut-off 65%) (Table 8). A sample collected at station LB (sample1 LB) was not grouped (Fig. 5); macrofauna of this sample was extremely poor (only two taxa in extremely low densities) (Table 8).

Table 7 Average densities (individuals m⁻²) of taxa which contributed >65% to the similarity in each sample group at Laki in September 1997 (SIMPER)

Table 8 Differences (< or >) in average densities (individuals m⁻²) of taxa which contributed >65% to the dissimilarity between the sample groups at Laki in September 1997 (SIMPER)

Four significantly distinct groups were formed by the monthly samples (ANOSIM, global R=0.815, P<0.001) corresponding to different periods of the year. The first group included the samples of February 1998 to May, the second those of June and July, the third those of August to October, and the fourth those of November to February 1999. The MDS ordination plot of the monthly samples at station LA₁ is shown in Fig. 6.

Fig. 6

Multidimensional scaling (MDS) ordination plot of the monthly samples at station LA₁; FEB February 1998; FEBR February 1999

The average densities of the species which contributed >65% to the similarity in each monthly sample group (SIMPER) are given in Table 9; it can be shown that almost the same taxa contributed to the similarity in each group. The average densities of the species which were responsible for the differences between the successive monthly sample groups (cut-off 65%; SIMPER) are given in Table 10. The differences between late winter/spring samples and those of summer were mainly caused by the disappearance of some taxa (Tubificidae, Cumacea), the decline in density of some molluscs (V. maritima, A. ovata) and the increase in density of the epibenthic amphipod G. aequicauda and the “enriched opportunist” polychaete S. shrubsolii in summer. The differences between summer and late summer/autumn samples were primarily due to the appearance of some taxa (Tubificidae, C. glaucum, Microprotopus maculatus), the increase in density of A. ovata and the decline in that of V. maritima during late summer/autumn. The differences between late summer/autumn samples and those of late autumn/winter were mainly caused by the appearance of some polychaetes (Polydora ciliata, G. tridactyla, A. cerrutii), the increase in density of some taxa (Tubificidae, S. shrubsolii) and the decrease in that of G. aequicauda, V. maritima and C. glaucum during late autumn/winter (Table 10).

Table 9 Average densities (individuals m⁻²) of taxa which contributed >65% to the similarity in each monthly sample group from February 1998 to February 1999 at station LA₁ (SIMPER)

Table 10 Differences (< or >) in average densities (individuals m⁻²) of taxa which contributed >65% to the dissimilarity between the monthly sample groups from February 1998 to February 1999 at station LA₁ (SIMPER)

Associated environmental variables

The highest values of the harmonic Spearman rank coefficient (ρ_w), deriving from the performance of the BIO–ENV analysis, are given in Table 11. From the examined abiotic variables (Table 1), the set of distance from the sea, depth, salinity and sediment organic matter correlated best with the distribution pattern of the macrobenthic community; the set of distance from the sea and salinity showed the next highest value of ρ_w (Table 11). From a total of 16 abiotic and biotic variables (Tables 2, 3), water temperature correlated best with the monthly variation in the macrobenthic community structure at station LA₁ (Table 11). The set of water temperature with the abundance of A. boyeri showed the next highest value of ρ_w (Table 11).

Table 11 Summary of the contributions of the environmental variables showing the highest values of the harmonic Spearman rank coefficient (ρ_w) with the distribution pattern of the macrobenthic community in Laki lagoon in September 1997, and with the monthly variation in the macrobenthic community structure at station LA₁ from February 1998 to February 1999

Discussion

The distribution of the macrobenthic fauna in Laki seems to conform to the “confinement” theory proposed by Guelorget and Perthuisot (1992). The outer areas of Laki lagoon, close to the lagoon mouths, showed a different faunal composition compared to the innermost part of the lagoon in early autumn at salinities of 32–35 psu, indicating a biological zonation.

Macrofauna in the outer part of Laki was characterized by a higher diversity and the occurrence of both brackish-water taxa (e.g. S. shrubsolii, tubificids, C. orientale, A. ovata) and marine species (e.g. A. cerrutii, G. tridactyla, C. minimus). Most of these marine species have also previously been found in brackish habitats (e.g. Picard 1965; Fèbvre 1968; Nicolaidou and Pitta 1986; Gravina et al. 1988; Nicolaidou and Papadopoulou 1989; Thiermann et al. 1997; Clarke and Warwick 1998; Gouvis et al. 1998). According to Cognetti and Maltagliati (2000) brackish habitats which communicate directly with the sea are usually colonized by euryhaline marine species, whose populations have developed a tolerance to unpredictable changes in habitat, and by stenohaline marine species which have developed euryhaline populations.

The macrobenthic assemblage in the innermost part of Laki was characterized by a higher total density and the occurrence of a few species which are exclusive inhabitants of coastal brackish habitats. These species have marine (e.g. S. shrubsolii, C. orientale, V. maritima, A. ovata, C. glaucum, H. diversicolor) or limnic origin (tubificid oligochaetes) (e.g. Diviacco 1983; Gravina et al. 1988; Kevrekidis and Koukouras 1988; Sarda and Martin 1993; Kevrekidis et al. 1996; Gouvis et al. 1998). Species composition was quite similar to that in other brackish habitats of the Evros Delta which have a moderate access to the sea (e.g. Kevrekidis 1997) and assemblage characteristics (a limited number of species, a strong dominance in abundance by a few of these species, relatively low diversity) were similar to those previously reported for other Mediterranean brackish habitats (e.g. Reizopoulou et al. 1996; Mistri et al. 2001). Cognetti and Maltagliati (2000) stated that the same euryhaline species usually constitute a typical brackish community, particularly in habitats in the same biogeographical region; their occurrence depends on basin area and on their individual tolerance to variations in some abiotic factors other than salinity (e.g. temperature, oxygen availability, substrate).

The set of distance from the sea, depth, salinity and organic matter of the sediment was found to correlate best with the spatial distribution pattern of the macrobenthic fauna in Laki. This suggests that confinement results from the combined effects of these variables. Koutsoubas et al. (2000) suggested that confinement is determined by a few variables, the extreme values of which produce a gradient along a lagoon and act as thresholds to species distribution.

Macrozoobenthic community structure in the innermost part of Laki lagoon showed seasonal variation. Its dynamics was mainly governed by water temperature. Temperature, along with photoperiod, provide a seasonal cue for the timing of reproduction of macrozoobenthic species and therefore for community changes (e.g. De March 1977; Whitlatch 1977; Moore 1981). The available information on the life history of some of the frequent species in the assemblage, such as A. ovata and G. aequicauda (e.g. Kevrekidis and Koukouras 1988, 1989, 1992), indicates that the observed seasonal variation in density of these species is mainly a result of the timing of their reproductive cycle.

The macrobenthic fauna in the innermost part of Laki lagoon was composed of only a few brackish-water species during the late winter/spring 1998 period at salinities of 0.1–2 psu. Extremely low salinities over a long period were possibly responsible for the decline in most community indices during spring. Species composition (C. orientale, S. shrubsolii, H. diversicolor, V. maritima, A. ovata, G. aequicauda, etc.) was quite similar to that in early autumn 1997 at 32–35 psu. This, along with the increase in total density during spring, reveals the euryhaline character of most community species. Most of the collected species have rarely been reported with this high level of density from habitats sustaining such low salinities for long periods (e.g. Cunha and Moreira 1995).

Some kind of stress most probably resulted in the impoverishment of the macrobenthic fauna during summer. That impoverishment was characterized by the disappearance of some taxa (Tubificidae, Cumacea) or the decline in density of some species (A. ovata, V. maritima); total density also decreased and community indices showed comparatively low values during summer. This disturbance was probably due to a dystrophic episode. The lowest dissolved oxygen concentrations during daytime (5.3 mg/l) were observed in June, probably due to a rapid decomposition of macroalgae. Moreover, anoxic/hypoxic conditions possibly occurred from late night to early morning during summer, when respiratory oxygen consumption by the increased number of macroalgae was probably not compensated for by photosynthetic oxygen production (e.g. Cunha and Moreira 1995; Lardicci et al. 1997; Tagliapietra et al. 1998).

The recovery of the macrofauna in the innermost part of Laki lagoon started in late summer. Furthermore, marine species (M. maculatus, P. ciliata, G. tridactyla, A. cerrutii) occurred in the innermost part of the lagoon during late summer, autumn and winter, probably as a result of immigration or transportation of larvae from the outer part of the lagoon or from marine environments surrounding the lagoon. That occurrence followed, with a time lag, the increase in salinity above 5 psu during summer and autumn, and was also reflected in an additional increase in community indices. Most of the aforementioned marine species were also previously found in brackish habitats (e.g. Gravina et al. 1988; Kevrekidis and Koukouras 1988; Nicolaidou and Papadopoulou 1989; Thiermann et al. 1997; Gouvis et al. 1998; Mistri et al. 2001), while the 5 psu salinity limit is considered as the critical threshold for the occurrence of marine species in these habitats (Cognetti and Maltagliati 2000).

The occurrence of marine species in the innermost part of Laki suggests a shift of the limits of the lagoonal zones during the recovery phase. Koutsoubas et al. (2000) also observed that the limits of the zones in another Mediterranean lagoon (Gialova Lagoon, Ionian Sea) show seasonal shifts; they suggested that the gradient produced by the extreme values of a few environmental variables along a lagoon disappears during the recovery phase.

Predation pressure by A. boyeri possibly contributed to the decline in density of some invertebrates (e.g. G. aequicauda, V. maritima, C. glaucum) during late autumn and winter and in total density in September and November as well. It is known that this small fish, which was collected only in September and November in Laki, feeds on macrobenthic invertebrates (Quignard and Pras 1986). Many authors have stressed the role of predation by fish and decapods (top-down control) in regulating the dynamics of macrobenthic fauna (e.g. Wiltse et al. 1984; Aarnio and Bonsdorff 1993; Sarda et al. 1996). However, the abundance of no other epibenthic decapod or small fish displayed a strong correlation with community structure, although the diet of most of them is known to include macrobenthic invertebrates (e.g. Muus 1967; Scalera Liaci et al. 1982; Cottiglia et al. 1983a, 1983b; Khoury 1984; Maugé 1986; Kevrekidis et al. 1990). The increased number of macroalgae during summer and autumn, when both epibenthic predators are most abundant and their feeding activity, especially in summer, is probably high, may offer macrobenthic invertebrates a refuge from predation (Arias and Drake 1994).

The abundance of the macrobenthic invertebrates in the innermost part of Laki lagoon displayed a seasonal variation with a peak in late spring and a lower one in mid autumn. As is typical of Mediterranean coastal lagoons, the macrofauna is exposed to extreme environmental conditions during summer, and macroinvertebrate abundance thus usually shows a seasonal variation with a peak in winter or spring (e.g. Guelorget and Michel 1979a) and sometimes a second one in autumn (e.g. Guelorget and Michel 1979a, 1979b; Gouvis et al. 1997; Mistri et al. 2001). A seasonal variation of total abundance with a peak in mid or late autumn has also been observed (e.g. Guelorget and Michel 1979a; Kevrekidis 1997).