- Open Access
Carbon stock under major land use/land cover types of Hades sub-watershed, eastern Ethiopia
Carbon Balance and Managementvolume 14, Article number: 7 (2019)
Developing land management scenarios that have the potential to sequester carbon and reduce greenhouse gasses (GHG) emission on a sustainable basis entails quantifying the current carbon stock under different land uses. In light of this, a study was conducted in Hades sub-watershed, eastern Ethiopia, to explore the carbon stock under four major land uses: natural forest, coffee agroforestry, grazing land and cropland, involving samples from four carbon pools: aboveground, belowground, litter, and soil. To this end, vegetation and soil samples were collected from the respective land uses following recommended procedures.
The results indicate that organic carbon concentration decreased with soil depth though substantial amount of carbon was found in the lower soil depths under land use with woody perennials. The mean total organic carbon stock ranged from 138.95 t ha−1 in the cropland to 496.26 t ha−1 in the natural forest. The soil organic carbon stock was found to be relatively higher than that of the vegetation carbon stock in the natural forest and coffee agroforestry land uses.
This study highlights the importance of assessing watershed level carbon stock for better and carbon friendly land use decision making. Land uses with woody perennials have high carbon stock than those without. Hence, conservation-based production systems with inclusion of woody perennials are options suggested to enhance carbon sequestration in the sub watershed.
The prominent factor deriving climate change is the increase in the concentration of greenhouse gases (GHGs) in the atmosphere. The ever increasing concentrations of carbon dioxide (CO2), methane (CH4), nitrous oxide (N2O), and other GHGs have distorted the balance between the incoming and the outgoing solar radiation emitted by the Sun and Earth, respectively . IPCC  estimated the increase in global temperature because of radiative forcing of GHGs in the atmosphere at 0.6 °C in the 1990s, while highlighting that this is expected to increase by 1.4 to 5.8 °C by 2100. Carbon dioxide, the second most important GHG, alone has increased from its pre-industrial 280 ppm to the current 397 ppm (International Energy Agency .
The two most common factors distorting the carbon cycle are land use change and combustion of fossil fuel . Practically, about 75% of the global CO2 emissions result from the combustion of fossil fuels in transportation, in building heating and cooling, and in the manufacture of cement and other goods . In Ethiopia, the emission from fossil fuel generated 2.3 million tons of CO2 in 1990, and the figure increased to 8.5 million tons of CO2 in 2013 .
Land use change, especially conversion of a natural system into a managed system, results in alteration of the carbon balance . There is also tremendous evidence showing the negative impact of agriculture on carbon stock [7,8,9,10,11]. However, agriculture is among the land use practices that emit as well as sequester CO2. It may lose soil organic matter due to intense decomposition following soil plowing, removal of above ground biomass during harvest, and severe soil erosion inherent in these activities . Nevertheless, if proper land use and management practices are put in place, agriculture can serve as an important sink [12,13,14,15].
In the tropics, deforestation is the second most important source of greenhouse gases (GHG) after fossil fuel combustion . Conversion of forest land to other land uses, such as agriculture, enhances decomposition and removal of carbon through harvest [8, 9, 17,18,19,20]. On the contrary, a significant increase (50%) of soil carbon was reported after the conversion of arable land into forest land [21, 22].
The positive effects of other land uses, such as agroforestry, on carbon balance were also highlighted. In soils that were previously under cultivation, agroforestry systems were found to be highly effective in restoring soil carbon . Similarly, higher aboveground carbon in coffee agroforestry (61.5 ± 25.0 t ha−1) than in woodland, pasture, and cropland, but slightly less than that in natural forest (82.0 ± 32.1 t ha−1) was reported in south western Ethiopia . The same study further indicated that about 59.5 t ha−1 organic carbon could have been lost if the coffee agroforestry had been converted into cropland.
Similar to the other land uses, rangelands and grazing areas are storing carbon above and below ground. In Europe, an increase in carbon stock following conversion of cropland into grass land was reported . Though there is scanty information on carbon stock of grazing lands in highland areas of Ethiopia, 128.39 t ha−1 belowground (soil and root) and 13.11 t ha−1 above ground, organic carbon was reported in the communally managed semi-arid rangelands in southern Ethiopia .
Hence, understanding the relationship between land-use systems and carbon stock is essential, as every land use system has either positive or negative impact on the carbon balance. Besides, considering the potential and constraints of a watershed in relation to carbon stock is vital to set recommendation on the maintenance and enhancement of carbon stock. In Ethiopia, few area-specific researches have been conducted on organic carbon stock and sequestration. The available researches are limited in their scope, mostly concentrating on soil carbon and giving much less emphasis to the carbon stock of the various carbon pools at watershed level. Furthermore, the impact, on carbon stock, of smallholder subsistence farming, mostly characterized by low level of management where limited or no inputs are used, under high population pressure has not been studied comprehensively. Under the looming climate change, such information on carbon emission and sequestration is essential for developing strategies that enhance productivity (through increasing carbon stock) and abate greenhouse gas emissions (through enhancing carbon sequestration). In view of this backdrop, this study assessed the current carbon stock under different land use types and carbon pools in Hades sub-watershed, eastern Ethiopia.
Description of study area
The study was conducted on major land uses of Hades sub-watershed in eastern Ethiopia. It is located at about 401 km away from Addis Ababa along the highway to Dire Dawa and Harar cities. The geographical location of the sub-watershed is 9°12′0″9°19′0′′ N and 41°15′0′′–41°17′0′′ E (Fig. 1). The altitude of the study area ranges from 1995 to 2775 m.a.s.l. The precipitation of the area is characterized by bi-modal distribution with mean annual rainfall of 930 mm. The mean annual temperature of the area ranges from 10.3 to 18.9 °C. The Hades forest used to cover most of the present crop and grass lands and coffee agroforestry sites. Currently, the forest is degraded due to continuous disturbance. The remaining dominant tree species are Afrocarpus gracilior, Croton macrostachyus, Hagenia abyssinica, Schefflera abyssinica, and Prunus africana. The farming practice in the area is mixed farming whereby the community cultivates annual and perennial crops and manages livestock. The dominant agricultural crops are sorghum (Sorghum bicolor), maize (Zea mays), wheat (Triticum aestivum), barley (Hordeum vulgare), and beans (Phaseolus vulgaris). Soil bunds are constructed on croplands and coffee agroforestry in order to abate soil erosion by water. According to small-scale studies carried out by , the dominant soil types are Leptosols.
Delineation of major land uses of the study area
Candidate land use/land cover (LULC) categories were identified using topographic map (1:50,000) obtained from Ethiopian Mapping Authority and Google Earth online imagery and Digital Elevation Model (DEM: 30 m × 30 m resolution). Following this, physical observation was made to confirm the basic information about the major land use types and topographic variations of the study area. The boundary of the sub-watershed was delineated using GPS recordings. The LULC map of the study period was produced in ArcGIS 10.5 software. Accordingly, four major land use types (cropland, grazing land, coffee agroforestry, and natural forest) were identified (Table 1).
Sampling sites for sample collection were the strata that were determined based on the land use/land cover types. Two transects were established along the slope at 500 m interval crossing cropland, grazing land, coffee agroforestry, and natural forest guided by compass. Besides, random sample points were established in areas not covered by the transect (Fig. 2). On crop and grass lands, sample plots of 10 m × 10 m were laid at 500 m interval. On coffee agroforestry, sample plots of 20 m × 20 m were used. A ‘nested’ sampling approach [28, 29] was used for collecting both vegetation and soil samples in the natural forest. To avoid boundary effect, the first plot was established 150 m inside the forest to the border line. The location of each sampling points were recorded using GPS.
In the sample plot, all living trees and shrubs with diameters at breast height (dbh) of ≥ 5 cm were inventorized. The dbh of the trees and shrubs were measured using diameter measuring tape . All the litter samples in a 0.5 m × 0.5 m quadrat from the four corners and the center of 40 m × 5 m plot were collected and a composite sample was made. The aboveground carbon stock of coffee shrubs was estimated from the diameter measured at 40 cm from the ground using caliper . As suggested by , the belowground biomass (root) of a plant was calculated from the aboveground biomass. The sites that were used for aboveground carbon stock measurement were also used for belowground carbon stock measurement in both the natural forest and the coffee agroforestry.
Soil samples were taken from the four corners and center of the 40 m × 5 m sub plot for the natural forest and the 10 m × 10 m plot of the crop land, grazing land, and coffee agroforestry. The samples were taken from three depths (0–20, 20–40, and 40–60 cm) using auger  taking into account the current standard depth of 30 cm recommended for C accounting purposes [34,35,36]. Hence, three composite samples per plot were collected for laboratory analysis of organic carbon, total nitrogen, soil texture, and pH. For bulk density determination, undisturbed soil samples were collected from the same depths using core sampler and from the center of the respective plots. All disturbed soil samples were air-dried, grinded, passed through a 2 mm mesh sieve for determination of the selected soil properties except organic carbon and total nitrogen in which case the samples were crashed further to pass through a 0.5 mm mesh sieve.
For litter carbon assessment, a sub sample of 100 g from a composite sample was oven dried at 105 °C until a constant weight was attained. Then, a sample of 4 g was ignited in Muffle furnace (Carbolite Aston-Lan. Hope, England) at 500 °C for 8 h . For all the samples, the determination was made in triplicate.
Soil texture was determined by the Bouyoucos hydrometer method as described in . The dry bulk density of the soils was determined using the core method as described in  in which case the core samples were dried in an oven set at a temperature of 105 °C to a constant weight. The bulk density was obtained by dividing the oven dry weight by the volume of the cores as is indicated in Equation . Soil pH in water was measured in soil–water (1:2.5) suspension using pH meter . The soil organic carbon content was determined following the Walkley–Black oxidation method . Total nitrogen was determined by the micro-Kjeldahl digestion, distillation, and titration method .
Biomass carbon estimation
Estimation of vegetation biomass and carbon
Species specific allometric equations were not available for all trees in the study area. Therefore, the dry weight biomass of each tree was calculated using the following formula :
where AGB is aboveground biomass in kg, D is diameter at breast height (dbh).
As described by , root biomass in ecosystems is often estimated from root to shoot ratios. The ratio ranges from 0.18 to 0.30, with tropical forests in the lower range. As suggested by , the belowground biomass (root) of a plant is close to 20 percent of its total aboveground biomass. Accordingly, the root biomass of trees was estimated using the following formula:
The corresponding carbon for above and below ground living biomass was obtained by multiplying the biomass with 0.5 .
Litter (undecomposed and unburned fallen leaves, twigs, and branches) was collected from the natural forest and coffee agroforestry, while undergrowth vegetation (all living seedling, sapling, and herbaceous species) was only collected from the natural forest. Litter and undergrowth vegetation samples were collected from 0.5 m × 0.5 m quadrants for natural forest and 10 m × 10 m plots for coffee agroforestry. The biomass and carbon stock in litter and undergrowth were calculated using the following equations :
Similarly, the organic matter content of the sample was derived from the following relationship:
where OM is organic matter, OC is organic carbon, and 1.724 is Van Bemmelen factor .
The following formula was used to calculate the biomass of the coffee shrub :
where Y is the biomass, d40 is the stump diameter at 40 cm height and b1 is the coefficient for the squared power equation for the above formula produced with the cross-validation ‘‘training’’ coffee plants in Ethiopia and whose value is 0.147 .
The carbon stock of the shade trees in coffee agroforestry was determined using the formula developed by :
Estimation of soil carbon stock
The organic carbon content of soils obtained from laboratory analysis was used to calculate carbon stock per unit area of land. The carbon stock for each layer per hectare was calculated as :
where Ci is C stock of the ith layer in t/ha, BDi is bulk density of the ith layer in kg/m3, CFi is coarse fragment content of the ith layer, OC is the soil content of organic C in %, di is thickness of the ith layer (m), 10 is conversion factor from kg/m2 to t/ha.
Coarse fraction was determined during sample preparation after repeated crushing of clods by hand, mechanical grinding, drying, and sieving until the sample was passed through a 2 mm sieve. Following this, the coarse fraction was weighed and its proportion was determined using the relationship :
The total carbon stock for the 0–60 cm depth was calculated as:
The dry bulk density was calculated using the following equation:
where MODS is the mass of the oven-dry soil (g); Vt is the total volume of the soil core calculated from:
where r is the internal radius of the cores measured using a caliber (cm), and h is height of the cores measured using a hand tape. 1000 is conversion factor from g/cm3 to kg/m3. π is a constant which is equal to 22/7.
The total carbon stock (ton/ha) of each land use of the sub-watershed was obtained from:
where CAG is aboveground carbon, CBG is belowground carbon, CLT is litter carbon, CUG is undergrowth vegetation carbon, and Csoil is soil carbon.
Calculation of deterioration index
Deterioration index (DI) of soils under coffee agroforestry, grazing land, and crop land was computed assuming that the level of organic carbon and total nitrogen contents under these land uses before conversion were once the same with soils under less influenced natural forest. The changes were averaged across the two dynamic soil properties to generate the overall soil deterioration index . The following relationship was employed to compute deterioration index for the individual soil properties:
where PSL is mean value of individual soil property (P) under specific land use (SL), PRL is mean value of individual soil property (P) under reference land use (RL), and DI is deterioration index.
The data were grouped and summarized by land use and soil depth. The standard error of mean was calculated for each parameter, and depths of the soil and vegetation. The data for individual depths were also subjected to ANOVA (analysis of variance). Least significant difference (LSD) was used to separate means that are significantly different from each other at p < 0.05. Besides, statistical differences were tested using two-way analysis of variance (ANOVA) to identify whether differences, in soil attributes, by land use and depth are significant following the general linear model (GLM) procedure of SPSS Version 20.0 for Windows.
Selected soil properties
The measured soil physical properties indicated significant (p < 0.05) differences in bulk density by land use and soil depth (Table 2). However, the interaction effect of land use by soil depth on bulk density was not significant (p ≥ 0.05) (Table 2). Across the three soil depths, significantly lower bulk density values were recorded in soils under the natural forest, whilst higher values were observed in soils under the crop land (Table 3). The other two land uses had intermediate values. Except in soils of the coffee agroforestry land use type, bulk density values increased down soil depth.
Sand, silt, and clay content, on the other hand, were significantly (p < 0.05) affected by land use alone (Table 2). Accordingly, significantly higher sand content was measured in soils of the natural forest across the three depths. On the contrary, significantly higher clay content was recorded at 0–20 cm soil depth of the cultivated and grazing lands (Table 3). Except in soils of the natural forest where it decreased, sand content did not show any consistent trend with soil depth in the other land use types. Silt content, on the other hand, exhibited a downward trend with soil depth in soils of the cultivated and grazing lands, showing no consistent variation with depth in the other land use types. Owing to some downward translocation, the clay content showed an upward trend with soil depth in all land uses but the grazing land, where it followed no consistent pattern.
The two way analysis of variance (Table 4) indicates that soil organic carbon and total nitrogen were significantly (p < 0.05) affected by soil depth and land use while pH was only significantly affected by land use (p < 0.05). These three soil properties, however, were not significantly affected by the interaction of land use and depth. On the contrary, C:N was neither affected by land use nor depth and the interaction of the two.
The pH of the soils under the four land use types exhibited significant (p < 0.05) variation across land uses but not within soil depth across the land uses (Table 5). At the surface layers, relatively higher pH value was recorded in the soils of the natural forest, whereas relatively lower value was recorded in the cultivated land. The soil organic carbon content was significantly different among land uses, and soil depths within a given land use (Table 5).
Across the three depths, significantly higher organic carbon content was found in the natural forest followed by the coffee agroforestry. The cultivated land, on the contrary, had comparatively lower values of soil organic carbon across its depths. In all the land uses, the soil organic carbon content significantly decreased with soil depth (Table 5). Similar to the soil organic carbon content, significantly high total nitrogen content was found in soils under natural forest followed by the coffee agroforestry (Table 5). This is also supported by the highly significant (p < 0.01) and positive correlation (r = 0.96) between total nitrogen and organic carbon content (Table 6). Furthermore, it also decreased with soil depth in all the land uses except the coffee agroforestry where it did not show any consistent trend down the soil depth. Similarly, sand, silt and clay correlated well with SOC (Table 6).
The C:N ratio, which is a derived parameter, did not follow the trend of soil organic carbon and total nitrogen with soil depth in all the land uses except the natural forest (Table 5). Comparatively high values of C:N ratio were recorded in soils of the cultivated and grazing lands as compared to the coffee agroforestry and natural forest soils. Across the land uses, it varied within a relatively narrow range of 11:1 to 17:1.
Soil and biomass carbon stock
Soil carbon stock
Comparing the four land uses, significantly higher soil organic carbon stock across the three soil depths was encountered in the natural forest. In coffee agroforestry, organic carbon stock was significantly higher compared with crop and grazing lands at 0–20 cm and 20–40 cm soil depths. However, there was no significant difference in SOC stock of the three land uses (NF; CL; GL) at 20–40 cm soil depth (Table 7). Similarly, no significant difference in SOC stock was observed between crop and grazing lands across soil depths. The decline in soil organic carbon was higher at 40–60 cm soil depth than it was at 20–40 cm in all the land uses except the natural forest where the decline was considerable at the 20–40 cm (Table 7). In all land uses, about 38–40% of SOC was found in the 0–20 cm soil layer while 68–60% was recorded in 20–40 cm and 40–60 cm soil depths respectively. In the natural forest and coffee agroforestry, respectively about 68% and 92% of the total carbon stock was found in the soil (Table 8).
The deterioration index clearly shows that the conversion of the once natural forest ecosystem into managed ecosystems in the form of grazing land, coffee agroforestry, and crop land has resulted in net degradation in organic carbon and total nitrogen levels. The decline in these soil attributes was more severe in soils that were under crop cultivation for the last many years (Fig. 3). This is followed by the grazing land. However, there was relatively lower deterioration index under the coffee agroforestry.
Biomass carbon stock
Due to absence of trees on sample plots under crop and grazing lands, biomass measurement was only made on natural forest and coffee agroforestry land uses. The natural forest was found to have significantly higher biomass carbon stock compared with the coffee agroforestry land use. In the natural forest, 81.5% share of the biomass carbon stock was attributed to the aboveground biomass. Particularly trees with dbh ≥ 30 cm had contributed the largest carbon. The contribution of the belowground (root), litter, and undergrowth vegetation in the natural forest was 18.5% of the total biomass carbon stock. In coffee agroforestry land use, only a small number of shade trees with small diameter were encountered and measured. Accordingly, the share of aboveground biomass (shade trees and coffee shrub) was 98.29% of the biomass carbon stock of coffee agroforestry land use. The undergrowth carbon stock recorded under the natural forest was 0.38 ± 0.08 t ha−1, which was lower than the other carbon pools (Table 8).
The presence of significantly more sand content in the natural forest than other land uses is probably due to the steep slope and high elevation where the forest is located. Nahusenay and Kibebew  also reported significantly higher sand content in soils under natural forest as compared to shrub, grazing, and cultivated lands in North-central highlands of Ethiopia. The presence of high clay content at 0–20 cm soil depth of the cultivated and grazing lands indicates the selective removal of the finer particles, such as silt and clay, by water erosion from the steep slopes where the natural forest land is located, and their subsequent accumulation in the gently sloping and low-lying parts of the study area where the grazing and cultivated lands are located. Ellerbrck and Gerke  pointed out that during erosion clay particles can be transported along hill slopes and form colluvic soil at the topographic depressions. In line with the findings of the current study,  reported accumulation of clay particles at lower soil depths.
The lower bulk density values in the natural forest could be associated with the high organic carbon content which encourages aggregation. Teshome et al. [52, 53] also reported significantly lower bulk density values in the natural forest as compared to grazing and cultivated lands in western and central highlands of Ethiopia. The inverse relationship between organic matter and bulk density was also reported by many researchers . The increase in soil bulk density with soil depth might be associated with the decline in organic matter, aggregation, and root penetration in addition to the expected compacting effect of the overlying soil mass. A number of studies conducted in Ethiopia have reported a general increase in bulk density value with soil depth [54, 55]. The bulk density values of the soils under the four land use types are not likely to restrict plant growth as a consequence of excessive compaction .
The lower pH value in soils of the cultivated land could be due to depletion of basic cations through crop harvest and continuous use of acid-forming fertilizers, such as di-ammonium phosphate (NH4)2HPO4), which produces strong acids when oxidized by soil microbes . The relatively higher pH values recorded at the bottom layers (40–60 cm) in some of the land uses could be related to the leaching of the basic cations from the upper layers and their subsequent deposition at the lower soil depths . Following soil pH rating suggested by , the pH values of the studied soils fall within the range of slightly acidic in the cultivated land to neutral in the other land use types. The pH values recorded in soils under the four land uses are within the range that is considered favorable for availability of most plant nutrients, growth of plants, and activity of microorganisms .
The higher total nitrogen level under the natural forest could be attributed to its better organic carbon content, which is the major source of total nitrogen (essentially organic nitrogen) . This is also supported by the highly significant (p < 0.01) and positive correlation (r = 0.96) between total nitrogen and organic carbon content (Table 6). Loss of considerable total nitrogen following conversion of land from forest to cultivated land was reported in many similar studies conducted elsewhere [49, 53, 61,62,63,64].
Some of the C:N ratios were outside the range that is considered normal (10–12) for arable mineral soils. As indicated by , ratios between 15 and 25 indicate a slow decomposition process due to scarcity of nitrogen. Accordingly, the ratios recorded in some surface and subsurface layers of the cultivated, grazing, and natural forest lands indicate the presence of slow decomposition process probably due to the presence of resistant structures and organic compounds. Strong and Mason  pointed out that organic matter with a high C:N ratio (> 20) locks up nitrogen as it decomposes, decreasing available nitrogen for the crop.
Soil carbon stock
The presence of high SOC in the upper 20 cm of the soil indicates the large amount of carbon to be lost if the upper 20 cm layer of the soil is disturbed, particularly in uncultivated lands. Besides the presence of high SOC in the upper 40 cm of coffee agroforestry, grazing, and crop land uses indicate the depth to that level is as important as the surface layers for storing carbon. On the other hand, the presence of 68–60% of the organic carbon in the lower (0.2–0.4 m and 0.4–0.6 m) soil layers indicates that the deeper layers of the soil are also important pools in terms of preserving soil organic carbon for a long time. Several other studies have also reported the presence of relatively higher organic carbon values in the upper layers of soils as compared with the lower or deeper layers [34, 67,68,69,70]. The lower SOC in the crop lands could be due to repeated cultivation before sowing , removal of crop residues [7, 14, 71], and burning of crop residues during land preparation . The current results indicate that the subsistence farming practices that have been followed by the smallholder farmers in the study area are exploitative in nature, and thus, have taken much of the soil organic carbon as compared to the natural forest. This loss of organic carbon is likely to increase emission from the soil system. A land use type that enhances the organic carbon content and minimizes emission has to be put in place in order to fight climate change and its unfavorable impacts. It is well understood that management practices have impact on the balance between input and output in a given system as well as on the rate of decomposition of organic matter . Therefore, appropriate land use and management practices [12, 14] such as conservation tillage , no tillage and fertilization/intensive agriculture , and agroforestry  should be employed to enhance the carbon sequestration potential of croplands under subsistence low input farming.
Even though the result of the study revealed higher SOC in grazing land compared with the national average of 26.2 t ha−1 , it is believed that excessive grazing pressure and the cut and carry system commonly used in the sub watershed have resulted in lower organic matter input into the system [75, 76]. The result also agreed with the findings of  who indicated the effect of grazing on plant biomass through removal of palatable plants which are high quality energy sources for decomposers. Enhancing carbon stock in grazing lands can be achieved through improved grazing management such as optimizing stock number, rotational grazing, and fertilization . However, the rotational grazing option is not feasible in the current context of the study area mainly due to population pressure and higher demand for agricultural land expansion. Hence, implementation of other options such as introduction of fodder species on marginal and other lands, and establishing fodder grasses on bunds are essential to minimize the grazing pressure and enhance the C sequestration potential of the grazing land.
The higher soil organic carbon stock in the natural forest might be attributed to the frequent addition of litter, the presence of network of roots, and modified microclimate, which retard decomposition rate of organic matter. The finding agrees with the findings of [78, 79] who reported higher SOC stock under forest than other land use types. Woldeamlak and Stroosnijder  also found high soil organic carbon in natural forest compared with grazing land, land under cultivation, and Eucalyptus plantation. Similarly,  reported high SOC in native forest compared with cultivated land and found a non-significant difference between native forest and grazing land. Noordwijk et al. [78, 79] also reported high SOC in disturbed forests compared with coffee based agroforestry.
According to information obtained from local agriculture office, most of the land that represents the different levels of the variable ‘land use’ used to be under natural forest cover many decades ago. Which was indeed attested by the presence of some remnants of naturally grown trees in the land under different land uses today. Hence, conversion from the natural system to human managed system has resulted in the loss of significant amount of SOC [9, 17,18,19,20, 75, 81, 82]. The declining trend in SOC could be due to repeated disturbance of the soil after conversion of forest land into other land uses [9, 11, 83].
Frequent addition of partially decomposed organic matter, presence of woody perennials, and well-built conservation structures have made coffee agroforestry to have higher SOC compared with grazing and crop lands. Reports  also witnessed the importance of management and biotic factors in influencing soil organic carbon at local level. The contribution of SOC to the total carbon stock of coffee agroforestry was 94.5%, which agrees with the findings of [78, 79, 85]. Nair et al.  indicated that time sequence studies on soil carbon are rare in agroforestry systems; hence the current study will serve as a bench mark for the study area.
The presence of SOC in lower soil depths indicates how these layers of the soil are important pools in terms of preserving soil organic carbon for a long time . Unlike studies elsewhere in Ethiopia [87,88,89,90,91], the current study yielded high SOC in the natural forest.
The relatively lower deterioration index under the coffee agroforestry could be associated with the relatively better soil management. Studies conducted elsewhere in Ethiopia (e.g., ) also indicated a high deterioration index under cultivated lands as compared with other land uses. These, suggesting that most of the smallholder subsistence farming practices in Ethiopia are highly exploitative and undermine the potential of the agriculture sector to sequester carbon and help in mitigating climate change. The results also imply that other alternative land uses with appropriate management strategies (e.g. climate-smart agriculture), which enhance the carbon stock and carbon sequestration potential of the lands while minimizing emissions, should be put in place to boost the subsistence farmers’ adaptive capacity against climate change.
Biomass carbon stock
The average AGC (116.46 t C ha−1) of the current finding is greater than tropical (45.45 t C ha−1) and regional (58.9 t C ha−1) estimates [43, 92]. On the other hand, all recently assessed forests in Ethiopia exhibited higher aboveground, root, and litter carbon stock compared to the finding of the current study [87,88,89,90,91]. The status of the forest and its composition of species, along with climatic and edaphic factors, might be responsible for the differences in aboveground, root, and litter carbon stock. Besides, the lower litter carbon in the study area could be related to the collection of twigs and branches for fuelwood by the local community. The lower carbon stock in the undergrowth vegetation of the natural forest could be due to the presence of large trees, though very small in number, whose shading effect may inhibit the growth of undergrowth vegetation. As put by , biomass of the understory declines as the forest age increases. The absence of undergrowth vegetation in coffee agroforestry sites is primarily due to intensive weeding and soil management works.
The coffee agroforestry had lower vegetation carbon stock compared with similar studies in Guatimala  and Costa Rica . The lower carbon stock in the current study could be attributed to the presence of a very small number of trees with small diameter. The aboveground carbon stock in agroforestry relate with the density of coffee shrub and shade trees , type of shade trees used [34, 94], and size of the plant (current study). Hence, the combination of these three criteria provides reasonable outcome. Nonetheless, it is also important to increase diversity of species, most importantly woody perennials, to enhance the service  and productive roles [95, 96] of coffee agroforestry in the study area.
The study demonstrates that even under poor management, soil has higher carbon stock compared with aboveground biomass. The subsistence farming practice followed by the smallholder farmers in the study area is exploitative in nature. As a result, it has taken much of the soil organic carbon, which was found to be relatively higher in the natural forest. This loss of organic carbon is likely to increase emission from the soil system. Also, the climax reaching vegetation and its impacts on undergrowth plants render the natural forest low carbon sequestration potential. Hence, it is necessary to implement conservation-based production systems with inclusion of woody perennials. Although there is scarcity of information on carbon stock under different land uses in Ethiopia, the present study highlights the importance of watershed level carbon stock assessment for better and carbon friendly land use decision making. The study, however, did not attempt to assess the carbon gain and loss as a function of each land use, so further study is needed to clearly map carbon sequestration potential of the sub-watershed.
above ground biomass
green house gas
soil organic carbon
Lal R. Managing soils and ecosystems for mitigating anthropogenic carbon emissions and advancing global food security. Bioscience. 2010;60:708–21.
IPCC. Agriculture in climate change mitigation. Working group III contribution to the fourth assessment report of the international panel on climate change. New York: Cambridge University Press; 2007.
International Energy Agency (IEA). CO2 emission from fuel combustion, highlight; 2015. http://www.iea.org/t&c. Accessed 23 Nov 2016.
Lal R. Challenges and opportunities in soil organic matter research. Eur J Soil Sci. 2009;60:158–69.
Steen M. Greenhouse gas emissions from fossil fuel fired power generation systems. Europian Union Joint Research Center, EUR 19754 EN; 2000. http://www.jrc.nl. Accessed 14 Feb 2018.
Yihenew GS, Getachew A. Effects of different land use systems on selected physicochemical properties of soils in North western Ethiopia. J Agric Sci. 2013;5(4):112–20.
Berry N. Whole farm Carbon accounting by smallholder, lesson from plan Vivo project presentation at the small-scale mitigation whole farm and landscape accounting workshop. FAO, Rome, 27–28th October, 2011. www.fao.org/climatechnage/micca/72531/en/.
Girmay G, Singh BR, Mitiku H, Borresen T, Lal R. Carbon stocks in Ethiopian soils in relation to land use and soil management. J Land Degrad Dev. 2008;19(4):351–67.
Lemenih M, Karltun E, Olsson M. Soil organic matter dynamics after deforestation along a farm field chronosequence in southern highlands of Ethiopia. Agr Ecosyst Environ. 2005;109:9–19.
Melero S, López BR, López BL, Muñoz RV, Moreno F, Murillo JM. Long-term effect of tillage, rotation and nitrogen fertilizer on soil quality in a mediterranean vertisol. Soil Tillage Res. 2011;114:97–107.
Yeshanew A, Zech W, Guggenberger G, Tekalign M. Soil aggregation, and total and particulate organic matter following conversion of native forests to continuous cultivation in Ethiopia. Soil Tillage Res. 2007;94:101–8.
Lal R. Enhancing crop yields in the developing countries through restoration of the soil organic carbon pool in agricultural lands. Land Degrad Dev. 2006;17:197–209.
Lipper L, Neves B, Wilkes A, Tennigketi T, Gerber P, Henderson B, Branca G, Mann W. Climate change mitigation finance for smallholder agriculture: a guide book to harvesting soil carbon sequestration benefit. Rome: FAO; 2011.
Pan G, Smith P, Pan W. The role of soil organic matter in maintaining the productivity and yield stability of cereals in China. Agric Ecosyst Environ. 2009;129:344–8.
Wright AL, Hons FM. Carbon and nitrogen sequestration and soil aggregation under sorghum cropping sequences. Biol Fertil Soils. 2005;41:95–100.
Don A, Schumacher J, Freibauer A. Impact of tropical land-use change on soil organic carbon stocks a meta-analysis. Glob Change Biol. 2011;17:1658–70.
International Panel for Climate Change (IPCC). Climate change. The physical science basis; 2013. http://www.climatechange2013.org/images/report/WG1AR5_ALL_FINAL.pdf. Accessed 18 Feb 2016.
Janzen HH. Carbon cycling in earth systems—a soil science perspective. Agr Ecosyst Environ. 2004;104:399–417.
Lasco RD. Forest carbon budgets in Southeast Asia following harvesting and land cover change. Sci China (Series C). 2002;45:55–64.
Vagen TG, Lal R, Singh BR. Soil carbon sequestration in sub-Saharan Africa: a review. Land Degrad Dev. 2005;16:53–71.
Dawson JJ, Smith P. Carbon losses from soil and its consequence for land-use management. Sci Total Environ. 1997;382:165–90.
Lal R. Forest soils and carbon sequestration. For Ecol Manag. 2005;220:242–58.
Guibin W, Welham C, Feng C, C L, Cao Fuliang. Enhanced soil carbon storage under agroforestry and afforestation in Subtropical China. Forests. 2015;6:2307–23.
Denu D, Platts PJ, Kelbessa E, Gole TW, Marchant R. The role of traditional coffee management in forest conservation and carbon storage in the Jimma Highlands, Ethiopia. For Trees Livelihoods. 2016;25(4):226–38.
Freibauer A, Rounsevell M, Smith P, Verhagen A. Carbon sequestration in European agricultural soils. In: Pete Smith, 2004. Carbon sequestration in croplands: the potential in Europe and the global. Eur J Agronomy. 2004;20:229–36.
Negasa B, Zewdu T, Gessa A. Carbon sequestration potentials of semi-arid rangelands under traditional management practices in Borana, Southern Ethiopia. Agr Ecosyst Environ. 2016;223:108–14.
Food and Agriculture Organization of the United Nations/International Institute for Applied Systems Analysis/International Soil Reference and Information Centre/Institute of Soil Science-Chinese Academy of Sciences/Joint Research Centre of the European Commission (FAO/IIASA/ISRIC/ISSCAS/JRC). Harmonized World Soil Database (version 1.2), Rome; 2012.
Hairiah K, Sitompul SM, Noordwijk MV, Palm, C. Methods for sampling carbon stocks above and below ground; 2001. http://www.icraf.cgiar.org/sea.
Iqbala S, Tiwari SC. Altitudinal variation of soil organic carbon stock in Achanakmar. Int J Curr Res. 2015;7(05):15885–90.
Pearson TR, Brown SL, Birdsey RA. Measurement guidelines for the sequestration of forest carbon. Washington, D.C: USDA Forest Service; 2007. https://doi.org/10.2737/NRS-GTR-18.
Negash M, Starr M, Kanninen M, Berhe L. Allometric equations for estimating aboveground biomass of Coffea arabica L. grown in the Rift Valley escarpment of Ethiopia. Agrofor Syst. 2013;87:953–66.
Santantonio D, Hermann RK, Overton WS. Root biomass studies in forest ecosystems. Pedobiologia. 1997;17:1–31.
Mulugeta L, Itanna F. Soil carbon stocks and turnovers in various vegetation types and arable lands along an elevation gradient in southern Ethiopia. Geoderma. 2004;123:177–88.
Ehrenbergerová L, Cienciala E, Kučera A, Guy L, Habrová H. Carbon stock in agroforestry coffee plantations with different shade trees in Villa Rica, Peru. Agrofor Syst. 2016;90:433–45.
Abera G, Wolde-Meskel E. Soil properties, and soil organic carbon stocks of tropical andosol under different land uses. Open J Soil Sci. 2013;3:153–62.
Noponen MR, Haggar JP, Edwards-Jones G, Healey JR. Intensification of coffee systems can increase the effectiveness of REDD mechanisms. Agric Syst. 2013;119:1–9.
Campbell CR, Plank CO. Preparation of plant tissue for laboratory analysis. In: Kalra YP, editor. Handbook of reference methods for plant analysis. Boca Raton: CRC Press; 1998. p. 37–49.
Van Reeuwijk LP. Procedures for soil analysis. Technical paper 9, 6th edition. International Soil Reference and Research Center; 2002.
Blake GR, Hartge KH. Bulk density. In: Klute A, editor. Methods of soil analysis. Part 1: physical and mineralogical properties. 2nd ed. Wisconsin: American Society of Agronomy; 1986. p. 363–82.
McLean E. Soil pH and lime requirement. In: Page AL, et al., editors. Methods of soil analysis. Part 2: chemical and microbiological properties. 2nd ed. Wisconsin: Soil Science Society of America Inc.; 1982. p. 199–224.
Walkley A, Black IA. An examination of the Degtjareff method for determining soil organic matter and a proposed modification of the chronic acid titration method. Soil Sci. 1934;37:29–38.
Bremner JM, Mulvaney CS. Nitrogen total in methods of soil analysis. In: Page AL, editor. Method of soil analysis part 2. Madison: Soil Science Society of America; 1982.
Brown S. Estimating biomass and biomass change of tropical forests: a primer, vol. 134. Rome: Food and Agriculture Organization of the United Nations; 1997. p. 20–55.
Cairns MA, Brown S, Helmer E, Baumgardener GA. Root biomass allocation in the world’s upland forests. Oecologia. 1997;111:1–11.
Doetterl S, Kearsley E, Bauters M, Hufkens K, Lisingo J, Baert G, et al. Aboveground vs belowground carbon stocks in african tropical lowland rainforest: drivers and implications. PLoS ONE. 2015;10(11):e0143209. https://doi.org/10.1371/journal.pone.0143209.
Armeein RB, Gabon FM. Biomass, organic carbon and mineral matter contents of abaca (Musa textilis Nee) at different stages of growth. Ind Crops Prod. 2008;28:340–5.
Zhang Y, Zhao YC, Shi XZ, Lu XX, Yu DS, Wang HJ, Sun WX, Darilek JL. Variation of soil organic carbon estimates in mountain regions: a case study form Southwest China. Geoderma. 2008;146:449–56.
Adejuwon JO, Ekanade O. A comparison of soil properties under different land use types in a part of the Nigerian cocoa belt. Catena. 1988;15:319–31.
Abate N, Kibret K. Effects of land use, soil depth and topography on soil physicochemical properties along the toposequence at the Wadla Delanta Massif, Northcentral Highlands of Ethiopia. Environ Pollut. 2016;5:2.
Ellerbrock R, Gerke H. Characterization of soil organic matter composition in top and sub soil samples from colluvial and eroded sites. ZALF, Inst. of soil landscape Research, Muncheberg, Germany. Geophysical Research Abstracts. 2013; 15.
Fantaw Y, Ledin S, A Abdu. Changes in soil organic carbon and total nitrogen contents in three adjacent land use types in the Bale Mountains, south-eastern highlands of Ethiopia. For Ecol Manag. 2007;242:337–42.
Nigussie T. Land use/land cover dynamics, soil characterization, and GIS based multi-criteria decision method land suitability evaluation in Wachacha Watershed, central highland of Ethiopia. Doctoral dissertation, Haramaya University; 2017. p. 169.
Gebrekidan H, Yitbarek T, Kibret K, Beyene S. Impacts of land use on selected physicochemical properties of soils of Abobo area, western Ethiopia. Agric For Fish. 2013;2(5):177–83.
Abate N, Kibret K, Gebrekidan H, Esayas A. Characterization and classification of soils along the toposequence at the Wadla Delanta Massif, north central highlands of Ethiopia. J Ecol Nat Environ. 2014;6(9):304–20.
Feyissa S. Characterization, classification, and spatial variability of soil fertility status under different land use/land cover in Fugnan Dimo Watershed, Western Hararghe of Oromia Region. Doctoral dissertation, Haramaya University; 2017. p. 186.
Jones CA. Effect of soil texture on critical bulk densities for root growth. Soil Sci Soc Am J. 1983;47:1208–11.
Emiru Nega, Gebrekidan Heluf. Effect of land use changes and soil depth on soil organic matter, total nitrogen and available phosphorus contents of soils in Senbat watershed, Western Ethiopia. ARPN J of Agric Biol Sci. 2013;8(3):206–12.
Soto B, Diazfierroz F. Interaction between plant ash leachates and soil. Int J Wildland Fire. 1993;3:207–16. https://doi.org/10.1071/WF9930207.
Tekalign T. Soil, plant, water, fertilizer, animal manure and compost analysis. Working Document No. 13. Addis Ababa: International Livestock Research Center for Africa; 1991.
Landon JR, editor. Booker tropical soil manual: a handbook for soil survey and agricultural land evaluation in the tropics and subtropics. Long man scientific and technical. New York: Booker Tate Ltd., Wiley; 2014.
Eyayu M, Heluf G, Tekaliign M, Mohammed A. Effects of land use change on selected soilproperties in the TeraGedam Catchment and adjacent agro-ecosystems, Northwest Ethiopia. Ethiop J Nat Resour. 2009;11(1):35–62.
Mojiri A, Aziz HA, Ramaji A. Potential decline in soil quality attributes as a result of land use change in a hill slope in Lordegan, Western Iran. Afr J Agric Res. 2012;7(4):577–82.
Mulugeta L, Karltun E, Olsson M. Assessing soil chemical and physical property responses to deforestation and subsequent cultivation in smallholders farming system in Ethiopia. Agric Ecosyst Environ. 2005;105:373–86.
Taye K. Chemical properties of wild coffee forest soils in Ethiopia and management implications. Agric Sci. 2011;2(4):443–50.
Hazelton P, Murphy B. Interpreting soil test results: what do all the numbers mean?. 2nd ed. Sydney: CSIRO Publishing; 2007.
Strong WM, Mason MG. Nitrogen. In: Peverill KI, Sparrow LA, Reuter DJ, editors. Soil analysis: an interpretation manual. Melbourne: CSIRO Publishing; 1999.
Grüneberg E, Schöning I, Kalko EK, Weisser WW. Regional organic carbon stock variability: a comparison between depth increments and soil horizons. Geoderma. 2010;155:426–33.
Sharma V, Hussain S, Sharma KR, Arya VM. Labile carbon pools and soil organic carbon stocks in the foothill Himalayas under different land use systems. Geoderma. 2014;232–234:81–7.
Wang Y, Fu B, Lü Y, Song C, Luan Y. Local-scale spatial variability of soil organic carbon and its stock in the hilly area of the Loess Plateau, China. Quat Res. 2010;73:70–6.
Zádorová T, Penížek V, Vašát R, Žížala D, Chuman T, Vaněk A. Colluvial soils as a soil organic carbon pool in different soil regions. Geoderma. 2015;253–254:122–34.
World Bank. Sustainable land management project, project appraisal document (PAD), Ethiopia/Report No. 42927-ET, Project I.D P107139; 2008. http://www-wds.worldbank.org/external/projects/. Accessed 14 Feb 2018.
Lal R, Kimble JM, Follet R. Land use and soil carbon pools in terrestrial ecosystems. In: Lal R, Kimble JM, Follet R, editors. Management of carbon sequestration in soils. New York: CRC Press; 1997.
Nair PKR, Nair VD, Kumar BM, Solomon GH. Soil carbon sequestration in tropical agroforestry systems: a feasibility appraisal. Environ Sci Policy. 2009;12:1099–111.
Federal Democratic Republic of Ethiopia (FDRE). Ethiopia—land degradation neutrality national report. Addis Ababa; 2015.
Nyssen J, Temesgen H, Mulugeta L, Zenebe A, Haregeweyn N, Mitiku H. Spatial and temporal variation of soil organic carbon stocks in a lake retreat area of the Ethiopian Rift Valley. Geoderma. 2008;146:261–8.
Riedo M, Gyalistras D, Fuhrer J. Net primary production and carbon stock in differently managed grasslands: simulation of site-specification sensitivity to an increase in temperature CO2 and to climate change. Ecol Model. 2000;134(2–3):207–27.
Guo LB, Gifford RM. Soil carbon stocks and land use change: a meta-analysis. Glob Change Biol. 2002;8:345–60.
Noordwijk M, Rahayu S, Hairiah K, Wulan YC, Farida A, Verbist B. Carbon stock assessment for a forest-to coffee conversion landscape in Sumber-Jaya (Lampung, Indonesia): from allometric equations to land use change analysis. Sci China. 2002;45:75–86.
Schmitt-Harsh M, Castellanos EJ, Evans TP, Randolph JC. Carbon stock in coffee agroforests and mixed dry tropical forests in the western highlands of Guatemala. Agrofor Syst. 2012;86:141–57.
Bewket W, Stroosnijder L. Effects of agroecological land use succession on soil properties in Chemoga watershed, Blue Nile basin, Ethiopia. Geoderma. 2003;111:85–98.
Houghton RA, Hackler JL. Changes in terrestrial carbon storage in the United States: the roles of agriculture and forestry. Glob Ecol Biogeogr. 2000;9:125–44.
Solomon D, Fritzche F, Lehmann J, Tekalign M, Zech W. Soil organic matter dynamics in the sub-humid agroecosystems of the Ethiopian highlands: evidence from natural 13C abundance and particle-size fraction. Soil Sci Soc Am J. 2002;66:969–78.
Shrestha RK, Ladha JK, Gami SK. Total and organic soil carbon in cropping systems of Nepal. Nutr Cycl Agroecosyst. 2006;75:257–69.
Allen DE, Pringle MJ, Page KL, Dalal RC. A review of sampling designs for the measurement of soil organic carbon in Australian grazing lands. Rangelands J. 2010;32:227–46.
Häger A. The effects of management and plant diversity on carbon storage in coffee agroforestry systems in Costa Rica. Agrofor Syst. 2012;86:159–74.
Harper RJ, Tibbett M. The hidden organic carbon in deep mineral soils. Plant Soil. 2013;368:641–8.
Feyisa A, Sormessa T, Argaw M. Forest carbon stocks and variation along altitudinal gradients in Egdu forest: implication of managing forest for climate change mitigation. Sci Technol Arts Res J. 2013;2(4):40–6.
Yohannes H, Soromessa T, Argaw M. Carbon stock analysis along altitudinal gradient in Gedo forest: implications for forest management and climate change mitigation. Am J Environ Prot. 2015;4(5):237–44.
Gedefaw M, Soromessa T, Belliethathan S. Forest carbon stock in woody plants of Tara Gedam forest: implication for climate change mitigation. Sci Technol Arts Res J. 2014;3(1):101–7.
Nega M, Sormessa T, Bayable E. Above and below ground reserved carbon in Danaba community forest of Oromiya Region, Ethiopia: implication for CO2 emission balance. Am J Environ Prot. 2015;4(2):75–82.
Tura T, Argaw M, Eshetu Z. Estimation of carbon stock in church forests: implications for managing church forest to help with carbon emission reduction. In: Leal Filho W, Mannke F, Mohee R, Schulte V, Surroop D, editors. Climate-smart technologies. Berlin: Climate Change Management; 2013. p. 403–14.
Gibbs HK, Brown S. Geographical distribution of woody biomass carbon stocks in tropical Africa: an updated database for 2000; 2007. https://doi.org/10.3334/cdiac/lue.ndp055.2007.
Pibumrung P, Gajaseni N, Popan A. Profiles of carbon stocks in forest, reforestation and agricultural land, Northern Thailand. J For Res. 2008;19(1):11–8.
Demessie A, Singh BR, Lal R, Børresen T. Leaf litter fall and litter decomposition under Eucalyptus and coniferous plantations in Gambo District, Southern Ethiopia. Acta Agric Scand Sect B Plant Soil Sci. 2012;62:467–76.
Rice RA, Ward J. Agricultural intensification within agroforestry: the case of coffee and wood products. Agric Ecosyst Environ. 2008;128:212–8.
Schroth G, D’Angelo SA, Geraldes W, Haag D, Lieberei R, et al. Conversion of secondary forest into agroforestry and monoculture plantations in Amazonia: consequences for biomass litter and soil carbon stocks after 7 years. For Ecol Manag. 2002;163:131–50.
This research was funded by Haramaya University, Vice president Office for research affairs.
The fund for this research was obtained from Haramaya University.
Availability of data and materials
Data is available in Springer nature.
Consent for publication
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.