Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Promoting Sustainable Development of Coal Mines: CNN Model Optimization for Identification of Microseismic Signals Induced by Hydraulic Fracturing in Coal Seams
Previous Article in Journal
How Environmental Policy Perception and Social Media Use Impact Pro-Environmental Behavior: A Moderated Mediation Model Based on the Theory of Planned Behavior
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effects of Land Use and Landform Transformation on the Vertical Distribution of Soil Nitrogen in Small Catchments

1
Academy of Data Science and Informetrics, Hangzhou Dianzi University, Hangzhou 310018, China
2
Chinese Academy of Science and Education Evaluation, Hangzhou Dianzi University, Hangzhou 310018, China
3
State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences, Xi’an 710061, China
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(17), 7590; https://doi.org/10.3390/su16177590
Submission received: 25 July 2024 / Revised: 25 August 2024 / Accepted: 30 August 2024 / Published: 2 September 2024

Abstract

:
The diversity of land use and consolidation is fundamental to ensuring sustainable development. However, the impact of diverse land uses and consolidation on the well-known shallow accumulation pattern of soil nitrogen (N) remains unclear. This existence of this knowledge gap severely constrains the sustainable production of newly created farmland. Therefore, the objective of this study was to investigate the effects of land use and gully land transformation on the vertical distribution of soil N in agricultural and nature catchments. Methodologically, soil nitrate (NO3), ammonium (NH4+) and total nitrogen (TN) were measured to a depth of 100 cm in the hillslope forestland, grassland and gully cropland areas of the treated (gully landform reshaping) and untreated (natural gully) catchments on the Chinese Loess Plateau (CLP). The results indicated that soil N in the hillslope forestland and grassland exhibited a shallow accumulation pattern, while the vertical distribution of soil N in the gully cropland areas displayed a homogeneous, random or deep accumulation pattern. In the hillslope areas, vegetable cover was the dominant factor controlling N variation in the topsoil. In contrast, in the gully areas, the interaction of landform transformation and hydrology was the primary factor influencing the distribution of soil N. In the treated catchment, soil NO3 exhibited a significant deep accumulation pattern in the newly created farmland through gully landform reshaping. In the untreated catchment, soil NH4+ showed a significant deep accumulation pattern in the undisturbed natural gully. This study provides valuable insights into how land use and gully landform transformation affect the soil N profile. This information is crucial for the sustainable development and scientific management of valley agriculture at the catchment scale.

1. Introduction

Land consolidation plays a crucial role in sustainable agricultural development by optimizing land use and enhancing productivity [1]. It contributes to the achievement of the Sustainable Development Goals (SDGs) such as ending poverty (SDG 1), ending hunger (SDG 2) and Life on Land (SDG 15) through increasing the economic viability of farming, agricultural productivity and food security, sustainable land management practices that prevent land degradation, and encouraging the preservation of ecosystems [2]. Land consolidation involves reorganizing fragmented land parcels into more efficient and larger plots, facilitating modern farming practices and the use of advanced agricultural technology [3]. Consequently, the landform and land uses are significantly changed after large-scale land consolidation, leading to substantial variations in the distribution of nutrient elements [4]. Soil nitrogen (N) is a vital element for maintaining soil fertility and enhancing ecosystem productivity [5,6]. Understanding the distribution pattern of soil N is crucial for effective field management. Over the past few decades, the vertical patterns of soil N at regional and global scales have been extensively investigated [7,8,9]. However, soil N cycling at reshaped catchments has received less attention, despite being a significant component of the terrestrial N cycling processes. Therefore, strengthening the research on soil N cycling in reshaped catchments is essential to enhance our understanding of soil N dynamics and management in the ecosystem.
Soil N concentration generally reflects the combined effects of climate, vegetation and topography, and its vertical distribution can provide insights into N inputs, outputs and cycling processes [10,11,12]. The vertical distribution of soil N is primarily controlled by four critical processes: atmospheric deposition, parent material weathering, ecosystem biological cycling and physical leaching [13]. Among these, biological cycling (i.e., plant uptake, root function and leaf litter return) is considered the most important factor regulating soil N cycling and profile distribution in terrestrial ecosystems [14,15,16]. Many studies have demonstrated that due to biological N cycling, the soil N profile generally exhibits a shallow accumulation pattern [10,17,18]. Land use is one of the most significant factors affecting biological N cycling and the vertical distribution of soil N [19]. However, the effects of land use are generally limited to the soil surface [20,21,22]. For instance, Wei et al. [23] found that soil total nitrogen (TN) in the 0–40 cm soil layer differed among various land use types, whereas it tended to be similar in the subsoil layer (40–100 cm).
Topography is also an important factor affecting the redistribution of soil N through hydrological and erosional processes [24,25,26]. For example, Kosmas et al. [27] demonstrated that soil nutrients can accumulate in gully areas due to runoff and sediment concentration. Additionally, Xue et al. [28] found that soil nitrate concentration in valley areas was significantly higher than that in top-of-mound crest and mound slope areas due to nutrient transportation. At the catchment scale, hydrology plays a crucial role in controlling the spatial and vertical distribution of soil N [26]. Hydrology often interacts with topography to influence the distribution and cycling of soil N. Different slope positions typically have distinct water regimes, leading to varying patterns of soil nutrient accumulation [29,30]. As mentioned earlier, lower slope positions (i.e., valley or gully, foot slope) generally exhibit higher soil N accumulation compared to upper slope positions (i.e., slope, summit or shoulder) in the watershed due to runoff and sediment concentration [12]. Therefore, soil N distribution at the catchment scale is generally controlled by the combined effects of land use, topography and hydrology.
Since 2000, the Chinese government has initiated the “Grain for Green Program” (GGP) on the Chinese Loess Plateau (CLP) to control hillslope erosion. After nearly 20 years of restoration, the vegetation cover on the CLP doubled, increasing from 32% in 1999 to 60% in 2013 [31]. However, the large-scale vegetation rehabilitation on the CLP has resulted in a substantial decrease in hillslope farmlands. To offset the loss of farmlands and reduce land reclamation on the slopes, the CLP launched a megaproject called “Gully Land Consolidation” (GLC) in 2011 [32,33,34,35]. The main approach of GLC includes reshaping the creek valleys by incising foot slopes, filling stream channels and ditches, constructing drainage canals, dams and reservoirs, and creating flat farmlands. This project improves the sustainable use of gully agricultural land and water resources in the Loess Plateau, which is vital for achieving the SDGs [36]. At the catchment scale, such large-scale modification of the landform can remarkably change the gully topography and, consequently, the hydrological pathways and material circulation. However, how such landform transformation affects the soil N distribution and how it differs from undisturbed catchments remains unknown.
To address the above questions, two adjacent small catchments on the CLP were selected: one catchment with hillslope revegetation and gully landform transformation and one catchment with only hillslope revegetation. The objectives of this study were to (1) characterize the vertical distribution patterns of soil NO3, NH4+ and TN under different land uses, (2) determine whether land use significantly affects deep soil N variation at the catchment scale and (3) explore the effect of gully landform transformation on the vertical distribution and sustainable management of soil N.

