Abstract

Understanding the role of soil carbon (C) dynamics and quantitative changes as affected by various land use patterns is very critical given the significance of carbon sequestration. In this context, the current study was conducted in the Lal Bakaiya watershed in Makawanpur District, Nepal, to assess the variation of soil organic carbon (SOC) and nitrogen (N) stocks in three different land use types, namely, natural forest, grassland, and cultivated land. Incremental soil depths method (i.e., 0–15 cm, 16–30 cm, and 31–45 cm) was applied to collect soil samples in bulk from each of the land use under the study to estimate SOC and N stocks in laboratory. A total of 90 soil samples were collected from three soil layers down the soil profile up to 45 cm for each land uses. The results show that both SOC and N contents decreased with soil depths; however, substantial amount of SOC and N stocks were reported in lower soil depths under land use with natural forest. Both SOC and N contents were found relatively higher at 0–15 cm depth in natural forest soil (1.40 ± 0.20% and 0.26 ± 0.04%) than those in grassland and cultivated land, respectively. The mean total SOC stock and N stock ranged from 46.3 ± 4.24 t ha−1 and 7.11 ± 1.86 t ha−1 in cultivated land to 62.05 ± 9.17 t ha−1 and 11.40 ± 1.92 t ha−1 in the land use with natural forest, respectively. Furthermore, the mean total carbon and nitrogen ratio (C/N ratio) of the soil was found to be higher in cultivated land (7.07 ± 1.93) than that in natural forest (5.75 ± 1.47) and grassland (5.62 ± 1.49), respectively. Two-way analysis of variance results showed that both land use type and soil depth have significantly () affected the SOC and N stocks in the study. From the results, it is suggested that well-managed land use can contribute significantly in offsetting global carbon emission.

1. Introduction

Soil comprises the largest active terrestrial carbon (C) pool globally [1, 2]. About 1200–1800 gigatons of C are estimated to be stored in soil worldwide. Soil C pool stores more than three times C than the amount stored in atmosphere and 3.8 times more than in biotic pool [1, 3]. Therefore, estimating the amount of C in soils can provide a significant opportunity to mitigate global warming [4, 5]. Thus, improving the capture and storage of atmospheric C through improved land use systems can be a good strategy to lower its concentration while also improving the quality of soil.

Soil plays a crucial role in the ecosystem functioning and provisioning of its services [6, 7]. The function and viability of soil depends on the dynamic symmetry among its biophysical and chemical properties [8]. Soil organic carbon (SOC) and nitrogen (N) are fundamental for soil health, sustainable biological productivity, and environmental quality. In addition, the SOC and N contents in soil reveal the long-term equilibrium between additions and losses of organic C from various processes and pathways [9]. Different land use practices often influence the fluxes of SOC and N stocks and have been reported to differ with the change in land use types due to the combined effect of biophysical and chemical processes over time [10, 11]. SOC and N stocks can either increase or decrease depending on a number of factors including climate, soil type, topography, and soil management techniques. Therefore, the land use system is one of the most significant factors that control SOC and N stocks build up because of the presence of SOC and N in the soil. Studies have shown that changes in land use, particularly from natural vegetation to others, can cause SOC losses of 20–50% [1214]. The SOC stock reduced in significant amount after using fire for land preparation and conversion from natural forest to maize field and change from maize field to other crops or grassland in northern Thailand [15, 16]. Similarly, converting natural vegetation cover to large-scale farmland resulted in significant loss of soil organic carbon and total nitrogen in Western Ethiopia [17]. Other studies also reported that a strong decrease of SOC and N after conversion of forest land to other land [18, 19], and this has provided a way to raise concern about the long-term sustainability of those land use types in the tropics [20].

