Corps de l’article

1. Introduction

Estuaries and coastal wetlands are critical transition zones linking terrestrial and marine habitats. Salt marshes are highly dynamic environments and are among the most productive ecosystems in the world (WOODWARD and WUI, 2001). They are characterized by high primary production rates and by intense remineralization processes within the surface sediment (BILLEN and LANCELOT, 1988; POULIN et al., 2007). In addition, they act as natural filters of allochtonous pollutants and actively transform nutrient loads received from the watershed before their introduction into the coastal environment (BALDWIN and MITCHELL, 2000). The environmental significance of these particular ecosystems is now widely recognized and protection efforts have been deployed to reduce impacts of human activities (CLOERN, 2001).

Sediment-dwelling microalgae and phototrophic bacteria (collectively named microphytobenthos hereafter, MPB) play a central role in the functioning of coastal salt marshes through their contribution to stabilize newly deposed sediment (CONSALVEY, 2002) and to support wetland consumers such as benthic macrofauna, meiofauna, fish and birds (GUARINI et al., 2002). Although cyanobacteria are often abundant, diatoms are considered as the most important primary producers in these systems (SULLIVAN and CURRIN, 2000). These diatoms can be divided in two distinct groups: the epipsammic forms attached to particles, and the epipelic free-living forms (ROUND, 1971, 1979). Epipelic diatoms are generally the most important component of the autotrophic biomass (FORSTER et al., 2006) found in these areas. They are affected by water movements and can migrate within the first top centimeters of the sediment (SABUROVA and POLIKARPOV, 2003) allowing the establishment and maintenance of epipelic forms within sedimentary environments (CONSALVEY et al., 2004).

There is a lack of knowledge concerning the density and diversity of microphytobenthic diatoms in northern coastal areas, which limits our understanding of the resilience of these communities exposed to environmental perturbations induced by climate changes and anthropogenic insults. This paper provides initial results on the density and diversity of diatom communities inhabiting two salt marshes located along the south shore of the St. Lawrence Estuary (Eastern Canada) during summer time. The structure of the diatom communities was examined in relation to pigment concentrations, total bacteria abundance, extracellular polymeric substance concentration, sediment grain size distribution and elementary organic matter composition in an attempt to determine the main environmental factors affecting MPB species distribution in such cold environments.

2. Material and methods

2.1 Study sites

Two marshes from the south shore of the Lower St. Lawrence Estuary (Quebec, Canada) were selected for this study: the Pointe-aux-Épinettes (PE) and the Pointe-au-Père (PP) salt marshes (Figure 1). Both located in the same geographical area (Rimouski, Qc), these marshes are affected by similar large-scale physical forcing such as tides, winds, and waves. Lower St. Lawrence Estuary marshes are exposed to important seasonal variations affecting surface water temperature (typically ranging from 0 to 25oC), salinity (from 10 to 30 psu) and nutrient concentrations ([NO3-] from 1 to 25 µM) (POULIN, 2008). At both studied sites, low and high marshes areas are included between -1 and 0 m and 1 and 2 m, respectively (related to the mean sea level).

Figure 1

Map of the St. Lawrence Estuary showing the location of the Pointe-aux-Épinettes (PE) and the Pointe-au-Père (PP) salt marshes. The sampling stations were located in the low marsh zone (PELM, PPLM) and high marsh zone (PEHM, PPHM) of each site.

Carte de l'estuaire du Saint-Laurent présentant la localisation des marais de Pointe-aux-Épinette (PE) et de Pointe-au-Père (PP). Les différentes stations d'échantillonnage sont localisées dans le marais inférieur (PELM, PPLM) et le marais supérieur (PEHM, PPHM) de chacun des sites.

Map of the St. Lawrence Estuary showing the location of the Pointe-aux-Épinettes (PE) and the Pointe-au-Père (PP) salt marshes. The sampling stations were located in the low marsh zone (PELM, PPLM) and high marsh zone (PEHM, PPHM) of each site.

-> Voir la liste des figures

