Translator Disclaimer
27 October 2020 Effects of Salinity and Suspended Solids on Tropical Phytoplankton Mesocosm Communities
Genevieve Sew, Peter Todd
Author Affiliations +

Phytoplankton play a fundamental role in marine food webs but are affected by both natural and anthropogenic fluctuations in environmental conditions. Here, to simulate a dynamic coastal environment, we used mesocosms to examine how different salinity levels and suspended solids concentrations (SSCs) impact a natural phytoplankton assemblage collected from a tropical estuary in Singapore. Significant differences in the phytoplankton composition between the baseline and treatments with medium and high SSC were found, but not among the three salinities tested. Differences can be attributed to nutrient limitation (particularly silicate) and the use of kaolinite for the suspended sediment. Silicate limitation is likely to have caused the observed switch in dominant genus from Skeletonema sp. to Chaetoceros sp. and the occurrence of weakly silicified genera such as Cylindrotheca. Kaolinite affected phytoplankton abundance through effects such as shading, flocculation, and nutrient adsorption. These results demonstrated how the combination of various physicochemical effects of suspended solids can influence tropical phytoplankton communities. Furthermore, as suspended solids such as kaolin can be found in the natural environment, this study showed that their potential effects should be evaluated beyond just their concentration.

As primary producers, phytoplankton form the base of most marine food webs, providing approximately 50% of carbon directly and indirectly for almost every aquatic organism (Field et al., 1998). Light is an essential resource for phytoplankton but, due to the dynamic nature of the ocean environment, its intensity varies in space and time. For example, rapid coastline development and associated activities such as dredging and land reclamation frequently result in elevated suspended solid concentrations (SSCs) and turbid conditions (Field et al., 1998). Even though phytoplankton possess the capacity to acclimatize to varying light conditions through mechanisms such as altering their photosynthetic pigment content and ratio (Dubinsky & Stambler, 2009; Nymark et al., 2009) and adjusting the number and size of their photosynthetic units (Falkowski & Raven, 2007), the shading effect created by turbidity can still impact phytoplankton primary productivity (Grobbelaar, 1985), especially within short time scales.

Changes in coastal water salinity occur due to natural causes such as heavy rainfall and terrestrial runoff and/or anthropogenic sources such as at outfalls of power plants. Changes in salinity can affect aspects of the phytoplankton community, such as species composition, due to interspecific differences in salinity optima (D’ors et al., 2016; Lionard et al., 2005; Qasim et al., 1972). Rapid and large salinity changes subject algal cells to osmotic stress (Lionard et al., 2005), disturbing cellular homeostasis. Suboptimal salinities can also lead to the suppression of photosynthesis, carbon fixation and cell division (D’ors et al., 2016), as well as affecting nitrogen metabolism (Dohler & Biermann, 1985), all of which are tied to cell growth rate. As salinity affects growth rate, it plays a role in phytoplankton competition that can lead to changes in community composition (Muylaert et al., 2000). These effects are generally more pronounced in coastal than in oceanic ecosystems, where the larger body of water buffers against salinity fluctuations (Smayda, 1958).

Studies investigating the interacting effects of SSCs and changes in salinity on phytoplankton are generally field-based surveys conducted in estuaries. Estuaries have a natural salinity gradient due to the freshwater input from upriver and saltwater intrusion at the river mouth and usually have concurrent high sediment levels. However, as a consequence of data being collected longitudinally (e.g., Lueangthuwapranit et al., 2011; Muylaert et al., 2009; Zhang et al., 2014), the phytoplankton community sampled analyzed along estuaries are often not independent. The effects of salinity and SSC are also potentially confounded by other factors such water flow (Lueangthuwapranit et al., 2011). Furthermore, the great majority of these studies have been carried out in temperate countries (Cloern et al., 1989; Kromkamp & Peene, 1995; Lionard et al., 2005; Muylaert et al., 2009; Zhang et al., 2014) with very few performed in the tropics (Lueangthuwapranit et al., 2011; Palleyi et al., 2008; Rochelle-Newall et al., 2011). To better understand the interaction effects of salinity and SSC on coastal phytoplankton communities, and minimizing influence from other factors, it is important to conduct manipulative experiments under controlled conditions.

The aim of our study was to determine the effects of SSC and salinity on a phytoplankton community collected from an estuary in tropical Singapore. This was achieved through testing three SSCs and three different salinity levels in combination using mesocosms. We hypothesized that the treatments with the highest concentration of SSC and with a change in salinity level from ambient would have the lowest phytoplankton density and species richness.


Experimental Design

The experiment was conducted outdoors, at ambient temperature and light intensities. The phytoplankton community was collected from Pandan River (1°18′9.29″N, 103°45′15.33″E) an estuarine water body with mean salinity of ∼30.0 ‰ at the collection site. A total of three salinities (Pandan River salinity and ±3‰) and three sediment concentrations (5 mg l−1, 20 mg l−1, and 80 mg l−1), were tested in a fully crossed design (n = 3), bringing the total number of experimental units (mesocosms) to 27 (Figure 1). Kaolin powder (KM20, Kaolin (Malaysia) Sdn Bhd, Appendix, Table A1) was used to simulate the physical effects of SSCs within the experimental units.

Figure 1.

