Spatial and temporal variation of the macrophyte assemblage in Santo Tomás, a wetland in the Caribbean Colombian

: Aim: The associated flood events to floodplain of many rivers in the world affect the composition and structure of aquatic biota due the water levels variation can ensure ecological integrity of associated wetlands. This study describes the spatial and temporal variation in the macrophyte assemblage in the Santo Tomás wetland during seasonally flood-pulsed of the Magdalena River (North of Colombia). Methods: For eight months, between 2017 and 2018 samplings were carried out in three stations for composition and abundance of aquatic macrophytes. Spatio-temporal pattern of richness was estimated with Chao methodology and vegetation cover was calculated using Kruskal-Wallis and U Mann-Whitney tests. A range abundance curve was used for species dominance between samplings and stations. A non-metric multidimensional scaling (nMDS) was used to analyze the Spatio-temporal distribution, Canonical Correspondence Analysis (ACC) was performed to relate the physicochemical variables to the species composition. Results: A total of 24 species of aquatic plants distributed in 23 genera and 15 families were registered in the transects. The most abundant species (% of vegetation cover) were Ipomoea aquatica (19%), followed by Ludwigia helminthorrhiza (14%) and Eichhornia azurea (13%). The most frequent life form was free-floating, followed by the emergent one. The Spatio-temporal changes and the highest values of richness and vegetation cover in the Santo Tomás wetlands occurred during the filling and high waters period. L. helminthorrhiza , E. azurea , Pistia stratiotes , Neptunia oleracea , I. aquatica , Salvinia auriculata, and Hymenachne amplexicaulis were the most dominant species. The quality water of Santo Tomás Wetlands showed spatial and temporal variations during flooding pulse and some physicochemical variables such as organic matter (COD, BOD 5 ), pH, depth, ammonia, and fecal coliforms were related to macrophyte community composition Conclusions: Spatial and temporal changes of aquatic plants in Santo Tomás wetlands were related to the flooding pulse of Magdalena River.

The environmental conditions in aquatic ecosystems regulate the life forms of macrophytes since emergent, rooted submerged, and freefloating use light, nutrients, and other resources in a specific way and also differed in their response to environmental changes (Schneider et al., 2018).According to Regmi et al. (2021) in sub-tropical floodplain wetlands of Nepal, excessive growth of floating macrophytes in particular invasive species such as Eicchornia crassipes and Typha angustifolia tend to increase biomass during high water levels periods, while with water level drawdown, species emergents like Alternantheria spp.and Cyperus spp increase and establishing.Additionally, several authors have been concluded that submerged aquatic vegetation (Example: Ceratophyllum sp., Egeria spp., Lemna spp, Marsilea spp) is an environment indicator it exists where there is a better water quality condition (Rameshkumar et al., 2019) Tropical floodplains can show a progression from grasslands to emergent aquatic plants
depending on the depth and frequency of flooding (Murray-Hudson et al., 2014).Thus, floods can act as a homogenizing factor that promotes changes n of aquatic macrophyte species, whether they increase or decrease during flood cycles (Catian et al., 2018;Delatorre et al., 2019).Consequently, this study aimed to examine macrophytes composition and their Spatial-temporal dynamic during different climatic periods in relationship with the water physicochemical variables in Santo Tomás wetland on flood plain Magdalena River (Colombian Caribbean).

Study area
Santo Tomás wetland (STW) is a lentic ecosystem located in the floodplain area of the Magdalena River (Arias, 1985).STW is in the Atlantic department between 10° 45' N and 74° 44' O in Santo Tomas municipality in the Caribbean Colombian.It is connected directly to the Magdalena River by a 2.12 km long canal entitled "El Canal" (The Channel) (Figure 1).
The total area of STW is 650 ha at 8 meters above sea level.Temperature ranges from 26 to 30°C with a relative humidity between 84 and 95% and a total annual precipitation of 1200 mm (Escolar, 2007).There are two dominant tropical seasons in the region.In the dry season (December-April), the surface water is significantly reduced (approximately 67 ha), with maximum depths of 1 m, while in the rainy season (May-October), the surface water can cover the total wetland area, 650 ha surface and reach 3 m in depth.However, El Niño Southern Oscillation (ENSO) can modify drastically the season in the Caribbean region (USAID, 2016).In fact, in 2010 the STW suffered severe flooding (Alvarado, 2016), and in 2016 caused a drought reduced surface water area by 90% (direct observations).Additionally, STW is affected by discharges of treated wastewater from lagoon stabilization ponds.

