# The Effect of Land Cover Change on Soil Properties around Kibale National Park in South Western Uganda

Authors: J. G. M. Majaliwa; R. Twongyirwe; R. Nyenje; M. Oluka; B. Ongom; J. Sirike; D. Mfitumukiza; E. Azanga; R. Natumanya; R. Mwerera; B. Barasa
Journal: Applied and Environmental Soil Science
 (2010)

---

## Abstract

The change from natural forest cover to tea and Eucalyptus is rampant in protected areas of western Uganda. The objectives were; to examine the trend in land-use /cover change and determine the effect of these changes on the physico-chemical properties of soils around Kibale National Park. The trend in land use/cover change was assessed by analyzing a series of landsat images. Focused group discussions and key informant interviews were used for land-use/cover reconstruction.  Three major land uses were included; woodlot  (Eucalyptus grandis; 5 years old) ), tea (57 years old) and natural forest used as a control. Each of these land-uses were selected at two different North facing landscape positions and were replicated three times.  A total of 36 composite soil samples were taken at 0–15 and 15–30 cm depth from natural forest, Tea plantation and eucalyptus on three ridges. Results showed that small scale farming, tea and eucalyptus plantation and built up area have increased over time, to the expense of woodlot and forest cover.  Tea and Eucalyptus have induced changes in: exchangeable Mg and Ca, available  P, SOM, pH, and bulk density of sub soil  (P<.05). Landscape positions within land use also significantly influenced most soil properties (P<.05). Similar findings were observed by Wang et al. (2006) in commercial tea plantations in China that received  nitrogen fertilizers.


---

## Body

## 1. Introduction

Forest degradation in Uganda continues to be an issue of serious concern threatening to accelerate land degradation [1]. This forest loss has been mostly attributed to anthropogenic activities other than climatic forces [2]. Areas affected by the degradation include forests within and around Kibale National Park.Kibale National Park is located in Western Uganda near the foothills of the Rwenzori Mountains and contains moist evergreen forest [3]. This natural forest is punctuated with patches of grassland which are mainly around the summits of the hill tops. Earlier, people within and around the forests coexisted sustainably and had access to the forest for basic needs such as water, firewood, vegetables, honey, and medicine for domestic use [4]. However, in the 1970s and 1980s the forest experienced serious degradation mainly due to anthropogenic activities. Despite the forest providing important social and environmental benefits, it is increasingly threatened by accelerating rates of forest conversion and degradation [5]. As the area experiences crop destruction due to wild animal raids especially elephants, people have resorted to tree plantation establishment notably, eucalyptus and pines [6]. Gradually these plantations have expanded horizontally forming new major vegetation covers.Land use practices affect the distribution and supply of soil nutrients by directly altering soil properties and by influencing biological transformations in the rooting zone [7]. For instance, cultivation of forests diminishes the soil carbon (C) within a few years of initial conversion [7] and substantially lowers mineralisable of Nitrogen (N) [8]. The conversion of forest to crop land has been associated with reduction in organic matter content of the top soil [9, 10] and subsequently, decline in productivity, since organic matter content is responsible for the productivity in soils [11, 12]. Islam and Weil [13] reported an increase in bulk density and a reduction in porosity and aggregate stability following the conversion of forest land to crop land. Similar findings were reported by Motavalli and McConnell [14] and Riezebos and Loerts [15]. Such changes in soil properties predispose the soil to soil erosion, a major soil degradation process in Uganda [16]. This study was therefore carried out to determine the trend in land use/cover change and assess how these trends have impacted the soil properties in and around Kibale National Park.


## 2. Materials and Methods

### 2.1. Description of the Study Area

The study was undertaken in Burahya county in Kabarole district (Uganda). The area is located between latitude0° 13′–0° 41′ N and longitude 30° 19′–30° 32′ E in Western Uganda near the foothills of the Rwenzori Mountains and covers approximately 10.16 km2. It lies between 1390 m and 1625 m above sea level. Currently, moist evergreen forest, tea plantation, small-scale agriculture, and woodland plantations are the dominant land use/cover in the area [3, 17]. The national park consists of mature, mid-altitude, moist semideciduous and evergreen forest (57%), grassland (15%), woodland (4%), lakes and wetlands (2%), colonizing forest (19%), and plantations of exotic trees (1%) primarily Cupressus lusitanica, Pinus patula, P.  caribaea, and Eucalyptus species, [3]. The soil is classified as Lixic Ferralsols.  The area receives bimodal rainfall occurring from March to May and from September to November. The mean annual rainfall in the region is 1750 mm whereas the mean daily minimum and maximum temperatures are 15.5°C and 23.7°C, respectively.


### 2.2. Research Approach

#### 2.2.1. Trends in Land Use/Cover Change