Randomized Layout of Experimental Units and Treatments. PAR Loggers Were Immersed Within Selected Experimental Units (Shaded Gray). SSC Treatments →L = Low; M = Medium; H = High; Salinity Treatment→ D = Decrease; N = No Change; I = Increase. Note That the Units Were Placed in a Single Line on Site.


Each experimental unit was filled with 20 L of ultraviolet radiated, filtered (0.2 µm) natural seawater from Pandan River. Other than introduced kaolin, the chemical components (e.g., TN, TP) were assumed to remain unchanged. The desired salinity (+3‰, no change and −3‰) treatments were achieved through the addition of Dionised water and aquarium salt (H2Ocean pro+) to the designated units.The treatments of ±3 psu were chosen with the aim of simulating a scenario of a heavy downpour buffered by an adjacent ocean body, a common occurrence in the coastal environments of Singapore. The experimental units were homogenized via manual stirring with a paddle and allowed one day to stabilize, after which the salinity was checked again. Natural phytoplankton populations were collected from Pandan River at 0.5 m depth during high tide (12:00 hr), using a 5 L niskin water sampler (General Oceanics) in triplicates. Immediately, subsamples were taken for the analysis of SSC, after which the remaining water was filtered through a 35 µm mesh to remove larger zooplankton. Further subsampling was then carried out for water quality analysis of other parameters (total nitrogen, total phosphorus, nitrates, nitrites, ammoniacal nitrogen, silicate, and chlorophyll a) and phytoplankton composition and abundance. Water quality samples were processed and analyzed according to the respective methods of the American Public Health Association (APHA, 2012; Appendix, Table A2). After subsampling, 100 ml of the remaining water was then inoculated into each experimental unit. The first replicate collected from Pandan River was used to inoculate the first replicate for each of the respective treatment combinations, and the second was used for second replicate and so forth. In total, 900 ml was taken from each of the 5 L replicates. From the time when samples were taken from Pandan River for inoculation, the whole experiment was established within an hour. The experimental units were distributed randomly using a random number table, so as to minimize the effect of location, and sheltered with clear plastic sheeting to prevent rain affecting salinities. Five loggers measuring photosynthetic active radiation (PAR; Odyssey, Dataflow Systems Ltd) were distributed among the 27 units to capture light distribution across the whole setup and to obtain PAR measurements from the different SSC treatments: two were placed in the low suspended solid treatment and one in the medium and two in the high SSC treatment. Temperature loggers (HOBO, Onset Computer Corporation) were placed within each unit. The experiment ran for a total of seven days, during which experimental units were gently stirred hourly between 09:00 and 17:00 hr with separate paddles to minimize cross contamination. The duration of seven days was chosen to identify the short-term changes that are characteristic of a dynamic nearshore environment. The salinity, temperature, and dissolved oxygen were checked at the beginning of each day using a portable YSI multimeter (YSI, Professional Plus). The experiment was terminated on the seventh day, and water samples were taken for analysis of nutrients (total nitrogen, total phosphorus, and silicate), chlorophyll a, and phytoplankton composition and abundance. Data from the PAR and temperature loggers were then extracted. Samples for phytoplankton composition and abundance were preserved with acidic lugol. All water samples were taken back to the laboratory (DHI Singapore Pte Ltd) in cooler boxes to minimize the influence of temperature changes. Water quality samples were analyzed as described earlier. Phytoplankton samples were enumerated and identified to the lowest possible taxon using the Utermöhl method (Utermöhl, 1931, 1958) with an inverted microscope (Olympus, CKX41) at 200× magnification and assessed for phytoplankton composition.

Statistical Analysis

Statistical analyses were carried out in R (R Core Team; Pinheiro et al., 2017). All data were checked for normality and variance. If the data were found to be normally distributed, with no significant difference in variance, an analysis of variance (ANOVA) was performed. Pairwise multiple comparisons were made using Tukey’s honest significant differences (Tukey’s HSD). If the data did not fulfill either of the assumptions, a linear mixed effects model (lme) (Pinheiro et al., 2017) or Kruskal–Wallis rank sum test was applied. For LME, all models were checked for heteroscedasticity and, if present, variance structures were input and the model with the lowest Akaike information criterion (AIC) was chosen. Pairwise multiple comparisons were then carried out using the “glht” function from the “multcomp” package (Hothorn et al., 2008).

Relationships between phytoplankton and environmental data were analyzed in Primer (7 with PERMANOVA+, Plymouth, UK) using canonical analysis of principle coordinates (CAP). For the CAP, the abundance of individuals within each genus was used as the input variable for each replicate. Due to the large range in plankton counts and the presence of a few highly abundant genera (e.g., Chaetoceros sp.), the data set was square root transformed to reduce the influence/bias of these relatively abundant genera when carrying out Bray–Curtis dissimilarities (M. J. Anderson & Willis, 2003). Following this, a distance matrix was created from Bray–Curtis similarity values. The CAP was then carried out on this distance matrix. The CAP plot was then overlaid with the phytoplankton genera and the water quality parameters (correlation >0.4) separately.


