Abstract

Different land use systems can have different soil properties. It is important to study the soil properties for wise use and sustainable management of land resources. This article reports the findings of a research conducted in Makwanpur district of Nepal, to determine and compare the selected physicochemical and hydrological properties of soil in forest, rainfed agriculture land, and grassland. These forest, agriculture land, and grassland represent tropical land use systems. Soil samples were collected from 0 to 30 cm depths of soil profile from nine randomly distributed pits dug in forest, grassland, and rainfed agriculture land in 2019. Soil samples were analyzed in the laboratory to determine the soil properties using standard methods. Bulk density, porosity, moisture content and infiltration, pH, total nitrogen, available potassium, and available phosphorus were quantified from the soils samples. It was found that the highest BD was found in the grassland (1.29 g/cm³) followed by the forest (1.23 g/cm³) and rainfed agriculture land (1.18 g/cm³). The highest porosity was found in rainfed agriculture land (55.50%) followed by the forest (53.74%) and grassland (51.63%). The highest MC was found in the grassland (26.94%) followed by the forest (10.17%) and rainfed agriculture land (9.92%). The mean cumulative infiltration amount was highest in the rainfed agriculture land (39.27 cm) followed by the forest (33.47 cm) and grassland (8.4 cm). The highest soil pH was found in the grassland (7.91) and the lowest pH (5.70) in the forest. The highest level of total nitrogen was found in rainfed agriculture land (0.121%), followed by the forest (0.106%) and grassland (0.096%). The highest level of available phosphorous was found in rainfed agriculture land (84.94 ppm), followed by the forest (67.76 ppm) and grassland (6.69). The highest level of available potassium was found in rainfed agriculture land (154.24 ppm), followed by the forest (84.49 ppm) and the grassland (44.71 ppm). Bulk density, porosity, and total nitrogen were not found to be significantly different and other soil properties were found to be significantly different between different land use systems. The contribution of farmers in maintaining soil properties on the farmlands is clearly reflected in the results, so their knowledge on soil management needs to be explored and adapted for wise use and sustainable management of other land use systems.

1. Introduction

Soil is the substrate of all living organisms in the terrestrial ecosystems, which contains nutrient reserves and supports many biological processes in vegetation development. Forest biodiversity directly depends on abiotic factors like forest soil and microclimate variables, which have vital roles to sustain forest biodiversity [1]. Land use involves the management and modification of natural environment or wilderness into built environment such as settlement and seminatural habitats such as arable fields, pastures, and managed woods.

Different land use systems exist in the world with variation in soil physical, chemical, and biological properties of soil. Soil properties vary due to different land use and management practices [2]. Soil properties also differ according to the land management practices, so even the similar land use system can have different soil properties. Land use and soil management practices influence the soil nutrients and related soil processes [3]. As a result, they can modify the processes of transport and redistribution of nutrients in the soils [4].

Land degradation resulting from land use change has been a major global challenge since the 20th century and will be focused as a major topic on the international agenda in the 21st century [5]. National Land Use Policy of Nepal has categorized lands of Nepal into seven classes: forestry, grassland, agricultural land, residential, commercial, archeological, and water bodies [6]. Rapid change in land use and land cover in Nepal has been felt due to the high population growth. Nepal is facing a severe problem of soil quality deterioration and of subsequent low crop yields because of different human activities and land use changes which cause lowering of organic matter, nutrients, and their hydrological parameters [7]. Dry soft rocks are the major sources of soils in this area. Variegated mudstone, siltstone, and shale with smaller amounts of fine-grained sandstone are found in this area [8]. The area is characterized by ancient river terraces, depositional basin, and moderately sloping terrain [9].

The evaluation of soil properties under the different land use systems is important to understand the role of different land uses and their effects on soil differentiation [10]. Soil management systems and changes in land use strongly influence soil physical and chemical properties [11]. The variation of soil properties should be monitored and quantified to understand the effects of land use and management systems on soils [12] so that this knowledge can be used in soil management. Though knowledge on variation of soil properties is essential in devising proper land use planning, previous researches devoted to this topic are absent in Hetauda Submetropolitan area. The lack of such researches could be due to low priority of the developing nations in allocating research funds across wider geographic areas over time. In this context, this study was carried out in topical land use systems to determine and compare the selected soil chemical and hydrological properties in forest soil, agriculture soil, and grassland soil in tropical region of Nepal.