Land use/cover change trends in the area were determined by analyzing series of satellite images (landsat images) of 1974 and 2000 covering the study area. In addition, the land use/cover map of 1986 was obtained from the National Biomass Study [18]. In Addition, a land use map of the two subcounties was clipped from the USGS (2005) land use/cover map of Africa. The major land use /cover classes included natural forest, grassland, small-scale cultivation and large-scale agriculture and woodland. Large scale agriculture was mainly Tea plantations; while woodland mainly consisted of eucalyptus in the study area. Preliminary maps were produced using Integrated Land and Water Information System (ILWIS version 3.1). The 1974 land use/cover of the area was reconstructed using targeted interviews. The preliminary map produced from the analysed landsat image of 2000 and field observations were used to update the USGS (2005) map of the study area and generate a land use/cover map of the year 2009.


#### 2.2.2. Effect of Land Use Change on Soil Properties

For selection of the sampling sites, two transects of 5 km cutting across the dominant land use/cover were drawn in Arc View 3.1. A reconnaissance was carried out in the field and three major land uses were selected based on their history, occurrence at different landscape positions, and possibility of having a minimum of three replications. The selected sampling sites were those originally under forest and have been converted to the current land use/cover. Two land uses satisfied the set criteria namely, woodlot (Eucalyptus grandis;5 years old) and tea (57 years old). The natural forest was used as a control. Each of these land uses was selected at two different North facing landscape positions (Back slope, foot slope) [19] and were replicated three times.For each replicate, a composite sample was obtained from three different points at two depths 0–15 cm and 15–30 cm. The soil depth of 0–30 cm was chosen because it is the agricultural layer [20]. Soil core samples were also taken one from each depth (0–15 and 15–30 cm) for determination of bulk density and hydraulic conductivity using the core method and falling head method, respectively [21]. Soil samples were air-dried and analyzed for pH, organic matter (SOM), total N, available P, exchangeable bases (Ca, Mg, and K) following standard methods [22]. Soil pH was measured using pH a meter (1 : 2.5 soil: water); SOM and total N were determined using Walkley and Black method and Kjedhal method, respectively. Soil texture was determined using Bouyoucos hydrometer method [23] and their textural classes according to the FAO classification [24].


### 2.3. Statistical Analysis

The trend in each of the land use/cover was estimated using regression techniques. Separation of means was done in Genstat Discovery version using the LSD. Values of soil properties that differed atP≤  .05 were considered significant trends.


## 2.1. Description of the Study Area

The study was undertaken in Burahya county in Kabarole district (Uganda). The area is located between latitude0° 13′–0° 41′ N and longitude 30° 19′–30° 32′ E in Western Uganda near the foothills of the Rwenzori Mountains and covers approximately 10.16 km2. It lies between 1390 m and 1625 m above sea level. Currently, moist evergreen forest, tea plantation, small-scale agriculture, and woodland plantations are the dominant land use/cover in the area [3, 17]. The national park consists of mature, mid-altitude, moist semideciduous and evergreen forest (57%), grassland (15%), woodland (4%), lakes and wetlands (2%), colonizing forest (19%), and plantations of exotic trees (1%) primarily Cupressus lusitanica, Pinus patula, P.  caribaea, and Eucalyptus species, [3]. The soil is classified as Lixic Ferralsols.  The area receives bimodal rainfall occurring from March to May and from September to November. The mean annual rainfall in the region is 1750 mm whereas the mean daily minimum and maximum temperatures are 15.5°C and 23.7°C, respectively.


## 2.2. Research Approach

### 2.2.1. Trends in Land Use/Cover Change

Land use/cover change trends in the area were determined by analyzing series of satellite images (landsat images) of 1974 and 2000 covering the study area. In addition, the land use/cover map of 1986 was obtained from the National Biomass Study [18]. In Addition, a land use map of the two subcounties was clipped from the USGS (2005) land use/cover map of Africa. The major land use /cover classes included natural forest, grassland, small-scale cultivation and large-scale agriculture and woodland. Large scale agriculture was mainly Tea plantations; while woodland mainly consisted of eucalyptus in the study area. Preliminary maps were produced using Integrated Land and Water Information System (ILWIS version 3.1). The 1974 land use/cover of the area was reconstructed using targeted interviews. The preliminary map produced from the analysed landsat image of 2000 and field observations were used to update the USGS (2005) map of the study area and generate a land use/cover map of the year 2009.


### 2.2.2. Effect of Land Use Change on Soil Properties

For selection of the sampling sites, two transects of 5 km cutting across the dominant land use/cover were drawn in Arc View 3.1. A reconnaissance was carried out in the field and three major land uses were selected based on their history, occurrence at different landscape positions, and possibility of having a minimum of three replications. The selected sampling sites were those originally under forest and have been converted to the current land use/cover. Two land uses satisfied the set criteria namely, woodlot (Eucalyptus grandis;5 years old) and tea (57 years old). The natural forest was used as a control. Each of these land uses was selected at two different North facing landscape positions (Back slope, foot slope) [19] and were replicated three times.For each replicate, a composite sample was obtained from three different points at two depths 0–15 cm and 15–30 cm. The soil depth of 0–30 cm was chosen because it is the agricultural layer [20]. Soil core samples were also taken one from each depth (0–15 and 15–30 cm) for determination of bulk density and hydraulic conductivity using the core method and falling head method, respectively [21]. Soil samples were air-dried and analyzed for pH, organic matter (SOM), total N, available P, exchangeable bases (Ca, Mg, and K) following standard methods [22]. Soil pH was measured using pH a meter (1 : 2.5 soil: water); SOM and total N were determined using Walkley and Black method and Kjedhal method, respectively. Soil texture was determined using Bouyoucos hydrometer method [23] and their textural classes according to the FAO classification [24].