Before the start of the experiment, samples from the Pandan River (PR; hereafter referred to as the baseline) were analyzed to assess water conditions. The mean salinity (±SE) was 30.06 ±0.03‰, SSC concentration was 14.31 ±0.66 mg l−1 and nutrients averaged 0.026 ±0.002 for µM for total nitrogen, 0.0006 ±0.00005 µM for nitrate, nitrite was below the detection limit (<0.0002 µM or <0.01 mg l−1) in all samples, 0.0069 ±0.0003 µM for ammoniacal nitrogen, 0.00179 ±0.00008 µM for total phosphorus, and 0.00025 ±0.00001 µM for phosphate. Chlorophyll a had a mean concentration of 3.60 ±0.30 µg L−1 and silicate 0.0130 ±0.0004 µM (Table 1).

Table 1.

Baseline (Pandan River, PR) Water Quality Parameters (n = 3).


Baseline (i.e., initial) phytoplankton samples had a total of 27 genera which were grouped into four categories (Diatoms, Dinoflagellates, Cyanobacteria, and Others; Table 2). The most abundant group of organisms was Diatoms followed by Dinoflagellates. The mean cell abundance was 231 cells ml−1, of which, the dominant genus was Skeletonema sp. with a mean abundance of 181 cells ml−1, followed by Pseudo-nitzschia sp. at 14 cells ml−1. The genera with the lowest mean abundance were from Coscinodiscus sp., Karenia sp., and Licmophora sp. at 0.1 cells ml−1 each. The mean number of genera collected in baseline samples was 17.3 (Figure 2).

Table 2.

Mean Phytoplankton Composition (Cells ml−1) for the Baseline and the Respective SSC Treatments.


Figure 2.

Mean (± SE) Genera Richness Recorded at the Baseline (Before) and in the Three Different SSC Treatments at the End of the Experiment (Seven Days Later). TSS = Total Suspended Solids.


The experiment ran for seven days, April 7–14, 2015. On the second day, four experimental units from different treatments were removed due to disturbance (building maintenance). At the end of the experiment, water temperature across the days had a mean (± SE) of 27.6 ±0.05°C, with a range of less than 2°C. Dissolved oxygen had a mean concentration of 4.17 ±0.44 mg l−1 and daytime PAR levels ranged from a mean of 968 to 1,490 µmol photon m−2 s−1. The total number of observed phytoplankton genera had increased from 27 (baseline data) to 35. The dominant group remained diatoms (Table 2). Among the SSC treatments, the medium SSC concentration had the highest mean (± SE) number of genera (16.1 ±1.1), whereas the high SSC treatment had the lowest mean number of genera (13.9 ±0.5). There was no significant difference in the genera richness across all treatments (F = 2.011, df = 3, p = 0 .142). Among the genera documented, Chaetoceros sp. had the highest mean abundance at 6,695 cells ml−1 followed by Skeletonema sp. with 637 cells ml−1, both from the treatment with low SSC. The least abundant genera were Cyclotella sp. from the medium SSC treatment (0.1 cells ml−1) and Oscillatoria sp. from the low SSC treatment (0.1 cells ml−1). Comparing across the SSC treatments, the low SSC treatment had the highest mean abundance with 7,579 cells ml−1, whereas the medium SSC treatment had the lowest mean abundance at 626 cells ml−1. Phytoplankton composition between SSC treatments was significantly different, whereas there was no significant difference in phytoplankton composition between salinity and no significant two-way interaction between SSC and salinity (Table 3). The baseline phytoplankton composition was significantly different from the medium (t = 1.976, p = 0.005) and high (t = 1.923, p = 0.015) SSC concentration treatments. For the low SSC treatment, the high abundance was largely due to Chaetoceros sp. In the medium SSC treatment, Leptocylindrus sp. and Asterionellopsis sp. were the two most dominant genera. In the high SSC treatment, Chaetoceros sp. was the most dominant genus. Also, all genera that were more abundant in the medium than in the low SSC treatment declined in the high SSC treatment (Asterionellopsis: 85 to 8 cells ml−1; Leptocylindrus: 203 to 133 cells ml−1; Navicula: 25 to 20 cells ml−1; Karlodinium: 30 to 19 cells ml−1; Protoperidinium: 10 to 4 cells ml−1), except for Cylindrotheca sp. As phytoplankton composition was not significantly different among salinity treatments, further analysis was conducted on the pooled data.

Table 3.

PERMANOVA Results Testing for the Effects of TSS and Salinity for Phytoplankton Genus Composition.


The CAP analysis produced 25 PCO axes, of which the first explained 79.8% of the variability and achieved 42.3% correct allocations (Figure 3). The phytoplankton composition had two distinct groupings, that is, the baseline and the experimental units, which was associated with silicate concentration. Among the experimental units, the spread of points to higher values on Axis 2 was due to high chlorophyll a concentration, particularly in treatments with low SSC which had a greater abundance of Chaetoceros sp.

Figure 3.

CAP Plot of Phytoplankton Composition for the Baseline Samples and Experimental Samples at Different SSC Levels at the End of Experiment. Overlayed With Phytoplankton Genera (Dotted Line) and Water Quality Parameters (Solid Lines, Bold Text) Which Have a Correlation of >.40. The Overlays Were Generated Separately but Combined in This Figure for Easy Viewing. TSS = Total Suspended Solids; CAP = Canonical Analysis of Principle Coordinates.



