The Role of Macrophytes in Biogenic Silica Storage in Ivory Coast Lagoons

Lagoons are shallow aquatic environments that typically show a broad variety in colonization by macrophytes. We present the biogenic silica (BSi) data obtained from 11 macrophyte species randomly collected in three small lagoons (Ono, Kodjoboue


INTRODUCTION
Lagoons are transitional water systems between land and ocean, characterized by a shallow depth (<5 m) and occasional water exchanges with the adjacent ocean (Boynton et al., 1996).They are subjected worldwide to increased nutrient inputs due to anthropogenic activities such as land use change, effluent disposal and aquaculture (Caumette et al., 1996).The resulting eutrophication is characterized by enhanced phytoplankton blooms and increased macrophyte biomass (Sidinei et al., 2001).Macrophytes provide food and shelter for aquatic organisms and regulate nutrient and element cycling dynamics within these systems (Rejmankova, 2011).The distribution of different aquatic macrophyte species and growth forms in lagoons is related in various ways to the physical characteristics of the system (depth, bottom slope, type of sediments, wind regime. . .), to riverine inputs of nutrients, and to the frequency of salt ocean water intrusions (Kjerfve and Magil, 1989;Aloo et al., 2013).
Aquatic macrophytes are also known to store significant amounts of silica (SiO 2 ), thus potentially affecting silica storage, fluxes and turnover rates in freshwater aquatic ecosystems (Schoelynck and Struyf, 2016).Silica in plants is present in an amorphous and hydrated form (SiO 2 •nH 2 O), usually referred to as biogenic silica (BSi) and mostly deposited as phytoliths (siliceous plant bodies) or associated with cell wall molecules (Broadley et al., 2012).Studies have demonstrated the significant roles of plant silicification in general (i.e., not only aquatic plants, but also terrestrial vegetation) on the silica cycle.It is well-established that 55-113 Tmol Si yr −1 is fixed globally in the terrestrial biosphere (Avg.84 ± 29 Tmol Si yr −1 ; Carey and Fulweiler, 2012).This annual ecosystem biogeochemical uptake exceeds the annual export from continents to the ocean through physicochemical weathering of silicates by two orders of magnitude (Conley, 2002), and thus exerts a strong control on global land-ocean interactions in the silica cycle.Specifically for freshwater aquatic ecosystems, studies have pointed to a significant effect of phytoliths on the silica cycle, and strong potential storage of phytolith-silica in soils and sediments (Cary et al., 2005;Struyf et al., 2007).The carbon occluded in the phytoliths of wetland plant species potentially even forms an important sink to consider in the carbon cycle (Li et al., 2013).Exemplary to this is the strong (permanent) silica sink in the Okavango Delta sediments (Botswana) that can be associated with a dominance of silica-rich tropical giant-grasses such as Cyperus papyrus L. (Struyf et al., 2015;Mosimane et al., 2017;Schoelynck et al., 2017).BSi accumulation in macrophytes, and thus potentially in sediments, varies strongly between species.Aquatic vegetation shows apparent plasticity regarding silica uptake, adaptive to silica availability, water and wind dynamics, light interception, herbivory and nutrient stress (Schoelynck and Struyf, 2016).High silica uptake results in high BSi concentrations in plant litter (Struyf et al., 2005;Carey and Fulweiler, 2014), which can impact on aquatic decomposition processes (Schaller and Struyf, 2013) and thus delays in and/or reductions of the transfer of dissolved silica (DSi) from land to ocean.Besides macrophytes, diatoms (Bacillariophyta, unicellular algae), testate amoebae (polyphyletic assemblage of protozoa with a shell) and sponges (Porifera) also contribute significantly to the BSi storage in wetland ecosystems [see Puppe et al. (2015) for a classification of BSi pools in terrestrial ecosystems].All diatom species and a large number of the testate amoebae synthesize amorphous siliceous shells (Aoki et al., 2007;Smol and Stoermer, 2010) or have an endoskeleton stiffened by (sponge) spicules (Maldonado et al., 2010).
It has been shown that wetlands, as is the case for many other nutrients, play an important role in the biogeochemical silica cycle (Struyf and Conley, 2009).In temperate and subarctic climates, both riparian (e.g., Struyf et al., 2009Struyf et al., , 2010) ) and tidal wetlands (e.g., Struyf et al., 2006;Müller et al., 2013) play an essential role in Si cycling, but tropical wetlands have not yet received similar attention.Tropical rivers deliver about 70-80% of the global DSi load into the ocean (Beusen et al., 2009), implying it is crucial to assess environmental factors that can influence its transport.Studies have indicated the potential influence of Si uptake by giant grasses and sedges (e.g., Ding et al., 2008;Cardinal et al., 2010) on Si-isotope biogeochemistry in tropical rivers, but did not consider the importance of the wetlands in large scale Si balances.
In this study, we hypothesized that macrophyte phytoliths play an important role in silica storage in three tropical lagoons of Ivory Coast.The observed plant species are typical for eutrophic waters of numerous tropical rivers, streams, reservoirs, and natural lakes.After vegetation senescence, plant BSi (phytoliths) may contribute to the BSi stock in lagoon sediments.We quantified this contribution by counting the relative abundance of plant phytoliths, diatom frustules and sponge spicules using a light microscope.After dissolution, DSi may again be taken up by vegetation or be transported downstream.Therefore, we made conservative estimates of the storage of BSi in the vegetation and sediments and compare it between the different lagoons.

