Next Article in Journal
Practical Use of Materials of Natural Origin as Loose-Fill Insulations in Open-Diffusion Constructions—Observation and Numerical Simulation
Next Article in Special Issue
Comparison of Carbon Storage in Forested and Non-Forested Soils in Tropical Wetlands of Caimanera, Colombia, and Llano, Mexico
Previous Article in Journal
Assessment of Sustainable Mobility Patterns of University Students: Case of Cameroon
Previous Article in Special Issue
Hydrolases Control Soil Carbon Sequestration in Alpine Grasslands in the Tibetan Plateau
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of Allometric Equations to Determine the Biomass of Plant Components and the Total Storage of Carbon Dioxide in Young Mediterranean Argan Trees

by
Assma Oumasst
1,2,
Fatima Ezzahra Tiouidji
1,3,
Salma Tabi
1,3,
Abdelaziz Zahidi
1,2,4,
Abdelhamid El Mousadik
2,
Aissam El Finti
2,
Abderrahmane Aitlhaj
5 and
Jamal Hallam
1,*
1
Regional Center of Agricultural Research of Agadir, National Institute of Agricultural Research, Avenue Ennasr, BP 415 Rabat Principale, Rabat 10090, Morocco
2
Laboratory of Biotechnology and Valorization of Natural Resources (LBVRN), Faculty of Sciences Ibn Zohr University, Agadir 80000, Morocco
3
Laboratory of Environmental, Ecological and Agro-Industrial Engineering (LGEEAI), Faculty of Science and Technology of Beni Mellal, Sultane Molay Slimane University, Beni Mellal 23000, Morocco
4
Laboratory of Entomology, Complex Horticol, Agronomic and Veterinary Institute Hassan II, Agadir 80000, Morocco
5
National Agency for the Development of Oasis and Argan Zones, Agadir 80000, Morocco
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(11), 4592; https://doi.org/10.3390/su16114592
Submission received: 16 February 2024 / Revised: 16 April 2024 / Accepted: 18 April 2024 / Published: 28 May 2024
(This article belongs to the Special Issue Soil Carbon Sequestration and Greenhouse Gas Emission)

Abstract

:
Morocco’s argan trees, native to forests, are now cultivated in large orchards within the Argan Biosphere Reserve, transforming “Arganiculture”. These orchards are anticipated to bolster carbon storage, but their precise contribution to carbon storage remains unclear. This study introduces for the first time allometric equations for estimating biomass in different components of argan plants within orchards. A total of 89 plant individuals, aged 2 to 6 years, were collected. Their diameter and total height were measured. The biomass, carbon content, and biomass carbon stock of each component were determined. The best-fit allometric equation incorporates diameter, height, growing years, and root-to-shoot ratio to estimate total biomass (R2 = 0.95). The estimated total carbon biomass stock ranged from 0.01 to 0.82 t CO2 ha−1 for plants, at a density of 200 plants ha−1. Between 2021 and 2023, the average annual carbon sequestration was 0.20 t CO2 ha−1 year−1. This model offers valuable tools for use when species-specific equations during the establishment growing stage are unavailable, enhancing carbon sequestration quantification for more reliable results and informing climate change mitigation strategies. The allometric parameters serve as benchmarks for trees resembling the argan tree. The methodology could be adapted for other forest plants undergoing conversion to orchard cultivation.

1. Introduction

Biomass estimation is one of the most crucial steps in estimating carbon sequestration in different biomass components [1,2]. It is important to determine the future projections of carbon dioxide (CO2) stored by biomass over the years [3]. This information is useful in ecology as it provides insights into ecosystem structure and function, serves as an indicator for sustainability, and helps in assessing ecosystem productivity and the global carbon budget [4,5].
Biomass estimation can be achieved using different methods, including the direct application of allometric equations, the multiplication of stand volume by wood density, or the use of the biomass expansion factor, which is a measure of volume-to-biomass ratio [4,5,6]. These methodologies have a common requirement for generalized allometric equations, although they differ in terms of predictor variables. The former method is often used for younger plants, while the latter is more appropriate for mature forest trees due to its rarity and adherence to ethical guidelines and regulations, aiming to mitigate any potential ecological impact on ecosystems [4,7].
There has been recent interest in the development of the allometric equation method for biomass estimation, which is a cost-effective technique and uses easily measurable variables such as diameter and height. Additionally, it has the potential to adapt to specific species, taking into account their unique characteristics and growth patterns [8]. These specialized models can also include site-specific factors, such as soil type, climate, and topography, which can enhance the accuracy, efficiency, and conciseness of biomass estimates and reduce errors in carbon biomass stock assessments [9,10,11,12].
This study focuses on creating allometric equations for estimating biomass in different components of the argan plant (Argania spinosa (L) Skeels), a vital species in Morocco supporting local economies and ecosystems. Argan trees, predominantly found in the 830,000-hectare Argan Biosphere Reserve (ABR), face challenges due to the increased demand for argan oil, impacting forest density and sustainability. To address this, Morocco has initiated the domestication of forest-origin species by establishing argan orchard perimeters in vulnerable areas [13]. This innovative project, based around the Arganiculture concept, aims to meet national commitments for environmental protection and sustainability outlined in the nationally determined contributions (NDCs). Arganiculture is expected to enhance carbon storage, mitigate soil erosion, and sustain the agroforestry system, supporting both mitigation and adaptation efforts. The development of allometric equations to estimate Argania biomass holds immense potential for quantifying the carbon sequestration capabilities within established orchards. These data would significantly contribute to policy formulation, enabling informed decisions on land use, conservation strategies, and carbon offset initiatives [14].
In earlier studies, Benzyane [15] and Belghazi [16] employed a semi-destructive method to estimate the biomass of argan trees in the Haha plateau (Essaouira region). Benzyane [15] generalized the equation for mature trees considered as part of a high forest, while Belghazi [16] focused on young stands classified as coppice pushed on old and already-well-established root systems. Determining the age of Argania trees in forest stands posed challenges; however, these studies relied on parameters like tree height and trunk circumference. The published biomass equations focus on large argan trees with diameters > 10 cm. No attempt has been made to explore the biomass estimation of argan trees at the establishment growing stage; these are considerably smaller and contribute little to the total stand biomass. In our current research, allometric equations focus on argan plants of known ages and are derived from nursery seedlings subsequently transplanted into the newly established Arganiculture orchards. By deriving allometric equations from nursery seedlings, this research aims to provide a more comprehensive understanding of the biomass of argan trees at the established growth stage, which can be identified as young argan plants. This approach contributes to a more accurate assessment of total stand biomass in argan forests. It enables more informed forest management and conservation strategies, improves carbon sequestration assessments, and contributes to better modeling of forest growth and development [14].
Accordingly, the present study attempts to (1) develop allometric equations for young argan plants and their components (e.g., leaves, stem, small branches, roots) to estimate their biomass, (2) apply these equations to estimate the total biomass in selected studied orchards, and (3) convert the estimated biomass to total carbon dioxide storage.

