Vertical distribution of soil chemical properties under different land use systems of Biligirirangana Hills, Karnataka

Govinda. K , Krishna Murthy, R , Bhavya., N , Uday Kumar S. N , Pushpa K

All India Coordinated Research Project on Soil Test Crop Response University of Agricultural Sciences, Bangalore

Corresponding Author Email: stcruasbangalore@gmail.com

DOI : https://doi.org/10.61739/TBF.2023.12.2.260

Keywords

available nitrogen, Soil organic carbon, Soil secondary and micronutrients

Download this article as:

Abstract

Globally, efforts are being made to identify land use types that could potentially improve carbon sequestration and soil fertility to ensure sustainable agriculture. Secondary nutrients and micronutrients in soil are crucial for improved crop output and require a detailed understanding of their status and sources of variation. Therefore, in order to evaluate the changes in soil chemical properties of soils, an experiment was conducted in Biligirirangana hills with soil samples from three soil depths i.e. 0-15, 15–30, and 30–45 cm of six prominent land-use systems viz. natural forest, a forest plantation, pure coffee system, coffee multi-storeyed system, agriculture mono, and mixed cropping systems. Among different land use systems (LUSs), significantly higher mean soil organic carbon and available nitrogen content were recorded in surface soils of natural forest (37.09 g kg-1) and pure coffee system (414.99 kg ha-1), respectively. The cation exchange capacity and percent base saturation were significantly higher in natural forest (28.22 cmol (p+) kg-1 soil and 76.56 % respectively) and lower in agriculture mono-cropping (13.83 cmol (p+) kg-1 soil and 39.50 %, respectively). The exchangeable calcium, magnesium, available sulfur and DTPA extractable micronutrients were higher in natural forest and tree-based LUSss compare to agriculture LUSss. The correlation analysis revealed that soil SOC had a positive effect on the availability of available secondary and DTPA extractable micronutrients in soils.

Introduction

The increasing population and its corresponding demand for food in the midst of changing climate require the identification of land use types and management practices that can promote carbon sequestration and retention for agricultural sustainability and food security (Okolo et al, 2019). Globally, land-use changes associated with forestry and agriculture, such as deforestation and subsequent conversion to agriculture, are estimated to be responsible for 25 percent of the accumulated GHG (greenhouse gas) emissions (Achard et al., 2014). Despite 37% of total land globally being used for agriculture (Chen et al 2018, and World Bank 2019), the world still grapples with food production on a sustainable basis due to soil and land degradation, rising sea levels, climate change, and increasing population growth (F.A.O, 2016, and Trendov, et al, 2019). This could affect the soil properties of a given land when it is changed from natural to converted LUSs or vice versa.

In the Indian context, around 275 million people rely on forests for their livelihood or means of subsistence (World Bank 2006). Of these, around 89 million people belong to the marginalized communities known as Adivasis or Scheduled tribes. Soligas, an Adivasi community residing in the Biligirirangana Hills located in South-Western Karnataka are the unique extension of the Western Ghats. Shifting cultivation and more recently settled agriculture have historically been important sources of subsistence for these communities (Mundoli et al., 2016). The Soliga community relocated to a stable agricultural area where the same plot of land was maintained and farmed continuously when this place was designated as a Sanctuary in 1974. One of the main causes of changes in land use and livelihood in the tropical forests is due to government policies that supported settled agriculture, combined with an ever-increasing demand for forest products from international markets, a changing climate, limited access to markets, and rising wildlife depredation, have increased rates of conversion of forests to plantation crops (Mundoli et al., 2016).

            The land use type could be determined by the need of the producer, the environmental condition (soil, climate, rainfall, altitude, etc.), socioeconomic status (landlord, tiller, or peasant), and political (tenure, land policy, and ownership) and cultural manners (beliefs, norms, and bylaws on the land) of the given area (DeFries et al., 2004). Land use changes affect basic processes such as erosion, soil structure, nutrient recycling, and carbon sequestration (Yoseph et al. 2017) and it greatly influences physicochemical properties (Paz-Kagan, et al. 2014) and affects the nutrient dynamics and supply (Agniva et al. 2018). Secondary (Ca, Mg, and S) and  Micronutrients (Fe, Mn, Zn, and Cu) plays a vital role in the growth and production of crops through their involvement in regulating enzyme, cell wall formation, synthesis of chlorophyll, proteins, lipids, nucleic acids, carbohydrates metabolism and tolerance to stress, etc (Barker and Pilbeam, 2015). Soil secondary nutrients and micronutrient cycling and their distribution in soils have been known to be varied with a shift in the land use from forest to cultivated systems and differences in land use patterns and depth of soil (Han, et al, 2017, Paul et al., 2018). The status of available micronutrients in soils depends on soil properties such as pH, EC, OC, clay and calcium carbonate (CaCO3) content, land use patterns, and soil depths (Shukla et al., 2014).

            Information soil secondary and micronutrient status when natural forest converted to different LUSss at different soil depths.  We get the chance to acquire a deeper understanding of the aforementioned subject in forest-dwelling Adivasi farmers situated within a protected tiger reserve (Biligirirangana Hills). Our major objectives of the current study were to investigate 1) soil available secondary and micronutrients for different land use systems during the conversion of natural forest to agricultural fields.

2. Materials and methods

2.1. Description of the study site

Biligirirangana Hills is located in South-Western, Karnataka and it is the unique extension of the Western Ghats. The region constitutes a live bridge between the Eastern and Western Ghats. The region extends between 11⁰ 40′ – 12⁰ 06′ N latitude and 76⁰ 24′ – 77⁰ 46′ E longitude and the altitude ranges from 1400 – 1800 m above MSL. A location map of the study area is presented in Figure. 1.  The Biligiris are Charnockite hills, covered with tropical dry broadleaf forest, part of the South Deccan Plateau dry deciduous forests ecoregion. The forests range from scrub forests at lower elevations, degraded by over-use, to the tall deciduous forests typical of the ecoregion, to stunted shola forests and montane grasslands at the highest elevations, which exceed 1800 meters. The region receives an annual rainfall of 1099 mm with the highest precipitation occurring in the month of September. The maximum and minimum temperature varied from 13.6 to 21.3⁰ C and 28.3 to 36.3 ⁰C, respectively (Figure 2).

2.2. Selection of land use types