Field sampling design
Field samplings were carried out, one by month, for eight months in the STW during the rainy and dry seasons.We sampled abiotic and biotic data from April to September 2017 (rainy season) and from February to April 2018 (dry season).Consequently, two hydrology periods were studied: rising water level (high water level) and decreasing water level (low water level) (Figure 2).
We established three permanent stations concerning impact degree: S1, The Canal's direct communication with the Magdalena River, S2 the area near the domestic wastewater discharge, and S3, Santo Tomás municipality population.
Three permanent stations registered with a GPS (Garmin Etrex 20) were established for monitoring.Site S1 was El Canal which communicates with the Magdalena River, site S2 was located in the area near the discharge of domestic wastewater of Santo Tomas municipality of Santo Tomás, and site S3 was located near the urban area.At each site, every 20 m, three transects of 10 m 2 each were established.These were located in parallel in the areas most exposed to aquatic vegetation.Acta Limnologica Brasiliensia, 2022, vol.34, e22

Macrophyte community
We sampled all macrophytes we found inside the quadrats.Aquatic macrophyte taxa were identified according to Velásquez (1994) The vegetation cover (%) was estimated in the climatic year according to the Magdalena River flood pulse, following the methodology by Ramos-Montaño et al. (2013).Through a small boat by direct observation and using a 1m2 quadrant which one was moving to cover the total transect area (10m 2 ).Three replicates at an interval of 20 m2 in each transect by sampling site were reached.Once the data was obtained, was calculated as a sum of relative frequency and relative abundance of macrophytes and 100% of each transect was weighted.Statistical and diversity analyses by site and period were performed with these data.
To facilitate the estimation of the vegetation cover in the transects, a general recognition of the aquatic plant species of the STW was carried out.For this purpose, a reference collection was made and deposited in the Armando Dugand Gnecco Herbarium of the Universidad del Atlántico.For their identification, they were compared with herbarium collections and with the help of specialized literature (Velásquez, 1994;Posada & López, 2011, Cirujano et al., 2014and Madriñán et al., 2017).According to Gyosheva et al. (2020), they were classified into three life forms: floaters (plants on the water surface), rooted floaters (plants rooted in sediment with leaves floating on the water surface), and emergent plants (the aerial part of the plant emerges above the waterline).

Environmental data
Water physicochemical parameters such as pH, ORP (EH), and dissolved oxygen (DO) were measured in situ using a multiparameter probe (Hanna HI 9829), both at the beginning and at the end of each transect to calculate the average of each parameter.To minimize fluctuation in water quality parameters due to temperature, measurements were carried out between 9-11 am.Transparency was measured with a Secchi disk.Water samples were collected, stored at 4 °C, and transported to Laboratories of Universidad del Atlántico for organic matter, nutrients, and coliforms analysis following APHA-AWWA-WEF (APHA, 2012) methods.

Data analysis
Spatial and temporal species composition was evaluated with Chao et al. (2014) methodology.Hill numbers were used to measure species effectiveness using the EstimateS 9.2 program (Colwell, 2009).True Diversity, species richness ( 0 D), Shannon Index exponential diversity ( 1 D), and Simpson Index inverse ( 2 D) were compared using the superposition of confidence intervals at 95% in R program version 3.1.Spatial and temporal variations were evaluated with the Kruskal-Wallis test, supported by the U Mann-Whitney paired comparisons test.To compare the dominance and uniformity of species between samplings and stations, rangeabundance curves were used (Cultid-Medina et al., 2012).
A community Principal Component Analysis (PCA) was used to establish the dominant macrophyte species.Last, Canonical Correspondence Analysis (CCA) was used to investigate the relationship of abundant species spatially and temporally.Nonmetric multidimensional scaling (nMDS) was performed to investigate the variations between periods and sampling months (Arellano et al., 2013).ANOSIM similarity analysis was used to detect differences in community composition.For the multivariate analysis, the dominant species were first determined using the Olmstead-Tukey test (Sokal & Rohlf, 1995), considering the frequency of occurrence and abundance of each taxon.In addition, before the construction of the matrices for the analyses, the physicochemical variables whose values between the samplings and study sites were below the detection limit of the respective methods and those that were collinear or correlated with each other were excluded, performing a collinearity and Spearman's correlation analysis.Canonical Correspondence Analysis (CCA) was then selected to establish the spatial and temporal relationship between the dominant species with the set of physicochemical variables.These analyzes were realized using the Past 2.17 program (Hammer et al., 2001).

