Skip to main content
Log in

Shade canopy density variables in cocoa and coffee agroforestry systems

  • Published:
Agroforestry Systems Aims and scope Submit manuscript

Abstract

An estimated 3.4 million hectares of cocoa and 9.7 million hectares of coffee are cultivated, globally, under shade trees, i.e. in agroforestry systems. Shade canopies are characterized in terms of tree density (N, trees ha−1), tree basal area (G, m2 ha−1) and percent canopy cover (%Cov). N, G and %Cov are named shade canopy density variables (SCDV). The use of these SCDV has two important limitations: (1) different combinations of values of the three SCDV variables generate very different shade tree stands (hence very different shading levels), and (2) Additional factors modify shading under shade canopies with constant SCDV values. This article uses the software ShadeMotion (www.shademotion.net) to show how 24 different, simple, even-sized, mono-layered, Cordia alliodora shade canopies with constant N, G and %Cov display significantly different shade levels and temporal patterns of shading depending on tree stem and crown diameter ratios, tree height, spatial planting configurations (square, random and alleys) and leaf fall patterns. A minimum set of variables capable of providing a more accurate description of the shading characteristics of a cocoa or coffee shade canopy is proposed. Our findings can shed light on the current debate on the pros and cons of the definitions of cocoa agroforestry used by chocolate and certification companies, governments, non-governmental organizations, and donors, especially in West and Central Africa. In this article, emphasis is given to cocoa, but the analysis, results and conclusions are equally applicable to coffee agroforestry systems.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6

Similar content being viewed by others

References

Download references

Acknowledgements

This research was funded by the University of Montpellier’s MAK’IT (Montpellier Advanced Knowledge Institute on Transitions) visiting scientist program, CATIE (Centro Agronómico Tropical de Investigación y Enseñanza), and CIRAD. Sergio Vílchez provided valuable orientation and help in the statistical analysis of the data. Two anonymous reviewers provided suggestions and comments that helped the authors to improve the quality of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

ES designed the study, ran the simulations with the software ShadeMotion, and wrote the various drafts of the manuscript. AS analyzed the bibliographic database and wrote the natural language processing algorithm. All authors carefully reviewed and commented on various versions of the manuscript and contributed to the interpretation of the results. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Eduardo Somarriba.

Ethics declarations

Conflict of interest

The authors have no relevant financial or non-financial interests to disclose, no competing interests, and no proprietary interests in any material discussed in this article.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary file1 (DOCX 275 KB)

Appendix 1: The mathematical relationships between N, G, %Cov, R, d, and k

Appendix 1: The mathematical relationships between N, G, %Cov, R, d, and k

With trees of same stem diameter (d), basal area per hectare (G) is obtained by multiplying the sectional stem area of one tree (g, m2 tree−1) by the population density represented in number of trees per hectare (N). Symbolically:

\({\text{G}} = {\text{g}}*{\text{N}}\)

Assuming tree stems are cylindrical, the sectional area can be estimated by the area of a circle of diameter d (in meters i.e. d in cm divided by 100),

$${\text{g }} = \, \left( {{\text{d}}/{1}00} \right)^{{2}} *\left( {\uppi /{4}} \right) \, = \, \left( {\uppi *{\text{d}}^{{2}} } \right)/{4}0000$$

Hence

$${\text{G }} = \, \left( {{\text{N}}*\uppi *{\text{d}}^{{2}} } \right)/{4}0000$$
(1)

Re-arranging Eq. 1

$${\text{N }} = \, \left( {{\text{G}}*{4}0000} \right)/\left( {\uppi *{\text{d}}^{{2}} } \right)$$
(2)

And

$${\text{d }} = \, \left[ {\left( {{\text{G}}*{4}0000} \right)/{\text{N}}*\uppi )} \right]^{\raise.5ex\hbox{$\scriptstyle 1$}\kern-.1em/ \kern-.15em\lower.25ex\hbox{$\scriptstyle 2$} }$$
(3)

Considering that

$$\begin{aligned} {\text{R }} = & {\text{ k}}/{\text{d}} \\ {\text{K }} = & {\text{ R}}*{\text{d}} \\ {\text{d }} = & {\text{ k}}/{\text{R}} \\ \end{aligned}$$

For a given R

$${\text{G }} = \, \left( {{\text{N}}*{\text{k}}^{{2}} *\uppi } \right)/\left( {{4}0000*{\text{R}}^{{2}} } \right)$$
(4)
$${\text{N }} = \, \left( {{4}0000*{\text{G}}*{\text{R}}^{{2}} } \right)/\left( {\uppi *{\text{k}}^{{2}} } \right)$$
(5)
$${\text{K }} = \, \left[ {\left( {{4}0000*{\text{G}}*{\text{R}}^{{2}} } \right)/\left( {{\text{N}}*\uppi } \right)} \right]^{\raise.5ex\hbox{$\scriptstyle 1$}\kern-.1em/ \kern-.15em\lower.25ex\hbox{$\scriptstyle 2$} }$$
(6)

Percent canopy cover (%Cov) is obtained by multiplying the opacity-adjusted crown projection area per tree (z) by the number of trees per hectare, divided by 10,000 m2 in one hectare and later multiply it by 100 to express the ratio in percent. Assuming that vertical crown projection is circular,

$${\text{z }} = {\text{ k}}^{{2}} *\left( {\uppi /{4}} \right)*{\text{p }} = \, \left( {{\text{R}}*{\text{d}}} \right)^{{2}} *\left( {\uppi /{4}} \right)*{\text{p }} = {\text{ R}}^{{2}} *{\text{d}}^{{2}} *\left( {\uppi /{4}} \right)*{\text{p}}$$

Hence

$$\% {\text{Cov }} = { 1}00*\left( {{\text{N}}*{\text{z}}*{\text{p}}} \right)/{1}0000 \, = \, \left( {{\text{N}}*{\text{R}}^{{2}} *{\text{d}}^{{2}} *{\text{p}}*\uppi } \right)/{4}00$$
(7)

Substituting N by Eq. (1) and rearranging terms,

$$\% {\text{Cov }} = { 1}00*{\text{G}}*{\text{p}}*{\text{R}}^{{2}}$$
(8)

Also

$$\% {\text{Cov }} = \, ({\text{N}}*{\text{p}}*\uppi *{\text{k}}^{{2}} )/{4}00$$
(9)

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Somarriba, E., Saj, S., Orozco-Aguilar, L. et al. Shade canopy density variables in cocoa and coffee agroforestry systems. Agroforest Syst 98, 585–601 (2024). https://doi.org/10.1007/s10457-023-00931-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10457-023-00931-2

Keywords

Navigation