Plankton are known to be sensitive to physicochemical changes in their marine environment (Hays et al., 2005). Our study used mesocosms to examine the independent effects of three concentrations of suspended solids and three salinity levels on a natural phytoplankton community. The greatest difference in phytoplankton composition was found between the baseline and the experimental units with medium and high SSC concentrations. Among the three SSC treatments, phytoplankton abundance was the highest in the low SSC treatment, potentially due to shading effects in the medium and high treatment. Also, the change in the number of phytoplankton genera documented was likely due to their low initial density, resulting in them not being captured during the first analysis (i.e., baseline samples). Salinity did not appear to have any effect on plankton numbers or composition.

The different SSC concentrations across the experimental units affected phytoplankton abundance, but not genera richness. The SSC levels in the low treatment were lower than the baseline (i.e., 5 mg l−1 vs. 14.31 mg l−1) determined from the initial samples, whereas the SSC levels for the medium and high treatments (20 mg l−1 and 80 mg l−1) were higher than the baseline. The higher SSC in the medium and high treatments probably created a shading effect (Kirk, 1985), reducing light penetration through the water in the experimental units. Light availability is a critical control of phytoplankton dynamics in estuarine phytoplankton (Cloern, 1987), for example, Dokulil et al. (1994) found that increasing SSC in the River Danube resulted in an exponential decrease in chlorophyll a (biomass). Shading can also result in reduced photosynthetic rates, impacting primary productivity and phytoplankton growth (Bilotta & Brazier, 2008; Cloern, 1987).

In addition to causing a shading effect, the use of kaolin in the experiments could have had some unintended influence on the phytoplankton density and nutrient availability within the experimental units. Natural clays such as kaolin have been used in various countries such as China, Japan, and South Korea as a cost-effective method for controlling harmful algal blooms (Table 1 in Sengco & Andersen, 2004). The clays are applied to the water column and through co-flocculation, sediment algal cells out of the water column ( Anderson, 1997; Yu et al., 1995). Cell removal efficiencies as high as 95% to 99% have been reported (Yu et al., 1995). Regular homogenization of the units during the experimental period may have encouraged flocculation, forming larger particles. These larger particles would have caused cells to quickly settle out from the water column before they could be collected. Thus, the higher concentrations of kaolin (medium and high SSC) would have not only have created a turbidity regime but may have caused cells to sediment out of the water column, possibly explaining the high phytoplankton abundance in the low SSC treatment. Kaolin has also been found to be capable of adsorbing metal ions (Duursma & Eisma, 1973; Noriki et al., 1985) and nutrients (Lu et al., 2015), some of which (e.g., iron, total nitrogen, total phosphate) are essential to the phytoplankton growth (Sunda, 1989). Trace metals were not tested for in this study; however, treatments with medium and high suspended sediment concentration had lower TN and TP concentrations than the treatment with low suspended sediment indicating decreased availability of nutrient uptake for the phytoplankton (Appendix, Table A3).This decreased availability of nutrients could then have further stressed the phytoplankton community. Unexpectedly, the average phytoplankton abundance was lowest in the medium SSC treatment, instead of the high SSC treatment. The kaolinite which was used to simulate the SSC is mostly made up of silicate and aluminium, thus although the higher concentration of SSC caused turbid conditions and may have adsorbed some nutrients, the silicate and aluminium within the kaolinite could have contributed toward phytoplankton growth.

The decrease in silicate concentration during the experiment (while other environmental factors such as temperature, dissolved oxygen, and pH levels varied minimally) was a strong factor separating the phytoplankton composition of the baseline and experimental units. By the end of the experiment, the mesocosms had a marked decrease (on average 20 times lower) in silicate concentration compared to the baseline, which may have resulted in the observed differences in phytoplankton composition. By Day 7, the experimental units were silicate limited (final: TN:Si = 100–5.6:1), whereas water quality measured in baseline samples were nitrogen limited (baseline: TN:Si = 0.46:1). The dominant phytoplankton genus in the baseline samples was Skeletonema sp. However, at the end of the experiment, the dominant genus in the experimental units for the treatments with low and high SSC concentration was Chaetoceros sp. Harrison and Davis (1979) observed a similar switch in dominant species with outdoor continuous cultures of a natural phytoplankton assemblage. Skeletonema costatum dominated under nitrogen limitation, but this changed to being dominated by Chaetoceros spp. with silicate limitation (Harrison & Davis, 1979). Silicon-limiting concentrations may also result in physiological changes such as the formation of thinner frustules on new cells (Paasche, 1973, 1975). This was not verified directly, but we did observe a higher mean abundance of Cylindrotheca sp., which is a weakly silicified diatom (Harrison & Davis, 1979), in all the SSC treatments.

