- Research article
- Open Access
The biomass–density relationship in seagrasses and its use as an ecological indicator
BMC Ecologyvolume 18, Article number: 44 (2018)
Biomass–density relations have been at the centre of a search for an index which describes the health of seagrass meadows. However, this search has been complicated by the intricacy of seagrass demographics and their complex biomass–density relations, a consequence mainly of their modular growth and clonality. Concomitantly, biomass–density upper boundaries have been determined for terrestrial plants and algae, reflecting their asymptotic maximum efficiencies of space occupation. Each stand’s distance to its respective biomass–density upper boundary reflects its effective efficiency in packing biomass, which has proved a reliable ecological indicator in order to discriminate between taxonomic groups, functional groups and clonal vs. non-clonal growth.
We gathered data from 32 studies on 10 seagrass species distributed worldwide and demonstrated that seagrasses are limited by their own boundary line, placed below the boundaries previously determined for algae and terrestrial plants. Then, we applied a new metric—dgrass: each stand’s perpendicular distance to the seagrass boundary—and used this parameter to review fundamental aspects such as clonal growth patterns, depth distribution, seasonality, interspecific competition, and the effects of light, temperature and nutrients.
Seagrasses occupy space less efficiently than algae and terrestrial plants. Using only their biomass and density data we established a new and efficient tool to describe space occupation by seagrasses. This was used with success to evaluate their meadows as an ecological indicator for the health of coastal ecosystems.
When individuals in an even-aged monospecific plant stand undergo active growth increasing their biomass, competitive stress may induce mortality. The consequent elimination of the weaker plants releases resources (space, light and nutrients) facilitating the further growth of survivors. This dynamic, commonly known as self-thinning, also reflects the efficiency of space occupation as more efficient stands (or species) exhibit higher biomasses under similar stand densities (numbers of individuals per unit area). The first studies of self-thinning, dating from the 1950s onwards [1,2,3,4], established a relationship between density (D) and mean plant mass (w) given by w = kD−3/2 or equivalently log10w = log10k − 1.5 log10D. Here, w refers to above-ground biomass and k is an allometric constant. In addition to many reported self-thinning slopes in plants ranging from small herbs to trees being close to -3/2, the law has also been applied to some mixed species stands  and across species regressions of different sized plants [6, 7]. Later the relationship evolved into an equivalent derived from stand biomass per unit area (B) and density: B = kD−1/2 or equivalently log10B = log10k − 0.5log10D . This new relation solved two problems: (i) auto-correlation, as the former w-D relation required the number of individuals to estimate the quantities on both sides of the equation, and (ii) mean biomass increasing without actual growth but just because smaller individuals died [8, 9]. This was the first improvement aiming at establishing a set of best-practices for the assessment of biomass–density relations and involving aspects such as the data quality and choice of regression methods [9, 10]. More recently, the numerical aspects were further improved with the development and application of model 2 quantile regression and careful data screening [11, 12]. These adjustments brought about the realization that the slope of the biomass–density regression should be reset at values closer to − 0.33 [10, 13]. Despite these advances, substantial debate followed as to what the self-thinning “law” is as well as controversy as to whether a law really existed at all [10, 14, 15]. Henceforth, biomass–density relationships were categorized into three different aspects (Fig. 1):
The (intraspecific) dynamic self-thinning line is the straight line that is approached, then followed by the time trajectory of a crowded monospecific stand as it grows [9, 10, 13,14,15]. Stands may have different dynamic thinning lines depending on the environmental conditions (= carrying capacity). Changes in the slope and intercept of this dynamic thinning line for stands of the same species usually relate to the allometry of a plant species  as well as to resource or temperature limitation [17,18,19]. In the later cases, flatter slopes associated with lower intercepts reflected smaller carrying capacities of the environment.
The species (or higher ranked taxon) boundary line—the upper boundary of possible biomass–density combinations for a given taxon from the plant kingdom . This line is fit to the most extreme density-biomass combinations of hundreds of stands of the same taxon. Theoretically, the y-intercept of a taxon provides information about its maximum capacity to pack biomass above-ground.
The interspecific biomass–density relationship (IBDR ) and its static upper boundary characterizing the maximum biomass–density limit for all the species which make up the plant kingdom. Weller  analysed plant data setting the boundary at log10B = 3.91 − 0.33 log 10D. However, this boundary was estimated from Ordinary Least Squares (OLS) and thus, although applied to the most extreme stands, it still determined a central tendency. Consequently, Scrosati  re-analysed the data from Weller  setting the plant boundary at log10B = 4.87 − 0.33 log10D and dubbing it ultimate biomass–density Line (UBDL) . More recently, Creed et al.  termed this boundary the interspecific boundary line (IBL).
Since the 1980s phycologists have tested the biomass–density relations in seaweeds, generally obtaining results similar to those obtained from land plants [22,23,24,25,26,27]. However, a recent analysis of data revealed that algae occupy space more efficiently . With a higher intercept and steeper slope (log10B = 6.69 − 0.67 log10D), the algal IBL was placed significantly higher than the plant IBL (Fig. 1), thus representing the ultimate boundary known for any life-type on the planet. The higher efficiency of algae was considered to be a consequence of the aquatic environment which facilitates the algal morphology and ecophysiology relative to that of plants. Compared to land plants, algae are neither limited by water availability nor bare costs of transporting water for transpiration; they also acquire nutrients and photosynthesize through the whole of their surfaces, and do not require (non-photosynthetic) deep roots or rigid tissues for support as buoyancy keeps them upright.
Seagrasses grow by the iteration of modules which may remain physiologically integrated or separate to form clones. Previously it has been shown in plants and algae that the intraspecific dynamic biomass–density relationship (self-thinning) often does not apply to clonal species. These may have clonal integration because modules (ramets) belong to the same genetic individual (genet) and are physically interconnected, which may allow the sharing of acquired resources between ramets and offset competition . Consequently, an increase of stand biomass may arise from the increase in plant (ramet) sizes, increase in plant (ramet) density, or both. Hutchings  identified several non-thinning biomass–density dynamics typical of clonal plants. On the other hand, Westoby  and de Kroon and Kalliola  identified clonal plants where the specificities of their life-histories may (or not) lead them to self-thin, depending on additional factors. The occurrence of self-thinning is also variable in clonal algae and seems to depend on their specific life-history, morphological characteristics and habitat [20, 31,32,33,34]. Although not necessarily self-thinning, terrestrial clonal plants and clonal algae were nevertheless demonstrated to be limited by an IBL [12, 28]. Thus, it was still possible to use their stand’s distance to their IBL as an estimator of their efficiency of space occupation. Having proved applicable to all sorts of algae, both clonal and non-clonal, this new method provided exciting new evidence about algal growth and ecology . In this work we tested whether the same applies to seagrasses. Taxonomically, seagrasses are flowering plants (i.e., angiosperms), sharing with their terrestrial counterparts the features that make angiosperms the most complex taxa among autotrophic life-forms: as well as having flowers they also possess rigid tissues, extensive roots and a primary circulatory system carrying water and nutrients. On the other hand, they spend a significant part (if not all) of their life underwater, where they can also benefit from the water environment, as do macroalgae. Furthermore, in contrast to terrestrial angiosperms, seagrasses acquire nutrients both through their leaves and roots, a capacity that can even become essential for their metabolic rates . Given this seagrass duality, the first question arising is whether seagrasses pack above-ground biomass similarly to algae and/or terrestrial plants?
The dynamic self-thinning line is thought to depend intrinsically on the plant taxon being considered, besides being dependent extrinsically on the environmental carrying-capacity (proposed by Weller [9, 13, 16]). Biologically more complex plants, appearing later in Life’s evolutionary history and generally comprehending larger individuals (namely trees), showed significantly steeper slopes along with lower intercepts (Fig. 1). Backed by elementary biostatistics, ecological theory proposed the intercept as a standard to estimate the biomass carrying capacity under equivalent densities. Hence, Weller’s results suggested that simpler plants occupy space (i.e., pack above-ground biomass) more efficiently. Weller’s  finding was refuted by Lonsdale  on alleged methodological grounds. Recently, new evidence emerged corroborating Weller’s findings. Creed et al.  also found that the algal efficiency of space occupation varies with algal traits, namely with taxonomic group, functional group and clonality. And similarly to Weller’s analysis simpler taxa occupy space more efficiently. Concurrently, the ability of plant communities to pack biomass into the volume effectively exploited has been demonstrated to depend on biodiversity and on the efficiency of resource exploitation . Creed et al.  obtained their results using a more straightforward and robust methodology than the one used by Weller  with each stand being measured by its perpendicular distance to the algal IBL (dalgal). Following the results obtained by Weller , Proulx et al.  and Creed et al. , in this study we use the perpendicular distance to the seagrass IBL (dgrass) to test whether different species show different efficiencies of space occupation in seagrasses too.
As seagrasses occur at the interface of land and sea they suffer diverse stressors, one of which is eutrophication . Increasing nutrient loadings are potentially toxic and also promote blooms of opportunistic macrophytes . Both effects have a negative impact on biomass or density in seagrass stands that may extend to the whole ecosystem [37,38,39,40,41,42]. Seagrasses more impacted by algal competitors tended to be smaller and/or less abundant (see Fig. 2 in Thomsen et al. ). As seagrasses are foundation species the dynamics of their stands has been used as a proxy for ecosystem health [39, 43,44,45,46,47]. The recent methodology by Creed et al. , which congregates the biomass and density data into a single metric irrespective of their specific correlation, estimates the stands’ efficiency of space occupation. We test here whether this metric applied to seagrasses is a meaningful ecological indicator. Do healthier seagrass stands better optimize their above-ground space occupation by falling closer to a biomass–density boundary line?
Biomass–density data retrieval
We gathered data of shoot density and above-ground biomass from 32 studies on 10 seagrass species distributed worldwide (Table 1). Most of the Halodule wrightii data was provided by one author (J.C.C.). The data from Plus et al.  was provided by Dr Martin Plus. The remaining data was retrieved from the respective publications using appropriate software. The compilation of data was carried out during years 2017 and 2018, and used the Google search engine as well as the search engines in the webpages of all cited publications. The search keywords included ‘biomass’, ‘density’, ‘seagrass’ and the species scientific denominations. We also searched the publication listings of the most cited authors in the subject and the reference lists of the cited works.
Biomass–density data analysis
From the 2954 biomass-shoot density observations available we selected the 500 closer to the top right corner of the seagrass biomass–density plot. With these 500 observations we estimated a seagrass-specific IBL following the procedure described by Creed et al.  and choosing the 99.9% quantile for the linefit. Once this boundary was known, we estimated the stands’ perpendicular distances to this line (the dgrass presented in “The dgrass metric” section below) in order to test:
How seagrasses worldwide have different efficiencies of space occupation. However, the data comprised species with widely different numbers of observations, which would inevitably bias comparisons based on properties of the samples’ distributions. In an attempt to overcome this problem and simultaneously compare among seagrasses at their maximum efficiency, for this task we selected from each taxon the five observations with the smallest dgrass;
How Z. marina and Z. japonica in Dadae Bay show different efficiencies of space occupation while inhabiting the same shoreline;
How the efficiency of space occupation varies seasonally and spatially;
How the efficiency of space occupation varies with light, temperature and nutrient concentrations.
The dgrass metric
The estimation of each stand’s perpendicular distance to the seagrass IBL (dgrass) requires the linear coefficients of this boundary line. With the general IBL equation corresponding to log10B = β0 + β1 log10D, the coefficients for the seagrass IBL were β0 = 4.569 and β1 = − 0.438 (as presented in “Results” and “Discussion” sections). The angle θ between the dgrass vector (oblique in the log10B-to-log10D orthogonal plane) and the log10B vertical axis is the same angle between the seagrass IBL and the log10D horizontal axis (Fig. 2). Hence, θ = arctg(|β1|), which in the seagrass case corresponds to θ = 0.413. The perpendicular distance is the cosine of θ, in the case of seagrass cosθ = 0.916, multiplied by the B vertical distance: dgrass = (log10Ḃ − log10B)∙cosθ. The vertical distances require the observed log10B and the estimated log10Ḃ = β0 + β1∙log10D.
The biomass–density data was correlated with additional biotic and abiotic data to test which factors may determine seagrasses demography. These additional data were taken from the studies from which the biomass–density data was taken (Table 1). The full dataset is provided as additional files.
Lee et al.  compared the demography of the Zostera marina and Zostera japonica native to the South Korean coast. The Zostera marina occupied the subtidal whereas the Zostera japonica occupied the intertidal. Sampling occurred from July 2001 to July 2002, and retrieved a set of biotic variables related to the populational, morphometric and ecophysiological properties of the seagrasses, and another set of abiotic variables characterizing the water properties, sediment properties and light regime.
Plus et al.  studied factors determining the primary production of Zostera noltii in the Thao lagoon from July 1996 to August 1998. The biotic variables retrieved included biomass and density. The abiotic variables included water temperature, salinity, photosynthetic active radiation (PAR) at the surface and on the bottom, light extinction coefficient, and the concentrations of NH4+, NO2−, NO3− and PO43−.
Stands of Zostera noltii were followed in the Ria Formosa lagoon system in southern Portugal at four locations along a pollution gradient . Site 1 was the closest to the wastewater treatment plant (more polluted—270 m away). At the opposite extreme, site 4 was the farthest from the wastewater treatment plant (little polluted—1500 m away and located in the main channel). All sites were simultaneously assessed during July 2001 (summer), November 2001 (autumn), February 2002 (winter) and May 2002 (spring). To this data we added data from the Z. noltii studied by Plus et al.  in the Thau Lagoon.
Aiming for an indicator of water quality and ecosystem health, Garcia-Marín et al.  developed the ZoNI index from compiled data for intertidal Z. noltii stands from the southern Iberian Peninsula. Five stands were sampled in the summer of 2010 in Ria Formosa, south Portugal, three of them assumed under normal conditions (R1, R2 and R3) and two of them impacted by extreme pollution from urban effluents (I1 and I2). Three stands were sampled both in Huelva (H1, H2 and H3) and in Cádiz (C1, C2 and C3), south Spain, during the summers of 2009–2011. Biotic variables included biomass, density, and other populational, morphometric and biochemical characteristics. Abiotic variables included ammonium, nitrate and phosphate concentrations, and the first principal component (PC1) extracted from a principal components analysis (PCA) synthesized the nutrients general dynamics.
Similarly to algae and terrestrial plants, seagrasses also had their specific biomass–density boundary (Fig. 3). With a slope of –0.438 and an intercept of 4.569, the seagrass boundary (IBL) placed far below that of the plant and algae boundaries. This lower placement was mainly due to the lower B-axis intercept. Below their own IBL, seagrasses often showed a positive correlation between their biomass and shoot density, while differentiating among themselves by occupying different bands of the biomass-shoot density spectra (Figs. 3 and 4). The five best efficiencies observed for each species showed that different seagrasses have different maximum efficiencies of space occupation (Fig. 5a). Permutation tests (a class of non-parametric ANOVA where the null hypothesis is simulated by randomly redistributing the data) with 10,000 replications showed that these maximum efficiencies were significantly different (p = 0.0001, d.f.error = 40, d.f.groups = 7).
Stands of Z. marina and Z. japonica at Dadae Bay showed similar biomasses, though Z. japonica achieved this under much higher densities of smaller shoots (shown in Fig. 3), resulting in a much better space occupation efficiency (Fig. 5b). The Dadae Bay situation seems to be case-specific, suggesting the occurrence of something particular, as the worldwide trend is for both species to show similar efficiencies of space occupation (Fig. 5a).
The efficiency of space occupation by Zostera marina in the Baltic and by Posidonia sinuosa in Western Australia decreased with depth (Fig. 6). For the P. sinuosa in Western Australia a seasonal trend was also observed where the efficiency of space occupation increased in the winter and decreased in the summer (Fig. 6). This seasonal pattern opposes the pattern observed for the other species and is presented as follows.
The dgrass metric confirmed the general seagrass seasonality but also highlighted regional differences (Fig. 7). The seagrasses from higher latitudes within the northern hemisphere, namely Cymodocea nodosa and the Zostera spp., peaked their efficiency of space occupation during the spring–summer and reached the low point during the winter. Located closer to the equator, the Halodule wrightii stands in the northern hemisphere peaked during June–September (summer). One of the stands located in the southern hemisphere peaked during November-January (summer) whereas the other slightly delayed its peak to the summer-autumn. Concomitantly with the seasonal dynamics, the dgrass also discriminated between stands at different locations experiencing different environments. This was particularly evident in the C. nodosa and the H. wrightii populations.
The effects of ammonium and phosphate on the efficiency of space occupation peaked at intermediate concentrations, mimicking the dynamics generally reported for the effects of these nutrients on the growth rates of seagrasses (Fig. 8). The optimal nutrient concentrations for Z. noltii in Ria Formosa were much higher than in the Thau Lagoon. The Z. noltii’s most efficient space occupation at Ria Formosa were found at the sampling sites located at intermediate distances from the wastewater treatment plant. The concentrations of other nutrients showed no effects over the Z. noltii stands in the Thau Lagoon and in Ria Formosa (Fig. 8).
In the Thau lagoon, temperature and PAR at the bottom (i.e., at stand level) were the fundamental factors governing Z. noltii’s efficiency of space occupation. The dgrass was better fit to the congregation of a quadratic function for temperature and another quadratic function for PAR (Eq. 1), although the resulting forth order term was disregarded from the numerical estimation of coefficients.
The efficiency of space occupation (i.e., dgrass) mimicked the photosynthetic and growth responses to irradiance as modelled by the classical photosynthesis-irradiance (P–I) curves. Hence, given a fixed temperature, there was an optimal PAR above and below which the space occupation deteriorated i.e., dgrass increased (Fig. 9). Concomitantly, given a fixed PAR, there was an optimal temperature above and below which the space occupation deteriorated (Fig. 9). The optimal temperature and optimal PAR were negatively correlated—i.e., as temperature increased its correspondent optimal PAR decreased—suggesting that temperature had a positive effect on the photosynthetic efficiency.
The ZoNI index correlated remarkably well with the dgrass (r = 0.92, see Additional file 1: Fig. S1). Correlating the dgrass with the abiotic data from Garcia-Marin et al.  (Fig. 10), it was found that Huelva’s stand 1 (H1) during 2011 showed an extremely weak space occupation as a consequence of some unaccounted factor, and hence was disregarded. The analysis on the remaining stands revealed that the fundamental constrain to Z. noltii space occupation was phosphate pollution above 0.5 μM. The polluted stands (I1 and I2) correlated phosphate concentrations above this threshold with an inefficient occupation of space (i.e., very large dgrass). All other stands were below this phosphate concentration threshold and showed much smaller and similar dgrass. Ammonium and nitrate had no influence on the inefficient space occupation by the polluted sites. In fact, these sites were even within the group of stands with lower nitrate concentrations. The correlation between dgrass and the nutrients PC1 was largely contaminated by the noise introduced in the nutrients PC1 by the presence of the ammonium and nitrate variables.
The lower placement of the seagrass IBL—due to its lower B-axis intercept—shows that seagrasses are poor occupiers of space, both compared to terrestrial plants, with whom they share taxonomic similarities, and when compared to algae, with whom they share the marine environment. The classical biomass–density models [1, 13] determine that one of the causes for a lower intercept is less volume exploited per unit plant stand surface (the other being less biomass per used volume). There are three possible causes for this event:
Wasted available volume due to short shoot heights, and thus also smaller height-to-width ratio. This constitutes vertical waste;
Wasted available volume due to internode lengths larger than shoot widths. This constitutes horizontal waste that is fixed along the vertical axis;
The scarcity of data with the required metrics did not allow generalization as to which aspects of seagrass allometry mentioned above affect their efficiency of space occupation. Nevertheless, the hypothesis (i) seems the most likely as seagrass canopies are relatively short when compared to their terrestrial counterparts. In permanently immersed stands, the shoots usually do not extend to fully occupy the available height (to sea-surface). Hence, it is unlikely that it is the sea-surface that usually acts as a barrier for vertical growth of seagrasses. A possible explanation for the height limitation comes from the fact that seagrasses generally lack secondary metabolites with anti-fouling properties frequently found in algae and hence need to avoid build-up of epiphytes which reduce light and nutrient uptake. They do so by using a basal meristem and changing leaves regularly, which puts a limit on both leaf length and biomass accumulation. The hypothesis (ii) about shoots narrower than internode length is unlikely, particularly if the stands’ have intermingled genets and/or branched rhizomes. Relative to hypothesis (iii) land crown shapes are less efficient at occupying space due to the structural adaptations of terrestrial plants for gravity. In water, less constrained by gravity, algae can use more efficient shapes, ultimately leading to an algae IBL placed above that of plants . Immersed seagrasses are similarly less constrained by gravity, but shoot shape is constrained by leaf number per shoot and leaf width, which are relatively non-plastic seagrass variables, so space occupation depends on leaves elongating and presenting themselves horizontally. It should also be noted that seagrasses may not benefit as much from the support offered by the aqueous medium as seaweeds do because they still present some limiting structural tissue of their terrestrial ancestors (see Creed et al.  on this topic). Compared to their terrestrial counterparts, in these marine plants there may be a trade-off between the benefit of reduced investment in structural tissue and disadvantage of increased light attenuation. Despite their inefficient space occupation, seagrasses can dominate the intertidal and shallow subtidal in estuarine and lagoon systems, leading to the conclusion that the efficiency of space occupation is not the characteristic that concedes their ecological advantage over their algal competitors. As most algae need hard substrate to attach to, the fact that seagrasses can colonize soft sandy bottoms is the most likely explanation for their ecological advantage.
Different seagrasses showed different maximum efficiencies of space occupation. However, these estimates were very sensitive to the placement of the seagrass IBL, which in turn was estimated based only on ten species for which data are available. Furthermore, in the case of Z. japonica it is possible that its maximum efficiency, estimated lower than many other species, results from insufficient sampling. On the contrary, H. wrightii was sampled at 25 stations scattered along both tropical bands in the northern and southern hemisphere totalling 1005 observations, whereas Z. marina was sampled at 23 stations scattered worldwide and totalling 908 observations. Hence, although we should be cautious analysing these results, it seems safe to say that Z. noltii, P. sinuosa, H. wrightii, and C. nodosa are able to occupy space more efficiently than T. testudinum, Z. japonica and Z. marina.
Different seagrasses occupy different bands of the biomass-shoot density spectra, suggesting conditional differentiation of co-occurring seagrass species. As an example, Z. marina co-occurs with Z. noltii in the Atlantic and with Z. japonica in the Pacific Ocean; where it occupies a biomass–density band conspicuously different from those of its competitors. Within their respective bands, seagrass species tend to show a positive biomass–density correlation contrasting with the negative correlation typical of non-clonal plants and algae. This positive correlation typical of seagrasses is a consequence of their clonality  and suggestive of the growth-form plasticity already demonstrated in clonal terrestrial plants [51,52,53] and clonal algae : the high shoot densities match the phalanx growth-form suited to dominate favourable environments . On the other hand, when the environment is not favourable, the low shoot densities match the guerrilla growth-form allowing for faster dispersion in the search for better locations . By mimicking the positive biomass–density correlation of seagrasses, the perpendicular distance to the seagrass IBL (i.e., the dgrass), besides being an index for the efficiency of space occupation, is also an index for the trade-off in clonal-growth-form plasticity: small dgrass corresponding to high efficiencies and phalanx growth, and large dgrass corresponding to low efficiencies and guerrilla growth. This dynamic may seem paradoxical at first sight: when the environment is sub-optimal, adopting the guerrilla growth form gets the stand further away from the efficient space occupation, and thus more susceptible to competitors. On the other hand, when the environment is favourable, the phalanx growth form allows for the stand to reach the maximum efficiency more quickly, but at the cost of decreasing the stand biomass potential. It intuitively seems advantageous when competing with other seagrasses but not when competing with algae that have much higher space occupation efficiencies (see Creed et al.  for the algal IBL). Regarding this unorthodox phalanx strategy, we highlight three aspects:
Seagrasses avoid competition from algae by colonizing a substrate to which most algae cannot attach and grow;
Thomsen et al.  found that small seagrasses (consequently, with lower biomasses and higher densities) are more vulnerable to their algae competitors;
The advantage of reaching the boundary more quickly may also be related to complementary aspects. For instance, reaching the maximum biomass quicker may release resources for other aspects of seagrass development such as sexual reproduction, rhizome storage or growth at the edges of the stand.
The Z. marina and Z. japonica stands at Dadae Bay  are a good example of how seagrasses make use of phenotypic plasticity to adapt to the abiotic and biotic conditions, and how this affects their efficiency of space occupation. Both species had equally long leaves but the Z. marina leaves were wider. Furthermore, Z. marina showed a higher leaf production rate, leading to a higher number of leaves per shoot . The consequent larger biomass per shoot of Z. marina relative to Z. japonica is a global attribute. However, at Dadae Bay it was fully compensated by the higher shoot density of Z. japonica, a consequence of its shorter internode length  and/or intermingled genets, and which led Z. japonica to occupy space more efficiently. The Z. japonica in Dadae Bay occupied the intertidal, which is its preferred environment [49, 55]. Hence, this stand adopted a phalanx growth-form whose high density and closeness to the IBL disables the colonization by its seagrass competitors. As for Z. marina, the application of the dgrass to the data taken from the studies by Krause-Jensen et al.  and Kim et al. , demonstrated that, worldwide, Z. marina is more efficient in the lower intertidal. However, in Dadae Bay the lower intertidal was already occupied by Z. japonica. Hence, Z. marina was relegated to the shallow subtidal where it adopted a guerrilla growth-form whose focus is on spreading to its preferred neighbouring lower intertidal. Contrasting the phalanx grow form of Z. japonica with the guerrilla grow form of Z. marina, the Z. marina internode length was on geometric average 2.5 times longer than that of Z. japonica while its reproductive effort was approximately 3 times higher (Additional file 2: data and Lee et al. ).
The efficiency of space occupation revealed itself to be an outstanding tool for understanding aspects of seagrass ecology. Based on this efficiency, we developed the dgrass metric, whose application to seagrasses worldwide enabled us to generalize regarding seasonality, their depth profiles, and their response to nutrients, light and temperature. We have already done the same for algae, developing the dalgal metric and applying it with success . In the case of seagrasses, the dgrass seasonality matched the generalized summer peaks of seagrass growth and Photosynthesis–Irradiance (P–I) curve parameters . The effects of ammonium and phosphate concentrations on dgrass mimicked the dynamics generally reported for the effects of these nutrients over the growth rates: very low concentrations were detrimental as autotrophs need nutrients to survive and grow. However, at high concentrations these molecules became toxic, shifting their role from nutrients to pollutants [39,40,41, 43,44,45,46]. The optimal nutrient concentrations for Z. noltii in Ria Formosa higher than in the Thau Lagoon suggested that some additional factor was influencing the optimal points. The Z. noltii’s most efficient space occupations in Ria Formosa located at intermediate distances from the wastewater treatment plant revealed that the wastewater often had a beneficial effect over those patches by raising the ammonium and phosphate concentrations to optimal values. The Z. noltii’s dgrass response to irradiance mimicked a P–I curve with photo-inhibition above the optimal irradiance, as has been reported for seagrasses  and for Z. noltii in particular . Concomitantly, the Z. noltii’s dgrass response to temperature mimicked the previously reported parabola-type relation between seagrass growth and temperature .
We take the remarkable correlation of the ZoNI index with the dgrass as evidence that both are good indicators of the health of seagrass meadows. Monitoring seagrass health has been proposed for coastal water quality assessment [43,44,45,46,47]. The development of an ecological indicator relying exclusively on the biomass–density relation of seagrasses has been previously attempted, unsuccessfully [47, 59]. At that time, the existence of a seagrass IBL and the fact that the efficiency of space occupation is a fundamental aspect of seagrass ecology were unknown. Hence, back then it was impossible for the respective authors to succeed in their objective. Nevertheless, our new findings demonstrate that this is a useful concept and tool, which we put back on track reformulated around the efficiency of space occupation. We present three strong arguments for its future application:
Its simplicity The dgrass only requires two variables—the stand’s above-ground biomass and shoot density—which are also the easiest to retrieve and describe well habitat availability and complexity;
Its generality It is applicable to any seagrass (estimating dgrass), any terrestrial plant (estimating dplant) and any algae (estimating dalgal), on any location on the planet, and even enables comparisons within and among these three groups. As examples, while the ZoNI is restricted to Z. noltii  and the POMI to P. oceanica , the dgrass allows comparisons among all seagrasses. Furthermore, it is possible to compare among seagrasses, algae and plants, estimating their distances to their respective IBLs, namely the dgrass, dalgal and dplant;
Its background The dgrass is grounded on the biomass–density relations of plants and algae, thoroughly studied since the 1950s with the self-thinning law accolade of “the only generalization worthy of the name law in plant ecology” . It has already performed well when applied to other subjects. Its previous application to worldwide algal data has determined that the efficiency of space occupation in the algae depends on taxonomic group, functional group, clonality and latitude . Its current application to seagrass worldwide has unveiled how it can be used to compare stands and interpret efficiency in light of intrinsic (population dynamics) and extrinsic (eco-physiological) factors.
The placement of seagrass meadows in the biomass–density plot is limited by their interspecific boundary line (IBL) setting a maximum efficiency of space occupation. Furthermore, species tend to differentiate the bands each occupies in this scatter-plot, which is evidence of their conditional differentiation. The efficiency of space occupation by seagrasses, requiring only the biomass and shoot-density of their stands and measured by their perpendicular distance to the seagrass IBL, revealed a highly useful indicator of their ecological condition. It identified the summer as their most favourable season and the lower intertidal as their preferred depth. It discriminated among locations. It identified which nutrients were in excess, thus acting as pollutants, and beyond which concentrations did these have a deleterious effect over seagrasses. The dgrass was revealed to be a most efficient ecological index with general application and comparable with similar indexes developed for seaweeds and terrestrial plants.
- dalgae :
perpendicular distance to the algae IBL
- dgrass :
perpendicular distance to the seagrass IBL
- Dplant :
perpendicular distance to the plant IBL
interspecific biomass–density relationship
interspecific boundary line
photosynthetic active radiation
principal components analysis
principal component #1
Posidonia oceanica multivariate index
ultimate biomass–density Line
individual plant weight
Zostera noltii multi-metric index
Yoda K, Kira T, Ogawa H, Hozumi K. Self-thinning in overcrowded pure stands under cultivated and natural conditions (Intraspecific competition among higher plants. J Biol Osaka City Univ. 1963;14:107–29.
Hutchings MJ. Ecology’s law in search of a theory. New Scientist. 1983;98:765–7.
Cousens R, Hutchings MJ. The relationship between density and mean frond weight in monospecific seaweed stands. Nature. 1983;301:240–1.
Lonsdale WM, Watkinson AR. Plant geometry and self-thinning. J Ecol. 1983;71(285–29):7.
White J, Harper JL. Correlated change in plant size and number in plant populations. J Ecol. 1970;58:467–85.
Gorham E. Shoot height, weight and standing crop in relation to density of monospecific plant stands. Nature. 1979;279:148–50.
Westoby M, Howell J. Influence of population structure on self-thinning of plant populations. J Ecol. 1986;74:343–59.
Westoby M. The place of the self-thinning rule in population dynamics. Am Nat. 1981;118:581–7.
Weller DE. A reevaluation of the -3/2 power rule of plant self-thinning. Ecol Monogr. 1987;57:23–43.
Lonsdale WM. The self-thinning rule: dead or alive? Ecology. 1990;71(4):1373–88.
Zhang L, Bi H, Gove JH, Heath LS. A comparison of alternative methods for estimating the self-thinning boundary line. Can J For Res. 2005;35:1507–14.
Creed J, Norton TA, Caetano D, Vieira VMNCS. Seaweeds surpass plants setting life-on-Earth’s limit for biomass concentration. BMC Ecol (submitted).
Weller DE. The interspecific size-density relationship among crowded plant stands and its implications for the -3/2 power rule of self-thinning. Am Nat. 1989;133:20–41.
Osawa A, Sugita S. The self-thinning rule: another interpretation of Weller’s results. Ecology. 1989;70:279–83.
Weller DE. The self-thinning rule: dead or unsupported? A reply to Lonsdale. Ecology. 1991;72:747–50.
Weller DE. Self-thinning exponent correlated with allometric measures of plant geometry. Ecology. 1987;68:813–21.
Morris EC, Myerscough PJ. Self-thinning and competition intensity over a gradient of nutrient availability. J Ecol. 1991;79:903–23.
Morris EC. Effect of localized placement of nutrients on root competition in self-thinning populations. Ann Bot. 1996;78:353–64.
Morris EC. How does fertility of the substrate affect intraspecific competition? Evidence and synthesis from self-thinning. Ecol Res. 2003;18(3):287–305.
Scrosati RA. The interspecific biomass–density relationship for terrestrial plants: where do clonal red seaweeds stand and why? Ecol Lett. 2000;3:191–7.
Scrosati R. Review of studies on biomass–density relationships (including self-thinning lines) in seaweeds: main contributions and persisting misconceptions. Phycol Res. 2005;53:224–33.
Cousens R, Hutchings MJ. The relationship between density and mean frond weight in monospecific seaweed stands. Nature. 1983;301:240–1.
Creed JC. Spatial dynamics of a Himanthalia elongate (Fucales, Phaeophyta) population. J Phycol. 1995;31:851–9.
Flores-Moya A, Fernández JA, Niell FX. Growth pattern, reproduction, and self-thinning in seaweeds: a re-evaluation in reply to Scrosati. J Phycol. 1997;33:1080–1.
Creed JC, Kain JM, Norton TA. An experimental evaluation of density and plant size in two large brown seaweeds. J Phycol. 1998;34:39–52.
Arenas F, Fernández C. Size structure and dynamics in a population of Sargassum muticum (Phaeophyceae). J Phycol. 2000;36:1012–20.
Steen H, Scrosati R. Intraspecific competition in Fucus serratus and F. evanescens (Phaeophyceae: Fucales) germlings: effects of settlement density, nutrient concentration, and temperature. Mar Biol. 2004;144:61–70.
Hutchings MJ. Weight–density relationships in ramet populations of clonal perennial herbs, with special reference to the 3/2 power law. J Ecol. 1979;67:21–33.
Westoby M. The self-thinning rule. Adv Ecol Res. 1984;14:167–225.
de Kroon H, Kalliola R. Shoot dynamics of the giant grass Gynerium sagittatum in Peruvian Amazon floodplains, a clonal plant that does show self-thinning. Oecologia. 1995;101:124–31.
Lazo ML, Chapman ARO. Components of crowding in a modular seaweed: sorting through the contradictions. Mar Ecol Prog Ser. 1998;174:257–67.
Scrosati R, Servière-Zaragoza E. Ramet dynamics for the clonal seaweed Pterocladiella capillacea (Rhodophyta): a comparison with Chondrus crispus and with Mazzaella cornucopiae (Gigartinales). J Phycol. 2000;36:1061–8. https://doi.org/10.1046/j.1529-8817.2000.00041.x.
Rivera M, Scrosati R. Self-thinning and size inequality dynamics in a clonal seaweed (Sargassum lapazeanum, Phaeophyceae). J Phycol. 2008;44:45–9.
Scrosati R. Crowding in clonal seaweeds: does self-thinning occur in Mastocarpus papillatus shortly before stand biomass peaks? Aquat Bot. 2006;84:233–8.
Stapel J, Aarts TL, van Duynhoven BHM, de Groot JD, van den Hoogen PHW, Hemminga MA. Nutrient uptake by leaves and roots of the seagrass Thalassia hemprichii in the Spermonde Archipelago, Indonesia. Mar Ecol Prog Ser. 1996;134:195–206.
Proulx R, Rheault G, Bonin L, Roca IT, Martin CA, Desrochers L, Seiferling I. How much biomass do plant communities pack per unit volume? PeerJ. 2015;3:e849.
Burkholder JM, Tomasko DA, Touchette BW. Seagrasses and eutrophication. J Exp Mar Biol Ecol. 2007;350:46–72.
Valiela I, McClelland J, Hauxwell J, Behr PJ, Hersh D, Foreman K. Macroalgal blooms in shallow estuaries: controls and ecophysiological and ecosystem consequences. Limnol Oceanogr. 1997;42:1105–18.
Lee K-S, Park SR, Kim YK. Effects of irradiance, temperature, and nutrients on growth dynamics of seagrasses: a review. J Exp Mar Biol Ecol. 2007;350:144–75.
Thomsen MS, Wernberg T, Engelen AH, Tuya F, Vanderklift MA, Holmer M, et al. A meta-analysis of seaweed impacts on seagrasses: generalities and knowledge gaps. PLoS ONE. 2012;7(1):e28595.
Brun FG, Olivé I, Malta EJ, Vergara JJ, Hernández I, Pérez-Lloréns J. Increased vulnerability of Zostera noltii to stress caused by low light and elevated ammonium levels under phosphate deficiency. Mar Ecol Prog Ser. 2008;365:67–75.
Pergent G, Boudouresque C-F, Dumay O, Pergent-Martini C, Wyllie-Echeverria S. Competition between the invasive macrophyte Caulerpa taxifolia and the seagrass Posidonia oceanica: contrasting strategies. BMC Ecol. 2008;8:20. https://doi.org/10.1186/1472-6785-8-20.
Cabaço S, Machás R, Santos R. Biomass–density relationships of the seagrass Zostera noltii: a tool for monitoring anthropogenic nutrient disturbance. Estuar Coast Shelf Sci. 2007;2007(74):557–64.
Romero J, Martínez-Crego B, Alcoverro T, Pérez M. A multivariate index based on the seagrass Posidonia oceanica (POMI) to assess ecological status of coastal waters under the water framework directive (WFD). Mar Pollut Bull. 2007;55:196–204.
Cabaço S, Machás R, Vieira V, Santos R. Impacts of urban wastewater discharge on seagrass meadows (Zostera noltii). Estuar Coast Shelf Sci. 2008;78:1–13.
García-Marín P, Cabaço S, Hernández I, Vergara JJ, Silva J, Santos R. Multi-metric index based on the seagrass Zostera noltii (ZoNI) for ecological quality assessment of coastal and estuarine systems in SW Iberian Peninsula. Mar Pollut Bull. 2013;68:46–54.
Cabaço S, Apostolaki ET, Garcıa-Marín P, Gruber R, Hernandez I, Martınez-Crego B, et al. Effects of nutrient enrichment on seagrass population dynamics: evidence and synthesis from the biomass–density relationships. J Ecol. 2013;101:1552–62.
Plus M, Deslous-Paoli J-M, Dagault F. Factors influencing primary production of seagrass beds (Zostera noltii Hornem.) in the Thau lagoon (French Mediterranean coast). J Exp Mar Biol Ecol. 2001;259:63–84.
Lee SY, Kim JB, Lee SM. Temporal dynamics of subtidal Zostera marina and intertidal Zostera japonica on the southern coast of Korea. Mar Ecol. 2006;27:133–44.
Purves DW, Lichstein JW, Pacala SW. Crown plasticity and competition for canopy space: a new spatially implicit model parameterized for 250 North American tree species. PLoS ONE. 2007;2(9):e870. https://doi.org/10.1371/journal.pone.0000870.
Doust LL. Population dynamics and local specialization in a clonal perennial (Rannunculus repens): 1. The dynamics of ramets in contrasting habitats. J Ecol. 1981;69(3):743–55.
Chen X-S, Xie Y-H, Deng Z-M, Hou Z-Y. A change from phalanx to guerrilla growth form is an effective strategy to acclimate to sedimentation in a wetland sedge species Carex brevicuspis (Cyperaceae). Flora Morphol Distrib Funct Ecol Plants. 2011;206(4):347–50.
Ye X-H, Yu F-H, Dong M. A Trade-off between Guerrilla and Phalanx growth forms in Leymus secalinus under different nutrient supplies. Ann Bot. 2006;98(1):187–91.
Monro K, Poore AGB. The evolvability of growth form in a clonal seaweed. Evolution. 2009;63:3147–57. https://doi.org/10.1111/j.1558-5646.2009.00802.x.
Ruesink JL, Hong J-S, Wisehart L, Hacker SD, Dumbauld BR, Hessing-Lewis M, Trimble AC. Congener comparison of native (Zostera marina) and introduced (Z. japonica) eelgrass at multiple scales within a Pacific Northwest estuary. Biol Invasions. 2009. https://doi.org/10.1007/s10530-009-9588-z.
Krause-Jensen D, Middelboe AL, Sand-Jensen K, Christensen PB. Eelgrass, Zostera marina, growth along depth gradients: upper boundaries of the variation as a powerful predictive tool. Oikos. 2000;91:233–44.
Kim JB, Lee W-C, Lee K-S, Park J-I. Growth dynamics of eelgrass, Zostera marina, in the intertidal zone of seomjin estuary, Korea. Ocean Sci J. 2013;48(3):239–50.
Peralta G, Pérez-Lloréns JL, Hernández I, Vergara JJ. Effects of light availability on growth, architecture and nutrient content of the seagrass Zostera noltii Hornem. J Exp Mar Biol Ecol. 2002;269:9–26.
Vieira VM, Leitão F, et al. Mateus M.Biomass–density data analysis: a comment on Cabaço (2013). J Ecol. 2015;103:537–40.
Hall MO, Durako MJ, Fourqurean JW, Zieman JC. Decadal changes in seagrass distribution and abundance in Florida Bay. Estuaries. 1999;22(2B):445–59.
Möller T, Kotta J, Martin G. Spatiotemporal variability in the eelgrass Zostera marina L. in the north-eastern Baltic Sea: canopy structure and associated macrophyte and invertebrate communities. Estonian J Ecol. 2014;63(2):90–108.
Olesen B, Sand-Jensen K. Biomass–density patterns in the temperate seagrass Zostera marina. Mar Ecol Prog Ser. 1994;109:283–91.
Plus M, Chapelle A, Ménesguen A, Deslous-Paoli J-M, Auby I. Modelling seasonal dynamics of biomasses and nitrogen contents in a seagrass meadow (Zostera noltii Hornem.): application to the Thau lagoon (French Mediterranean coast). Ecol Model. 2003;161:213–38.
Tomasko DA, Hall MO. Productivity and Biomass of the Seagrass Thalassia testudinum along a Gradient of Freshwater Influence in Charlotte Harbor, Florida. Estuaries. 1999;22(3):592–602.
Gallegos ME, Merino M, Marbá N, Duarte CM. Biomass and dynamics of Thalassia testudinum in the Mexican Caribbean: elucidating rhizome growth. Mar Ecol Prog Ser. 1993;95:185–92.
Kaldy JE, Dunton KH. Above-and below-ground production, biomass and reproductive ecology of Thalassia testudinum (turtlegrass) in a subtropical coastal lagoon. Mar Ecol Prog Ser. 2000;193:271–83.
Enríquez S, Pantoya-Reyes NI. Form-function analysis of the effect of canopy morphology on leaf self-shading in the seagrass Thalassia testudinum. Oecologia. 2005;145:235–43.
Paynter CK, Cortés J, Engels M. Biomass, productivity and density of the seagrass Thalassia testudinum at three sites in Cahuita National Park, Costa Rica. Rev Biol Trop. 2001;49(2):265–72.
Medina-Gómez I, Madden CJ, Herrera-Silveira J, Kjerfve B. Response of Thalassia Testudinum morphometry and distribution to environmental drivers in a pristine tropical Lagoon. PLoS ONE. 2016;11(10):e0164014.
Larsson S. The production of the seagrass Thalassia hemprichii in relation to epiphytic biomass. Degree project for Master of Science in Biology, Department of Marine Ecology, University of Gothenburg, 2009.
Duarte CM, Sand-Jensen K. Seagrass colonization: biomass development and shoot demography in Cymodocea nodosa patches. Mar Ecol Prog Ser. 1990;67:97–103.
Agostini S, Pergent G, Marchand B. Growth and primary production of Cymodocea nodosa in a coastal lagoon. Aquat Bot. 2003;76:185–93.
Cunha AH, Duarte CM. Biomass and leaf dynamics of Cymodocea nodosa in the Ria Formosa lagoon, South Portugal. Bot Mar. 2007;50:1–7.
Peduzzi P, Vukovic A. Primary production of Cymodocea nodosa in the Gulf of Trieste (Northern Adriatic Sea): a comparison of methods. Mar Ecol Prog Ser. 1990;64:197–207.
Sghaier YR, Zakhama-Sraieb R, Charfi F. Spatio-temporal dynamics and biomass of Cymodocea nodosa in Bekalta (Tunisia, Southern Mediterranean Sea). Mar Ecol. 2017. https://doi.org/10.1111/maec.12383.
Terrados J, Medina-Pons FJ. Inter-annual variation of shoot density and biomass, nitrogen and phosphorus content of the leaves, and epiphyte load of the seagrass Posidonia oceanica (L.) Delile off Mallorca, western Mediterranean. Sci Mar. 2011;75(1):61–70.
Keulen M. Water flow in seagrass ecosystems. PhD thesis, Murdoch University, 1998.
Collier CJ, Lavery PS, Masini RJ, Ralph PJ. Morpjhological, growth and meadow characteristics of the seagrass Posidonia sinuosa along a depth-related gradient of light availability. Mar Ecol Prog Ser. 2007;337:103–15.
Fraser MW, Kendrick GA. Belowground stressors and long-term seagrass declines in a historically degraded seagrass ecosystem after improved water quality. Sci Rep. 2017. https://doi.org/10.1038/s41598-017-14044-1.
Jones BL, Cullen-Unsworth LC, Unsworth RFK. Tracking nitrogen source using δ15N reveals human and agricultural drivers of seagrass degradation across the British Isles. Front Plant Sci. 2018;9:133. https://doi.org/10.3389/fpls.2018.00133.
Vieira VM, Engelen AH, Huanel OR, Guillemin M-L. Linear-in-the-parameters Oblique Least Squares: a case study with the estimation of density-dependent survival in algae with isomorphic biphasic life-cycles. PLoS ONE. 2016;11(12):e0167418. https://doi.org/10.1371/journal.pone.0167418.
VV collected the data, developed the models and software, designed and performed the data analysis, interpreted the results and wrote the article. JCC collected the data, interpreted the results and reviewed the article. IEL collected the data and reviewed the article. All authors read and approved the final manuscript.
Our gratitude to the researchers who shared their data with us. Their contributions were fundamental for the success of our study.
The authors declare that they have no competing interests.
Availability of data and materials
The dataset supporting the conclusions of this article is included within the Additional file 2.
Consent to publish
Ethics approval and consent to participate
VV and IEL were funded by ERDF Funds of the Competitiveness Factors Operational Programme-COMPETE and national funds of the FCT—Foundation for Science and Technology under the project UID/EEA/50009/2013. This work was supported by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (JCC, Ciências do Mar 1137/2010); Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (JCC, FAPERJ-E-26/111.574/2014 and E26/201.286/2014); and Conselho Nacional de Desenvolvimento Científico e Tecnológico (JCC, CNPq- 307117/2014-6). The funders took no part in the design of the study, in the collection, analysis, and interpretation of data, and in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.