Water quality of the Santo Tomás Wetland (STW)
The variation ranges of the physicochemical parameters between the hydrological periods and between the sampling stations indicated slight transitions in most of the variables (Table 1).Electrical conductivity and COD values increased during low water levels while nitrates and total coliforms decreased under the same water level conditions.The concentrations of organic matter (COD, BOD 5 ) and nutrients (Ammonium, Nitrate, and Orthophosphate) were relatively similar in rising and falling waters, as well as in each sampling station.High concentrations of nitrate were recorded concerning ammonium, probably due to the effect of runoff from nearby agricultural areas.The transparency of the water was around 30 cm at low and high water levels.The pH was kept in neutral conditions and dissolved oxygen presented relatively low values throughout the study.

Aquatic plant: composition and abundance
A total of 41 aquatic plant species were registered and distributed in 36 genera and 22 families (Figure 3 and Table 2).Transects registered 24 species being Ipomoea aquatica was the most abundant of the species with 19% of cover, followed by Ludwigia helminthorrhiza (14%) and Eichhornia azurea (14%).
Ipomoea aquatica, Ludwigia helminthorriza, and Eichhornia azurea showed the highest vegetation cover in both stations 1 and 2 (Figure 5).These species were also present at station 3 but to a lesser extent.

Aquatic macrophyte richness: spatial and temporal variation
The greatest species richness occurred in April-May 2017 (raising water), while the lowest values occurred in February-April 2018 (falling water).During the high-water months (Apr-May-Jul/17), species richness ( 0 D), common ( 1 D), and abundant species ( 2 D) were high with values of 24, 10, and 8 respectively, compared to the recorded in the low water months (Feb-Apr/18) with values of 15, 5 and 4 respectively (Figure 6).The S2 and S3 stations were statically different from S1 and showed values of 24, 10, and 8 richness, common and abundant species, respectively.S1 presented lower values of 21, 8, and 5 for species richness, common species, and abundant species, respectively.Similarly, aquatic plant coverage percentage (Figure 4) was higher in April and May 2017 (13 and 14%, respectively) and significantly different (P <0.05), in February and April 2018 (11 and 9%, respectively).Spatially, vegetation cover was higher and significantly different in S1 (38%) compared to S2, where the lowest value was recorded (31%).
I. aquatica, E. azurea, and L. helminthorrhiza, were more abundant in S1, S2, and S3 (Figure 5).P. stratiotes, N. oleracea, E. crassipes, S. auricular I. aquatica, L. helminthorrhiza, and E. azurea were observed in all months of sampling.During April, May, and Jun 2017 during the high water period, E. azurea, L. helminthorrhiza, and N. oleracea were the dominant species.In August and September 2017, when the water level was at their highest P. stratiotes and S. auriculata were abundant, and February-March 2018 when l water levels were at their lowest, I. aquatica and E. crassipes, were mainly dominants.
According to nMDS analysis, the species composition during the sampling period did not vary considerably and there were no statistical differences (P=0.181) between plant groups shown in the level water variations periods (Figure 7B).nMDS analysis shows a weak relationship between  species composition and sampling periods, contrary to ANOSIM similarity results (Table 4) showed that spatially a plant group of S1 station is separated from S2 and S3 stations with statistical differences (P<0.001).

Environmental variables and aquatic macrophyte composition
Differences between water variables, hydrological periods, station sampling, and its relation with macrophytes composition were registered in Canonical Correspondence Analysis (CCA).CCA results show 62.31% of the variation of the data (Figure 8) in the first two axes with statistical significance (p <0.02).I. aquatica and L. subintegra showed the highest vegetation cover in the S1 station (36.38% of the variance).Particularly, I. aquatica showed a significant increase in vegetation cover during the dry period (Aug/17 and Apr/18), and then appear spatial-temporally differentiated.N. oleracea, O. cubense, and L. helminthorrhiza were associated with the negative quadrant of the first axis.This was related to nitrate concentrations,    ), and ammonium (NH 4 ).M. polycarpa, H. amplexicaulis, and S. auriculata species showed high vegetation cover percentages in all sampling periods.M. polycarpa had a distribution restricted to the S1 station in samples Jun, Aug, Feb, and Apr.