Diatoms (Phylum: Ochrophyta) remained dominant in all the units despite silicate limitation. Usually, silicate limitation selects for nonsiliceous flora such as flagellates (Egge & Aksnes, 1992; Weston et al., 2008) and cyanobacteria (Rocha et al., 2002). Egge and Aksnes (1992), using floating outdoor sea enclosures, identified a minimum silicate concentration of 2 µM for diatom dominance. Both Egge and Aksnes (1992) and Weston et al. (2008) reported an increase in the flagellate Phaeocystis sp. in silicon-limited environments. With the exception of one unit from the treatment LI (3.49 µM), silicate concentration in the experimental units ranged from 0.17 to 1.45 µM. Besides limiting silicate concentrations, other criteria such as nutrient ratios and light have been reported in the literature for this diatom-flagellate switch. Zhou et al. (2017) found that, aside from a low Si:N ratio (< 1), a high N:P ratio (> 50) accelerated the bloom of the dinoflagellate, Prorocentrum donghaiense. Nitrogen concentrations in our experimental units remained low, and TN:TP ratios in the experimental units only reached a maximum of 27.3. In terms of light, flagellates have a higher requirement than diatoms (Peperzak et al., 1998). Weston et al. (2008) found that, in the Thames estuary, UK, the switch to flagellate dominance was brought about by a combination of limiting silicon concentration and higher irradiance during the spring. The (dynamic) natural light that the mesocosms were exposed to may have slowed the competitive exclusion as compared to a constant light regime. Fluctuations in light cause recurring interruptions that might reverse successional processes. These interruptions create opportunities for different phytoplankton species to occur, leading to a higher diversity and impeding competitive exclusion (Flöder et al., 2002). For example, cloud cover due to rainy weather on the last two days of our experiment would have caused a reduction in available light. The light-dependent uptake of nutrients such as nitrogen and phosphate require different optimal nutrient uptake strategies under light and dark conditions (Litchman, 2003; Litchman et al., 2004). At the time of final sampling, the experimental units may have been in a time-lag occurring before the (expected) sharp increase of flagellate abundance (Escaravage & Prins, 2002).

Regular stirring of the experimental units to maintain the kaolin in suspension may have negatively affected the growth of dinoflagellates. Among the various phytoplankton groups, dinoflagellates are the most sensitive to turbulence (Thomas & Gibson, 1990). Agitation of the water column has been found to have negative effects on dinoflagellates, some of which include disturbance of cell division, morphological changes, cell disruption and behavior interference (Estrada & Bardalet, 1997). Water column stability is needed for dinoflagellate growth (Peperzak et al., 1998), but diatoms can survive in well-mixed environments (Escaravage & Prins, 2002).In calm environments, diatoms sediment out of the surface layer, whereas more motile organisms such as dinoflagellates are able to move vertically and remain close to the surface. Mixing resuspends nonmotile cells such as diatoms and also disrupts gradients (e.g., light and nutrients) creating a homogenous environment (Hinder et al., 2012). Diatoms have a higher growth rate, higher photosynthetic rate, and are capable of surviving at low nitrogen concentrations (Ross & Sharples, 2007); hence, they may have been able to outcompete dinoflagellates under the experimental conditions we created.

The ability of species to acclimatize to salinity changes could explain the lack of significant difference among the phytoplankton communities subjected to the different salinity treatments. Often, variation in the phytoplankton community along salinity gradients can be attributed to the fact that most species present are stenohaline (Brand, 1984; Larson & Belovsky, 2013). Sudden changes in salinity can cause osmotic stress and stenohaline organisms are be particularly vulnerable toward fluctuations in salinity (Brand, 1984). The baseline samples were collected at a river mouth that opens into the Singapore Straits. Freshwater input from drainage upstream due to heavy downpours during monsoon seasons along with overflow from the adjacent Pandan Reservoir will cause a decrease in salinity. Conversely, seawater intrusion due to tidal forces can potentially cause localized increases in salinity. Thus, phytoplankton species living in this estuarine environment are likely to have had some level of salinity tolerance (Brand, 1984). Skeletonema sp. was the most abundant genus recorded in the baseline samples, and a study on 10 strains of Skeletonema showed that all displayed euryhaline responses, with those from estuarine populations possessing the greatest range of tolerance (Balzano et al., 2011).

This experiment was carried out at a mesocosm scale, in a simplified environment, and in the absence of larger predators (>35 µm). In the natural environment, currents, tides, and runoff are potential avenues by which nutrients for coastal phytoplankton populations can be replenished. Also, connectivity to other aquatic bodies allows for the renewal of phytoplankton cells that may have been lost to processes such as mortality or sedimentation. The exclusion of predators removed the effect of grazing pressure on the phytoplankton community. Microzooplankton (<200 µm) have been shown to be key consumers of primary production (Calbet & Landry, 2004), and their grazing can significantly affect the phytoplankton community, the number of phytoplankton genera, and their abundance.

Using a manipulative mesocosm-based approach, this study tested the influence of salinity stress and SSC concentrations in combination on a natural phytoplankton community. Overall, salinity had little effect on the response variables, however, the SSC treatments probably resulted in poor nutrient and light conditions, such as suboptimal nutrient ratios and low irradiance, impacting the phytoplankton abundance and affecting community composition. Higher concentrations of kaolinite increased turbidity, leading to low irradiance thus affecting the phytoplankton primary productivity. Furthermore, the use of kaolinite likely had some additional impacts such as flocculation and nutrient absorption. Among the nutrients analyzed, silicate limitation was a prominent factor affecting the phytoplankton community, distinguishing the baseline and experimental units.

Implications for Conservation

Even though phytoplankton perform vital functions in marine ecosystems, these roles are often overlooked when designing coastal conservation strategies. Phytoplankton are sensitive to changes in their environment, but estuaries and nearshore waters are often impacted by high levels of sediment pollution (Barton et al., 2016; Hays et al., 2005). This is likely to be compounded by climate change, especially extreme weather events that can lead to lowered salinity and increased turbidity from terrestrial run off. Although salinity in this study had no significant effect, larger changes have the potential to alter phytoplankton communities (Brand, 1984) and suboptimal salinities for resident species may allow opportunistic species to thrive (Guinder et al., 2010).

