Close Home
iForest - Biogeosciences and Forestry
vol. 9, pp. 626-634
Copyright © 2016 by the Italian Society of Silviculture and Forest Ecology
doi: 10.3832/ifor1866-008

Collection: IUFRO RG7.01.00 - Nice (France 2015)
“Global Challenges of Air Pollution and Climate Change to Forests”
Guest Editors: Elena Paoletti, Pierre Sicard

Research Articles

Carbon storage and soil property changes following afforestation in mountain ecosystems of the Western Rhodopes, Bulgaria

Miglena Zhiyanski (1)Corresponding author, Maria Glushkova (2), Angel Ferezliev (3), Lorenzo Menichetti (4), Jens Leifeld (4)


Mountain ecosystems are important for their distinct ecological features and the ecosystem services they offer such as regulating atmospheric composition and providing climate change mitigation via carbon (C) sequestration ([42], [64]). Under contemporary global change with widespread abandonment of agricultural lands and massive afforestation activities mostly in the last century, as well as the expected increase in summer temperatures and shifts in precipitation patterns, C storage patterns likely will change in various compartments of mountain ecosystems. Spatial variability of organic matter stocks and turnover in different types of ecosystem depend on many factors (air temperature, aboveground vegetation composition, soil type, soil moisture regime etc. - [27], [37]). Mountain ecosystems are highly sensitive to environmental changes ([32]). A recent report in Europe highlights the fragility and vulnerability of mountain areas ([13]), and projects that both climate and land-use changes (LUC) will have profound effects on their services and functioning. In addition, climate change impacts interact with LUC in mountains ([57], [59]). These changes are expected to affect many services that mountains provide, including C sequestration ([51], [36], [37]).

Environmental degradation, including soil C loss caused by inappropriate land use, is a worldwide problem that has fostered research into sustainable production systems ([75], [19], [62], [29]). Land-use changes influence the balance of soil organic carbon (SOC) and hence may cause CO2 emissions or C sequestration ([47]). In recent years, numerous studies have focused on SOC changes in various ecosystems ([50], [30], [67]). Nevertheless, Poeplau & Don ([47]), showed that in Europe different forms of land use and LUC can lead to either C loss or C accumulation in soils, and emphasized the lack of studies covering all major LUC and biomes with respect to qualitative and quantitative LUC effects.

Soil quality plays an important role in sustainable land management and for environmental integrity. During the recent decades, soil quality concepts have emerged and have been used to assess land or soil quality for various ecosystems ([28], [41]). Changes in land use and land management have great effects on the direction and degree of changes in soil physical and chemical properties ([22]). Land-use change obviously affects the cycles of C ([21], [72]) and nitrogen (N - [49]) but also of soil potassium and phosphorus. Nitrogen is a key factor that influences the primary production of ecosystems ([45]) but, on the other hand, is also controlled by land use and management practices ([69], [68]). With continuous cultivation, physical soil properties may change and organic matter (OM) content, an important carrier of soil quality and storage compound for soil N, may decline ([38]). An important factor for changes in OM content is the flux of plant residues, which is smaller in cropland compared with grassland or forest ([2], [54]). Organic matter input rates also differ among land-use types owing to the different share of root-derived residues ([26], [66]). Many studies have reported significant effects of LUC from cropland to forest and vice versa on soil C ([48], [18]). Conversion of a native ecosystem to cultivated land may cause changes in some soil properties, and often, long-term cropland management induces decreasing soil C stocks. The latter are caused not only by smaller residue inputs but also by less protection of soil C owing to physical disintegration of soil ([9], [5]).

Soil bulk density is the major discriminating factor in soil properties in relation to land use and management ([60], [70], [44]). Soil bulk density integrates effects of management and OM decline; these factors in conjunction result in smaller pore volumes and thus higher bulk densities. Celik ([10]), for instance, showed that soil bulk density increased with cultivation relative to that of adjacent soils in forests and pastures in the southern Mediterranean highlands of Turkey.

In Europe, grasslands at relatively low elevations and flat topography are more and more intensively managed, whereas less fertile pastures at higher altitudes are increasingly abandoned or actively converted into forests ([12]). The area and management of mountain pastures are changing rapidly at national scale ([15]). Despite this rapid change, and in contrast to the vast literature on agricultural soils, only few studies have examined the effects of the long-term conversion from grassland or cropland on soil properties in mountain forest ecosystems ([56]). Historical review on land-use in the mountain zones in Bulgaria showed that past forests have been threatened by excessive and improper use by the local population, which has led to replacement of typical tree species by others ([43], [1]). In Bulgaria, afforestation is of high importance for mountain ecosystems, and substantial land conversion took place during the last decades to contrast erosion, to increase the treeline zone and to use low-productive areas.

Our objective was to characterize soil properties and characteristics of aboveground vegetation in ecosystems of the Western Rhodope Mountains, Bulgaria. We studied the effect of typical regional afforestation practices on soil physical and chemical properties and organic C stocks to evaluate which land-use types are most valuable for ecosystem C storage and, hence, climate change mitigation.

Materials and methods 

Study area

The study area was located in the temperate climatic zone and Trace forest region, sub-zone of the Western Rhodopes ([52]). The average annual precipitation in the studied region was 630.5 mm with peaks from May to July, and the mean annual temperature was 7.2 °C. A specific feature of the regional climate is the period of increasing temperatures, which starts from early October to mid-November. The studied sites were located at 860-1270 m a.s.l. in the mountain climatic zone in Bulgaria ([52]).

The study was conducted in 2013/2014 at four experimental sites, representing the main forest-related conversions in this region. The four land-use types included in the study were Rd1, a Douglas fir (Pseudotsuga menziesii [Mirb.] Franco) plantation created on former cropland RdA1; Rd2, a mixed forest plantation dominated by black pine (Pinus nigra Arn.) in association with Scots pine (Pinus sylvestris L.), created on abandoned agricultural land RdA2. At both forest sites, one experimental plot with a total area of 0.1 ha (25 × 40 m) was randomly selected near the adjacent site with maintained previous land use.