2. Materials and Methods

2.1. Study Area

The study was conducted in argan orchards across the marginal and vulnerable areas within the regions of Souss Massa, Marrakech-Safi, and Guelmim Oued Noun, all located within the ABR in the North Africa region. Figure 1 presents the orchard’s locations in the ABR. These areas are within three different climatic zones, including the semi-arid climatic zone, the arid climatic zone, and the pre-arid or Saharan climatic zone. The mean annual temperature is 18.05 °C, with average monthly temperatures ranging between 9.4 °C (December, January) and 26 °C (July, August). The mean annual precipitation is 301.6 mm, with the highest rainfall occurring from October to April, with extremely low precipitation occurring between June and August [17]. The argan soils in these areas are mostly lithosols and regosols, associated with fluvisols and saline soils on lowlands [18]. The argan plants studied exhibit different growing stages, ranging from 2 to 6 years.

2.2. Sampling and Field Measurement

At the beginning of the study in 2021, six individual argan plant samples were randomly collected from 11 experimental plots (each 0.5 ha) across 6 distinct argan orchard perimeters within the ABR, following the experimental design (Figure 2). In total, 66 young argan plants were collected. Plants in newly transplanted orchards were only 2, 3, and 6 years old, creating a data gap. The biomass trend of plants aged between 3 and 6 years could be difficult to interpret without data on intermediate ages. To avoid the risk of a false biomass trend due to missing age groups, additional plants were collected later in 2023, when they had reached 4 and 5 years of age. A total of 36 plants were collected from 6 plots, distributed across 3 distinct perimeters within the ABR. The potential seasonal differences between the initial and later data collection were considered. The plots maintained a density of 200 argan plants per hectare. The total height (H) and diameter (D) of the sampled plants were measured using a measuring tape and calliper. The diameter was measured at the base of the stem. Degraded plants without leaves or with broken stems were removed, resulting in a sample size of 89 young plants.

2.3. Laboratory Analysis

The leaves, stems, small branches, and roots were separated from each harvested plant. In the case of very small individuals, the stem and small branches were considered as one part [19]. The fresh biomass of each component was weighed separately. Roots were washed and cleaned of adhering soil, drained, and air-dried, and their weight were recorded. A sub-sample were then oven-dried at 70 °C (leaves) and 105 °C (stems, branches, and roots) until a constant weight was achieved. Root-to-shoot ratio was calculated in argan plants by dividing the root weight by the shoot weight. The carbon concentration of each component was measured using the ignition loss method [20]. The above- and below-ground carbon stock was calculated by multiplying the tree’s biomass by its carbon content [21]. To convert carbon stock to CO2 emissions, we multiplied it by 3.67.

2.4. Allometric Model Development and Statistical Analysis

The study analysed the relationship between leaf, stem, and root biomass as dependent variable, and specific independent variables including diameter (D), height (H), growing years (age), and root-to-shoot ratio [22]. A linear regression was applied (Equation (1)) and the data were checked for linear relationship, independence, homoscedasticity, and normality following the protocol described by Zuur et al. [23]. To enhance linearity and homoscedasticity, the regression models underwent the natural logarithm transformation of the variables (Equation (2)) [4,24]. Possible sources of error were checked for, and outliers were identified. Only data that were correctly identified were retained for this study. To remove bias when converting natural logarithmic values to biomass, a correction factor (CF) was derived for each equation (Equation (3)), following Sprugel’s [25] method. The correction factor is determined by the standard error of estimate (SEE) of the regression.
B i o m a s s = a X
L n B i o m a s s = l n a + b L n X
B i o m a s s = e a X b C F
where X is the independent variable (e.g., D), a and b are empirically estimated scaling factors, and CF is the correction factor.
Several studies have referred to the allometric models frequently used to estimate biomass [9,21,26,27,28]. The best allometric models were evaluated and selected based on the primary criteria as outlined by Chave et al. [4]. Table 1 summarises the primary criteria and their descriptions and interpretations. The allometric models was validated using the leave-one-out cross-validation (LOOCV) technique, in which one data point at a time is systematically left out and the model is trained on the remaining data points [29,30]. This iterative process is repeated for each data point in the dataset, enabling multiple rounds of training and validation to assess the models’ performance. The LOOCV technique assesses its predictive accuracy and reliability. The mean squared error (MSE) was adopted as robust metric for evaluating model performance within the LOOCV framework. The statistical modelling was performed using the lm() function from the R package, and the statistical comparisons were made using the basic R package under version R 3.2.3.
The study analyzed the effect of growing years and argan plant components on carbon content using a two-way analysis of variance (ANOVA) after checking for normality and homogeneity of variance. The statistical analysis was conducted using the SPSS program (IBM Corp. Released in 2019, Version 26), and significance was determined at p < 0.05 level.

3. Results

3.1. Biomass and Carbon Content

Table 2 shows that the mean D, H, and total dry biomass for selected young argan plants were 1.61 ± 1.28 cm, 0.68 ± 0.25 m, and 0.48 ± 0.85 kg, respectively. The biomass distribution was concentrated mainly in the stems, which represented 49% of the total dry biomass, with an average of 0.22 ± 0.39 kg. Roots represented 27% of the total dry biomass, with an average of 0.12 ± 0.26 kg, while the leaves represented 24% of the total dry biomass, with an average of 0.14 ± 0.22 kg. The mean root-to-shoot ratio for the study was 0.64, with a range of 2.67 for the smallest diameter classes and 0.69 for the largest.
Figure 3 illustrates the relationship between diameter and measured dry biomass, alongside height and measured biomass for different components of young argan plants. The analysis indicates that biomass is directly affected by both diameter and height. However, this effect is non-linear.
Figure 4 shows the mean carbon content among different components of young plants. The carbon content varied significantly (p < 0.001) based on the tree component and its age (Table 3). The root had a lower carbon concentration (33–55%), while the stem and leaf of the argan tree had a significantly higher carbon content (52–56%).
The study found a medium, yet significant, correlation (0.5 < r < 0.9) between the biomass of young argan plant components and growing years, which were the dependent and the independent variables, respectively, as shown in Figure 5. Additionally, a large, significant correlation was observed between root and stem biomass (r = 0.95) and between diameter and both stem and leaf biomass (r = 0.90 and r = 0.92, respectively). The results also indicated a negative correlation between the shoot-to-root ratio and years of growth (r = −0.30).