## 2.2.1. Trends in Land Use/Cover Change

Land use/cover change trends in the area were determined by analyzing series of satellite images (landsat images) of 1974 and 2000 covering the study area. In addition, the land use/cover map of 1986 was obtained from the National Biomass Study [18]. In Addition, a land use map of the two subcounties was clipped from the USGS (2005) land use/cover map of Africa. The major land use /cover classes included natural forest, grassland, small-scale cultivation and large-scale agriculture and woodland. Large scale agriculture was mainly Tea plantations; while woodland mainly consisted of eucalyptus in the study area. Preliminary maps were produced using Integrated Land and Water Information System (ILWIS version 3.1). The 1974 land use/cover of the area was reconstructed using targeted interviews. The preliminary map produced from the analysed landsat image of 2000 and field observations were used to update the USGS (2005) map of the study area and generate a land use/cover map of the year 2009.


## 2.2.2. Effect of Land Use Change on Soil Properties

For selection of the sampling sites, two transects of 5 km cutting across the dominant land use/cover were drawn in Arc View 3.1. A reconnaissance was carried out in the field and three major land uses were selected based on their history, occurrence at different landscape positions, and possibility of having a minimum of three replications. The selected sampling sites were those originally under forest and have been converted to the current land use/cover. Two land uses satisfied the set criteria namely, woodlot (Eucalyptus grandis;5 years old) and tea (57 years old). The natural forest was used as a control. Each of these land uses was selected at two different North facing landscape positions (Back slope, foot slope) [19] and were replicated three times.For each replicate, a composite sample was obtained from three different points at two depths 0–15 cm and 15–30 cm. The soil depth of 0–30 cm was chosen because it is the agricultural layer [20]. Soil core samples were also taken one from each depth (0–15 and 15–30 cm) for determination of bulk density and hydraulic conductivity using the core method and falling head method, respectively [21]. Soil samples were air-dried and analyzed for pH, organic matter (SOM), total N, available P, exchangeable bases (Ca, Mg, and K) following standard methods [22]. Soil pH was measured using pH a meter (1 : 2.5 soil: water); SOM and total N were determined using Walkley and Black method and Kjedhal method, respectively. Soil texture was determined using Bouyoucos hydrometer method [23] and their textural classes according to the FAO classification [24].


## 2.3. Statistical Analysis

The trend in each of the land use/cover was estimated using regression techniques. Separation of means was done in Genstat Discovery version using the LSD. Values of soil properties that differed atP≤  .05 were considered significant trends.


## 3. Results

### 3.1. Land Use/Cover Change Maps in 1973, 1986, and 2009

Figure2 and Table 1 show the land use/cover for three time series. Generally small-scale farming, tea plantation, and built up area have increased over time, and woodlot has decreased with time. Previous studies have shown this continuing trend [25].Table 1
showing land use and cover changes.
1973%1987%2009%Built up areas90.5181281.5Grasslands151493361.9Small-scale farming76041125969130871.4Tea plantations1514831789.7Tropical forest828453852125413.8Woodlands20711744291.6Figure 1
Bunyanganbu and Burahya sub-counties in Fort Portal district in Uganda.land use/cover for three time series, 1973, 1987, 2009.
(a)(b)(c)From the figures above, the highest percentage of land in 1973 was tropical forest (45%) compared to small-scale farming in 1987 (69%) and 2009 (71.4%). The smallest percentage of land in 1973 (0.5%) was built up areas compared to tea plantations and grasslands in 1987 (3%) and 2009 (1.9%), respectively. The changes in land use could be attributed to agricultural expansion that took place in the areas adjacent to Kibale National Park [25] as well as an increase in energy demand (firewood). Such expansion was associated with arrival of new immigrants into the areas.


#### 3.1.1. Effect of Land Use on Selected Soil Properties