The PE salt marsh is located in the Bic Provincial Park, a protected area created in 1984. The marsh is located at the far end of a protected bay (Anse à l’Orignal) and is protected by small islands, between two steep hills composed of Cambrian claystone/sandstone. The erosion of these geologic formations has promoted the deposition of a large sandy low marsh area. The marsh is dominated by an homogenous emergent macrophyte community (Spartina alterniflora, Spartina pattens and other halophytes in its upper part) surrounded by a creek network free of vegetation. The PE marsh watershed is part of a forested area without agricultural activity and is entirely protected from anthropogenic stressors.

The PP salt marsh is part of the Pointe-au-Père National Wildlife area, a protected area since 2002. The marsh is confined in a shallow elongated basin parallel to the main St. Lawrence Estuary axis and sheltered from wave actions by a low, partially submerged ridge made of Palaeozoic claystone. The lower marsh is dominated by a Spartina alterniflora community growing between numerous boulders buried in muddy sediment. The upper marsh is nearly a flat platform which terminates seaward with a steep slope or a micro-cliff of 0.5 m. Unlike PE, the PP marsh has been subjected for decades to various anthropogenic pressures (urban sewage and agricultural runoff, landfills and road construction) and is still receiving dissolved nutrient inputs from small tributaries (POULIN, 2008). Erosion structures (such as micro-cliff, micro-ravine, high marsh collapse) occur along drainage channel edges. A preliminary survey conducted in May and July 2004 indicated that surface sediment of the PP salt marsh was not heavily impacted with organic pollutants such as chlorinated compounds but aliphatic hydrocarbons and some light polycyclic aromatic hydrocarbons (PAHs) were detected and their presence was attributed to nearby road traffic (PELETIER, E. unpublished results).

2.2 Sampling and analyses

Sediments were sampled in summer 2005 (7-11 July) from the Pointe-aux-Epinettes (48o 21.3’N; 68o 46.5’W) low marsh zone (PELM) and high marsh zone (PEHM), and from the Pointe-au-Père (48o 30.3’N; 68o 28.0’W) low marsh zone (PPLM) and high marsh zone (PPHM) as illustrated in Figure 1. Six cores were randomly collected (according to the benchmark sampling technique) using truncated 10 cm3 sterile syringes in a predetermined square of 10 m2 (benchmark site) and the first two top centimeters were kept for analysis.

Surface sediment samples were processed for diatom analysis following a modification of the Schrader method (1974). A 0.5 g sub-sample of dried sediment was treated with hydrogen peroxide and hydrochloric acid. The samples were washed with distilled water and centrifuged until no acid remained in the samples. An aliquot of 0.5 mL was placed on a cover slip and air dried. Permanent slides were prepared using Hyrax as the mounting medium. A Leitz Wetzlar Orthoplan® microscope, equipped with phase contrast and magnification up to x1000 was used for diatom identification. A minimum of 400 cells were counted in random transects for each sample. Diatom abundance was expressed as number of cells•g-1. The Shannon index H’ (MAGURRAN, 1988) was used for diversity measures. For statistical analysis, tychoplanktonic, periphytic, epiphytic and epilithic cells were included in the epipsammic group whereas planktonic species were included in the epipelic group. Total bacteria abundances (TB) were determined on formaldehyde fixed samples by epifluorescence microscopy using a 4’,6-diamidino-2-phenylindole (DAPI) staining procedure adapted from KUWAE and HOSOKAWA (1999).

Chlorophyll-a (Chl-a) and phaeopigments (Pheo) were first extracted from sediment using 90% acetone solution and analyzed by fluorescence detection (STRICKLAND and PARSONS, 1972). Extracellular polymeric substances (EPS) concentrations were determined by a phenol-sulphuric acid dosage technique (UNDERWOOD et al., 1995) which measures total carbohydrate concentrations. Surface sediment elementary analysis for carbon (Corg) and nitrogen (Ntot) was performed with a ECS 4010 Costech® instrument equipped with a zero blank auto-sampler, using a certified sediment sample (n = 10) NIST-1941b as a reference, whereas grain size analyses were done using a Coulter® counter model LS13320C.

