Hostname: page-component-745bb68f8f-grxwn Total loading time: 0 Render date: 2025-02-06T06:44:50.652Z Has data issue: false hasContentIssue false

Standing stock variations, growth and CaCO3 production by the calcareous green alga Halimeda opuntia

Published online by Cambridge University Press:  30 August 2016

Pedro Bastos De Macedo Carneiro*
Affiliation:
Instituto de Ciências do Mar, Universidade Federal do Ceará, Av. Abolição, 3207, Meireles. CEP 60.165-081. Fortaleza, Ceará, Brazil
Jamile Ulisses Pereira
Affiliation:
Instituto de Ciências do Mar, Universidade Federal do Ceará, Av. Abolição, 3207, Meireles. CEP 60.165-081. Fortaleza, Ceará, Brazil
Helena Matthews-Cascon
Affiliation:
Centro de Ciências, Universidade Federal do Ceará, Campus do Pici. CEP 60440-554. Fortaleza, Ceará, Brazil
*
Correspondence should be addressed to:P.B.M. Carneiro, Instituto de Ciências do Mar, Universidade Federal do Ceará, Av. Abolição, 3207, Meireles. CEP 60.165-081. Fortaleza, Ceará, Brazil email: pedrocarneiro@ufc.br
Rights & Permissions [Opens in a new window]

Abstract

The present paper investigates standing stock variations of Halimeda opuntia on a sandstone reef of the South-west Atlantic Ocean, in order to better understand the role of this seaweed as a CaCO3 producer. The study was conducted over two 3-month periods, using photo quadrats to analyse the coverage area, and destructive sampling to quantify area-specific biomass and CaCO3 percentage. The alga occupied 2.4% of the substrate (4464 m2), growing as clumps with an average biomass of 1.59 kg m−2, resulting in a standing stock of 7097.8 kg of alga. This standing stock varied with environmental conditions, particularly wind speed. Assuming an exponential model for these variations, H. opuntia produced at least 13,050.14 kg (54.37 g m−2 day−1) of carbonate sediments. There was a positive correlation between changes in standing stock and coverage, but not with area-specific biomass. This suggests that net algal growth results in the occupation of new spaces, with minimal increases in height or segment density. Therefore monitoring coverage should complement traditional individual-based methods for estimating Halimeda growth and production. Combined, these approaches should result in more accurate models of the role of this alga on marine carbonate budgets.

Type
Research Article
Copyright
Copyright © Marine Biological Association of the United Kingdom 2016 

INTRODUCTION

Among marine calcium carbonate (CaCO3) producers, calcareous green algae of the genus Halimeda J.V. Lamouroux are one of the major contributors of CaCO3-rich sediments in tropical regions (Freile et al., Reference Freile, Milliman and Hillis1995; Rees et al., Reference Rees, Opdyke, Wilson and Henstock2006). These algae have widespread distribution, usually attain high coverage, and have brittle segmented thalli with high percentages of calcium carbonate (Hillis-Colinvaux, Reference Hillis-Colinvaux, Blaxter, Russel and Younge1980; Kooistra et al., Reference Kooistra, Coppejans and Payri2002). In many tropical continental shelves, large fractions of calcareous sediment deposits are composed by Halimeda dead fragments (Milliman, Reference Milliman and Flugel1977; Alexandersson & Milliman, Reference Alexandersson and Milliman1981). Due to its importance, the genus has been the subject of many ecological studies, and direct estimates of growth rates and CaCO3 production are available for many of its species (Multer, Reference Multer1988; Vroom et al., Reference Vroom, Smith, Coyer, Walters, Hunter, Beach and Smith2003; Mayakun et al., Reference Mayakun, Bunruk and Ko2014).

Within the genus, the rock-dwelling Halimeda opuntia (L.) J. V. Lamouroux is one of the most common species worldwide (Kooistra et al., Reference Kooistra, Coppejans and Payri2002; Bandeira-Pedrosa et al., Reference Bandeira-Pedrosa, Pereira and Oliveira2004). Due to its abundance, this alga may play a key ecological role on marine hard-bottom communities (Hillis-Colinvaux, Reference Hillis-Colinvaux, Blaxter, Russel and Younge1980; Bandeira-Pedrosa et al., Reference Bandeira-Pedrosa, Pereira and Oliveira2004). Nevertheless, direct measures of CaCO3 production by this species are rare (e.g. Drew, Reference Drew1983; Hudson, Reference Hudson, Toomey and Nitecki1985; Multer & Clavijo, Reference Multer and Clavijo2004), and the recognition of its importance on sediment dynamics is mostly based on indirect data, derived from sedimentological studies of carbonate build-ups (Enos & Perkins, Reference Enos and Perkins1977; Hine et al., Reference Hine, Hallock, Harris, Mullins, Belknap and Jaap1988; Johns & Moore, Reference Johns and Moore1988).

The paucity of studies on H. opuntia is partially due to its complex morphology, which differs from many of its congeners by the tendency of assuming sprawling forms, with indeterminate growth (Verbruggen & Kooistra, Reference Verbruggen and Kooistra2004; Yñiguez et al., Reference Yñiguez, McManus and Deangelis2008). In such cases it is often not feasible to isolate individuals within a population, hindering quantitative measurements of CaCO3 production, since these estimates are commonly based on monitoring mean individual growth (e.g. Bach, Reference Bach1979; Payri, Reference Payri1988; Garrigue, Reference Garrigue1991; Freile & Hillis, Reference Freile, Hillis, Lessios and MacIntyre1997; Vroom et al., Reference Vroom, Smith, Coyer, Walters, Hunter, Beach and Smith2003). The scarcity of quantitative studies in many locations prevents more accurate analyses on the role of H. opuntia on global and regional sediment budgets.