2. Materials and Methods

2.1. Research Sites

The research was conducted in three land use systems: forest, grassland, and rainfed agriculture in Makawanpur district, a tropical region of Nepal (27°25′ N latitude and 85°02′ E longitude; 300–390 m above the sea level) in 2019. The soils of the Makawanpur district are Fluvisols, Cambisols, and Regosols [9, 13].

A tropical forest in ward number 10 of Hetauda Sub-Metropolitan with total area 94 hectares was selected for data collection (Figure 1).

The forest represented a mature tropical forest with Shorea robusta as dominant species. The associated species in the forest are Dalbergia sissoo, Acacia catechu, Trewia nudiflora, Bombax ceiba, etc. (own observation). According to Rautiainen [14], the estimated age of the stand with average diameter above 48 cm is 120 years. Hence, the study assumed the age of the stand at around 120 years observing the maturity level of the dominant Shorea robusta (own observation), though scientific assessment of stand age, is yet to be conducted.

A tropical grassland adjacent to the forest was selected for data collection, with a total area of around 3 hectares (Figure 1).

Small grasses like Cynodon dactylon, Artemisia vulgaris, and Imperata cylindrica (own observation). were major species on the grassland, with usual cattle grazing.

Rainfed agricultural land nearby the forest and grassland was selected for data collection, with a total area around 4 hectares (Figure 1).

2.2. Soil Sampling

A stratified random sampling design was used for collecting soil samples in the three representatives land use systems purposively selected: Forest, grassland and rainfed agriculture land. Soil samples, each around half kilogram, were collected from 0 to 30 cm depths of soil profile from nine randomly distributed pits dug in forest, grassland, and rainfed agriculture land. One additional soil sample (undisturbed soil sample) from each pit was collected using the soil core sampler to determine the bulk density, porosity, and moisture content of the soil.

Soil formation is predominantly a vertical process related to climate, parent material, relief, organisms, time, and spatial position [15]. The inputs from weathering and humus are strongest at the soil surface and decrease with depth [16]. The organic matter in soil tends to concentrate in the upper horizons of the soil because of high input from the vegetation on soil surface. According to IPCC [17], roughly half of the soil organic carbon of the top 100 cm of mineral soil is held in the upper 30 cm layer. It is good practice to include a minimum depth of at least 30 cm in nutrient stock assessment. The carbon held in the upper profile is often the most chemically decomposable and the most directly exposed to natural and anthropogenic disturbances. So, we focused our soil sampling upto 30 cm soil depth.

The collected soil samples were packed, labelled, air-dried, and taken to the Soil Test Laboratory of a Non-Governmental Organization (NGO)—Local Initiatives for Biodiversity, Research and Development (LIBIRD), Kaski, for analysis. The overall field measurement methods were based on guidelines of Asia Network for Sustainable Agriculture and Bio resource ANSAB [18].

2.3. Infiltration Test

The infiltration test was carried out in the forest, grassland, and rainfed agriculture land using double ringed infiltrometer [19].

2.4. Soil Analysis

The percentage of moisture content (MC) in the soil was obtained bywhere M1 = fresh weight of soil.

Bulk density (BD) of the soil was determined using the soil core sampler of volume 245.22 cm3 (V) [20].

Porosity (P) of the soil was determined using the equation taking where the particle density (PD) equals 2.66 g/cm3:

The infiltration rate (IR) was obtained bywhere t is the time taken to drop the water from initial level to the final level in the inner ring in minutes. Cumulative infiltration was obtained by adding each infiltration amount in each time recorded.

Soil pH was determined using the distilled water method as indicated by Davis and Freitas [21].

Nitrogen was determined by using the semi-Kjeldhal Dri-Block method [22].

Available phosphorus was determined by using the Bray-p1 method [23].

Potassium in soil was determined by the flame photometer with 1M ammonium acetate extracting solution [24].

2.5. Statistical Analysis

Statistical Package for Social Science (SPSS) software (IBM SPSS Statistics, IBM Corporation, Version 20) was used for statistical analysis of the data. Analysis of variance (ANOVA) and least significance difference (LSD) tests were conducted to test the significance difference of the variables.