Three concentrations of kaolinite were used in this study, and changes in the phytoplankton community at higher concentrations were possibly due to various factors working in combination (e.g., nutrient limitation, cell adhesion). As shorelines undergo development and coastal cities expand, excess suspended solids due to sediment disturbance and terrestrial runoff are introduced into the marine environment. Clays similar to kaolinite occur naturally and elevated concentrations in the water column could lead to observations similar to this study. Being the foundation of most aquatic food webs, potential impacts to the phytoplankton community (e.g., changes in species composition) should not be overlooked, as these can have repercussions for higher trophic levels (Butler, 1995). Overall, this study has shown that the impacts of SSC are multi-faceted with co-varying factors that can make general predictions of sediment effects on phytoplankton challenging to formulate.


Table A1.

Kaolin Particle Size Distribution Used in the Experiments.


Table A2.

Schedule of Samples Taken for Each Water Quality Parameter.


Table A3.

Mean (± SE) Concentration of Water Quality Parameters Tested From the Experimental Units at the End of Seven Days for Each SSC Treatment Against the Baseline.



The authors thank DHI Singapore Pte Ltd for the use of their facilities and also members of the Experimental Marine Ecology Laboratory at the National University of Singapore for their input and assistance.

Declaration of Conflicting Interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.


The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was supported by the Economic Development Board—Singapore and DHI Water Environment Pte Ltd though the EDB-IPP scholarship, grant number: R-154–000-662–592 and the National Research Foundation, Prime Minister’s Office, Singapore under its Marine Science Research and Development Programme (Award No. MSRDP-05).



American Public Health Association. (2012). Standard methods for the examination of water and wastewater (22nd ed.). Port City Press. Google Scholar


Anderson, D. M. (1997). Turning back the harmful red tide. Nature, 388(6642), 513–514. Google Scholar


Anderson, M. J., Willis, T. J. (2003). Canonical analysis of principal coordinates: A useful method of constrained ordination for ecology. Ecology, 84(2), 511–525 Google Scholar


Balzano, S., Sarno, D., Kooistra, W. H. C. F. (2011). Effects of salinity on the growth rate and morphology of ten Skeletonema strains. Journal of Plankton Research, 33(6), 937–945. Google Scholar


Barton, A. D., Irwin, A., J., Finkel, Z., Stock, C. A. (2016). Anthropogenic climate change drives shift and shuffle in North Atlantic phytoplankton communities. Proceedings of the National Academy of Sciences of the United States of America, 113(11), 2964–2969. Google Scholar


Bilotta, G. S., Brazier, R. E. (2008). Understanding the influence of suspended solids on water quality and aquatic biota. Water Research, 42(12), 2849–2861. Google Scholar


Brand, L. E. (1984). The salinity tolerance of forty-six marine phytoplankton isolates. Estuarin, Coastal and Shelf Science, 18(5), 543–556. Google Scholar


Butler, N. M. (1995). Effects of sediment loading on food perception and ingestion by freshwater copepods. Marine and Freshwater Behaviour and Physiology, 26(2–4), 287–294. Google Scholar


Calbet, A., Landry, M. R. (2004). Phytoplankton growth, microzooplankton grazing, and carbon cycling in marine systems. Limnology and Oceanography, 49(1), 51–57. Google Scholar


Cloern, J. E. (1987) Turbidity as a control on phytoplankton biomass and productivity in estuaries. Continental Shelf Research, 7(11/12), 1367–1381. Google Scholar


Cloern, J. E., Powell, T. M., Huzzey, L. M. (1989). Spatial and temporal variability in South San Francisco Bay (USA). 11. Temporal changes in salinity, suspended sediments, and phytoplankton biomass and productivity over tidal time scales. Estuarine, Coastal and Shelf Science, 28(6), 599–613. Google Scholar


D’ors, A., Bartolomé, M. C., Sánchez-Fortún, S. (2016). Repercussions of salinity changes and osmotic stress in marine phytoplankton species. Estuarine, Coastal and Shelf Science, 175, 169–175. Google Scholar


Dohler, G., Biermann, I. (1985). Effect of salinity on 15N-ammonia and 15N-nitrate assimilation of bellerochea yucatanensis and thalassiosira rotula. Biochemie Und Physiologie Der Pflanzen, 180(8), 589–598. Google Scholar


Dokulil, M. T. (1994). Environmental control of phytoplankton productivity in turbulent turbid systems. Hydrobiologia, 289(1–3), 65–72. Google Scholar


Dubinsky, Z., Stambler, N. (2009). Photoacclimation processes in phytoplankton: Mechanisms, consequences, and applications. Aquatic Microbial Ecology, 56, 163–176. Google Scholar


Duursma, E. K., Eisma, D. (1973). Theoretical, experimental and field studies concerning reactions of radioisotopes with sediments and suspended particles of the sea. Part C: Applications to field studies. Netherlands Journal of Sea Research, 6(3), 265–324. Google Scholar


Egge, J. K., Aksnes, D. L. (1992). Silicate as regulating nutrient in phytoplankton competition. Marine Ecology Progress Series, 83, 281–289. Google Scholar