Description of the Study Area
Kodjoboue, Ono and Hebe are three small lagoons located in the South-East of Ivory Coast with a joined surface area of about 1,150 ha (Figure 1).These lagoons are mainly fed by ground water seepage and flow into the Comoé River before reaching the Atlantic Ocean through different small channels.The catchment area of Ono lagoon has been dominated by pineapple cultivation since 1960 while the watersheds of Hebe and Kodjoboue lagoons are influenced by palm and coconut crops.Today, the three lagoons behave as freshwater systems due to the absence of seawater input since the opening of the Vridi channel in the near of Abidjan in 1950, that discharges most of the flow.The lagoons became gradually dominated by freshwater macrophytes following the reduced salt water inputs.The development of plants is particularly high in Ono lagoon where it starts hampering water transport and fishery activities.Between 1986 and 1989, only Salvinia molesta D.S.Mitch.and Eichhornia crassipes (Mart.)Solms were reported occasionally along the shores of Ono lagoon (Amon et al., 1991;Guiral and N'da, 1991).Today, Ono lagoon is 66% covered by different macrophyte species while the two other lagoons remain more or less in their original state with (emergent) plants limited to the littoral zones.The climate of the area is sub-equatorial with bimodal character.Two rainy seasons, a long rainy season (April-July) and a short rainy season (October to November), are intercepted by two dry periods: a long dry season (December-March) and a short dry season (August-September).The average annual precipitation is 1,704 mm (referred to the years 1970-2014); average annual temperature is 26.3

Sampling and Analytical Techniques
The following physicochemical parameters were measured in situ with a probe (Ap-5000, AQUAREAD Limited, Broadstairs, UK): depth (m), pH, temperature ( • C), and electric conductivity (µS cm −1 ).Subsurface waters (depth ∼30 cm) were sampled with a 1.7 L Niskin bottle.Water samples for DSi analyses were filtered through Sartorius cellulose acetate filters and refiltered through 0.2 µm pore size polysulfone filters.DSi concentrations (mg L −1 ) were determined with a spectrophotometer (DR 6000) according to standard techniques (Grasshoff et al., 1983).DSi samples were collected only in Ono and Kodjoboue lagoons during eight sampling campaigns at three different stations in each lagoon (Figure 1).For Ono lagoon, station 1 was close to dense Echinochloa pyramidalis (Lam.)Hitchc.and Chase vegetation, station 2 was close to Hydrilla verticillata (L.f.) Royle vegetation, and station 3 was devoid of macrophytes and located in the main channel connecting Ono lagoon to the Comoé River.In Kodjoboue lagoon, station 1 was located in the center of the lagoon, station 2 was further to the west and close to the main channel, and station 3 was to the east.All three stations were devoid of macrophytes because plants are limited to the littoral zones.
In parallel, sediment samples were randomly collected at three occasions (April, July, and September) in Ono and Kodjoboue lagoons (same locations as for water samples) and one time (April) in Hebe lagoon (sampling station devoid of macrophytes because plants are limited to the littoral zones).Sediment samples were taken in the organic rich top layer of the sediment using a sediment corer of 3.3 cm diameter and 15 cm length.Sediment cores were subsampled at three different depths (0-5, 5-10, and 10-15 cm).Sediment samples were dried at 75 • C for 5 days and manually homogenized afterwards.Sediment subsamples of 30 mg were mixed with 25 mL of NaOH solution (0.5 M) and incubated in a water bath maintained at 85 • C for 5 h.Subsamples taken at 3, 4, and 5 h were filtered through 0.45 µm nitrocellulose Chromafil syringe filters (A-45/25) and analyzed for silica using the spectrophotometric molybdate-blue method (Grasshoff et al., 1983) on a colorimetric segmented flow analyser (SAN++, Skalar, Breda, The Netherlands).Sediment BSi concentration is derived from the intercept of the linear regression between sampling time and measured concentrations, as adapted from DeMaster (1981).
Macrophytes were collected randomly on a monthly base from October 2015 to September 2016 in the three lagoons (12 sampling campaigns).Free floating macrophytes were collected with a net by sweeping a demarcated quadrant (surface area: 0.315 m 2 ).The other plants were cut manually from a fixed area of 1 m 2 .After collection, macrophytes were thoroughly rinsed to remove sediments, algae and macroinvertebrates, sorted per species and total fresh mass was weighted.A subsample of ±100 g fresh mass of each plant species was oven dried at 75 • C for 5 days and weighed again.This conversion factor between fresh and dry mass of the subsample enabled to calculate total macrophyte dry mass per unit surface area.The dried samples were then ground to 300 µm.BSi was extracted from 25-30 mg of dried plant material by incubation in a 0.5 M NaOH solution at 80 • C for 5 h (DeMaster, 1981), and analyzed for silica using the spectrophotometric molybdate-blue method (Grasshoff et al., 1983) on a colorimetric segmented flow analyser (SAN++, Skalar, Breda, The Netherlands).

Diatom, Sponge Spicule, and Phytolith Fixation Protocol
The presence/absence of BSi particles (primarily diatom frustules, sponge spicules and phytoliths) was determined in 12 sediment samples and 6 (non-rinsed) dominant plant samples (see Table 1 for species list) across all three lagoons (0-15 cm depth).Samples were prepared for light microscopy following the method described in Van der Werf (1955).Small parts of the sample were cleaned by adding 37% H 2 O 2 and heating to 80 • C for about 1 h.The reaction was completed by addition of saturated KMnO 4 .Following digestion and centrifugation (three times 10 min at 3,700 rpm), cleaned material was diluted with distilled water to avoid excessive concentrations of BSi particles on the object slides.Cleaned material was mounted in Naphrax.The object slides were analyzed using an Olympus BX53 microscope, equipped with Differential Interference Contrast (Nomarski) optics and the Olympus UC30 Imaging System.The variation in BSi particles in both sediment and plant samples was visualized and BSi particles were quantified (up to 100 particles) in sediment samples on random transects at 400x magnification and the relative abundance of each group (diatoms, sponge spicules, and phytoliths) was determined.Diatom species were identified up to genus level using Lange-Bertalot et al. ( 2017) and phytolith morphotypes were described according to Madella et al. (2005).

Stock Calculations
To calculate the BSi stocks in macrophyte biomass, we selected the dominant macrophyte species from each lagoon (see Table 1 for species list) and multiplied the average calculated dry biomass of each plant species (g m −2 ) with their respective average BSi concentration (mg g −1 ).These values were converted to t of BSi by multiplying with the respective vegetation coverage to obtain the total stock in vegetation in each lagoon.For Ono lagoon, the area covered by macrophytes was obtained by a satellite image of 2016 (Landsat 8-OLI/TIRS) using ArcGIS.For the two other lagoons, we did not get the images but field observations (Dr.Egnankou from the SOS Forêts, pers.obs.) show that the area covered by macrophytes is small and represents only about 1 and 5% for Hebe and Kodjoboue lagoon, respectively.To calculate the BSi stocks in the sediment, BSi concentrations in sediments were converted to t of BSi by multiplying the average concentrations (mg g −1 ) by the density of the sediment (g m −3 ), a depth of 10 cm (ca.rooting depth of macrophytes) and by the total surface area of the respective lagoon.

Statistical Analyses
By means of an Anova One Way test, average values of macrophyte BSi were statistically compared between species within the same lagoon, and average values of water DSi were statistically compared between months within the same lagoon [using Prism 5.00 (GraphPad)].P-values are not explicitly mentioned hereafter but "significant(ly)" refers to p < 0.05 and "not significant(ly)" refers to p ≥ 0.05.Using linear regression, we did not find relationships between plant BSi, sediment BSi and water DSi, hence these results were excluded from the text.

Water Quality Parameters and DSi Concentrations
In general, pH was low with average values below 7 (Table 2).
Conductivity was also low with average values ranging from 13.5 to 36.7 µS cm −1 showing the highest values in Hebe lagoon and the lowest in Kodjoboue lagoon.Both parameters are characteristic for freshwater systems and demonstrate the present-day absence of marine influence.Temperature varied slightly between the three lagoons with average values ranging from 27.5 to 29.9 • C. DSi showed a strong temporal variability with a similar pattern in Ono and Kodjoboue lagoons (Figure 2).The highest values were observed between March and April and the lowest in September.When comparing the two lagoons, the highest DSi values were generally observed in Ono lagoon.Overall, DSi concentrations ranged from 1.7 to 10.8 mg L −1 .

Biogenic Silica Accumulation in Macrophytes
Eight plant species were found in Ono lagoon, seven in Kodjoboue lagoon and only three in Hebe lagoon (Table 1).The most dominant species were Echinocloa pyramidalis followed by Hydrilla verticillata and Eichhornia crassipes for Ono lagoon, E. pyramidalis, Panicum parvifolium, and Acroceras zizanoides for Kodjoboue lagoon, and Salvinia molesta, E. pyramidalis, and Rhynchospora corymbosa for Hebe lagoon.All other species were less abundant.BSi concentrations ranged from 0.8 to 54.9 mg g −1 with the lowest concentration found in E. crassipes in Ono lagoon and the highest found in S. molesta in Kodjoboue lagoon (Figure 3).Generally, free floating species had low BSi concentrations and much lower than submerged and emergent species.Among the emergent species, Brachiaria villosa and A. zizanioides showed the highest values while the lowest were found in Jussiaea repens.BSi concentrations of B. villosa were significantly higher than those of the other plants collected in Ono lagoon.In contrast, the BSi concentrations of A. zizanioides in Kodjoboue lagoon did not differ significantly than those of S. molesta (a freefloating species), P. parvifolium and B. villosa.For the single submerged macrophyte (H.verticillata) found in Ono lagoon, the concentrations were relatively high (5.1-19.8mg g −1 ) and similar to those found in emergent species.In Hebe lagoon, BSi concentrations found in R. corymbosa were significantly higher than those obtained in S. molesta and E. pyramidalis.
Ono lagoon stored less BSi per m² macrophyte biomass than the two other lagoons.Average values for the three lagoons were 0.054, 0.102, and 0.288 t BSi ha −1 for Ono, Hebe and Kodjoboue, respectively.By extrapolating these values to the total surface covered by macrophytes in each lagoon, the BSi stored in macrophyte biomass was 0.2, 6.1, and 17.4 t for Hebe, Kodjoboue and Ono, respectively.

Biogenic Silica Accumulation in Sediment
The concentrations of BSi generally decreased with depth (Table 3): the highest values were found near the subsurface (0-5 cm) while the lowest occurred at the bottom of the sediment core.The lowest values were observed in July in Ono lagoon.Overall, average BSi concentrations in sediments ranged from 49.6 to 89.4 mg g −1 .The BSi stock in the sediments over 10 cm depth was higher in Kodjoboue lagoon (8.2 t BSi ha −1 ) in comparison to those observed in Hebe (4.2 t BSi ha −1 ) and Ono (5.5 t BSi ha −1 ).Extrapolating these values to the total surface of each lagoon, the BSi stock in the sediments was 1,025, 2,651, and 3,477 t for Hebe, Ono and Kodjoboue, respectively.

Characterization of the BSi Particles in Sediment and Plant Samples
All sediment samples were dominated by diatom frustules, mainly belonging to the genera Aulacoseira (>90%), Pinnularia (2-3%), Eunotia (2-3%), and Diadesmis (<1%).In the samples from Kodjoboue lagoon, several sponge spicules were observed (Figure 4a, arrow).Only megascleres were found.Gemmoscleres and microscleres were never observed.In Ono and Hebe lagoons, the amount of spicules was extremely low and spicules were  usually not included in the 100 particle-counts.Phytoliths were only very rarely observed (Figure 4b).In the samples from Ono lagoon, no phytoliths were observed, even when scanning an entire object slide.In Hebe lagoon and Kodjoboue lagoon, occasionally, only fragments of phytoliths were recorded after scanning entire object slides but never during counting.The observed phytolith fragments were not identifiable.Intact shells of testate amoebae were not found, but very rarely idiosomes (siliceous platelets synthesized by testate amoebae for shell construction) were observed (Figures 4c,d).

The Role of Macrophytes and Diatoms in Lagoon BSi Storage
BSi accumulation in vegetation varied per lagoon and was dependent on plant species dominance [see Prychid et al. (2003) and Hodson et al. (2005) for BSi contents in plants related to their phylogenetic position].The BSi stock in macrophytes ranged from 0.054 to 0.288 t BSi ha −1 , which is situated between values reported for submerged macrophytes in the Okavango Delta channels (0.008 t BSi ha −1 ; Schoelynck et al., 2017) and values reported for emergent wetland vegetation in the Okavango Delta (ranging from 0.010 to 1.600 t BSi ha −1 ; Struyf et al., 2015).The highest BSi concentrations were found in macrophytes collected from Kodjoboue lagoon and were related to the dominance of emergent species (A.zizanioides, B. villosa, R. corymbosa, E. pyramidalis, P. lanceolatus, and P. parvifolium) belonging to the Poaceae and Cyperaceae families.Species of these two families are often characterized by high BSi contents (Ma and Takahashi, 2002) and were also identified as a macrophyte growth form with high BSi content in Schoelynck and Struyf (2016).The generally lower BSi concentrations found in Ono lagoon were counterbalanced by the high macrophyte coverage giving this lagoon the highest BSi stock in vegetation of all three lagoons: 17.4 t.We did not find a relationship between plant BSi and sediment BSi.Other studies (e.g. in the Okavango Delta, Struyf et al., 2015) identified an impact of vegetation on phytolith input to the soils and thus on BSi content, but results are not always consistent.In saltmarshes, de Bakker et al. (1999) found no relation between plant BSi content and sediment BSi content.Carey and Fulweiler (2014) concluded that silica uptake in saltmarsh plants is not always directly linked to sediment BSi and sediment DSi, since active uptake mechanisms stimulated by multiple abiotic factors (e.g., hydrological stress) can strongly impact plant silica uptake.
Light microscopy of the lagoon sediments revealed that the majority of BSi particles in the sediment of all lagoons comprise almost entirely of diatom frustules, occasionally sponge spicules and almost never phytolith particles.Tychoplanktonic diatoms such as Aulacoseira species form long filaments of large, heavily silicified, cylindrical cells, and are known to bloom mainly in spring and autumn (Siver and Kling, 1997).The light microscopy investigated samples in the current study were taken in November 2017 and April 2018, most likely in the middle of large Aulacoseira blooms.During this sampling period, stimulated by the nutrient rich character of these waters, the shallow depths and the abundant sunlight, Aulacoseira blooms may have produced large amounts of frustules (Conley et al., 1993).Additionally, non-rinsed plant samples were also dominated by large amounts of epiphytic diatoms, mainly belonging to the genus Eunotia, typically found in acidic environments such as these lagoons (Van de Vijver, pers.obs.).It is thus clear that diatoms could have depleted the plantavailable DSi in the water rapidly, before the plants could take up the necessary silica to produce phytoliths.This can be observed in Figure 2 where the amount of DSi strongly decreased between April and May suggesting either a dilution effect in the rainy season (resulting in lower concentrations; Meunier et al., 2011), or an increased uptake by organisms compared to other months (Rabosky and Sorhannus, 2009).
The fact that we do find BSi in plants but rarely phytoliths, may be explained by different deposition mechanisms.Rapid deposition of silica results in lumen and cell wall phytoliths, whereas slow deposition results in silica deposition in intercellular spaces or in an extracellular (cuticular) layer (reviewed by Prychid et al., 2003;Hodson, 2019).These nonhomogeneous masses can be i.a. a silica ring formed around the periphery of the cell, deposited on a dispersed organic matrix, or laid down as deposits within the cell wall or between the cellulose wall and the plasma membrane or in cortical intercellular spaces within the cell (Prychid et al., 2003).In general, these fragile phytogenic silica structures as well as very small phytoliths (<5 µm) are difficult to trace in sediments, but they potentially play an important role in silica cycling (Meunier et al., 2017;Puppe et al., 2017).Additionally, given the very large number of epiphytic diatoms on non-rinsed plant samples, it cannot be excluded that a fraction of plant-silica in our results is actually diatom-silica, even though plants were washed thoroughly before analysis, according to protocol.Furthermore, it could be hypothesized that plant organic matter may float away through the outlets of the lagoons, or accumulate on the shores, hence not contributing to the sediment silica storage (cf.Struyf et al., 2015).Alternatively, Puppe et al. (2017) suggest the formation of a layer of coarse organic matter on top of sediments from which phytoliths cannot be released easily and hence are missing in underneath sediment layers.We have no arguments supporting either of these hypotheses, none of our sediment cores contained such a distinct organic top layer and samples were visually homogeneous organic rich (fine material) over the entire depth of the core (15 cm).

Sediments as a BSi Sink
Overall, the BSi stocks in sediments of the three lagoons ranged from 4.2 × 10 3 kg ha −1 to 8.2 × 10 3 kg ha −1 and were much higher than those found in the vegetation, suggesting that lagoon sediments are important sites for long-term BSi storage.Similar studies on African wetland BSi accumulation are scarce.The best studied ecosystem in this regard is the Okavango Delta (Botswana) where similar values were found: between 3.8 × 10 3 kg ha −1 and 36 × 10 3 kg ha −1 in the riparian wetlands [depending on flooding frequency and in top 5 cm of sediment; Struyf et al. (2015)] and between 0.2 × 10 3 kg ha −1 and 1.5 × 10 3 kg ha −1 in the main channels [depending on vegetation cover and in top 10 cm of sediment; Schoelynck et al. (2017)].Since phytoplankton production in the Delta is limited by low ambient nutrient concentrations, Okavango sediment BSi particles mainly originate from dead terrestrial and aquatic plant fragments (Struyf et al., 2015).In case of the analyzed groundwater-fed Ivory Coast lagoons, most of the BSi is produced within the lagoons and not brought in by any river.Overland runoff and topsoil erosion as in the Nyong basin in Cameroon (Cary et al., 2005) can potentially introduce allochthonous BSi particles, especially by transformation of an area from pristine to agricultural land use (cf.Smis et al., 2011).However, we have no indication at the moment that this is a significant process in the analyzed lagoons and no phytoliths with confirmed terrestrial origin were found in corresponding sediment samples.The investigated lagoons sediment BSi pool is thus dominated by diatom frustules, occasional sponge spicules and phytolith fragments, and an unquantified amount of fragile phytogenic silica structures and small phytoliths (<5 µm).Entire shells of testate amoebae were not observed probably because the organic cement, which glues the idiosomes (i.e., the building blocks of the shell) together, was destroyed during the H 2 O 2 sample preparation step.But also single idiosomes were very rare in our samples, although other studies showed that they can be quite abundant in sediment samples in general (cf.Douglas and Smol, 1987;Cary et al., 2005).

Implications for the Ecosystem Silica Cycling
The in situ produced BSi is deposited locally due to limited exchange with the adjacent coastal ocean (Boynton et al., 1996).This accumulation strongly influences BSi cycling and storage in the sediments.In general, plant and phytoplankton BSi is returned in detritus to soil and sediments, where decomposition processes and pedogenic transformations can result in either stored BSi or re-dissolved DSi (Frings et al., 2014).As the analyzed lagoons have an acidic pH (<6.4), the preservation of diatom frustules and phytoliths can be assumed to be quite good (Flower and Ryves, 2009) making silica only very slowly available again for uptake.The average DSi concentration for the Comoé river draining these lagoons is with 2.86 mg L −1 rather low (Koné et al., 2008).Invertebrate fauna might be important contributors to remineralisation processes in general as was reported for other nutrients in Taabo Lake, which is similar to the studied lagoons (Kouamé et al., 2011).Invertebrate fauna, especially burrowers, are known to amplify BSi cycling in shallow environments like coastal lagoons and estuaries (Viaroli et al., 2013).The presumed fragile phytogenic silica in the plants (potentially explaining the relatively high plant BSi concentrations) may be very important in such conditions as they may be able to dissolve more easily than homogeneous masses (i.e., phytoliths and diatom frustules) and can hence contribute to local silica cycling.Plants may then be an important source of silica availability in the lagoons, although this warrants further investigation.
Generally, BSi stored in wetland sediments is several orders of magnitude more soluble than mineral silicates (Farmer et al., 2005) and its temporal or permanent storage vs. recycling or downstream transport exhibits an important control on silica export toward the ocean (e.g., Derry et al., 2005;Fulweiler and Nixon, 2005;Sommer et al., 2006), which is crucial for the silica cycling in these (sub) tropical regions (Beusen et al., 2009).Decreasing availability of silica compared to other nutrients can negate any competitive advantage the diatoms have and can lead to nuisance and toxic blooms of green and blue green algae in coastal ecosystems (Cloern, 2001).Diatoms need an optimal nutrient ratio of C:Si:N:P = 106:15:16:1, and diatom growth will cease when DSi supplies are depleted, allowing other phytoplankton classes to proliferate using any excess of N and P (Anderson and Burkholder, 2002).Because diatoms are primary producers at the base of the food-web, any change in that production might have knock-on effects to the entire ecosystem.Phytoplankton assemblages at Ivory Coast coastal environment show that, although diatoms are the most species-rich, cyanobacteria are the most abundant (up to 90% of the phytoplankton biomass), and silica was identified as predominant abiotic factor controlling phytoplankton dynamics (Osemwegie et al., 2016).Expected human population rise in the greater Abidjan area however, will probably put more pressure on local water bodies and nutrient levels may increase to levels which could escalate to dramatic levels of eutrophication [projections made for the nearby (and larger) Ebrié Lagoon which receives water from the Comoé River; (Scheeren et al., 2004)].

CONCLUSION
It can be concluded that relatively recent colonization of the lagoons by aquatic vegetation had no major effect yet on the lagoons' sediment BSi storage that is still dominated by phytoplankton, entirely comparable to any other open water system lacking plants.Perhaps it is just a matter of time before plants start having an impact here.Values reach up to 8.2 t BSi ha −1 , which is similar to other tropical wetlands and may exhibit an important control on fluxes toward the tropical coastal zone.Relatively high plant BSi concentrations are likely explained by non-homogeneous fragile phytogenic silica structures that may be more able to dissolve and contribute to local silica cycling, rather than the more homogeneous diatom frustules.Using other techniques are advised measuring or visualizing these fragile phytogenic silica structures to prove this hypothesis.

AUTHOR CONTRIBUTIONS
YK and JS conceived the program research.YK collected the samples.BV analyzed the BSi particles in sediment and plant samples.YK, JS, and BV wrote the paper.JS and BV read and approved the final manuscript and we are pretty confident that YK would have agreed with our corrections after revision.

FIGURE 1 |
FIGURE 1 | Map showing the location of the different lagoons with a photo impression of their macrophyte cover (A, Ono; B, Kodjoboue; C, Hebe).

FIGURE 2 |
FIGURE 2 | DSi concentrations (mg L −1 ) in the water column of Ono and Kodjoboue lagoons.Data represent averages of different sample points (n = 3) and error bars indicate standard deviation.Different letters indicate a significant difference (p < 0.05).

FIGURE 3 |
FIGURE 3 | Box-plots showing average BSi concentrations in the different macrophytes for the three lagoons.The vertical bars refer to the minimum/maximum values and the different letters indicate a significant difference (p < 0.05).

FIGURE 4 |
FIGURE 4 | Silica particles in the sediment and plant samples.Sediment samples originate from Kodjoboue lagoon, and are exemplary to all examined samples from the three lagoons.(a) Large number of diatom valves, mainly belonging to the genus Aulacoseira and one sponge spicule (megasclere, see arrow).(b) Shows a few broken diatom valves and one phytolith (see arrow).Plant samples originate from dominant plant species from across all three lagoons.(c,d) Show testate amoeba idiosomes.(e-h) Show different phytoliths found in Echinochloa pyramidalis samples: short cell bilobate (e,f), cylindrical polylobate (g), and globular echinate (h) phytoliths were observed.(i-k) Show different phytoliths found in a Panicum parviflorum sample: short cell cross (i), short cell bilobate (j), and long cell echinate (k) phytoliths were observed.Scale bar represents 20 µm for panels (a,b) and 10 µm for panels (c-k).

TABLE 2 |
Physicochemical characteristics (average ± SD) of the three lagoons.

TABLE 1 |
Characteristics and occurrence of macrophyte species collected in the three lagoons (×, presence of the species; #, species used for phytolith visualization; *, species' values used in BSi stock calculation).

TABLE 3 |
BSi concentrations (mg g −1 ) in sediments of the three lagoons (average ± SD; nd, no data).