3. Results and Discussion

3.1. Comparison of Physical Properties on Land Use Systems
3.1.1. Bulk Density (BD)

The highest BD was found in the grassland (1.29 g/cm³) followed by the forest (1.23 g/cm³) and rainfed agriculture land (1.18 g/cm³), similar to the results of Muche et al. [25], who reported that the highest average value of BD (1.2 gm/cm3) was obtained under the grazing land. The ANOVA test shows that there was no significant difference of BD () between different land use systems at 5% level of significance (Table 1). The lowest BD of the rainfed agriculture land might be due to the loosening of soil due to frequent nonmechanized tillage operations. The highest BD of the grassland might be due to the regular grazing and movement of the humans on the ground which causes soil compaction [26]. Well-aggregated loose and porous soil has high organic matter and lower bulk density while poorly aggregated soil has low organic matter content and high bulk density that make total pore spaces become greater [27].

3.1.2. Soil Porosity (SP)

The highest porosity was found in rainfed agriculture land (55.50%) followed by the forest (53.74%) and grassland (51.63%). The ANOVA test shows that there was no significant difference of porosity () between different land use systems at 5% level of significance (Table 1). Porosity was highest in the rainfed agriculture land which may be due to the frequent addition of manure and nonmechanized tillage. The bulk density might have played a role in the percentage of the porosity which has reported that the decrease in bulk density with organic amendments and nutrients is directly related to the increased porosity, which is related to the improved soil aggregation [28]. Cattle grazing and human movement might have caused substantial reduction of the soil pore space causing lowest porosity in the grassland [26].

3.1.3. Moisture Content (MC)

The highest MC was found in the grassland (26.94%) followed by the forest (10.17%) and rainfed agriculture land (9.92%). The ANOVA test shows that there was significant difference of MC () between different land use systems at 5% level of significance (Table 1). The LSD test shows that there was significant difference of MC between grassland and forest and between grassland and rainfed agriculture land but no significant difference between forest and rainfed agriculture land at 5% level of significance (Table 2). A similar result was found by Tuffour et al. [26] in which the moisture content of the soil was lower in the ungrazed field (17%) in relation to the grazed field (24%). It was stated that this result could be due to lower porosity resulting from the higher bulk density in the grazed field. Reduction of the pores in the grassland may have caused a restriction in the transmission of water with the soil becoming frequently too wet [29] resulting in an increased amount of water available on the soil surface.

3.1.4. Infiltration Parameter

The mean cumulative infiltration amount was highest in the rainfed agriculture land (39.27 cm) followed by the forest (33.47 cm) and grassland (8.4 cm) (Table 1).

For the first minute, the mean infiltration rate of the rainfed agriculture land, forest, and grassland was 1.9 cm/min, 1.53 cm/min, and 0.567 cm/min, respectively.

The ANOVA test shows that there was significant difference of cumulative infiltration amount () between different land use systems at 10% level of significance (Table 1). The LSD test at 10% level of significance revealed that there is significance difference in cumulative infiltration amount in between grassland and rainfed agriculture land () and forest and grassland but no significant difference between forest and rainfed agriculture land (Table 2).

Higher infiltration was found in the rainfed agriculture land which might be due to the higher porosity (55.40%). The lower porosity (51.63%) in the grassland caused decrease in infiltration due to few pores and the compaction of soil. Also, the moisture content of the land might also be the reason for the difference in the infiltration amount. Similar study was done by the Tuffour et al. [26], and they found a lower infiltration amount in the grazed compacted field with respect to the ungrazed field where the terminal infiltration amount of spot 1 on the ungrazed field was 1320 mm compared to 510 mm of the grazed field.

3.2. Comparison of Chemical Properties on Land Use Systems
3.2.1. Soil pH

The highest soil pH was found in the grassland (7.91) and the lowest pH (5.70) in the forest. The ANOVA test shows that there was significant difference of soil pH () between different land use systems at 5% level of significance (Table 1). The LSD comparison at 5% level of significance between different land use systems for the soil pH revealed that there is a significance difference between all land use systems (Table 2).