The dominant land use types selected were natural forest, a forest plantation, pure coffee system, coffee multi-storeyed system, agriculture mono, and mixed cropping systems. The details of the land use types have been described in Table 1.

2.3. Land demarcation and soil sampling

Garmin eTrex 30x (global positioning system (GPS) device was used to record the boundary coordinates of the entire study site and each land use type. A stratified random sampling technique was used for soil sampling where GPS coordinates were recorded at each sampling point. In each LUS, ten mini pits were randomly dug and marked at 0–15 cm, 15–30 cm, and 30–45 cm depths. At each fixed soil depth, 1 soil sample were randomly collected resulting in a total soil sample of 180 (1 sample x 3 soil depths x 10 pits × 6 land use types) and then kept in labeled transparent plastic zip lock bags and sealed for laboratory analysis. These samples collected across the entire field were used for the determination of the spatial variation and vertical distribution of soil chemical properties.

2.4. Laboratory analysis of soil samples

Plant materials and clods were removed from the air-dried soil samples, crushed using a mortar with a pestle, and sieved through a 2 mm mesh sieve. The chemical properties analyzed in the laboratory soil organic carbon, cation exchange capacity, secondary and micronutrient contents. The determination of soil organic carbon was carried out by wet oxidation method (Walkley and Black, 1934), and available nitrogen by Subbaiah and Asija (1956). The neutral normal ammonium acetate extract was used to determine exchangeable calcium and magnesium content (Page et al., 1982). Sulfur was extracted using a 0.15 percent calcium chloride solution, and the amount of sulfur in the extract was determined using turbidimetry as described by Black (1965). The available micronutrients in soil were extracted with DTPA-extractant (0.005 M Diethylene Triamine Penta Acetic acid + 0.01 M CaCl2 + 0.1 M Triethanol Amine) at a 1:2 soil-to-extractant ratio as described by Lindsay and Norvell (1978). The concentrations of zinc, iron, manganese, and copper in the extract were determined by an Atomic Absorption Spectrophotometer (AAS) under suitable measuring conditions (Page et al., 1982). The cation exchange capacity of the soil was estimated by the neutral normal ammonium acetate method Jackson (1973).

2.5 Statistical analysis

The statistical analysis of the experimental data gathered during the study was performed using the suitable analysis of variance (ANOVA) approach to examine the significance of the overall differences across treatments using the “F” test as described by Gomez and Gomez (1984). Using SPSS software, the data was statistically examined for different soil properties using a factorial approach (version 1.2.5042). When the “F” value was found to be significant, the Tukey HSD process was used to compare the means of different land uses and soil depths (Steel and Torrie, 1960). At the (P 0.05) level of significance, differences in land use mean were judged statistically significant.

3 Result and discussion

3.1 Chemical properties of soils under different land use systems

3.1.1 Soil organic carbon

The results of soil organic carbon (SOC) content of soil at different depths under different LUSs are presented in Table 2. The soil organic carbon content differed significantly with LUSs and depth. The soil organic carbon content of the surface soil was significantly higher (29.84 g kg-1) than the lower soil layer (17.33 g kg-1). Among different LUSs, the mean SOC content under the natural forest system (37.09 g kg-1) was significantly high followed by forest plantation (31.82 g kg-1). The mean SOC content in soils under a coffee multi-storeyed system (26.03 g kg-1) was statistically on par with the pure coffee system (25.55 g kg-1) but significantly higher than agricultural LUSs (Table 2). The interaction between the land use system and depth was found to be significant. A significantly higher SOC content of 45.9 g kg-1 was recorded in the surface soil layer under natural forest type than the values found in the surface soil layer of other systems. Similarly, a significantly lower SOC content of 10.24 g kg-1 was found in agriculture mono-cropping systems.

Higher SOC content in natural forest, forest plantation, and coffee LUSs was mainly attributed to addition of larger amounts of litter biomass on surface soil and low oxidation process as the soils are less disturbed. Similar, higher SOC content in forest soils has been reported by Nagaraja (1997) and Yao et al. (2010).  Sainepo et al. (2018) revealed different land use types had significantly different SOC where shrubland exhibited high levels of SOC which could be attributed to higher litter deposition, soil moisture content, C: N ratio and grazing management.

            The higher SOC content in soils under coffee could be attributed to the addition of litter and pruned plant materials from shade trees and coffee plants and other management practices involving the addition of organic manures and fertilizers over the years helps to maintain the higher organic carbon level than other agriculture LUSs. Whereas, significantly lower SOC content was recorded in soils under the agriculture mono-cropping system (8.5 g kg-1) than all other systems. The variation in SOC content among agriculture LUSs compared to other land uses could be attributed to variations in litter input and amount of FYM used, root growth, the extent of surface cover and soil erosion, type and time of tillage operation and addition of biomass and intensive cultivation practices (Doran 2002). The accumulation of organic matter in soils is the result of a balance between residue input rates and net decomposition rates as affected by soil temperature and moisture conditions (Murage et al., 2007). These results are corroborated by the findings of Shivakumar et al. (2020). Islam and Well (2000) reported that lower SOC in soil under paddy might be attributed to an imbalance in C input and output in the system. A decrease in soil organic carbon is observed as the impact of conversion from natural forest to coffee plantation (Karthikaet al. 2020;; Dariah et al., 2008).

The variation in soil organic carbon content may be attributed to variation in potential carbon returned to soil through the litter. This was evident in the significant positive correlation observed in the present investigation among litter C content and SOC content which suggest that as the plant carbon in litter increases the soil organic carbon.   Takele et al. (2014) reported that the SOC decreases with increasing soil depth, with more accumulation on the upper surface soil layer. Toru and Kibret (2019) found that soil organic carbon and total nitrogen were significantly (p<0.05) affected by soil depth and land use under four major land uses viz., natural forest, coffee agroforestry, grazing land, and cropland in Hades sub-watershed, eastern Ethiopia.

3.1.2 Available nitrogen

The data on the available nitrogen status of soils as influenced by different LUSs at different depths are presented in Table 2. The difference in available nitrogen content with depth was significant. The soil available nitrogen content in surface soil was significantly higher (424.41 kg ha-1) than that of subsurface soil (275.97 kg ha-1). Among LUSs, the available nitrogen content was significantly higher in the pure coffee system (414.99 kg ha-1) followed by coffee multi-storeyed system (410.81 kg ha-1), natural forest (393.04 kg ha-1) than other systems. Similarly, the lowest available N content of 263.42 kg ha-1 was recorded in the agriculture mono-cropping system and it was statistically on par with agriculture mixed cropping system (270.74 kg ha-1). The interaction between depth and land use was significant. Thus, the maximum value of 495.49 kg ha-1 was recorded in surface soil of the coffee multi-storeyed system and the minimum value (185.02 kg ha-1) was recorded in the subsurface layer (30-45 cm) of the agriculture mono-cropping system.

Soil organic matter (SOM) acts as a storehouse and supplier of N to plant roots and microorganisms, almost 95 percent of the total soil N is closely associated with SOM. The higher available N content in manmade systems especially under coffee plantations might be attributed to the application of very high dose of N fertilizers along with very high level of FYM. Similarly, the higher available N content under natural forest and forest plantation system might be attributed to the variation in quality and quantity of litter input and rate of mineralization. Secondly, the leguminous tree species in the natural forest might have fixed atmospheric N. Thus, the higher available N content was observed under natural forests. The amount so fixed along N derived from mineralization might be sufficient to replace the N removed through the litter. These results are in accordance with Shivakumar et al. (2020) who reported that higher available N content under natural systems namely evergreen and semi-evergreen compared to forest plantation systems viz., teak and acacia might be attributed to the variation in quality and quantity of litter as it is derived from mixed vegetation type under natural system.

            Further, the removal of vegetation by livestock grazing and exposure of soil surface may impact on the surface runoff due to direct drop of rainfall under agricultural land use system, which can remove the animal and plant residues from the surface soil layer thereby causing nitrogen depletion. Similarly, Nigussie and Kissi (2012), Ufot et al. (2016), Tufa et al. (2019) and Chemada et al. (2017) stated that the higher total N was obtained under forest land compared to the adjacent grazing and cultivated lands. 

            The soil available nitrogen content in surface soil was significantly higher (424.41 kg ha-1) than that of subsurface soil (275.97 kg ha-1) which may be attributed to the accumulation of more organic matter on the surface layer and its subsequent decomposition than the subsurface soil layer. These results are in conformity with the results obtained by Shivakumar (2020) explaining the presence of higher available nitrogen content in the surface layer than that of the subsurface layer of different LUSs. The contents of available nitrogen were strongly associated (r =0.875*) with soil organic carbon in all the depths and decreased consistently with increasing soil depth under all LUSs (Table 6). The results of the present study agree with the findings of many other workers Malo et al. (2005) and Heluf and Wakene (2006) reported that the decrease in available nitrogen with increasing depth was due to a decline in humus content with depth.

3.1.3 Cation exchange capacity

The data presented in Table 3 indicates that among different LUSs mean cation exchange capacity (CEC) values of soil were higher in natural forest (28.22 cmol (p+) kg-1 soil)  and lower in agriculture mono – cropping (13.83 cmol (p+) kg-1 soil), while CEC of pure coffee and coffee multi-storeyed system was on par with each other. Among different soil depths, cation exchange capacity of surface soil depth (0-15 cm) were significantly higher (23.38 cmol (p+) kg-1 soil) than lower soil depth (30-45 cm) (18.72 cmol (p+) kg-1 soil). The CEC of soil differed significantly with land use (Table 3). Higher CEC of soil was recorded in the surface depth of natural forest and tree-based LUSs compare to agricultural LUSs, which might be due to high exchangeable bases and organic carbon resultin in  the dissociation  of   its   functional   groups  which providing   more negative charges (Das  et al., 2023: Karthika et al. (2020)). Shivakumar et al. (2020) reported that CEC of soil under the coffee system was significantly higher as compared to paddy LUSs. Das et al., 2023 and Shivakumar et al. (2020) observed a decline of CEC with soil depth along with a decline in SOM but unchanged clay content and composition.

NS: Non significant, 

3.1.4 Base saturation (%)

  The percent base saturation (PBS) values of soils were significantly affected by land use types but not among soil depths (Table 3). The PBS values range from  39.50 to 76.56 % from all the LUSs. The study site is situated in high rainfall area and these soils are basically derived from granite and granite gneiss which are silica-saturated igneous and metamorphic rocks, as a result the soils shows lower PBS (Anonymous, 2005). The highest mean value of PBS was found in natural forest (76.56 %) which was on par with forest and coffee plantations, and the lowest was observed in agriculture mono-cropping system (39.50 %). In all the LUSs, the PBS is significantly increased with increasing depths and the highest was recorded in surface soil depth (60.25 %). Natural forest and tree-based LUSs’ surface soil depths recorded greater percentages of base saturation; these results were comparable to those of the soil’s CEC and soil organic carbon content.  Bouajila, et al. (2023) reported that FYM and Sulla amendment resulted in PBS increased up to 73.3% and 63.8% respectively compared to the soil control (50.6%). Relatively lower percent base saturation was recorded in agricultural LUSs because of nutrient removal resulting from forest clearance and less organic input with intensive tillage operation (Chimdie and Gurmessa, 2023).

3.1.5 Calcium and magnesium saturation (%)

Conversion of natural forest to different LUSs had a significant influence on the calcium and magnesium saturation in soil. The highest calcium and magnesium saturation were recorded in the surface depth of pure coffee (61.61%) and coffee multi storeyed systems (18.04%) when compared to other LUSs (Figure 3 and 4). Coffee and forest plantations, natural forest recorded higher Ca and Mg saturation values compare to agricultural LUSs, which might be due to the organic carbon of soil which increase cation holding capacity by increasing the CEC of soil  (Chimdie and Gurmessa, 2023).

3.2 Secondary nutrients

 3.2.1 Exchangeable calcium

The data on exchangeable-Ca content in soils under different LUSs are presented in Table 4. Calcium content in soils differed significantly due to the land use system. The significantly higher calcium content (10.62 cmol (p+) kg-1) was recorded in natural forest soils and which was on par with the coffee multi-storeyed system (10.38 cmol (p+) kg-1). The highest exchangeable calcium content (10.62 cmol (p+) kg-1 was recorded in natural forests (Table 5) compare to agriculture, pure coffee system may be attributed to higher turnover of organic matter having higher Ca, and its mineralization was the cause for the higher available Ca in soils. The recycling of litter having a higher amount of Ca keeps this element in continuous circulation in the system. These results are in accordance with the findings of Geetha (2005) reported that the exchangeable calcium content in surface soils under high rainfall areas was highest in evergreen forest vegetation followed by moist deciduous forest. However, the quantity of calcium content in soils under plantations such as teak systems is comparable with each other. In coffee plantations, the exchangeable calcium content was higher may be attributed to the recycling of higher litter and the application of agricultural lime and dolomite (CaCO3MgCO3) in coffee plantations might have enhanced exchangeable (Shivakumar et al. 2020) and Karthika et al. (2020) found that coffee soils recorded relatively exchangeable Ca and Mg than forest soils in all the profiles. The calcium content in soil was positively correlated with SOC (r = 0.296), CEC (r = 0.418) and PBS (r = 0.629) in all the soil depths (Table 6).

The lowest (3.88 cmol (p+) kg-1) value exchangeable-Ca was recorded in the agriculture mixed cropping system which was on par with the agriculture mono-cropping system (4.46 cmol (p+) kg-1) which might be due to erosion loss and variation in litter input and rate of decay and leaching. Secondly, the soils in the hill zone are predominantly kaolinitic and consequently had low negative charge density to hold cation in exchange surface. Further, these cations are not adsorbed; they are susceptible for leaching in high-rainfall areas. Besides, the variation may also be attributed to management practices adopted to grow these crops. Lima et al. (2022b), analyzing the forest-to-pasture conversion, also found low levels of Ca2+, Mg2+ and K+ in the first layers of soil in cultivated areas. Prabhudev (2011) has also opined that leaching is the main cause of low availability.  

Calcium content in soils differed significantly with depth. Surface (0-15 cm) soil had higher NH4OAC – Ca content (9.34 cmol (p+) kg-1) than subsurface (30-45 cm) soil (7.15 cmol (p+) kg-1). The variation in depth could be attributed to low-activity clay and organic matter content in the subsurface soil. Arevalo et al. (2009) stated that the increase in exchangeable Ca and Mg with soil depth were likely to have been brought about by weathering and leaching. Constant mixing of the soil by ploughing and disking of the soil surface release more exchangeable Ca and Mg, which are leached to subsurface layers with precipitation and infiltration. Salve et al. (2018) reported that exchangeable Ca (5.52 mg/100g) was significantly higher in agri–the horti-silviculture system and at 0–15 cm depth (6.03 mg/100g).    Interaction between depth and land use showed a significant difference and the highest value of 11.90 cmol (p+) kg-1 was recorded in the surface soil of natural forest followed by coffee multi-storeyed system (11.79 cmol (p+) kg-1).

3.2.2 Exchangeable magnesium

            The data relating to exchangeable Mg in soils under different LUSs is presented in Table 4. The concentration of NH4OAC-Mg varied significantly between LUSs, ranging from 3.11 to 5.54 cmol (p+) kg-1. Soils under the forest plantation had a higher Mg content of 5.54 cmol (p+) kg-1, which was significantly higher than soils under natural forest (3.94 cmol (p+) kg-1), agriculture mono-cropping (3.28 cmol (p+) kg-1), and agriculture mixed cropping (3.11 cmol (p+) kg-1), and other systems were comparable. There was no significant variation in exchangeableMg with depth, but surface soils (4.62 cmol (p+) kg-1) contained higher exchangeableMg than subsurface (30-45 cm) samples (4.12 cmol (p+) kg-1). Furthermore, there was no significant variation between land use and depth.

Higher exchangeableMg content was recorded in forest and coffee plantations compared to agricultural LUSs. This kind of variation was mainly attributed to management practices being practiced and the quantity and quality of litter input added to the particular land use. Application of dolomite (CaCO3MgCO3), for instance, in coffee plantations might have enhanced exchangeableMg content in soil under these crops. Similar observations were also reported by Prabhudev (2011). While agriculture’s land use system has a lower amount of organic matter, intensive cultivation, leaching and erosion loss of exchangeable Mg through percolating water might have reduced Mg content (Shivakumar et al., 2020: Lima et al. (2022b).

3.2.3 Available sulphur

            The data on depth-wise changes in the available sulphur status of soils under different LUSs are presented in Table.4.

The available sulfur content differed significantly with depth, so the available sulphur value (49.62 mg kg-1) of surface soil was higher than that of subsurface (30-45 cm) layer soil (25.73 mg kg-1). Similarly, the available sulfur content differed significantly among LUSs. The available sulfur content in soils under multi-storeyed coffee soils was significantly higher (49.31 mg kg-1) than in pure coffee system (46.79 mg kg-1), natural forest (43.62 mg kg-1), forest plantation (39.12 mg kg-1), agriculture mixed system (25.21 mg kg-1) and agriculture mono-cropping system (18.81 mg kg-1).

In the interaction of land use and depth, a maximum available sulfur content of 62.96 mg kg-1 was observed in the surface soil layer of natural and a minimum available sulphur content of 12.26 mg kg-1 was observed in the subsurface (30-45 cm) layer of an agriculture mono-cropping system.

When compared to agricultural mono and mixed systems, the available sulfur level in soils under coffee plantations was much higher, followed by natural forest (Table 4). The reason for this could be due to increasing soil organic carbon content (Table 4), which reduces sulfate ion leaching. Second, the use of S-containing fertilizers (single super phosphate or ammonium

sulfate) and CuSO4 (bordeaux mixture) could be contributing to the elevated S levels in that system. In coffee plantations, Shivakumar et al. (2020) and Prabhudev (2011) observed comparable results in terms of S availability. Herojith Singh et al. (2007), with the higher available sulfur content ascribed to higher organic matter content.

3.4 DTPA extractable micronutrients (mg kg-1)

3.4.1 DTPA-Fe

DTPA extractable Fe content in soils differed significantly with depth and LUSs (Table 5). The surface soil layer had significantly higher DTPA-Fe levels (33.38 mg kg-1) than the subsurface soil layer (30-45 cm) (24.24 mg kg-1). While the amount of DTPA-Fe was significantly higher in soils under forest plantation (42.32 mg kg-1) than in soils under natural forest (36.80 mg kg-1), coffee multi-storeyed system (31.98 mg kg-1), pure coffee system (28.85 mg kg-1), agriculture mixed cropping (15.92 mg kg-1). The amount of DTPA-Fe was significantly lower in agriculture mono-cropping (15.13 mg kg-1) compare to other land uses. Kumar et al. (2009) showed that DTPA–extractable Fe and Mn (but not Zn and Cu) in soil declined from their respective initial values as a result of intensive cropping for more than three decades in long‐term experiments in New Delhi, India. The interaction effect of depth x land use was non-significant, thus the higher value of 46.88 mg kg-1 DTPA-Fe was recorded in surface soils of forest plantation and lower value was recorded in subsurface soils under agriculture mono-cropping (11.22 mg kg-1).

3.4.2 DTPA-Mn

The data presented in Table 5 revealed that DTPA-Mn content in soils differed significantly with depth and LUSs. However, the effect of the season was non-significant.

Significantly lower (20.17 mg kg-1) DTPA-Mn content was recorded in the subsurface (30–45 cm) soil layer as compared to the surface soil layer (25.76 mg kg-1). However, the amount of DTPA-Mn was significantly higher in soils under forest plantation (28.33 mg kg-1), which was on par with natural forest (27.98 mg kg-1) and pure coffee systems (26.12 mg kg-1) followed by the coffee multi-storeyed system (22.99 mg kg-1). A significantly lower amount of 14.86 mg kg-1 was recorded in soils under the agriculture mono-cropping system, which was on par with the agriculture mixed cropping system (16.27 mg kg-1).

The interaction effect of depth and land use was non-significant. However, a numerically higher value of 31.68 mg kg-1 DTPA-Mn was recorded in surface soils of pure coffee plantations and lower value was recorded in subsurface (30-45 cm) soils under agriculture mono-cropping (12.04 mg kg-1).

 3.4.3 DTPA-Zn

The data on the DTPA-Zn content of soils under different LUSs are presented in Table 5. The soil DTPA-Zn showed a significant difference among LUSs. The highest mean DTPA-Zn content was found in the coffee multi-storeyed system (2.71 mg kg-1), followed by the pure coffee system (2.60 mg kg-1), and forest plantation (2.43 mg kg-1). A significantly lower amount of mean DTPA-Zn was found in the soil of agriculture mixed cropping (1.18 mg kg-1), which was on par with agriculture mono-cropping (1.31 mg kg-1). Similarly, the available zinc content decreased with depth. A significantly higher amount of DTPA-Zn was found in surface soil (2.69 mg kg-1) as compared to the subsurface (30–45 cm) soil layer (1.47 mg kg-1). The interaction effect of depth x land use was significant, thus the higher value of 3.52 mg kg-1 DTPA-Mn was recorded in surface soils of pure coffee plantations and the lower value was recorded in subsurface (30-45 cm) soils under agriculture mixed cropping (0.85 mg kg-1). These results were in line with Shivakumar,(2013) who reported that the highest DTPA zinc availability was found in the coffee system (1.81mg kg-1) followed by paddy (1.51 mg kg-1), areca nut (1.36 mg kg-1) and banana (1.27 mg kg-1). Significantly, a lower amount of available Zn was found in the soil of the grassland system (0.93 mg kg-1).

 3.4.4 DTPA-Cu

 The data on DTPA-Cu content in soils, as influenced by depth and LUSs, are presented in Table 5. The mean DTPA-Cu content of soil varied significantly among LUSs and depths. Soils under pure coffee (1.74 mg kg-1) had the highest mean DTPA-Cu content, which was on par with natural forest (1.72 mg kg-1), coffee multi-storeyed system (1.69 mg kg-1) and forest plantation (1.62 mg kg-1). In soils under agricultural mono-cropping systems, significantly lower levels of DTPA-Cu (1.16 mg kg-1). With increasing depth, the mean DTPA-Cu content reduced significantly, and a higher amount of DTPA-Cu was found in surface soil (1.63 mg kg-1) as compared to the subsurface (30–45 cm) soil layer (1.45 mg kg-1).

DTPA-extractable micronutrients (Fe, Mn, Zn, and Cu) were found to be above critical levels in both surface and subsurface soils under a different of LUSs (Tables 5).  Higher content over the threshold level could be due to low soil pH, as the solubility of certain micronutrients increases as pH falls (Brady and Weil, 1996). Second, increased intensity of OM in forest soils developed through higher litterfall and root biomass that elevated the aeration status of soils, prevent oxidation and precipitation of micro-nutrients in bound forms along with supplementation of chelating agents that enhanced solubility and availability of micronutrients in soils. (Dhaliwal, et al. 2023). Lindsay and Norwell (1978) found that soil reactivity and organic matter concentration were the two most critical factors controlling micronutrient availability. Soil organic carbon showed a positive and significant correlation with all the micronutrient cations (Zn, Fe, Mn, and Cu), being the highest correlation (r = 0.887) with DTPA Cu (Table 7).

The soil DTPA-Fe and Mn content was found to be higher in forest plantations, followed by natural forest and coffee plantations. The increased concentration could be related to increased litter recycling and chelation processes in forest soils (Dhaliwal, et al. 2023:). Coffee soils contained higher levels of DTPA-Fe and Mn among coffee crops (22.74 and 8.29 mg kg-1). Higher litter turnover and the application of external organic inputs such as FYM, leftover pruning materials, coffee processing waste, and compost could be ascribed to this (Shivakumar et al., 2020). Extractable DTPA Cu and Zn followed the same trends as Fe and Mn, but their content was maximum in the coffee multistoried system (Zn) and pure coffee system (Cu). Copper and Zn concentrations were substantially higher in a coffee-based land-use system. Higher Cu and Zn content in coffee plantation systems due to the frequent application of Bordeaux mixture, Zn, and Cu-fungicides, as well as the high organic carbon content. DTPA extractable micronutrient concentration significantly decreases with increases in soil depth from 0–45 cm, which might be due to decreasing organic carbon content and increasing soil pH with increasing soil depth.

The increased rates of agricultural and tillage practices that lead to a faster rate of OM decomposition, loss of soil nutrients, and increased nutritional demands with lower use efficiencies by plants may be associated with a lower content of DTPA-extractable micronutrients under agricultural LUSs (Dhaliwal, et al. 2023). The decreased input of organic carbon to soils and the ongoing removal of nutrients through nutrient agglomeration by crops in cultivated soils may also be responsible for the decline in micronutrients observed with the conversion of natural forest to agricultural land.

Correlation and regression analysis

Irrespective of different LUSs and soil depths, available N, secondary, and micronutrients of soil were found significantly well-correlated (Table 6) with SOC (P=0.01). The available sulfur showed the highest positive correlation with SOC (r = 0.909, P = 0.01) followed by Cu (r=0.887, P= 0.05) available nitrogen (r=0.876, P= 0.01), Mn (r=0.757, P= 0.01), Fe (r= 0.743, P= 0.01) and exchangeable Ca exhibited least (r= 0.296, P= 0.01). In all six LUSs, available N, S (Figure 5c), and micronutrients Zn, Cu (Figure 5a) Fe, Mn (Figure 5b) exhibited a positive linear relationship with SOC. The regression analysis showed a trend of nitrogen (R = 0.87) > Fe (R = 0.84) > Cu (R = 0.0.81) > Mn (R = 0.79) > S (R = 0.0.73) > Zn (R = 0.67) regarding their relationship with SOC under different LUSs. The aggradation of SOC under different LUSs is reflected by increasing the concentration of all soil nutrients.

Conclusion

            LUSs play an important role in the distribution of carbon stock, nitrogen stock, secondary and micronutrients in soils. The study concluded that among different LUSs viz. natural forest, coffee plantations, and agricultural land, the forest land use was the most eco-friendly and sustainable system followed by the coffee land use system. Natural forest conversion to coffee-based LUSs results in losses of 23-32, 30-37, and 30-39 percent carbon stock at soil depths of 0-15, 15-30 and 30-45 cm, respectively. Similarly, 58-72, 60-74, and 70-71 percent loss of C stock from 0-15, 15-30, and 30-45 cm soil depths, respectively after conversion from natural forest to agriculture LUSs. The natural forests had the highest levels of soil CEC, PBS, secondary nutrients, and micronutrients followed by coffee land use and the least in agricultural LUSs. The correlation study showed that SOC had a favorable impact on the availability of secondary nutrients and DTPA-extractable micronutrients in soils. The present study’s findings had a key impact in soil and crop productivity, which is useful for farmers in implementing appropriate land use selection and nutrient management techniques.

Funding statement: This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Data availability statement: Data will be made available on request.

Declaration of competing interest: The authors declare no conflict of interest.

Reference

  1. Achard, F., Beuchle, R., Mayaux, P., Stibig, H. J., Bodart, C., Brink, A., Carboni, S., Desclée, B., Donnay, F., Eva, H.D. and Lupi, A., 2014. Determination of tropical deforestation rates and related carbon losses from 1990 to 2010. Global change biology, 20(8), 2540-2554.
  2. Mandal, A., Toor, A.S. and Dhaliwal, S.S., 2018. Effect of land-uses on physico-chemical properties and nutrient status of surface (0–15 cm) and sub-surface (15–30 cm) layers in soils of South-Western Punjab, India. Int. J. Curr. Microbiol. Appl. Sci, 7: 2659-2671.
  3. Anand, M.R., Shiva Kumar, H.D., Kommireddy, P. and Kalyana Murthy, K.N., 2019. Secondary and Micronutrient Management Practices in Organic Farming-An Overview. Current Agriculture Research Journal, 7(1).
  4. Anonymous. 2005. Watershed atlas of Chikmagalore (sub, mini and micro watersheds) report sponsored by Indian Space Research Organization and Government of Karnataka, Delineated and published by Karnataka State Remote Sensing Application Center (KSRSAC), Bangalore. (34-38).
  5. Arevalo, C.B., Bhatti, J.S., Chang, S.X. and Sidders, D., 2009. Ecosystem carbon stocks and distribution under different land-uses in north central Alberta, Canada. Forest Ecology and Management, 257(8): 1776-1785.
  6. Barker, A.V. and Pilbeam, D.J. eds., 2015. Handbook of plant nutrition. CRC press.
  7. Black, C.A., 1965. Method of soil analysis part 2. Chemical and microbiological properties, 9, pp.1387-1388.
  8. Bouajila, K., Hechmi, S., Mechri, M., Jeddi, F.B. and Jedidi, N., 2023. Short-term effects of Sulla residues and farmyard manure amendments on soil properties: cation exchange capacity (CEC), base cations (BC), and percentage base saturation (PBS). Arabian Journal of Geosciences, 16(7): 410.
  9. Brady, N.C., Weil, R.R. and Weil, R.R., 2008. The nature and properties of soils (Vol. 13, pp. 662-710). Upper Saddle River, NJ: Prentice Hall.
  10. Chandel, S., Hadda, M.S. and Mahal, A.K., 2018. Soil quality assessment through minimum data set under different land uses of submontane Punjab. Communications in Soil Science and Plant Analysis, 49(6): 658-674.
  11. Chemeda, M., Kibret, K. and Fite, T., 2017. Influence of different land use types and soil depths on selected soil properties related to soil fertility in Warandhab area, Horo Guduru Wallaga Zone, Oromiya, Ethiopia. International Journal of Environmental Sciences and Natural Resources, 4(2): 555634.
  12. Chen, B., Han, M.Y., Peng, K., Zhou, S.L., Shao, L., Wu, X.F., Wei, W.D., Liu, S.Y., Li, Z., Li, J.S. and Chen, G.Q., 2018. Global land-water nexus: Agricultural land and freshwater use embodied in worldwide supply chains. Science of the Total Environment, 613: 931-943.
  13. Chimdi, A. and Gurmessa, E., 2023. Effects of land use on selected soil properties and carbon sequestration in the Gurra watershed: The case of Chalia District, West Shoa Zone, Oromia, Ethiopia. South African Journal of Botany, 156: 21-28.
  14. Dariah, A., F. Agus, and Maswar. 2008. Soil quality of the land under coffee-based farming system (Case study at Sumberjaya, West Lampung). Indonesian Soil Climate Journal1: 48-57.
  15. Das, S., Liptzin, D. and Maharjan, B., 2023. Long-term manure application improves soil health and stabilizes carbon in continuous maize production system. Geoderma, 430, :116-338.
  16. DeFries, R.S., Foley, J.A. and Asner, G.P., 2004. Land‐use choices: Balancing human needs and ecosystem function. Frontiers in Ecology and the Environment, 2(5): 249-257.
  17. Dessalegn, D., Beyene, S., Ram, N., Walley, F. and Gala, T.S., 2014. Effects of topography and land use on soil characteristics along the toposequence of Ele watershed in southern Ethiopia. Catena, 115: 47-54.
  18. Dhaliwal, S.S., Sharma, V., Shukla, A.K., Kaur, J., Gupta, R.K., Verma, V., Kaur, M., Singh, V.K. and Singh, P., 2023. Interactive effect of land use systems on depth-wise soil properties and micronutrients minerals in North-Western, India. Heliyon, 9(2).
  19. Doran, J.W., 2002. Soil health and global sustainability: translating science into practice. Agriculture, ecosystems & environment, 88(2): 119-127.
  20. Duguma, L.A., Hager, H.E. and Sieghardt, M.O., 2010. Effects of land use types on soil chemicalproperties in smallholder farmers of central highland Ethiopia. Ekológia (Bratislava. Ekologia, 29(1): 1-14.
  21. Forests, F.A.O., 2016. agriculture: Land-use challenges and opportunities. State of the World’s Forests, pp.1545-2050.
  22. Geetha, K., 2005. Quantification of litter production, nutrient content in litter and in soil under different vegetation in both high rainfall area and low rain fall area of north Karnataka M.Sc (Agri). Thesis, submitted to Univ. Agric. Sci., Dharwad.
  23. Gomez, K. A. and A. A. Gomez, 1984, Statistical procedures for agric. Res. 2nd Ed. Jhonwilley and sons, New York.
  24. Han, F.X., Kingery, W.L., Hargreaves, J.E. and Walker, T.W., 2007. Effects of land uses on solid-phase distribution of micronutrients in selected vertisols of the Mississippi River Delta. Geoderma, 142(1-2): 96-103.
  25. Heluf, G. and N. Wakene, 2006. Impact of land use and management practices on chemical properties of some soils of Bako area, western Ethiopia. Ethiopian Journal of Natural Resources 8(2): 177-197.
  26. Herojit Singh AS, Singh RK, Singh LN, Singh NG, Chongtham N, Singh AK. 2007, Sulphur status and forms in acid soils of Manipur. Indian Journal of Agricultural Research. ;41(3):205-9.
  27. Ipcc, C. C. 2007. The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment report of the Intergovernmental Panel on Climate Change, vol. 996, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 113–119.
  28. Islam, K.R. and Weil, R.R., 2000. Land use effects on soil quality in a tropical forest ecosystem of Bangladesh. Agriculture, Ecosystems & Environment, 79(1), 9-16.
  29. Jackson, M.L., 1958. Soil chemical analysis prentice Hall. Inc., Englewood Cliffs, NJ, 498, pp.183-204.
  30. Karthika, K. S., K. A. Kumar, K. M. Nair, M. V. D’souza, J. S. Nagaraj, M. Lalitha, A. Koyal, M. Jayaramiah, S. Parvathy, K. Sujatha, and V. Ramamurthy. 2020. Land quality and soil organic carbon stocks consequent to land use change from natural forest to coffee plantation in a hot humid tropical ecosystem of Western Ghats, South India: Is it restorative enough? Authorea Preprints.
  31. Kaur, T., Sehgal, S.K., Singh, S., Sharma, S., Dhaliwal, S.S. and Sharma, V., 2021. Assessment of seasonal variability in soil nutrients and its impact on soil quality under different land use systems of lower Shiwalik foothills of Himalaya, India. Sustainability, 13(3), p.1398.
  32. Kumar Behera, S., Singh, D. and Swaroop Dwivedi, B., 2009. Changes in fractions of iron, manganese, copper, and zinc in soil under continuous cropping for more than three decades. Communications in soil science and plant analysis, 40(9-10): 1380-1407.
  33. Lambin, E.F., Turner, B.L., Geist, H.J., Agbola, S.B., Angelsen, A., Bruce, J.W., Coomes, O.T., Dirzo, R., Fischer, G., Folke, C. and George, P., 2001. The causes of land-use and land-cover change: moving beyond the myths. Global environmental change, 11(4), pp.261-269.
  34. Lima, A.F.L., Campos, M.C.C., Martins, T.S., Silva, G.A., Brito, W.B.M., Dos Santos, L.A.C., de Oliveira, I.A. and da Cunha, J.M., 2022. Soil chemical attributes in areas under conversion from forest to pasture in southern Brazilian Amazon. Scientific Reports, 12(1), p.22555.
  35. Lindsay, W. L. 1979. Chemical equilibria in soils. John Wiley and Sons Ltd..
  36. Lindsay, W. L. and W. A. Norwell. 1978. Development of a DTPA soil test for Zn, Fe, Mn and Cu. Soil Science Society of America Journal 42: 421-428.
  37. Maini, A., Sharma, V. and Sharma, S., 2020. Assessment of soil carbon and biochemical indicators of soil quality under rainfed land use systems in North Eastern region of Punjab, India. Carbon Management, 11(2): 169-182.
  38. Malo, D.D., Schumacher, T.E. and Doolittle, J.J., 2005. Long-term cultivation impacts on selected soil properties in the northern Great Plains. Soil and tillage research, 81(2): 277-291.
  39. Mandal, A. and Dhaliwal, S.S., 2018, November. Distribution of micronutrients in soils under different land-uses in south-western plains of Punjab. In 83rd Annual Convention of ISSS (pp. 27-30).
  40. Maqbool, M., Akhtar, F., Najar, G.R., Peer, F.A., Baba, Z.A., Dar, K.R. and Ramzan, S., 2018. Altitudinal and depth-wise micro nutrient indexing of high density apple orchards in North Kashmir, India. Journal of Pharmacognosy and Phytochemistry, 7(3): 1905-1911.
  41. Mundoli, S., Joseph, G. and Setty, S., 2016. “Shifting agriculture”: the changing dynamics of Adivasi farming in the forest-fringes of a tiger reserve in south India. Agroecology and Sustainable Food Systems, 40(8), pp.759-782.
  42. Murage, E.W., Voroney, P.R., Kay, B.D., Deen, B. and Beyaert, R.P., 2007. Dynamics and turnover of soil organic matter as affected by tillage. Soil Science Society of America Journal, 71(4), pp.1363-1370.
  43. Nigussie, A. and Kissi, E., 2012. The contribution of coffee agroecosystem to soil fertility in Southwestern Ethiopia. African journal of agricultural research, 7(1), pp.74-81.
  44. Okolo, C.C., Gebresamuel, G., Retta, A.N., Zenebe, A. and Haile, M., 2019. Advances in quantifying soil organic carbon under different land uses in Ethiopia: a review and synthesis. Bulletin of the National Research Centre, 43, pp.1-24.
  45. Page, A.L., Miller, R.H. and Keeny, D.R., 1982. Methods of Soil Analysis, Part II. 2nd (ed.) Agronomy Monogr. ASA and SSSA, Medison, WI.
  46. Paul, O.O., Sekhon, B.S. and Sharma, S., 2019. Spatial variability and simulation of soil organic carbon under different land use systems: geostatistical approach. Agroforestry Systems, 93, pp.1389-1398.
  47. Paz-Kagan, T., Shachak, M., Zaady, E. and Karnieli, A., 2014. A spectral soil quality index (SSQI) for characterizing soil function in areas of changed land use. Geoderma, 230, pp.171-184.
  48. Pingali, P., 2023. Are the lessons from the green revolution relevant for agricultural growth and food security in the twenty-first century. Agricultural Development in Asia and Africa; Springer: Singapore, pp.21-32.
  49. Prabhudev Dhumogod, 2011. Studies on organic carbon distribution pattern under selected cropping systems in hill zone acid soils of Karnataka. M.Sc (Agri). Thesis, submitted Univ. Agric. Sci., Bangalore.
  50. Sainepo, B.M., Gachene, C.K. and Karuma, A., 2018. Assessment of soil organic carbon fractions and carbon management index under different land use types in Olesharo Catchment, Narok County, Kenya. Carbon balance and management, 13(1), pp.1-9.
  51. Salve, A., Bhardwaj, D.R. and Tahkur, C.L., 2018. Soil Nutrient study in different agroforestry systems in north western Himalayas. Bulletin of Environment, Pharmacology and Life Sciences, 7(2), pp.63-72.
  52. Shivakumar, K.M., Prakash, S.S., Nagaraja, M.S., Kumar, C.V. and Dhumgond, P., 2020. Effect of Different Land Use Systems on Major Nutrient Status in Soils of Westernghat-Chikamagalur, Karnataka, India. International Journal of Current Microbiology and Applied Sciences, 9(11), pp.3502-3510.
  53. Shivakumar, K.M. 2013. Studies on carbon and nutrient dynamics in different land use systems of a watershed in the Western Ghat. Ph.D. Thesis, Univ. Agric. Sci., Bangalore.
  54. Shukla, A.K., Tiwari, P.K. and Prakash, C., 2014. Micronutrients deficiencies vis-a-vis food and nutritional security of India. Indian J. Fert, 10(12), pp.94-112.
  55. Shukla, A. K. and S. K. Behera. 2014. Total and extractable manganese and iron in some cultivated acid soils of India: status, distribution and relationship with some soil properties, Pedosphere 24:196–208.
  56. Si, B.C., 2008. Spatial scaling analyses of soil physical properties: A review of spectral and wavelet methods. Vadose Zone Journal, 7(2), pp.547-562.
  57. Singh, M.V., 2008. Micronutrients deficiency in Indian soils and crops. Micronutrient deficiencies in global crop production, pp.93-125.
  58. Steel, R.G.D. and Torrie, J.H., 1960. Principles and procedures of statistics. Principles and procedures of statistics.
  59. Stępień, A., Wojtkowiak, K., Pietrusewicz, M., Skłodowski, M. and Pietrzak-Fiećko, R., 2016. The yield and grain quality of winter rye (Secale cereale L.) under the conditions of foliar fertilization with micronutrients (Cu, Zn and Mn). Polish Journal of Natural Sciences, 31(1), pp.33-46.
  60. Subbiah, B.V. and Asija, G.L., 1956. A rapid procedure for the estimation of available nitrogen in soils. Current science, 25(8), pp.259-260.
  61. Takele, L., Chimdi, A. and Abebaw, A., 2014. Dynamics of soil fertility as influenced by different land use systems and soil depth in West Showa Zone, Gindeberet District, Ethiopia. Agriculture, Forestry and Fisheries, 3(6), pp.489-494.
  62. Toru, T. and Kibret, K., 2019. Carbon stock under major land use/land cover types of Hades sub-watershed, eastern Ethiopia. Carbon balance and management, 14, pp.1-14.
  63. Trendov, M., Varas, S. and Zeng, M., 2019. Digital technologies in agriculture and rural areas: status report. Digital technologies in agriculture and rural areas: status report.
  64. Tufa, M., Melese, A. and W. Tena, 2019. Effects of land use types on selected soil physical and chemical properties: the case of Kuyu District, Ethiopia. Eurasian Journal of Soil Science 8(2): 94-109.
  65. Turner, B. L. E. F. Lambin, and Reenberg. 2007. The emergence of land change science for global environmental change and sustainability, Proceedings of the National Academy of Science USA 104: 20666–20671, https://doi.org/10.1073/pnas.0704119104.
  66. Ufot, U.O., Iren, O.B. and Chikere Njoku, C.U., 2016. Effects of land use on soil physical and chemical properties in Akokwa area of Imo State, Nigeria. International Journal of Life Sciences Scientific Research, 2(3), pp.273-278.
  67. Van der Werf, G.R., Morton, D.C., DeFries, R.S., Olivier, J.G., Kasibhatla, P.S., Jackson, R.B., Collatz, G.J. and Randerson, J.T., 2009. CO2 emissions from forest loss. Nature Geoscience, 2(11), pp.737-738.
  68. Walkley, A. and Black, I.A., 1934. An examination of the Degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method. Soil science, 37(1), pp.29-38.
  69. Walkley, A. and Black, I.A., 1934. An examination of the Degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method. Soil science, 37(1), pp.29-38.
  70. World Bank. South Asia Regional Office. Agriculture and Rural Development Sector Unit, 2006. India: Unlocking Opportunities for Forest-Dependent People. Oxford University Press, USA.
  71. Yao, M.K., Angui, P.K., Konaté, S., Tondoh, J.E., Tano, Y., Abbadie, L. and Benest, D., 2010. Effects of land use types on soil organic carbon and nitrogen dynamics in Mid-West Cote d’Ivoire. European Journal of Scientific Research, 40(2), pp.211-222.
  72. Yoseph, T. D., P. Witoon, B. Amila, Y. Birru, W. Tesfaye, G. Hans, and L.G. Douglas. 2017. Changes in land use alter soil quality and aggregate stability in the highlands of northern Ethiopia. Scientific Reports 7: 13602.

Similar Posts

Leave a Reply

Your email address will not be published.