Soil physical and chemical properties in natural forests, tea plantations, and Eucalyptus plantations are presented in Tables2 and 3 for soil depths 0–15 cm and 15–30 cm, respectively. Top soil (0–15 cm) properties except total N, exchangeable potassium, and sodium differed significantly (P<  .05) across the land uses. Textural analysis indicates that the fraction of clay highly significantly (P<  .001) differed between LUs. Clay contents were higher in Eucalyptus than in tea plantations and natural forests. Bulk density was lower in the forest compared to Eucalyptus and Tea. Soil pH (mean = 3.53 ± 0.41) in tea plantations was statistically lower than in natural forests and in Eucalyptus plantations, respectively. Similarly, mean values of SOM were statistically higher in the respective LUs. Total N in natural forests and Eucalyptus were roughly equal.  Available phosphorus in natural forests was about 6 and 3 times lower than in tea plantations and in Eucalyptus plantations, respectively. The highest level of exchangeable Mg was observed in natural forest followed by Eucalyptus plantations and tea plantations, respectively.Table 2
Mean± SD of selected soil (0–15 cm) physical and chemical properties in the three land uses.
Soil propertiesLand use/coverP<  .05Natural forestTeaEucalyptuspH4.72±0.55a3.53±0.41b5.13±0.61c<  .001SOM (%)6.48±0.56a7.20±1.15b4.97±0.45c.002Total N (%)0.26± 0.070.31± 0.050.24± 0.04NsAv. P (ppm)3.32±0.99a20.50±9.52b11.30±10.03a.004Ca (cmol(p+)/kg)9.70±4.18a2.30±1.90b9.30±5.20a.010Mg (cmol(p+)/kg)2.26±0.80a0.66±0.74b2.00±0.55a<  .001K (cmol(p+)/kg)0.21± 0.150.29± 0.190.34± 0.20NsNa (cmol(p+)/Kg)0.04± 0.010.05± 0.010.04± 0.01NsBulk density (g/cm3)0.72±0.05a1.18±0.08b1.22±0.08b.018Ksat (cm/hr)21.90± 24.3519.03± 14.7814.86± 17.84NsClay (%)15.67±4.46a23.33±7.55b36.00±3.35c<  .001Silt (%)13.67± 2.9413.00± 3.0316.00± 2.53NsTextural classSandy LoamSandy clay LoamSandy ClayValues followed by the same letters in a row are not significantly different atP<  .05; n.s.: not significant.Table 3
Mean± SD of selected soil (15–30 cm) physical and chemical properties in the three land use/covers.
Soil propertiesNatural forestTeaEucalyptusP<  .05pH4.70±0.69a3.72±0.24b4.88±0.85a0.018SOM (%)5.92± 0.846.14± 1.495.26± 1.95NsTotal Nitrogen (%)0.26± 0.060.26± 0.040.26± 0.08NsAv. P (ppm)3.32±0.99a19.40±10.04b12.00±8.70c.002Ca (cmol(p+)/kg)7.63±2.87a3.55±1.98b5.66±2.97a.014Mg (cmol(p+)/kg)2.20±0.38a1.01±0.46b1.57±0.59c<  .001K (cmol(p+)/kg)0.14± 0.030.22± 0.100.33± 0.17NsNa (cmol(p+)/kg)0.042± 0.020.04± 0.010.12± 0.19NsBulk density (g/cm3)0.84±0.03a1.48±0.15b1.40±0.09b.05Ksat (cm/h)15.27± 12.387.00± 12.243.99± 3.21NsClay (%)23.33± 6.5327.00± 7.6731.00± 7.56NsSilt (%)15.00±1.67a16.67±7.12a24.67±10.25bNsTextural classSandy clay loamSandy clay loamSandy clay loamValues followed by the same letters in a row are not significantly different atP<  .05; n.s.: not significant.Likewise, soil properties in the lower depths (15–30 cm) differed significantly (P<  .05) across the land uses. The properties that were similar were clay and silt fractions, composition, total N, SOM, exchangeable sodium, fraction of clay and silt. Comparing soil properties across the LUs, pH, available P and bulk density differed considerably. For instance, soil pH in tea plantations and saturated hydraulic conductivity in natural forests were much lower and higher, respectively than in the other two LUs. As for the top soil (0–15 cm), bulk density was lower in the forest compared to Eucalyptus and Tea.


#### 3.1.2. Interactions Between Land Use and Landscape Positions on Selected Properties in the Top and Subsoil Layers