The lowest pH value of the forest soil reveals that the forest soil is acidic which might be due to the presence of the dominant Shorea robusta tree in the forest [30]. The highest pH value of the grassland reveals that the grassland soil is basic which might be due to the accumulation of the salt concentration because of the nearer distance from the river. The pH value of the rainfed agriculture land is around neutral (6.55) which might be due to the ongoing soil management practices, for example, use of compost and proper irrigation. Similar results were found by the Islam and Weil [31], where forests and reforested soils contributed significantly to lower pH values than cultivated and grassland soils in Bangladesh.

3.2.2. Total Nitrogen

The highest level of total nitrogen was found in rainfed agriculture land (0.121%), followed by the forest (0.106%) and grassland (0.096%). The ANOVA test shows that there was no significant difference of total nitrogen () between different land use systems at 5% level of significance (Table 1). The highest level of total nitrogen in rainfed agriculture land might be due to the cultivation of the leguminous crops, mulching of such crops, and the application of the compost as well as inorganic fertilizers. The lowest level of total nitrogen in the grassland might be due to the sparse vegetation due to human and livestock disturbance.

3.2.3. Available Phosphorous

The highest level of available phosphorous was found in rainfed agriculture land (84.94 ppm), followed by the forest (67.76 ppm) and grassland (6.69). The ANOVA test shows that there was significant difference of available phosphorus () between different land use systems at 5% level of significance (Table 1). The LSD comparison at 5% level of significance between different land use systems for the available phosphorus revealed that there is significant difference between all land use systems (Table 2). Highest amount of available phosphorous in the rainfed agriculture land might be due to the ongoing soil management practices, for example, application of manure and fertilizers. The lowest available phosphorus content in the forest soil might be due to forest vegetation with availability of different plant species which absorbs a larger amount of available phosphorus causing to deplete it in the soil [32]. Leaching might be the cause for lowest level of available phosphorus in the soil of the grassland. Kharal et al. [33] found a similar result where available phosphorus was the highest in the vegetable farm (41.07 mg kg−1) and was the lowest in grazing land (2.89 mg kg−1). Bizuhoraho et al. [34] and Ivezic et al. [35] also found the higher available phosphorus in cultivated land than in the forest.

3.2.4. Available Potassium (ppm)

The highest level of available potassium was found in rainfed agriculture land (154.24 ppm), followed by the forest (84.49 ppm) and the grassland (44.71 ppm). The ANOVA test shows that there was significant difference of available potassium () between different land use systems at 5% level of significance (Table 1). The LSD comparison at 5% level of significance between different land use systems for the available potassium revealed that there is significance difference between grassland and rainfed agriculture land (Table 2). The highest level of available potassium in the rainfed agriculture land might be due to the litter mixing with soil and undergoing recycling to release available potassium and applying ashes of the burned crop residues [36] and application of the cow dung, poultry manure and compost manure in the field. The lowest amount of available potassium in the grassland might be due to leaching of the available potassium from soil surface and low external input [36] and might be due to the disturbance to the soil from trampling by livestock, human movement, and soil compaction. Kharal et al. [33] found similar result where available potassium was the highest in the vegetable farm (130.2 mg kg−1) and lowest in grazing land (36.8 mg kg−1). Ivezic et al. [35] also found the higher available potassium in cultivated land than in the forest.

4. Conclusion

This article indicates the first reporting of hydrological properties in soils of tropical forest, rainfed agriculture land, and grassland in Nepal. Results of this research show evidence of variation of soil physicochemical and hydrological properties in different land use systems of tropical region of Nepal. The soils with external disturbances showed low levels of important physiochemical and hydrological properties of the soils. The control of rampant grazing on the grassland will improve the hydrological properties of the grassland such as porosity and infiltration. The rainfed agriculture land shows superiority in soil properties, so existing farming practices should be promoted, and associated indigenous knowledge should be documented and disseminated. How soils have changed over time under specific management scenarios needs to be studied in-depth.

Data Availability

The data supporting the conclusions of this study are available from the first author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

The data collection for this research was supported by the Faculty of Forestry, Agriculture and Forestry University, Nepal. The authors are highly grateful to the Faculty of Forestry, Agriculture and Forestry University, for providing lab equipment and library facility.