Introduction
Cocoa (Theobroma cacao L.) is a major cash crop for millions of farmers, particularly in developing countries. Indonesia is the third-largest cocoa-producing country in the world, after Ghana and Ivory Coast, with about 50 to 75% of its output coming from Sulawesi. From 1990 to 2010, the land allocated for cocoa production in Indonesia increased 10-fold, reaching a plateau of approximately 17 000 km2 in 2010 (FAOSTAT, 2020). However, the gross production of cocoa has decreased since 2010 because of declining yields per hectare that started in the early 1990s (FAOSTAT, Reference Abdulai, Hoffmann, Jassogne, Asare, Graefe, Tao, Muilerman, Vaast, Van Asten, Läderach and Rötter2020). Because of the limited availability of land, productivity per hectare will need to increase if Indonesia is to contribute to the increasing global demand for cocoa.
The decline of productivity of cocoa plantations has been related to inadequate management leading to problems such as lower soil fertility and increased levels of pests and disease (Asare et al., Reference Asare, Afari-Sefa, Muilerman and Anim-Kwapong2018). Due to a lack of resources to improve the fertility and health of existing farms, low profits have often pushed farmers to move cultivation to forested areas (Ruf, Reference Ruf2001). To prevent deforestation, many public and private stakeholders in the cocoa industry are seeking to improve crop productivity on existing farms (Carodenuto, Reference Carodenuto2019; Weber, Reference Weber2017). The revitalization of cocoa farms is being implemented through improved planting material, the integrated control of pests and diseases, optimal shade control, and the long-term management of soil fertility (Asare et al., Reference Asare, Afari-Sefa, Muilerman and Anim-Kwapong2018; Wood & Lass, Reference Wood and Lass2008).
Soil fertility often deteriorates on cocoa farms (Adeniyi et al., Reference Adeniyi, de Clercq and van Niekerk2017). Typically, cocoa trees planted on a freshly cleared forest initially benefit from high fertility due to high organic matter levels and well-developed soil structure. However, the subsequent removal of the harvested pods and beans can reduce nutrient levels (Boyer, Reference Boyer1973; Fassbender et al., Reference Fassbender, Alpízar, Heuveldop, Fölster and Enríquez1988; Hartemink, Reference Hartemink2005; Thong & Ng, Reference Thong and Ng1978; van Vliet et al., Reference van Vliet, Slingerland and Giller2015), and soil fertility declines if they are not replenished with organic or mineral/inorganic fertilizers (Aikpokpodion, Reference Aikpokpodion2010; Hartemink, Reference Hartemink2005). Degradation can also occur due to the acidification of the soil from the use of acidifying fertilizers, like urea, organic matter decay, the removal and leaching of basic cations (Goulding, Reference Goulding2016), and an increased availability of toxic elements (Lal et al., Reference Lal, Hall and Miller1989).
To maintain soil fertility, farmers typically apply amendments and fertilizers to replenish nutrients stocks and correct soil acidity. Few peer-reviewed studies have evaluated their effects both on cocoa productivity and soil properties, and where there is research, it often focuses on short-term effects on seedlings or young cocoa trees (Ahenkorah et al., Reference Ahenkorah, Halm, Appiah, Akrofi and Yirenkyi1987; van Vliet et al., Reference van Vliet, Slingerland and Giller2015; Wessel, Reference Wessel1971). Verlière (Reference Verlière1981) reported that only a few fertilizer experiments with cocoa had provided significant results, and there was a need to determine the interactions between shade management, cocoa nutritional needs, and productivity. Fewer fertilizer studies have been conducted in South-East Asia than in West Africa, and the work completed in Sulawesi is scarce (Mulia et al., Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019). In addition, there have been few studies examining the combined effects of fertilizers, organic inputs (e.g., compost or manure), and other amendments (e.g., lime dolomite).
To address this lack of information for Indonesia, Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019) reported the responses for the first 4 years of a cocoa field trial established in 2011 with various combinations of fertilizers and amendments (organic and inorganic). Because the treatments were continued after the first 4 years, this study aims to evaluate the effects of the treatments on soil properties, growth, and yield on more mature cocoa plants in order to develop fertilizer recommendations based on organic amendments.
Methodology
Experimental site
The experimental area has been described in detail by Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019), but the main points are repeated here for clarity. The plot is located in Bone-Bone, South Sulawesi, Indonesia (2.605833°S, 120.612333°E). The principal activities carried out between 2011 and 2018 on the experiment are presented in Table S1, Table S2, and Table S3. Cocoa plants of clone PBC123, also known as Sulawesi 1, were planted at a 3 m x 3 m spacing in December 2011, that is, 1111 trees per hectare. Shade trees of Gliricidia sepium were planted at the same time as the cocoa trees, as well as existing coconuts. In January 2016, after the completion of Mulia’s study, new cocoa trees were planted to replace the ones that died during the first phase of the experiment. However, after this replanting, productivity was not recorded on those trees. Monthly precipitation for the 2011–2018 period (Table S4 in supplementary material) was obtained from Mars’s Cocoa Research Station, located in Tarengge, approximately 20 km from the site (South Sulawesi, Indonesia).
Treatments
Initially, the experiment followed a randomized block design repeated four times, with 16 cocoa trees for each repetition (Table 1). As described in Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019), there was a control treatment with no fertilizers or amendments (Treatment A), and treatments comprising either the application of NPK fertilizer + urea (Treatment B), compost (Treatment C), and dolomite (Treatment D). Then four possible combinations of the three primary inputs were implemented with cumulated application rates (Treatments E to H). There were eight treatments in total, coded from A to H (one control + three with individual inputs + four with combinations of inputs). The treatments described in Table 1 were split into two applications applied at a 6-month interval. The inorganic fertilizer applied to each tree was a mix of 374 g of Phonska (15% N, 15% P2O5, and 15% K2O with traces of S) and 250 g of urea. The compost was locally made of 60% cow manure, 15% empty oil palm bunches, 10% rice straw, 10% diverse leaves (banana, grass, Gliricidia, and maize), and 5% cocoa pod husks. A microorganism mix (EM4) was also applied. The compost treatments comprised the application of 5 kg of compost to each tree twice a year (10 kg per year, distributed in six small pits located at 1 m from the trunk, Table S2), and the full chemical composition is described by Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019). The dolomite amendment comprised 18–22% of MgO, but the content for other constituents such as carbonates or calcium oxide CaO was unknown, but typical contents are around 22% Ca, 13% Mg, 13% C, and 52% O in elemental terms (Mineralogical Society of America, 2003).
Table 1. Breakdown of the treatments applied between 2012 and 2018 (adapted from Mulia et al., Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019)
Soil amendments and fertilizers were applied twice per year and per tree to provide total quantities as follows: 374 g NPK (‘Phonska’) and 250 g urea (mineral fertilizer), 5 kg dolomite, and 10 kg compost. Combinations (Treatments E–H) were additive. The values represent the total quantities of elements (g) provided per tree each year in each treatment. Phonska is a subsidized compound fertilizer made from three raw materials: urea, DAP (diammonium phosphate), rock phosphate, MOP (potassium chloride), and ‘other macronutrients’ according to the manufacturer (https://www.pupukkaltim.com/en/distribution-product-product-knowledge). At the time of planting, mineral fertilizer, 100 g NPK (Phonska) and 150 g triple superphosphate (36%), was added to each tree in equal amounts to provide adequate and uniform nutrient conditions for the establishment of all plants in the first few months after planting out (Mulia et al., Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019).
Sampling and analyses
Soil
In each plot, one soil sample was collected, at the center of the plot, 1 m from a cocoa trunk, at a depth of 0–20 cm, below the scraped surface litter, using an auger. The same day, samples were air-dried before being ground and sieved to 2 mm. The core ring method was used to collect bulk density samples on each plot, next to the soil sample, at a 0–5 cm depth. These were later air-dried at 60°C for 48 hours before weighing. In total, 32 soil samples were collected for soil analyses and bulk density measurement, corresponding to eight treatments with four repetitions each. The samples were then split with one sample sent to Asian-Agri Laboratory in North Sumatra, Indonesia, and the second sample was sent to Cranfield University in the United Kingdom. The analytical methods corresponding to each soil property are listed in Table 2 (Asian-Agri analyses) and Table 3 (Cranfield University analyses). A difference must be noted between the term ‘extractable’, referring to Asian-Agri analyses (using HCl at 25% v/v as an extractant), and the term ‘total’ which refers to Cranfield analyses (using microwave-assisted aqua regia digestion).
Table 2. Soil bulk density, pH, organic carbon, and nitrogen content and extractable and exchangeable nutrient contents in December 2018, 6.5 years after Treatments A–H were first applied
Soil treatments were applied twice per year, beginning in May 2012. Soil properties were determined in the soil laboratory of the Asian-Agri Laboratory in North Sumatra, Indonesia (samples collected in December 2018). Soil treatments are: A, control; B, mineral; C, compost; D, dolomite; E, mineral/compost; F; mineral/dolomite; G, compost/dolomite; H, all amendments. Means are calculated on four samples and given with three significant figures. Numbers in brackets are standard deviations. Treatments with the same letter are not statistically different (P > 0.05). K–W attached to a statistical significance refers to the p-value of Kruskal–Wallis test, used if ANOVA assumptions were not met. Methods: ♠After Fahmy (1977); ○Core ring method; ∗pH determined by AIAT Soil Laboratories, Maros; †Walkley–Black method; ‡Kjeldahl method; §25% HCl extraction; ¶Bray-I method; ††Ammonium acetate (pH 7) extraction; and ‡‡KCl (1 N) extraction. ‘na’ stands for not applicable (calculated data). ¥ B.S. refers to base saturation.
Table 3. Total contents of some soil elements of the experiment in December 2018
Soil treatments were applied twice per year, beginning in May 2012. These soil properties were determined at Cranfield University laboratory (samples collected in Indonesia in December 2018). Soil treatments are: A, control; B, mineral; C, compost; D, dolomite; E, mineral/compost; F; mineral/dolomite; G, compost/dolomite; and H, all amendments. Means are calculated on four samples and given with three significant figures. Numbers in brackets are standard deviations. Treatments with the same letter are not statistically different (P > 0.05). K–W attached to a statistical significance refers to the p-value of Kruskal–Wallis test, used if ANOVA assumptions were not met. ‡ The total soil contents of these elements were determined after dry combustion (SOP based on ISO 10694:1995), while the total contents marked by ♠ were analyzed by ICP-MS after using an HCl/HNO3 extractant (SOP based on ISO 11047:1998).
Growth of cocoa
Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019) reported the height of each cocoa tree height. In this study, the cocoa growth was determined from the circumference of the trunks (30 cm from the soil surface) because the bushes had recently been pruned. These measurements were converted to mean basal areas per treatment (formula and results available in Table S5). Both the soil sampling and the trunk measuring occurred at the same time in December 2018. The number of dead cocoa trees was used to calculate average survival rates per treatment and is presented in the supplementary material (Table S6).
Productivity and bean quality
From January 2015 to December 2018, cocoa tree productivity was assessed by counting the total number of pods produced per tree over a year, and measuring the weight of pods and the annual yield of cocoa beans per plot, as described in Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019). A pod index (PI, number of pods required to produce 1 kg of dry beans) for each treatment was derived from the pod counts and the mean pod weights (PI = 1000 ÷ average yield per pod in grams) (Table S7). The number of healthy and infected pods was also recorded. The pods were categorized as being uninfested or infested following the method previously described in McMahon et al. (Reference McMahon, Purung, Lambert, Mulia, Susilo, Sulistyowati, Sukamto, Israel, Saftar, Amir, Purwantara, Iswanto, Guest and Keane2015). A sample of the harvested beans was also collected in November 2017 to determine the waste fraction and the average weight of the fresh beans. Harvest quality results were obtained from Mars Laboratory in Makassar. The mean pod count per tree, yield per pod, pod index, and proportion of infected pods for each year and treatment are presented as supplementary material (Table S7). In Mulia’s study, productivity data covered the second semester of 2014 and the first of 2015, whereas this study recorded the annual yields from January 2015 to December 2018.
Statistical analyses
Basal areas, dry bean yields, and pod counts were averaged per surviving tree (≤ 16) and then extrapolated to 1000 trees to provide a per hectare value (as in Mulia et al., Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019). Additionally, the average yield per planted tree was also calculated, this time dividing the yield per plot by the number of initial trees, 16, and then multiplying by 1000 to derive the adjusted yield per hectare. Averaging per surviving tree minimizes treatment differences due to mortality rates. However, it can lead to the confounding effect in that trees with lower competition (i.e., where the mortality is high) could show greater growth rates and higher yields on a per tree basis.
Means and standard deviations of the four replicates were calculated for each variable (growth, productivity, and the soil analyses) and each treatment. Statistical differences between the treatments were evaluated by submitting the individual observations to an ANOVA followed by a Tukey HSD test in R 3.6.0 using the package ‘agricolae’ (at P = 0.05). If ANOVA assumptions were not met, the Kruskal–Wallis test was applied using the ‘kruskal’ function (at P = 0.05), and multiple comparisons between the treatments were obtained by using Fischer’s least significant difference post hoc test (with a level of significance at 0.05). Because of the lack of dispersion data in 2014, we used Welch one-sample t-tests to estimate the statistical difference between the two soil sampling periods, 2014 and 2018. In this case, we had to assume that the 2014 result was equal to the true mean (Table S8).
Results
Tree basal areas and survival rates
Tree basal areas and survival rates in December 2018, 7 years after planting, are presented in Figure 1a. Compared to the control (Treatment A) and the mineral fertilizer-only plots (Treatment B), we found significantly higher mean basal areas where compost and dolomite together (Treatment G) and the full combination were applied (Treatment H). The observed basal areas (of surviving trees only) were relatively heterogeneous for all treatments, with coefficients of variation ranging from 25 to 41% (Table S5). Treatment B had the lowest survival rate (41%), followed by the control (67%) and dolomite Treatment D (72%). All the other treatments had a survival rate higher than 80%.
Figure 1. Effect of soil amendments on (a) mean basal area per cocoa tree and survival rates in December 2018 and (b) annual yield of dry beans in each year from 2015 to 2018. In December 2018, cocoa trees were 7.4 years old (89 months). For each variable, means correspond to an average of the four replicates. Error bars represent standard errors of the mean based on surviving trees (light bars). Upper case letters on the X-axis refers to each treatment, namely: A is control; B is mineral fertilizer; C is compost; D is dolomite; E is mineral fertilizer and compost; F is mineral fertilizer and dolomite; G is compost and dolomite; and H is the full combination. Treatments with the same lower-case letter are not statistically significant (P > 0.05). For (b), the dark bars represent the harvest averaged for 16 initially planted trees, while the addition of light bars shows the average yield calculated per surviving tree (i.e., excluding dead trees, excluding mortality). From 2016, one replicate of Treatment B was excluded from the average calculation, as all trees were dead, leaving three replicates instead of four.
Yields and harvest quality
Dry bean yield
Dry bean yields for each treatment from 2015 to 2018 are shown in Figure 1b and detailed in Table S9. Productivity was very low in 2015, with a pattern similar to Mulia et al.’s results for 2014–2015 (low for A, B, D, and F; high for C, E, G, and H). In 2015, the treatments receiving no compost (Treatments A, B, D, and F) provided the lowest annual yields of only 110 kg ha−1 on average. In 2016, the mean level of yields increased, but the treatments without compost (A, B, D, and F) again produced less than the composted treatments (C, E, G, and H). The following year, differences between treatments were less marked. Finally, in 2018, only the yield in the mineral fertilizer treatment (B) was greater than the yield in the control treatment (A); the other treatments’ yields were similar. The coefficients of variation associated with a given year declined from 69% in 2015 to 27% in 2018 (Table S9). Over the 4 years, the compost-only treatment (C), the fertilizer and compost treatment (E), and the compost and dolomite treatment (G) produced the highest cumulative yields (2600–2750 kg ha−1; Table S9).
Pod count and yield per pod
The pattern for pod counts in 2015 and 2016 (Table S7) was similar to dry bean yields: treatments with compost (C, E, G, and H) had high pod counts (27–47), the treatments without compost (A, B, D, and F) had low counts (5–16). From 2017, the pod counts between treatments seemed to converge in a similar way as the yields.
The mean yield of dry beans per harvested pod was highly variable across treatments over the years (Table S7), ranging from a minimum of 10.8 g attained by Treatment D in 2015 to a maximum of 30.3 g (for pod-producing replicates) attained by the mineral fertilizer Treatment B in 2017. The low yields obtained in 2015 were associated with the smallest mean yield per pod of 12.4 g across all the treatments. Mean yields of dry beans per pod increased to 23.2 g in 2016 and 26.8 g in 2017, before declining to 19.1 g in 2018.
Yield index
To partially account for the differences in competition caused by the replanting in 2016, a yield index was calculated for each treatment (Table S10). It was determined as the ratio of the yield of dry beans per hectare (kg) divided by the basal area (cm²). On average, the mineral fertilizer treatment had the highest yield index of around 10 g cm−2, while the value for the other treatments ranged from 3.8 to 5.6 g cm−2.
Pod index
The dolomite treatment attained the highest pod index in 2015 (Table S10), reaching 96 pods kg−1, equating to a very small mean pod size of 10.8 g. The lowest pod index of 33 pods kg−1 was attained by the mineral fertilizer treatment (B) in 2017 when the mean pod weight reached 30.3 g. In 2016 and 2017, the pod indices were uniform across treatments with a mean value of 44 and 38 pods kg−1. Even though yields were more uniform in 2018, pod indices differed, ranging from 43 pods kg−1 (Treatment B) to 70 pods kg−1 (Treatment G).
Proportion of infected pods, dry bean weight, and fraction of waste beans
The proportion of infected pods was high, ranging from 0.62 to 0.97 for individual treatments per year (Table S7), and there was no consistent treatment effect. The analysis of a production sample, collected in November 2017, revealed differences between average dry bean weights (Table S11). The smallest beans were obtained for Treatment B (1.25 g), which were lower than the control (1.30 g). The largest beans were found for Treatment D (1.59 g). Compared to the mineral fertilizer Treatment B, Treatments C to H produced beans that were, on average, 18% heavier, while those in the control were only 65% of the weight of those in the mineral fertilizer treatment. The fraction of waste beans ranged from 8% (Treatment E) to 18% (Treatment F), with the others between 10 and 14%.
Soil properties
The results of the soil analyses are presented in Table 2 and Table 3, and the soil changes that occurred in 4 years are shown in Figure 2 (arranged in terms of soil property) and Figure S1 (arranged in terms of the treatments). The mean annual rates of change for each element in kg ha−1 year−1 are presented in Table S12.
Figure 2. Summary of the differences observed in the soil properties between 2014 and 2018. Each difference was calculated by taking 2014 as the initial value and 2018 as the final one. Stars rating correspond to the following rule, calculated after a Welch one-sample t-test: P ≤ 0.001, ***; P ≤ 0.01, **; P ≤ 0.05, *. Bars with no stars indicate no statistical difference between the 2 years (p = 0.05). Please refer to Table S8 in the supplementary material for the exact p-values). For extractable Ca and Mg, F had the largest effect, but the difference is not significant, only because the variability was very high. The mean annual rate of change for each element in kg ha−1 year−1 is shown in Table S12.
Bulk density and pH
The mean surface bulk density of the soil was statistically similar for each treatment, with a mean value of 1.09 g cm−3. For all treatments, soil pH values increased between 2014 and 2018 (Figure 2). Treatments A, B, and C, which did not receive dolomite, were the most acidic in 2018 with a value of 5.21–5.36 (Table 2). The pH was about neutral (6.75–6.87) for D, F, and G, where dolomite was applied. The pH of E and H were intermediate (respectively 5.79 and 6.25).
Carbon (C)
Despite the addition of compost, the organic C contents across the treatments in 2018 were statistically similar, with a mean value of 1.25% (Table 2). The organic C content in Treatments C, G, and H appeared to have significantly declined since 2014 (Figure 2). Subsequent analyses of the total C after dry combustion (following ISO 10694:1995, in Table 3) determined with an Elementar Analyzer at Cranfield University on sample duplicates also showed no statistical difference between treatments. However, the mean total carbon content of 1.65% was about 30% higher than the organic carbon values. Possible reasons for the higher reading could be the presence of inorganic forms of C (included in the Cranfield measurement) or a systematic difference between the Elementar and Walkley–Black methods.
Nitrogen (N)
Despite the different treatments, there was no statistical difference between treatments in the level of soil nitrogen, determined within either the Kjeldahl method (Table 2) or the Elementar Analyzer (Table 3). The nitrogen values reported in 2018 of 0.123–0.153% were broadly similar to measurements made in 2014 and reported by Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019). In 2018, the C:N ratios based on the measurements with the Asian-Agri’s dataset were approximately around 9:1. The C:N ratio determined using Cranfield University’s data also indicated similar results between treatments, but the higher carbon measurement resulted in a higher mean C:N ratio of 12:1.
Phosphorus (P)
In 2018, the extractable P content in the control Treatment A (18.4 ppm) was lower than that (68.4 ppm) in Treatment E receiving mineral fertilizer and compost (Table 2). The other treatments were not statistically different. The total P contents also followed the same pattern, with the value (283 ppm) for Treatment A being less than that (366 ppm) for Treatment E; the rest were statistically similar (Table 3). The extractable P contents measured in 2018 (18–68 ppm) were substantially lower than those measured in 2014 (227–471 ppm; Figure 2).
In 2018, the mean value for available P (Bray 1) for each treatment ranged from 1.9 to 12.6 ppm (Table 2). These values were lower than those (10–40 ppm) in 2014 and those (13 ppm) in 2011 prior to planting. The fraction of available P relative to total P ranged from 5% (Treatment D) to 25% (Treatment C). The higher pH observed for Treatment D was not associated with a noticeable increase in P availability (lowest available P at 1.9 ppm), nor did the addition of compost affect P availability (higher for Treatments C and E, 10.3 ppm on average, but not G and H, 4.2 ppm on average).
Potassium (K)
There were no statistical differences for exchangeable, extractable, and total K between the treatments in 2018 (Table 2). However, the measured exchangeable K decreased in all treatments from 2014 to 2018 (Figure 2). In 2018, the mean value across all treatments was below adequate levels (between 117 and 235 ppm, according to Nelson et al., Reference Nelson, Webb, Berthelsen, Curry, Yinil, Fidelis, Fisher, Oberthur and Oberthür2011). By contrast, between 2014 and 2018, extractable K increased with all the treatments, except in Treatment B, where it decreased, and in the control for which the change was negative but not significant (Figure 2).
Calcium (Ca)
In 2018, there were treatment differences in the level of extractable and total Ca (Table 2 and Table 3). Treatments D and F, which received dolomite, had the highest concentrations. The lowest levels were found in the control (Treatment A) and the plots receiving only mineral fertilizer (Treatment B). Extractable Ca increased significantly in the limed treatments, with a peak at 2 g kg−1 for Treatment F. The treatment effects on exchangeable and total Ca were similar to those for extractable Ca. The highest value occurred where dolomite was applied. Again, there was potentially an outlier in Treatment F, where a reading reached 5.7 g kg−1, perhaps due to the presence of incompletely dissolved dolomite. After 2014, the Ca contents of all the treatments which did not receive dolomite decreased to low levels.
Magnesium (Mg)
In 2018, there was no treatment effect on the level of extractable Mg (Table 2), but the exchangeable (Table 2) and total Mg (Table 3) showed a similar treatment response to Ca. The highest contents corresponded to the treatments with dolomite application; the lowest values were observed for the control and mineral fertilizer-only plots. There was again a peak for Treatment F, possibly due to an outlier, while the other observations were closer to 900 ppm. The lowest values were found for Treatments A and B, slightly exceeded by Treatments C and E. For exchangeable Mg, a pattern similar to Ca was found: highest where dolomite was applied, peaks for Treatments F and D and minimums for the control and the mineral fertilizer Treatments A and B. As for Ca, one very high Mg measurement was found for one of the samples, suggesting the presence of high concentrations of dolomite. As with Ca, the Mg contents have decreased since 2014 in the plots without dolomite, while it increased where it was applied.
Sodium (Na)
In 2018, no statistical differences were found between the treatments for exchangeable Na contents (Table 2). However, between 2014 and 2018, the Na contents decreased for all treatments (Figure 2). Sodium adsorption ratios were all very low (<1), suggesting no degradation risks of soil structure (USDA NRCS, 2017).
Aluminum (Al)
Distinctively, the lowest exchangeable Al contents were recorded in the treatment with the application of dolomite, which also decreased compared to 2014. Overall, total Al contents were statistically the same for all treatments, around 2%.
CEC and base saturation
The treatments did not appear to influence the cation exchange capacity (CEC), which was low in all the treatments. The treatments showed distinct base saturations. Treatment F attained 125% of saturation, possibly due to very high Ca and Mg found in the samples, themselves probably coming from the dolomite.
Discussion
Recontextualization
The aim of extending the study by Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019) was to determine if the yield responses of 7-year-old cocoa trees to different soil amendments were similar to those reported for 3-year-old trees and if the differences in responses could be related to soil chemical properties. This information was then used to develop fertilizer recommendations.
The initial study highlighted the beneficial effects of compost in increasing cocoa trees’ height, increasing flowering and yields, and pod quality. There were no obvious treatment effects on leaf nutrient contents and pest and disease incidence. The initial study also highlighted that mineral fertilizer and dolomite application to the young cocoa trees were associated with high mortality rates. Soil-wise, it was suggested that the high yields in the compost treatments were associated with increased nutrient availability and uptake. Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019) also suggested that the mineral fertilizer application could have resulted in nutritional deficiencies in Ca and Mg. This discussion examines how an additional 4 years of measurements improves our understanding of the medium-term effects of the different soil amendments.
Soil response to treatments
The effect of the treatments on soil properties is considered in terms of the effect of the compost, the mineral fertilizer, and the dolomite.
Effect of the compost
Despite the addition of almost 1 Mg of organic carbon (OC) per hectare per year from 2012 to 2018, compost inputs did not result in significantly higher surface bulk density, CEC, and soil C contents than the control. Several factors could explain this.
First, although the mean value of OC in the treatments ranged from 1.17% to 1.40%, these differences were not statistically different. A difference of 0.33% OC over one hectare to a depth of 20 cm and assuming a bulk density of 1.09 g cm−3 is equivalent to 7 Mg per hectare. Hence, the level of replication (four per treatment) described in this paper would be insufficient to identify a statistically significant difference from the addition of approximately 6 Mg of OC per hectare over 6 years. Such analysis highlights the very high levels of replication needed to detect soil changes because of the innate spatial variability in soil properties (Upson et al., Reference Upson, Burgess and Morison2016).
A second factor relates to the mode of compost application. Compost was applied in six small pits, surrounding the cocoa trees (Table S2). However, even though we sampled soil at the same distance from the trunk, the location of our samples was randomly positioned. Moreover, since we only collected one sample per treatment, it is easy to miss one of those little pits and therefore miss the local effect of compost on those spots.
A third factor is that the compost had a low C:N ratio (7:1), suggesting that it was prone to rapid mineralization instead of being turned into more ‘passive’ and recalcitrant forms of soil organic matter through humification (Nicolardot et al., Reference Nicolardot, Recous and Mary2001; Tian et al., Reference Tian, Kang and Brussaard1992; Weil & Brady, Reference Weil and Brady2017). By contrast, the use of more recalcitrant organic matter (e.g., higher C:N or lignin:N ratios) could be a more promising method to raise soil organic matter levels (Talbot & Treseder, Reference Talbot and Treseder2012).
Various authors have argued that a key benefit of compost addition is the increase in soil water holding capacity (Adugna, Reference Adugna2016; Blanco-Canqui & Lal, Reference Blanco-Canqui and Lal2007; Nguyen, Reference Nguyen2013; Smith, Reference Smith2018; Zemánek, Reference Zemánek2011). This is useful as cocoa yield is strongly affected by soil water regimes (Abdulai et al., Reference Abdulai, Hoffmann, Jassogne, Asare, Graefe, Tao, Muilerman, Vaast, Van Asten, Läderach and Rötter2020; Dada, Reference Dada2018; Kotei, Reference Kotei2019). Improving soil water storage, it may have served as a water reserve which may have resulted in different levels of drought stress between the composted and non-composted treatments.
Effect of dolomite
Over the 4 years, the dolomite (without compost) plots received about 4950 kg Ca ha−1 and 2900 kg Mg ha−1; the changes in measured soil contents were equivalent to 1820–4300 kg Ca ha−1 and 2006–3074 kg Mg ha−1. Seven years after the experiment’s start, the pH in the top 20 cm of the four treatments receiving dolomite was 6.25–6.83, compared to 5.21–5.79 in the treatments receiving no dolomite. Although the addition of dolomite alone increased yields relative to the control (+40%; presumably because of the higher pH increasing nutrient availability), the increase was less than that for the addition of compost alone (+195%) and the addition of mineral fertilizer alone (+118%; Figure 1b). The effect of dolomite on soil pH may also have contributed to reduce the availability of Al (Figure 2).
Effect of mineral fertilizer
In the initial study, mineral fertilizer alone did not have a beneficial effect on yield, but 7 years after planting, the mineral fertilizer treatment was the only treatment producing a statistically higher yield than the control plot. However, if yields including mortality rates are considered (dark bars in Figure 1), this improved performance is no longer apparent. Mineral fertilizer addition to young tree crops can be problematic, and high rates can lead to toxicities and root scorching (Gauthier et al., Reference Gauthier, Kaiser, Dockery and Goodin2014). This is less of a problem once a tree has established a robust root system. Also, regularly harvesting cocoa pods removes nitrogen, phosphorus, and potassium and can create deficiencies for these nutrients. For example, a yield of approximately 1000 kg ha−1 described in Treatment B would represent the removal of approximately 40 kg of nitrogen, 6 kg of phosphorus, and 62 kg of potassium per hectare (Singh et al., Reference Singh, Sanderson, Field, Fidelis and Yinil2019). Compared to the control plot, the mineral fertilizer treatment led to higher yields in 2016, 2017, and 2018. Despite this, there was no significant measurable difference in the soil properties between the control and the mineral fertilizer plot.
Calcium and Magnesium availability
Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019) reported that the poor initial yields in the fertilizer treatment were due to the low availability of soil Ca and Mg. By contrast, in 2018, the soil analyses detected low Ca and Mg contents (both total and exchangeable) for the compost (alone) and the mineral fertilizer + compost treatments (C and E). The Treatments C and E were among the most productive plots, while their soil Ca, and Mg contents were similar to controls. These new results suggest that point measurements of soil Ca and Mg contents may not be good indicators of crop productivity because low-nutrient contents may also be a result of increased plant uptake. Presumably, by promoting the activity of beneficial microorganisms (Bünemann et al., Reference Bünemann, Schwenke and Van Zwieten2006), compost could enhance the availability of nutrients like Ca. This could have resulted in similar soil Ca and Mg contents between the composted and non-composted plots. The difference would simply be that plants absorbed more nutrients under the influence of compost, therefore resulting in similar levels to other treatments. Determining the tissue contents of these elements would be necessary to confirm this hypothesis.
Exchangeable Aluminum
Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019) also associated low initial yields with high concentrations of exchangeable Al and low pH. The exchangeable Al contents at year 7 ranged from 21 to 30 ppm for Treatments D, F, and G, which received dolomite, up to 705–743 ppm for Treatments A, B, and C, which did not. In both cases, the soil Al concentrations are substantially below the levels reported by Shamshuddin et al. (Reference Shamshuddin, Muhrizal, Fauziah and Husni2004). Our results showed that high productivity was obtained for the compost (Treatment C) and the compost + mineral fertilizer (Treatment E) plots, which featured high exchangeable Al contents (720 and 678 ppm; also reported for Treatment C by Mulia et al., Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019), comparable to the control and mineral fertilizer-only treatments (A and B; 705 and 743 ppm). On the other hand, we found low exchangeable Al concentrations for two low-yield treatments, dolomite-only and dolomite + mineral fertilizer (D and F; Table 3). These observations imply that exchangeable Al, even though it may affect yields, was not a determining factor. The presence of organic matter may still offset Al toxicity, but the soil analyses revealed no difference in C contents across all treatments. Nevertheless, with the same C contents, differences in the microbiological profiles and the type of organic compounds present could play a role in mitigating Al toxicity.
Cocoa response to treatments
Seven years after field planting, the highest cumulative yields were still achieved from those treatments that had received compost (e.g., Treatments C, E, G, and H). The most effective single treatment to improve cocoa productivity was compost (C). The treatment including compost alone (C) led to higher survival rates, basal areas, and dry bean yields than the control (A), mineral fertilizer (B), and dolomite (D) alone treatments (except in 2018 for dry bean yields). Of the combined treatments, the treatment without compost (F) resulted in lower yields than those with compost (E, G, and H).
From year 4 to 7, the yield benefits of compost tended to decline relative to the other treatments. In 2015, which was a particularly dry year (Table S4), the benefit from compost was particularly strong. However, by 2018, the cocoa yields had become similar within all the treatments receiving additions, although they remained greater than the control treatment yields. Hence, as the cocoa bushes became more established, the benefits of adding organic material declined, and the benefits of supplying specific nutrients such as nitrogen (removed during harvesting) became more important.
The results indicate that the basal area of 7-year-old cocoa plants was not a good predictor of yield, since the largest cocoa trees were not the most productive (Figure 1). Verlière (Reference Verlière1981) suggested that the positive relationship between basal area and yield is only significant for younger trees (Jones & Maliphant, Reference Jones and Maliphant1958; Longworth & Freeman, Reference Longworth and Freeman1963), and this was observed in the initial phase of this experiment (Mulia et al., Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019).
Despite the various soil measurements, no immediately obvious relationship was found between soil nutrient levels and yields. For example, although yields varied between treatments in 2018, soil nitrogen levels across the treatments were similar. By contrast, there was evidence of the yield benefits of organic matter and an increase in pH. Nevertheless, it appeared that yield benefits from compost and/or dolomite addition had declined by years 6 and 7. The benefits of compost may have been particularly pronounced in 2015 because of the low rainfall in that year (Table S4). Another factor related to compost application could be the quantity of organic P provided by the compost (50% more than the fertilizer input). For example, a strong positive correlation (r = 0.85) between organic P and cocoa yields was reported for Southern Nigeria (Omotoso, Reference Omotoso1971).
One difficulty limiting the results’ interpretation is that soil characteristics were only measured at the beginning of year 4 and the end of year 7, whereas the yields were measured continuously. In addition, the lack of significant differences between treatments in many chemical soil properties also constrained the identification of relationships. One possible way of addressing this is to construct models describing the inputs and outputs of nutrients.
Over the first 7 years, one of the lowest cumulative yields was achieved in the plots where only mineral fertilizer was applied. These results do not mean that mineral fertilizers should be avoided. Instead, one of the largest cumulated yields was obtained from the compost + fertilizer (E; three times the control), and fertilizer alone led to double the yield of that from the control treatment. Furthermore, the rate at which productivity was increasing during the last 3 years for the fertilizer-only treatment (B) suggests that productivity could continue to rise substantially as the cocoa matures. Conversely, the average yield of the four composted treatments all declined during the last 3 years, while the non-composted increased. Another argument in favor of fertilizer is the high-yield index reached for this treatment, surpassing the others.
As described before, it appears that the yield response was different between the composted and non-composted treatments. In 2015 (dry year), the response to treatments was strong for the composted treatment and less for the non-composted ones. In the following years, the yields slowly declined for the composted treatment, while the non-composted gradually increased. In 2018, productivity was equivalent for all, except for the control (lowest average yield) and the fertilizer treatment (highest average yield). We could hypothesize that in the long run, fertilizer can be particularly helpful in maintaining productivity, while compost is useful in encouraging high yields during establishment and particularly during drought. Continuing the analysis during more years as the cocoa further matures could support or contradict these trends.
It is also important to note that we did not observe additional effects on yields when combining the treatments. For example, the yields of the treatment combinations with compost were approximately at the same level as compost alone. In Treatment E, there was no benefit from using the fertilizer as compared to compost alone. The fertilizer + dolomite treatment was slightly lower than fertilizer alone, suggesting that the addition of dolomite reduced the benefits of fertilizer application. The full combination did not give the highest yields. However, we should add that, if we consider the productivity including mortality rates (averaging per 16 trees), dolomite + fertilizer showed an improvement. The cumulated productivity of Treatment F (dolomite + fertilizer) was also 1.27 times Treatment B (fertilizer alone). In addition, average survival rates were also increased by adding dolomite. The survival rate of Treatment F (88%) was more than double that of Treatment B (41%) and showed an improvement as compared to the control. Treatment A had a survival rate of 67%, which means that adding dolomite approximately halved the mortality rate of the control treatment (33% for A and 12% for F).
A secondary cause suggested by Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019) for the response to the treatments was the young age of the cocoa trees. We can presume that young cocoa trees are more sensitive to environmental factors than older trees. In our study, cocoa age seems to be an essential factor since yields increased over time but also became almost homogenous the last year. The explanation behind the harmonization of yields by 2018 across the trial is uncertain. The positive effects of amendments may have been only useful in the first years of cocoa growth and development, but yield disparities between treatments had narrowed after 7 years. One reason could be as simple as the design of the trial itself. A block size of 16 trees separated spaced by 3 m may not be enough. Belowground, root systems may now be accessing adjacent plots and therefore blending the responses. Expansion of the root zone, possibly resulting in access to soil nutrition sources in adjacent plots, may have been influential.
Developing fertilization and amendment recommendations
The mean annual rate of change of soil nutrients can provide insights into the quantity of each element that is either taken up or lost by leaching or volatilization (Table S12) and into appropriate cocoa fertilization rates. There were major differences in the contents of certain soil nutrients between those reported for 2014 by Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019) and those reported here for 2018 (Figure 2 and Figure S1). For example, available P and exchangeable K declined drastically in 4 years. This could result from substantial leaching associated with the low-nutrient retention capacity of the soil. However, total N concentrations were maintained over 4 years. Another explanatory factor for the decline in nutrients is the uptake by cocoa trees and the associated shade trees at a rate higher than the rate of inputs. Based on those averages, the stock of extractable K, Ca, and Mg, declined, even in the control plot, suggesting that the supply of these nutrients may be insufficient for this plot. Considering organic C, the statistically significant changes that occurred over 4.5 years for Treatments C, G, and H correspond respectively to annual rates of change of 1671, 1986, and 1853 kg of C ha−1 yr−1, while 930, 1580, and 1580 kg ha−1 of C were applied yearly. About 896 kg of C ha−1 yr−1 were lost on average in the control plot, almost equivalent to the compost treatment inputs. However, for the composted treatments, the loss of C was 1.5 to 2 times that in the control treatment. This could suggest that the application of compost with a low C:N ratio was associated with increased decomposition of the preexisting soil organic matter (Kuzyakov, Reference Kuzyakov2010), whereas the incorporation of more recalcitrant forms of organic inputs could have resulted in greater organic matter stability. The reported rate of change of extractable P was extremely large, and initially assuming a technical error, the Asian-Agri laboratory manager double-checked the 2018 results, and the analysis of total P at Cranfield University was repeated three times. In theory, ‘extractable’ (Asian-Agri) (Table 2) and ‘total’ (Cranfield) contents (Table 3) should be similar (i.e., ‘pseudo-total’), but on this occasion, the quantities varied widely. Unfortunately, the 2014 and 2018 soil samples were not stored by Asian-Agri, preventing additional analyses. Because Cranfield’s results were consistent, it is possible that there is a technical error in the values of extractable P reported from Asian-Agri in 2018. This issue raises the question of the consistency and comparability of chemical analyses between different laboratories, methods, and years. The soil analyses presented by Mulia et al. (Reference Mulia, McMahon, Purwantara, Bin Purung, Djufry, Lambert, Keane and Guest2019) were produced by ICCRI (Indonesian Coffee and Cocoa Research Institute) and AIAT (Assessment Institute of Agricultural Technology), while ours came from Asian-Agri and Cranfield University laboratories.
The above analysis demonstrates the difficulty in developing appropriate fertilizer and amendment recommendations based on soil nutrient measurements alone. However, soil measurements can provide insights within an experiment examining yield responses. Two other tools for developing fertilizer recommendations are analysis of plant tissue and nutrient budgets. A nutrient budget requires considering the flows in terms of inputs and the quantity of nutrients being removed either as harvested pods or stored in the cocoa and shade trees’ accumulated biomass. Nitrogen can illustrate this issue since additional inputs may have increased yields but did not lead to a measurable change in 2018.
Conclusions
The effects of the treatments on soil properties were variable, with distinct changes for some variables (such as organic C, pH, and extractable nutrients) and little to no response for others (such as total N). Composted treatments resulted in the highest cocoa cumulative yields (on average, 2.8 times the control). In contrast, the addition of mineral fertilizer or dolomite without compost provided yields that were 1.6 times the control. The composted treatments yielded significantly more pods than the other treatments. The relative benefits from compost (with a low C:N ratio) compared to fertilizer applications were greatest in the initial years of establishment, gradually declining as the cocoa matured. Furthermore, soil C contents were similar between treatments despite the inputs and were not adequate to raise soil C levels. This issue raises questions about the feasibility of improving soil carbon storage in cocoa systems. The effect of altering the C:N ratio of the compost could be an area for further study, as well as experimenting with other combinations of organic inputs. The results also demonstrate that developing a site-specific soil fertility management strategy cannot be based on soil nutrient analysis alone, but soil nutrient analysis can be useful when integrated with experimental yield results and the analysis of nutrient flows. While it seemed that adding compost was sufficient to support cocoa productivity, the applications may be unattainable to many farmers. For this reason, future research should evaluate other combinations of compost + fertilizer + dolomite, with lower compost application rates than our experiment, combined with other organic inputs, to determine which ones are the most cost-effective to meet cocoa’s needs.
Acknowledgments
This research project has been made possible thanks to the support of the British Council Newton Institutional Link Program with Indonesia (Grant number: 217494906), Cocoa Research United Kingdom (CRUK), Cranfield University Industrial Partnership, and Mars Wrigley. We would like to commend the Cocoa Research Centre staff in Tarengge (Sulawesi, Indonesia) for their precious help, including the provision of a team of technicians and a translator, local insights, in-country laboratory analysis as well as their logistical support. Data underlying this study can be accessed through the Cranfied University repository at https://doi.org/10.17862/cranfield.rd.14885847.
Supplementary materials
For supplementary material for this article, please visit https://doi.org/10.1017/S0014479721000132
Open access
