Tropical agroforestry has an enormous potential to sequester carbon while simultaneously producing agricultural yields and tree products. The amount of soil organic carbon (SOC) sequestered is influenced by the type of the agroforestry system established, the soil and climatic conditions, and management. In this regional-scale study, we utilized a chronosequence approach to investigate how SOC stocks changed when the original forests are converted to agriculture, and then subsequently to four different agroforestry systems (AFSs): home garden, coffee, coconut and mango. In total we established 224 plots in 56 plot clusters across 4 climate zones in southern India. Each plot cluster consisted of four plots: a natural forest reference, an agriculture reference and two of the same AFS types of two ages (30–60 years and > 60 years). The conversion of forest to agriculture resulted in a large loss the original SOC stock (50–61 %) in the top meter of soil depending on the climate zone. The establishment of home garden and coffee AFSs on agriculture land caused SOC stocks to rebound to near forest levels, while in mango and coconut AFSs the SOC stock increased only slightly above the agriculture SOC stock. The most important variable regulating SOC stocks and its changes was tree basal area, possibly indicative of organic matter inputs. Furthermore, climatic variables such as temperature and precipitation, and soil variables such as clay fraction and soil pH were likewise all important regulators of SOC and SOC stock changes. Lastly, we found a strong correlation between tree species diversity in home garden and coffee AFSs and SOC stocks, highlighting possibilities to increase carbon stocks by proper tree species assemblies.
Land-use changes in the tropics are responsible for approximately 10 % of the human-induced greenhouse gas emissions and are expected to remain the second largest source of carbon (C) emissions in the near future (Achard et al., 2014). Considering that tropical forest soils store a similar amount of organic carbon (692 Gt in the top 3 m; Jobbágy and Jackson, 2000) as the atmosphere (589 Gt C; Ciais, et al., 2013), and that tropical climates foster rapid organic matter decomposition, land-use changes can result in strong carbon fluxes into or out of the soil. The conversion of tropical forests to agriculture causes a release of stored soil organic carbon, often in the form of carbon dioxide, but it also results in a decline in soil productivity. To reduce carbon emissions from agriculture while simultaneously maintaining agricultural productivity it is necessary to identify and implement simple and cost-effective measures to store and capture carbon. In this context agroforestry practices, which integrate trees into agricultural systems, offer a unique opportunity to sequester atmospheric carbon while also growing food, diversifying incomes (e.g., from sale of wood, fruit and staple foods), and simultaneously providing numerous environmental benefits. These include mitigating soil erosion (Montagnini and Nair, 2004), improving soil structure (Lal, 2007), pumping up nutrients from the subsoil (Das and Chaturvedi, 2008) and sequestering atmospheric carbon (Lal, 2007; Nair et al., 2009). Agroforestry systems (AFSs) have higher soil organic carbon (SOC) sequestration rates than conventional agricultural systems (Nair et al., 2009) as the trees have comparatively higher litter inputs and are capable of inserting carbon deep in the soil with their root systems (Montagnini and Nair, 2004).
Furthermore, tree species diversity in AFSs can have a large impact on organic matter turnover as diverse species mixtures can add different qualities of organic matter which correspondingly influence soil microbial communities and decomposition processes (Six et al., 2002; Acker et al., 2002). Niche differentiation and resource partitioning may lead to a better use of space and nutrient uptake and thus increase ecosystem carbon inputs (Thakur, et al., 2015) . Although it is recognized that AFSs have many benefits, their C-sequestration potential, especially below ground, remains largely unexplored (Montagnini and Nair, 2004). Concentration and SOC turnover rates in AFSs vary significantly with biophysical site properties such as climate (Liu et al., 2011), vegetation, land-use types (Cadotte, 2013; Saha et al., 2010), soil type and texture (Six et al., 2002; Chaplot et al., 2010), land management (Hevia et al., 2003) and their interactions (Powers and Schlesinger, 2002). It is generally recognized that temperature and precipitation are the most important variables regulating SOC (Chaplot et al., 2010; Liu et al., 2011), since both affect the type of vegetation cover, the quantity of biomass production and the rate of soil organic matter turnover (Hevia et al., 2003).
Environmental variables and structural characteristics of
agroforestry systems in the study area (mean
It is estimated that there are approximately 25.3 million hectares (Dhyani et al., 2013) of AFSs established across India, whereby the type of AFS established depends on the biophysical site conditions and the socioeconomic status of the owners. Despite this, regional-scale studies evaluating the impacts the establishment these land-use types have on SOC stocks remain relatively scarce (Albrecht and Kandji, 2003; Mutuo et al., 2005; Saha et al., 2010). Here in this study, we quantified SOC changes associated with the conversion of forest to agriculture and subsequently from agriculture to four different AFS types (home garden, coffee, coconut and mango). Our plots were established across southern India in a broad range of biophysical conditions, ranging from semi-arid to humid climates and in soils with low- and high-activity clays. The objectives of this study were as follows: to quantify SOC stocks and changes to SOC stocks along a forest–agriculture–AFS trajectory, and to determine the biophysical drivers regulating SOC stocks and its changes.
The study was conducted across three states in southern India (Kerala,
Karnataka and Andhra Pradesh; Fig. 1), in an area that has a long-standing
history of diverse agroforestry-based land-use practices. The region has a
tropical monsoon climate, with a rainy season from May to October and a
distinct dry season from November to April. The mean annual precipitation
(MAP) ranges from 627 to 3422 mm, and mean annual temperature (MAT) ranges
between 21.9
Geographical location of the plot clusters in southern India. Each point represents a cluster of four plots: natural forest, agriculture and the AFSs at two ages (medium and old).
In this study we investigated how SOC stocks changed when forests are
converted to agriculture and subsequently to agroforestry systems in four
different climatic zones along a precipitation gradient (humid, moist
sub-humid, dry sub-humid and semi-arid; based on classification by ICAR, 1984).
Each AFS type was sampled in two of four climatic zones based on
its relative importance in terms of area, production and income in the
region (Table 1). Using a chronosequence approach we established 56 plot
clusters (Fig. 1). A plot cluster consisted of four land-use types: a
natural forest reference and an agriculture reference and two AFS plots at
different ages: medium (30–60 years) and old (> 60 years). Seven
plot clusters were established for each AFS type in two climatic zones. In
total 224 plots were set up (4 AFSs
The plot clusters were always centered around the natural forest reference
plot and located within a maximum distance of 1 km of each other. To reduce
edge effects, the forest plots were selected at least 40 m from the forest
edge. Plot pairs (either forest with agriculture or agriculture with AFS)
were carefully selected to ensure soil and climate conditions were similar.
(1) All plots were located on similar landscape positions on flat to gently
sloping terrain (average: 3 %; maximum: 7 %). (2) Using a feel test we
compared the subsoil texture of all possible candidate sites and chose only
those with comparable soil textures. An a posteriori texture analysis revealed that
there were small differences in surface clay percent at 10–30 cm. For the
forest to agriculture comparison, the clay content difference was
Each AFS type consisted of a unique combination of trees and crops (Table 1) and originated from former agricultural land with the exception of 19 coffee plots and 1 home garden plot which replaced forests directly. Pictures of the investigated land-use types in the different climate zones are found in Fig. S1 in the Supplement.
The home garden AFS has a multilayered canopy, consisting of a diverse tree admixture of different ages and sizes. The multipurpose trees grown here are found in association with shrubs and herbaceous species (Kumar et al., 1994). The majority of households in the humid and moist sub-humid climate zones manage home garden AFS to (partially) satisfy their fruit, spice and vegetable needs (such as cassava, banana and ginger). For management, farmers loosen the soil once a year during vegetable cultivation (Fig. S1).
Coffee AFSs are mainly grown in the humid and moist sub-humid region of the Western Ghats. Coffee is grown in the understory of both native shade tree species and planted trees and is often inter-cropped with spices such as pepper, cardamom, cinnamon, clove and nutmeg (Fig. S1). Litterfall from shade trees and pruning products mean that this system receives substantial organic matter inputs. The soils in coffee AFS are typically hand tilled once every 2–3 years.
Coconut is primarily grown by smallholder farmers in southern India. The four southern states Karnataka, Andhra Pradesh, Tamil Nadu and Kerala together produce 92 % of India's total coconut production. In the humid region of Kerala state, coconut is extensively grown with vegetables, whereas in the dry sub-humid region it is grown together with food grains (such as maize, turmeric and finger millet) (Fig. S1). These crops are cultivated at different stages of coconut plantation development depending on the amount of incoming light. Since most of the coconut plant parts (leaves, stems) are useful to the farmers for various household uses, little organic matter is left on site. In the dry sub-humid zone, farmers plough the land two to three times a year for agriculture crop production.
Mango is an important commercial fruit crop in India which ranks first among world's mango producing countries, accounting for about 40 % of the world production (Sekhar et al., 2013). In Karnataka and Andhra Pradesh states, mango is predominantly grown in drier climate (dry sub-humid and semi-arid) and is predominately cultivated by smallholder farmers. Many farmers adopt wide row spacing and grow agriculture crops between the rows (Fig. S1). Since it normally takes 12–15 years to establish a closed canopy, farmers utilize this duration to cultivate agriculture crops, such as finger millet and maize. In later stages, those crops are no longer profitable and farmers switch to grow fodder and short rotation crops for the household consumption.
Agriculture under humid and moist sub-humid climate zones is only practiced by smallholder farmers. Here, staple foods like tubers and vegetables for subsistence use are grown, with only organic matter inputs as nutrient supplement sources (Fig. S1). In the dry sub-humid and semi-arid climate zones commercial agriculture crops are grown with a high input of both fertilizers and pesticides. All agricultural fields are typically ploughed two to four times a year, often with a small tractor. During agriculture establishment the forests were cleared using hand tools, and most aboveground biomass was removed for domestic use. Sites were not burnt.
The forest plots of our study were most often community managed and considered “sacred groves”. For religious reasons people do not remove any wood from there. The remaining forests which we sampled were government owned; these too were relatively undisturbed. All forests were located in or within the vicinity of the village and were highly protected. Evergreen forests were found in humid climate, moist deciduous forest in moist sub-humid climate zone, dry deciduous forest in dry sub-humid and scrub forest in semi-arid climate (Fig. S1).
In each 20
Total soil carbon and nitrogen concentrations were analyzed using a CN
Elemental Analyzer (Vario EL III, Elementar, Hanau, Germany). As the entire
region is underlain by high-grade metamorphic rocks and granites of the
Indian Shield, no carbonates were expected in these soils and no attempts
were made to remove them. The soil carbon stock (Mg C ha
At each plot, we measured tree basal area for all trees with a diameter at breast height greater than 10 cm. These tree species were identified to the species level. Furthermore, we recorded information on slope, elevation and geographical coordinates of each plot. Through informal interviews with the land owners we got information on current and past land uses and their management practices. Meteorological data such as mean annual temperature and mean annual precipitation for the selected plots were retrieved from the WorldClim database (Hijmans et al., 2005).
To verify that plots were satisfactorily selected in the field and that the
soils of the plot pairs were inherently similar we did an a posteriori comparison of soil
clay percentages in the subsoil (10–30 cm) of the plot pairs using a
paired
The residuals of all models were checked for normality with Q–Q plots; models
were considered significant at the
Among the land-use types investigated, SOC concentrations (data not shown) and SOC stocks were highest in natural forest, lowest in agriculture and intermediate for the different AFSs (Fig. 2). Deforestation for agriculture resulted in a strong decrease in SOC stocks across all climate zones, creating a loss of 50 to 61 % of the original SOC stock. The resulting SOC stocks of AFSs however were dependent on AFS type. While SOC stocks in home garden and coffee AFSs rebounded to near forest levels, the SOC stocks in both coconut and mango AFSs increased only marginally compared to the agricultural reference (Figs. 2 and 3).
Absolute SOC stocks (0–100 cm) in the four different land-use
systems. Each graph shows the respective forest and agriculture references
and the two AFS age categories (medium and old) in the two different
climatic zones. Point represents the mean of seven plots, while error bars
indicate the 95 % confidence intervals based on the Student's
Of the AFSs studied, SOC stocks in the top meter of soil were highest in
coffee (156
Relative change in SOC stock in the 1 m soil profile from
agriculture to agroforestry systems: home garden (
The climate zone had a marginal influence on overall SOC stocks of a given AFS where only coconut and mango AFSs showed a tendency towards higher SOC stocks in the respective wetter climate zone (Fig. 2). Also the SOC stocks of the AFSs did not differ among the two age categories sampled (Fig. 2).
In undisturbed forest ecosystems 82 % of the variance in SOC stocks could
be significantly explained by five variables (MAP, MAT, basal area, clay
fraction and soil pH) using stepwise multivariate regression analyses
(Table 2). While tree basal area and clay fraction exhibited positive linear
correlations with SOC stock, both soil pH and MAP exhibited parabolic and
inverse parabolic relationships respectively (Fig. 4a–d). For pH, SOC stocks
were lowest at near-neutral conditions, while SOC stocks peaked between
2000 and 3000 mm yr
Multivariate regression models predicting SOC stocks in different
land-use systems and the relative changes in SOC stock (%) from either
forest or agriculture references in the 100 cm soil profile. MAP is mean
annual precipitation (mm), MAT is mean annual temperature (
Scatterplots
Basal area was the single best predictor of forest SOC stocks. This is evident in both the scatterplot graphs in Fig. 4b and from its dominant position in the multivariate regressions for all land uses investigated (Table 2), except in agriculture which had no trees. Likewise, soil clay fraction was also an important predictor of SOC and was present in all AFS regression equations except home gardens. Although MAT was an influential predictor of forest SOC stocks in the stepwise regression (found in three of the six prediction models; Table 2), yet taken as single predictor its influence on SOC stocks remains insignificant. SOC stocks in both home garden and coffee AFSs are further positively correlated to the Shannon–Wiener species diversity index of trees (Fig. 5). Coconut and mango AFSs however were not included as they are monocultures in terms of tree or palm species admixture.
Scatterplot showing the relationship SOC stock of home garden and
coffee AFSs exhibit with plot Shannon–Wiener species diversity index:
home garden in the humid zone (
The SOC stock losses attributed to the conversion of forests to agriculture could be predicted by two variables: clay fraction and MAT (Table 2). Thereafter, when agriculture plots were converted to AFSs, SOC stock changes could be predicted by MAP, MAT, basal area, clay fraction and soil pH in varying importance for the individual AFSs and for all AFSs combined (Table 2). While MAP, basal area, and clay fraction all exhibited positive linear correlations with SOC stock change, soil pH exhibited a parabolic relationship with SOC change. For the latter, SOC stock losses were highest in acidic soils, lowest near neutral pH, and again higher in slightly alkaline soils (Fig. 4e–h).
The SOC stock is a balance of incoming carbon from organic matter and carbon losses either through decomposition processes or dissolved organic carbon (DOC) leaching (Davidson and Janssens, 2006; Raich et al., 2006). Both the carbon inputs and outputs however are strongly affected by ecosystem productivity, vegetation type, climate, clay mineralogy, soil pH, nutrient availability, soil aggregates and texture (Lal, 2004; Six et al., 2002; Chaplot et al., 2010; Don et al., 2011).
Although we did not measure organic matter inputs directly, we found a very strong correlation between SOC stocks and plot basal area, which could be indicative of organic matter inputs (Fig. 4b; Lebret et al., 2001). Once organic matter enters the soil, the soil's physical characteristics and biochemical environment play an important role in how SOC is either stabilized, mineralized or leached (Raich and Schlesinger, 1992; Six et al., 2002; Zinn et al., 2007a). The clay mineralogy and soil clay fractions play a critical role in how much SOC can be stored by the soil, its residence times and its susceptibility to land-use change (Zinn et al., 2007b). Clays can stabilize organic matter particles through clay–humus complexation, or physically protect organic matter from further mineralization by trapping them in clay macroaggregates (Sollins et al., 1996). The strong positive correlation between SOC stocks and clay fraction we measured reflect the importance of these processes on SOC storage (Table 2, Fig. 4c). Next, the parabolic relationship we observed between soil pH and forest SOC stocks (Fig. 4d) reflect how the soil biochemical environment is critically important for soil microbial communities that decompose organic matter (Motavalli et al., 1995). The higher SOC stocks found in the acidic and basic soils in this study indicate that the unfavorable biochemical environment retards microbial communities that decompose the organic matter. At a near-neutral soil pH, conditions were ideal for microbial communities and accordingly decomposition rates were high, resulting in little SOC accumulation.
As expected, regions of higher rainfall stored more SOC than drier zones
(Fig. 2). This is likely because the natural forests in the humid zone have
higher net primary production given the favorable year-round water
availability compared to forests in the drier regions. Unexpectedly, however,
at very high levels of precipitation (> 3000 mm yr
Although MAT is a significant predictor of forest SOC stocks, its influence can only be evaluated in interaction with other variables (Fig. S3a). Although one would expect a strong correlation between MAT and MAP, it was not present here. While MAT is mainly driven by altitude, MAP depends on more complex weather phenomena that differ among climate zones. This might also explain why there is no clear relationship between MAT and SOC when other variables are not included in the analysis.
In comparison to a large pool of studies on this land-use conversion conducted in the tropics, the SOC losses detected here for the different climate zones were at the higher end (50–61%). Results from two meta-analyses report SOC stock decreases on average between 18 % (Powers et al., 2011) and 25 % (Don et al., 2011) for this land-use conversion in the tropics. Powers et al. (2011) however report that large SOC losses are possible (to a maximum of 76 %), but the losses depend on the clay mineralogy and precipitation regimes. Nevertheless, in our study we measured significant decreases in SOC stocks irrespective of soil type and precipitation regimes. We primarily attribute these large SOC losses to the frequent tilling and low organic matter inputs. Tillage exposes SOC to microbial activity through the destruction of aggregates which as a result makes SOC complexes vulnerable to decomposition (Six et al., 2002; Mangalassery et al., 2013). Furthermore, fine-grained particles (such as clay) and associated organic matter can also be lost through soil erosion, runoff and leaching (Gonzalez and Laird, 2003) due to a lack of soil protection measures especially following plowing at the onset of the monsoon (Dourte et al., 2012). Associated with the erosional clay losses is a corresponding reduction of the soil's carbon storage potential. However, since all our plots were established on flat to gentle slopes, and because we were only on site for 1 day, we did not measure soil erosion.
As previously reported by van Straaten et al. (2015), we also found that SOC stock losses were proportional to the initial forest SOC stock, whereby the higher the SOC stock was initially, the larger the corresponding SOC stock loss when converted to agriculture (Fig. S4).
Increases in SOC stocks resulting from agroforestry establishment are ultimately attributed to higher organic matter inputs from above- and belowground sources (leaves, wood, roots, fungi, animals, etc. (Montagnini and Nair, 2004) and a reduction of SOC losses from decomposition and leaching. Similar to the carbon stocks in forests, MAP, basal area, clay fraction and pH also control carbon stock changes when agricultural land is converted to agroforestry (Table 2, Fig. 4e–h). Furthermore, the accumulation of SOC depends highly on the quality of incoming litter (Lemma et al., 2006) and is reflected in the significantly higher soil C : N ratios in home garden, coconut and mango AFSs (Fig. S5). In comparison to litter from agricultural crops, which generally have low C : N ratios and decompose rapidly, organic matter inputs from trees are generally of poorer quality (higher C : N ratios) because of the higher lignin and polyphenolic contents, which in turn results in slower decomposition rates and more SOC accumulation (Davidson and Janssens, 2006).
The large SOC stock increases in the subsoil (Fig. 3) indicate that belowground carbon inputs from roots and/or leaching of organic acids and soluble humus fractions to deeper layers are important processes for SOC accumulation. Considering the trees of the four AFSs have deeper rooting profiles than agricultural crops and that often more than half of the carbon assimilated by trees is transported belowground for root production (Montagnini and Nair, 2004; Poeplau and Don, 2013), it is no surprise that SOC stocks also increased substantially at depth. Furthermore, tillage activities will have mixed soils in the top 30 cm and therein homogenized soil carbon concentrations to a certain extent (Yang and Kay, 2001). However, the size of the SOC stock change hinges on the type of AFS established and its management practices. While all AFSs gained carbon compared to agriculture, the amounts gained varied strongly between the different types (Fig. 2). While coconut and mango SOC stocks increased just marginally above the agriculture reference, home garden and coffee SOC stocks rebounded to forest levels. Clearly, the carbon cycling dynamics of both home garden and coffee AFSs resemble that of natural forests since both AFSs support many different tree species of different ages, have varied stand structures and have high basal areas.
However, not only carbon input, but also losses are a function of AFS type, especially in terms of the plantation management schemes. In coconut and mango AFSs (especially in the dry regions), the low SOC stock increases are linked to the removal of crop residues (including leaves) from the site, which are used as fodder or fuel. In coconut, even the palm fronds are removed and utilized.
The effect of climate zone on SOC stock changes can only be quantified within each AFS type. While the climate zone did not affect SOC stock changes in either coffee or home garden AFSs, only a slight difference was found for coconut and mango (Fig. 2). This is however more likely attributed to the implementation of different management practices and cannot be disentangled from climate itself. For instance, in humid climates coconut farmers utilize a “planted fallow system” which includes a fallow period in the cropping cycle where organic matter is reintegrated into the soil. Such systems have been shown to improve soil fertility and maintain SOC stocks (Salako et al., 1999). In contrast, in the dry sub-humid climate farmers use a continuous cropping system which inevitably results in lower organic matter inputs and therein lower SOC stocks. The age of each respective AFS also had no significant effect on SOC stocks (Fig. 2), indicating that soil carbon had already reached a new equilibrium within the first 30–60 years. This is consistent with literature which reports that a new SOC equilibrium can be attained in 20 to 40 years following land-use conversion (Detwiler, 1986; de Blécourt et al., 2013; Chiti et al., 2014).
Lastly, tree species diversity (Shannon–Wiener index) correlated strongly positively with SOC stocks in home garden and coffee AFSs (Fig. 5), highlighting the role that species and resource complementarity have in maximizing biomass production (Cadotte, 2013). Interestingly, both home garden and coffee AFSs investigated in this study had similar tree diversities as natural forests. Since plant diversity was also positively correlated with basal area (data not shown) it is possible that the favorable climatic conditions where these AFSs exist can allow both high species diversity and high ecosystem productivity. Literature has shown that plant diversity is integrally linked to ecosystem productivity (Cadotte, 2013) and ecosystem resource utilization (Tilman et al., 2012), which both affect SOC storage potential (Thakur et al., 2015). In contrast, monocultures have been shown to have lower organic inputs than species-diverse systems (Cardinale et al., 2007)
Agroforestry systems provide a unique opportunity to produce food and tree products, while also improving livelihoods, protecting and improving soils and, as discussed here, to sequester carbon. Nevertheless, not all AFSs provide the same benefits. Soils in home garden and coffee plantations for instance can sequester much more carbon than coconut or mango AFSs. Additionally, the soil carbon sequestration potential of AFSs can be maximized by cultivating a broad range of different tree species, minimizing tillage activities and leaving crop residue on site.
We kindly acknowledge financial support from the International Fellowship Program of the Indian Council of Agricultural Research (ICAR), Government of India. We also thank N. B. Prakash at the University of Agricultural Sciences in Bangalore, for his support and guidance and our field assistants M. R. Jagadish and B. S. Sathya. Lastly, we also thank the two anonymous reviewers for their critical yet constructive reviews. This open-access publication was fundedby the University of Göttingen. Edited by: A. A. Berhe