3.2. Biomass Model Generation

Table 4 presents models of allometric equations and their corresponding coefficients used to estimate the biomass of leaf, stem, and root. The statistical parameters of these allometric equation models for each young plant’s components are compared in Table 5. This study reported higher values of R2 and adjusted R2 (0.85 to 0.96 and 0.86 to 0.93, respectively) for leaf and root parts, while lower values (0.82 to 0.88) were recorded for stem parts of young plants.
For leaf biomass, model 9 was identified as the best-fitting regression equation, explaining 95% of the variance. It had the lowest values of RSE, AIC, and RSME compared to other models (RSE = 0.51, AIC= −116.45, RSME = 4.35) and no potential for heteroskedasticity. For stem biomass, only models 1, 2, 10, and 11 showed no heteroskedasticity of variance. Model 11 was selected as the best allometric model with the highest R2 and R2-adjusted value of 0.86 and the lowest RSE and AIC (RSE = 0.70, AIC= −60.06, RSME = 3.44). For root biomass, model 4 with three variables (D, H, and age) met all the selection criteria with an R2 of 0.90. However, when specific variables such as the root-to-shoot ratio (model 13) and aboveground biomass (ABG) (model 14) were added, the R2 increased to 0.93 and 0.92, respectively, while RSE decreased from 0.51 to 0.47 and 0.49, respectively. Furthermore, the AIC, RMSE, and MSE parameters decreased. Model 13 was selected as the most reliable and strong model for predicting root biomass.
The residuals of the three selected models were considered independent, as indicated by a Durbin–Watson statistic of 2 (DW = 2). Model 9 had an MSE of 0.27 for leaf biomass, Model 11 had an MSE of 0.51 for stem biomass, and Model 13 had an MSE of 0.24, representing the lowest values among the models presented in Table 5. The selected allometric models for leaf and roots were more reliable and efficient compared to the stem model (Table 5).
In all fitted models for all three biomass components, the residuals did not reveal heterogeneity of variance, as shown in Figure 6. Q-Q plots demonstrate a close alignment of residuals with the theoretical normal distribution, which supports model assumptions. Predicted vs. observed values plots show a strong linear relationship, indicating good model-prediction accuracy.

3.3. Biomass Carbon Stock

The average biomass per plant of different components varied across growing years, with a total biomass range of 0.03 to 2.34 kg per plant (Figure 7). Table 6 presents the total biomass of leaves, stems, roots of argan plants in the studied orchards. The total carbon biomass of 200 argan plants per hectare at the ages of 2 to 6 years ranges from 0.003 to 0.22 t C ha−1, corresponding to 0.01 to 0.82 t CO2 ha−1 (Figure 8). From 2021 to 2023, the average annual carbon sequestration was 0.20 tonnes of CO2 per hectare per year.

4. Discussion