2. Materials and Methods

2.1. Study Site

This study was conducted in two small catchments locating in the Gutun watershed of Ganguyi Town, Yan’an City, Shaanxi Province, China (latitude: 36°45.32′–36°50.60′ N, longitude: 109°46.16′–109°51.30′ E) (Figure 1). The mean annual temperature in the area is 10.3 °C, and the mean annual rainfall is 495 mm, of which approximately 70% falls from June to September. The dominant soil type is loessial soil, which corresponds to Clacis Cambisols in the FAO soil taxonomy [37]. The dominant tree, shrub and grass species are Robinia pseudoacacia, Sophoraviciifolia Hippophae and Artemisia sacrorum, respectively.
The areas of the treated (S1) and untreated (S2) catchments are 44.7 and 47.0 ha, respectively, with similar elevations ranging from 972.8 m to 1188.3 m a.s.l. In 2012, the GLC project was implemented in one of the two catchments (treated, S1); therefore, the land use type of the gully area in the catchment was converted to flat farmland. The untreated catchment (S2) represents a natural catchment, and the gully area in this catchment has not been reshaped by the GLC project. A vegetation rehabilitation project (Green for Grain Project) has been implemented in both catchments since 1999. The hilly slopes and ridges in the two catchments currently consist of forestland, shrubland and grassland. The two catchments present marked differences in land use types (Figure 2). Forestland is the dominant land use type in the treated catchment (67.9%), followed by grassland (15.0%), gully farmland (8.4%), shrubland (5.6%) and other land use types (3.1%) (Table 1). The untreated catchment is composed of 33.2% forestland, 34.0% shrubland, 30.2% grassland and 2.6% other land use types (Table 1).
Because of the same dominant soil type and the relatively small area in the two catchments, we assumed that the effects of atmospheric deposition and parent material weathering on the vertical distribution of soil N in the two catchments are equivalent; so, these effects were not explored and discussed in this study [38].

2.2. Sampling

Figure 3 is the flowchart regarding sampling, laboratory procedures and data analysis. Based on an 80 m × 80 m grid on Google maps in the slope area and a 30 m sampling interval in the gully area, 81 and 72 sampling sites were selected in the treated (S1) and untreated (S2) catchments, respectively, ensuring that the sampling sites were uniformly distributed and represented all land uses in each watershed. In each grid cell, one representative site where the dominant land use appeared was selected to collect soil samples. At each sampling site, the soil sampling depth was 100 cm, and soil samples at the layers of 0–20, 20–40, 40–60, 60–80 and 80–100 cm were collected using soil augers. A total of 725 soil samples were obtained from the two catchments, comprising 385 soil samples from the treated catchment and 340 soil samples from the untreated catchment. During the process of soil sampling, each soil sample was divided into three parts, which were used for measuring soil moisture, inorganic N (NO3 and NH4+) and soil TN.

2.3. Laboratory Analysis

The first part of each soil sample was stored in an aluminum box (diameter, 40 mm; depth, 23.5 mm) to determine the soil moisture content, and the sample was oven-dried at 105 °C to a constant mass [39]. The second part of the soil sample was passed through a 1 mm mesh and then stored at 4 °C to prevent microbial activity that could alter NO3 and NH4+ concentrations. We placed 5 g soil of each sample in an extraction vessel, added 50 mL of 1 M KCl solution to the vessel (at a soil/water ratio of 1:10), then shook the mixture on a mechanical shaker for 1 h to ensure thorough mixing and extraction of NO3 and NH4+, and subsequently centrifuged the mixture and filtered the supernatant through Whatman No. 42 filter paper. The clear extract was used to measure NO3 and NH4+ using the colorimetric method, and this was analyzed by using an automatic continuous flow analyzer (AutAnalyel, Bran + Luebbe GmbH, Krefeld, Germany) [40]. The third part of the soil sample was passed through a 0.25 mm mesh and was used to determine the soil TN concentration using the Kjeldahl digestion procedure in a Kjeltec auto analyzer (Kjeltec 8400, FOSS, Hillerød, Denmark) [41].

2.4. Statistical Analysis

In this study, the number of soil sample sites was limited in the shrublands because they are located in the deep hillslopes and soil sampling is very difficult to carry out in that area (Table 1). Therefore, data analysis of soil N in the shrublands was excluded. In this paper, means and standard errors of soil N are calculated, and significant differences are obtained. The Mann–Whitney (M-W) U test is a nonparametric statistical test used to compare differences between two independent groups. Unlike parametric tests (the t-test), the M-W U test does not assume that the data follow a normal distribution, making it particularly useful when dealing with non-normally distributed continuous data. Therefore, the nonparametric Mann–Whitney U test was used to compare the differences in soil N among the different land uses, soil layers and catchments because most data did not follow a normal distribution [42]. The most common method, Bonferroni correction, was applied for multiple comparisons. The alpha level of the significant difference test was 0.05; that is to say, if the p-value is less than 0.05, we reject the null hypothesis, concluding that there is a statistically significant difference between groups. The results of the M-W U test can reveal the effect of land use and landform transformation on soil N concentration to some extent. Moreover, correlation analysis was used to determine the relationships between soil N and soil moisture. The results of the correlation analysis can provide valuable information about how hydrology affects soil N changes. All calculations and statistical analyses were conducted using R 4.4.0 software (R Development Core Team, 2024).

3. Results

3.1. Vertical Distribution Patterns of Soil N under Different Land Use Types