C sequestration is the long-term process of capturing and storing of CO2 from the atmosphere to reverse atmospheric CO2 pollution and to reduce global climate change [21, 22]. With ever-increasing global warming issues, the topic of soil C sequestration is something that we cannot afford to pass over [1, 22]. Globally, soils are viable sinks for atmospheric C and they are an important part of terrestrial C pool [1, 23, 24]. About 1500 pentagrams of C is reserved in soils in the form of organic matter that accounts approximately twice the terrestrial C pool [25], and Nepalese soil holds 308.78 million tons [26]. Conversion of land use particularly of a natural into a managed system can change soil C pools and can exert critical impact on balance between soil properties and the atmosphere [27]. In tropics, about 12 to 20% anthropogenic greenhouse gases (GHGs) is due to land use changes which is reported as the second largest source of GHGs [19, 21]. Also, this scenario is likely to remain also for the future if we do not address the situation in time. In Nepal, total greenhouse gas emissions in 2014 were 44.06 million metric tons of CO2, and the figure has increased to 51.24 million metric tons in 2017 [28]. Estimating SOC and N stock in various land uses has become very important because it will assist policy makers to work out approaches for managing land use systems sustainably as well as preventing loss of SOC and N stocks. In this backdrop, the objective of this study was to evaluate the variability in the SOC and N stocks with respect to various land use types in the northern part of Lal Bakaiya watershed in Bagmati Province, Nepal. Specifically, this study intended to address two concerns: (i) variation of SOC and N contents within different soil layers and (ii) impact of different land use types on SOC and N stocks.

2. Materials and Methods

2.1. Study Site

This study was carried out in three major land use types such as natural forest, grassland, and cultivated land in the northern part of the Lal Bakaiya watershed located between 27°30′00″ N and 85°20′00″ E in Makawanpur District, Nepal. It is located in the central southern part of the country and covers about 868 km2 of the total basin area (Figure 1). Out of the total watershed area, 42% lies in Makawanpur District, 20% in Bara, and 38% in Rautahat District [29]. The study focuses only on Makawanpur District which covers the geography from Terai to Mahabharata ranges. Furthermore, the watershed area can be partitioned into three major geological units: the Terai and Bhabar in the south, Chure/Siwalik in the middle, and the Mahabharata range in the north. The length of the river is 134 km and flows through elevations ranging from 2,000 to 71 m above MSL. The elevation of the study site ranges from 435 to 2000 m above MSL. The mean annual precipitation of the area is 1,434 mm (at Rautahat) to 2,306 mm (at Makawanpur) [29]. The soil texture in this area is sandy loam, and the area is environmentally vulnerable due to frequent landslide, soil erosion, forest encroachment, and deforestation.

For land use land cover classification, Landsat images were extracted for free from Earth Explorer of USGS, and object-based image analysis (OBIA) technique was employed to classify Landsat images to acquire land use map for the particular site [30, 31]. After that, on-field observation was carried out to validate the basic information topography of the study site, following that the boundary of the watershed was identified by using ArcGIS 10.2. Accordingly, three major land use types as natural forest, grassland, and cultivated land were defined for allocation of various samples.(1)Natural forest: It includes land having at least 0.5 ha area and a minimum width/length of 20 m with trees higher than 5 m at maturity and a crown cover of more than 10%, not being utilized principally other than forestry purposes (Figure 2) [32].(2)Grassland: It includes land in which the vegetation is dominated by a nearly continuous cover of grasses. They are commonly used for open grazing, maintained by grazing animals and sometimes for cut and carry system (Figure 3) [32].(3)Cultivated land: It includes land primarily used for agricultural purposes such as production of crops and raring of livestock (Figure 4) [33].

2.2. Sampling Techniques

The field visit and soil sampling carried out starting from 13 March to 11 May 2018 for two months. Stratified random sampling technique was adopted for the study, and land use types were identified as the strata for soil sample collections. Upstream and downstream approach was followed for sampling. Two transect lines were delineated along the slopes; one in the downstream (at 500 m elevation) and the other in the upstream (at 1500 m) part of the watershed crossing natural forest, grassland, and cultivated land. The possible sampling points through transect lines were pointed out in Google earth, and when all likely points were being defined, a random selection was done. A total of 15 pairs of sample sites were identified (15 in each transect line); total 30 sample sites with 90 soil samples were collected for the study. At every sampling site, a pit of 30 cm by 50 cm was dug, and undisturbed soil core samples were collected by using a cylindrical core sampler (5.5 cm diameter and 5 cm height) from the 0–15 cm, 16–30 cm, and 31–45 cm soil depths for the determination of bulk density (BD), SOC content, and N content. The soil samples collected were analyzed at the Provincial Soil Laboratory, Hetauda, Nepal.