Site Rd1 was representative of Douglas fir plantations involving Scots pine and was, together with RdA1, located in the area of St. Petka, Velingrad. The plantation Rd1 was established in 1968 through afforestation in rows on existing cropland. Site RdA1 represented the cropland adjacent to Rd1, and was used for producing annual crops (vegetables and cereals including potatoes, beans, peas, corn, fodder beet, beetroot, rye and oats).

Site Rd2 was part of the area afforested with Scots and black pine plantations before the Second World War in the region of Velingrad for prevention of ongoing soil erosion ([55]). It extended over 2700 ha. The experimental area was located in the locality Baba and Diado on a southern slope. Site Rd2 was the LUC of RdA2. Site RdA2 represented abandoned land, located near Rd2 in the upper part of the south-facing slope. Currently, the area is used occasionally for grazing animals (horses, sheep and goats). Until 80 years ago, RdA2 was managed as a cropland.

Details of the selected sites are given in Tab. 1. Afforestation at both forest sites was performed with 2-year-old seedlings with open root systems and planted with the “Sward of Kolesov”.

Tab. 1 - Characteristics of the experimental sites. (FL): forest land; (CL): cropland; (AL): abandoned land. (*): Age of stand and land-use history (2014); (**): IUSS WG ([25]).

Soil sampling and laboratory measurements

Sites were sampled during the growing season between April and October in 2013 and 2014. Soils were sampled randomly in five repetitions within a 0.1 ha experimental plot for all studied land-use types (Tab. 1). All these sites presented similar slope and soil type within the sampled site. Soils were sampled by a hand auger (4 cm in diameter) in four soil layers: 0-10, 10-20, 20-30 and 30-40 cm (20-40 cm depth was analyzed as an aggregated sample of 20-30 and 30-40 cm). In total, 80 soil samples were collected.

Soil properties were determined in accordance with the standardized methods in the Laboratory of Forest Soil Science at the Forest Research Institute - BAS ([11]). Recorded data included morphological field description, bulk density, particle size distribution (wet sieving, including coarse fragments, i.e., particles with a diameter of 2-10 mm), textural classes, pH in H2O ([24]), N content and organic C content. The bulk density was measured per each soil layer of undisturbed soil with cylinders with known volume in two repetitions. Fresh soil samples were individually weighed and dried for 48 h at 105 °C. Each soil sample was sieved in a 2 mm mesh to remove coarse sands while soil aggregates were broken. In the next step, the sample was ground and the total organic C content was measured by the modified method of Tjurin ([3]). Total N content was measured according to the method of Kjeldhal ([7]). The C content, combined with bulk density, was used to estimate the amount of carbon per unit area.

Forest floor (at afforested sites Rd1 and Rd2) and dead litter (RdA1 and RdA2) were sampled within each experimental plot using a 25 × 25 cm frame in three repetitions. Aboveground biomass at RdA1 and RdA2 formed by grasses and crop species was harvested once in the growing season, and biomass was measured by 1 × 1 m frames in three repetitions. All grass species were collected and their fresh weight was measured. Then the grasses were dried in an oven for 48 h at 105 °C and weighed. The samples from forest floor layers and grass species were separated and their dry mass (after 48 h at 105 °C) was measured individually; then each sample was ground. The C content in each layer (L, F and H) was estimated by multiplying the dry weight by the C concentration in dry matter and converted to t C ha-1.

Tree stand characteristics were measured using Bulgarian biometrical methods ([34]). At both forest sites, number of trees, mean diameter at breast height, mean height and growing stock were determined and calculated per hectare. The volumes of biomass were calculated based on tree species tables for trees with bark for the region of the Western Rhodopes ([33], [46]) (Tab. 2), adjusted for stand density. For the present study, we applied IPCC reference values ([23]) for wood dry matter density and C content (1 m3 of wood contains a dry matter of 500 kg with a C fraction of dry matter of 0.5).

Tab. 2 - Biometric characteristics of trees at the studied forest sites (mean ± SE).

The organic C stocks in soil and litter were quantified according to the IPCC Good Practice Guidance for Land Use, Land Use Change and Forestry ([23]). The SOC stock at each site was based on equivalent soil masses ([14], [65]).

Statistical analysis

Statistics were run using the SPSS 16.0 program. Departure from normal distribution of each soil variable was tested by χ2 normality test. An analysis of variance (ANOVA) was carried out, and differences of means of studied soil properties among depth layers per site and between study sites were verified by Tukey’s test (α = 0.05). Pearson’s linear correlation coefficients among studied parameters were calculated ([58]) and used to reveal the direction of relationships between selected soil properties and among study sites and soil depths.

Results and discussion 

Soil bulk density

Bulk density at 0-10 cm depth in the studied land-use types ranged from 0.82 to 0.93 g cm-3. Land use had no significant effect (P = 0.546), although higher densities occurred in cropland (RdA1) and abandoned land (RdA2), and lower densities in the forest stands. In contrast, at 20-40 cm soil depth, bulk densities were slightly higher in the forest stands (0.93 and 0.92 g cm-3) than in abandoned land (0.90 g cm-3 - Tab. 3). Although these values were not significantly different, the relatively higher bulk density in the upper layer of the cropland soil requires additional studies concerning probable indication for soil compaction due to more intensive cultivation by traditional farming and poor crop rotation practices in the region, i.e., a dominance of vegetables and missing cover crops. Some previous studies also indicated that soil cultivation significantly increases soil bulk density ([40]).

Tab. 3 - Main soil and litter/forest floor parameters at the experimental sites (mean ± SE). (DW): Dry weight of forest floor and dead litter.

Soil coarse fragment distribution