Escaravage, V., Prins, T. C. (2002). Silicate bioavailability, vertical mixing and grazing control of phytoplankton blooms in mesocosms. Hydrobiologia, 484(1/3), 33–48. Google Scholar


Estrada, M., Bardalet, E. (1997). Phytoplankton in a turbulent world. Scientia Marina, 61, 125–140. Google Scholar


Falkowski, P. G., Raven, J. A. (2007). Aquatic photosynthesis (2nd ed., p. 484). Princeton University Press.Falkowski, P. G., Raven, J. A. (2007). Aquatic photosynthesis (2nd ed., p. 484). Princeton University Press.


Field, C. B., Behrenfeld, M. J., Randerson, J. T., Falkowski, P. (1998). Primary production of the biosphere: Intergrating terrestrial and oceanic components. Science (New York, N.Y.), 281(5374), 237–240. Google Scholar


Flöder, S., Urabe, J., Kawabata, Z.-I. (2002). The influence of fluctuating light intensities on species composition and diversity of natural phytoplankton communities. Oecologia, 133(3), 395–401. Google Scholar


Grobbelaar, J. U. (1985). Phytoplankton productivity in turbid waters. Journal of Plankton Research, 7(5), 653–663. Google Scholar


Guinder, V. A., , Popovich, C. A., , Molinero, J. C., & , Perillo, G. M. E. (2010). Long-Term Changes in Phytoplankton Phenology and Community Structure in the Bahía Blanca Estuary, Argentina. Marine Biology, 157(12), 2703–2716. Google Scholar


Harrison, P. J., Davis, C. O. (1979). The use of outdoor phytoplankton continuous cultures to analyse factors influencing species selection. Journal of Experimental Marine Biology and Ecology, 41(1), 9–23. Google Scholar


Hays, G. C., Richardson, A. J., Robinson, C. (2005). Climate change and marine plankton. TRENDS in Ecology and Evolution, 20(6), 337–344. Google Scholar


Hinder, S. L., Hays, G. C., Edwards, M., Roberts, E. C., Walne, A. W., Gravenor, M. B. (2012). Changes in marine dinoflagellate and diatom abundance under climate change. Nature Climate Change, 2(4), 271–275. Google Scholar


Hothorn, T., Bretz, F., Westfall, P. (2008). Simultaneous inference in general parametric models. Biometrical Journal. Biometrische Zeitschrift, 50(3), 346–363. Google Scholar


Kirk, J. T. (1985). Effects of suspensoids (turbidity) on penetration of solar radiation in aquatic ecosystems. Hydrobiologia, 125(1), 195–208. Google Scholar


Kromkamp, J., Peene, J. (1995). Possibility of net phytoplankton primary production in the turbid Schelde estuary (SW Netherlands). Marine Ecology Progress Series, 121, 249–259. Google Scholar


Larson, C. A., Belovsky, G. E. (2013). Salinity and nutrients influence species richness and evenness of phytoplankton communities in microcosm experiments from grate salt lake, Utah, USA. Journal of Plankton Research, 35(5), 1154–1166. Google Scholar


Lionard, M., Muylaert, K., Gansbeke, D. V., Vyverman, W. (2005). Influence of changes in salinity and light intensity on growth of phytoplankton communities from the Schelde river and estuary (Belgium/The Netherlands). Hydrobiologia, 540(1–3), 105–115. Google Scholar


Litchman, E. (2003). Competition and coexistence of phytoplankton under fluctuating light: Experiments with two cyanobacteria. Aquatic Microbial Ecology, 31, 241–248. Google Scholar


Litchman, E., Klausmeier, C. A., Bossard, P. (2004). Phytoplankton nutrient competition under dynamic light regimes. Limnology and Oceanography, 49(4part2), 1457–1462. Google Scholar


Lu, G., Song, X., Yu, Z., Cao, X., Yuan, T. (2015). Effects of modified clay flocculation on major nutrients and diatom aggregation during Skeletonema costatum blooms in the laboratory. Chinese Journal of Oceanology and Limnology, 33(4), 1007–1019. Google Scholar


Lueangthuwapranit, C., Sampantarak, U., Wongsai, S. (2011). Distribution and abundance of phytoplankton: Influence of salinity and turbidity gradients in the Na Thap river, Songkhla province, Thailand. Journal of Coastal Research, 27(3), 585–594. Google Scholar


Muylaert, K., Sabbe, K., Vyverman, W. (2000). Spatial and temporal dynamics of phytoplankton communities in a freshwater tidal estuary (Schelde, Belgium). Estuarine, Coastal and Shelf Science, 50(5), 673–687. Google Scholar


Muylaert, K., Sabbe, K., Vyverman, W. (2009). Changes in phytoplankton diversity and community composition along the salinity gradient of the Schelde estuary (Belgium/The Netherlands). Esturaine, Coastal and Shelf Science, 82(2), 335–340. Google Scholar


Noriki, S., Ishimori, N., Harada, K., Tsunogai, S. (1985). Removal of trace metals from seawater during a phytoplankton bloom as studied with sediment traps in Funka Bay, Japan. Marine Chemistry, 17(1), 75–89. Google Scholar


Nymark, M., Valle, K C., Brembu, T., Hancke, K., Winge, P., Andresen, K., Johnsen, G., Bones, A M. (2009). An integrated analysis of molecular acclimation to high light in the marine diatom, Phaeodactylum tricornutum. PLoS One, 4(11), e7743. Google Scholar