Allometric models were developed to estimate the biomass of different components. including stems, leaves, roots, and the total biomass of argan plants within recently established argan orchards perimeters across different regions within the Argan Biosphere Reserve. The studied plants were 2 to 6 years old, corresponding to the establishment growing stage (encompassing the lag phase and the very early part of the exponential growth phase). They exhibited diameters ranging from 0.33 to 7.91 cm and a height between 0.29 and 1.42 m. The mean root-to-shoot ratio obtained was 0.64. The root-to-shoot ratio indicates that there are on average 0.64 units of root per 1 unit of aboveground biomass across the plants in the study. The variation in the root-to-shoot ratio across diameter classes indicates that smaller plants allocate a larger proportion of their resources to belowground structures compared to larger plants.
Young plants with smaller stems may require larger root systems for stability and support in arid and semi-arid climates. This could be explained by the need to acquire resources such as water and nutrients in the soil [33]. The ratio varies with environmental conditions such as precipitation, soil moisture, texture, and fertility and stand conditions such as age, height, forest type, or origin [34,35].
The study of carbon content in plant components identified differences in the mean carbon content of roots based on the age of the argan plants. In particular, the carbon content of the roots of argan plants aged between 4 and 5 years old was higher than that of other growing years. This variation could be attributed to natural fluctuations in carbon content throughout the growing season and across years, which are influenced by environmental factors such as temperature, rainfall, and nutrient availability. Harvesting at different times (2021 vs. 2023) may capture these variations. Furthermore, differences in management practices may also be a contributing factor. If there had been changes in management practices, such as fertilization or irrigation, between 2021 and 2023, then this could have had a greater impact on the root C content.
The study revealed that the roots of argan trees have a lower carbon concentration (33–55%) compared to the stems and leaves (52–56%). This could be explained by argan plant roots concentrating on acquiring resources and water directing less carbon towards the belowground parts of the plant. Additionally, age-related variations in carbon content can occur, with younger plants investing more carbon in rapid growth and development [33,36]. Mahmood et al. [26] found that the carbon content was between 42% and 47% in leaves and between 40 and 49% in branches, bark, and stems in tropical evergreen forests. Thomas et al. [37] reported a wide variation in carbon concentration of 45.7–60.7% in subtropical/Mediterranean species, and Ma et al. [38] found that the carbon content in crop roots was 38.20 ± 5.23%. In the global assessments, carbon content has been assumed to be 50% of tree biomass; however, recent studies indicate that this assumption is not accurate, with substantial variation in carbon content between tree species, as well as between tissue types [37].
The study results show a non-linear relationship between diameter and height, which suggests that a simple regression (linear model) might not accurately represent the data. These involve using non-linear equations (e.g., power functions, exponential functions) or applying transformations to the data to linearize the relationship. Among the transformation methods, log–log models were most often reported in many studies of different ecosystems [4,26,39].
Several studies in the literature recommend that biomass equations based on both diameter (D) and height (H) are more accurate than equations based only on D [4,40,41]. Also, the precision and the accuracy of allometric equations increased with the level of specificity of equations as found by Paul et al. [42]. This study included diameter, height, age, and root-to-shoot ratio for the root biomass model as specific independent variables. The results revealed that the best model for leaves included D, H, and age; D and age for stem biomass, and D, H, age, and root–shoot ratio for root biomass. Therefore, the fit model for total biomass included all the mentioned variables with the highest biomass variation (R2 = 0.95). On the other hand, the accuracy of biomass estimation largely depends on the appropriate selection of allometric models [26,43]. In the study, both models 3 and 9 performed well in predicting leaf biomass, with high R-squared adjusted values of 0.95. However, model 1 showed a lower AIC score than model 9, suggesting that it may capture the relationship between predictors and biomass slightly better. Model 3, on the other hand, shows a potential for heteroskedasticity, which means that the error variances may not be consistent across the data. This could have an impact on the reliability of the conclusions drawn from the model. Our findings agree with those of Mokany et al. [35]. They emphasize the need for reliable root–shoot ratios across various vegetation types to enhance the accuracy of root biomass estimates, particularly for purposes such as National Greenhouse Gas Inventories and carbon accounting and for studies of ecosystem dynamics.
Henry et al. [14] highlighted the limited consideration of root biomass in current carbon (C) inventories under an arid climate (only 1.3% of equations). They reported the challenges in measuring root biomass as the high cost and time required for the direct measurement of root biomass and the difficulty of extracting the entire root [14]. In this study, the roots of 5- and 6-year-old argan plants were found to have deeper roots and develop further into the soil compared to the above-ground parts. This limited the extraction of these root samples. This suggests that we should consider alternative methods for estimating root biomass for argan, such as fractal geometry, minirhizotrons, and ground-penetrating radar. The importance of accurate root biomass data is emphasized, as they can contribute to improving carbon inventories and enhancing our understanding of ecosystem dynamics.
In summary, this study formulated biomass allometric equations for each component of young argan plants, including the leaf (Equation (4)), stem (Equation (5)), root (Equation (6)), and total biomass (Equation (7)). The total biomass is calculated by summing the biomass of each component.
B i o m a s s L e a f = e 7.21 × ( D 2 H ) 0.6 × a g e 3.2 × 1.28
B i o m a s s S t e m = e ( 8.15 + 2.20 × D + 1.24 × a g e ( 0.35 × D × a g e ) ) × 1.41
B i o m a s s R o o t = e 5 × D 1.48 × H 0.4 × a g e 1.38 × R o o t   t o   s h o o t   r a t i o 0.31 × 1.26
B i o m a s s T o t a l   B i o m a s s = e 4.2 × D 1.36 × H 0.57 × a g e 1.67 × R a t i o   s h o o t   t o   r o o t 0.3 × 1.23
The study demonstrates that carbon biomass storage increases with argan plant age, ranging from 0.01 to 0.82 t CO2 ha−1. This highlights the significance of conserving older argan trees to improve carbon sequestration and inform conservation strategies. The total biomass of argan plants in the orchards studied was between 0.006 ± 0.002 and 0.47 ± 0.06 tonnes of dry biomass (DB) per hectare. The findings of the current study differ from those of the Belghazi [16] study in terms of reported biomass values. Belghazi [16] found that the total dry aboveground biomass of a young argan coppice aged 5 years old was 1.62 t ha−1 for a density of 75 coppice ha−1, corresponding to 21.6 kg per coppice [16]. This value is approximately 21 times higher than the current study’s value (1 kg per plant). The observed difference could be attributed to the disparity in stand characteristics. While Belghazi [16] studied coppice derived from well-established root systems, our study focused on plants grown from nursery seedlings subsequently transplanted into the field. In coppicing practice, trees are trimmed close to the ground, stimulating vigorous new growth from the well-established root system. Variations in site and plant density, competition for water availability, and nutrient limitations may also affect growth.
Between 2021 and 2023, the average annual carbon sequestration within the studied age range of argan plants was 0.20 t CO2 per hectare per year, representing the establishment phase of growth. The low rate observed in our findings could be attributed to the species-specific growth phase examined, in contrast to previous studies integrating the exponential growth phase. Environmental conditions, including soil properties, genetics, and arid and semi-arid climates with prolonged dry seasons, particularly those experienced in Morocco during the study period, could account for these differences. The observed carbon sequestration observed is lower than the reported values for some broader agroforestry systems. These systems have been found to have CO2 removal rates ranging from 0.8 to 15.6 t CO2 ha−1 year−1 during the first 20 years of growth [44]. Further research is required to investigate the long-term carbon sequestration potential of argan trees as they mature.

5. Conclusions

Allometric models were developed to estimate the biomass of different components including leaves, stems, roots, and total biomass of young argan plants (age ≤ 6 years). Additionally, the total carbon dioxide (CO2) storage during this establishment growth phase was calculated for the different components. In this study, we found that the performance of allometric equations varied among plants biomass components and the independent variable also varied. The highest performance of the allometric models was found in total biomass estimates with high explained variance (R2 = 0.95). Allometric equations provide biomass estimates based on diameter, height, age, and root-to-shoot ratio. This research has raised many questions in need of further investigation. Further work is needed to determine whether locally specific allometric equations should be developed to account for seasonal dynamics in estimating leaf biomass and to identify areas where fractal geometry can be used to estimate root biomass. This contribution can enhance the accuracy of biomass estimation models, leading to improved carbon inventories and a better understanding of ecosystem dynamics. The implementation of the new concept of Arganiculture through the recently established argan orchards in various regions of the ABR in Morocco presents an opportunity for research concerning argan plants. Sustained data collection and the monitoring of tree establishment, especially during the accelerated or exponential growth phase, will be crucial. This continued monitoring will aid in developing equations specific to various growth phases, encompassing exponential, linear, and deceleration phases. Such an approach will be essential to the comprehensive modelling and analysis of tree growth dynamics.

Author Contributions

Methodology and field work, A.O. and J.H.; data curation, A.O. and F.E.T. and S.T.; writing—original draft preparation, A.O.; writing—review and editing, A.O. and J.H.; visualization, A.E.M. and A.E.F. and A.A.; supervision, J.H. and A.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Green Climate Fund (GCF) grant number FP0022/Grant/DARED 2021-23 and the APC was funded by the Green Climate Fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are part of an ongoing PhD project and can be provided by the corresponding author upon request, due to restrictions concerning the originality of the PhD work.

Acknowledgments