Composition and abundance of aquatic macrophytes
The composition of aquatic macrophytes in Santo Tomás Wetland (STW) was related to the water level in the floodplain of the Magdalena River.During the falling and high-water period, the highest values of richness and total vegetation cover occurred.L. helminthorrhiza, E. azurea, P. stratiotes, N. oleracea, I. aquatica, S. auriculata, and H. amplexicaulis were the most dominant species.The highest richness and diversity were recorded in the first two months of sampling when rising water occurred.Similar behavior was registered by (Agostinho et al., 2004) in flood pulse Paraná river (Brazil).During rising water levels, depth influences richness since macrophytes is commonly distributed along with the water mirror and the littoral zones of lowland wetlands (Santos & Thomaz, 2007).
The increase in the water level creates biotic gradients (Vieira et al., 2019).Richness can be explained from the littoral zone until the habitat is off with greater depth.
More than 170 aquatic plant species have been reported in 34 shallow freshwater wetlands in the Colombian Caribbean region (Cortés-Castillo, 2017).Several studies affirm that the cosmopolitan nature of many species of aquatic plants confers a similar composition (Álvarez-León et al., 2004;Rangel-Ch, 2010;Cataño-Vergara et al., 2008).However, some species in this study were not included in those studies, indicating that diverse factors can affect species composition (Rial, 2006).E. azurea and E. crassipes species, typical in tropical wetlands, spread on the surface of the water and persist throughout the flood cycle (Cortés-Castillo & Rangel-Ch, 2013).Regularly, dominant at the end of the rainy period and the beginning of the dry period (Cortés-Castillo & Rangel-Ch, 2013).E. azurea formed scattered abundant patches in the water body and on the banks of the STW (S2 and S3 stations).This behavior has also been observed in L. helminthorrhiza, especially on the banks and in shallow water sectors.Both E. azurea and L. helminthoriza species present rapid growth, high reproductive potential, latency, and adaptability that facilitate high and rapid proliferation (Cortés-Castillo & Rangel-Ch, 2013).During the dry period, I. aquatica presented greater vegetation cover, mainly on the shores of the STW.In America, this species is considered an invasive aquatic weed with high reproductive efficiencies (more than 200 seeds per plant) and a high colonizer potential of calm and shallow water (Staples, 1996).
Shallower areas are dominated by emerging species rooted while deeper zone, favor the development of floating, epiphytic, or submerged aquatic plants.According to (Santos & Thomaz, 2007), the depth dynamic in flood pulse causes the coexistence of various forms of life due to replacing species and consequently, increasing richness.Similarly, the light incidence in the water column also influences the structure of macrophyte communities, especially for submerged species.Low radiation can seed germination restricted, and avoid the establishment of propagules (Kettenring et al., 2006;Roldán & Ramírez, 2008).The highest richness values recorded at the S2 station, likely coincided, with the highest depth and transparency recorded at that station.Wallsten & Forsgren (1989) and Gordon (1998) state that in shallow tropical wetlands the abundance increases during the flood period but richness has a relative decrease.In STW, this situation was presented during April-Jun/2017, in rising and high water, when the highest abundance and decreased richness were recorded.Rial (2006) reported that the diversity and abundance of aquatic plants change with water level variation.Periodic fluctuations in the water level allow the establishment of plant communities with changes in the dominance of certain species (Delatorre et al., 2019).During the dry period, when the water level decreases, some species establish taking advantage of the supply of nutrients from the decomposition of other plants (Pérez-Vásquez et al., 2015).Functional diversity also varies according to the flood cycle, due to macrophytes with life forms predominant defining the community characteristics (Catian et al., 2018).
Therefore, exclusively aquatic plants such as E. azurea, L. helminthorrhiza, N. oleracea, P. stratiotes, and S. auriculata, dominated in rising period water (April-May-Jun/17).Emergent, amphibians or rooted plants such as I. aquatica, E. crassipes, and H. amplexicaulis, were more abundant in the low water level (February -April/18).L. helminthorrhiza, N. oleracea, P. stratiotes, and S. auriculata, (freefloating) present stems or leaves with aerenchyma and rhizomes for propagation.At a high water level, these forms of life have a competitive advantage over other species and are environmentally selected (Catian et al., 2018;van der Valk, 1981).During the descent water period, the decomposition of the plants that were submerged in the flood phase occurs.This condition provides nutrients to the water, electrical conductivity increases, and floating and amphibian species optimize the nutrients available (Almeida & Melo, 2009).In this environment species such as I. aquatica and H. amplexicaulis, may have high growth due to not depending completely on the water to survive, making the best use of nutrients available, and the electrical conductivity of water conditions.
The aquatic plant's diversity responds to the evolutionary process and adaptations in landwater transition (Sculthorpe, 1967).Progressive adaptations to aquatic life suggested different reactions to fluctuations in the water level (Thomaz & Bini, 2003).Consequently, the heterogeneous environment imposes selection pressure on aquatic macrophytes that results in considerable Spatiotemporal diversity of life forms (Eckert et al., 2016).The clear temporal differentiation of macrophytes in the STW could be determined by the flood pulse and connectivity with the Magdalena River.The aquatic plant groupings correspond with water levels were contrasting.The cover, richness, and diversity changed throughout the study period as a response to the depth and extent of water body changes.This behavior has been reported in similar wetlands in the Colombian Caribbean (Cortés-Castillo & Rangel-Ch, 2013).
Extensive communities of L. helminthorrhiza, N. oleracea, and P. repens, increase the richness of species because, regularly during the rainy period, come accompanied by other species (Cortés-Castillo & Rangel-Ch, 2013;Gordon, 2000;Rangel-Ch, 2010).During the flood period, free-floating macrophytes proliferate in favorable conditions of nutrients, currents, and wind action (Vieira et al., 2019).Contrary, during the dry period, numerous species lose aerial tissues and die, decreasing the richness of species (Mereles et al., 2004).L. aequinoctalis and U. gibba (free-floating) and Aeschynomene cf.indica (aquatic-rooted), disappeared completely during the dry period.These variations in the composition and vegetation cover during the hydroclimatic regime have also been reported in Venezuela wetlands (Rial, 2006).
Regarding spatial similarity, the nMDS analysis showed that S1 stations were separate from S2 and S3 stations.Direct connection with the Magdalena River through the Channel caused a different dynamic.Probably, the current flow prevented the permanent establishment of some free-standing plants.However, during the water lowering period, the sediment accumulation favored slopes or small islands formation where species I. Spatial and temporal variation of the macrophyte assemblage... Brasiliensia, 2022, vol. 34, e22 aquatica was particularly abundant.This aquatic plant and M. polycarpa occupied a big extension from the litoral zone until the open water zone, similar to those reported in the Zaptosa wetland (Cesar, department) by Cortés-Castillo & Rangel-Ch, (2013).M. polycarpa is usually associated with shallow water and muddy substrates, while I. aquatica is related to shallow flooded sand soils (Cortés-Castillo & Rangel-Ch, 2013).S2 and S3 stations, areas far from the Magdalena River channel, allowed species establishments like N. arguta, a floating submerged macrophyte that grows in the most stable habitat and less current disturbance (Delatorre et al., 2019).In these stations, aquatic plant community composition was homogenous, likely the result of ecological functioning in the floodplain (Cortés-Castillo & Rangel-Ch, 2013).