The coarse fragment distribution ranged from 13.2 to 28.6% in different layers of the four land-use types (Tab. 3), but the variability was not related to land use (P = 0.077). The coarse fraction tended to be highest in soil of abandoned land RdA2 (23.0-28.6%), probably because abandonment can result in increased surface stoniness ([31]). This effect could also be related to erosion events and preferential loss of particles <2 mm during previous intensive land management. Under forest vegetation, coarse fragment contents at 0-10 cm depth tended to be low (13.2 and 21.7%). This implies abundant water-stable aggregates that may reduce the soil vulnerability to erosion. For cropland RdA1, the fraction of coarse fragments of between 18.9 and 26.8% indicated that soil in the upper layer was relatively little influenced by erosion but part of the silt fraction may have been lost. Generally, the coarse fragment content in the studied area was low, suggesting good conditions for seedling emergence.

Sand, silt and clay contents

Mean sand content differed significantly between sites (F = 6.050; P = 0.019). Both the highest and the lowest sand contents occurred in the forests (Tab. 3). Considering the three soil depths among sites, sand content was highest at 20-40 cm depth. For all land-use types except Rd2, sand content increased with soil depth. Soil erosion can modify soil properties by reducing soil depth, changing soil texture and inducing loss of nutrients and OM ([16]). Soil silt and clay contents differed significantly between sites (F = 5.943, P = 0.020). This could be due to differences in soil formation but may also be caused by a preferential loss of fine fractions upon erosion ([74]). Sand content increases upon soil erosion after forest conversion and, hence, can be used as an indicator for evaluating soil degradation ([4]). The textural data tentatively suggested that erosion processes at RdA1 were still ongoing. Opposite to sand content, the fraction formed by silt and clay was largest in the superficial soil with highest values for Rd1 and RdA2 (Tab. 3). Plantation Rd2 had the lowest content of medium and small particles (silt and clay). Among sites, sand content was highest in the surface layer of the mixed forest plantation whereas clay content was highest in the surface layer of the RdA2. Unlike the clay fractions in forest sites and abandoned land, that in cropland RdA1 did not change with depth. Based on the data of textural composition, the studied soils were classified as sandy-loams.

Soil pH

Soil pH values were significantly affected by both land use (F = 83.516, P < 0.001) and soil depth (F = 57.693, P < 0.001). Generally, the pH values at the four sites were within the range of moderately acidic to neutral soil reactions as indicated by Foth & Ellis ([17]). Considering the three soil depths, pH was highest in surface layers of RdA1, followed by RdA2 (Fig. 1). By contrast, forest soils were acidified in their upper layers, most likely due to differences in the quality of litter inputs and the acidifying effect of coniferous materials ([6]). Moreover, soil pH at agricultural sites could be increased by management practices. Together, differences in land use came along with systematic and significant differences in soil physical and chemical properties, especially in the superficial soil layer.

Fig. 1 - Soil pH in soil layers. Different letters indicate significant differences among means after Tukey’s test (P ≤ 0.05). Error bars represent the standard error of the mean values per layer for each land use. Land uses compared: (Rd1) Douglas fir forest plantation; (RdA1) Arable cropland; (Rd2) Mixed P. nigra and P. sylvestris forest plantation; (RdA2) Abandoned land.

Carbon content

The average organic C content in mineral soil layers ranged from 0.35 to 4.14% and differed significantly between sites (F = 9.157, P = 0.006). At 0-10 cm depth, the highest C content was found at Rd1 and the lowest at RdA2 (Fig. 2). The highest content in the forest was most likely related to the presence of large biomass and deep-rooting plants that produced root input and prevented easy removal of nutrients from the soil. On the other hand, with conversion from cropland (RdA1) to forest (Rd1), soil C in the 0-10 cm soil layer increased by 3.7% after 40 years. The low C content in cropland RdA1 likely was related to the continuous cultivation of selected crops with high C exports, a factor known to be important for OM decline ([54]). Sites Rd2 (mixed forest land) and RdA2 (abandoned land) had even lower C contents than site RdA1, probably due to previous long-term arable cultivation (>100 years). The forest vegetation at Rd2 caused a small increase in C content at 0-10 cm depth relative to the abandoned land site. An SOC inventory in southeast Germany revealed only slightly lower total SOC stocks in cropland soils compared with forest soils when both top- and subsoils were considered ([67]). In the present study, ANOVA indicated significant differences in C content between the three soil layers (F = 11.046, P < 0.001) for the three land-use types with permanent vegetation, but not for arable cropland (RdA1).

Fig. 2 - Carbon content (%) in soil layers at the study sites. Different letters indicate significant differences among means after Tukey’s test (P ≤ 0.05). Error bars represent the standard error of the mean values per layer for each land use. Land uses compared: (Rd1) Douglas fir forest plantation; (RdA1) Arable cropland; (Rd2) Mixed P. nigra and P. sylvestris forest plantation; (RdA2) Abandoned land.

Nitrogen content

The N content varied between 0.04 and 0.79%, detected at Rd2 (20-40 cm) and Rd1 (0-10 cm), respectively (Fig. 3). Both land use (F = 59.676, P < 0.001) and soil depth (F = 8.225, P = 0.003) were significant factors affecting soil N. The highest N content at 0-10 cm depth was found at Rd1 and the lowest at RdA2. Average N contents at 0-40 cm depth differed between sites Rd1 and RdA1 and between sites Rd2 and RdA2 (P ≤ 0.05 - Tab. 4). It is important to consider that cropland soil generally contains little OM, and thus little organically bound nitrogen, because OM is removed during harvesting. Furthermore, N fertilization of croplands helps to maintain a higher level of soil N as compared with unfertilized fields ([35]), although much of the added nitrogen is exported by harvest. The occasional inputs of fertilizers at study site RdA1 likely did not change the soil total N content. Long-term cultivation and harvesting at this site resulted in smaller SOC and N contents compared with Rd1, despite additional organic fertilization. Moreover, surrounding forest vegetation significantly affected soil N content, with higher values obtained at Douglas fir plantation Rd1 and mixed pine plantation Rd2 than at cropland RdA1 and abandoned land RdA2. This result confirms the findings of Holubík et al. ([20]) at sites in the Czech Republic, where soil development resulted in a higher production of C and N pools in the forest stands compared with agricultural land. Hence, we interpret the different soil N contents largely as being a result of continuous N removal from cropland that was not fully compensated by fertilizer additions.