Within land uses, only soil pH and exchangeable Mg in the top soil depth differed significantly with landscape positions (P<  .05). Soil pH (P=  .014) and exchangeable Mg (P=  .033) were significantly lower in the back slope than foot slope positions. Nonetheless, the other properties particularly SOM (6.08 versus 6.35), available P (15.8 versus 7.6), K (0.33 versus 0.22), and saturated hydraulic conductivity though not significant were lower in the respective slope positions (Table 4 and Figure 3, resp.).Table 4
Effect of land use/cover change on selected physicochemical properties.
Depth (cm)Land useLandscape positionClaySiltOMNAv.PBDCaKMgNapH%ppmMg/m3Cmol (p+)/kg0–15ForestBack slope14.6712.676.420.273.701.696.430.261.670.054.30Foot slope16.6714.676.540.243.001.6512.970.162.850.045.13TeaBack slope17.3314.007.070.2928.101.551.230.400.340.053.30Foot slope29.3312.007.330.3212.901.393.370.180.990.043.77EucalyptusBack slope34.6716.004.750.2515.501.5311.530.342.320.045.53Foot slope37.3316.005.190.227.001.267.130.341.680.054.7315–30 cmForestBack slope19.3314.005.960.302.901.366.000.122.240.034.30Foot slope27.3316.005.890.233.701.529.270.162.160.055.10TeaBack slope24.0013.336.620.2626.601.193.770.280.900.043.63Foot slope30.0020.005.660.2512.201.173.330.161.130.043.80EucalyptusBack slope30.0025.335.810.2516.601.277.900.322.060.045.37Foot slope32.0024.004.710.277.401.173.200.331.080.204.40Land use (LSD0.05)5.004.951.08ns5.310.152.41ns0.43ns0.41Landscape (LSD0.05)4.08nsnsns4.33nsnsnsnsnsnsDepth (LSD0.05)ns4.04nsnsns0.12nsnsnsnsnsLand use*LP(LSD0.05)nsnsnsns7.51ns3.40ns0.60ns0.58Figure 3
Interaction effect of landscape position and depth on saturated hydraulic conductivity in the different land uses.Significant interactions between land use and landscape positions in the lower soil depth were observed for Ca (P=  .016) and exchangeable Mg (P=  .036). The mean values of the two parameters in the back slope were lower than in the foot slope positions within a land use.In addition, Magnesium, Calcium, and Phosphorus were significantly different at the different landscape positions (P<  .05). Also, clay content was significant at the two depths (P<  .05) with depth 15–30 cm generally having higher clay content. Eucalyptus had the highest clay content in both depths although depth 15–30 cm had a slightly lower content than 0–15 cm.


## 3.1. Land Use/Cover Change Maps in 1973, 1986, and 2009

Figure2 and Table 1 show the land use/cover for three time series. Generally small-scale farming, tea plantation, and built up area have increased over time, and woodlot has decreased with time. Previous studies have shown this continuing trend [25].Table 1
showing land use and cover changes.
1973%1987%2009%Built up areas90.5181281.5Grasslands151493361.9Small-scale farming76041125969130871.4Tea plantations1514831789.7Tropical forest828453852125413.8Woodlands20711744291.6Figure 1
Bunyanganbu and Burahya sub-counties in Fort Portal district in Uganda.land use/cover for three time series, 1973, 1987, 2009.
(a)(b)(c)From the figures above, the highest percentage of land in 1973 was tropical forest (45%) compared to small-scale farming in 1987 (69%) and 2009 (71.4%). The smallest percentage of land in 1973 (0.5%) was built up areas compared to tea plantations and grasslands in 1987 (3%) and 2009 (1.9%), respectively. The changes in land use could be attributed to agricultural expansion that took place in the areas adjacent to Kibale National Park [25] as well as an increase in energy demand (firewood). Such expansion was associated with arrival of new immigrants into the areas.


### 3.1.1. Effect of Land Use on Selected Soil Properties