2.3 Statistical analyses

Spearman’s rank correlations were firstly used to test the correlation between diatom community characteristics (i.e. epipelic cells abundance, epipsammic cells abundance, total cells abundance), total bacteria abundances and environmental variables (i.e. Chl-a, Pheo, EPS, Corg, Ntot, Clay, Silt Fine Sand, Medium Sand and Coarse Sand). Thereafter, multiple linear regression analyses were used to link these environmental parameters to the variations of diatom community characteristics considering correlation significance at a p < 0.05 level. Normality and homoscedasticity were confirmed by the examination of the residuals. Analysis of variance (One Way ANOVA on Ranks) was performed to test for significant differences in diatom density, epipelic and epipsammic species as well as for other environmental variables among different sampling sites (PELM, PEHM, PPLM and PPHM). Prior to ANOVA, normality and homoscedasticity were confirmed by the examination of the residuals. To isolate groups that differed from the others, Holm-Sidak and Tukey post-hoc multiple comparison procedures were performed. Differences between groups were considered significant at the p < 0.05 level. Comparisons of means were performed with Sigma Stat® Software Inc.

3. Results

3.1 Environmental variables

Mean (AV ± SD) variable values are presented in Tables 1a and 1b. According to the ANOVAs, sampling sites show significant differences for several variables [Corg (p < 0.001); Ntot (p = 0.002); C/N (p = 0.002); EPS (p < 0.001); Chl-a  p = 0.002); Pheo (p = 0.017), and the five sediment types (p < 0.001)] while mean Chl-a/Pheo ratios did not show any significant inter-site differences (p = 0.611). The a posteriori test showed that Corg, Chl-a, Pheo as well as C/N ratios remained statistically higher in PPHM compared to PELM and PEHM with no significant difference between PPLM and other sampling sites (PPHM, PELM and PEHM). The Ntot concentration was statistically higher at the PPHM and PPLM sites compared to the PELM site, whereas EPS concentrations remained statistically higher at the PPHM site compared to the other sites, except for PEHM. Sediment grain size distribution showed important spatial variability with a significant higher proportion of clay and silt at the PPHM site, a significant higher proportion of fine sand at PPLM and a significant higher proportion of medium and coarse sand at PELM.

3.2 Diatom composition and bacterial abundance

TB abundances ranged from 1.57 x108 cell•g-1 at PELM to 2.48 x109 cell•g-1 at PPHM. TB remained statistically higher at PPHM compared to PELM and PEHM with no significant differences between PPLM and other sampling sites (PPHM, PELM and PEHM).

Table 1a

Mean values (AV) and standard deviations (SD) of sediment characteristics analyzed at stations PELM (Pointe-aux-Épinettes low marsh), PEHM (Pointe-aux-Épinettes high marsh), PPLM (Pointe-au-Père low marsh) and PPHM (Pointe-au-Père high marsh).

Valeurs moyennes (AV) et écarts-types (SD) des caractéristiques géochimiques mesurées aux stations PELM (bas marais de Pointe-aux-Épinettes), PEHM (haut marais de Pointe-aux-Épinettes), PPLM (bas marais de Pointe-au-Père) et PPHM (haut marais de Pointe-au-Père).

Mean values (AV) and standard deviations (SD) of sediment characteristics analyzed at stations PELM (Pointe-aux-Épinettes low marsh), PEHM (Pointe-aux-Épinettes high marsh), PPLM (Pointe-au-Père low marsh) and PPHM (Pointe-au-Père high marsh).

-> Voir la liste des tableaux

Table 1b

Mean values (AV) and standard deviations (SD) of biological variables analyzed at stations PELM (Pointe-aux-Épinettes low marsh), PEHM (Pointe-aux-Épinettes high marsh), PPLM (Pointe-au-Père low marsh) and PPHM (Pointe-au-Père high marsh).