Figure 4 illustrates the vertical distribution patterns of soil NO3, NH4+ and TN across different land use types in the two catchments. In the forestland and grassland, soil NO3 was depleted from 0 to 60 cm, and then stabilized at greater depths (60–100 cm) (Figure 4a,d). The distribution of soil NH4+ showed a homogeneous or random pattern in the same areas of both catchments (Figure 4b,e). Additionally, soil TN exhibited a typical shallow distribution pattern in the forestland and grassland. However, the vertical distribution patterns of soil NO3, NH4+ and TN in the gully land were significantly different from those in the forestland and grassland (Figure 4). In the gully farmland of the treated catchment, soil NO3 showed a significant accumulation pattern below 60 cm, increasing by 3.97 mg kg−1 from 60 to 100 cm—much greater than that in the natural gully of the untreated catchment (0.61 mg kg−1). In contrast, the vertical distribution of soil NH4+ in the gully of the untreated catchment revealed significant accumulation below 60 cm, with NH4+ increasing by 4.50 mg kg−1, which was much greater than the increase observed in the treated catchment’s gully land of 0.67 mg kg−1.

3.2. Soil N Variation among Different Land Use Types in the Same Soil Layer

Figure 5 illustrates the variations in soil N concentration among the different land use types in the same soil layer. Land use significantly affected soil NO3 levels in the various soil layers of the two catchments, except for the soil layer of 40–60 cm in the treated catchment. In both catchments, soil NO3 concentrations in the soil layers of 0–20 and 20–40 cm were significantly higher in the forestland compared to the grassland and gully land. However, in the 60–80 and 80–100 cm soil layers, soil NO3 concentrations were much higher in the gully land than in the forestland and grassland. Significant differences in soil NH4+ levels among the different land uses were observed in each soil layer of the two catchments, except for the 0–20 and 20–40 cm soil layers in the treated catchment. The highest soil NH4+ concentrations occurred in the natural catchment’s gully. In both catchments, topsoil TN concentrations (0–20 and 20–40 cm soil layers) in the gully land were significantly lower than those in the forestland and grassland. However, in the untreated catchment, subsoil TN concentration (80–100 cm soil layer) in the gully land was considerably higher than the concentrations in the forestland and grassland.

3.3. Soil N Variation between Treated and Untreated Catchments

Figure 6 illustrates the differences in soil N concentrations between the treated and untreated catchments within the same land use type and soil layer. In the forestland and grassland, surface soil NO3 and both surface and subsoil NH4+ concentrations were significantly higher in the treated catchment compared to the untreated catchment, while no significant difference was observed for TN in these two land uses. However, there were significant differences in soil TN in the gully land between the two catchments, but no significance was found for soil inorganic N. Although subsoil NO3 concentrations showed no significance across all land uses, they were higher in the treated catchment than in the untreated catchment, particularly in the gully land. Conversely, subsoil NH4+ concentrations in the gully land were higher in the untreated catchment compared to the treated catchment.

3.4. Vertical Distribution of Soil Moisture and Its Relationship with Soil N

Figure 7 illustrates the average soil water contents across the different land use types in the two catchments. Mean soil moisture in the forest-dominated catchment (treated, S1) was lower than that in the mixed land use catchment (untreated, S2) (Figure 7a). Additionally, soil moisture in the gully land was significantly higher than that in the forestland and grassland (Figure 7b). To investigate the impact of soil moisture on soil N across different soil layers, correlations between soil N and soil moisture were determined (Figure 8). The relationship between soil NO3 and moisture shifted from negative to positive with increasing soil depth in both catchments. This suggested that less soil NO3 migrates from the topsoil to the subsoil when soil moisture is low. Whereas when soil moisture is high, a significant amount of soil NO3 leaches down to the subsoil and accumulates in the deep layers. The same pattern was also observed in the relationship between soil TN and moisture. However, the correlation between soil NH4+ and moisture was positive across all soil layers in both catchments, with the exception of the 20–40 cm soil layer in the treated catchment.

4. Discussion

4.1. Effects of Land Use on Vertical Distribution of Soil N

Our results showed that soil NO3 and TN in the forestland and grassland located along the slope and ridge areas of the two catchments presented shallow distribution pattern, and significant differences in soil N between the forestland and grassland were limited to the soil layers spanning 0–60 cm in depth in both the treated and untreated catchments. These findings are consistent with many studies [20,23,43]. This corroborated the finding of Carter et al. (1997) [17], who found a similar vertical distribution pattern of soil TN in the 16 different agroecosystem sites, in which approximately 30, 28, 25 and 16% of the soil TN mass was distributed in the 0–10, 10–20, 20–40 and 40–60 cm soil layers, respectively. Previous studies have demonstrated that the shallow distribution of soil N is mostly controlled by biological N cycling under different land uses and land cover [44,45,46,47,48]. Jobbágy and Jackson (2001) [10] demonstrated that the shallow distributions of the most limiting soil nutrients for plants, including N, were dominated by plant cycling (i.e., litter return and root uptake and excretion). This mechanism was supported by Francaviglia et al. (2014) [49], who found crop residues and root turnover to be important factors explaining the higher soil TN concentrations and stocks in the topsoil. Many studies demonstrated that the heterogeneity of soil properties would also significantly affect biological N cycling, especially in the topsoil, which showed higher heterogeneity of soil organic matter, root mass, soil moisture, temperature and microbial activity than in the subsoil and thus showed more significant differences in topsoil N variation [20,21,50].
In this study, the area of gully land showed a different distribution pattern of soil N when compared with the forestland and grassland, which was probably attributed to the concentrated soil nutrients from the surrounding hillslopes. Numerous studies have demonstrated that soil nutrients, including soil N, concentrate in the area of gullies due to the concentration effects of runoff and sediments from the hillslopes [27,51]. Moreover, soil moisture and groundwater tables in the area of gullies are high and thus can result in intensive leaching of soil N, especially nitrate. Tan et al. [52] found that large amounts of soil N could be transported down to the deep soil layer by strong leaching processes in the gully land. In addition, we found that the content of soil N in the subsoil was positively correlated with the content of soil moisture (Figure 8), which further demonstrates that higher soil moisture can lead to greater soil N leaching. Therefore, the area of gully land showed higher subsoil N than that in the forestland and grassland due to high leaching effects (Figure 5).
The comparison of soil N between the treated and untreated catchments revealed that, with the exception of subsoil NH4+ in the gully land, the concentrations of soil inorganic N across all land uses were higher in the treated catchment than in the untreated catchment. Given the similar elevations, soil types and vegetation in the two catchments, these differences can be attributed to the distinct slope aspects between them. Similar findings were reported by Yu et al. [53]. Numerous studies have highlighted that slope aspect is an important topographical factor influencing the spatial variation of soil N [21]. Our results showed that while no significance was found for soil inorganic N in the gully land, notable differences in subsoil NO3 and NH4+ were observed between the treated and untreated catchments. Additionally, soil TN concentrations in the gully land of the untreated catchment were higher than those of the treated catchment. These findings suggest that changes in land use within gully land might alter the patterns of soil N accumulation.