2.3. Data Analysis

Soil BD was estimated using the core sampling method [29]. BD was room dried for 7 days, and then it was oven dried at constant temperature of 105°C for 24 h. Oven-dried weight was taken and soil was processed through 2 mm sieve to differentiate stones. Volume of the stone was taken for stone correction by putting it in a cylindrical jar with water. BD was then determined by the following formula:where ρ = bulk density in g cm−3, Md = oven dry weight of soil in g, Vt = volume of the soil in cm3, and volume of the soil = volume of the core − volume of the stone.

The SOC content in the soil samples was estimated by using the Walkley–Black wet oxidation method [34]. Accordingly, soil N content was determined by the Kjeldahl method [31]. Finally, SOC and N stocks for each depth of land uses were determined by the following equation suggested by Pearson et al. [32].where SOC = soil organic carbon stock per unit area in t ha−1, ρ = soil bulk density in g cm−3, D = depth of soil horizon at which the sample was taken in cm, and %C = organic carbon concentration in %.where NS = soil nitrogen stock per unit area in t ha−1, ρ = soil bulk density in g cm−3, D = depth of soil horizon at which the sample was taken in cm, and %C = nitrogen concentration in %.

The collected data were then categorized and summarized by land use types and soil depths. The two-way ANOVA model was applied to test the effect of land use types and soil depth on SOC and soil N. Furthermore, the correlation analysis was applied to determine relationship between the studied variables. Final data were analyzed by using Statistical Package for the Social Sciences (SPSS) software (version 20.0).

3. Results

3.1. Bulk Density

The study found that there was a significant variation in the BD with reference to soil depths in all three land uses. Higher variation was recorded in 0–15 cm depth, whereas BD values in 31–45 cm depths were clustered. The BD value was ranged from 0.96 g cm−3 to 1.54 g cm−3. Significantly, lower BD values were observed in soil under natural forest while higher values were found in land under cultivation over three studied soil depths (Table 1), whereas in the grassland, the BD values were found intermediate between forest and cultivated land. BD values were found to increase with incremental soil depths in all three land use types. Furthermore, the result from two-way ANOVA showed that there are significant () differences in BD by land use types and soil depths (Table 2). Even so, the interaction effect of land use types and soil depths on BD was not found significant ().

3.2. Organic C, Total Nitrogen Content, and C/N

In this study, higher variation of the SOC content was found in 0–15 cm depth than in 16–30 cm and 31–45 cm depths in all land uses. SOC contents were found from 0.57% to 1.70%. Significantly, higher SOC contents were found in natural forest than grassland and cultivated land over the three soil depths (Table 3). In all three land uses, SOC contents were significantly decreased with increase in depths; however, it did not show any consistent trend in the case of grassland (Table 3). In the case of grassland, organic carbon contents reported higher in 16–30 cm depth than in 0–15 cm depth. Accordingly, higher variations of the total nitrogen content were observed in 0–15 cm and 16–30 cm depth than in 31–45 cm in all three land uses. N content values were reported from 0.05% to 0.33%. Significantly, a higher N content was reported in natural forest followed by grassland and cultivated land, respectively. Furthermore, the N content is also decreased with increase in soil depths except the grassland in which it did not indicate any consistent trend down the soil profile (Table 3).

In this study, the C/N ratio did not follow likely trend of SOC and N contents with soil depths in all three land use types except in the case of cultivated land that showed the increasing trend with increasing soil depths (Table 3). In comparison, higher values of the C/N ratio were observed in cultivated land compared to forest and grasslands. The two-way ANOVA showed that both SOC content and N content were significantly () affected by land use types and soil depths (Table 4). Furthermore, both soil properties were significantly () affected by the interaction effect of land use and soil depth (Table 4). In addition, the C/N ratio was also affected by land use and soil depths but not affected by interaction of land use and soil depths (Table 4).

3.3. SOC and Nitrogen Stocks

