Using a Macroalgal Functional Form Approach to Assess the Level of Disturbance of Seagrass Meadows in Bahía of Nuevitas, Cuba (2000-2002)

A study on the spatial and seasonal variations of the associate macroalgae and epiphytes of Thalassia testudinum was carried out in Bahía de Nuevitas. Sixty-two species were identified: 3 Cyanophyta, 25 Chlorophyta, 8 Phaeophyceae, 23 Rhodophyta and 3 Magnoliophyta, with two new records for Cuba and 43 for the area. The differences in the specific composition of the macroalgae communities are determined by a space component related to the type of affectation in each area. The morpho-functional groups of macroalgae in the station with more nutrient influence were mainly foliaceous and filamentous. In the stations far from the city, the predominant morpho-types were the leathery and articulate calcareous indicators of lower nitrification levels. The abundance and diversity of macroalgae in the site affected by fisheries were lower due to the damage by bottom trawls. Seasonal variations were found in the relative abundance of the species, not in the diversity, which makes evident seasonal changes in the structure of the seaweeds, where some species replace others in the community.


Introduction
Seagrasses form highly productive ecosystems comparable to the most produc-tive biomes on Earth [1]. Seagrass ecosystems operate based on food chains underpinned by detritivores and herbivores, where seagrass, as well as their accompanying epiphytes and associated macroalgae, is consumed [2]. Seagrass systems have been recognized as an important resource with its functions listed by [2] [3] [4] [5]. Similar to other marine ecosystems, seagrass meadows have diminished their coverage in part driven by warming sea surface temperatures but also by the direct effect of human intervention [6]. The most obvious anthropogenic-driven impacts on seagrasses occur in coastal areas adjacent to human settlements, where physical disturbance, organic enrichment, and overfishing activities result in the greatest negative impacts [7].
Seagrasses constitute one of the largest biotopes in Cuba, with Thalassia testudinum dominating meadows (1 -4 m) on the island shelf. The total area of the Earth covered by seagrass beds is poor known, but recent estimates are between 300,000 and 600,000 km 2 [8] [9]. Thalassia is associated with a great diversity of organisms, among which macroalgae stand out. Macroalgae are more efficient than seagrasses in their ability to uptake nutrients, vulnerable to mechanical damage, and their abundance and species composition can be used to determine the degree of disturbance in the ecosystem [10] [11], refer to different morpho-functional groups in macroalgae and describe their relationship to the level of environmental disturbance. A sufficient understanding of these categories allows simplifying the interpretation of patterns in environmental conditions, whether biotic or abiotic, by including a high number of species [12]. In fact, macroalgae have been used as indicators of environmental quality, particularly when linked to pollution [13] [14].
Bahía of Nuevitas, Cuba, in one of the sites on the island where seagrass meadws are disturbed directly by wastewaters and dragnet fishing, increasing stress on marine biodiversity of the area [7]. To date, there are no records of any studies using the morpho-functional approach described by [10] as a tool for evaluating the degree of disturbance. There are several records of a survey conducted by [15] [16] [17] [18], identifying macroalgae and phanerogams, as well as their biomass. As such, in this study we report on the analysis of spatial and temporal variations of the specific composition and abundance of the accompanying macroalgae and epiphytes of Thalassia testudinum forming beds in three areas of Bahía of Nuevitas. These sites differ in the degree of anthropogenic influence, which will help to evaluate the sensitivity of using the morpho-functional approach. It will also contribute to expanding the knowledge of the marine flora in the Caribbean.