4.2. Effects of Landform Transformation on Accumulation of Soil N

In this study, we found that the most significant differences in the vertical patterns of soil N between the treated and untreated catchments were located at the gully area and these differences were mostly due to gully landform transformation through the GLC project. In the treated catchment that carried out the GLC project, the gully channels had been flattened, and thus, the streamflow was blocked, which could lead to the surface runoff and sediments being concentrated in the channel areas (created flat farmland). Moreover, the organic and inorganic N fertilizers applied in the farmland would also leach down to the subsoil and accumulate, likely further inducing water pollution [54]. Sun [55] found that the GLC project could lead to a 27–45% decrease in surface runoff and a 55–73% increase in soil moisture in the reshaped creek valley. Zhao et al. [38] found that the reshaped catchment could significantly intercept rainfall runoff and increase local infiltration after the implementation of the GLC project. Sun et al. [56] also simulated the effects of GLC on runoff–sediment–nitrogen transportation and found that nitrate output in the surface runoff decreased by 31–48% after GLC. Therefore, gully landform transformation is the primary factor that leads to subsoil NO3 accumulation in the reshaped catchment.
In the catchment without conducting the GLC project, there are no blocking effects and surface runoff is intense during rainstorm events [57]. Therefore, larger amounts of soil nitrate are lost through surface runoff. However, compared with the treated catchment, the area of gully in the untreated catchment exhibited a significant accumulation of NH4+ in the subsoil, which was possibly due to a relatively higher ammonification rate and lower nitrification rate under the high soil moisture environment [58]. Rui and Wienhold [59] also found that when soil moisture levels were greater than 80% (v/v), NO3 concentration declined rapidly and NH4+ concentration increased due to the existence of anaerobic conditions. Moreover, Jackson-Blake et al. [60] found that the high soil NH4+ level under the high soil moisture environment resulted from the high net mineralization input. In this study, we also found that the content of soil NH4+ was positively correlated with soil moisture.

4.3. The Implications of Land Consolidation on Sustainable Agricultural N Management

Gully land consolidation has been implemented in many mountainous agricultural areas worldwide [61,62]. This process involves reorganizing the fragmented and scattered plots into larger, more contiguous units, leading to various environmental and agronomic changes [63]. For example, the potential for increased mechanization and intensification of agriculture could result in higher N inputs [64]. However, few studies have addressed the effect of gully land consolidation on soil N distribution in the Chinese Loess Plateau and how to manage N fertilizers through scientific strategies. This study is a pioneering attempt to explore the impacts of gully land consolidation on the vertical levels of soil N. Our findings showed that agricultural land consolidation alters the vertical distribution pattern of soil NO3, resulting in subsoil accumulation. The deep pattern of soil NO3 increases the risk of groundwater contamination, which poses a threat to the drinking water source of local villages and towns. Similar results were found by Yu et al. [54]. Soil NH4+ and TN concentrations of surface soil were relatively low in the gully farmland.
Therefore, actions must be taken to prevent drinking water contamination from soil NO3 pollution. First, precision agriculture technologies such as GPS-guided equipment, remote sensing and site-specific soil testing can be easily applied to the uniform and larger-sized consolidated farmland. These technologies ensure that N fertilizers are applied in the right amounts, at the right time and in the right place [65,66]. For example, since corn is the main crop in the study areas, we could monitor maize growth and soil N dynamics using these technologies. Consequently, precision fertilization and optimal use of residual soil N could be implemented to increase yield, reduce nitrogen loss and minimize environmental pollution. Second, consolidated farmland facilitates the logistics of collecting, storing and applying organic amendments, making integrated nutrient management more feasible and effective. This approach can enhance soil fertility and microbial activity, improving the efficiency of N utilization [67,68]. In general, gully land consolidation, coupled with scientific farmland management, could promote more balanced and sustainable agricultural development.

5. Conclusions

Soil NO3, NH4+ and TN concentrations at a depth of 0–60 cm were significantly higher than those at a depth of 60–100 cm in the forestland and grassland, which exhibited a shallow distribution pattern in the treated and untreated catchments due to biological cycling. However, the gully land, characterized by superior hydrological conditions, displayed different vertical distributions of soil N in the two catchments. In the treated catchment, where valley reshaping had occurred, higher NO3 levels were observed in the 60–100 soil layer, showing a pattern of subsoil accumulation. Conversely, the untreated catchment, which features a natural gully landform, showed significantly higher subsoil NH4+ concentrations compared to the surface soil, which displayed a pattern of subsoil accumulation. The GLC project changed the gully landform and, consequently, the hydrological pathways within the catchment, intensifying soil NO3 leaching and leading to subsoil NO3 accumulation. Overall, land use and landform changes significantly affect the vertical distribution of soil N. Precision agriculture technologies and integrated nutrient management can address the subsoil NO3 problem and support sustainable agricultural development in the reshaped farmland.

Author Contributions