This research was undertaken as part of the Moroccan Implementation Project for Research on Argan Trees, a component of the Arganiculture Development Project in Vulnerable Areas (“DARED”). Noureddine Ibnezzyn, along with his dedicated teams from the National Agency for the Development of Oasis and Argan Zones (ANDZOA), played an active role in this study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Saglam, S.; Ozdemir, E.; Ozkan, U.Y.; Demirel, T.; Makineci, E. Biomass Estimation of Aboveground Tree Components for Turkey Oak (Quercus Cerris L.) in South-Eastern Turkey. Environ. Monit. Assess. 2020, 192, 418. [Google Scholar] [CrossRef]
  2. Sebrala, H.; Abich, A.; Negash, M.; Asrat, Z.; Lojka, B. Tree Allometric Equations for Estimating Biomass and Volume of Ethiopian Forests and Establishing a Database: Review. Trees For. People 2022, 9, 100314. [Google Scholar] [CrossRef]
  3. Koornneef, J.; van Breevoort, P.; Hamelinck, C.; Hendriks, C.; Hoogwijk, M.; Koop, K.; Koper, M.; Dixon, T.; Camps, A. Global Potential for Biomass and Carbon Dioxide Capture, Transport and Storage up to 2050. Int. J. Greenh. Gas Control 2012, 11, 117–132. [Google Scholar] [CrossRef]
  4. Chave, J.; Andalo, C.; Brown, S.; Cairns, M.A.; Chambers, J.Q.; Eamus, D.; Fölster, H.; Fromard, F.; Higuchi, N.; Kira, T. Tree Allometry and Improved Estimation of Carbon Stocks and Balance in Tropical Forests. Oecologia 2005, 145, 87–99. [Google Scholar] [CrossRef] [PubMed]
  5. Yang, M.; Zhou, X.; Peng, C.; Li, T.; Chen, K.; Liu, Z.; Li, P.; Zhang, C.; Tang, J.; Zou, Z. Developing Allometric Equations to Estimate Forest Biomass for Tree Species Categories Based on Phylogenetic Relationships. For. Ecosyst. 2023, 10, 100130. [Google Scholar] [CrossRef]
  6. Forrester, D.I.; Tachauer, I.H.H.; Annighoefer, P.; Barbeito, I.; Pretzsch, H.; Ruiz-Peinado, R.; Stark, H.; Vacchiano, G.; Zlatanov, T.; Chakraborty, T. Generalized Biomass and Leaf Area Allometric Equations for European Tree Species Incorporating Stand Structure, Tree Age and Climate. For. Ecol. Manag. 2017, 396, 160–175. [Google Scholar] [CrossRef]
  7. Joosten, R.; Schumacher, J.; Wirth, C.; Schulte, A. Evaluating Tree Carbon Predictions for Beech (Fagus Sylvatica L.) in Western Germany. For. Ecol. Manag. 2004, 189, 87–96. [Google Scholar] [CrossRef]
  8. Picard, N.; Saint-Andre, L.; Henry, M. Manual for Building Tree Volume and Biomass Allometric Equations: From Field Measurement to Prediction; FAO: Rome, Italy, 2012. [Google Scholar]
  9. Chave, J.; Réjou-Méchain, M.; Búrquez, A.; Chidumayo, E.; Colgan, M.S.; Delitti, W.B.; Duque, A.; Eid, T.; Fearnside, P.M.; Goodman, R.C. Improved Allometric Models to Estimate the Aboveground Biomass of Tropical Trees. Glob. Change Biol. 2014, 20, 3177–3190. [Google Scholar] [CrossRef] [PubMed]
  10. Kettler, T.A.; Doran, J.W.; Gilbert, T.L. Simplified Method for Soil Particle-Size Determination to Accompany Soil-Quality Analyses. Soil Sci. Soc. Am. J. 2001, 65, 849–852. [Google Scholar] [CrossRef]
  11. Luo, Y.; Wang, X.; Ouyang, Z.; Lu, F.; Feng, L.; Tao, J. A Review of Biomass Equations for China’s Tree Species. Earth Syst. Sci. Data 2020, 12, 21–40. [Google Scholar] [CrossRef]
  12. Maulana, S.I.; Wibisono, Y.; Utomo, S. Development of Local Allometric Equation to Estimate Total Aboveground Biomass in Papua Tropical Forest. Indones. J. For. Res. 2016, 3, 107–118. [Google Scholar] [CrossRef]
  13. Green Climate Fund. FP022: Development of Arganiculture Orchards in Degraded Environment (DARED). Available online: https://www.greenclimate.fund/project/fp022 (accessed on 15 January 2024).
  14. Henry, M.; Picard, N.; Trotta, C.; Manlay, R.; Valentini, R.; Bernoux, M.; Saint-André, L. Estimating Tree Biomass of Sub-Saharan African Forests: A Review of Available Allometric Equations. Silva Fenn. 2011, 45, 477–569. [Google Scholar] [CrossRef]
  15. Benzyane, M. Estimation de La Biomasse et Étude de La Croissance de l’arganier (Argania Spinosa L. Skeels) Dans Le Plateau de Haha (Essaouira); I.A.V. Hassan II: Rabat, Morocco, 1989. [Google Scholar]
  16. Belghazi, T. Typologie de l’arganeraie et Estimation de La Biomasse Aérienne Des Taillis d’arganier (Plateau Des Haha-Essaouira, Maroc). Ph.D. Thesis, UCL-Université Catholique de Louvain, Ottignies-Louvain-la-Neuve, Belgium, 2013. [Google Scholar]
  17. CCKP WB Climate Change Knowledge Portal. Available online: https://climateknowledgeportal.worldbank.org/country/ morocco/climate-data-historical (accessed on 9 December 2023).
  18. Lahou, O. Morocco—Global Yield Gap Atlas. Available online: https://www.yieldgap.org/Morocco (accessed on 9 December 2023).
  19. Nyamukuru, A.; Whitney, C.; Tabuti, J.R.S.; Esaete, J.; Low, M. Allometric Models for Aboveground Biomass Estimation of Small Trees and Shrubs in African Savanna Ecosystems. Trees For. People 2023, 11, 100377. [Google Scholar] [CrossRef]
  20. Allen, R.G. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements-FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998. [Google Scholar]
  21. Romero, F.M.B.; Jacovine, L.A.G.; Torres, C.M.M.E.; Ribeiro, S.C.; Da Rocha, S.J.S.S.; Novais, T.D.N.O.; Gaspar, R.D.O.; Da Silva, L.F.; Vidal, E.; Leite, H.G.; et al. Aboveground Biomass Allometric Models for Large Trees in Southwestern Amazonia. Trees For. People 2022, 9, 100317. [Google Scholar] [CrossRef]
  22. Kenzo, T.; Himmapan, W.; Yoneda, R.; Tedsorn, N.; Vacharangkura, T.; Hitsuma, G.; Noda, I. General Estimation Models for Above- and below-Ground Biomass of Teak (Tectona Grandis) Plantations in Thailand. For. Ecol. Manag. 2020, 457, 117701. [Google Scholar] [CrossRef]
  23. Zuur, A.F.; Ieno, E.N.; Elphick, C.S. A Protocol for Data Exploration to Avoid Common Statistical Problems. Methods Ecol. Evol. 2010, 1, 3–14. [Google Scholar] [CrossRef]
  24. Brown, S.; Gillespie, A.J.R.; Lugo, A.E. Biomass Estimation Methods for Tropical Forests with Applications to Forest Inventory Data. For. Sci. 1989, 35, 881–902. [Google Scholar] [CrossRef]
  25. Sprugel, D.G. Correcting for Bias in Log-Transformed Allometric Equations. Ecology 1983, 64, 209–210. [Google Scholar] [CrossRef]
  26. Mahmood, H.; Siddique, M.R.H.; Islam, S.M.Z.; Abdullah, S.M.R.; Matieu, H.; Iqbal, M.Z.; Akhter, M. Applicability of Semi-Destructive Method to Derive Allometric Model for Estimating Aboveground Biomass and Carbon Stock in the Hill Zone of Bangladesh. J. For. Res. 2020, 31, 1235–1245. [Google Scholar] [CrossRef]
  27. Saha, C.; Mahmood, H.; Nayan, S.N.S.; Siddique, M.R.H.; Abdullah, S.M.R.; Islam, S.M.Z.; Iqbal, M.Z.; Akhter, M. Allometric Biomass Models for the Most Abundant Fruit Tree Species of Bangladesh: A Non-Destructive Approach. Environ. Chall. 2021, 3, 100047. [Google Scholar] [CrossRef]
  28. Ispikoudis, S.; Zianis, D.; Tziolas, E.; Damianidis, C.; Rapti, D.; Tsiros, E.; Michalakis, D.; Karteris, A. Assessment of Forest Biomass and Carbon Storage in Habitat 9340 Quercus Ilex L. to Support Management Decisions for Climate Change Mitigation. Sustainability 2024, 16, 1403. [Google Scholar] [CrossRef]
  29. Nath, A.J.; Tiwari, B.K.; Sileshi, G.W.; Sahoo, U.K.; Brahma, B.; Deb, S.; Devi, N.B.; Das, A.K.; Reang, D.; Chaturvedi, S.S. Allometric Models for Estimation of Forest Biomass in North East India. Forests 2019, 10, 103. [Google Scholar] [CrossRef]
  30. Sileshi, G.W. A Critical Review of Forest Biomass Estimation Models, Common Mistakes and Corrective Measures. For. Ecol. Manag. 2014, 329, 237–254. [Google Scholar] [CrossRef]
  31. Akaike, H. A New Look at the Statistical Model Identification. IEEE Trans. Automat. Contr. 1974, 19, 716–723. [Google Scholar] [CrossRef]
  32. Durbin, J.; Watson, G.S. Testing for Serial Correlation in Least Squares Regression. II. Biometrika 1951, 38, 159. [Google Scholar] [CrossRef]
  33. Dijkstra, F.A.; Zhu, B.; Cheng, W. Root Effects on Soil Organic Carbon: A Double-Edged Sword. New Phytol. 2021, 230, 60–65. [Google Scholar] [CrossRef] [PubMed]
  34. He, H.; Zhang, C.; Zhao, X.; Fousseni, F.; Wang, J.; Dai, H.; Yang, S.; Zuo, Q. Allometric Biomass Equations for 12 Tree Species in Coniferous and Broadleaved Mixed Forests, Northeastern China. PLoS ONE 2018, 13, e0186226. [Google Scholar] [CrossRef] [PubMed]
  35. Mokany, K.; Raison, R.J.; Prokushkin, A.S. Critical Analysis of Root: Shoot Ratios in Terrestrial Biomes. Glob. Change Biol. 2006, 12, 84–96. [Google Scholar] [CrossRef]
  36. Xia, J.; Yuan, W.; Wang, Y.-P.; Zhang, Q. Adaptive Carbon Allocation by Plants Enhances the Terrestrial Carbon Sink. Sci. Rep. 2017, 7, 3341. [Google Scholar] [CrossRef] [PubMed]
  37. Thomas, S.C.; Martin, A.R. Carbon Content of Tree Tissues: A Synthesis. Forests 2012, 3, 332–352. [Google Scholar] [CrossRef]
  38. Ma, S.; He, F.; Tian, D.; Zou, D.; Yan, Z.; Yang, Y.; Zhou, T.; Huang, K.; Shen, H.; Fang, J. Variations and Determinants of Carbon Content in Plants: A Global Synthesis. Biogeosciences 2018, 15, 693–702. [Google Scholar] [CrossRef]
  39. Mensah, S.; Kakaï, R.G.; Seifert, T. Patterns of Biomass Allocation between Foliage and Woody Structure: The Effects of Tree Size and Specific Functional Traits. Ann. For. Res. 2016, 59, 49–60. [Google Scholar] [CrossRef]
  40. Kafuti, C.; Van den Bulcke, J.; Beeckman, H.; Van Acker, J.; Hubau, W.; De Mil, T.; Hatakiwe, H.; Djiofack, B.; Fayolle, A.; Loubota Panzou, G.J.; et al. Height-Diameter Allometric Equations of an Emergent Tree Species from the Congo Basin. For. Ecol. Manag. 2022, 504, 119822. [Google Scholar] [CrossRef]
  41. Payandeh, B. Choosing Regression Models for Biomass Prediction Equations. For. Chron. 1981, 57, 229–232. [Google Scholar] [CrossRef]
  42. Paul, K.I.; Roxburgh, S.H.; England, J.R.; Ritson, P.; Hobbs, T.; Brooksbank, K.; John Raison, R.; Larmour, J.S.; Murphy, S.; Norris, J.; et al. Development and Testing of Allometric Equations for Estimating Above-Ground Biomass of Mixed-Species Environmental Plantings. For. Ecol. Manag. 2013, 310, 483–494. [Google Scholar] [CrossRef]
  43. Nam, V.T.; van Kuijk, M.; Anten, N.P.R. Allometric Equations for Aboveground and Belowground Biomass Estimations in an Evergreen Forest in Vietnam. PLoS ONE 2016, 11, e0156827. [Google Scholar] [CrossRef]
  44. Bernal, B.; Murray, L.T.; Pearson, T.R.H. Global Carbon Dioxide Removal Rates from Forest Landscape Restoration Activities. Carbon Balance Manag. 2018, 13, 22. [Google Scholar] [CrossRef]
Figure 1. Location of the orchards studied in the Argan Biosphere Reserve in North Africa. Two experimental study plots were installed within each perimeter.
Figure 1. Location of the orchards studied in the Argan Biosphere Reserve in North Africa. Two experimental study plots were installed within each perimeter.
Sustainability 16 04592 g001
Figure 2. The applied experimental design for the sampled argan plants situated within the argan orchard.
Figure 2. The applied experimental design for the sampled argan plants situated within the argan orchard.
Sustainability 16 04592 g002
Figure 3. Relationship in logarithmic scale (base 10) between measured dry biomass (Biomass) and diameter (D) (a) and height (H) (b) for different components of young argan plants.
Figure 3. Relationship in logarithmic scale (base 10) between measured dry biomass (Biomass) and diameter (D) (a) and height (H) (b) for different components of young argan plants.
Sustainability 16 04592 g003
Figure 4. Mean carbon content within the components of young argan plants (2 to 6 years). Error bars refer to standard deviation.
Figure 4. Mean carbon content within the components of young argan plants (2 to 6 years). Error bars refer to standard deviation.
Sustainability 16 04592 g004
Figure 5. Pearson’s correlation for young argan plant components.
Figure 5. Pearson’s correlation for young argan plant components.
Sustainability 16 04592 g005
Figure 6. Standardized residuals against the fitted values plots, Q-Q plots, and predicted vs. observed values plots for leaf biomass (a), stem biomass (b), and root biomass (c), using the selected models.
Figure 6. Standardized residuals against the fitted values plots, Q-Q plots, and predicted vs. observed values plots for leaf biomass (a), stem biomass (b), and root biomass (c), using the selected models.
Sustainability 16 04592 g006
Figure 7. Mean biomass of leaves, stems, roots and total biomass across growing years. n, number of young plants sampled. The inset graph is a rescaling of the values for growing years with low biomass / tree value. Error bars refers to standard deviation.
Figure 7. Mean biomass of leaves, stems, roots and total biomass across growing years. n, number of young plants sampled. The inset graph is a rescaling of the values for growing years with low biomass / tree value. Error bars refers to standard deviation.
Sustainability 16 04592 g007
Figure 8. Total carbon stock and carbon dioxide (CO2) storage in the biomass of argan plants in the orchards studied. Error bars refer to standard deviation.
Figure 8. Total carbon stock and carbon dioxide (CO2) storage in the biomass of argan plants in the orchards studied. Error bars refer to standard deviation.
Sustainability 16 04592 g008
Table 1. Goodness-of-fit indicators and model-selection criteria.
Table 1. Goodness-of-fit indicators and model-selection criteria.
CriteriaDescriptionBest SelectionReference
Coefficient of determination (R2)Measures the proportion of the variance in the dependent variable (R2) that is explained by the independent variables in the model
  • Highest value
[4,19,30]
Adjusted coefficient of determination (R2)Adjusts R2 by acknowledging the number of independent variables in the model
  • Highest value
Residual standard deviation (RSE)Describes the difference between the standard deviations of the observed values compared with the predicted values
  • Lowest value
Root-mean-square error (RMSE)Measures the average distance between the predicted values and the observed values
  • Lowest value
Akaike information criteria (AIC)Minimizes issues associated with overfitting
  • Lowest value
[31]
Durbin–Watson (DW)Tests the autocorrelation of residuals in a linear regression model and reflect the independence of the residuals
  • 0 < DW < 2 indicates a positive correlation between successive residues
  • 2 < DW < 4 indicates a negative correlation between residues.
  • DW = 2 indicates no correlation
[32]
Mean squared error (MSE)Measure the average difference between the predicted and observed values in LOOCV data
  • Lowest value
[29]
Table 2. Descriptive statistics of young plants’ dry biomass components (kg), diameter (D), and total height (H), collected from argan orchards located within the Argan Biosphere Reserve in North Africa.
Table 2. Descriptive statistics of young plants’ dry biomass components (kg), diameter (D), and total height (H), collected from argan orchards located within the Argan Biosphere Reserve in North Africa.
Dry Biomass (kg/Plant)D (cm)H (m)
LeafStemRootTotalShoot-to-Root Ratio
Min.0.00080.0020.00060.0060.080.340.29
Max.1.022.051.424.482.677.911.42
Mean0.140.220.120.480.641.610.68
SD0.220.390.260.850.511.280.25
n *89898989898989
* Min, minimal; Max, maximal; SD, standard deviation; n, number of young plants sampled.
Table 3. Two-way ANOVA results on the effect of growing years and argan plant components on carbon content.
Table 3. Two-way ANOVA results on the effect of growing years and argan plant components on carbon content.
SourceType III Sum of SquaresdfMean SquareFSig.
Growing years1413.214353.3026.890.000
Plant components3383.9221691.96128.800.000
Growing years × plant components2365.678295.7122.510.000
Table 4. Regression models of allometric equation and associated coefficients used for the estimation of dry biomass in different components of argan plants including leaf, stem, root, and total biomass. n is the number of young plants studied.
Table 4. Regression models of allometric equation and associated coefficients used for the estimation of dry biomass in different components of argan plants including leaf, stem, root, and total biomass. n is the number of young plants studied.
NoteAllometric Equation ModelnCoefficients
abcde
Leaf1Ln b = a + b Ln (D)89−4.35 *2.89 *---
2Ln b = a + b Ln (D) + c Ln (H)89−3.57 *2.27 *1.44 *--
3Ln b = a + b Ln (D) + c Ln (H) + d Ln (age)89−6.95 *0.88 *1.21 *3.29 *-
4Ln b = a + b Ln (D × age)89−6.39 *2.04 *---
5Ln b = a + b Ln (D × H)89−3.26 *2.00 *---
6Ln b = a + b Ln (D × H × age)89−5.07 *1.56 *---
7Ln b = a + b Ln (D2 × H)89−3.71 *1.20 *---
8Ln b = a + b Ln (D2 × H × age)89−4.83 *1.02 *---
9Ln b = a + b Ln (D2 × H) + c ln (age)89−7.21 *0.60 *3.21 *--
10Ln b = a + b (D) × c (H)89−9.92 *2.84 *6.18 *−1.90 *-
11Ln b = a + b (D) × c (age)89−10.83 *2.79 *1.81 *−0.50 *-
12Ln b = a + b (D) × c (H) × d (age)89−10.08 *1.380.431.33 *-
Stem1Ln b = a + b Ln (D)89−3.40 *2.39 *---
2Ln b = a + b Ln (D) + c Ln (H)89−3.02 *2.08 *0.71 *--
3Ln b = a + b Ln (D) + c Ln (H) + d Ln (age)89−4.90 *1.30 *0.58--
4Ln b = a + b Ln (D × age)89−5.04 *1.65 *---
5Ln b = a + b Ln (D × H)89−2.50 *1.63 *---
6Ln b = a + b Ln (D × H × age)89−3.97 *1.26 *---
7Ln b = a + b Ln (D2 × H)89−2.87 *0.98 *---
8Ln b = a + b Ln (D2 × H × age)89−3.78 *0.83 *---
9Ln b = a + b Ln (D2 × H) + ln (age)89−4.88 *0.63 *1.84 *--
10Ln b = a + b (D × H)89−7.59 *2.36 *4.21 *−1.45 *-
11Ln b = a + b (D × age)89−8.15 *2.20 *1.25 *−0.35 *-
12Ln b = a + b (D × H × age)89−7.09 *1.11−0.170.63-
Root1Ln b = a + b Ln (D)89−4.06 *2.23 *---
2Ln b = a + b Ln (D) + c Ln (H)89−3.84 *2.06 *0.40 *--
3Ln b = a + b Ln (D) + c Ln (H) + d Ln (age)89−5.05 *1.56 *0.32--
4Ln b = a + b Ln (D × age)89−5.56 *1.52 *---
5Ln b = a + b Ln (D × H)89−3.22 *1.52 *---
6Ln b = a + b Ln (D × H × age)89−4.57 *1.16 *---
7Ln b = a + b Ln (D2 × H)89−3.56 *0.91 *---
8Ln b = a + b Ln (D2 × H × age)89−4.40 *0.77 *---
9Ln b = a + b Ln (D2 × H) + ln (age)89−4.89 *0.68 *1.22 *--
10Ln b = a + b (D × H)89−7.56 *2.01 *3.30 *--
11Ln b = a + b (D × age)89−7.95 *1.78 *0.97 *--
12Ln b = a + b (D × H × age)89−7.78 *3.93 *−3.08--
13Ln b = a + b Ln (D) + c Ln (H) + d Ln (age) + e Ln (root-to-shoot ratio)89−5.00 *1.48 *0.40 *1.38 *0.31 *
14Ln b = a + b Ln (D) + c Ln (H) + d Ln (ABG)89−2.81 *1.20 *0.010.41 *-
* Significant difference of 0.05.
Table 5. Statistical properties of regression models to estimate dry biomass in different components of argan plants including leaf, stem, root, and total biomass. Initials refer to statistical parameters mentioned in Table 1; CF stands for the correction factor.
Table 5. Statistical properties of regression models to estimate dry biomass in different components of argan plants including leaf, stem, root, and total biomass. Initials refer to statistical parameters mentioned in Table 1; CF stands for the correction factor.
NoteR2R2 AdjustedRSEAICRSMEDWCFMSE
Leaf10.850.850.85−26.314.2911.530.74
20.870.870.80−35.534.3021.490.67
30.950.950.50−120.24.3521.280.26
40.930.930.61−85.984.3321.360.38
50.870.870.81−35.704.3021.500.67
60.930.920.60−90.014.3321.340.36
70.870.870.80−37.184.3021.490.66
80.910.910.66−71.434.3321.390.45
90.950.950.51−116.454.3521.280.27
100.870.870.82−32.414.3021.500.69
110.940.940.54−106.194.3521.300.30
120.960.950.48−127.844.3521.260.23
Stem10.840.830.75−46.993.4221.460.59
20.840.840.74−48.313.4221.450.58
30.880.870.66−68.883.4421.390.47
40.870.870.66−69.633.4421.400.45
50.830.830.76−44.633.4221.460.60
60.870.860.67−67.073.4321.400.47
70.840.840.74−49.833.4221.450.57
80.870.860.68−65.273.4421.410.47
90.880.870.66−70.863.4411.390.45
100.830.820.78−39.263.4211.470.64
110.860.860.70−60.063.4421.410.51
120.880.860.68−59.363.4411.390.51
Root10.890.890.56−99.523.9521.320.33
20.890.890.56−99.423.9521.320.33
30.910.900.52−113.363.9621.290.28
40.900.900.52−114.613.9621.290.27
50.870.870.60−87.383.9521.350.37
60.900.900.54−107.613.9621.310.30
70.890.890.57−98.383.9521.330.33
80.900.900.52−112.873.9621.300.28
90.910.900.52−113.413.9621.300.28
100.860.860.63−77.583.9421.360.44
110.890.890.55−100.263.9521.320.34
120.910.910.52−108.493.9621.280.29
130.930.920.47−127.63.9721.260.24
140.920.910.49−124.043.9621.270.25
Total biomass 0.950.940.42−1482.7821.230.19
Table 6. Total biomass of leaves, stems, and roots of argan plants in the orchards studied.
Table 6. Total biomass of leaves, stems, and roots of argan plants in the orchards studied.
Growing YearsBiomass (t ha−1)
LeafStemRootTotal
20.001 ± 0.0010.003 ± 0.0020.002 ± 0.0020.006 ± 0.002
30.003 ± 0.0020.009 ± 0.0040.004 ± 0.0030.017 ± 0.003
40.013 ± 0.0110.032 ± 0.0240.01 ± 0.0090.055 ± 0.015
50.082 ± 0.0630.125 ± 0.0650.046 ± 0.0360.253 ± 0.055
60.201 ± 0.1060.185 ± 0.0230.083 ± 0.0510.469 ± 0.060
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Oumasst, A.; Tiouidji, F.E.; Tabi, S.; Zahidi, A.; El Mousadik, A.; El Finti, A.; Aitlhaj, A.; Hallam, J. Development of Allometric Equations to Determine the Biomass of Plant Components and the Total Storage of Carbon Dioxide in Young Mediterranean Argan Trees. Sustainability 2024, 16, 4592. https://doi.org/10.3390/su16114592