Soil physical and chemical properties in natural forests, tea plantations, and Eucalyptus plantations are presented in Tables2 and 3 for soil depths 0–15 cm and 15–30 cm, respectively. Top soil (0–15 cm) properties except total N, exchangeable potassium, and sodium differed significantly (P<  .05) across the land uses. Textural analysis indicates that the fraction of clay highly significantly (P<  .001) differed between LUs. Clay contents were higher in Eucalyptus than in tea plantations and natural forests. Bulk density was lower in the forest compared to Eucalyptus and Tea. Soil pH (mean = 3.53 ± 0.41) in tea plantations was statistically lower than in natural forests and in Eucalyptus plantations, respectively. Similarly, mean values of SOM were statistically higher in the respective LUs. Total N in natural forests and Eucalyptus were roughly equal.  Available phosphorus in natural forests was about 6 and 3 times lower than in tea plantations and in Eucalyptus plantations, respectively. The highest level of exchangeable Mg was observed in natural forest followed by Eucalyptus plantations and tea plantations, respectively.Table 2
Mean± SD of selected soil (0–15 cm) physical and chemical properties in the three land uses.
Soil propertiesLand use/coverP<  .05Natural forestTeaEucalyptuspH4.72±0.55a3.53±0.41b5.13±0.61c<  .001SOM (%)6.48±0.56a7.20±1.15b4.97±0.45c.002Total N (%)0.26± 0.070.31± 0.050.24± 0.04NsAv. P (ppm)3.32±0.99a20.50±9.52b11.30±10.03a.004Ca (cmol(p+)/kg)9.70±4.18a2.30±1.90b9.30±5.20a.010Mg (cmol(p+)/kg)2.26±0.80a0.66±0.74b2.00±0.55a<  .001K (cmol(p+)/kg)0.21± 0.150.29± 0.190.34± 0.20NsNa (cmol(p+)/Kg)0.04± 0.010.05± 0.010.04± 0.01NsBulk density (g/cm3)0.72±0.05a1.18±0.08b1.22±0.08b.018Ksat (cm/hr)21.90± 24.3519.03± 14.7814.86± 17.84NsClay (%)15.67±4.46a23.33±7.55b36.00±3.35c<  .001Silt (%)13.67± 2.9413.00± 3.0316.00± 2.53NsTextural classSandy LoamSandy clay LoamSandy ClayValues followed by the same letters in a row are not significantly different atP<  .05; n.s.: not significant.Table 3
Mean± SD of selected soil (15–30 cm) physical and chemical properties in the three land use/covers.
Soil propertiesNatural forestTeaEucalyptusP<  .05pH4.70±0.69a3.72±0.24b4.88±0.85a0.018SOM (%)5.92± 0.846.14± 1.495.26± 1.95NsTotal Nitrogen (%)0.26± 0.060.26± 0.040.26± 0.08NsAv. P (ppm)3.32±0.99a19.40±10.04b12.00±8.70c.002Ca (cmol(p+)/kg)7.63±2.87a3.55±1.98b5.66±2.97a.014Mg (cmol(p+)/kg)2.20±0.38a1.01±0.46b1.57±0.59c<  .001K (cmol(p+)/kg)0.14± 0.030.22± 0.100.33± 0.17NsNa (cmol(p+)/kg)0.042± 0.020.04± 0.010.12± 0.19NsBulk density (g/cm3)0.84±0.03a1.48±0.15b1.40±0.09b.05Ksat (cm/h)15.27± 12.387.00± 12.243.99± 3.21NsClay (%)23.33± 6.5327.00± 7.6731.00± 7.56NsSilt (%)15.00±1.67a16.67±7.12a24.67±10.25bNsTextural classSandy clay loamSandy clay loamSandy clay loamValues followed by the same letters in a row are not significantly different atP<  .05; n.s.: not significant.Likewise, soil properties in the lower depths (15–30 cm) differed significantly (P<  .05) across the land uses. The properties that were similar were clay and silt fractions, composition, total N, SOM, exchangeable sodium, fraction of clay and silt. Comparing soil properties across the LUs, pH, available P and bulk density differed considerably. For instance, soil pH in tea plantations and saturated hydraulic conductivity in natural forests were much lower and higher, respectively than in the other two LUs. As for the top soil (0–15 cm), bulk density was lower in the forest compared to Eucalyptus and Tea.


### 3.1.2. Interactions Between Land Use and Landscape Positions on Selected Properties in the Top and Subsoil Layers

Within land uses, only soil pH and exchangeable Mg in the top soil depth differed significantly with landscape positions (P<  .05). Soil pH (P=  .014) and exchangeable Mg (P=  .033) were significantly lower in the back slope than foot slope positions. Nonetheless, the other properties particularly SOM (6.08 versus 6.35), available P (15.8 versus 7.6), K (0.33 versus 0.22), and saturated hydraulic conductivity though not significant were lower in the respective slope positions (Table 4 and Figure 3, resp.).Table 4
Effect of land use/cover change on selected physicochemical properties.
Depth (cm)Land useLandscape positionClaySiltOMNAv.PBDCaKMgNapH%ppmMg/m3Cmol (p+)/kg0–15ForestBack slope14.6712.676.420.273.701.696.430.261.670.054.30Foot slope16.6714.676.540.243.001.6512.970.162.850.045.13TeaBack slope17.3314.007.070.2928.101.551.230.400.340.053.30Foot slope29.3312.007.330.3212.901.393.370.180.990.043.77EucalyptusBack slope34.6716.004.750.2515.501.5311.530.342.320.045.53Foot slope37.3316.005.190.227.001.267.130.341.680.054.7315–30 cmForestBack slope19.3314.005.960.302.901.366.000.122.240.034.30Foot slope27.3316.005.890.233.701.529.270.162.160.055.10TeaBack slope24.0013.336.620.2626.601.193.770.280.900.043.63Foot slope30.0020.005.660.2512.201.173.330.161.130.043.80EucalyptusBack slope30.0025.335.810.2516.601.277.900.322.060.045.37Foot slope32.0024.004.710.277.401.173.200.331.080.204.40Land use (LSD0.05)5.004.951.08ns5.310.152.41ns0.43ns0.41Landscape (LSD0.05)4.08nsnsns4.33nsnsnsnsnsnsDepth (LSD0.05)ns4.04nsnsns0.12nsnsnsnsnsLand use*LP(LSD0.05)nsnsnsns7.51ns3.40ns0.60ns0.58Figure 3
Interaction effect of landscape position and depth on saturated hydraulic conductivity in the different land uses.Significant interactions between land use and landscape positions in the lower soil depth were observed for Ca (P=  .016) and exchangeable Mg (P=  .036). The mean values of the two parameters in the back slope were lower than in the foot slope positions within a land use.In addition, Magnesium, Calcium, and Phosphorus were significantly different at the different landscape positions (P<  .05). Also, clay content was significant at the two depths (P<  .05) with depth 15–30 cm generally having higher clay content. Eucalyptus had the highest clay content in both depths although depth 15–30 cm had a slightly lower content than 0–15 cm.


## 3.1.1. Effect of Land Use on Selected Soil Properties

Soil physical and chemical properties in natural forests, tea plantations, and Eucalyptus plantations are presented in Tables2 and 3 for soil depths 0–15 cm and 15–30 cm, respectively. Top soil (0–15 cm) properties except total N, exchangeable potassium, and sodium differed significantly (P<  .05) across the land uses. Textural analysis indicates that the fraction of clay highly significantly (P<  .001) differed between LUs. Clay contents were higher in Eucalyptus than in tea plantations and natural forests. Bulk density was lower in the forest compared to Eucalyptus and Tea. Soil pH (mean = 3.53 ± 0.41) in tea plantations was statistically lower than in natural forests and in Eucalyptus plantations, respectively. Similarly, mean values of SOM were statistically higher in the respective LUs. Total N in natural forests and Eucalyptus were roughly equal.  Available phosphorus in natural forests was about 6 and 3 times lower than in tea plantations and in Eucalyptus plantations, respectively. The highest level of exchangeable Mg was observed in natural forest followed by Eucalyptus plantations and tea plantations, respectively.Table 2
Mean± SD of selected soil (0–15 cm) physical and chemical properties in the three land uses.
Soil propertiesLand use/coverP<  .05Natural forestTeaEucalyptuspH4.72±0.55a3.53±0.41b5.13±0.61c<  .001SOM (%)6.48±0.56a7.20±1.15b4.97±0.45c.002Total N (%)0.26± 0.070.31± 0.050.24± 0.04NsAv. P (ppm)3.32±0.99a20.50±9.52b11.30±10.03a.004Ca (cmol(p+)/kg)9.70±4.18a2.30±1.90b9.30±5.20a.010Mg (cmol(p+)/kg)2.26±0.80a0.66±0.74b2.00±0.55a<  .001K (cmol(p+)/kg)0.21± 0.150.29± 0.190.34± 0.20NsNa (cmol(p+)/Kg)0.04± 0.010.05± 0.010.04± 0.01NsBulk density (g/cm3)0.72±0.05a1.18±0.08b1.22±0.08b.018Ksat (cm/hr)21.90± 24.3519.03± 14.7814.86± 17.84NsClay (%)15.67±4.46a23.33±7.55b36.00±3.35c<  .001Silt (%)13.67± 2.9413.00± 3.0316.00± 2.53NsTextural classSandy LoamSandy clay LoamSandy ClayValues followed by the same letters in a row are not significantly different atP<  .05; n.s.: not significant.Table 3
Mean± SD of selected soil (15–30 cm) physical and chemical properties in the three land use/covers.
Soil propertiesNatural forestTeaEucalyptusP<  .05pH4.70±0.69a3.72±0.24b4.88±0.85a0.018SOM (%)5.92± 0.846.14± 1.495.26± 1.95NsTotal Nitrogen (%)0.26± 0.060.26± 0.040.26± 0.08NsAv. P (ppm)3.32±0.99a19.40±10.04b12.00±8.70c.002Ca (cmol(p+)/kg)7.63±2.87a3.55±1.98b5.66±2.97a.014Mg (cmol(p+)/kg)2.20±0.38a1.01±0.46b1.57±0.59c<  .001K (cmol(p+)/kg)0.14± 0.030.22± 0.100.33± 0.17NsNa (cmol(p+)/kg)0.042± 0.020.04± 0.010.12± 0.19NsBulk density (g/cm3)0.84±0.03a1.48±0.15b1.40±0.09b.05Ksat (cm/h)15.27± 12.387.00± 12.243.99± 3.21NsClay (%)23.33± 6.5327.00± 7.6731.00± 7.56NsSilt (%)15.00±1.67a16.67±7.12a24.67±10.25bNsTextural classSandy clay loamSandy clay loamSandy clay loamValues followed by the same letters in a row are not significantly different atP<  .05; n.s.: not significant.Likewise, soil properties in the lower depths (15–30 cm) differed significantly (P<  .05) across the land uses. The properties that were similar were clay and silt fractions, composition, total N, SOM, exchangeable sodium, fraction of clay and silt. Comparing soil properties across the LUs, pH, available P and bulk density differed considerably. For instance, soil pH in tea plantations and saturated hydraulic conductivity in natural forests were much lower and higher, respectively than in the other two LUs. As for the top soil (0–15 cm), bulk density was lower in the forest compared to Eucalyptus and Tea.


## 3.1.2. Interactions Between Land Use and Landscape Positions on Selected Properties in the Top and Subsoil Layers

Within land uses, only soil pH and exchangeable Mg in the top soil depth differed significantly with landscape positions (P<  .05). Soil pH (P=  .014) and exchangeable Mg (P=  .033) were significantly lower in the back slope than foot slope positions. Nonetheless, the other properties particularly SOM (6.08 versus 6.35), available P (15.8 versus 7.6), K (0.33 versus 0.22), and saturated hydraulic conductivity though not significant were lower in the respective slope positions (Table 4 and Figure 3, resp.).Table 4
Effect of land use/cover change on selected physicochemical properties.
Depth (cm)Land useLandscape positionClaySiltOMNAv.PBDCaKMgNapH%ppmMg/m3Cmol (p+)/kg0–15ForestBack slope14.6712.676.420.273.701.696.430.261.670.054.30Foot slope16.6714.676.540.243.001.6512.970.162.850.045.13TeaBack slope17.3314.007.070.2928.101.551.230.400.340.053.30Foot slope29.3312.007.330.3212.901.393.370.180.990.043.77EucalyptusBack slope34.6716.004.750.2515.501.5311.530.342.320.045.53Foot slope37.3316.005.190.227.001.267.130.341.680.054.7315–30 cmForestBack slope19.3314.005.960.302.901.366.000.122.240.034.30Foot slope27.3316.005.890.233.701.529.270.162.160.055.10TeaBack slope24.0013.336.620.2626.601.193.770.280.900.043.63Foot slope30.0020.005.660.2512.201.173.330.161.130.043.80EucalyptusBack slope30.0025.335.810.2516.601.277.900.322.060.045.37Foot slope32.0024.004.710.277.401.173.200.331.080.204.40Land use (LSD0.05)5.004.951.08ns5.310.152.41ns0.43ns0.41Landscape (LSD0.05)4.08nsnsns4.33nsnsnsnsnsnsDepth (LSD0.05)ns4.04nsnsns0.12nsnsnsnsnsLand use*LP(LSD0.05)nsnsnsns7.51ns3.40ns0.60ns0.58Figure 3
Interaction effect of landscape position and depth on saturated hydraulic conductivity in the different land uses.Significant interactions between land use and landscape positions in the lower soil depth were observed for Ca (P=  .016) and exchangeable Mg (P=  .036). The mean values of the two parameters in the back slope were lower than in the foot slope positions within a land use.In addition, Magnesium, Calcium, and Phosphorus were significantly different at the different landscape positions (P<  .05). Also, clay content was significant at the two depths (P<  .05) with depth 15–30 cm generally having higher clay content. Eucalyptus had the highest clay content in both depths although depth 15–30 cm had a slightly lower content than 0–15 cm.


## 4. Discussion

Although mean soil pH in all the LUs is lower than the critical value of 5.5 [26], soils from only tea plantations are strongly acidic, pH less than 4.5 [27]. The lowest pH, though corresponds to high SOM contents in the latter LU, could partly be attributed to the long-term effect of fertilizer application as it was also observed in analogous; total N in tea plantations was much above the critical value of 0.2%. Similar findings were observed by Wang et al. [28] in commercial tea plantations in China that received mostly nitrogen fertilizers.Soil organic matter in all the LUs is well above 3.0% cited by Okalebo et al. [22] as the critical value. The observed differences in SOM seem to relate to the level of input of organic materials which was higher (thickness of debris/litter was qualitatively higher) in tea plantations than natural forests.  Residues from regular pruning and leaf fall were common in tea plantations and combined with application of inorganic fertilizer-enhanced decomposition and mineralization of plant residues compared to natural conditions [29]. Arguably, OM turnover is higher in tea plantations than in natural forests. Consequently, available P was in only tea plantations higher than the critical of 15 ppm [22]. According to Brady [27], the bulk of phosphorus exists in organic form and P is very stable [30] in the soil, thus a good indicator of organic matter accumulation and quality; Exchangeable Ca was below the critical value of 4.0% only in tea plantations. Except in natural forests, the mean value of exchangeable Mg in other land uses is slightly below the critical values of 0.5 cmol(p+)/kg.Interpretations of values of hydraulic conductivity with ranges established by Landon [26] indicate that water flow in soils of all the LUs is rapid (Ksat > 15 cm/hr). According to Lichner et al. [31], roots can affect saturated hydraulic conductivity in topsoil directly by inducing the preferred flow through macro-pores or indirectly by their influence on soil structure. Soil texture under the three different land use covers was dominantly sandy-soil type which is characterized by high-saturated hydraulic conductivity.The lower values of soil properties in back slope than foot slope within land uses could be related to deposition of eroded materials which occurs in the latter landscape [19]. Steeper landscapes are generally more prone to detachment/loss of soil materials by water erosion [32].


## 5. Conclusions

Small-scale farming, tea and eucalyptus plantation, and built up area have increased over time, at the expense of the woodlot and forest cover.The change from natural forest cover to tea and Eucalyptus induced changes in top soil properties namely, exchangeable Mg and Ca, available P, SOM, pH, and bulk density of sub soil. Landscape positions within land use also significantly influenced most of the soil properties denoting interactions. There is a need therefore of assessing the spatiotemporal consequence of these changes at watershed level, carrying out an economic analysis of these activities at household level, and identifying the most beneficial Kibale National Park trees for afforestation.



---

*Source: 185689-2010-06-20.xml*