Materials and Methods
Assessments were conducted in three different stations within Thalassia testudinum beds growing in Nuevitas Bay. Station 1 was located in Cayo Puto (21˚31'06''N; 077˚16'15''W). This station is continuously affected by wastewater discharges from an industrial and urban origin that mix with its muddy substrate full of debris derived from red mangrove trees. Station 2 in Playa Bagá (21˚28'43''N; 077˚15'84''W) is frequented used for dragnet fishing, while station 3 in Júcaro which was used as reference site (21˚32'89''N; 077˚08'0''W), is a more pristine location deprived of direct influence of contaminants or fishing. This site has a muddy-sandy substrate with a variety of shells and rocks available for colonization ( Figure 1 were collected randomly to evaluate epiphytic macroalgae. After collection, a value of relative abundance was assigned to each species of macroalgae identified [19]. These values were used to calculate Shannon-Weaver diversity index [20] following a Jackknife resampling method. To evaluate the appropriateness of our sample size (i.e. minimum area), we created cumulative curves of the number of species vs. UM [21] [22].
Many studies had been done by other authors, which was to our advantage because classification was an important part of this paper [23]- [28]. Moreover, all macroalgal species were organized according to their morpho-functional grouping [9], season, and month when they were collected. Epiphytes were not grouped because of the majority of species corresponded to early stages of succession and Thalassia showing a high rate of leaves renewal, thus increasing the challenge of interpreting temporal patterns.
A fixed-effect two-factor ANOVA was used to analyze all data. Transformations (x 0.25 ) were performed as required to fulfill the assumptions of normality and homogeneity of variances per factor and level [29]. For all data, where sig- Moreover, the percentage dissimilarity coefficient [21] was used as an affinity index with the unprocessed abundance data. This included only species that accumulated 95% of the total number of relative abundance [30]. A spatial representation of the morpho-functional per across months and sampling stations grouping was obtained using a Non-Metric Multidimensional Scaling (NMMDS) approach computed with the same affinity index. These analyses were performed on the MVSP [31] and STATISTICA v. 7.0 [32] programs. The index of constancy and nodal fidelity were used to assess the correspondence between the grouping results by season and by species [21].

Results and Discussion
The curves of the number of species vs. UM showed that. The cumulative curves show an asymptotic tendency that stabilizes at 12 and 9 UMs, for the epiphytic and associated macroalgae, respectively. This confirms that the sample size used was optimal in defining the structure of the macroalgal communities ( Figure 2). [33], argues that the best dimension for UM should be the smallest possible and that random sampling ensures greater accuracy in the determination of the mean and variance.
Our results show 59 epiphytic or associated macroalgae species. These included 3 Cyanophyta, 25 Chlorophyta, eight Phaeophyceae, 23 Rhodophyta and 3 Magnoliophyta (Table 1). Moreover, the inverse dissimilarity coefficient (by species) allowed defining three distinct morpho-functional groups ( Figure 3, Table 1). All of them present in different proportions at the three sampling stations and represented with high fidelity according to the nodal analysis ( Table  2). The model of morpho-functional groups was developed to synthesize the variations in the morphology of the thallus as a response to different ecological conditions [11]. Analyzing the distribution of macroalgae within this context allows grouping species with greater fidelity per station, thus enabling to distinguish the dominance of specific morphotypes as a function of the environmental conditions per site.
The analysis of relative abundance shows significant differences between stations (F × 21.83, p × 0.001**) and months (F × 2.55, p < 0.038*), while the  Shannon-Weaver diversity index only shows differences between stations (F × 33.50, p < 0.001**) ( Figure 4). Cayo Puto station shows the highest abundance and diversity of macroalgae (Figure 4(a), Figure 4(b)) presumably because the site is sheltered from strong oceanic bulk water flow. These results coincide with observations by [34], who found that the number of species increases in areas with lower mechanical disturbance and rich in nutrients compared to areas with higher disturbance and low nutrient availability. Other authors, however, consider that low diversity is typical of areas with heavy nutrient loads [35]. Contrary to Cayo Puto, Playa Bagá station shows the lowest abundance and diversity of macroalgae (Figure 4(a), Figure 4(b)). This could be associated with the mechanical damage caused by trawls and dragnets increasing physical stress in the meadow [7]. The spatial variations were best explained by the differences nutrients and substratum between sampling sites [17].
On the other hand, Júcaro station shows the greatest number of associated species, while the epiphyte component was lower than in Cayo Puto. Moreover, October showed the lowest values of abundance across sites (Figure 4(c)). Temporal variations were only evident in terms of the relative abundance but not diversity, showing seasonal changes in the qualitative structure of macroalgal communities, where some species replace others across seasons. Similar results observed by [36].
The NMMDS constructed with percentage dissimilarity coefficient shows a clear separation of macroalgal communities as a function of sampling stations R. Cabrera et al.       Figure 6). These results coincide with macroalgal morphologies frequent in areas subjected to organic enrichment [35]. Similarly, [38] found that the genera Spyridia, Ulva, Bryopsis, and Dictyota predominate in nutrient-enriched areas, such as Cayo Puto. Furthermore, Gracilaria caudata and G. damaecornis, which were also present at site, are also typical of areas with high nutrient loads [38] [39]. These results coincide with who compared zones with different levels of degradation.
Our results highlight the option of using morpho-functional grouping as a tool to analyze the overall macroalgae community structure across sites [11] [12] [36]. Results also confirm that the spatial variation of the macroalgal community across seagrass meadows at Nuevitas is much more noticeable and decisive than the temporal variation. These variations appear to be determined, at least to some extent, by the degree of anthropogenic impact (i.e. nutrient enrichment and physical disturbance be trawling and dragnet fishing) but also due to abiotic factors such as water movement, as described in [7]. In Cuba, [41] recognized that seagrass ecosystems are very important for sustaining fisheries, along with coral reefs and mangroves.
Sites showing the highest physical disturbance have the lowest relative abundance and diversity of macroalgal, while a greater abundance and diversity was observed in stations with higher levels of nutrients.
Differences in the specific composition of macroalgal communities are determined by a spatial rather than a temporal component.
Macroalgae within the foliose and filamentous groups dominated in the station with the highest levels of nutrients. On the other hand, articulated-calcified and leathery macroalgae dominated stations with lower nutrients loads and physical disturbance.