Valeurs moyennes (AV) et écarts-types (SD) des caractéristiques biologiques mesurées aux stations PELM (bas marais de Pointe-aux-Épinettes), PEHM (haut marais de Pointe-aux-Épinettes), PPLM (bas marais de Pointe-au-Père) et PPHM (haut marais de Pointe-au-Père).

Mean values (AV) and standard deviations (SD) of biological variables analyzed at stations PELM (Pointe-aux-Épinettes low marsh), PEHM (Pointe-aux-Épinettes high marsh), PPLM (Pointe-au-Père low marsh) and PPHM (Pointe-au-Père high marsh).

-> Voir la liste des tableaux

A total of 53 diatom species were identified during this study (Table 2). Most of them are benthic species with only six planktonic species observed: Aulacoseira subarctica, Coscinodiscus marginatus, Cyclostephanos dubius, Ondotalla aurita, Thalassiosira eccentrica and T. proschkinae. Among the 47 benthic taxa observed, 24 belong to the epipelic group, 12 to the epipsammic group, 4 are tychoplanktonic, 4 are periphytic, 2 are epiphytic and 1 is epilithic. Except for PPHM, a great similarity in species composition and distribution was observed between stations PELM, PEHM and PPLM. The lowest number of species (26 taxa) was observed at PPLM whereas the highest number (31) was observed at PEHM stations. In terms of diatom density, the lowest values (1.22 x106 cells•g-1) were observed at PEHM whereas the highest one (2.18 x106 cells•g-1) was observed at PELM. The three dominant species were: Achnanthes delicatula (abundance ≥ 50%), Navicula cryptocephala (abundance 10-30%) and Cocconeis disculus (abundance 5-10%). Diatom density was significantly higher at the PPHM station (p < 0.001) with an appreciable difference in the community composition: 44 species were identified and, even though A. delicatula, M. criptocephala and C. disculus remained the most dominant species, each represented less than 30% of total microphytobenthos community. The a posteriori test showed that total cell density was statistically higher at the PPHM site compared to the PELM and PEHM areas, with no significant difference between PPLM and other sampling sites.

Table 2

List of microphytobenthic diatom species and their relative abundance (%) at stations PELM (Pointe-aux-Épinettes low marsh), PEHM (Pointe-aux-Épinettes high marsh), PPLM (Pointe-au-Père low marsh) and PPHM (Pointe-au-Père high marsh) in July 2005.

Liste des espèces de diatomées microphytobenthiques et les abondances relatives (%) dénombrées aux stations PELM (bas marais de Pointe-aux-Épinettes), PEHM (haut marais de Pointe-aux-Épinettes), PPLM (bas marais de Pointe-au-Père) et PPHM (haut marais de Pointe-au-Père).

List of microphytobenthic diatom species and their relative abundance (%) at stations PELM (Pointe-aux-Épinettes low marsh), PEHM (Pointe-aux-Épinettes high marsh), PPLM (Pointe-au-Père low marsh) and PPHM (Pointe-au-Père high marsh) in July 2005.

-> Voir la liste des tableaux

On the basis of their habitat, diatom species were grouped into epipelic and epipsammic forms. Comparison between sampling sites showed significant differences in total diatom density (p < 0.001), epipelic (p < 0.001) and epipsammic (p < 0.001) species with significantly higher total abundance in PPHM and no statistical difference between other sampling sites. The a posteriori test showed the dominance of epipelic forms in PPHM (relative abundance > 57%) compared to PPLM and PEHM whereas epipsammic forms are significantly higher in PPLM compared to all other sites (relative abundance > 87%). Epipelic species are mainly represented by Diploneis interrupta, D. smithii, Gyrosigma fasciola, G. spenceri, Navicula cryptocephala, N. peregrina, N. radiosa var. tenella, N. salinarum, Nitzschia clausii and N. tryblionella while the epipsammic species are mainly represented by Achnanthes delicatula and Cocconeis disculus.

Multiple linear regressions showed that among all sampling sites, the variance of the dependent variable diatom density was mainly affected by a linear combination of surface sediment nutrient concentrations (Table 3). About 72% of the variance in diatom cell density is explained by the combination of Corg and Ntot. Similarly, 81% of the variance of epipelic species abundance could be explained by the combination of Corg and Ntot whereas the best fit model showed that only 33% of the variance of the epipsammic species could be explained by coarse sand relative concentration variability.