AMA Style

Oumasst A, Tiouidji FE, Tabi S, Zahidi A, El Mousadik A, El Finti A, Aitlhaj A, Hallam J. Development of Allometric Equations to Determine the Biomass of Plant Components and the Total Storage of Carbon Dioxide in Young Mediterranean Argan Trees. Sustainability. 2024; 16(11):4592. https://doi.org/10.3390/su16114592

Chicago/Turabian Style

Oumasst, Assma, Fatima Ezzahra Tiouidji, Salma Tabi, Abdelaziz Zahidi, Abdelhamid El Mousadik, Aissam El Finti, Abderrahmane Aitlhaj, and Jamal Hallam. 2024. "Development of Allometric Equations to Determine the Biomass of Plant Components and the Total Storage of Carbon Dioxide in Young Mediterranean Argan Trees" Sustainability 16, no. 11: 4592. https://doi.org/10.3390/su16114592

APA Style

Oumasst, A., Tiouidji, F. E., Tabi, S., Zahidi, A., El Mousadik, A., El Finti, A., Aitlhaj, A., & Hallam, J. (2024). Development of Allometric Equations to Determine the Biomass of Plant Components and the Total Storage of Carbon Dioxide in Young Mediterranean Argan Trees. Sustainability, 16(11), 4592. https://doi.org/10.3390/su16114592

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop