Skip to main content

The effect of parameter variability in the allometric projection of leaf growth rates for eelgrass (Zostera marina L.) II: the importance of data quality control procedures in bias reduction



Eelgrass grants important ecological benefits including a nursery for waterfowl and fish species, shoreline stabilization, nutrient recycling and carbon sequestration. Upon the exacerbation of deleterious anthropogenic influences, re-establishment of eelgrass beds has mainly depended on transplantation. Productivity estimations provide valuable information for the appraisal of the restoration of ecological functions of natural populations. Assessments over early stages of transplants should preferably be nondestructive. Allometric scaling of eelgrass leaf biomass in terms of matching length provides a proxy that reduces leaf biomass and productivity estimations to simple measurements of leaf length and its elongation over a period. We examine how parameter variability impacts the accuracy of the considered proxy and the extent on what data quality and sample size influence the uncertainties of the involved allometric parameters.


We adapted a Median Absolute Deviation data quality control procedure to remove inconsistencies in the crude data. For evaluating the effect of parametric uncertainty we performed both a formal exploration and an analysis of the sensitivity of the allometric projection method to parameter changes. We used parameter estimates obtained by means of nonlinear regression from crude as well as processed data.


We obtained reference leaf growth rates by allometric projection using parameter estimates produced by the crude data, and then considered changes in fitted parameters bounded by the modulus of the vector of the linked standard errors, we found absolute deviations up to 10 % of reference values. After data quality control, the equivalent maximum deviation was under 7 % of corresponding reference rates. Therefore, the addressed allometric method is robust. Even the smaller sized samples in the quality controlled dataset produced better accuracy levels than the whole set of crude data.


We propose quality control of data as a highly recommended step in the overall procedure that leads to reliable allometric surrogates of eelgrass leaf growth rates. The proliferation of inconsistent replicates in the crude data points towards the importance of discarding incomplete leaves. We also recommend avoiding errors in estimating the biomass of small leaves for which precision of the used analytical scale might be an issue.


Eelgrass is a relevant seagrass species that distributes worldwide in estuaries and nearshore environments. Eelgrass meadows provide habitat and foraging grounds for marine animals, buffers the shoreline from erosion, filter the water, and oxygenate the sediments. Recently at a global scale, deleterious effects derived from anthropogenic influences, have been exacerbated to such an extent that in the absence of remediation efforts, the important ecological services resulting from the permanence of eelgrass meadows could be irreversibly lost. Due to its ecological relevance eelgrass has been the subject of intense research and conservation efforts that are mainly carried out by means of transplanting plots. The measurement of biomass and productivity, provide key information for the evaluation of the overall status of a given eelgrass population. The noticeable growth form of this species makes the average rate of leaf growth per shoot-day measured over a growing interval of length ∆t determinant of overall productivity. In what follows we denote the biomass of an individual eelgrass leaf at time t through the symbol w(t) and its corresponding length by means of l(t). Similarly, we will symbolize the observed values of the biomass of leaves in shoots by means of w s (t) and the average rate of leaf growth per shoot-day through L g (t, ∆t). Conventional techniques for the assessment of L g (t, ∆t), require intensive sampling that involves the removal of shoots, and then tedious, time consuming dry weight measurement procedures in the laboratory. Even though the elimination of eelgrass shoots linked to a typical evaluation activity does not infringe damage to natural populations, the effects of shoot removal could be severe for transplanted plots. Therefore, in an overall program that leads to eelgrass conservation it is important to include nondestructive assessment methods [1, 2].

Bivariate allometric scaling relationships between measured quantities X and Y expressed as a power function of the form Y = βX α, appear in a myriad of research problems in physics, biology, and earth and planetary sciences [311]. The parameter α is named allometric exponent and β termed normalization constant. Particularly, Echavarria-Heras et al. [1] and Echavarria-Heras et al. [2], stressed the significance of the model

$$ w(t)=\beta l{(t)}^{\alpha }, $$

in the adaptation of an allometric proxy for L g (t, ∆t). This allometric surrogate will be here denoted by means of the symbol L ga (α, β, t, ∆t), and its explicit formulae derived in Appendix A.

The appropriateness of this device to provide truly nondestructive assessments depends on the time invariance of the parameters α and β, because estimates previously fitted at a site could be used to readily obtain surrogates for currently observed values of L g (t, ∆t) Several studies show that within our geographical region the parameters α and β can be considered as time invariant ([1, 1114]). But even though α and β are statistically invariant within a given region, environmental influences are expected to induce a relative extent of variability on local estimates of α and β [11] and as it was stated by Echavarria-Heras et al. [2] this could propagate significant uncertainties on L ga (α, β, t, ∆t) values. Indeed the results of Appendix B corroborate that these proxies could produce significantly accurate assessments for L g (t, ∆t) only in case the fitting of equation (1) yields highly precise estimates of the parameters α and β.

One important factor contributing to the uncertainty of the estimates of the parameters α and β, in equation (1) associates to the set of biological influences that could have been disregarded when assuming that eelgrass leaf biomass depends solely on linked length. Another important source of uncertainty relates to the failure of the model of equation (1) to handle environmental effects on a proper way. But the high values of determination coefficients as well as results of residual analysis corresponding to fittings of the model of equation (1) to independent data sets collected in our geographical region ([1, 2, 11, 13, 15, 16]) reveal that on spite of its simplicity, equation (1) is a paradigm that provides highly consistent estimates of the involved parameters.

But beyond uncertainties that biological or environmental influences could convey on estimates of the parameters α and β in equation (1), there could be other factors affecting their precision. Moreover, the results of Savage et al. [17], Hui and Jackson [18], Packard and Birchard [19], illustrate the relevance of factors like measurement error, data quality, sample size, and regression method on the precision of estimates of the parameters in a bivariate allometric scaling relationship. Moreover, Hui and Jackson [18], Packard and Birchard [19] and Packard and Boardman [20] concluded that in a bivariate allometric scaling nonlinear regression provides better estimates of the allometric exponent than these commonly estimated by linear regressions. Moreover, nonlinear regression was the analysis method used to obtain the significant estimates of the parameters α and β reported in all the eelgrass studies performed at our geographical region. Therefore, in the present settings we also considered nonlinear regression as an analysis method. However, in comparison with previously reported results ([2, 11, 15]) applying this analysis method to the present data set, resulted in an smaller value for the determination coefficient. Besides, in the current crude data, we detected a proliferation of inconsistent replicates, in the whole showing a noticeable discrepancy to the dispersion pattern that is normally expected for a scaling relationship of the form given by equation (1). Then, we assumed that on the light of the proven consistency of this model, the unusual spreading observed in the present data could on its own; provide an explanation for the reduction in determination coefficient obtained here. This view was strengthened by the fact that the discussed anomalous spreading was not observed for the data sets used in previous analysis of the fitting of a scaling relationship like equation (1) at our study site [1116]. And, it was perhaps a lack of standardization in processing tasks that led to a relatively larger incidence of measurement error in the present data set. Therefore, we considered that the exploration of data quality effects on the precision of estimates of the parameters α and β in equation (1) is a fundamental step in the overall procedure that leads to reliable allometric surrogates for the assessment of eelgrass leaf growth rates. Since, the aforesaid study had not been yet done; we made here an attempt to fill this gap.

Eelgrass leaf biomass is lognormally distributed, and the sizes of the different groups of replicates that make up our data are reduced. Therefore, we removed inconsistent replicates from the set of crude data by using a robust median absolute deviations cleaning procedure. We also relied on a sensitivity analysis study in order to assess the improvement in the accuracy of the L ga (α, β, t, ∆t) proxy that may possibly be associated to whatever reduction on the uncertainties of α and β could be derived from an improvement in data quality. Besides, we explored the corresponding effects of sample size, by randomly drawing differently sized samples out of both the crude and the processed data sets and then comparing the values of the determination coefficients for the fittings of equation (1) to the selected samples. In the Methods section we further substantiate the data cleaning approach, the procedure for estimating the effects of sample size as well as the steps of the performed sensitivity analysis. In the Results and Discussion section we elaborate on the relevance and limitations of this study. In the Conclusions section we presented the summary and potential implications of our findings.


Raw data

For the aims of the present research we assembled an extensive data set containing measurements taken on a total of 10412 individual eelgrass leaves collected in San Quintin Bay Baja California as previously reported in ([1, 2, 1114]). Crude data includes measurements of length (mm), width (mm) and dry weight (gr).

Data quality control procedures

Short [21] and Gaeckle and Short [22] estimated leaf biomass in eelgrass by using an isometric weight to length ratio, and the appropriateness of allometric methods in eelgrass research has been validated for a number of independent data sets (e.g. [1, 2, 1116] and references therein). Particularly, Echavarria-Heras et al. [2] formally demonstrated that the conspicuous leaf architecture and growth form of eelgrass makes leaf length a reliable allometric descriptor of the associated biomass. Indeed, results show that for independent data sets collected in our geographical region the model of equation (1) always produced reliable fits ([1, 2, 1113, 15, 16]). Moreover, the results of Solana-Arellano et al. [11] stating that the parameters associated to the scaling relationship of equation (1) can be considered invariant within a given geographical region, endorse that this model is highly consistent. This makes it reasonable to assume that the true relationship linking w(t) and l(t) should conform to a power function like dependency and consequently we could expect a typical dispersion for the linked data. In other words, the inherent spreading pattern should fit in a dominant power function trend masked by expected stochastic variability induced by intrinsic biological factors and environmental forcing. Particularly, at our study site the fitting of equation (1) produced a determination coefficient of R 2 = 0.90 for the data set addressed in [2]. Meanwhile, for data set in [11] as well as for that reported in [15] the corresponding value was R 2 = 0.92, with consistency of residuals verified for all these fits. Since in comparison with these results, for the present set of crude data the fitting of the model of equation (1) produced a reduced determination coefficient (R 2 = 0.81) we concluded that differences in determination coefficients could be explained by inconsistencies in the present data that could be linked to factors beyond biological influences or local environmental forcing. Moreover, by using a direct survey of the dispersion pattern shown by the present crude data (Fig. 1) we observed that departure of a given replicate from an expected power function-like trend was more pronounced for leaves with smaller sizes. Since, this effect was not observed in related data sets collected in our region, we judged that these inconsistencies could be in a fundamental way tied to the considerably larger extension of the present data set, which required processing performed by several technicians and it was perhaps a lack of standardization in these tasks that explains the present proliferation of inconsistencies that should be removed in order to prevent their permanence from exerting a reduction in the precision of estimates for the parameters in the scaling relationship of equation (1). Notwithstanding, the judgment to remove an outlier or inconsistent measurement in a data set, it is necessary to be able to identify its occurrence. Data cleaning procedures tasks are commonly performed by using a mean plus or minus three standard deviations method. This technique is based on the property of a normal distribution for which 99.87 % of the data appear within this range [23], but its application in the present settings presents difficulties. It firstly assumes that the distribution of data is normal while eelgrass leaf biomass is lognormally distributed [2]. Secondly, both the mean and standard deviation are strongly impacted by outliers [24]. Thirdly, our data is composed by groups of measurements that include a reduced number of leaf biomass replicates, and the mean plus or minus three standard deviations method is very unlikely to detect outliers in small samples [25]. This fact made it also inconvenient to use alternative data cleaning techniques like studentized residuals, the hat elements, Cook’s distance, or the Mahalanobis distance procedure [26]. Alternatively, the median of a group of data is totally immune to the sample size and a robust estimator of scale, reason why we adapted a Median Absolute Deviation (MAD) method [27] which provided the present criteria for the removal of inconsistencies in groups of more than 10 replicates of biomass values associated to a given leaf length. By eliciting the consistency of the allometric scaling of equation (1), we assumed that the addressed length-to-weight variation pattern should conform to a power function-like trend Therefore, individual leaf biomass measurements that in an initial data analysis were found to severely deviate from this pattern in sets under 10 replicates were removed from the set of unprocessed data. In summary, by analyzing the spreading of leaf biomass values in a leaf length-to-weight plot, we detected replicates that we supposed represented unduly deviations from an expected power function-like variation pattern. By applying the explained data cleaning procedures we removed these from the set of crude data, because in light of the proven consistency of the scaling relationship of equation (1) these inconsistent replicates could have been caused by errors due to a lack of standardization in leaf length or biomass measurements, or perhaps, explained by errors owed to faulty equipment for dry weight assessment or even explained by incorrect recordings.

Fig. 1
figure 1

Observed leaf length and biomass values. A plot of leaf length versus leaf biomass data reveals a power function-like trend masked by a great variability in observed leaf weight replicates for a given leaf size. Departure from an expected power function-like trend is more pronounced for smaller and larger leaves. The dispersion pattern shown suggests that data processing errors might have induced unduly biased leaf biomass values

In order to formalize the present MAD data cleaning procedure, we arranged the crude data in different groups G(l) = {w G1(l), …, w Gn (l)} formed by an observed leaf length l and associated leaf biomass replicates w G1(l), …, w Gn (l). For almost all groups of replicates G(l) we observed several leaf biomass values that parted from the expected variability pattern and were considered as inconsistencies. For each group G(l) we firstly obtained its median denoted by means of the symbol MED{w G1(l), …., w Gn (l)} or simply by means of M(G(l)) for short. Then, for each replicate w Gj (l) in (G(l) we calculated its absolute deviation from the group median δ Gj (l), that is, δ Gj (l) = |w Gj (l) − M(G(l))|. Similarly, we obtained the median of the set of absolute deviations denoted by the symbol MED{ δ G1(l), …, δ Gn (l)}. Following, Huber [28] and also recalling that eelgrass leaf biomass values are log-normally distributed [2], we obtained the Median Absolute Deviation of a group G(l) denoted here through MAD(G(l)) and given by

$$ MAD\left(G(l)\right)= bMED\left\{\ {\delta}_{G1}(l),\dots, {\delta}_{Gn}(l)\right\}, $$

where b = 1/Q (0.75), being Q (0.75) the 0.75 quantile of the lognormal distribution.

For the removal of inconsistent replicates in a group G(l) we used the decision criterion

$$ M\left(G(l)\right)-T\cdot MAD\left(G(l)\right)<{w}_j(l)<M\left(G(l)\right)\kern0.5em +T\cdot MAD\left(G(l)\right), $$

where T is the rejection threshold that following Miller [24], we set at a value of T = 3. For groups G(l) under ten replicates we applied a direct data cleaning procedure by removing replicates that we considered were severely deviated from the central power function-like trend.

Sensitivity analysis

In order to evaluate the influences that the uncertainties of the parameters α and β convey in the performance of L ga (α, β, t, ∆t), we used both a formal study and also simulation runs. The analytical exploration is presented in Appendix B. For the simulation study, initially, at equally spaced iteration (i.e., sampling) times t, we drawn uniformly distributed random numbers representing a number of shoots retrieved NS (t, ∆t) and the number nl(s) of leaves that each one of these shoots holds. Afterwards, the lengths l(t + ∆t) of each one the leaves in a shoot were simulated according to a fitted lognormal distribution, (p > 0.05), with a mean of 207.24 mm and a standard error of 3.8 [2]. Similarly, the corresponding leaf length increments, ∆l, were also extracted from a fitted lognormal distribution (p > 0.05) with a mean of 140.04 mm and a standard error of 7.3 [2]. The values for these increments coincided with the length of the whole leaf, as it occurs in newly produced leaves. Subsequently, we utilized the simulated leaf lengths and estimates \( \widehat{\alpha} \) and \( \widehat{\beta} \) of the parameters α and β in equation (1), in order to produce the proxy growth rates L ga (α, β, t, ∆t). In order to produce estimates \( \widehat{\alpha} \) and \( \widehat{\beta} \) for the parameters α and β we followed Hui and Jackson [18], Packard and Birchard [19] and Packard and Boardman [20] and fitted equation (1) to the present data sets using an iterative nonlinear least-squares method rather than the traditional approach of linearizing the equation through a logarithmic transformation of data. All fittings were performed using the Matlab Statistics Toolbox.

To study the sensitivity of L ga (α, β, t, ∆t) to changes in the parameters α and β we selected as reference values the estimates \( \widehat{\alpha} \) and \( \widehat{\beta} \) and produced changing values α ep for \( \widehat{\alpha} \) and β eq for \( \widehat{\beta} \) namely,

$$ {\alpha}_{ep} = \widehat{\alpha}+\varDelta {\alpha}_p, $$
$$ {\beta}_{eq}=\widehat{\beta}+\varDelta {\beta}_q, $$

with the values ∆α p and ∆α q satisfying

$$ \left|\varDelta {\alpha}_p\right| = p\cdot stde\left(\widehat{\alpha}\right) $$


$$ \left|\varDelta {\beta}_q\right| = q\cdot stde\left(\widehat{\beta}\right), $$

where \( stde\left(\widehat{\alpha}\right) \) and \( stde\left(\widehat{\beta}\right) \) stand for the standard errors of \( \widehat{\alpha} \) and \( \widehat{\beta} \) respectively and p and q are numbers satisfying, 0 < p ≤ 1 and 0 < q ≤ 1.

Therefore, the fluctuating values α ep and β eq for \( \widehat{\alpha} \) and \( \widehat{\beta} \) are scaled through proportions p and q of their standard errors. In our analysis we also relied on a parameter change index denoted by means of ρ(p, q), and defined through

$$ \rho \left(p,\ q\right) = \sqrt{{\left(\widehat{\beta}-{\beta}_{eq}\right)}^2+{\left(\widehat{\alpha}-{\alpha}_{ep}\right)}^2}. $$

Particularly, the maximum value that ρ(p, q) can attain as p and q vary is denoted through the symbol ρ max (p, q) and is given by,

$$ {\rho}_{max}\left(p,q\right) = \sqrt{ste{\left(\widehat{\alpha}\right)}^2+ste{\left(\widehat{\beta}\right)}^2\kern0.75em } $$

In what follows we will also use the mean value of ρ(p, q) denoted here by ρ av (p, q) and given by

$$ {\rho}_{av}\left(p,q\right)=\sqrt{stde{\left(\alpha \right)}^2+ stde{\left(\beta \right)}^2/2}. $$

The simulated leave lengths l(t + ∆t), corresponding increments ∆l the NS(t, ∆t) and nl(s) values, and the parameter estimates \( \widehat{\alpha} \) and \( \widehat{\beta} \) produced by means of equation (17) a reference trajectory \( {L}_{ga}\left(\widehat{\alpha}, \widehat{\beta},t,\varDelta, t\right) \). It turns out that every picked value of ρ(p, q) yields different pairs (α ep , β eq ), each one associated to a couple (∆α p , ∆β q ) that comply with the condition of equation (8). For each pair(α ep , β eq ), the simulated leave and shoot data and equation (17) shape a trajectory L ga (α ep , β eq , t, ∆t). The average of the values that these trajectories produce at each sampling time t is then calculated and its value denoted through 〈L ga (α ep , β eq , t, Δt)〉 ρ . Following the procedure, requires to calculate the deviations between 〈L ga (α ep , β eq , t, Δt)〉 ρ and \( {L}_{ga}\left(\widehat{\alpha}, \widehat{\beta},t,\varDelta t\right) \), that are denoted by means of the symbol δL gaρ (α ep , β eq , t, Δt) and obtained through,

$$ \delta {L}_{ga\rho}\left({\alpha}_{ep},{\beta}_{eq},,,t,,\varDelta t\right)={\left\langle {L}_{ga}\left({\alpha}_{ep},{\beta}_{eq},,t,\varDelta t\right)\right\rangle}_{\rho }-{L}_{ga}\left(\widehat{\alpha},\ \widehat{\beta},t,\varDelta t\right). $$

We also need to compute the average through time of the δL gaρ (α ep , β eq , t, Δt) deviations whose output is represented here through the symbol 〈δL gaρ (α ep , β eq , t, Δt)〉 t . Finally, it is necessary to obtain the time average of the values of the reference trajectory \( {L}_{ga}\left(\widehat{\alpha}, \widehat{\beta},t,\varDelta t\right) \). The resulting record is denoted using the symbol \( {\left\langle {L}_{ga}\left(\widehat{\alpha}, \widehat{\beta},t,\varDelta t\right)\right\rangle}_t \). These statistics are then used to produce the relative deviation index ϑ(Δα p Δβ q ) ρ given by

$$ \vartheta {\left(\varDelta {\alpha}_p,\ \varDelta {\beta}_q\ \right)}_{\rho }=\frac{\left|{\left\langle \delta {L}_{ga\rho}\left({\alpha}_{ep},\ {\beta}_{eq},t,\varDelta t\right)\right\rangle}_t\right|\kern0.5em }{{\left\langle {L}_{ga}\left(\widehat{\alpha},\ \widehat{\beta},t,\varDelta t\right)\right\rangle}_t}. $$

The value of ϑ(Δα p Δβ q ) ρ provides a measure of the sensitivity of the reference trajectory \( {L}_{ga}\left(\widehat{\alpha}, \widehat{\beta},t,\varDelta t\right) \) to a change of tolerance ρ(p, q) on the pair (\( \widehat{\alpha} \), \( \widehat{\beta} \)). Moreover, by calculating ϑ(Δα p , Δβ q ) max defined through,

$$ \vartheta {\left(\varDelta {\alpha}_p,\ \varDelta {\beta}_q\right)}_{max} = {max}_{\rho}\left\{\vartheta {\left(\varDelta {\alpha}_p,\ \varDelta {\beta}_q\right)}_{\rho}\right\}, $$

we can get in what percentage of \( {\left\langle {L}_{ga}\left(\widehat{\alpha}, \widehat{\beta},t,\varDelta t\right)\right\rangle}_t \) the maximum absolute deviation between L ga (α ep , β eq , t, Δt) and \( {L}_{ga}\left(\widehat{\alpha}, \widehat{\beta},t,\varDelta t\right) \) amounts.

Sample size effects

In order to explore the extent of sample size influences on \( {L}_{ga}\left(\widehat{\alpha}, \widehat{\beta},t,\varDelta t\right) \) we performed the fitting of equation (1), using differently sized samples of leaf length and biomass (l(t),w(t)) data which were randomly drawn from the available datasets. For each sample of size n, we obtained the determination coefficient (r 2 n ) as well as the values of the estimators of the parameters \( {\widehat{\alpha}}_n \) and \( {\widehat{\beta}}_n \), their respective standard errors \( std{e}_n\left(\widehat{\alpha}\right) \) and \( std{e}_n\left(\widehat{\beta}\right) \) and the values of the upper bound ρ max (p, q) n and ϑ(Δα p Δβ q  ) ρmax .

Results and discussion

In seagrass research, allometric equations have provided useful empirical models aimed to the representation of biologically relevant traits Y in terms of an easily measured variable X (e.g. [2933]). This is readily exemplified by equation (1) that provides a convenient representation of eelgrass leaf biomass w(t) in terms of associated length l(t). But the realm of allometric methods in eelgrass research is not only limited to empirical formulations that provide convenient proxies for leaf biomass estimations. Indeed from a theoretical standpoint, allometric approaches have provided for instance, a framework for the evaluation of the suitability of the leaf-biomass-to length ratio as a tool for nondestructive leaf biomass evaluations [15] and have also endorsed the substantiation of the plastochrone method of leaf growth assessments [16]. Moreover, equation (1) is the foundation of the derivation of the L ga (α, β, t, Δt) proxy of equation (17) a paradigm aimed to the indirect assessment of eelgrass leaf growth rates ([1, 2, 12]).

The effectiveness of the L ga (α, β, t, ∆t) construct for providing truly accurate and non-destructive assessments for the observed average rate of leaf growth per shoot-day L g (t, ∆t) depends on both, the time invariability of the parameters α and β, and on the accuracy of their estimates. From a general standpoint, the variability of the allometric exponent \( \alpha \) has been the focus of theoretical and empirical studies because it often seems to have a constant value specific to a particular biological relationship (e.g. [3439]). On the other hand, several studies have provided evidence that support certain variability in the exponent of allometric scaling laws (e.g. [4045]). Accordingly, the value of the normalization constant \( \beta \) is thought to be characteristic of species or populations [46]. And the variability observed in the normalization constant, explained as a differential response to environmental conditions [34, 35, 4749]. Particularly, for eelgrass, Solana-Arellano et al. [11] analyzed independent data set collected in different geographical regions to conclude that no universal values can be found for the allometric parameter α in equation (1), suggesting as well that this scaling relationship might be considered static, thus implying that local factors determine the extent of the variability of the actual values of the parameters α and β. On the presence of this variability associated uncertainties are expected to spread inaccuracies in eelgrass leaf growth rates produced by means of equation (17).

Figure 1 displays the variability of observed leaf lengths and linked biomasses. The survey of the crude data reveals a notorious proliferation of replicates associated to an observed leaf length. Data was arranged into 755 different groups G(l) = {w G1(l),…,w Gn (l)}. Figure 2 displays the variation of the number of replicates for a given leaf size. We can observe that the number of replicates in groups G(l) decreases with leaf size. This is probably due to the fact that separation from shoots by drag forces is more pronounced for longer and older leaves than it is for small and younger leaves. Figure 3 shows examples of the dispersion within G(l) groups. Moderated differences in replicated leaf biomass values associating to a given length can be explained by the inherent stochastic variability along with normal systematic errors introduced by data processing. Nevertheless, an exploration of the dispersion pattern of leaf biomass values in the present crude data reveals replicates representing unduly deviations from the expected power function-like trend. These inconsistent measurements are more visible for leaves under 100 mm long and also for those longer than 800 mm. For groups with 10 or more replicates and according to the present MAD data quality control criterion, whenever an individual leaf biomass replicate w Gj (l) differed in magnitude from the group median by more than 3 times the MAD it was considered as an inconsistency and it was rejected [50]. Similarly, the power function-like trend data cleaning criteria described in the methods section permitted the removal of inconsistent leaf biomass values in groups with less than 10 replicates. Figure 4 presents a plot of the leaf biomass values that remained after applying data quality control. This procedure retained a representative number of replicates showing the normally expected variability of eelgrass leaf biomass values ([1, 2, 1116]). Moreover, in order to exemplify differences in quality among equivalent data sets we applied the present data cleaning procedure to equivalent data reported in Echavarria-Heras et al. [12]. For easy of presentation we chose leaf length class intervals covering the extent of observed leaf lengths in both data sets. For each one of the lengths of in these class intervals we applied the present MAD procedure and calculated the percentage of discarded replicates. Then we obtained the average percentage of discarded leaves that associate to each length class interval. Figure 5 compares the average percentages of discarded replicates per leaf length class for the present data and for the Echavarria-Heras et al. [12] data set. For leaf lengths lying in the first class interval (0 ≤ l(t) < 50) the average percentage of rejected replicates in the present data set was 30 % while the Echavarria-Heras et al. [12] data reported no inconsistent replicates for the same interval. And generally, the percentages of discarded replicates in the present data set were markedly bigger than those associated to the Echavarria-Heras et al. [12] data. This comparison can explain why the fitting the model of equation (1) to the present crude data set produced an smaller determination coefficient than that reported in Echavarria-Heras et al. [12].

Fig. 2
figure 2

Leaf lengths and numbers of associated leaf biomass replicates. The distribution of the number of replicates of leaf biomasses for a given length decreases with leaf length. Maximum number of replicates was 107 and corresponded to a leaf length of 3 mm

Fig. 3
figure 3

Examples of the distribution of leaf biomass replicates for a given leaf length. a Leaf biomass replicates linked to l = 50 mm, (b) Leaf biomass replicates corresponding to l = 101 mm, and (c) Leaf biomass replicates associated to l = 200 mm. In these plots red lines stand for the group median: M(G(l)), yellow lines represent the M(G(l)) - 3∙MAD(G(l)) threshold and green lines show the M(G(l)) + 3∙MAD(G(l) threshold

Fig. 4
figure 4

Plot of quality controlled data. This plot shows the distribution of replicates of eelgrass leaf biomass that remained after applying the addressed quality control procedures. About 30 % of replicates in the raw data set were found inconsistent

Fig. 5
figure 5

Comparison of results of present data cleaning procedures applied to the current and an independent data set. This figure exhibits the results of applying the present data cleaning procedures applied to the current data and another equivalent data set reported in [12]. For both data sets we formed leaf length class intervals covering the extent of observed values, with the extension of each class interval set at 50 mm. The first length class interval is labeled by 1 and represents the set of leaf length lying within interval 0 < l(t) ≤ 50. Class 2 associates to the interval 50 < l(t) ≤ 100 and so one (horizontal axis). For each set of replicates associated to a given length in a class interval we obtained the percentage of removed inconsistencies, then for each class interval we obtained the average of these percentages (vertical axis). For all of the considered length class intervals we observed relatively greater average percentages of data inconsistencies in the present crude data

Observations show that the order relationship 0 < β < α holds [1, 2, 11, 12] and that it is also reasonable to consider the variation ranges of ∆α and ∆β satisfying: |∆α| < α and |∆β| < β. By embracing these assumptions in Appendix B we derived explicit formulae for the deviations δL ga (ΔαΔβ, t, Δt) that changed values of the parameters α and β having the form α e  = α + ∆α and β e  = β + ∆β induce in the reference values L ga (α, β, t, ∆t). These are given by equations (21). By assuming that the observed order relationship for the parameters α and β holds and also by considering the aforesaid expected variation ranges for ∆α and ∆β, we found that the proxies L ga (α, β, t, ∆t) for the average rate of leaf growth per shoot-day, will be overestimated by L ga (α e , β e , t, ∆t) values, primarily when the ordered pair (∆α, ∆β) lies on the domain ∆α > 0 and ∆β > 0. Nevertheless, as it is explained in Appendix B, each positive value of ∆β, can be associated to a set of negative values of ∆α for which L ga (α, β, t, ∆t) will be also overestimated by L ga (α e , β e , t, ∆t). Correspondingly, L ga (α e , β e , t, ∆t) underestimates L ga (α, β, t, ∆t) whenever the pair (∆α, ∆β) is placed inside the region -α < ∆α < 0 and -β < ∆β < 0, but this time, each value of ∆β satisfying -β < ∆β < 0 can be associated to a set of positive values of ∆β for which L ga (α, β, t, ∆t) will be also underestimated by L ga (α e, β e, t, ∆t).

By fitting of equation (1) to the crude leaf length and weight data we obtained r 2 = 0.81 and estimates, \( \widehat{\alpha}=1.32367 \) with \( stde\left(\widehat{\alpha}\right) = 0.0143 \) for α and of \( \widehat{\beta}=0.000014 \) with \( stde\left(\widehat{\beta}\right)=1.4e-6 \) for β. Therefore, according to equation (8) we have 0 ≤ ρ(p, q) ≤ 0.0143. Using the \( \widehat{\alpha} \) and \( \widehat{\beta} \) estimates we obtained ρ max (p, q) = 0.0143 (cf. Eq. (9)) and ρ av (p, q) = 0.0072 (cf. Eq. (10)). By means of the simulated leave and shoot data and equation (17) we estimated the \( {L}_{ga}\left(\widehat{\alpha},\ \widehat{\beta},t,\varDelta t\right) \) reference trajectory. Afterwards for fixed values of the parameter change index ρ(p, q) we used equations (6) and (7) to produce values for the ∆α p and ∆β p increments complying with the condition of equation (8). This allowed the characterization of the associated L ga (α ep , β eq , t, ∆t) trajectories. Figure 6 displays trajectories obtained for the case in which ρ(p, q) takes its average value ρ av (p, q). These simulation outputs are consistent with the results of the appendix setting domains of variation of the pair (∆α, ∆β) where L ga (α, β, t, ∆t) values are overestimated or underestimated. Similarly, Fig. 7 shows the time variation of the average deviations 〈δL ga (Δα p Δβ q , t, Δt)〉 ρ obtained for ρ(p, q) = ρ av (p, q). Moreover, Fig. 8 shows that whenever 0 ≤ ρ(p, q) ≤ ρ max (p, q) the values of the relative deviation index ϑ(Δα p Δβ q ) ρ are bounded above by a value ϑ(Δα p Δβ q max  = 0.01, that is, the maximum absolute deviation between L ga (α ep , β eq , t, ∆t) and \( {L}_{ga}\left(\widehat{\alpha},\ \widehat{\beta},t,\varDelta t\right) \) amounts to 10 % of \( {\left\langle {L}_{ga}\left(\widehat{\alpha},\ \widehat{\beta},t,\varDelta t\right)\right\rangle}_t \). In turn, when fitting of equation (1) to the quality controlled data we obtained r 2 = 0.91 and estimates \( \widehat{\alpha}=1.403 \) with \( stde\left(\widehat{\alpha}\right)=0.0120 \) for α and of \( \widehat{\beta}=8.462e-006 \) with \( stde\left(\widehat{\beta}\right) = 6.4300e-007 \), for β. Therefore, according to equation (9) this time we have ρ max (p,q) = 0.01. In turn, performing the sensitivity study of equations (4) through (13) we obtained ϑ(Δα p Δβ q ) max  = 0.07, which amounts to 7 % of \( {\left\langle {L}_{ga}\left(\widehat{\alpha},\ \widehat{\beta},t,\varDelta t\right)\right\rangle}_t \) (see Fig. 9).

Fig. 6
figure 6

Examples of deviations between allometrically projected leaf growth rates. We portrait several trajectories of allometrically projected leaf growth rates L ga (α ep , β eq ,t, ∆t) in units of mg. leaf −−1. These trajectories were produced by variations α ep and β ep in fitted parameters \( \widehat{\alpha} \) and \( \widehat{\beta} \), complying with the condition ρ(p, q)2 = stde(\( \widehat{\alpha} \))2 + stde(\( \widehat{\beta} \))2/4. We show the reference trajectory \( {L}_{ga}\left(\widehat{\alpha},\ \widehat{\beta},t,\varDelta t\right) \) in red. The number of days elapsed between sampling times is 15

Fig. 7
figure 7

Example of the trajectory of deviations δL gaρ (α ep , β eq , t, Δt). We show an example of the behavior of the trajectory of deviations δL gaρ (α ep , β eq , t, Δt). This is obtained by subtracting at each time t, the value of the reference trajectory \( {L}_{ga}\left(\widehat{\alpha},\ \widehat{\beta},t,\varDelta t\right) \) from the average of the values of the different trajectories L ga (α ep , β eq , t, Δt) produced by a fixed value of the parameter change index ρ(p, q). This average is denoted by means of the symbol 〈L ga (α ep , β eq , t, Δt)〉 ρ (cf. Equation (11)). The number of days between sampling times is 15

Fig. 8
figure 8

The behavior of the relative deviation index ϑ(Δα p Δβ q  ) ρ . This plot shows the behavior of the relative deviation index for different values of ρ(p, q). This index ϑ(Δα p Δβ q  ) ρ is obtained by dividing the absolute value of the average through time of the deviations between the 〈L ga (α ep , β eq , t, Δt)〉 ρ and \( {L}_{ga}\left(\widehat{\alpha},\ \widehat{\beta},t,\varDelta t\right) \) rates, by the time average of the \( {L}_{ga}\left(\widehat{\alpha},\ \widehat{\beta},t,\varDelta t\right) \) values, (cf. Equation (12)). The relative deviation index ϑ(Δα p Δβ q  ) ρ provides a measure of the sensitivity of the allometric method L ga (α, β, t, Δt) to changes in the parameters α and β. The values of the relative deviation index ϑ(Δα p Δβ q  ) ρ are bounded above by a value ϑ(Δα p Δβ q  )  max  = 0.01. (cf. Equation (13))

Fig. 9
figure 9

Comparison of relative deviation index for crude and quality controlled data. This plot compares the values for the relative deviation index ϑ(Δα p Δβ q  ) ρ calculated as the parameter change index ρ(p, q) varies. Values of ϑ(Δα p Δβ q  ) ρ produced by crude data are shown using blue lines and those corresponding to quality controlled data in red lines

This sensitivity exploration, sets the accuracy of the allometric method of equation (17) to be mainly dependent on the extent of the upper bound ρ max (p, q) of the parameter change index ρ(p, q), This means that an improvement in the quality of the fit of equation (1) reducing the magnitude of the upper bound for ρ(p, q) will certainly increase the accuracy of the proxies L ga (α, β, t, ∆t). This study also demonstrates the importance of data quality control in the reduction of ρ max (p, q). Therefore, within the realm of parametric uncertainties as determined by of the values of \( stde\left(\widehat{\alpha}\right) \) and \( stde\left(\widehat{\beta}\right) \) obtained for the quality controlled data set, the addressed projection method can be considered robust relative to numerical differences in the estimators \( \widehat{\alpha} \) and \( \widehat{\beta} \).

Figure 10 presents the variation of ϑ(Δα p Δβ q )  max depending on sample size for both the set of raw data and that resulting after quality control procedures. As expected Fig. 10 shows that for both data sets ϑ(Δα p Δβ q )  max decreases as sample size increases, but these plots also reveal that the effect of sample size on the accuracy of the L ga (α, β, t, Δt) proxy is more pronounced for the raw data set. Moreover, even the smallest sized sample in the quality controlled data set (n = 1000) produces a smaller value for ϑ(Δα p Δβ q ) max than that induced by the largest sized sample in the raw data set (n = 10412). This adds on to the aforesaid on the prominent role of data quality control in improving the accuracy of eelgrass leaf growth rate assessments obtained by means of L ga (α, β, t, ∆t).

Fig. 10
figure 10

The effects of sample size on relative deviation index values. For a sample of a given size, the maximum value that the relative deviation index ϑ(Δα p Δβ q  ) ρ attains is smaller for the quality controlled data. Even the smallest sized sample in the quality controlled data produces a higher precision that the largest sized one in the crude data


The results of the present study highlight the important role that the precision of parameter estimates linked to the allometric model of equation (1) plays in the overall suitability of the allometric proxy for eelgrass leaf growth rates L ga (α, β, t, ∆t) given by equation (17). The basic allometric model of equation (1) has been consistently identified using different data sets, and a property of invariance for the involved parameters statistically verified. But on spite of model consistency, the evaluation of the extent on what factors such as data quality sampling size, and analysis method can influence parameter estimates is an important entry in the comprehensive analysis of an allometric scaling relationship like equation (1). Moreover, the recommendation by Hui and Jackson [12] stating that for variables that are measured with errors, it is important to obtain accurate estimates with repeated measurements on similar individuals grown under similar conditions highlights the relevance of controlling factors that could affect data quality in our settings. Since the present results were obtained following the recommendation of Hui and Jackson [18], Packard and Birchard [19] and Packard and Boardman [20] concerning analysis method, and since sample size in our study was optimal, it is reasonable to assume that errors in data processing could explain a drop in the determination coefficient for the fitting of equation (1) to crude data collected at our study site from R 2 = 92 obtained by Solana-Arellano et al. [11] and Echavarria-Heras et al. [15] to a value of R 2 = 81 produced by the present data set. And the fact that after the data cleaning procedures performed on the present data we obtained for the resulting determination coefficient a value of R 2 = 91 seems to endorse our judgment. Therefore, our results seem to confirm that data quality can be raised as a crucial factor explaining precision of estimates for the parameters α and β in equation (1). Moreover, the present data quality control procedure reveals the remarkable influence of inconsistencies in reducing precision of allometric projections. Indeed, as it is shown in Fig. 9 the accuracy of the L ga (α, β, t, ∆t) proxy improved after processing data irregularities. Our contribution provides guidelines that could prevent the proliferation of undesirable data quality related effects on the accuracy of allometric projections of eelgrass leaf growth rates. It was perhaps a lack of standardization that could explain a relatively reduced quality in the present data set. In other words the unduly deviations from an expected power function--like spreading could be probably tied to errors in leaf length estimation, faulty equipment for dry weight assessment, or even due to the incorrect recording of measurements. Nevertheless, the spread displayed in Fig. 4, shows notorious remnant variability for leaves longer than 500 mm, even after the application of data quality control procedures. This suggests that a data quality approach must be mandatory not only during the processing stage but even since the gathering step of forming a data set. Indeed for leaves with the smaller sizes, that normally yield reduced biomasses we might expect estimation errors that are imputable to the precision of the analytical scale used to obtain leaf dry weight estimations, along with the reading and/or recording of the output. But longest leaves remain exposed longer to environmental effects such as drag forces or herbivore. These influences could remove sections of the leaves rendering underestimated biomasses for a given leaf length. Therefore, special care must be also taken when processing longer an older leaves which are more often damaged of trimmed to such an extent that the biomass they report do not necessarily associate to the true length to weight relationship. We anticipate that the overlooking of these important issues as well as other data processing mistakes like rounding off errors or even incorrect data recording could explain the anomalous proliferation of inconsistencies (about 30 %) found in the present crude data set. Our findings also reveal that a data quality approach could optimize sampling endeavors, this because as shown in Fig. 10, even the smallest sized sample of the quality controlled data could entail better model identification results than the whole unprocessed data set. In summary, the present findings show how data quality control can be conceived as a fundamental step in the overall procedure that leads to reliable nondestructive allometric proxies of eelgrass leaf growth rates. And, it is also worth to point out, that it can be misleading to implicitly assume that a given set of data should conform to an idealized allometric relationship and then use a data cleaning procedure to ensure the fit. It can be misleading because the dispersion pattern itself could act as an indicator of eelgrass population status. But, we must emphasize that in accomplishing the present approach we are not neglecting the fact that there could be other unknown influences explaining the dispersion observed in the present crude data, but we are in fact acknowledging that the proven consistency and invariance of the addressed allometric relationship for eelgrass leaf biomass, is expected to provide a dependable account of pertinent causal influences, being these linked to the biology of the plant or associated to environmental forcing, thus highlighting the relevance of data quality influences in determining an unexpected variation pattern in the data as the one detected here.


  1. Echavarría-Heras HA, Lee KS, Solana-Arellano ME, Franco-Vizcaino E. Formal analysis and evaluation of allometric methods for estimating above-ground biomass of eelgrass. Ann Appl Biol. 2011;159(3):503–15.

    Article  Google Scholar 

  2. Echavarría-Heras HA, Solana-Arellano ME, Leal-Ramírez C, Franco-Vizcaíno E. An allometric method for measuring leaf growth in eelgrass, Zostera marina, using leaf lenght data. Bot Mar. 2013;56(3):275–86.

    Article  Google Scholar 

  3. Newman MEJ. Power laws, Pareto distributions and Zipf’s law. Contemporary Phys. 2005;46:323–51.

    Article  Google Scholar 

  4. Marquet PA, Quiñones RA, Abades S, Labra F, Tognelli M. Scaling and power-laws in ecological systems. J Exp Biol. 2005;208:1749–69.

    Article  PubMed  Google Scholar 

  5. West GB, Brown JH. The origin of allometric scaling laws in biology from genomes to ecosystems: towards a quantitative unifying theory of biological structure and organization. J Exp Biol. 2005;208:1575–92.

    Article  PubMed  Google Scholar 

  6. Harris LA, Duarte CM, Nixon SW. Allometric laws and prediction in estuarine and coastal ecology. Estuar Coast. 2006;29:343–7.

    Article  Google Scholar 

  7. Filgueira R, Labarta U, Fernández-Reiriz MJ. Effect of condition index on allometric relationships of clearance rate in Mytilus galloprovincialis Lamarck. 1819. Rev Biol Mar Oceanogr. 2008;43(2):391–8.

    Article  Google Scholar 

  8. Kaitaniemi P. How to derive biological information from the value of the normalization constant in allometric equations. PLoS One. 2008;3(4):e1932.

    PubMed Central  Article  PubMed  Google Scholar 

  9. Martin RD, Genoud M, Hemelrijk CK. Problems of allometric scaling analysis: examples from mammalian reproductive Biology. J Exp Biol. 2005;208:1731–47.

    Article  PubMed  Google Scholar 

  10. González-Trilla G, Borro MM, Morandeira NS, Schivo F, Kandus P, Marcovecchio J. Allometric scaling of dry weight and leaf area for spartina densiflora and spartina alterniflora in two southwest Atlantic Saltmarshes. J Coast Res. 2013;29(6):1373–81.

    Article  Google Scholar 

  11. Solana-Arellano ME, Echavarría-Heras HA, Leal-Ramírez C, Lee KS. The effect of parameter variability in the allometric projection of leaf growth rates for eelgrass (Zostera marina L.). Lat Am J Aquat Res. 2014;42(5):1099–108.

    Article  Google Scholar 

  12. Echavarría-Heras HA, Solana-Arellano ME, Franco-Vizcaino E. An allometric method for the projection of eelgrass leaf biomass production rates. Math Biosci. 2010;223(1):58–65.

    Article  PubMed  Google Scholar 

  13. Solana-Arellano E, Echavarría-Heras H, Díaz-Castañeda V, Flores-Uzeta O. Shoot biomass assessments of the marine phanerogam Zostera marina for two methods of data gathering. Am J Plant Sci. 2012;3:1541–5.

    Article  Google Scholar 

  14. Solana-Arellano E, Borbón-González DJ, Echavarria-Heras H. A general allometric model for blade production in Zostera narina L. Bull South Cal Acad Sci. 1998;97:39–49.

    Google Scholar 

  15. Echavarria-Heras H, Solana-Arellano E, Lee KS, Hosokawa S, Franco-Vizcaino E. An evaluation of leaf biomass: length ratio as a tool for nondestructive assessment in eelgrass (Zostera marina L.). Sci World J. 2012;2012:543730.

    Article  Google Scholar 

  16. Echavarría-Heras H, Solana-Arellano E, Leal-Ramírez C, Castillo O. Using allometric procedures to substantiate the plastochrone method for eelgrass leaf growth assessments Theoretical Biology and Medical Modelling. 2013; 10:34,

  17. Savage VM, Gillooly JF, Woodruff WH, West GB, Allen AP. The predominance of quarter-power scaling in biology. Funct Ecol. 2004;18:257–82.

    Article  Google Scholar 

  18. Hui D, Jackson RB. Uncertainty in allometric exponent estimation: a case study in scaling metabolic rate with body mass. J Theor Biol. 2007;249:168–77.

    Article  PubMed  Google Scholar 

  19. Packard GC, Birchard GF. Traditional allometric analysis fails to provide a valid predictive model for mammalian metabolic rates. J Exp Biol. 2008;211:3581–7.

    Article  PubMed  Google Scholar 

  20. Packard GC, Boardman TJ. Model selection and logarithmic transformation in allometric analysis. Physiol Biochem Zool. 2008;81:496–507.

    Article  PubMed  Google Scholar 

  21. Short FT. Effects of sediment nutrients on seagrasses: literature review and mesocosm experiment. Aquat Bot. 1987;27(1):41–57.

    CAS  Article  Google Scholar 

  22. Gaeckle JL, Short FT. Plastochrone method for measuring leaf growth in eelgrass, Zostera marina L. Bull Mar Sci. 2002;71(3):1237–46.

    Google Scholar 

  23. Howell DC. Statistical methods in human sciences. New York: Wadsworth; 1998.

    Google Scholar 

  24. Miller J. Reaction time analysis with outlier exclusion: bias varies with sample. Q J Exp Psychol. 1991;43(4):907–12.

    CAS  Article  Google Scholar 

  25. Cousineau D, Chartier S. Outliers detection and treatment: a review. Int J Psychol Res. 2010;3(1):58–67.

    Google Scholar 

  26. Stevens JP. Outliers and Influential data points in regression analysis. Psychol Bull. 1984;95(2):334–44.

    Article  Google Scholar 

  27. Leys C, Ley C, Klein O, Bernard P, Licata L. Detecting outliers: do not use standard deviation around the mean, use absolute deviation around the median. J Exp Soc Psychol. 2013;49(4):764–6.

    Article  Google Scholar 

  28. Huber PJ. Robust statistics. New York: John Wiley; 1981.

    Book  Google Scholar 

  29. McRoy CP. Standing stock and other features of eelgrass (Zostera marina) populations on the coast of Alaska. J. Fish. Res. Bd. Can. 1970; 27:1811–1821.

  30. Patriquin DG. Estimation of growth rate, production and age of the marine angiosperm, Thalassia testudinum. Mar Biol. 1973;15(1):35–46.

    Article  Google Scholar 

  31. Jacobs RPWM. Distribution and aspects of the production and biomass of eelgrass, Zostera marina L. at Roscoff France. Aquat Bot. 1979;7:151–72.

    Article  Google Scholar 

  32. Hamburg SP, Homann PS. Utilization of growth parameters of eelgrass, Zostera marina, for productivity estimates under laboratory and in situ conditions. Mar Biol. 1986;93:299–303.

    Article  Google Scholar 

  33. Duarte CM. Allometric scaling of seagrass form and productivity. Mar Ecol Prog Ser. 1991;77:289–300.

    Article  Google Scholar 

  34. Winter JE. A critical review on some aspects of filter feeding in lamellibranchiate bivalves. Haliotis. 1976;7:71–87.

    Google Scholar 

  35. Bernard FR. Physiology and mariculture of some northeastern Pacific bivalve molluscs. Can Spec Pub Fish Aquat Sci. 1983;63:1–24.

    Google Scholar 

  36. Jones HD, Richards OG, Southern TA. Gill dimensions, water pumping rate and body size in the mussel Mytilus edulis L. J Exp Mar Biol Ecol. 1992;155:213–37.

    Article  Google Scholar 

  37. Enquist BJ, Brown JH, West GB. Allometric scaling of plant energetic and population density. Nature. 1998;395:163–5.

    CAS  Article  Google Scholar 

  38. Niklas KJ, Enquist BL. Invariant scaling relationships for interspecific plant biomass production rates and body size. Proc Natl Acad Sci. 2001;98:2922–7.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  39. West GB, Woodruff VH, Brown HJ. Allometric scaling of metabolic rate from molecules and mitochondria to cells and mammals. Proc Natl Acad Sci. 2002;99:2473–8.

    PubMed Central  Article  PubMed  Google Scholar 

  40. Riisgård HU. No foundation of a ‘3/4 power scaling law’ for respiration in biology. Ecol Lett. 1998;1:71–3.

    Article  Google Scholar 

  41. Atanasov AT, Dimitrov DB. Changes of the power coefficient in the ‘metabolism-mass’ relationship in the evolutionary process of animals. Biosystems. 2002;66:65–71.

    Article  PubMed  Google Scholar 

  42. Bokma F. Evidence against universal metabolic allometry. Funct Ecol. 2004;18:184–7.

    Article  Google Scholar 

  43. Muller-Landau HC, Condit RS, Chave J, Thomas SC, Bohlman SA, Bunyavejchewin S, et al. Testing metabolic ecology theory for allometric scaling of tree size, growth and mortality in tropical forests. Ecol Lett. 2006;9:575–88.

    Article  PubMed  Google Scholar 

  44. Reich PB, Tjoelker MG, Machado JL, Oleksyn J. Universal scaling of respiratory metabolism, size and nitrogen in plants. Nature. 2006;439:457–61.

    CAS  Article  PubMed  Google Scholar 

  45. White CR, Phillips NF, Seymour RS. The scaling and temperature of vertebrate metabolism. Biol Lett. 2006;2:125–7.

    PubMed Central  Article  PubMed  Google Scholar 

  46. Niklas KL. Plant allometry: the scaling of form and process. Chicago: University of Chicago; 1994.

    Google Scholar 

  47. West GB, Brown JH, Enquist BJ. A general model for the origin of allometric scaling laws in biology. Science. 1997;276:122–6.

    CAS  Article  PubMed  Google Scholar 

  48. Gillooly JF, Brown JH, West GB, Savage VM, Charnov EL. Effects of size and temperature on metabolic rate. Science. 2001;293:2248–51.

    CAS  Article  PubMed  Google Scholar 

  49. Brown JH, Gupta L, Li BL, Milne BT, Restrepo C, West GB. The fractal nature of nature: power laws, ecological complexity and biodiversity. Phil Trans R Soc Lond B. 2002;357:619–26.

    Article  Google Scholar 

  50. Miller JN. Rouges and suspects: how to tackle outliers. Amc Thechnical Briefs No. 39, 2009.

Download references

Author information

Authors and Affiliations


Corresponding author

Correspondence to Héctor Echavarría-Heras.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

H conceived, designed and performed the research. C performed formal derivations and simulation and numerical tasks. E performed the statistical analysis and interpretation of data. N supervised the whole research and revised the manuscript critically both at statistical and formal levels. All authors read and approved the final manuscript.


Appendix A

In what follows a subscript s will be used to label a generic eelgrass shoot, holding a number nl(s) of leaves with combined biomass w s (t). In time, if ∆w l (t, ∆t) stand for the increment in biomass that is gained by an individual leaf over a growing interval [t, t + ∆t], then denoting by means of L sg (t, ∆t) the resulting average growth rate of the leaves on the shoot s we will have,

$$ {L}_{sg}\left(t,\varDelta t\right)=\frac{{\displaystyle {\sum}_{nl(s)\ }}\varDelta {w}_l\left(t,\varDelta t\right)\ }{\varDelta t}, $$

and consequently if L g (t, ∆t) denotes the linked average rate of leaf growth per shoot-day, we then have,

$$ {L}_g\left(t,\varDelta t\right)=\frac{{\displaystyle {\sum}_{NS\left(t,\varDelta t\right)}}{L}_{sg}\left(t,\varDelta t\right)}{NS\left(t,\varDelta t\right)}, $$

where ∑ NS(t,Δt) indicates summation of the shoots collected over the marking interval [t, t + ∆t] and being NS(t, ∆t) their number. Moreover, as it is thoroughly explained by Echavarría-Heras et al. [2], we can use equation (1) in order to derive an allometric approximation for L sg (t,∆t), which we here denote through the symbol L sga (α, β, t, Δt) and formally express by

$$ {L}_{sga}\left(\alpha, \beta, t,\varDelta t\right)=\frac{{\displaystyle {\sum}_{nl(s)}}\varDelta {w}_{la}\left(t,\varDelta t\right)}{\varDelta t}, $$

where Δw la (t, Δt) is an allometric proxy for Δw l (t, Δt).

It turns out that equation (1) yields,

$$ \varDelta {w}_{la}\left(t,\varDelta t\right) = \beta l{\left(t+\varDelta t\right)}^{\alpha } - \beta l{(t)}^{\alpha } $$

then, factoring l(t + Δt)α and taking into account that l(t + Δt) = l(t) + Δl where Δl stands for the increment in leaf length gained over the interval [t, t + Δt] we get,

$$ \varDelta {w}_{la}\left(t,\varDelta t\right) = \beta l{\left(t+\varDelta t\right)}^{\alpha}\left(1-{\left(1 - {\rho}_l\ \left(t,\varDelta t\right)\right)}^{\alpha}\right) $$


$$ {\rho}_l\ \left(t,\varDelta t\right)=\frac{\varDelta l}{l\left(t+\varDelta t\right)}. $$

And by letting

$$ \delta \left(t,\varDelta t\right)=\left(1-{\left(1-{\rho}_l\ \left(t,\varDelta t\right)\right)}^{\alpha}\right) $$

we equivalently have

$$ \varDelta {w}_{la}\left(t,\varDelta t\right) = \beta l{\left(t+\varDelta t\right)}^{\alpha}\delta \left(t,\varDelta t\right). $$

Therefore, from (15) we obtain,

$$ {L}_{sga}\left(\alpha, \beta, t,\varDelta t\right) = \frac{{\displaystyle {\sum}_{nl(s)}}\beta l{\left(t+\varDelta t\right)}^{\alpha}\delta \left(t,\varDelta t\right)}{\varDelta t} $$

Similarly, if L ga (α, β, t, Δt) denotes the allometric proxy for L ga (t, Δt), this is given by

$$ {L}_{ga}\left(\alpha, \beta, t,\varDelta t\right)=\frac{{\displaystyle {\sum}_{NS\left(t,\varDelta t\right)}}{L}_{sga}\left(\alpha, \beta, t,\varDelta t\right)\ }{NS\left(t,\varDelta t\right)}, $$

and in turn we have

$$ {L}_g\left(t,\varDelta t\right) = {L}_{ga}\left(\alpha, \beta, t,\varDelta t\right)+{\in}_{ga}, $$

with the term ϵ ga standing for the involved approximation error. It is worth to point out that the derivation of the results of equation (18) could have been done using leaf area a(t) in place of l(t).

Appendix B

In this appendix we firstly derive explicit formulae for the deviations δL ga (∆α, ∆β, t, ∆t) that a change ∆α and ∆β in α and β will induce in L ga (α, β, t, ∆t). Using the derived forms we obtain domains of variation of ∆α and ∆β where L ga (α, β, t, ∆t) are underestimated or overestimated. Without loss of generality, we will assume through that the inequality 0 < β < α holds, as it regularly occurs in fittings of equation (1) to eelgrass leaf biomass and length data. Also, for the sake of mathematical tractability, without sacrificing pertinence of the analysis, regarding the increments ∆α and ∆β we will assume that |∆α| ≤ α and |∆β| < β.

By definition above we have,

$$ \delta {L}_{sga}\left(\varDelta \alpha,\ \varDelta \beta, t,\varDelta t\right)={L}_{sga}\left(\alpha +\varDelta \alpha,\ \beta +\varDelta \beta, t,\varDelta t\right)-{L}_{sga}\left(\varDelta \alpha,\ \varDelta \beta, t,\varDelta t\right) $$

And from equation (17) one gets

$$ \delta {L}_{ga}\left(\varDelta \alpha, \varDelta \beta, t,\varDelta t\right)=\frac{{\displaystyle {\sum}_{NS\left(t,\varDelta t\right)}}\delta {L}_{sga}\left(\varDelta \alpha,\ \varDelta \beta, t,\varDelta t\right)\ }{NS\left(t,\varDelta t\right)} $$

Now equation (16) yields,

$$ \delta {L}_{sga}\left(\varDelta \alpha,\ \varDelta \beta, t,\varDelta t\right)=\frac{{\displaystyle {\sum}_{nl(s)}}\beta l{\left(t+\varDelta t\right)}^{\alpha}\delta \left(t,\varDelta t\right){\mu}_{Lg}\left(\varDelta \alpha, \varDelta \beta, l\left(t+\varDelta t\right)\right)}{\varDelta t}, $$


$$ {\mu}_{Lg}\left(\varDelta \alpha, \varDelta \beta, l(t)\right)=\left(\left(\frac{\beta +\varDelta \beta }{\beta}\right)\left(\left(1-{\left(1-\frac{\varDelta l}{l\left(t+\varDelta t\right)}\right)}^{\alpha +\varDelta \alpha}\right)\ \right)\ \delta {\left(t,\varDelta t\right)}^{-1}l{\left(t+\varDelta t\right)}^{\varDelta \alpha } - 1\right), $$

which can be equivalently written as

$$ {\mu}_{Lg}\left(\varDelta \alpha, \varDelta \beta, l(t)\right)=\frac{\varphi \left(\varDelta \alpha \right)}{\theta \left(\varDelta \beta \right)}-1, $$


$$ \theta \left(\varDelta \beta \right)=\kern0.5em \frac{\beta }{\beta +\varDelta \beta }, $$


$$ \varphi \left(\varDelta \alpha \right)=\kern0.5em \frac{\left(1-{\left(\frac{l(t)}{l\left(t+\varDelta t\right)}\right)}^{\alpha +\varDelta \alpha}\right)l{\left(t+\varDelta t\right)}^{\varDelta \alpha }}{\left(1-{\left(\frac{l(t)}{l\left(t+\varDelta t\right)}\right)}^{\alpha}\right)}. $$

Now by assumption 0 < β < α and |∆α | ≤ α, and since observations show that the statement that l(t) being time increasing is consistent, for ∆t > 0 we necessarily have l(t + ∆t) > l(t) and therefore, the factor φ(∆α) as given above increases in its entire domain and also remains positive. Again the case ∆α = ∆β = 0 will set θ(∆β) = φ(∆α) = 1 (Fig. 11). Moreover, L ga (α, β, t, ∆t) will be overestimated by L ga (α + ∆α,β + ∆β,t) whenever the inequality

Fig. 11
figure 11

The behavior of the auxiliary factors θ(∆β) and φ(∆α). This plot presents the behavior of the auxiliary factors θ(∆β) and φ(∆α) defined by equations 24 and 25 respectively. Both ∆α and ∆β vary in the horizontal axis. Similarly the variations of θ(∆β) and φ(∆α) are projected in the vertical axis. We assumed that the order relationships 0 < β < α, and |∆β| < β, hold. Therefore, θ(∆β) remains positive and decreases and in turn we have φ(−α) = 0, and also that φ(α) remains positive and increasing for ∆α >−α. Particularly, for ∆α = ∆β = 0, we have φ(∆α) = θ(∆β) = 1

$$ {\mu}_{Lg}\left(\varDelta \alpha, \varDelta \beta, l(t)\right)>0 $$

holds, but since we also assumed that |Δβ | < β, as it has been already pointed out, the factor θβ) remains positive, by virtue of (23) we have that the statement in (26) implies

$$ \varphi \left(\varDelta \alpha \right) > \theta \left(\varDelta \beta \right). $$

Now for ∆α ≥ 0, φ(∆α) satisfies φ(∆α) ≥ 1 and increases, while for ∆β ≥ 0 the factor θ(∆β) satisfies θ(∆β) ≤ 1 and decreases, therefore the statement (27) will clearly hold in the domain ∆α > 0 and ∆β > 0. Notwithstanding, for ∆β = ∆β * fixed, and such that 0 < ∆β * < β, we have θ(Δβ *) ≥ 1/2, then since for -α ≤ ∆α < 0, we have 0 ≤ φ(Δα) ≤ 1, and φ(Δα) is continuous there, we can exhibit a value ∆α *, that satisfies, − α ≤ Δα * < 0, for which we have φ(Δα * ) = θ(Δβ *), and since φ(Δα) increases whenever ∆α * ≤ ∆α ≤ 0, we will have φ(Δα) ≥ θ(Δβ * ). Therefore, to each positive value of ∆β, we can associate a set of negative values of ∆α for which inequality (27) also holds.

Now, L ga (α + ∆α, β + ∆β, t) will underestimate L ga (α, β, t, ∆t), whenever

$$ {\mu}_{Lg}\left(\varDelta \alpha, \varDelta \beta, l(t)\right)<0 $$

is satisfied, but again since |∆β| < β, from (23) the statement above implies

$$ \varphi \left(\varDelta \alpha \right)<\theta \left(\varDelta \beta \right) $$

Now for -α ≤ ∆α ≤ 0, the factor φ(∆α) vanishes at ∆α = −α and then increases monotonically towards a value φ(0) = 1. Meanwhile, as it has been discussed before for -β < ∆β ≤0 the term θ(∆β) decreases monotonically towards a value θ(0) = 1. Therefore, inequality (29) will particularly hold in the domain ∆α < 0 and ∆β < 0. This time, for ∆β = ∆β τ fixed and satisfying -β < ∆β τ <0, we have θβ τ ) > 1 and since φ(∆α) is continuous, there exists ∆α τ positive, for which we have φα τ  ) = θβ τ ) and since φ(∆α) increases and changes continuously over the domain |∆α| < α, particularly for 0 ≤ ∆α < ∆α τ we will have, 1 ≤ φα) < θβ τ ). Therefore, each negative value of ∆β can be associated to a set of positive values of ∆α for which inequality (29) also holds.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Echavarría-Heras, H., Leal-Ramírez, C., Villa-Diharce, E. et al. The effect of parameter variability in the allometric projection of leaf growth rates for eelgrass (Zostera marina L.) II: the importance of data quality control procedures in bias reduction. Theor Biol Med Model 12, 30 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Eelgrass leaf growth rates
  • Allometric estimation
  • Parametric uncertainty effects
  • Data quality control