Paasche, E. (1973). Silicon and the ecology of marine plankton diatoms. I. Thalassiosira pseudonana (cyclotella nana) grown in a chemostat with silicate as limiting nutrient. Marine Biology, 19(2), 117–126. Google Scholar


Paasche, E. (1975). Growth of the plankton diatom Thalassiosira nordenskioeldii cleve at low silicate concentrations. Journal of Experimental Marine Biology and Ecology, 18(2), 173–183. Google Scholar


Palleyi, S., Kar, R N., Panda, C R. (2008). Seasonal variability of phytoplankton population in the Brahmani estuary of Orissa. India. Journal of Applied Sciences and Environmental Management, 12(3), 19–23. Google Scholar


Peperzak, L., Colijn, F., Gieskes, W. W. C., Peeters, J. C. H. (1998). Development of the diatom-Phaeocystis spring bloom in the Dutch coastal zone of the North Sea: The silicon depletion versus the daily irradiance threshold hypothesis. Journal of Plankton Research, 20(3), 517–537. Google Scholar


Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D., & R Core Team (2017). _nlme: Linear and Nonlinear Mixed Effects Models_. R package version 3.1-131., J., Bates, D., DebRoy, S., Sarkar, D., & R Core Team (2017). _nlme: Linear and Nonlinear Mixed Effects Models_. R package version 3.1-131.


Qasim, S. Z., Bhattathiri, P. M. A., Devassy, V. P. (1972). The influence of salinity on the rate of photosynthesis and abundance of some tropical phytoplankton. Marine Biology, 12(3), 200–206. Google Scholar


Rocha, C., Galvao, H., Barbosa, A. (2002). Role of transient silicon limitation in the development of cyanobacteria blooms in the Guadiana estuary, South-Western Iberia. Marine Ecology Progress Series, 228, 35–45. Google Scholar


Rochelle-Newall, E. J., Chu, V. T., Pringault, O., Amouroux, D., Arfi, R., Bettarel, Y., Bouvier, T., Bouvier, C., Got, P., Nguyen, T. M. H., Mari, X., Navarro, P., Duong, T. N., Cao, T. T. T., Pham, T. T., Ouillon, S., Torréton, J.-P. (2011). Phytoplankton distribution and productivity in a highly turbid, tropical coastal system (Bach Dang Estuary, Vietnam). Marine Pollution Bulletin, 62(11), 2317–2329. Google Scholar


Ross, O. N., Sharples, J. (2007). Phytoplankton motility and the competition for nutrients in the thermocline. Marine Ecology Progress Series, 347, 21–38. Google Scholar


Sengco, M. R., Anderson, D. M. (2004). Controlling harmful algal blooms through clay flocculation. The Journal of Eukaryotic Microbiology, 51(2), 169–172. Google Scholar


Smayda, T. J. (1958). Biogeographical studies of marine phytoplankton. Oikos, 9(2), 158–191. Google Scholar


Sunda, W. G. (1989). Trace metal interactions with marine phytoplankton. Biological Oceanography, 6(5–6), 411–442. Google Scholar


Thomas, W. H., Gibson, C. H. (1990). Effects of small-scale turbulence on microalgae. Journal of Applied Phycology, 2(1), 71–77. Google Scholar


Utermöhl, H. (1931). Neue wege in der quantitativen erfassung des planktons. (Mit besonderer berücksichtigun des ultraplanktons.) [New ways of quantifying plankton. (With special consideration of the ultraplankton.)]. Verhandlungen der Internationalen Vereinigung für Theoretische und Angewandte Limnologie, 5(2), 567–596. Google Scholar


Utermöhl, H. (1958). Zur vervollkommnung der quantitativen Phytoplankton-Methodik [To perfect the quantitative phytoplankton methodology]. Mitteilungen Internationale Vereinigung Für Theoretische Und Angewandte Limnologie, 9(1), 1–38. Google Scholar


Weston, K., Greenwood, N., Fernand, L., Pearce, D. J., Sivyer, D. B. (2008). Environmental controls on phytoplankton community composition in the Thames plume, U. K. Journal of Sea Research, 60(4), 246–254. Google Scholar


Yu, Z., Zou, J., Ma, X. (1995). Application of clays to removal of red tide organisms III. The coagulation of kaolin on red tide organisms. Chinese Journal of Oceanology and Limnology, 13(1), 62–70. Google Scholar


Zhang, X., Zhang, J., Huang, X., Huang, L. (2014). Phytoplankton assemblage structure shaped by key environment variables in the Pearl river estuary. Journal of Ocean University of China, 13(1), 73–82. Google Scholar


Zhou, Y., Zhang, Y., Li, F., Tan, L., Wang, J. (2017). Nutrients structure changes impact the competition and succession between diatom and dinoflagellate in the East China Sea. Science of the Total Environment, 574, 499–508. Google Scholar
© The Author(s) 2020 Creative Commons Non Commercial CC BY-NC: This article is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License ( which permits non-commercial use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (
Genevieve Sew and Peter Todd "Effects of Salinity and Suspended Solids on Tropical Phytoplankton Mesocosm Communities," Tropical Conservation Science 13(1), (27 October 2020).
Received: 29 November 2019; Accepted: 11 June 2020; Published: 27 October 2020

Back to Top