From the study, significantly higher amount of SOC stock was found in all three depths of natural forest (Table 5). Furthermore, a declining trend in the SOC stock was reported in all three land uses with incremental soil depths. However, the contrasting trend was found in grassland as higher SOC stock was found in 16–30 cm than in 0–15 cm. Similarly, the same kind of pattern was found in the N stock within all three land uses (Table 5). Furthermore, the two-way ANOVA result shows that both SOC and N stocks were significantly () affected by land use types and soil depth and further by the interaction effect of land use types and soil depths (Table 6).

4. Discussion

4.1. Variation in BD

In this study, significantly lower BD values were observed in natural forest compared to other two land uses. Higher organic inputs and lower losses might be reasons for lower BD values in the natural forest. Other similar studies have also highlighted significantly lower BD values in the natural forest soil compared to cultivated land, barren land, and grassland in central midland of Nepal [33, 35, 36]. Accordingly, BD values were found to be increased with the increasing soil depth for all three land uses. The increase in BD with incremental soil depth could be related with the decrease in OM source, root penetration, and the compaction pressure of the overlying soil mass. The findings are in line with those in [37, 38] who reported lower OM due to fewer root networks and less aggregation, and also the compaction caused by overlying soil layers increases the BD of soil. Many researchers in Nepal and around the world also reported a general increase in BD values with increasing soil depths [33, 36, 3941]. The result is further supported by highly significant () and negative correlation (r = −0.59) between BD and soil C contents.

4.2. Variation in SOC, N Contents, and C/N within Different Land Uses

This study found that soil under natural forest had higher organic C and N contents’ level than grassland and cultivated land. Higher amount of organic C content in forest soil could be associated to improved accumulation of higher OM and decomposition of fine roots and also to frequent tillage in cultivated soil. As reported in [33, 36, 42], a higher level of organic C content is observed in forest soil than in agricultural and pastures land in central and western Nepal. Similarly, the authors in [4144] reported that the organic C content in the cultivated soils is less protected than in undisturbed soils due to removal of large quantities of the biomass during clearing and land preparation. Accordingly, the higher N content in the natural forest could be associated with its higher organic C concentration, which is the prime source of soil N [45, 46]. This finding is further assisted by significant () and positive correlation (r = 0.69) between organic C content and N content. Many similar studies carried out in Nepal and elsewhere also reported a higher level of organic C content and N content in forested land than in agricultural and other land uses [38, 47, 48].

The average C/N ratio of soils in forest, grassland, and cultivated land were found to be 5.57 ± 1.33, 5.26 ± 1.49, and 7.06 ± 1.93, respectively. However, the observed values of C/N ratios were out of the normal range of international soil which is about 10 to 12 [49]. As reported in [50], the OM matter with lesser C/N ratio (<10) opens up the N content leading to rapid release of N into the soil and increasing availability of N for the crop. The two-way ANOVA showed that both SOC content and N content were significantly () affected by land use types and soil depths (Table 4). Kafle [51] also reported a significant difference of C/N with incremental soil depths in a community forest of Chitwan District, Nepal, within 1 m soil profile.

4.3. Variation in SOC and N Stocks within Soil Depths of Different Land Uses

Environmental factors such as vegetation cover, climate, and management practices have a significant impact on the soil’s capacity to hold or release C [52]. The result found that both SOC stock and N stock were significantly higher in natural forest in comparison to grassland and cultivated land (Table 5). The result is further supported by a highly significant () and positive correlation (r = 0.60) between SOC stock and N stock. The results show that natural forest land has a higher stock of SOC and soil N than grassland and cultivated land. Higher amount of SOC and N stocks in the natural forest could be associated with a higher organic C content due to the accumulation of aboveground leaves, ground litter, and underground root litter [53]. The findings are in line with those reported in [36, 5457] which reported higher SOC stock than other land use types in Nepal and around the world.