Environmental variables vs composition and structure of aquatic plants
A direct relationship between organic matter and nutrient concentrations and macrophyte species development has been reported by several authors (Endut et al., 2016;Gezie et al., 2018;Ondiba et al., 2018;Rameshkumar et al., 2019).For example, E. crassipes have high phenotypic plasticity and tolerance to level water variations (Gezie et al., 2018).This plant can fast grow in low contaminate waters or hypertrophic wetlands, even in wastewater with high contents of heavy metals (Rameshkumar et al., 2019).In the rainy period, floating associations of L. helminthorrizha, E. crassipes, and N. oleracea, have been reported in Zapatosa wetlands and Ciénaga Grande de Santa Marta in the Colombian Caribbean (Cataño-Vergara et al., 2008;Cortés-Castillo & Rangel-Ch, 2013;Rivera-Díaz et al., 2013).In STW these associations were present in April-May/17 (rainy period) and coincided with high transparency and nitrate concentrations.Probably, the fertilization by flood pulse influence during the rainy period favors the fast growth and assemblage diversity of aquatic plants (Gómez-Ródriguez et al., 2017).Floating macrophytes such as N. oleracea, L. helminthorrhiza, E. crassipes, E. azurea, and P. stratiotes, intervene directly in water turbidity decrease and transparency increase (Martelo & Borrero, 2012;Meerhoff & Mazzeo, 2004).
On the other hand, aquatic plant velocity growth is affected by temperature.N. arguta was positively related to temperature during April and May/17, while, L. subintegra, and I. aquatica were negatively related to depth.Probably, low water level exposes areas susceptible to being rapidly colonized, by floating and emerging rooted species (Rangel-Ch., 2010) and according to the results obtained in the ACC, it could be indicated that M. polycarpa, H. amplexicaulis, and S. auriculata proved to be tolerant to low oxygen concentrations or can certainly grow in eutrophic wetlands.