In the tropical South-west Atlantic Ocean, the continental margin is a mixed carbonate-siliclastic platform (Testa & Bosence, Reference Testa, Bosence, Wright and Burchette1998, Reference Testa and Bosence1999; Vital et al., Reference Vital, Stattegger, Amaro, Schwarzer, Frazão, Tabosa, Silveira, Hampson, Steel, Burgess and Dalrymple2008,  Reference Vital, Gomes, Tabosa, Frazão, Santos and Plácido2010; Bastos et al., Reference Bastos, Quaresma, Marangoni, D'Aagostini, Bourguignon, Cetto, Silva, Amado-Filho, Moura and Collins2015; Gomes et al., Reference Gomes, Vital, Eichlere and Sen-Gupta2015). In shallower waters, calcareous accumulations are frequently found near reef habitats (Segal & Castro, Reference Segal and Castro2011; Castro et al., Reference Castro, Segal, Negrão and Calderon2012; Silva et al., Reference Silva, Leão, Kikuchi, Costa and Souza2013). These carbonate deposits are mostly composed of skeletal debris of reefal CaCO3 producers, such as cnidarians, molluscs and algae, including Halimeda (Leão et al., Reference Leão, Dutra, Spanó, Dutra, Alen, Werner and McKenna2006; Pereira et al., Reference Pereira, Manso, Macedo, Dias and Silva2013). Understanding the ecology of the H. opuntia in this context is therefore paramount to comprehending regional sediment and carbonate dynamics, which so far have been poorly studied (e.g. Gherardi, Reference Gherardi2004; Araújo & Machado, Reference Araújo and Machado2008; Amado-Filho et al., Reference Amado-Filho, Moura, Bastos, Salgado, Sumida, Guth, Francini-Filho, Pereira-Filho, Abrantes, Brasileiro, Bahia, Leal, Kaufman, Kleypas, Farina and Thompson2012; Carneiro & Morais, Reference Carneiro and Morais2016).

The present study aims to analyse H. opuntia population dynamics on a SW Atlantic reef, by measuring temporal variations on coverage and area-specific biomass, as a means to quantify its importance as a space occupier. Furthermore, standing stock changes were used to estimate CaCO3 production in order to better characterize its role as a producer of carbonate sediments. To our knowledge, this is one of the first quantitative studies on biogenic production of sediments in South Atlantic reefs and the first to relate standing stock and carbonate productivity by Halimeda. This different approach results in net instead of gross CaCO3 production and seems to be suited to species with sprawling, indeterminate habits.

MATERIALS AND METHODS

Study area

The Parrachos de Pirangi are ferruginous sandstone outcrops, possibly of tertiary origin (Branner, Reference Branner1904), located ~1000 m, toward 90°E, off the mouth of the Pirangi river in the state of Rio Grande do Norte, Brazil (05°58.824′S 35°06.495′W).

The formation is composed of numerous flat sandstone elevations of varying sizes, which may or may not emerge during low tide, forming a discontinuous barrier parallel to the coastline (Figure 1). On previous visits to the study area, Halimeda opuntia was found growing only at the intertidal region of the largest of these platforms, an area of ~186,000 m2, that could be divided into three zones: a protected back zone with many submerged channels, a plateau with flat relief and a wave-exposed crest.

Fig. 1. Location of the study area. Dark grey region in the detail indicates the area of occurrence of Halimeda opuntia.

The climate of the region is tropical wet and dry (Aw in Köppen classification). The rainiest months are April to June and the driest October to December. Temperatures are usually high and with low amplitudes throughout the year. Average temperature varies from 24.5°C between June and August to 27°C in January–March. The wind speed is higher outside the rainy season, the average speed being the highest between August and October.

Field sampling

Three transects were established parallel to the major axis (roughly parallel to the coastline) of the reef where H. opuntia was found, one in each zone: back, plateau and crest. Each transect was 500 m long and separated from the adjacent transect by 75 m.

Six field campaigns were conducted to study the H. opuntia population, three in the dry season (August, September and October 2009) and three in the rainy season (March, April and May 2010).

On each campaign, the three transects were visited in periods of diurnal spring tide, when the area is mostly emerged. These transects were used to collect data on biomass and population dynamics (see below).

Area-specific biomass and calcium carbonate

In the present study, biomass was defined as the area-specific mass of the algae (in g cm−2), and was calculated dividing the algal dry weight (DW) by the area of the samples (100 cm2). In each campaign, 15 samples of H. opuntia were haphazardly collected, five per transect. Each sample consisted of algal clumps measuring at least 100 cm2, scraped from the substrate with the aid of 10 × 10 cm quadrats. This sample size was chosen to ensure sampling of mature individuals, which should cover the entire quadrat area to allow comparisons among months, and it is similar to the average clump size in the study area (see Results section). The samples were stored in plastic bags containing a solution of 4% formalin in seawater.

In the laboratory, samples were washed and cleaned to remove epibionts and dried at 50°C until their weight was constant. To measure CaCO3 content, dry samples were decalcified with 0.1M HCl. After the reaction, samples were carefully rinsed with distilled water and the liquid phase was discarded to remove any dissolved salts. The material was once again dried at 50°C until constant weight and the CaCO3 percentage was calculated as the difference between the DW of the fresh sample and that of the decalcified sample (Van Tussenbroek & Van Dijk, Reference Van Tussenbroek and Van Dijk2007).

Population dynamics

Coverage was defined as the average reef area covered by H. opuntia (in m2). To estimate this parameter, 30 quadrats measuring 50 × 50 cm were evenly distributed over the reef (one quadrat every 50 m, in total 10 per transect) and were photographed using a digital camera with the aid of mounted supports to avoid angle variations. This quadrat size was selected to ensure that no sample was smaller than the largest Halimeda clump observed on a pilot survey (~1200 cm2). To avoid bias in coverage estimates, only the Halimeda observed inside quadrats was considered in the analysis, even if the clumps extended beyond the delimited area. The digital images were analysed on ImageJ software, which was used to visually delimit and measure the area occupied by H. opuntia. Since our data did not show significant differences in area-specific biomass (see Results section), this 2D projected area was considered a good proxy to clump size and to the area occupied by the alga on the reef.

For the purpose of calculation, in campaigns where H. opuntia was observed on the reef but none of the 30 quadrats contained any clumps we used a percentage occupancy estimate that was 10% of the smallest coverage measured during the study.

Four environmental variables were used to explain eventual variations in population: air temperature (°C), photosynthetically active radiation (PAR) (μmol photons m−2 s−1), rainfall (mm) and wind speed (m s−1). These variables were chosen to represent those used by Yñiguez et al. (Reference Yñiguez, McManus and Deangelis2008), with rainfall as an estimator of the nutrient input to the system and wind speed as a measurement of disturbance by hydrodynamic forces, which seems reasonable for the study area (Barros & Rocha-Barreira, Reference Barros and Rocha-Barreira2014). To reduce bias, environmental variables were measured in the 30 days immediately prior to each campaign, therefore: temperature and wind speed were the average daily value along this period; PAR was the average value per hour of sunlight; and rainfall was the rain accumulated in these 30 days.

These data were obtained from the solarimetric station of the Laboratório de Variáveis Ambientais Tropicais (LAVAT) of the Instituto Nacional de Pesquisas Espaciais (INPE), located in Brazil.

Halimeda opuntia net growth

Standing stock was defined as the total amount of algae on the reefs (in kg), which was calculated by multiplying the coverage by the area-specific biomass estimated in each campaign. The net growth rate was the rate of variation of the standing stock between consecutive months.

We have assumed that standing stock variations were exponential. This assumption results in simple but useful models (Hastings, 1996; Yong et al., Reference Yong, Yong and Anton2013), that are widely used in seaweed studies (e.g. DeBoer et al., Reference DeBoer, Guigli, Israel and D'Elia1978; Littler & Littler, Reference Littler and Littler1985; Pedersen & Borum, Reference Pedersen and Borum1997; Martins et al., Reference Martins, Oliveira, Flindt and Marques1999; Barr et al., Reference Barr, Kloeppel, Rees, Scherer, Taylor and Wenzel2008; Hadley et al., Reference Hadley, Wild-Allen, Johnson and Macleod2015), particularly when density-dependent effects do not need to be considered, such as when population sizes are small or highly variable, which is the case in the present study (see Results section). A more complex approach would require estimates of the carrying capacity of tropical intertidal environments, which are not readily available and would be difficult to determine due to the complexity of this type of ecosystem (Edgar, Reference Edgar1993; Christensen & Pauly, Reference Christensen and Pauly1998; Monte-Luna et al., Reference Monte-Luna, Brook, Zetina-Rejón and Cruz-Escalona2004). Furthermore, other studies have shown that H. opuntia may cover 90–100% of the substrate in some locations (Hillis-Colinvaux, Reference Hillis-Colinvaux, Blaxter, Russel and Younge1980), which indicates that the species is able to flourish into very dense aggregates. As such, density-dependent effects should be an issue only at extreme coverages. Therefore we have considered that exponential models are reasonable approximations, sufficiently accurate to describe the Halimeda population under study.

According to this model, the standing stock and the coverage at a given time are calculated by the formulae:

(1) $$S_{\rm f}^{\,\rm s} = S_{\rm o}^{\,\rm s} \big(1 + r_{\rm s} \big)^{\rm t}$$
(2) $$S_{\rm f}^{\,\rm c} = S_{\rm o}^{\,\rm c} \big(1 + r_{\rm c} \big)^{\rm t}$$

where $S_{\rm f}^{\,\rm s} $  = final standing stock, $S_{\rm o}^{\,\rm s} $  = initial standing stock, r s = net growth rate, t = time in number of days, $S_{\rm f}^{\rm c} $  = final coverage, $S_{\rm o}^{\rm c} $  = initial coverage and r c = rate of variation in coverage (rate of area growth). Rearranging either equation (1) or (2), the net growth rate is:

(3) $$r = \big(S_{\rm f} /S_0 \big)^{{\rm 1/t}} {\rm -} 1$$

which is a reliable estimate of seaweed growth when compared with other formulae (Yong et al., Reference Yong, Yong and Anton2013). Finally, the average daily growth per square metre (χ) is:

(4) $$\chi = \displaystyle{1 \over t}\sum\limits_{i = 0}^t {\displaystyle{{S_i^{\,\rm s} r_{\rm s}} \over {S_i^{\,\rm c} }}}$$

which represents how much algal mass was lost per day, after adjusting for the respective variations in coverage.

Statistical analyses

Due to the non-normal skewed nature of abundance data, a Kruskal–Wallis test was used to compare reef zones in terms of H. opuntia coverage. An Analysis of Variance (ANOVA) followed by a Tukey HSD test in cases where the P-value was <0.05 was used to examine whether there was a significant variation in area-specific biomass or percentage of CaCO3. The normality of the data was analysed by Lilliefors test and a log transformation was applied in the cases where data seemed to be non-normal.

The relation of the environmental variables with area-specific biomass, % CaCO3, coverage and standing stock was examined via Pearson correlation analysis, with the data log-transformed (log10 + 1).

A linear regression analysis was used to understand which environmental variable could explain most of the variation in the standing stock. Considering that the environmental variables are often correlated, to prevent multicollinearity we performed a Principal Component Analysis (PCA) of these variables and used the Components with eigenvalue higher than one in the multiple regression.

All tests were performed using R 3.1.1 (R Foundation for Statistical Computing, Vienna, Austria).

RESULTS

Halimeda opuntia population

Throughout the study period, H. opuntia was the only species belonging to the Halimeda genus that was observed growing on the studied reefs. No evidence was found that these clusters were more common in any particular transect (Kruskal–Wallis χ2 = 3.09, P = 0.21). Therefore, it seems that Halimeda may grow over their entire area of occurrence, without preference for any specific reef zone (Figure 2).

Fig. 2. Area of occurrence of Halimeda opuntia. The black circles illustrate the average area of coverage at each point sampled indicating no differences between transects or reef areas.

Population parameters are summarized in Table 1. When present on the reef, H. opuntia appeared as aggregates with an average of 2.67 ± 0.4 clumps m−2 (mean ± SEM), each clump measuring 99.2 ± 18.0 cm2. The average biomass was 1.59 ± 0.13 kg m−2, most of it (84.06 ± 0.6%) being calcium carbonate.

Table 1. Summary statistics (mean ± SEM) of Halimeda opuntia J.V. Lamouroux on Brazilian sandstone reefs.

a Calculated by multiplying mean coverage and mean biomass.

b None of the quadrats had tufts in October.

c For calculation purposes, coverage in October was considered to be 10% of the smallest estimated coverage.

The percentage of CaCO3 differed between months (F = 24.99, P < 0.001). The Tukey HSD test did not show a clear pattern, grouping August and September (highest % CaCO3); October, March and May (lowest % CaCO3); and April (intermediate % CaCO3), every group differing significantly to each other.

Halimeda opuntia net growth and sediment production

Halimeda opuntia occupied an average of 2.4% of its area of occurrence (4464 m2). This area, however, varied between the months, from 2758.4 m2 in September to 8799.8 m2 in August (Figure 3). In October, as none of the 30 quadrats contained any clumps, we considered a coverage area of 275.8 m2 for calculation purposes.

Fig. 3. Area covered by Halimeda opuntia in each month of the study.

Multiplying average coverage by the area-specific biomass, it was estimated that there was on average 7097.8 kg DW of H. opuntia over the reef. Variations in this standing stock correlated positively with coverage (r = 0.90, P = 0.01), but not with biomass (r = 0.52, P = 0.10), which have remained relatively stable throughout the study period (Figure 4).

Fig. 4. Temporal variation in standing stock, coverage and biomass of Halimeda opuntia on Brazilian reefs.

Considering that during the study period the standing stock of H. opuntia decreased, the net growth rate estimated by equation (3) was −0.07% of the total amount of algae on the reefs per day. Applying equation (4), this equates to an average net loss of 1.0 gDW m−2 day−1.

This rate, however, was not constant throughout the study period. The maximum daily net growth observed was 75.2 g DW m−2 day−1 between March and April. On the other hand, between September and October, there was an average net loss of 72.2 g DW m−2 day−1 (Figure 5).

Fig. 5. Temporal variation of average growth rate and daily net growth of Halimeda opuntia. Daily net growth represents how much algal mass was lost per day after adjusting for the respective area variations.

If we consider only the periods when the population declined (August–September, September–October and April–May) and that 84.06% of the dry mass of H. opuntia was CaCO3, it was estimated that over the study period, this alga alone contributed toward at least 13,050.14 kg (54.37 g m−2 day−1) of carbonate sediments to the surrounding reef areas.

Environmental variables and population dynamics

As expected by the climate in the study area, the 30 days preceding the field campaigns between August and October had on average higher wind speeds and PAR, but lower rainfall intensities and air temperature than the 30 days leading collections between March and May (Figure 6).

Fig. 6. Average of six environmental variables in the 30 days prior to the campaigns to study Halimeda opuntia on Brazilian reefs.

The correlation analysis showed a strong positive relation between rainfall and both coverage (r = 0.90, P = 0.01) and standing stock (r = 0.88, P = 0.02). On the other hand, these two variables had a strong negative relation with wind (coverage: r = −0.83, P = 0.04; standing stock: r = −0.88, P = 0.02) (Figure 7). None of the other variables showed significant relationships, with the exception of a negative correlation between rainfall and PAR (r = −0.86, P = 0.03).

Fig. 7. Scatterplot showing a negative relationship between Halimeda opuntia standing stock and wind speed (r = −0.88, P = 0.02), but a positive relation with rainfall (r = 0.88, P = 0.02).

Only the first two components generated by the PCA had eigenvalue greater than one, and together represented 87.34% of the variance of the environmental variables. The plot of environmental variables in the space formed by the two components (Figure 8) showed that PC1 separated PAR and wind speed from rainfall, with minor influence from air temperature. Thus it can be understood as representing climatic variations between rainy and dry season in the study area, which result in increased agitation in shallow marine environments. Based on the same reasoning, PC2 could represent climatic variations of the transition from the coldest to the warmest months, which could lead to desiccation and heat stress.

Fig. 8. Plot of four environmental variables on the space formed by the first and second Principal Components Analysis. PAR, photosynthetically active radiation.

The regression analysis of the standing stock with the first two PCs generated a model with a low but non-significant P-value (F = 9.17, P = 0.052). However, considering that the t-test showed a significant coefficient for the first PC (t = −4.14, P = 0.02) and that the two components are by definition orthogonal (non-correlated), we interpreted this non-significant F-test as a consequence of the inclusion of the PC2, which seems to be unrelated to the standing stock (t = 1.06, P = 0.36) and, as so, could negatively affect the result of the model. Because of that, we excluded this second component from the model and repeated the regression only with PC1.

This second linear regression with only the first principal component and the standing stock was significant (F = 16.65, P = 0.02). Both the intercept (P < 0.001) and the linear coefficient (P = 0.01) were also significant. The model formula was:

$$S^{\,\rm s} = 3.65{\rm -} 0.36{\rm PC}$

where S s = standing stock and PC = first principal component. The 95% confidence limits for the intercept are from 3.321 to 3.990, and for the linear coefficient from −0.608 to −0.116.

DISCUSSION

The role of Halimeda on calcium carbonate production is well established, and the genus is considered an important player in the global carbonate budget (Freile et al., Reference Freile, Milliman and Hillis1995; Rees et al., Reference Rees, Opdyke, Wilson and Henstock2006). In accordance with the genus standards, the present study has detected variations in abundance of Halimeda opuntia which suggests that it is a fast growing species and an important sediment producer in tropical reefs. Even though it did not cover large areas of the substrate (i.e. we have estimated a coverage of 2.4% of the area), the variations in standing stock were responsible for at least 13,000 kg of CaCO3-rich sediments to the reef habitat and its surroundings.

Regarding the influence of environmental conditions on H. opuntia, our results agree with the model proposed by Yñiguez et al. (Reference Yñiguez, McManus and Deangelis2008), which have indicated that disturbance is one of the main factors structuring its populations. In the present study the population declined in times of stronger winds, which determine more intense waves and currents, and consequently result in increased disturbances (Vital et al., Reference Vital, Stattegger, Amaro, Schwarzer, Frazão, Tabosa, Silveira, Hampson, Steel, Burgess and Dalrymple2008, Reference Vital, Gomes, Tabosa, Frazão, Santos and Plácido2010). There was also a significant correlation with the first principal component, which was interpreted as summarizing agitation-mediated disturbances of abiotic origin. Such conditions are also known to negatively affect seagrass beds in the SW Atlantic (Barros & Rocha-Barreira, Reference Barros and Rocha-Barreira2014). Intense wind-driven disturbance seems to affect coverage by breaking algal clumps, while preventing fragment resettlement (Yñiguez et al., Reference Yñiguez, McManus and Deangelis2008).

Despite this likely importance of wind-driven disturbance in structuring H. opuntia populations, environmental variables are typically correlated in the study area, and it is not always easy to isolate the effects of a single factor on the population. The PCA analysis suggests that temperature may have had only a minor effect on Halimeda standing stock, whereas rainfall significantly benefited the population. During rainy seasons, the average wind speed usually decreases and this benefit could have been the result of a reduced disturbance. Nevertheless, considering that our data did not show a linear correlation between rain and wind, we believe that the positive effect of rainfall was mainly due to increased nutrient inputs. The CaCO3 percentage of thalli supports this interpretation, since algae tended to be less calcified in rainy months (Table 1), indicating an active production of new segments. Such a reduction in calcium carbonate would not be expected if the larger coverage was due mainly to a less disturbed season. Furthermore, there is a river near the reefs, which is a natural source of fresh water to the system, increasing nutrient supply during rainy seasons (Rodrigues et al., Reference Rodrigues, Knoppers, Souza and Santos2009; Seitzinger et al., Reference Seitzinger, Mayorga, Bouwman, Kroeze, Beusen, Billen, Van Drecht, Dumont, Fekete, Garnier and Harrison2010; Dias et al., Reference Dias, Castro and Lacerda2013).

Moreover, new studies are needed to assess the effects of PAR on H. opuntia. In our analysis, this variable did not have an influence on the population. Nevertheless high levels of radiation may have inhibited photosynthesis (Littler et al., Reference Littler, Littler and Lapointe1988; Payri, Reference Payri1988), reducing growth rates. This could have led to an accentuation of the destructive effects of disturbance, even if PAR did not affect the population directly.

Regarding H. opuntia growth, most of the previous data came either from in vitro analyses (Drew & Abel, Reference Drew, Abel, Harmelin-Vivien and Salvat1985; Payri, Reference Payri1988) or from monitoring individual colonies on shallow subtidal environments (Drew, Reference Drew1983; Hudson, Reference Hudson, Toomey and Nitecki1985; Multer & Clavijo, Reference Multer and Clavijo2004). These studies have shown that the species grows through spasmodic events, which may be strongly seasonal, wherein algae may produce segments luxuriantly, but unevenly over their thalli. The present study adds information on intertidal populations growing on reef flats, and it suggests temporal variations on coverage, but a relatively stable area-specific biomass of individual clumps. This indicates that most of the growth results in an increased coverage, the algae occupying new spaces with relatively minor increases in height or segment density.

Since H. opuntia may assume a laxer morphology in shallow subtidal habitats, exhibiting looser and longer branches (Hillis-Colinvaux, Reference Hillis-Colinvaux, Blaxter, Russel and Younge1980; Littler & Littler, Reference Littler and Littler2000), new studies are needed to estimate area-specific biomass changes in these environments. Furthermore, new analyses are needed to understand growth in younger and smaller individuals, since we have restricted our biomass observations to mature algae (with at least 100 cm2). Nevertheless, due to the limited height that is achievable by this species in field conditions, and based on descriptions by other authors (e.g. Hillis-Colinvaux, Reference Hillis-Colinvaux, Blaxter, Russel and Younge1980; Multer and Clavijo, Reference Multer and Clavijo2004), sprawling morphotypes seem to be the rule rather than the exception for the species. Hence our results should be extendable to other populations under similar environmental conditions, irrespective of their average clump sizes.

Furthermore, previous studies have also pointed to the same relationship between variations in standing stock and coverage in Halimeda, with a relatively minor importance of individual sizes or biomass (Van Tussenbroek & Van Dijk, Reference Van Tussenbroek and Van Dijk2007). This seems to be particularly true for populations under disturbance – which seems to be the case in the present study – where fragmentation may be a rapid source of new individuals (Walters & Smith, Reference Walters and Smith1994; Yñiguez et al., Reference Yñiguez, McManus and Deangelis2008), since segment resettlement seems to be an important mechanism for dispersion in Halimeda (Walters et al., Reference Walters, Coyer, Hunter, Beach and Vroom2002). Nevertheless, population growth by reproduction and space colonization remains poorly measured in this alga. Most analyses of production are based on counting the number of new segments produced by previously marked individuals over a given period of time (Bach, Reference Bach1979; Payri, Reference Payri1988; Garrigue, Reference Garrigue1991; Freile & Hillis, Reference Freile, Hillis, Lessios and MacIntyre1997; Vroom et al., Reference Vroom, Smith, Coyer, Walters, Hunter, Beach and Smith2003; Carneiro & Morais, Reference Carneiro and Morais2016). This approach, because it does not consider the effect of propagation, may have been underestimating Halimeda contributions to both CaCO3 and sediment production.

New studies are also necessary to assess the role of sexual reproduction in the process of space occupation and CaCO3 production by this alga. All Halimeda species undergoes holocarpic sexual reproduction, in which the cellular content is converted into gametes, resulting in algal death, total dismantling of thallus and segment release (Hillis-Colinvaux, Reference Hillis-Colinvaux, Blaxter, Russel and Younge1980). The literature on Halimeda suggests that only a small fraction of the population is actively reproducing at any given time (Clifton & Clifton, Reference Clifton and Clifton1999; Vroom et al., Reference Vroom, Smith, Coyer, Walters, Hunter, Beach and Smith2003). Nevertheless, some authors have pointed to the coincidence between periods of sexual reproduction and population declines (Clifton & Clifton, Reference Clifton and Clifton1999). During our field campaigns, no specimen was found with clear signs of sexual activity, such as visible gametangia. However, the fast population reduction, as well as the rapid subsequent recovery, observed around October may have been caused, at least in part, by episodes of sexual reproduction. As such, due to their potential to produce segment release, these events may be an important source of CaCO3 in reef environments.

The exponential function used in the analysis of standing stock variations is simple and may be further improved to better capture the complexity of this alga. However, it seems to be adequate to describe a typical population that is far from the ecosystem's carrying capacity (Weiner et al., Reference Weiner, Kinsman and Williams1998). This seems to be the case in the present study, due to the low mean coverage observed. Moreover, we have only aimed to quantify standing stock variations and test some hypotheses on their causes, and not to develop full models of H. opuntia population dynamics.

Additionally, estimating population net growth through coverage variations, besides being suited to H. opuntia and other sprawling species, is a way of detecting increases in population size due to reproduction and colonization, refining CaCO3 production estimates by these algae. Furthermore, this approach may overcome other limitations of individual methodologies, such as the impossibility of detecting population reductions, which may happen even if some algae are actively growing (Multer & Clavijo, Reference Multer and Clavijo2004). On the other hand, this methodology lacks the capacity to quantify total growth, which could, at least theoretically, be done with individual measurements of growth. Therefore, instead of an alternative methodology, estimates of coverage variation should be used to complement Halimeda growth measurements, resulting in more precise estimates of its CaCO3 and biomass production.

In conclusion, the quantitative analyses on the present study, based on standing stock measurements, support the view that H. opuntia has an active role in sediment production in tropical intertidal reef environments. Nevertheless, they also indicate that reproduction and dispersion should be better quantified in studies of calcium carbonate production. Since growth in such environments seems to be strongly related to substrate occupation, as opposed to increases in algae height or area-specific biomass, more accurate measurements of Halimeda standing stock and new models of population dynamics may be useful tools in estimating current CaCO3 production by this alga, improving the knowledge on their role on marine carbonate budgets.

ACKNOWLEDGEMENTS

The authors are deeply indebted to Aline Martinez, Antônio Xavier, Bruna Sarkis, Camila Juanes, Caroline Castro, Cristina Rocha-Barreira, Daniel da Silva, Eurico Neto, Hortência Gonzales, João Correia Jr. Lívia Shell, Luciana Faustino, Murillo Ribeiro, Nathane Freitas, Rafael Rosseto, Rhayan Ramalho, Sônia Pereira, Tatiane Garcia, Victor Azevedo and Wilson Franklin for their help with field activities or data processing. We would also like to thank two anonymous reviewers, whose comments and suggestions greatly improved this work.

References

REFERENCES

Alexandersson, E.T. and Milliman, J.D. (1981) Intragranular Mg-Calcite cement in Halimeda plates from the Brazilian Continental Shelf. Journal of Sedimentary Research 51, 13091314.Google Scholar
Amado-Filho, G.M., Moura, R.L., Bastos, A.C., Salgado, L.T., Sumida, P.Y., Guth, A.Z., Francini-Filho, R.B., Pereira-Filho, G.H., Abrantes, D.P., Brasileiro, P.S., Bahia, R.G., Leal, R.N., Kaufman, L., Kleypas, J.A., Farina, M. and Thompson, F.L. (2012) Rhodolith beds are major CaCO3 bio-factories in the tropical South West Atlantic. PLoS ONE 7, e35171.CrossRefGoogle ScholarPubMed
Araújo, H.A.B. and Machado, A.D.J. (2008) Benthic Foraminifera associated with the South Bahia Coral Reefs, Brazil. Journal of Foraminiferal Research 38, 2338.CrossRefGoogle Scholar
Bach, S.D. (1979) Standing crop, growth and production of calcareous Siphonales (Chlorophyta) in a South Florida Lagoon. Bulletin of Marine Science 29, 191201.Google Scholar
Bandeira-Pedrosa, M.E., Pereira, S.M.B. and Oliveira, E.C. (2004) Taxonomy and distribution of the green algal genus Halimeda (Bryopsidales, Chlorophyta) in Brazil. Revista Brasileira de Botânica 27, 363377.Google Scholar
Barr, N.G., Kloeppel, A., Rees, T.A.V., Scherer, C., Taylor, R.B. and Wenzel, A. (2008) Wave surge increases rates of growth and nutrient uptake in the green seaweed Ulva pertusa maintained at low bulk flow velocities. Aquatic Biology 3, 179186.CrossRefGoogle Scholar
Barros, K.V.S. and Rocha-Barreira, C.A. (2014) Influence of environmental factors on a Halodule wrightii Ascherson meadow in northeastern Brazil. Brazilian Journal of Aquatic Sciences and Technology 18, 3141.CrossRefGoogle Scholar
Bastos, A.C., Quaresma, V.S., Marangoni, M.B., D'Aagostini, D.P., Bourguignon, S.N., Cetto, P.H., Silva, A.E., Amado-Filho, G.M., Moura, R.L. and Collins, M. (2015) Shelf morphology as an indicator of sedimentary regimes: a synthesis from a mixed siliciclastic–carbonate shelf on the eastern Brazilian margin. Journal of South American Earth Sciences 63, 125136.CrossRefGoogle Scholar
Branner, J.C. (1904) The stone reefs of Brazil, their geological and geographical relations, with a chapter on the coral reefs. Bulletin of the Museum of Comparative Zoology 44, 1285.Google Scholar
Carneiro, P.B.M. and Morais, J.O. (2016) Carbonate sediment production in the equatorial continental shelf of South America: Quantifying Halimeda incrassata (Chlorophyta) contributions. Journal of South American Earth Sciences 72, 16.CrossRefGoogle Scholar
Castro, C.B., Segal, B., Negrão, F. and Calderon, E.N. (2012) Four-year monthly sediment deposition on turbid southwestern Atlantic coral reefs, with a comparison of benthic assemblages. Brazilian Journal of Oceanography 60, 4963.CrossRefGoogle Scholar
Christensen, V. and Pauly, D. (1998) Changes in models of aquatic ecosystems approaching carrying capacity. Ecological Applications 8, S104S109.CrossRefGoogle Scholar
Clifton, K.E. and Clifton, L.M. (1999) The phenology of sexual reproduction by green algae (Bryopsidales) on Caribbean coral reefs. Journal of Phycology 35, 2434.CrossRefGoogle Scholar
DeBoer, J.A., Guigli, H.J., Israel, T.L. and D'Elia, C.F. (1978) Nutritional studies of two red algae. 1. Growth rate as a function of nitrogen source and concentration. Journal of Phycology 14, 261266.CrossRefGoogle Scholar
Dias, F.J.S., Castro, B.M. and Lacerda, L.D. (2013) Continental shelf water masses off the Jaguaribe River (4S), northeastern Brazil. Continental Shelf Research 66, 123135.CrossRefGoogle Scholar
Drew, E.A. (1983) Halimeda biomass, growth rates and sediment generation on reefs in the central great barrier reef province. Coral Reefs 2, 101110.CrossRefGoogle Scholar
Drew, E.A. and Abel, K.M. (1985) Biology, sedimentology and geography of the vast inter-reefal Halimeda meadows within the Great Barrier Reef Province. In Harmelin-Vivien, M. and Salvat, B. (eds) Proceedings of the Fifth International Coral Reef Congress, Ecole Pratique des Hautes Etudes, Tahiti, 27 May-01 June 1985, Volume 5. Tahiti: Antenne Museum-Ephe, pp. 1520.Google Scholar
Edgar, G.J. (1993) Measurement of the carrying capacity of benthic habitats using a metabolic-rate based index. Oecologia 95, 115121.CrossRefGoogle ScholarPubMed
Enos, P. and Perkins, R.D. (1977) Quaternary sedimentation in south Florida. Boulder, CO: Geological Society of America.Google Scholar
Freile, D. and Hillis, L. (1997) Carbonate productivity by Halimeda incrassata in a land proximal lagoon, Pico Feo, San Blas, Panama. In Lessios, H.A. and MacIntyre, I.G. (eds) Proceedings Eighth International Coral Reef Symposium, University of Panama, Panama, 24–29 June 1996, Volume 1. Panama City: Smithsonian Tropical Research Institute, pp. 767771.Google Scholar
Freile, D., Milliman, J.D. and Hillis, L. (1995) Leeward bank margin Halimeda meadows and draperies and their sedimentary importance on the western Great Bahama Bank slope. Coral Reefs 14, 2733.CrossRefGoogle Scholar
Garrigue, C. (1991) Biomass and production of two Halimeda species in the southwest New Caledonian lagoon. Oceanologica Acta 14, 581588.Google Scholar
Gherardi, D.F.M. (2004) Community structure and carbonate production of a temperate rhodolith bank from Arvoredo Island, southern Brazil. Brazilian Journal of Oceanography 52, 207224.CrossRefGoogle Scholar
Gomes, M.P., Vital, H., Eichlere, P.P.B. and Sen-Gupta, B.K. (2015) The investigation of a mixed carbonate-siliciclastic shelf, NE Brazil: side-scan sonar imagery, underwater photography, and surface-sediment data. Italian Journal of Geosciences 134, 922.CrossRefGoogle Scholar
Hadley, S., Wild-Allen, K., Johnson, C. and Macleod, C. (2015) Modeling macroalgae growth and nutrient dynamics for integrated multi-trophic aquaculture. Journal of Applied Phycology 27, 901916.CrossRefGoogle Scholar
Hastings (1996) Population biology: concepts and models. New York, NY: Springer-Verlag.Google Scholar
Hillis-Colinvaux, L. (1980) Ecology and taxonomy of Halimeda: primary producer of coral reefs. In Blaxter, J.H.S., Russel, F.S. and Younge, M. (eds) Advances in marine biology. London: Academic Press.Google Scholar
Hine, A.C., Hallock, P., Harris, M.W., Mullins, H.T., Belknap, D.F. and Jaap, W.C. (1988) Halimeda bioherms along an open seaway: Miskito Channel, Nicaraguan Rise, SW Caribbean Sea. Coral Reefs 6, 173178.CrossRefGoogle Scholar
Hudson, J.H. (1985) Growth rate and carbonate production in Halimeda opuntia: Marquesas Keys, Florida. In Toomey, D. and Nitecki, M. (eds) Paleoalgology: contemporary research and applications. Berlin: Springer-Verlag, pp. 257263.CrossRefGoogle Scholar
Johns, H.D. and Moore, C.H. (1988) Reef to basin sediment transport using Halimeda as a sediment tracer, Grand Cayman Island, West Indies. Coral Reefs 6, 187193.CrossRefGoogle Scholar
Kooistra, W.H.C.F., Coppejans, E.G.G. and Payri, C. (2002) Molecular systematics, historical ecology, and phylogeography of Halimeda (Bryopsidales). Molecular Phylogenetics and Evolution 24, 121138.CrossRefGoogle ScholarPubMed
Leão, Z.M.A.N., Dutra, L.X.C. and Spanó, S. (2006) The characteristics of bottom sediments. In Dutra, G.F., Alen, G.R., Werner, T. and McKenna, S.A. (eds) A rapid marine biodiversity assessment of the Abrolhos Bank, Bahia, Brazil. RAP Bulletin of Biological Assessment 38. Washington, DC: Conservation International, pp. 7581.Google Scholar
Littler, D.S. and Littler, M.M. (2000) Caribbean reef plants: an identification guide to the reef plants of the Caribbean, Bahamas, Florida and Gulf of Mexico. Washington, DC: Offshore Graphics.Google Scholar
Littler, M.M. and Littler, D.S. (1985) Ecological field methods: macroalgae. New York, NY: Cambridge University Press.Google Scholar
Littler, M.M., Littler, D.S. and Lapointe, B.E. (1988) A comparison of nutrient- and light-limited photosynthesis in psammophytic versus epilithic forms of Halimeda (Caulerpales, Halimedaceae) from the Bahamas. Coral Reefs 6, 219225.CrossRefGoogle Scholar
Martins, I., Oliveira, J.M., Flindt, M.R. and Marques, J.C. (1999) The effect of salinity on the growth rate of the macroalgae Enteromorpha intestinalis (Chlorophyta) in the Mondego estuary (west Portugal). Acta Oecologica 20, 259265.CrossRefGoogle Scholar
Mayakun, J., Bunruk, P. and Ko, R. (2014) Growth rate and calcium carbonate accumulation of Halimeda macroloba Decaisne (Chlorophyta: Halimedaceae) in Thai waters. Songklanakarin Journal of Science and Technology 36, 419423.Google Scholar
Milliman, J.D. (1977) Role of calcareous algae in Atlantic continental margin sedimentation. In Flugel, E. (ed.) Fossil algae: recent results and developments. New York, NY: Springer, pp. 232247.CrossRefGoogle Scholar
Monte-Luna, P., Brook, B.W., Zetina-Rejón, M.J. and Cruz-Escalona, V.H. (2004) The carrying capacity of ecosystems. Global Ecology and Biogeography 13, 485495.CrossRefGoogle Scholar
Multer, H.G. (1988) Growth rate, ultrastructure and sediment contribution of Halimeda incrassata and Halimeda monile, Nonsuch and Falmouth Bays, Antigua, W.I. Coral Reefs 6, 179186.CrossRefGoogle Scholar
Multer, H.G. and Clavijo, I. (2004) Halimeda investigations: progress and problems. Miami, FL: NOAA/RSMAS.Google Scholar
Payri, C.E. (1988) Halimeda contribution to organic and inorganic production in a Tahitian reef system. Coral Reefs 6, 251262.CrossRefGoogle Scholar
Pedersen, M.F. and Borum, J. (1997) Nutrient control of estuarine macroalgae: growth strategy and the balance between nitrogen requirements and uptake. Marine Ecology Progress Series 161, 155163.CrossRefGoogle Scholar
Pereira, N.S., Manso, V.A.V., Macedo, R.J.A., Dias, J.M.A. and Silva, A.M.C. (2013) Detrital carbonate sedimentation of the Rocas Atoll, South Atlantic. Anais da Academia Brasileira de Ciências 85, 5772.CrossRefGoogle Scholar
Rees, S.A., Opdyke, B.N., Wilson, P.A. and Henstock, T.J. (2006) Significance of Halimeda bioherms to the global carbonate budget based on a geological sediment budget for the Northern Great Barrier Reef, Australia. Coral Reefs 26, 177188.CrossRefGoogle Scholar
Rodrigues, R.P., Knoppers, B.A., Souza, W.F.L. and Santos, E.S. (2009) Suspended matter and nutrient gradients of a small-scale river plume in Sepetiba Bay, SE-Brazil. Brazilian Archives of Biology and Technology 52, 503512.CrossRefGoogle Scholar
Segal, B. and Castro, C.B. (2011) Coral community structure and sedimentation at different distances from the coast of the Abrolhos Bank, Brazil. Brazilian Journal of Oceanography 59, 119129.CrossRefGoogle Scholar
Seitzinger, S.P., Mayorga, E., Bouwman, A.F., Kroeze, C., Beusen, A.H.W., Billen, G., Van Drecht, G., Dumont, E., Fekete, B.M., Garnier, J. and Harrison, J.A. (2010) Global river nutrient export: a scenario analysis of past and future trends. Global Biogeochemical Cycles 24, GB0A08.CrossRefGoogle Scholar
Silva, A.S., Leão, Z.M.A.N., Kikuchi, R.K.P., Costa, A.B. and Souza, J.R.B. (2013) Sedimentation in the coastal reefs of Abrolhos over the last decades. Continental Shelf Research 70, 159167.CrossRefGoogle Scholar
Testa, V. and Bosence, D.W.J. (1998) Carbonate-siliciclastic sedimentation on a high-energy, ocean-facing, tropical ramp, NE Brazil. In Wright, V.P. and Burchette, T.P. (eds) Carbonate ramps. London: The Geological Society, pp. 5571.Google Scholar
Testa, V. and Bosence, D.W.J. (1999) Physical and biological controls on the formation of carbonate and siliciclastic bedforms on the north-east Brazilian shelf. Sedimentology 46, 279301.CrossRefGoogle Scholar
Van Tussenbroek, B.I. and Van Dijk, J.K. (2007) Spatial and temporal variability in biomass and production of psammophytic Halimeda incrassata (Bryopsidales, Chlorophyta) in a Caribbean reef lagoon. Journal of Phycology 43, 6977.CrossRefGoogle Scholar
Verbruggen, H. and Kooistra, W.H.C.F. (2004) Morphological characterization of lineages within the calcified tropical seaweed genus Halimeda (Bryopsidales, Chlorophyta). European Journal of Phycology 39, 213228.CrossRefGoogle Scholar
Vital, H., Gomes, M.P., Tabosa, W.F., Frazão, E.P., Santos, C.L.A. and Plácido, J.S. Jr (2010) Characterization of the Brazilian continental shelf adjacent to Rio Grande do Norte state, NE Brazil. Brazilian Journal of Oceanography 58, 4354.CrossRefGoogle Scholar
Vital, H., Stattegger, K., Amaro, V.E., Schwarzer, K., Frazão, E.P., Tabosa, W.F. and Silveira, I.M. (2008) A modern high-energy siliciclastic-carbonate platform. In Hampson, G., Steel, R., Burgess, P. and Dalrymple, (eds) Recent advances in models of siliciclastic shallow-marine stratigraphy. Tulsa, OK: SEPM Society for Sedimentary Geology, pp. 177190.CrossRefGoogle Scholar
Vroom, P., Smith, C., Coyer, J., Walters, L., Hunter, C., Beach, K. and Smith, J. (2003) Field biology of Halimeda tuna (Bryopsidales, Chlorophyta) across a depth gradient: comparative growth, survivorship, recruitment and reproduction. Hydrobiologia 501, 149166.CrossRefGoogle Scholar
Walters, L.J., Coyer, J.A., Hunter, C.L., Beach, K.S. and Vroom, P.S. (2002) Assexual propagation in the coral reef macroalga Halimeda (Chlorophyta, Bryopsidales): production, dispersal and attachment of small fragments. Journal of Experimental Marine Biology and Ecology 278, 4765.CrossRefGoogle Scholar
Walters, L.J. and Smith, C.M. (1994) Rapid rhizoid production in Halimeda discoidea decaisne (Chlorophyta, Caulerpales) fragments: a mechanism for survival after separation from adult thalli. Journal of Experimental Marine Biology and Ecology 175, 105120.CrossRefGoogle Scholar
Weiner, J., Kinsman, S. and Williams, S. (1998) Modeling the growth of individuals in plant populations: local density variation in a stand population of Xanthium strumarium (Asteraceae). American Journal of Botany 85, 16381645.CrossRefGoogle Scholar
Yñiguez, A.T., McManus, J.W. and Deangelis, D.L. (2008) Allowing macroalgae growth forms to emerge: use of an agent-based model to understand the growth and spread of macroalgae in Florida coral reefs, with emphasis on Halimeda tuna . Ecological Modelling 216, 6074.CrossRefGoogle Scholar
Yong, Y.S., Yong, W.T.L. and Anton, A. (2013) Analysis of formulae for determination of seaweed growth rate. Journal of Applied Phycology 25, 18311834.CrossRefGoogle Scholar
Figure 0

Fig. 1. Location of the study area. Dark grey region in the detail indicates the area of occurrence of Halimeda opuntia.

Figure 1

Fig. 2. Area of occurrence of Halimeda opuntia. The black circles illustrate the average area of coverage at each point sampled indicating no differences between transects or reef areas.

Figure 2

Table 1. Summary statistics (mean ± SEM) of Halimeda opuntia J.V. Lamouroux on Brazilian sandstone reefs.

Figure 3

Fig. 3. Area covered by Halimeda opuntia in each month of the study.

Figure 4

Fig. 4. Temporal variation in standing stock, coverage and biomass of Halimeda opuntia on Brazilian reefs.

Figure 5

Fig. 5. Temporal variation of average growth rate and daily net growth of Halimeda opuntia. Daily net growth represents how much algal mass was lost per day after adjusting for the respective area variations.

Figure 6

Fig. 6. Average of six environmental variables in the 30 days prior to the campaigns to study Halimeda opuntia on Brazilian reefs.

Figure 7

Fig. 7. Scatterplot showing a negative relationship between Halimeda opuntia standing stock and wind speed (r = −0.88, P = 0.02), but a positive relation with rainfall (r = 0.88, P = 0.02).

Figure 8

Fig. 8. Plot of four environmental variables on the space formed by the first and second Principal Components Analysis. PAR, photosynthetically active radiation.