Fig. 3 - Nitrogen content (%) in soil layers at the study sites. Different letters indicate significant differences among means after Tukey’s test ( P ≤ 0.05). Error bars represent the standard error of the mean values per layer for each land use. Land uses compared: (Rd1) Douglas fir forest plantation; (RdA1) Arable cropland; (Rd2) Mixed P. nigra and P. sylvestris forest plantation; (RdA2) Abandoned land.
Tab. 4 - Carbon and nitrogen contents C/N ratios and C stock in forest floor/dead litter ± standard deviation. Different letters indicate significant differences among the means after Tukey’s test (p ≤ 0.05).

C/N ratio

Soil C/N ratios ranged from 3.4 (RdA2) to 15.2 (Rd2), detected at 20-40 cm depth. The difference between the average values along the soil profile at the four sites was statistically significant (F = 9.951, P < 0.001). In the deepest layer, forest soils tended to have higher C/N ratios than cropland and abandoned land (Fig. 4), probably owing to a larger root-derived input.

Fig. 4 - The C/N ratio in soil layers at the study sites. Different letters indicate significant differences among means after Tukey’s test (P ≤ 0.05). Error bars represent the standard error of the mean values per layer for each land use. Land uses compared: (Rd1) Douglas fir forest plantation; (RdA1) Arable cropland; (Rd2) Mixed P. nigra and P. sylvestris forest plantation; (RdA2) Abandoned land.

SOC stock

Cumulated SOC stock differed significantly among sites (F = 44.069, P < 0.001) and soil depths (F = 45.168, P < 0.001). The SOC stock ranged from 8.9 to 35.6 t ha-1 at 0-10 cm depth. It was highest at Rd1 and lowest at RdA2 (Fig. 5). The SOC stock had a C content distribution along the profiles with the highest values at 0-10 cm depth and a general decline with depth (0-10 > 10-20 > 20-30 > 30-40 cm) for all sites but RdA1. The total SOC stock over the entire 0-40 cm soil profile differed by a factor of up to five, ranging from 21 to 103 t ha-1. Abandoned land RdA2 stored the least and forest lands the most carbon in mineral soil. The lowest SOC stock in RdA2 was caused by the low C content but not by a different soil bulk density (Tab. 3). Similarly, Li et al. ([40]) reported that conversion of natural forests to continuous cultivation and grazing decreased both the concentration and stock of SOC. Zhang et al. ([73]) found that SOC stock increased when arable cropland was abandoned. In our case, the lowest SOC stock was found under abandonment (RdA2). In comparison, SOC stock in topsoil was higher at Rd2. If we assume a similar starting condition prior to LUC, there are two possible explanations for this finding. Either there was a smaller SOC accumulation under abandonment, after cropland conversion, as compared with afforestation, or the abandoned site actually lost carbon since conversion. In any case, our results challenge the hypothesis that LUC from cropland to abandoned land might be beneficial for C accumulation in soil.

Fig. 5 - The SOC stock (t ha-1) in mineral soil layers at the study sites. Different letters indicate significant differences among means after Tukey’s test (P ≤ 0.05). Error bars represent the standard error of the mean values per layer for each land use. Land uses compared: (Rd1) Douglas fir forest plantation; (RdA1) Arable cropland; (Rd2) Mixed P. nigra and P. sylvestris forest plantation; (RdA2) Abandoned land.

Our data also indicate that SOC stock increased with elevation, which is consistent with the results of Wang et al. ([63]) and Zhang & McGrath ([71]). However, because of the different starting conditions prior to LUC, we cannot say whether this increase was related to climate and associated gradients in OM decomposition ([61], [39]).

The C content, SOC stock and clay content of all samples correlated highly significantly and positively with land-use type (Tab. 5). Soil C correlated significantly and positively with both N content (r = 0.815, P < 0.01) and clay content (r = 0.500, P < 0.05). The N content was significantly and positively correlated with pH (r = 0.538, P < 0.05) and C content (r = 0.815, P < 0.01) and negatively with sand content (r = -0.518, P < 0.05). These results are in good agreement with findings of Kilic et al. ([30]), who indicated significant correlations between some physical and chemical soil properties under different land uses in Turkey. Our results showed a mutual dependency of OM content and soil physico-chemical properties for soils of the Western Rhodopes and confirmed the role of OM as an ecosystem property ([53]). Our results further showed that total SOC stocks were greater in the forests compared with other sites, and similar findings were reported by Kammer et al. ([27]).

Tab. 5 - Pearson’s correlation matrix for various soil physico-chemical parameters. (**): p<0.01 (2-tailed); (*): p<0.05 (2-tailed).

Forest floor/dead litter

The forest floor formed in both forests was classified as “Moder” type ([8]). The pH was acidic with variations between 4.8 and 4.9 and, according to Tukey’s test, did not differ between both coniferous plantations (Tab. 2). The pH of litter materials from the other land-use types differed significantly among each other and from the forest lands (P < 0.05). The dead litter pH value was highest at RdA2 (abandoned land), where pH was near to neutral in the turf layer.

The dry weight of forest floor/dead litter differed significantly (P < 0.05) among the four sites. It was highest under abandonment (RdA2: 178.21 g m-2) and lowest under cropland use (RdA1: 14.92 g m-2). Sites RdA1 and RdA2 had completely different depth profiles, and the formed turf layer in the abandoned grassland likely played an important role in formation of plant inputs towards the mineral soil. On the other hand, the regular harvest export in the cropland resulted in very small amounts of litter.