Table 3

Results of multiple linear regressions testing the effect of environmental variables on diatom density and abundance of epipelic and epipsammic forms at stations PELM (Pointe-aux-Épinettes low marsh), PEHM (Pointe-aux-Épinettes high marsh), PPLM (Pointe-au-Père low marsh) and PPHM (Pointe-au-Père high marsh).

Résultats des régressions linéaires multiples testant l'effet des variables environnementales sur la densité des diatomées et les abondances des formes épipéliques et épipsammiques aux stations aux stations PELM (bas marais de Pointe-aux-Épinettes), PEHM (haut marais de Pointe-aux-Épinettes), PPLM (bas marais de Pointe-au-Père) et PPHM (haut marais de Pointe-au-Père).

Results of multiple linear regressions testing the effect of environmental variables on diatom density and abundance of epipelic and epipsammic forms at stations PELM (Pointe-aux-Épinettes low marsh), PEHM (Pointe-aux-Épinettes high marsh), PPLM (Pointe-au-Père low marsh) and PPHM (Pointe-au-Père high marsh).

-> Voir la liste des tableaux

The determination of the Shannon and Weaver diversity measures (H’) revealed a higher diversity at the PPHM station (mean H’ = 2.58) compared to mean values of 1.46 bits•cells-1 at the PELM, PEHM and PPLM stations. MPB biomass was estimated using the calculation model proposed by de Jonge (1980). The conversion from µg Chl-a•g-1 dry sediment to mg Chl-a•m-2 was done by multiplying recorded Chl-a concentrations by 12.7, based on an average bulk density of 1.27 g•cm-3, whereas the conversion from Chl-a to carbon was obtained using the C/Chl-a ratio of 40 (DE JONGE, 1980). The estimated surface biomass at PE marsh ranged from 267 to 1762 mg Chl-a•m‑2 (11 to 71 g C•m-2) while it ranged from 576 to 12170 mg Chl-a•m-2 at the PP marsh with the highest values observed at the PPHM site. These values correspond to carbon concentrations ranging from 132 to 486 g C•m-2 at PPHM and from 24 to 65 g C•m-2 at PPLM.

4. Discussion

Our study is the first to describe the summer composition of diatom communities in two northern salt marshes of the Lower St. Lawrence Estuary. Diatom abundance and Chl-a concentrations reported here offer some similarities with those reported in the literature in temperate intertidal areas (PELETIER, 1996; SABBE, 1993). However, significant differences were observed among the four sampling stations, with higher diatom cell density and Chl-a concentration recorded at the PPHM station compared to both PE sampling sites. Although these differences may be related to spatial heterogeneity, the higher biomass registered at the PPHM station seems to result from favorable environmental conditions for MPB growth. Among important variables involved in diatom growth, nutrient load due to anthropogenic activities is generally suggested as the most probable factor leading to an important increase of MPB in coastal areas (HOWARTH and MARINO, 2006). The particular location of the PP saltmarsh, partly enclosed within an urbanized area, might be expected to contribute to the higher diatom abundance observed at the PPHM station as nutrient deposition may be higher in impacted areas compared to pristine ones. In addition, the higher proportion of clay and silt found at this site might promote nutrient retention compared to more sandy sites which favor nutrient runoff toward estuarine water.