Conclusions
Spatial and temporal variations in the macrophytes assemblage were related to the hydrological period.During rising and high-level water periods the richness, vegetation cover, and diversity were increased.In the dry period, richness and abundance were lower due to the water low level.Water quality and physicochemical variables such as depth, temperature, organic matter, and nutrients were related to the greater abundance of E. azurea, L. helminthorrhiza, N. oleracea, P. stratiotes, I. aquatica, L. subintegra, and E. crassipes species Some species of the aquatic plant reached highs cover percentage, such I. aquatica (19%) or E. azurea (14%) principally in S1 station (Channel with the Magdalena River) and S2 stations.The presence of opportunistic species such as I. aquatica and, E. azurea covered almost the entire area of these stations.

Figure 1 .
Figure 1.Santo Tomás municipality in Colombia and Santo Tomás wetland showing the sample sites.The black points indicate the transects number by the station.
; Posada & López, (2011), Cirujano et al. (2014) and Madriñán et al. (2017).When necessary, plants were compared in the herbarium reference: Armando Dugand Gnecco Herbarium from Universidad del Atlantico to confirm identification.We classified the macrophyte taxa into three life forms, based on Gyosheva et al. (2020) classification: free-floating (plants on the water surface), rooted-floating (rooted plants in sediment with floating leaves on the water surface), and emergent plants (the aboveground part of the plant emerges above the waterline).

Figure 3 .
Figure 3. Accumulation of aquatic plant curve using Caho 1 and Chao 2 methodology in Santo Tomás wetlands during the study.

Figure 4
Figure4shows the highest vegetation cover by month and their relation to the water level.In 2017 when the level of river water was around 60 cm, Eichhornia azurea, and Ludwigia helminthorriza were dominants, while in 2018, and low-level water, Ipomoea aquatica was dominant.High-level water favors the aquatic macrophytes' diversity.Ipomoea aquatica, Ludwigia helminthorriza, and Eichhornia azurea showed the highest vegetation cover in both stations 1 and 2 (Figure5).These species were also present at station 3 but to a lesser extent.

Figure 4 .
Figure 4. Vegetation cover during sampling months in the Santo Tomás wetland and water level variation.

Figure 5 .
Figure 5. Vegetation cover variation of principal aquatic macrophytes by station sampling in Santo Tomás wetlands during the study period.

Figure 6 .
Figure 6.Diversity (richness 0 D, common 1 D, and abundant species 2 D) analysis for station and month sampling in Santo Tomás wetlands.

Figure 7 .
Figure 7. Non-metric multidimensional scaling (nMDS) analysis of species composition in Santo Tomás wetland by stations (A) and months sampling (B).
coliform bacteria, and water transparency.This association was evident in the first three months of sampling (Aug.-Sep.-Oct.2017), corresponding to rising waters (Figure 2).A clear separation of the samplings carried out in S1 (negative quadrant of the second axis) is evidenced, mainly in the function of depth, (BOD, orthophosphates (PO 4 -3

Table 1 .
Variation ranges of physicochemical characteristics in high-and low-level water scenery and average values and standard deviation (in brackets) of the physicochemical characteristics of water in the sampling stations of the Santo Tomás wetlands.

Table 2 .
Aquatic macrophytes species registered in the Santo Tomás Wetlands, their life forms, and total vegetation cover during the study period.

Table 3 .
Richness and abundance variation in coverage percentage of the aquatic plants by stations (S) sampling in Santo Tomás wetlands.

Table 4 .
Station sampling comparisons with ANOSIM similarity test.