The C content in forest floor/dead litter was not significantly different among the four sites (between 42.7 and 44.7% - Tab. 4). The N content varied significantly between land-use types (P < 0.05). It was lowest in the coniferous plantations (0.63-0.68%) and highest under cropland and abandoned land use (0.83 and 0.94%, respectively). The differences in C/N ratios between sites were statistically significant (P < 0.05, range 46.1-70.8) and indicated higher litter qualities, i.e., better decomposability for litter in non-forests than in forests. The C/N ratios were highest in forest floor litter (Rd1 > Rd2) and lowest in abandoned land litter (RdA2).

Forest floor/dead litter C stocks differed significantly between the four sites (F = 69.132, P < 0.001). Overall mean values for forest plantations were 5.78 and 7.01 t ha-1 for Rd2 and Rd1, respectively (Tab. 4). The cropland had the lowest dead litter C stock (0.34 t ha-1). A high dead litter accumulation on the abandoned land suggested a reduced microbial activity at this site or higher input of dead organic materials. Statistically significant and positive correlation was found for pH and N content (r = 0.748, P < 0.01), whereas the C/N ratio correlated significantly negatively with pH (r = -0.910, P < 0.01) and indicated that processes of nitrification of organic materials were related with pH.

Aboveground vegetation

The C stock was higher in the aboveground biomass of forests than of cropland or abandoned land (Fig. 6). Site Rd1 stored the most carbon in vegetation (196.7 t ha-1), followed by Rd2 (160.8 t ha-1). A substantial amount of carbon was accumulated also in aboveground parts of abandoned land, but more studies are required in order to reveal the seasonal variation and to estimate the annual mean. Concerning the aboveground tree biomass and carbon stored in it, more detailed studies are required to obtain data for branches and needles. The latter two are not explicitly considered in the growing tables applied in national forestry, especially for Douglas fir stands.

Fig. 6 - Carbon stock distribution in ecosystem compartments of studied land-use types. Different letters indicate significant differences among means after Tukey’s test (P ≤ 0.05). Error bars represent the standard error of the mean per layer for each land use. Land uses compared: (Rd1) Douglas fir forest plantation; (RdA1) Arable cropland; (Rd2) Mixed P. nigra and P. sylvestris forest plantation; (RdA2) Abandoned land. Compartments included are: Soil, FF/litter - forest floor (F+H) and litter layer, Aboveground - tree and grass vegetation.

Carbon stock distribution in ecosystem compartments

The C stock distribution in different ecosystem compartments showed that in forest land and under abandonment, aboveground vegetation represented the largest C pool. The amount of carbon in aboveground vegetation was similar for both forest stands, although they consisted of different coniferous species and with different stand characteristics (Fig. 6). The forest floor was another important C pool at the forest sites. At the abandoned land site, both dead litter and standing biomass contributed the most to the overall ecosystem C storage whereas mineral soil C played a minor role. At the arable cropland site, carbon was mostly stored in the mineral soil whereas fresh litter was removed from the system via harvest, leaving little litter material in the ecosystem.


Land-use change in the Western Rhodope Mountains is one of the major drivers for modulating soil characteristics. We found significant changes and variations of soil C storage down to 40 cm of soil depth, although differences were most pronounced in the superficial 0-10 cm of organo-mineral soil. Afforestation with coniferous tree species caused changes in soil pH expressed as acidification of the superficial soil layers. Typical mountainous agricultural land use was related with changes in textural fractions suggesting ongoing processes of superficial erosion. Different soil N contents were most significant between Rd1 and RdA1 as a result of continuous N removal from arable cropland (RdA1) that was not fully compensated by fertilizer additions. It could be summarized that due to afforestation activities, the N content in upper soil layers increased.

Afforestation helped accumulate and preserve soil C, but even several decades after the conversion from arable cropland or abandoned land to forest land use, this effect was expressed only weakly and could not compensate the organic C loss due to previous active cultivation. The C content in soil and its profile distribution were shown to be very sensitive to land-use and coverage changes and could be proposed as appropriate indicators for studying the long-term effects of LUC on soil properties. Of all studied compartments of mountain ecosystems, mineral soil remained the main accumulator of SOC for all studied land-use types. Together with mineral soil, aboveground tree vegetation also accumulated a major part of carbon at the studied forest and abandoned sites.

Afforestation considerably increased organic C stock in the uppermost soil layers of different land-use types in the Western Rhodope Mountains. Afforestation could be identified as a good practice that could preserve soil C and improve the regulating and sustaining ecosystem services of mountain regions. Therefore, we recommend conserving mountain soils by preventing abandonment and by establishing appropriate forest plantations.


The work was a part of the project “Land-use and management impacts on carbon sequestration in mountain ecosystems” financed within the framework of the Bulgarian-Swiss Research Programme (BSRP), 2011-2016, between the Government of Bulgaria and the Swiss Federal Council, contract number IZEBZ0_142972.


Alexandrov AH (1994). The anthropogenic interference at restoring the upper forest limit. Forest Science 1: 3-7. [in Bulgarian]
::Google Scholar::
Angers DA, Caron J (1998). Plant-induced changes in soil structure: process and feedbacks. Biogeochemistry 42: 55-72.
::CrossRef::Google Scholar::
Arinushkina EV (1970). Manual in soil chemical analysis. Publishing House at Moscow State University, Moscow, pp. 487. [in Russian]
::Google Scholar::
Ayoubi S, Khormali F, Sahrawat KL, Rodrigues de Lima AC (2011). Assessing impacts of land use change on soil quality indicators in a loessial soil in Golestan Province, Iran. Journal of Agricultural Science and Technology 13: 727-742.
::Online::Google Scholar::
Beheshti A, Raiesi F, Golchin A (2012). Soil properties, C fractions and their dynamics in land use conversion from native forests to croplands in northern Iran. Agriculture, Ecosystems and Environment 148: 121-133.
::CrossRef::Google Scholar::
Binkley D (1995). The influence of tree species on forest soils: processes and patterns. In: Proceedings of the “Trees and Soil Workshop” (Mead DJ, Cornforth IS eds). Lincoln University (New Zealand) 28 Feb - 2 Mar 1994. Agronomy Society of New Zealand, Special Publication No. 10, Lincoln University Press, Canterbury, New Zealand, pp. 1-33.
::Google Scholar::
Bradstreet RB (1954). Kjeldahl method for organic nitrogen. Analytical Chemistry 26 (1): 185-187.
::CrossRef::Google Scholar::
Brêthes A, Brun JJ, Jabiol B, Ponge JF, Toutain F (1998). Types of humus forms in temperate forests. In: “A sound reference base for soils” (Baize D ed). AFES, INRA, Versailles, France, pp. 265-282.
::Google Scholar::
Caravaca F, Lax A, Albaladejo J (2004). Aggregate stability and carbon characteristics of particle-size fractions in cultivated and forested soils of semiarid Spain. Soil and Tillage Research 78: 83-90.
::CrossRef::Google Scholar::
Celik I (2005). Land-use effects on organic matter and physical properties of soil in a southern Mediterranean highland of Turkey. Soil and Tillage Research 83: 270-277.
::CrossRef::Google Scholar::
Donov V, Gentcheva S, Yorova K (1974). Guidance for exercises in forest soil science. Zemizdat, Sofia, Bulgaria, pp. 220. [on Bulgarian]
::Google Scholar::
EEA (2003). Environmental assessment report no. 10, Europe’s environment: the third assessment. European Environmental Agency, Copenhagen, Denmark, pp. 344.
::Online::Google Scholar::
EEA (2013). Annual report 2012 and environmental statement 2013. European Environmental Agency, Copenhagen, Denmark, pp. 100.
::Online::Google Scholar::
Ellert BH, Bettany JR (1995). Calculation of organic matter and nutrients stored in soils under contrasting management regimes. Canadian Journal of Soil Science 75: 529-538.
::CrossRef::Google Scholar::
ExAF (2010). Annual report of Executive Agency of Forestry for 2009. Ministry of Agriculture and Food, Sofia, Bulgaria, pp. 71. [in Bulgarian]
::Online::Google Scholar::
Foster GR (2001). Keynote: soil erosion prediction technology for conservation planning. In: Proceedings of the 10th ISCO Meeting “Sustaining the global farm” (Stott DE, Mohtar RH, Steinhardt GC eds). Purdue University (West Lafayette, IN, USA) 24-29 May 1999. International Soil Conservation Organization, Tuscon, AZ, USA, pp. 847-851.
::Online::Google Scholar::
Foth HD, Ellis BG (1997). Soil fertility (2nd edn). Lewis CRC Press LLC, USA, pp. 290.
::Online::Google Scholar::
Guo LB, Gifford RM (2002). Soil carbon stocks and land use change: a meta analysis. Global Change Biology 8: 345-360.
::CrossRef::Google Scholar::
Hebel A (1998). Soil degradation-diagnosis, appraisal and reversing measure. Introduction. In: “Toward sustainable land use” (Blume HP, Eger H, Fleischhauer E, Hebel A, Reij C, Steiner KG eds). Advances in Geoecology 31, Vol. I, Catena Verlag, Reiskirchen, Germany, pp. 1-2.
::Google Scholar::
Holubík O, Podrázský V, Vopravil J, Khel T, Remeš J (2014). Effect of agricultural lands afforestation and tree species composition on the soil reaction, total organic carbon and nitrogen content in the uppermost mineral soil profile. Soil and Water Resources 9: 192-200.
::Online::Google Scholar::
Houghton RA (1999). The annual net flux of carbon to the atmosphere from changes in land use 1850-1990. Tellus B 51: 298-313.
::CrossRef::Google Scholar::
Hulugalle NR, Lobry De Bruyn LA, Entwistle P (1997). Residual effects of tillage and crop rotation on soil properties, soil invertebrate numbers and nutrient uptake in an irrigated Vertisol sown to cotton. Applied Soil Ecology 7: 11-30.
::CrossRef::Google Scholar::
IPCC (2003). Good practice guidance for land use, land-use change and forestry. Intergovernmental Panel on Climate Change, Institute for Global Environmental Strategies, Kanagawa, Japan, pp. 590.
::Online::Google Scholar::
ISO-10390 (2005). Soil quality: determination of pH. International Organization of Standardization, Geneva, Switzerland.
::Online::Google Scholar::
IUSS WG WRB (2014). World reference base for soil resources 2014. International soil classification system for naming soils and creating legends for soil maps. World Soil Resources Reports No. 106, FAO, Rome, Italy, pp. 191.
::Google Scholar::
Johnson JMF, Allmaras RR, Reicosky DC (2006). Estimating source carbon from crop residues, roots and rhizodeposits using the national grain-yield database. Agronomy Journal 98: 622-636.
::CrossRef::Google Scholar::
Kammer A, Hagedorn F, Shevchenko I, Leifeld J, Guggenberger G, Goryacheva T, Rigling A, Moiseev P (2009). Treeline shifts in the Ural Mountains affect soil organic matter dynamics. Global Change Biology 15: 1570-1583.
::CrossRef::Google Scholar::
Karlen DL, Mausbach MJ, Doran JW, Cline RT, Harris RF, Schuman GE (1997). Soil quality: a concept definition and framework for evaluation. Soil Science Society of America Journal 90: 644-650.
::CrossRef::Google Scholar::
Khormali F, Nabiollahi K (2009). Degradation of Mollisols in Western Iran as affected by land use change. Journal of Agricultural Science and Technology 11: 363-374.
::Online::Google Scholar::
Kilic K, Kilic S, Kocyigit R (2012). Assessment of spatial variability of soil properties in areas under different land use. Bulgarian Journal of Agricultural Science 18: 722-732.
::Online::Google Scholar::
Koulouri M, Giourga C (2007). Land abandonment and slope gradient as key factors of soil erosion in Mediterranean terraced lands. Catena 69: 274-281.
::CrossRef::Google Scholar::
Körner C, Paulsen J (2004). A world-wide study of high altitude tree-line temperatures. Journal of Biogeography 31: 713-732.
::CrossRef::Google Scholar::
Krastanov K, Beliakov P, Shikov K (2004). Table for form heights of Scots pine plantation with bark. In: “Reference book in dendrobiometry” (Krastanov K, Raikov R, eds.). Bulprophor Publishing House, Sofia, Bulgaria, pp. 611. [in Bulgarian]
::Google Scholar::
Krastanov K, Raikov R (2004). Handbook of dendrobiometry. Bulprophor Publishing House, Sofia, Bulgaria, pp. 612. [in Bulgarian]
::Google Scholar::
Ladha JK, Reddy CK, Padre AT, Van Kessel C (2011). Role of nitrogen fertilization in sustaining organic matter in cultivated soils. Journal of Environmental Quality 40: 1756-1766.
::CrossRef::Google Scholar::
Laganière J, Angers DA, Pare D (2010). Carbon accumulation in agricultural soils after afforestation: a meta-analysis. Global Change Biology 6: 439-453.
::CrossRef::Google Scholar::
Laganière J, Paré D, Bergeron Y, Chen HH, Brassard B, Cavard X (2013). Stability of soil carbon stocks varies with forest composition in the Canadian boreal biome. Ecosystems 16: 852-865.
::CrossRef::Google Scholar::
Lal R (1986). Conversion of tropical rainforest: agronomic potential and ecological consequences. Advances in Agronomy 39: 173-263.
::CrossRef::Google Scholar::
Leifeld J, Bassin S, Fuhrer J (2005). Carbon stocks in Swiss agricultural soils predicted by land-use, soil characteristics and altitude. Agriculture, Ecosystems and Environment 105: 255-266.
::CrossRef::Google Scholar::
Li J, Edwards PC, Burghammer M, Villa C (2004). The effects of land use change on soil properties and organic carbon at Dagdam River catchment in Turkey. Journal of Environmental Biology 30 (5): 825-830.
::Google Scholar::
Lima AMN, Silva IR, Neves JCL, Novais RF, Barros NF, Mendonca ES, Smyth TJ, Moreira MS, Leite FP (2006). Soil organic carbon dynamics following afforestation of degraded pastures with eucalyptus in southeastern Brazil. Forest Ecology and Management 235: 219-223.
::CrossRef::Google Scholar::
MAES (2014). Mapping and assessment of ecosystems and their services. Technical Report 2014-080, European Commission, Directorate-General for the Environment, Brussels, Belgium, pp. 82.
::CrossRef::Google Scholar::
Masarov S (1940). Juniperus as a pioneer species of coniferous stands over pastures. Forestry Ideas 2: 123-136. [in Bulgarian]
::Google Scholar::
Matano AS, Kanangire CK, Anyona DN, Abuom PO, Gelder FB, Dida GO, Owuor PO, Ofulla AVO (2015). Effects of land use change on land degradation reflected by soil properties along Mara River, Kenya and Tanzania. Open Journal of Soil Science 5: 20-38.
::CrossRef::Google Scholar::
Mooney HA, Vitousek PV, Matson PA (1987). Exchange of materials between terrestrial ecosystems and the atmosphere. Science 238: 926-932.
::CrossRef::Google Scholar::
Nedyalkov S, Krastanov K, Tsakov H (2004). Table for form heights of black pine plantations with bark. In: “Reference book in dendrobiometry” (Krastanov K, Raikov R eds). Bulprophor Publishing, Sofia, Bulgaria, pp. 612. [in Bulgarian]
::Google Scholar::
Poeplau C, Don A (2013). Sensitivity of soil organic carbon stocks and fractions to different land-use changes across Europe. Geoderma 192: 189-201.
::CrossRef::Google Scholar::
Post WM, Kwon KC (2000). Soil carbon sequestration and land-use change: processes and potential. Global Change Biology 6: 317-327.
::CrossRef::Google Scholar::
Potter CS, Matson PA, Vitousek PM, Davidson E (1996). Process modeling of controls on nitrogen trace gas emissions from soils worldwide. Journal of Geophysical Research 101 (D1): 1361-1377.
::CrossRef::Google Scholar::
Powlson DS (2005). Will soil amplify climate change? Nature 433: 204-205.
::CrossRef::Google Scholar::
Preston CM, Schmidt MWI (2006). Black (pyrogenic) carbon: a synthesis of current knowledge and uncertainties with special consideration of boreal regions. Biogeosciences 3: 397-420.
::CrossRef::Google Scholar::
Sabev A, Stanev S (1963). Climatic districts in Bulgaria and their climate. Zemizdat Edition, Sofia, Bulgaria, pp. 180-185. [in Bulgarian]
::Google Scholar::
Schmidt MWI, Torn MS, Abiven S, Dittmar T, Guggenberger G, Janssens IA, Kleber M, Kögel-Knabner I, Lehmann J, Manning DAC, Nannipieri P, Rasse DP, Weiner S, Trumbore SE (2011). Persistence of soil organic matter as an ecosystem property. Nature 478: 49-56.
::CrossRef::Google Scholar::
Schulze ED, Luyssaert S, Ciais P, Freibauer A, Janssens IA, Soussana JF, Smith P, Grace J, Levin I, Thiruchittampalam B, Heimann M, Dolman AJ, Valentini R, Bousquet P, Peylin P, Peters W, Rodenbeck C, Etiope G, Vuichard N, Wattenbach M, Nabuurs GJ, Poussi Z, Nieschulze J, Gash JH (2009). Importance of methane and nitrous oxide for Europe’s terrestrial greenhouse-gas balance. Nature Geoscience 2: 842-850.
::CrossRef::Google Scholar::
Shikov K (1972). Growth and productivity of the Scots and black pine under the same conditions of habitat. Gorsko stopanstvo 8: 18-22. [In Bulgarian]
::Google Scholar::
Skousen J, Sencindiver J, Owens K, Hoover S. (1998). Physical properties of minesoils in West Virginia and their influence on wastewater treatment. Journal of Environmental Quality 27: 633-639.
::CrossRef::Google Scholar::
Spehn EM, Liberman M, Körner C (2006). Fire and grazing’a synthesis of human impacts on highland biodiversity. In: “Land use change and mountain biodiversity” (Spehn EM, Liberman M, Körner C eds). CRC Publishers, Boca Raton, FL, USA, pp. 337-347.
::Google Scholar::
Swan ARH, Sandilands M (1995). Introduction to geological data analysis. Blackwell, London, UK, pp. 446.
::Online::Google Scholar::
Tasser E, Tappeiner U (2009). Impact of land use changes on mountain vegetation. Applied Vegetation Science 5: 173-184.
::CrossRef::Google Scholar::
Tisdall JM, Oades JM (1982). Organic matter and water stable aggregates in soils. Journal of Soil Science 33: 141-163.
::CrossRef::Google Scholar::
Trumbore SE, Chadwick OA, Amundson R (1996). Rapid exchange between soil carbon and atmospheric carbon dioxide driven by temperature change. Science 272: 393-396.
::CrossRef::Google Scholar::
Vagen TG, Andrianorofanomezana MAA, Andrianorofanomezana S (2006). Deforestation and cultivation effects on characteristics of Oxisols in the Highlands of Madagascar. Geoderma 131: 190-200.
::CrossRef::Google Scholar::
Wang HQ, Hall CAS, Cornell JD, Hall MHP (2002). Spatial dependence and the relationship of soil organic carbon and soil moisture in the Luquillo Experimental Forest, Puerto Rico. Landscape Ecology 17: 671-684.
::CrossRef::Google Scholar::
Ward A, Dargusch P, Thomas S, Liu Y, Fulton EA (2014). A global estimate of carbon stored in the world’s mountain grasslands and shrublands, and the implications for climate policy. Global Environmental Change 28: 14-24.
::CrossRef::Google Scholar::
Wendt JW, Hauser S (2013). An equivalent soil mass procedure for monitoring soil organic carbon in multiple soil layers. European Journal of Soil Science 64: 58-65.
::CrossRef::Google Scholar::
Wichern F, Eberhardt E, Mayer J, Joergensen RG, Müller T (2008). Nitrogen rhizodeposition in agricultural crops: methods, estimates and future prospects. Soil Biology and Biochemistry 40: 30-48.
::CrossRef::Google Scholar::
Wiesmeier M, Spörlein P, Geuß U, Hangen E, Haug S, Reischl A, Schilling B, Von Lützow M, Kögel-Knabner I (2012). Soil organic carbon stocks in southeast Germany (Bavaria) as affected by land use: soil type and sampling depth. Global Change Biology 18: 2233-2245.
::CrossRef::Google Scholar::
Wiesmeier M, Hübner R, Barthold F, Spörlein P, Geuß U, Hangen E, Reischl A, Schilling B, Von Lützow M, Kögel-Knabner I (2013). Amount, distribution and driving factors of soil organic carbon and nitrogen in cropland and grassland soils of southeast Germany (Bavaria). Agriculture, Ecosystems and Environment 176: 39-52.
::CrossRef::Google Scholar::
Yang JC, Huang JH, Pan QM, Tang JW, Han XG (2004). Long-term impacts of land-use change on dynamics of tropical soil carbon and nitrogen pools. Journal of Environmental Science (China) 16: 256-261.
::Google Scholar::
Zenga X, Zhang W, Cao J, Liua X, Shen H, Zhao X (2014). Changes in soil organic carbon, nitrogen, phosphorus, and bulk density after afforestation of the “Beijing-Tianjin Sandstorm Source Control” program in China. Catena 118: 186-194.
::CrossRef::Google Scholar::
Zhang CS, McGrath D (2003). Geostatistical and GIS analyses on soil organic carbon concentrations in grassland of southeastern Ireland from two different periods. Geoderma 119: 261-275.
::CrossRef::Google Scholar::
Zhang GL (2010). Changes of soil labile organic carbon in different land uses in Sanjiang Plain, Heilongjiang Province. Chinese Geographical Science 20: 139-143.
::CrossRef::Google Scholar::
Zhang LH, Xie ZK, Zhao RF, Wang YJ (2012). The impact of land use change on soil organic carbon and labile organic carbon stocks in the Longzhong Region of Loess Plateau. Journal of Arid Land 4: 241-250.
::CrossRef::Google Scholar::
Zhao P, Shao MA, Omran W, Amer AM (2011). Effects of erosion and deposition on particle size distribution of deposited farmland soils on the Chinese Loess Plateau. Revista Brasileira de Ciência do Solo 35: 2135-2144.
::CrossRef::Google Scholar::
Zink JA, Farshad A (1995). Issues of sustainability and sustainable land management. Canadian Journal of Soil Science 75: 407-412.
::CrossRef::Google Scholar::


Paper Contents

Paper Sections

Paper Figures

Paper Tables



Zhiyanski M, Glushkova M, Ferezliev A, Menichetti L, Leifeld J (2016).
Carbon storage and soil property changes following afforestation in mountain ecosystems of the Western Rhodopes, Bulgaria
iForest - Biogeosciences and Forestry 9: 626-634. - doi: 10.3832/ifor1866-008
First Previous Next Last
© iForest

Download Reference

Paper ID# ifor1866-008
Title Carbon storage and soil property changes following afforestation in mountain ecosystems of the Western Rhodopes, Bulgaria
Authors Zhiyanski M, Glushkova M, Ferezliev A, Menichetti L, Leifeld J
Close Download