Due to their diversity and abundance in aquatic systems, diatoms have frequently been used as a proxy for aquatic ecosystem perturbations (CHRISTIE and SMOL, 1993; UNDERWOOD and KROMKAMP, 1999). One of the well-documented effects of nutrient inputs on diatom communities is the reduction in number of species along with an increase in the abundance of the remaining ones (COOPER, 1995; DE SÈVE and GOLDSTEIN, 1981). A high dominance of epipsammic species within a relatively low diversified community was observed at the PPLM, PEHM and PELM stations, whereas a slight dominance of epipelic forms within a much more diverse community was found at the PPHM station. The higher diversity observed at the PPHM station is attributable to an increase of epipelic species and the presence of some particular species (such as Diploneis interrupta and Gyrosigma acuminatum), suggesting the presence of eutrophic conditions at this station. The differences observed in Chl-a concentrations among different sites could be partially attributed to the proportion of epipelic versus epipsammic diatoms within the community since epipelic forms are generally larger and susceptible to contain more pigments than the epipsammic ones (Underwoodet al., 1995). Although the PP salt marsh can be considered as a human disturbed area, the high MPB diversity recorded at the PPHM station does not reflect the reduction of diatom diversity generally observed in temperate human impacted sites. The absence of toxic industrial wastes and the presence of labile clayey sediment in the PP system might explain the maintenance of high diatom diversity in this marsh.

Organic matter lability, grain size distribution and sediment stability are among marsh sediment characteristics influencing the distribution of microphytobenthic species. In our study, the average C/N ratio ranged from 6.4 to 34.7 among sites, but cannot be used as a relevant proxy to determine the lability of the organic matter due to the high relative carbon content (and low nitrogen) of the Spartina rhizosphere (Spartina C/N ~ 40 following Meyers,1994). For example, Spartina roots at the PPHM station may account for an important proportion of the observed high C/N ratio, but this value does not necessarily imply a refractory character of the interstitial sediment. Thus, high epipelic cell abundances at the PPHM site could be better explained by the microbial degradation potential within the Spartina rhizosphere rather than by the C/N ratio observed in the sediment. Degradation processes estimated through the Chl-a/phaeopigments ratio support this finding. A ratio of 10 is generally observed in growing populations but values below 3 usually indicate some cell degradation or important grazing by meiofauna (CARIOU-LEGALL and BLANCHARD, 1995). Elevated pheopigment concentrations, particularly at station PPHM, further suggest active degradation processes of vegetation. In the present study, the lowest mean Chl-a/phaeopigments ratio was recorded at the PPHM site with no significant differences among the sampling sites (Table 1b) showing strong pigment degradation.

The composition of diatom communities has been previously reported to be influenced by the size of the sediment particles (ZONG and HORTON, 1998). Previous studies in temperate areas showed that epipsammic taxa were mainly associated with sandy substrates whereas epipelic taxa were mainly found in muddy to clayey substrates (FORSTER et al., 2006). Our results in St. Lawrence marshes are in agreement with these findings since the abundance of epipsammic forms can be roughly predicted from coarse sand concentrations (Table 3). In addition, the muddy character of the sediments at the PPHM station (i.e. high clay and silt content; Table 1a) may have contributed to the settling of many epipelic species in this specific area.

The benthic carbon concentrations, derived from Chl-a measures at surface sediment in July 2005, ranged from 11 to 71 g C•m-2 in the PE marsh, and from 24 to 65 g C•m-2 at the PPLM site, with a peak at the PPHM station (>132 g C•m-2). According to these values, the annual potential marsh productivity can be estimated to range from 10 to 70 g C•m-2•y-1. This preliminary calculation may underestimate the marsh production since marsh biomass evaluation (using Chl-a concentrations) does not take in account bacterial and EPS carbon production which could represent up to 30% of the MPB production (HALL and FISHER, 1985). Diatom productivity observed in both studied salt marshes is of the same order of magnitude as the primary production (104 g C•m-2•y-1) recorded in the adjacent biologically rich St. Lawrence Estuary waters (THERRIAULT and LEVASSEUR, 1985).

In conclusion, this work provides the first results on the abundance and diversity of diatom communities inhabiting two northern salt marshes during summer time. Rich organic carbon clayey substrates seem to account for the high diatom counts at the PPHM station along with effluent enrichment, even if no direct links can been established yet. In addition, the relative abundance of epipelic and epipsammic forms is mainly related to the sediment grain size distribution. This study provides reference data from which diatom-based indices (SGRO et al., 2007) could be calculated and used in monitoring salt marsh environmental changes.