Compared with natural forest, lower SOC and N stock in grassland land is associated with low organic inputs due to removal of aboveground biomass for livestock feed (i.e., cut carry system) in grassland and cultivated land in the watershed area which is further assisted by grazing pressure particularly in grassland [36, 41, 46, 58]. Furthermore, soil disturbance during land preparation and tillage activities exposes OM to decomposition and causes rapid losses in SOC and N stocks in cultivated land due to loss of root biomass and removal of crop residues [42, 5961]. Due to the low input of organic matter from harvested farm residues, cultivated lands often have lower SOC stocks. In addition, the study also observed that farming systems (i.e., agriculture-based livelihood) practiced by local farmers in the study site are exploitative in nature. In this regard, a land use type that can improve SOC and soil N accumulation and protect against the loss of cations through leaching and biological processes should be promoted for better land management and to fight against the undesirable impacts of climate change as well [62, 63]. Thus, appropriate management intervention such as controlled or rotational grazing, optimizing livestock number, and addition of organic inputs should be applied to enhance C stock and N stock in grassland [41, 64]. Similarly, suitable practices such as conservation tillage, zero tillage, terrace farming, and agroforestry practices should be practiced to improve the carbon storage capacity of cultivated land [36, 59, 65, 66]. These findings are in line with the research of Leul et al. [17] on effects of land use dynamics on soil organic carbon and total nitrogen stock.

Depth-wise average SOC stock and N stock results are presented in Table 5. Significantly, a higher amount of SOC stock and N stock was reported in the top soil layer than in subsoil layers within the profile of the same land use. The highest SOC in the 0–15 cm depth was observed in the natural forest land followed by cultivated land and grassland (Table 5). Furthermore, the highest soil N stock in the topsoil was measured in the natural forest land followed by cultivated land and grassland (Table 5). The result demonstrates that both SOC and N stocks decrease with the increasing depths down the soil profile and it is attributed to the less accumulation of OM content at lower depths down the soil profile [41, 54]. Higher amount of SOC and N stocks in the 16–30 cm layer of grassland is attributed to the rapid leaching and greater accumulation of organic inputs from the top soil layer to the subsoil layer. Findings from other studies [33, 36, 6769] also suggest that relatively higher amounts of SOC and N stocks are present in the top soil layer than those of the deeper layers. High OM content, greater root biomass, and higher accumulation of vegetative residues are the major reasons for higher amount of SOC and N stocks in the upper soil layer. The losses from decomposition, leaching, and soil erosion are major reasons for the low SOC stock. Low SOC storage levels in the topsoil may also be a result of poor land management and postharvest grazing [70]. The findings of this study are consistent with those reported in [41,71] which reported that land use types had a significant impact on the SOC and N contents, and soils under natural vegetation had a higher amount of SOC stock compared to shrub land and cultivated soil in the middle mountain region of Nepal. Studies around the world highlights that land use and land use changes are important determinants of landscape C stocks and that careful management of existing natural and managed ecosystem is critical to global C stock [72]. Therefore, this study recommends protection of soil C sequestration through improved land use system is a better strategy to reduce increasing concentration of atmospheric CO2.

5. Conclusions

In order to identify and understand the influence of various land use and land cover on soil carbon storage, as well as how various land use and land management practices influence soil C storage in soils, a better understanding of the impact of various land use on soil carbon storage is necessary. This study highlights the first reporting of variations in SOC and N stocks from three distinct land use types in the Lal Bakaiya watershed, Makawanpur District, Nepal. The study found that both land use types and soil depth significantly affected the amount of SOC and N stocks in the soil. The mean total SOC stock observed in the order as natural forest > grassland > cultivated land. Accordingly, the mean total N stock followed the order as natural forest > grassland > cultivated land. Land uses with forest cover have higher SOC and N stocks than those with less vegetation. As the study found that land use types do have significant impact on SOC and N storage capacity of soil, we recommend that soil C sequestration through an improved land use system is a better strategy to reduce increasing concentration of atmospheric C. Though there is lack of proper knowledge about C trading in Nepal, the present study also indicates the significance of watershed-level SOC valuation for better and C-friendly land use decision making in the future.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

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

Authors’ Contributions

Pramod Ghimire and Uchita Lamichhane conceptualized the study, performed data curation, and investigated the data; Pramod Ghimire, Uchita Lamichhane, and Sandesh Bolakhe performed formal analysis, contributed to methodology, and prepared the original draft; Pramod Ghimire validated the data; Pramod Ghimire, Uchita Lamichhane, Sandesh Bolakhe, and Chun-Hung Jeb Lee reviewed and edited the manuscript.

Acknowledgments

The authors are grateful to Provincial Soil Laboratory, Bagmati Province, Hetauda, for extending laboratory facilities.