Conceptualization, Y.Y., S.W. and J.Q.; methodology, Y.Y.; software, Y.Y.; validation, Y.Y., S.W. and J.Q.; formal analysis, Y.Y.; investigation, Y.Y.; resources, Y.Y.; data curation, Y.Y., S.W. and J.Q.; writing—original draft preparation, Y.Y.; writing—review and editing, Y.Y., S.W. and J.Q.; visualization, Y.Y.; supervision, Y.Y.; project administration, S.W. and J.Q.; funding acquisition, Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the open-end funds of the State Key Laboratory of Loess and Quaternary Geology (SKLLQG2026), the Research Startup Fund of Hangzhou Dianzi University (KYS265624173), and the Basic Research Expenses for Provincial Colleges and Universities—Special Project for Young Teachers’ Research and Innovation (General Project, KYS268224002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the author Y.Y. at [email protected].

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Nilsson, P. The role of land use consolidation in improving crop yields among farm households in Rwanda. J. Dev. Stud. 2018, 55, 1726–1740. [Google Scholar] [CrossRef]
  2. Asiama, K.Q.; Voss, W.; Bennett, R.; Rubanje, I. Lan consolidation activities in Sub-Saharan Africa towards the agenda 2030: A tale of three countries. Land Use Policy 2021, 101, 105140. [Google Scholar] [CrossRef]
  3. Zeng, S.; Zhu, F.; Chen, F.; Yu, M.; Zhang, S.; Yang, Y. Assessing the impacts of land consolidation on agricultural technical efficiency of producers: A survey from Jiangsu Province, China. Sustainability 2018, 10, 2490. [Google Scholar] [CrossRef]
  4. Yu, Q.; Zeng, Q.; Yu, G. The influence of land consolidation on biomass and ecologiacl environment. Res. J. App. Sci. Eng. Technol. 2014, 7, 3656–3662. [Google Scholar] [CrossRef]
  5. Shah, Z.; Shah, S.H.; Peoples, M.B.; Schwenke, G.D.; Herridge, D.F. Crop residue and fertiliser N effects on nitrogen fixation and yields of legume–cereal rotations and soil organic fertility. Field Crop Res. 2003, 83, 1–11. [Google Scholar] [CrossRef]
  6. Dong, J.; Hengsdijk, H.; Dai, T.B.; Boer, W.D.; Qi, J.; Cao, W.X. Long-Term Effects of Manure and Inorganic Fertilizers on Yield and Soil Fertility for a Winter Wheat-Maize System in Jiangsu, China. Pedosphere 2006, 16, 25–32. [Google Scholar]
  7. Ding, F.; Hu, Y.-L.; Li, L.-J.; Li, A.; Shi, S.; Lian, P.-Y.; Zeng, D.-H. Changes in soil organic carbon and total nitrogen stocks after conversion of meadow to cropland in Northeast China. Plant Soil 2013, 373, 659–672. [Google Scholar]
  8. Li, D.; Niu, S.; Luo, Y. Global patterns of the dynamics of soil carbon and nitrogen stocks following afforestation: A meta-analysis. New Phytol. 2012, 195, 172–181. [Google Scholar]
  9. Batjes, N.H. Total carbon and nitrogen in the soils of the world. Eur. J. Soil Sci. 1996, 47, 151–163. [Google Scholar]
  10. Jobbágy, E.G.; Jackson, R.B. The Distribution of Soil Nutrients with Depth: Global Patterns and the Imprint of Plants. Biogeochemistry 2001, 53, 51–77. [Google Scholar]
  11. Tian, L.; Zhao, L.; Wu, X.; Fang, H.; Zhao, Y.; Yue, G.; Liu, G.; Chen, H. Vertical patterns and controls of soil nutrients in alpine grassland: Implications for nutrient uptake. Sci. Total Environ. 2017, 607–608, 855–864. [Google Scholar]
  12. Zhu, H.; Wu, J.; Guo, S.; Huang, D.; Zhu, Q.; Ge, T.; Lei, T. Land use and topographic position control soil organic C and N accumulation in eroded hilly watershed of the Loess Plateau. Catena 2014, 120, 64–72. [Google Scholar]
  13. Trudgill, S.T. Soil and Vegetation Systems; Oxford Universtiy Press: New York, NY, USA, 1988. [Google Scholar]
  14. Fu, X.; Shao, M.; Wei, X.; Horton, R. Soil organic carbon and total nitrogen as affected by vegetation types in Northern Loess Plateau of China. Geoderma 2010, 155, 31–35. [Google Scholar]
  15. Debela, N.; Gebrekidan, H. Effect of land use changes and soil depth on soil organic matter, total nitrogen and available phosphorus contents of soils in Senbat watershed, western Ethiopia. ARPN J. Agri. Biol. Sci. 2013, 8, 206–212. [Google Scholar]
  16. Mekonnen, K.; Buresh, R.J.; Coe, R.; Kipleting, K.M. Root length and nitrate under Sesbania sesban: Vertical and horizontal distribution and variability. Agroforest. Syst. 1998, 42, 265–282. [Google Scholar]
  17. Carter, M.R.; Angers, D.A.; Gregorich, E.G.; Bolinder, M.A. Organic carbon and nitrogen stocks and storage profiles in cool, humid soils of eastern Canada. Can. J. Soil Sci. 1997, 77, 205–210. [Google Scholar]
  18. Chai, H.; Yu, G.; He, N.; Wen, D.; Li, J.; Fang, J. Vertical distribution of soil carbon, nitrogen, and phosphorus in typical Chinese terrestrial ecosystems. Chin. Geogr. Sci. 2015, 25, 549–560. [Google Scholar]
  19. Jin, Z.; Zhu, Y.; Li, X.; Dong, Y.; An, Z. Soil N retention and nitrate leaching in three types of dunes in the Mu Us desert of China. Sci. Rep. 2015, 5, 207–216. [Google Scholar]
  20. Omonode, R.A.; Vyn, T.J. Vertical distribution of soil organic carbon and nitrogen under warm-season native grasses relative to croplands in west-central Indiana, USA. Agric. Ecosyst. Environ. 2006, 117, 159–170. [Google Scholar]
  21. Yimer, F.; Ledin, S.; Abdelkadir, A. Soil organic carbon and total nitrogen stocks as affected by topographic aspect and vegetation in the Bale Mountains, Ethiopia. Geoderma 2006, 135, 335–344. [Google Scholar]
  22. Zhou, X.; Dong, H.; Lan, Z.; Bacon, G.; Hao, Y.; Chen, C. Vertical distribution of soil extractable organic C and N contents and total C and N stocks in 78-year-old tree plantations in subtropical Australia. Environ. Sci. Pollut. Res. 2017, 24, 22312–22320. [Google Scholar]
  23. Wei, X.R.; Shao, M.G.; Fu, X.L.; Horton, R.; Li, Y.; Zhang, X.C. Distribution of soil organic C, N and P in three adjacent land use patterns in the northern Loess Plateau, China. Biogeochemistry 2009, 96, 149–162. [Google Scholar]
  24. Burt, T.P.; Butcher, D.P. Topographic controls of soil moisture distributions. Eur. J. Soil Sci. 2010, 36, 469–486. [Google Scholar]
  25. Perakis, S.S.; Tepley, A.J.; Compton, J.E. Disturbance and Topography Shape Nitrogen Availability and δ15N over Long-Term Forest Succession. Ecosystems 2015, 18, 573–588. [Google Scholar]
  26. Zhu, Q.; Schmidt, J.P.; Lin, H.S.; Sripada, R.P. Hydropedological processes and their implications for nitrogen availability to corn. Geoderma 2009, 154, 111–122. [Google Scholar]
  27. Kosmas, C.; Danalatos, N.; Cammeraat, L.H.; Chabart, M.; Diamantopoulos, J.; Farand, R.; Gutierrez, L.; Jacob, A.; Marques, H.; Martinez-Fernandez, J.; et al. The effect of land use on runoff and soil erosion rates under Mediterranean conditions. Catena 1997, 29, 45–59. [Google Scholar]
  28. Xue, Z.; Cheng, M.; An, S. Soil nitrogen distributions for different land uses and landscape positions in a small watershed on Loess Plateau, China. Ecol. Eng. 2013, 60, 204–213. [Google Scholar]
  29. Elwell, H.A.; Stocking, M.A. Vegetal cover to estimate soil erosion hazard in Rhodesia. Geoderma 1976, 15, 61–70. [Google Scholar]
  30. Francis, C.F.; Thornes, J.B. Runoff hydrographs from three Mediterranean vegetation cover types. Veg. Erosion: Process Environ. 1990, 363–384. [Google Scholar]
  31. Chen, Y.; Wang, K.; Lin, Y.; Shi, W.; Song, Y.; He, X. Balancing green and grain trade. Nat. Geosci. 2015, 8, 739–741. [Google Scholar]
  32. Liu, Q.; Wang, Y.; Zhang, J.; Chen, Y. Filling gullies to create farmland on the loess plateau. Environ. Sci. Technol. 2013, 47, 7589–7590. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, X.B.; Jin, Z. Gully land consolidation project in Yan’an is inheritance and development of wrap land dam project on the Loess Plateau. J. Earth Environ. 2015, 6, 251–264. (In Chinese) [Google Scholar]
  34. Li, X.; Jin, Z.; Zhang, X.B.; Zhou, W.J. Analysis of ecosystem management of the Loess Plateau during the past 60 years and suggestions for the future development. J. Earth Environ. 2015, 6, 248–254. (In Chinese) [Google Scholar]
  35. Liu, Y.; Guo, Y.; Li, Y.; Li, Y. GIS-based effect assessment of soil erosion before and after gully land consolidation: A case study of Wangjiagou project region, Loess Plateau. Chin. Geogr. Sci. 2015, 25, 137–146. [Google Scholar]
  36. Wu, G.; Fu, S.H.; Zhou, G.Y.; Yin, B. Effect of the gully land consolidation project on soil erosion and sediment yield on the Loess Plateau, China. Land Degrad. Dev. 2023, 34, 1464–1476. [Google Scholar]
  37. WRB, I.W.G. World Reference Base for Soil Resources, 2nd ed.; FAO: Rome, Italy, 2006; p. 127. [Google Scholar]
  38. Zhao, Y.; Wang, Y.; Wang, L.; Zhang, X.; Yu, Y.; Jin, Z.; Lin, H.; Chen, Y.; Zhou, W.; An, Z. Exploring the role of land restoration in the spatial patterns of deep soil water at watershed scales. Catena 2019, 172, 387–396. [Google Scholar] [CrossRef]
  39. Wang, Y.; Shao, M.a.; Liu, Z. Large-scale spatial variability of dried soil layers and related factors across the entire Loess Plateau of China. Geoderma 2010, 159, 99–108. [Google Scholar]
  40. Keeney, D.R.; Nelson, D.W. Nitrogen-Inorganic Forms, in Methods of Soil Analysis, 2nd ed.; ASA and SSSA: Madison, WI, USA, 1982; pp. 643–698. [Google Scholar]
  41. Bremner, J.M.; Tabatabai, M.A. Use of an ammonia electrode for determination of ammonium in Kjeldahl analysis of soils. Commun. Soil Sci. Plant Anal. 1972, 3, 159–165. [Google Scholar]
  42. Neilen, A.D.; Chen, C.R.; Parker, B.M.; Faggotter, S.J.; Burford, M.A. Differences in nitrate and phosphorus export between wooded and grassed riparian zones from farmland to receiving waterways under varying rainfall conditions. Sci. Total Environ. 2017, 598, 188–197. [Google Scholar]
  43. Yu, Z.; Li, C.; Wang, M. National-Scale Meta-Analysis of Soil Carbon and Nitrogen Accumulation Potential in China’s Grain for Green Program. Eur. Soil Sci. 2020, 53, 892–901. [Google Scholar]
  44. Chang, R.; Jin, T.; Lü, Y.; Liu, G.; Fu, B. Soil Carbon and Nitrogen Changes following Afforestation of Marginal Cropland across a Precipitation Gradient in Loess Plateau of China. PLoS ONE 2014, 9, e85426. [Google Scholar]
  45. Liu, W.; Chen, S.; Qin, X.; Baumann, F.; Scholten, T.; Zhou, Z.; Sun, W.; Zhang, T.; Ren, J.; Qin, D. Storage, patterns, and control of soil organic carbon and nitrogen in the northeastern margin of the Qinghai-Tibetan Plateau. Environ. Res. Lett. 2012, 7, 35401–35412. [Google Scholar]
  46. Owen, A.G.; Jones, D.L. Competition for amino acids between wheat roots and rhizosphere microorganisms and the role of amino acids in plant N acquisition. Soil Biol. Biochem. 2001, 33, 651–657. [Google Scholar]
  47. Gherardi, L.A.; Yahdjian, L. Preference for different inorganic nitrogen forms among plant functional types and species of the Patagonian steppe. Oecologia 2013, 173, 1075–1081. [Google Scholar] [PubMed]
  48. Makarov, M.I.; Malysheva, T.I.; Menyailo, O.V. Isotopic Composition of Nitrogen and Transformation of Nitrogen Compounds in Meadow-Alpine Soils. Eur. Soil Sci. 2019, 52, 1028–1037. [Google Scholar]
  49. Francaviglia, R.; Benedetti, A.; Doro, L.; Madrau, S.; Ledda, L. Influence of land use on soil quality and stratification ratios under agro-silvo-pastoral Mediterranean management systems. Agric. Ecosys. Environ. 2014, 183, 86–92. [Google Scholar]
  50. Mayes, M.T.; Marinspiotta, E.; Ozdogan, M.; Erdogan, M.A. A landscape-scale study of land use and parent material effects on soil organic carbon and total nitrogen in the Konya Basin, Turkey. In Proceedings of the American Geophysical Union, Fall Meeting 2011, San Francisco, CA, USA, 5–9 December 2011. [Google Scholar]
  51. Soon, Y.K.; Malhi, S.S. Soil nitrogen dynamics as affected by landscape position and nitrogen fertilizer. Can. J. Soil Sci. 2005, 85, 579–587. [Google Scholar]
  52. Tan, X.; Shao, D.; Liu, H.; Yang, F.; Xiao, C.; Yang, H. Effects of alternate wetting and drying irrigation on percolation and nitrogen leaching in paddy fields. Paddy Water Environ. 2013, 11, 381–395. [Google Scholar]
  53. Yu, Y.; Jin, Z.; Lin, H.; Wang, Y.; Zhao, Y.; Chu, G.; Zhang, J.; Song, Y.; Zheng, H. Spatial variation and soil nitrogen potential hotspots in a mixed land cover catchment on the Chinese Loess Plateau. J. Mt. Sci. 2019, 16, 1353–1366. [Google Scholar] [CrossRef]
  54. Yu, Y.; Jin, Z.; Chu, G.; Zhang, J.; Wang, Y.; Zhao, Y. Effects of valley reshaping and damming on surface and groundwater nitrate on the Chinese Loess Plateau. J. Hydrol. 2020, 584, 124702. [Google Scholar]
  55. Sun, P. Simulation of the Effects of Gully Land Consolidation on Precipitation Transformation; The Northwest A & F University: Yangling, China, 2017. (In Chinese) [Google Scholar]
  56. Sun, P.C.; Gao, J.E.; Han, S.Q.; Yin, Y.; Zhou, M.F.; Han, J.Q. Simulation study on the effects of typical gully land consolidation on runoff-sediment-nitrogen emissions in the loess hilly-gully region. J. Agro-Environ. Sci. 2017, 36, 1177–1185. (In Chinese) [Google Scholar]
  57. Lou, X.Y.; Gao, J.E.; Han, S.Q.; Guo, Z.H.; Yin, Y. Influence of Land Consolidation Engineering of Gully Channel on Watershed Runoff Yield and Concentration in Loess Hilly and Gully Region. Water Resour. Power. 2016, 34, 23–27. (In Chinese) [Google Scholar]
  58. Vernimmen, R.R.E.; Verhoef, H.A.; Verstraten, J.M.; Bruijnzeel, L.A.; Klomp, N.S.; Zoomer, H.R.; Wartenbergh, P.E. Nitrogen mineralization, nitrification and denitrification potential in contrasting lowland rain forest types in Central Kalimantan, Indonesia. Soil Biol. Biochem. 2007, 39, 2992–3003. [Google Scholar] [CrossRef]
  59. Rui, Z.; Wienhold, B.J. The effect of soil moisture on mineral nitrogen, soil electrical conductivity, and pH. Nutr. Cycl. Agroecosys. 2002, 63, 251–254. [Google Scholar]
  60. Jackson-Blake, L.; Helliwell, R.C.; Britton, A.J.; Gibbs, S.; Coull, M.C.; Dawson, L. Controls on soil solution nitrogen along an altitudinal gradient in the Scottish uplands. Sci. Total Environ. 2012, 431, 100–108. [Google Scholar] [CrossRef]
  61. Janus, J.; Łopacka, M.; John, E.J.G. Land consolidation in mountain areas. Case study from southern Poland. Geod. Cartog. 2017, 66, 241–251. [Google Scholar] [CrossRef]
  62. Janečková Molnárová, K.; Sklenička, P.; Bohnet, I.C.; Lowther-Harris, F.; van den Brink, A.; Movahhed Moghaddam, S.; Fanta, V.; Zástěra, V.; Azadi, H. Impacts of land consolidation on land degradation: A systematic review. J. Environ. Manag. 2023, 329, 117026. [Google Scholar] [CrossRef]
  63. Ramírez del Palacio, Ó.; Hernández-Navarro, S.; Sánchez-Sastre, L.F.; Fernández-Coppel, I.A.; Pando-Fernández, V. Assessment of Land Consolidation Processes from an Environmental Approach: Considerations Related to the Type of Intervention and the Structure of Farms. Agronomy 2022, 12, 1424. [Google Scholar] [CrossRef]
  64. Huddell, A.M.; Galford, G.L.; Tully, K.L.; Crowley, C.; Palm, C.A.; Neill, C.; Hickman, J.E.; Menge, D.N.L. Meta-analysis on the potential for increasing nitrogen losses from intensifying tropical agriculture. Glob. Chang. Biol. 2020, 26, 1668–1680. [Google Scholar] [CrossRef]
  65. Sishodia, R.P.; Ray, R.L.; Singh, S.K. Applications of Remote Sensing in Precision Agriculture: A Review. Remote Sens. 2020, 12, 3136. [Google Scholar] [CrossRef]
  66. Toriyama, K.; Sadaki, R.; Shibata, Y.; Sugimoto, M.; Chosa, T.; Omine, M.; Saito, J. Development of a site-specific nitrogen management system for paddy rice. Japan Agric. Res. Quat. 2003, 37, 213–218. [Google Scholar]
  67. Toda, M.; Walder, F.; van der Heijden, M.G.A. Organic management and soil health promote nutrient use efficiency. J. Sustain. Agric. Environ. 2023, 2, 215–224. [Google Scholar]
  68. Bargaz, A.; Lyamlouli, K.; Chtouki, M.; Zeroual, Y.; Dhiba, D. Soil microbial resources for improving fertilizers efficiency in an integrated plant nutrient management system. Front. Microbiol. 2018, 9, 1606. [Google Scholar]
Figure 1. The locations of the study area and soil sampling sites in the treated and untreated catchments. S1: the treated catchment; S2: the untreated catchment.
Figure 1. The locations of the study area and soil sampling sites in the treated and untreated catchments. S1: the treated catchment; S2: the untreated catchment.
Sustainability 16 07590 g001
Figure 2. Land use types in the treated and untreated catchments. (a) The treated catchment (S1); (b) the untreated catchment (S2).
Figure 2. Land use types in the treated and untreated catchments. (a) The treated catchment (S1); (b) the untreated catchment (S2).
Sustainability 16 07590 g002
Figure 3. The technology roadmap.
Figure 3. The technology roadmap.
Sustainability 16 07590 g003
Figure 4. Average distribution patterns of nitrate (NO3), ammonium (NH4+) and total nitrogen (TN) under different land use types in the soil profile of the treated and untreated catchments. S1: the treated catchment; S2: the untreated catchment.
Figure 4. Average distribution patterns of nitrate (NO3), ammonium (NH4+) and total nitrogen (TN) under different land use types in the soil profile of the treated and untreated catchments. S1: the treated catchment; S2: the untreated catchment.
Sustainability 16 07590 g004
Figure 5. The differences in soil nitrate (NO3), ammonium (NH4+) and total nitrogen (TN) among different land uses under the same soil layer in the treated and untreated catchments. Different letters in the same soil layer indicate the significant difference in soil N concentration under different land use types (p < 0.05). S1: the treated catchment; S2: the untreated catchment.
Figure 5. The differences in soil nitrate (NO3), ammonium (NH4+) and total nitrogen (TN) among different land uses under the same soil layer in the treated and untreated catchments. Different letters in the same soil layer indicate the significant difference in soil N concentration under different land use types (p < 0.05). S1: the treated catchment; S2: the untreated catchment.
Sustainability 16 07590 g005
Figure 6. The differences in soil nitrate (NO3), ammonium (NH4+) and total nitrogen (TN) between the treated and untreated catchments. The star (*) indicates a significant difference in soil N concentration between the two catchments (p < 0.05). S1: the treated catchment; S2: the untreated catchment.
Figure 6. The differences in soil nitrate (NO3), ammonium (NH4+) and total nitrogen (TN) between the treated and untreated catchments. The star (*) indicates a significant difference in soil N concentration between the two catchments (p < 0.05). S1: the treated catchment; S2: the untreated catchment.
Sustainability 16 07590 g006
Figure 7. Soil moisture in the soil profiles of the treated and untreated catchments. Different letters in the same soil layer indicate a significant difference in soil moisture under different land uses (p < 0.05). S1: the treated catchment; S2: the untreated catchment. a: Soil moisture of all samples in the two catchments. b: Soil moisture among different land uses.
Figure 7. Soil moisture in the soil profiles of the treated and untreated catchments. Different letters in the same soil layer indicate a significant difference in soil moisture under different land uses (p < 0.05). S1: the treated catchment; S2: the untreated catchment. a: Soil moisture of all samples in the two catchments. b: Soil moisture among different land uses.
Sustainability 16 07590 g007
Figure 8. Pearson correlation analysis of soil N and soil moisture in different soil layers in the treated and untreated catchments. NN: nitrate (NO3, mg kg−1), AN: ammonium (NH4+, mg kg−1), TN: total nitrogen (g kg−1) and SM: soil moisture (%). The Arabic numerals 1, 2, 3, 4 and 5 following the capital letters represent the 0–20, 20–40, 40–60, 60–80 and 80–100 cm soil layers, respectively. The color and direction of the symbols represent the negativity or positivity of the corresponding correlation coefficients; and the areas of the symbols mean the absolute values of the corresponding correlation coefficients (the larger the area, the smaller the value). (a) The treated catchment (S1); (b) the untreated catchment (S2).
Figure 8. Pearson correlation analysis of soil N and soil moisture in different soil layers in the treated and untreated catchments. NN: nitrate (NO3, mg kg−1), AN: ammonium (NH4+, mg kg−1), TN: total nitrogen (g kg−1) and SM: soil moisture (%). The Arabic numerals 1, 2, 3, 4 and 5 following the capital letters represent the 0–20, 20–40, 40–60, 60–80 and 80–100 cm soil layers, respectively. The color and direction of the symbols represent the negativity or positivity of the corresponding correlation coefficients; and the areas of the symbols mean the absolute values of the corresponding correlation coefficients (the larger the area, the smaller the value). (a) The treated catchment (S1); (b) the untreated catchment (S2).
Sustainability 16 07590 g008
Table 1. Land use types and number of soil sampling sites in treated and untreated catchments.
Table 1. Land use types and number of soil sampling sites in treated and untreated catchments.
CatchmentLand UseNumber of
Sample Sites
CatchmentLand UseNumber of
Sample Sites
Treated
(S1)
Forestland (67.9%)47Untreated
(S2)
Forestland (33.2%)24
Shrubland (5.6%)3 Shrubland (34.0%)4
Grassland (15.0%)13 Grassland (30.2%)32
Cropland (8.4%)17 Wasteland (1.2%)12
Others (3.1%)- Others (1.4%)-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yu, Y.; Wang, S.; Qiu, J. The Effects of Land Use and Landform Transformation on the Vertical Distribution of Soil Nitrogen in Small Catchments. Sustainability 2024, 16, 7590. https://doi.org/10.3390/su16177590

AMA Style

Yu Y, Wang S, Qiu J. The Effects of Land Use and Landform Transformation on the Vertical Distribution of Soil Nitrogen in Small Catchments. Sustainability. 2024; 16(17):7590. https://doi.org/10.3390/su16177590

Chicago/Turabian Style

Yu, Yunlong, Shanshan Wang, and Junping Qiu. 2024. "The Effects of Land Use and Landform Transformation on the Vertical Distribution of Soil Nitrogen in Small Catchments" Sustainability 16, no. 17: 7590. https://doi.org/10.3390/su16177590

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop