Next Article in Journal
Assessing the Quality and Efficiency of Education in Rural Schools through the Models of Integration and the Cooperation of Educational Institutions: A Case Study of Russia and Indonesia
Previous Article in Journal
Innovation and High-Quality Development of Enterprises—Also on the Effect of Innovation Driving the Transformation of China’s Economic Development Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Land Use Change on Rainfall Erosion in Luojiang River Basin, China

School of Water Resources, North China University of Water Resources and Electric Power, Zhengzhou 450011, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(14), 8441; https://doi.org/10.3390/su14148441
Submission received: 6 May 2022 / Revised: 28 June 2022 / Accepted: 5 July 2022 / Published: 11 July 2022

Abstract

:
This paper, based on daily rainfall erosivity model, ArcGIS, trend analysis and Kriging interpolation method, analyzed the spatial and temporal distribution characteristics of rainfall erosivity in the Luojiang River Basin of China, and then explored the influence relationship between land use change types and rainfall erosivity potential. The results showed the following: (1) from 1980 to 2019, the distribution range of multi-annual rainfall erosivity in the Luojiang River Basin was 14,674–15,227 MJ·mm/ (hm2·h), with an average value of 14,102 MJ·mm/(hm2·h), showing an overall increasing trend; (2) the spatial distribution of rainfall erosivity value tends to be consistent with the multi-year average rainfall, showing a decreasing trend from the middle to the periphery of the basin; (3) land use change is an important factor affecting the spatial and temporal distribution characteristic of rainfall erosivity value in the basin. The increase in rainfall erosivity will undoubtedly increase the potential of soil erosion. This study can provide theoretical reference for future basin land use planning and put forward preventive suggestions according to the distribution characteristics of rainfall erosivity.

1. Introduction

Rainfall is the power source of regional soil and water loss. Rainfall erosivity reflects the potential ability of rainfall to soil erosion, and is the main dynamic factor causing soil erosion. It is also one of the main parameters of soil erosion, sediment yield and water environment modeling [1,2]. In 1958, Wischmeier [3] and Smith et al. [4] proposed the concept of rainfall erosivity, firstly, which was further developed by Hudson [5], Wischmeier and Smith [6,7,8] and used the product of the kinetic energy of subrainfall (E) and the maximum 30 min rainfall intensity (I30) as the measurement index [9]. However, due to the difficulty in obtaining the data required for calculation and the complexity and time-consuming nature of calculation steps [10,11], a simple algorithm for rainfall erosivity came into being [12]. Richardson et al. [13] established a rainfall erosivity calculation model based on daily rainfall scale, which is cited by many scholars. Yu et al. [4,14] introduced cosine function on the basis of daily rainfall model and analyzed rainfall in New South Wales, Australia and obtained good verification results. Duulatov et al. [8] predicted rainfall erosivity in Central Asia based on four global climate models. The research on rainfall erosivity in China began in the 1980s. Jia Zhijun [15], Wang Wanzhong [16], Huang Yanhe [17], Zhang Xiankui [18] and other scholars summarized the calculation model of rainfall erosivity suitable for different regions in China [11]. Zhang Wenbo et al. [19,20] studied and established a simple model for daily rainfall estimation of rainfall erosivity, which has been well-used in many drainage basins in southern China. Gao Ge et al. [21] discussed the temporal and spatial variation of rainfall erosivity in the past and future of the Yellow River Basin based on the multi-year daily precipitation and future medium discharge scenarios. Wang Yousheng et al. [22] introduced seasonal variation into the Richardson model, estimated rainfall erosivity of Huangshan City, Anhui Province, China, and tested the seasonal performance of the Richardson model.
The soil loss equation RUSLE (Revised Universal Soil Loss Equation) [23] model considers soil erosion as a function of climate erosion (which is affected by rainfall and intensity), topography, soil erodibility and vegetation cover [23,24], and studies the variation law of soil erosion from rainfall, soil properties, topography, land use and vegetation cover [25]. Because soil and topography have a large but relatively stable influence on water erosion, scholars attribute the increase in soil erosion to heavy rain, inappropriate land use and degraded vegetation [26,27]. Vegetation controls soil erosion through its canopy, root and litter components, but soil erosion also affects the composition, structure and growth pattern of plant communities [28,29]. However, due to their interaction, the relationship between precipitation, vegetation and erosion is complex and uncertain. Therefore, the positive or negative correlation between precipitation and erosion mainly depends on land use and vegetation [30].
Land use mainly affects soil erosion characteristics by influencing soil properties and vegetation coverage [7,31,32]. Land use changes include the quantity and quality of land resources and the change of land use structure over time, as well as the change in the spatial structure and the combination of land use types [33]. However, with the increasing intensity of human activities, unreasonable land use changes the regional topographic conditions, deteriorates soil characteristics and destroys vegetation resources, which aggravates soil erosion, limits the way and structure of land use, leads to land productivity degradation, and further intensifies the contradiction between human and land. A large number of studies have shown that the impact of rainfall on the characteristics of slope erosion is significantly different under different land use types [34]. Slopes with good vegetation coverage can effectively resist rainstorm erosion [35,36]. The degree of soil erosion in forest, shrub and dense grassland is lower than that in woodlands and sparse grasslands [37]. The greater the vegetation density, the more significant the interception effect of raindrops, preventing them from directly hitting the soil surface [38] and thus reducing the erosivity of rainfall. In addition, dense vegetation also slows down the formation of surface runoff, the root system of trees and shrubs increases soil porosity and organic matter content, thereby increasing the infiltration rates and reducing erosion. Sparse grassland and cultivated land are more likely to form runoff, so soil particles are directly exposed to the driving effect of surface runoff, which is prone to rainfall erosion [39,40].
At present, most scholars have studied the spatial and temporal characteristics of rainfall erosivity in different regions from different perspectives. However, due to the influence of vegetation coverage, the potential of rainfall erosivity under different land use conditions is also very different. The upper reaches of the Luojiang River Basin are a hilly region with multi-weathered sandy loam, loose soil, single land use type and serious soil and water loss. The middle and lower reaches are gentle slopes and plains, with red loam and yellow loam; the soil fertility is good. However, due to the serious influence of human activities, vegetation cover varies greatly, land use types are complex and changeable, and soil erosion occurs easily. In view of this, on the basis of previous studies, this paper still takes the Luojiang River Basin as an example to calculate the value and spatial–temporal variation characteristics of rainfall erosivity, and then it explores the influence relationship of rainfall erosion potential under land use change in the basin. The results not only can provide important reference for agricultural management in the basin, but they also have practical prospects for rational development and utilization of water and soil resources and soil erosion prevention.

2. Materials and Methods

2.1. Study Area

The Luojiang River is located in the low latitude area south of the North Return Line (110°10′ E~110°50′ E, 21°30′ N~22°30′ N), most of which are located in Huazhou City, Maoming City, Guangdong Province, China. It is the largest tributary of the Jianjiang River (Figure 1). Luojiang River originates from Wangjiang, Xinyi City, with a total length of 143 km, a basin area of 2055.067 km2 and an average slope of 0.64‰. This region belongs to the south subtropical monsoon climate. The annual average temperature is 23.1 °C, and the rainfall is abundant. The annual average rainfall is 1890 mm, but the distribution of rainfall in space and time is extremely uneven. The rainfall in the middle of a year is mainly concentrated in May to September, accounting for 75% of the total rainfall of the year [41,42].

2.2. Data Source

The data required in this paper include land use data (2010, 2015 and 2018) and daily rainfall data of the Hejiang, Jialong and Baowei meteorological station from 1980 to 2019. The description, pretreatment and sources of each data are shown in Table 1. Land use data of the three calculated years are shown in Figure 2.

2.3. Calculation Model of Rainfall Erosivity

The rainfall erosivity factor (R) is a basic factor in the soil loss equation, which is mainly estimated by establishing the power function structure or cosine function structure of rainfall and rainfall intensity. This paper adopts the simplified model of calculating rainfall erosivity by daily rainfall proposed by Zhang Wenbo et al. [19]. The model combines rainfall data in half a month period, uses the product of daily rainfall and daily rainfall intensity combination as the index to measure rainfall erosivity. In addition, the regional model parameters are used to modify the model, which can be well applied to the red soil region with abundant rainfall in southern China.
R = α k = 1 n P k β
α = 21.586 β 7.1891
β = 0.8363 + 18.144 P d 12 + 24.455 P y 12
where R is rainfall erosivity within half a month, MJ·mm /(hm2·h); K is erosive rainfall days within half a month; P K is the daily rainfall, which greater than or equal to 12 mm on the kth day in half a month, mm; P d 12 refers to the average daily rainfall greater than or equal to 12 mm, the unit is mm; P y 12 is the annual average rainfall with daily rainfall greater than or equal to 12 mm, the unit is mm; α and β are model parameters, the calculation formula is shown in Formulas (2) and (3).

2.4. Spatio-Temporal Analysis Method

2.4.1. Mann–Kendall Trend Analysis

Mann–Kendall (MK) is a nonparametric test method recommended by the World Meteorological Organization and has been widely used. This method judges the trend change of the sequence to be tested by the positive and negative statistics Z, which is suitable for the trend analysis of meteorological and hydrological sequences [43]. In this paper, MK trend analysis is used to analyze the change trend of multi-year rainfall erosivity in the basin from the perspective of time, and the confidence level of explicitness test is set to be 95 % (|Z| > 1.96).

2.4.2. Kriging Interpolation Method

Kriging estimation method [44] based on geostatistics is a commonly used spatial interpolation method. It can better fit the spatial pattern of rainfall erosivity by establishing the covariance function or variation function of variables and interpolating the estimated points according to the spatial autocorrelation of variables. In this paper, ArcGIS software is used to realize Kriging interpolation to analyze the spatial distribution of multi-year average rainfall erosivity in the basin.

2.5. Land Use Transfer Matrix

Land use state transition matrix is a quantitative description between system state and state transition in system analysis, which can describe the structural characteristics of regional land use change, reflect the direction of land use change, and better reveal the evolution process of land use pattern [45]. Its mathematical expression is as follows:
S i j = [ S 11 S 12 S 1 n S 12 S 22 S 2 n S n 1 S n 2 S n n ]
where S is the area, n is the number of land use types before and after the transfer, i and j are the land use types at the beginning and end of the research period, respectively, and Sij is the area of land use conversion from land type i at the beginning of the research period to land type j at the end of the research period.

2.6. Sub-Basins Division

The division of sub-basins is an indispensable step in constructing hydrological models. The SWAT (Soil and Water Assessment Tool) model is a semi-distributed hydrological model constructed and developed by DILE et al. [46] based on ArcGIS software, which can be based on digital elevation model by setting the catchment area threshold, and the subwatershed in the simulated study area can be divided.

3. Results

3.1. Spatial and Temporal Distribution Characteristics of Rainfall Erosivity

Based on the daily observed rainfall data from 1980 to 2019 at meteorological stations in the Hejiang, Jialong and Baowei, the rainfall erosivity of each station is calculated by using Equations (1)–(3). The monthly rainfall erosivity, annual rainfall erosivity and annual average rainfall erosivity can be obtained through summary

3.1.1. Temporal Variation Characteristics of Rainfall Erosivity

The trend analysis of rainfall erosivity in the Luojiang Basin from 1980 to 2019 was conducted by MK test method, as shown in Figure 3.
It can be seen that from 1980 to 2019, the temporal pattern and trend of annual rainfall erosivity and annual rainfall basically tend to be consistent, showing an obvious trend of increase, with a growth rate of 29.048 MJ·mm/ (hm2·h·a). After the 1990s, the rainfall erosivity value showed an obvious upward trend, and the inter-annual fluctuation amplitude was significantly different. The maximum annual rainfall erosivity was 24,429.774 MJ·mm/(hm2·h) in 1993, and the minimum was 7171.761 MJ·mm/(hm2·h) in 1998, with a difference of 3.4 times.

3.1.2. Spatial Distribution Characteristics of Rainfall Erosivity

The ArcGIS software was used to interpolate the whole basin by Kriging interpolation method combined with the calculation results of rainfall erosivity at each station, and the spatial distribution map of rainfall erosivity in the Luojiang River Basin was obtained (Figure 4a). The rainfall erosivity of the Luojiang Basin ranges from 14,674 to 15,227 MJ·mm/(hm2·h), with an average of 14,102.67 MJ·mm/(hm2·h), and a standard deviation of 4325.206 MJ·mm/(hm2·h), indicating a large spatial difference. The rainfall erosivity showed a decreasing trend from the middle to the surrounding areas. The rainfall erosivity was the largest in the middle and most parts of the north-central basin, with the average annual rainfall erosivity exceeding 15,146 MJ·mm/(hm2·h). The rainfall erosivity value is the smallest in most areas of southwest of the basin, and the average rainfall erosivity is concentrated in 14,674–15,002 MJ·mm/(hm2·h). It can be seen from Figure 4b that the value of multi-year average rainfall erosivity is relatively consistent with the spatial pattern of multi-year average rainfall. The linear regression method is adopted to establish the correlation between the two (Figure 5); the results show that the correlation coefficient R of the multi-year average rainfall erosivity and the multi-year average rainfall is 0.974, and they are extremely significantly correlated.

3.2. Analysis of Land Use Change

IDRISI 17.0 Selva and ArcGIS 10.2 software developed by relevant departments in the United States were used, the land use data of the basin in 2010, 2015 and 2018 were statistically analyzed, such as Table 2; at the same time, the land use area transfer matrix from 2010 to 2015 and from 2015 to 2018 was produced, as shown in Table 3 and Table 4. It can be seen from Table 2 that in the three calculation years, the land use types are mainly garden land, accounting for more than 80% of the total basin area, while the grassland use type is relatively small, accounting for less than 1.5% of the total basin area. The order of land use types by area is Garden land > Cultivated land > Forest land > Commercial land > Grassland. As can be seen from Table 3, the grassland area remained unchanged and the garden land area decreased slightly in 2010 to 2015, mainly converted to forest land, commercial land and arable land, with forest land being the largest. As can be seen from Table 4, more complex two-way conversions occurred in all land use types from 2015 to 2018, with the area of grassland, forest land and garden land decreased by 6.144 km2, 9.198 km2 and 5.145 km2, respectively; the reduced area is mainly converted into cultivated land and commercial land.

3.3. Impact of Land Use Change on Rainfall Erosivity

The SWAT (Soil and Water Assessment Tool) plug-in in ArcGIS software was used to divide the basin into 27 sub-basins based on the digital elevation model. Combined with the rainfall erosivity model, the rainfall erosivity values of each sub-basin under land use scenarios in 2010, 2015 and 2018 were calculated successively (Figure 6).
The comparison shows that under the land use scenario in 2018, the rainfall erosion potential of the whole basin is minimal, with the rainfall erosivity value of each sub-basin is between 11,773–12,973 MJ·mm/(hm2·h). Under the land use scenario in 2010, the rainfall erosivity of the whole basin was relatively the largest, with the rainfall erosivity value of each sub-basin was between 17,889–18,209 MJ·mm/(hm2·h). Under the three land use conditions, the rainfall erosivity potential of each sub-basin in the middle and south is greater than that in the north of the basin, which is mainly because: the cultivated land use types in the north and northeast of the basin occupy more areas, while the grassland is mainly concentrated in the middle and south of the basin. Grassland can effectively reduce the direct effect of raindrops on soil particles, and has a stronger retention effect on surface runoff, which reduces the fluidity of soil particles. Moreover, according to the actual cultivated land types in the basin, most of them are low vegetation cover crops such as rice, peanut and vegetables; their interception effect on raindrops is not as good as that of grassland, so the potential of forming rainfall erosivity is greater.
Based on the distribution maps of rainfall erosivity in each year, the distribution maps of rainfall erosivity change rate from 2010 to 2015 and from 2015 to 2018 can be calculated (Figure 7). Compared with the land use scenario in 2010, the potential of rainfall erosion in the basin showed a decreasing trend in 2015, but the decrease was not significant, and the change rate was between 0.031 and 0.034. Among them, the decline of No.4, No.11, No.12 in the northeast and No.5, No.13 in the west of the basin is relatively obvious, with a change rate more than 0.032. Compared with the land use scenario in 2015, the potential of rainfall erosion in each sub-basin also showed a downward trend in 2018, and the change rate increased from the northern part of the basin to the surrounding areas, with a relatively large change range between 0.091 and 0.176. Among them, the southern and southwestern sub-basins of No. 17, No. 21, No. 26, have the most significant decline, with a change rate of more than 0.150.
The sub-basins No. 3, No. 17, No. 21 and No. 26, with the most significant rainfall erosion rate during 2010–2015 and 2015–2018, were selected as the main objects for land use change analysis (Table 5). According to the comparison of land use transfer matrix in the basin, compared with the land use conditions in 2010, the area of garden land in the basin decreased in 2015, while the area of forest land and commercial land increased, which led to a downward trend of rainfall erosivity in 2015. In 2018, compared with 2015, the area of grassland, forest and garden land decreased and was mainly converted to commercial land, followed by cultivated land. Although the increase in cultivated land area will promote the increase of rainfall erosivity to a certain extent, the large increase in commercial land area is the main reason for the decrease in rainfall erosivity potential under the land use scenario in 2018.
Through the analysis of sub-basins No. 13, No. 17, No. 21, No. 26, it is found that the conversion of garden land to commercial land makes the area of regional hardening increase greatly. Although the increase in commercial land area easily produces runoff and form surface runoff, the hardening layer protects the soil surface and avoids the direct exposure of soil particles to splash erosion of raindrops, thereby reducing the likelihood of rainfall erosion. The conversion of garden land to woodland land increases the vegetation coverage in the watershed, and the interception effect on rainfall is more significant, which weakens the kinetic energy of rainfall, thus reducing the erosion potential of rainfall and reducing the possibility of soil erosion. The conversion of grassland, forest land and garden land into cultivated land will reduce the vegetation coverage of the watershed, enhance the splash erosion of raindrops on soil, and increase the turbulence of surface runoff and rainfall erosion force. On the other hand, under the type of cultivated land use, the underground root amount and root distribution in the basin will also be reduced accordingly, the soil structure is relatively single, the porosity becomes smaller, and the infiltration performance is weakened. When encountering heavy rainfall, it is easy to form super-infiltrating runoff, resulting in serious soil erosion and water loss of cultivated land, and the potential of rainfall erosion also increases correspondingly.
Although the conversion area of the above land use types is small, it provides reference for the optimal allocation of land resources and the reduction of soil loss caused by land use conversion in the basin.

4. Discussion

4.1. Effects of Land Use Change on Rainfall Erosivity

The process of soil erosion includes the comprehensive effect of rainfall, vegetation, soil and human factors, among which rainfall and vegetation are the two main factors affecting the degree of soil erosion, and other factors can be reflected in vegetation, while rainfall erosivity is the main index reflecting erosive rainfall. Based on the meteorological data from 1980 to 2019 and the land use data of 2010, 2015 and 2018 in the Luojiang River Basin, combined with the rainfall erosivity model, this paper compared the response of rainfall erosivity in each sub-basin under different land use change conditions, and concluded that the main reason for the decreasing trend of rainfall erosion potential in 2010, 2015 and 2018 is the decrease of rainfall, but land use change is also an important influencing factor. As the annual rainfall decreases, the rainstorm intensity and surface runoff decrease correspondingly, so the value of rainfall erosivity decreases and the potential of rainfall erosion decreases. On the one hand, land use change can change the watershed vegetation coverage; on the other hand, it can affect the rate and magnitude of runoff production and confluence over the underlying surface, and then affect the magnitude of rainfall erosion potential.
However, this study only explores the relationship between land use change and rainfall erosivity in the Luojiang River Basin in general, and then it provides a theoretical basis for the management of soil erosion. The difference of rainfall erosivity caused by the change of vegetation type, vegetation coverage and farming pattern is not much involved, and more attention should be paid to this aspect in the future. In addition, it should be noted that the spatial distribution characteristics of rainfall and rainfall erosivity values in the paper were obtained by Kriging interpolation method, which may deviate from the measured values in different regions.

4.2. Recommendations

The potential of rainfall erosion is different with different land use patterns. The increase in rainfall erosion potential will undoubtedly lead to soil erosion and other problems. Only by optimizing the land use structure, selecting a reasonable way of land use, efficiently utilizing limited land resources and controlling the potential of rainfall erosion can we realize the sustainable development of agriculture. According to the research experience of land use and rainfall erosion at home and abroad, combined with the characteristics of land use change and the change trend of rainfall erosivity in the basin, it is suggested that the watershed should develop crops with high vegetation coverage in the central and northern regions where the potential of rainfall erosion is relatively large, and establish a multi-level rainfall erosion prediction model, carry out the quantitative evaluation of the relationship between land use dynamic change and rainfall erosion, and construct a distributed land use and rainfall erosion monitoring network based on modern space technology.

5. Conclusions

This paper is based on the simple model of daily rainfall erosivity, combined with Mann–Kendall trend analysis and Kriging interpolation method to analyze the temporal and spatial distribution characteristics of rainfall erosivity in the Luojiang Basin. On this basis, by comparing the response relationships of rainfall erosivity values of each sub-basin under different land use scenarios in 2010, 2015 and 2018, the following conclusions are drawn:
(1)
From 1980 to 2019, the distribution range of annual rainfall erosivity in the Luojiang River Basin was 14,674–15,227 MJ·mm/(hm2·h), and the time series showed an obvious increasing trend, with the growth rate of 29.048 MJ·mm/(hm2·h·a), which was significantly correlated with annual rainfall (correlation coefficient R = 0.974).
(2)
The average value of rainfall erosivity is 14,102 MJ·mm/(hm2·h), and the standard deviation is 4325.206 MJ·mm/(hm2·h). The spatial difference is large, and the overall trend is decreasing from the middle to the surroundings of the basin.
(3)
Land use change is an important factor affecting the spatial and temporal distribution characteristics of rainfall erosivity in the basin. It is suggested to improve the vegetation coverage in the middle and north of the basin and establish a multi-level rainfall erosion prediction model.

Author Contributions

J.H. was responsible for the original concept and writing the paper. Y.-R.W. and H.-T.C. processed the data and conducted the program design. S.-L.W. revised the manuscript and shared numerous comments and suggestions to improve the study quality. All authors have read and agreed to the published version of the manuscript.

Funding

The authors are grateful to the project of key science and technology of the Henan province (No: 222102320333).

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Not Applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, B.; Tao, H.; Song, C.; Guo, B.; Shi, Z.; Zhang, C.; Kong, B.; He, B. Temporal and spatial variations of rainfall erosivity in China during 1960–2009. Geogr. Res. 2013, 32, 245–256. [Google Scholar] [CrossRef]
  2. Yin, S.; Xue, X.; Yue, T.; Xie, Y.; Gao, G. Spatial-temporal distribution and return period of rainfall erosivity in China. Trans. CSAE 2019, 35, 105–113. [Google Scholar] [CrossRef]
  3. Wischmeier, W.H.; Smith, D.D. Rainfall energy and its relationship to soil loss. Eos Trans. Am. Geophys. Union 1958, 39, 285–291. [Google Scholar] [CrossRef]
  4. Wischmeier, W.H. A Rainfall Erosion Index for a Universal Soil-Loss Equation. Soil Sci. Soc. Am. J. 1959, 23, 246–249. [Google Scholar] [CrossRef]
  5. Hudson, N. Soil Conservation; Cornell University Press: Ithaca, NY, USA, 1971. [Google Scholar]
  6. Wischmeier, W.H.; Smith, D.D. Predicting Rainfall Erosion Losses: A Guide to Conservation Planning; Agriculture Handbook No. 537; U.S. Department of Agriculture: Washington, DC, USA, 1978. [Google Scholar]
  7. Wei, O.; Wu, Y.; Hao, Z.; Zhang, Q.; Bu, Q.; Gao, X. Combined impacts of land use and soil property changes on soil erosion in a mollisol area under long-term agricultural development. Sci. Total Environ. 2018, 613–614, 798–809. [Google Scholar] [CrossRef]
  8. Duulatov, E.; Chen, X.; Amanambu, A.C.; Ochege, F.U.; Orozbaev, R.; Issanova, G.; Omurakunova, G. Projected Rainfall Erosivity Over Central Asia Based on CMIP5 Climate Models. Water 2019, 11, 897. [Google Scholar] [CrossRef] [Green Version]
  9. Renard, K.G.; Foster, G.R.; Weesies, G.A.; McCool, D.K.; Yoder, D.C. Predicting Soil Erosion by Water: A Guide to Conservation Planning with the Revised Universal Soil Loss Equation (RUSLE); United States Department of Agriculture: Washington, DC, USA, 1997. [Google Scholar]
  10. Gu, C.; Mu, X.; Gao, P.; Zhao, G.; Sun, W.; Yu, Q. Rainfall erosivity and sediment load over the Poyang Lake Basin under variable climate and human activities since the 1960s. Theor. Appl. Climatol. 2019, 136, 15–30. [Google Scholar] [CrossRef]
  11. Li, Y.; Liu, H. Temporal and spatial distribution characteristics of rainfall erosivity in Pingjiang River Basin of upper Ganjiang River. Bull. Soil Water Conserv. 2020, 40, 23. [Google Scholar] [CrossRef]
  12. Gu, Z.; Feng, D.; Duan, X.; Gong, K.; Li, Y.; Yue, T. Spatial and Temporal Patterns of Rainfall Erosivity in the Tibetan Plateau. Water 2020, 12, 200. [Google Scholar] [CrossRef] [Green Version]
  13. Richardson, C.W.; Foster, G.R.; Wright, D.A. Estimation of Erosion Index from Daily Rainfall Amount. Trans. ASABE 1983, 26, 153–156. [Google Scholar] [CrossRef]
  14. Yu, B.; Rosewell, C.J. An assessment of a daily rainfall erosivity model for New South Wales. Soil Res. 1996, 34, 139–152. [Google Scholar] [CrossRef]
  15. Jia, Z.; Wang, X.; Li, J.; Zeng, B. Study on Rainfall Erosion in Loess Plateau of Western Shanxi. Soil Water Conserv. China 1991, 1, 45–48, 66. [Google Scholar] [CrossRef]
  16. Jiang, Y.; Wang, W.; Wang, Z.; Hu, G.; Lei, H. The Soil Erosion in the North of Shaanxi Province of the Loess Plateau and Its Synthetizing Harness. Sresearch Soil Water Conserv. 1999, 6, 174–180. [Google Scholar] [CrossRef]
  17. Huang, Y.; Lu, C.; Zheng, T.; Fu, Q.; Xu, J. Study on R value of rainfall erosivity index in Southeast Fujian. J. Soil Water Conserv. 1992, 4, 1–5. [Google Scholar]
  18. Zhang, X.; Xu, J.; Lu, X.; Deng, Y.; Gong, D. Study on Soil Loss Equation in Heilongjiang Province. Bull. Soil Water Conserv. 1992, 4, 18. [Google Scholar] [CrossRef]
  19. Zhang, W.; Xie, Y.; Liu, B. Study on Calculation of Rainfall Erosion by Daily Rainfall. Geogr. Sci. 2002, 6, 705–711. [Google Scholar] [CrossRef]
  20. Zhang, W. Spatial and temporal distribution of rainfall erosivity in agro-pastoral ecotone of North China. Prog. Nat. Sci. 2003, 13, 93–96. [Google Scholar] [CrossRef]
  21. Gao, S.; Han, Z.; Yin, S.; Huang, D.; Wang, W. Characteristics of rainfall erosivity and prediction of future changes in the Yellow River Basin from 1961 to 2017. J. Basic Sci. Eng. 2021, 29, 575–590. [Google Scholar] [CrossRef]
  22. Wang, Y.; Tan, S.; Liu, B.; Yang, Y. Estimating rainfall erosivity by incorporating seasonal variations in parameters into the Richardson model. J. Geogr. Sci. 2017, 27, 275–296. [Google Scholar] [CrossRef] [Green Version]
  23. Benavidez, R.; Jackson, B.; Maxwell, D.; Norton, K. A review of the (Revised) Universal Soil Loss Equation ((R)USLE): With a view to increasing its global applicability and improving soil loss estimates. Hydrol. Earth Syst. Sci. 2018, 22, 6059–6086. [Google Scholar] [CrossRef] [Green Version]
  24. Sun, W.; Shao, Q.; Liu, J. Soil erosion and its response to the changes of precipitation and vegetation cover on the Loess Plateau. J. Geogr. Sci. 2013, 23, 1091–1106. [Google Scholar] [CrossRef]
  25. Panagos, P.; Borrelli, P.; Poesen, J.; Ballabio, C.; Lugato, E.; Meusburger, K.; Montanarella, L.; Alewell, C. The new assessment of soil loss by water erosion in Europe. Environ. Sci. Policy 2015, 54, 438–447. [Google Scholar] [CrossRef]
  26. Alatorre, L.C.; Beguería, S.; Lana-Renault, N.; Navas, A.; García-Ruiz, J.M. Soil erosion and sediment delivery in a mountain catchment under scenarios of land use change using a spatially distributed numerical model. Hydrol. Earth Syst. Sci. 2012, 16, 1321–1334. [Google Scholar] [CrossRef] [Green Version]
  27. Mohammad, A.G.; Adam, M.A. The impact of vegetative cover type on runoff and soil erosion under different land uses. CATENA 2010, 81, 97–103. [Google Scholar] [CrossRef]
  28. Bakker, M.M.; Govers, G.; Kosmas, C.; Vanacker, V.; Oost, K.v.; Rounsevell, M. Soil erosion as a driver of land-use change. Agric. Ecosyst. Environ. 2005, 105, 467–481. [Google Scholar] [CrossRef]
  29. Gyssels, G.; Poesen, J.; Bochet, E.; Li, Y.J.P.i.p.g. Impact of plant roots on the resistance of soils to erosion by water: A review. Prog. Phys. Geogr. Earth Environ. 2005, 29, 189–217. [Google Scholar] [CrossRef] [Green Version]
  30. Xu, J. Precipitation–vegetation coupling and its influence on erosion on the Loess Plateau, China. CATENA 2005, 64, 103–116. [Google Scholar] [CrossRef]
  31. Vijith, H.; Hurmain, A.; Dodge-Wan, D. Impacts of land use changes and land cover alteration on soil erosion rates and vulnerability of tropical mountain ranges in Borneo. Remote Sens. Appl. Soc. Environ. 2018, 12, 57–69. [Google Scholar] [CrossRef]
  32. Li, H.; Yan, F.; Jia, J.; Tang, D.; Zhang, Y. Soil Water Availability and Holding Capacity of Different Vegetation Types in Hilly-Gullied Region of the Loess Plateau. Acta Ecol. Sin. 2018, 38, 3889–3898. [Google Scholar] [CrossRef]
  33. Liu, C.; Qi, S.; Shi, M. Research Progress on Relationship between Land Use Change and Soil Erosion. J. Soil Water Conserv. 2001, 6, 10–13, 17. [Google Scholar]
  34. Dai, C.; Liu, Y.; Wang, T.; Li, Z.; Zhou, Y. Exploring optimal measures to reduce soil erosion and nutrient losses in southern China. Agric. Water Manag. 2018, 210, 41–48. [Google Scholar] [CrossRef]
  35. Jiao, J.; Wang, Z.; Wei, Y.; Su, Y.; Cao, B.; Li, Y. Characteristics of erosion sediment yield with extreme rainstorms in Yanhe Watershed Based on field measurement. Trans. Chin. Soc. Agric. Eng. 2017, 33, 159–167. [Google Scholar] [CrossRef]
  36. Huang, C.; Yang, Q.; Cao, X.; Li, Y. Assessment of the Soil Erosion Response to Land Use and Slope in the Loess Plateau—A Case Study of Jiuyuangou. Water 2020, 12, 529. [Google Scholar] [CrossRef] [Green Version]
  37. Sharma, A.; Tiwari, K.N.; Bhadoria, P.B.S. Effect of land use land cover change on soil erosion potential in an agricultural watershed. Environ. Monit. Assess. 2011, 173, 789–801. [Google Scholar] [CrossRef]
  38. Wei, W.; Chen, L.; Fu, B.; Chen, J. Water erosion response to rainfall and land use in different drought-level years in a loess hilly area of China. CATENA 2010, 81, 24–31. [Google Scholar] [CrossRef]
  39. Chen, H.; Zhang, X.; Abla, M.; Lü, D.; Yan, R.; Ren, Q.; Ren, Z.; Yang, Y.; Zhao, W.; Lin, P.; et al. Effects of vegetation and rainfall types on surface runoff and soil erosion on steep slopes on the Loess Plateau, China. CATENA 2018, 170, 141–149. [Google Scholar] [CrossRef]
  40. Cerdan, O.; Bissonnais, Y.L.; Souchère, V.; Martin, P.; Lecomte, V. Sediment concentration in interrill flow: Interactions between soil surface conditions, vegetation and rainfall. Earth Surf. Processes Landf. 2002, 27, 193–205. [Google Scholar] [CrossRef]
  41. Dai, W.; Liu, Q.; Zhao, Y.; Li, W.; Zhou, B.; Dai, J.; Cheng, J.; Deng, N. Utilization and quality study of new cultivated land in Huazhou. Guangdong Agric. Sci. 2013, 40, 193–196. [Google Scholar] [CrossRef]
  42. He, J.; Wan, Y.-R.; Chen, H.-T.; Wang, W.-C. Study on the Impact of Land-Use Change on Runoff Variation Trend in Luojiang River Basin, China. Water 2021, 13, 3282. [Google Scholar] [CrossRef]
  43. Patakamuri, S.K.; Muthiah, K.; Sridhar, V. Long-Term Homogeneity, Trend, and Change-Point Analysis of Rainfall in the Arid District of Ananthapuramu, Andhra Pradesh State, India. Water 2020, 12, 211. [Google Scholar] [CrossRef] [Green Version]
  44. Cressie, N. Spatial prediction and ordinary kriging. Math. Geol. 1988, 20, 405–421. [Google Scholar] [CrossRef]
  45. Yang, C.; Wei, T.; Li, Y. Simulation and Spatio-Temporal Variation Characteristics of LULC in the Context of Urbanization Construction and Ecological Restoration in the Yellow River Basin. Sustainability 2022, 14, 789. [Google Scholar] [CrossRef]
  46. Dile, Y.T.; Daggupati, P.; George, C.; Srinivasan, R.; Arnold, J. Introducing a new open source GIS user interface for the SWAT model. Environ. Model. Softw. 2016, 85, 129–138. [Google Scholar] [CrossRef]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Sustainability 14 08441 g001
Figure 2. Land use type data in 2010 (a), 2015 (b), and 2018 (c).
Figure 2. Land use type data in 2010 (a), 2015 (b), and 2018 (c).
Sustainability 14 08441 g002
Figure 3. Trend analysis of rainfall erosivity.
Figure 3. Trend analysis of rainfall erosivity.
Sustainability 14 08441 g003
Figure 4. Annual average rainfall erosivity and annual average rainfall distribution map. (a) Annual average rainfall erosivity map. (b) Annual average rainfall distribution map.
Figure 4. Annual average rainfall erosivity and annual average rainfall distribution map. (a) Annual average rainfall erosivity map. (b) Annual average rainfall distribution map.
Sustainability 14 08441 g004
Figure 5. Correlation between multi-year average rainfall erosivity and multi-year average rainfall.
Figure 5. Correlation between multi-year average rainfall erosivity and multi-year average rainfall.
Sustainability 14 08441 g005
Figure 6. Rainfall erosivity distribution in each sub-basins.
Figure 6. Rainfall erosivity distribution in each sub-basins.
Sustainability 14 08441 g006
Figure 7. (a) Change rate of rainfall erosivity from 2010 to 2015; (b) Change rate of rainfall erosivity from 2015 to 2018.
Figure 7. (a) Change rate of rainfall erosivity from 2010 to 2015; (b) Change rate of rainfall erosivity from 2015 to 2018.
Sustainability 14 08441 g007
Table 1. Data description and Preprocessing.
Table 1. Data description and Preprocessing.
Data TypeData DescriptionData SourceData Processing
Land UseLand use type maps of 2010, 2015 and 2018 with a resolution of 1 km × 1 kmData Center for Resources and Environmental Sciences, Chinese Academy of Sciences http://www.resdc.cn (accessed on 30 March 2021)Projection conversion, splicing, reclassification, calibration with actual data
Meteorological dataDaily precipitation at Jialong, Baowei and Hejiang stations from 1980 to 2019Maoming Hydrology Sub-bureau of Guangdong Hydrology BureauData extraction, missing data filling, correlation analysis
Table 2. Statistical table of land-use types in different years.
Table 2. Statistical table of land-use types in different years.
Land Type201020152018
Area/km2Percentage/%Area/km2Percentage/%Area/km2Percentage/%
Arable land272.19513.40272.21213.40278.82113.77
Garden land1643.11280.881627.74880.121622.60380.10
Forest land53.9652.6662.3423.0753.1442.62
Grassland23.2131.1423.2131.1417.0690.84
Commercial land39.1261.9346.0972.2755.0602.72
Table 3. 2010–2015 Land use transfer matrix.
Table 3. 2010–2015 Land use transfer matrix.
Land Transfer2015
GrasslandArable LandForest LandCommercial LandGarden LandSum /km2
2010Grassland23.213 23.213
Arable land 272.1860.001 0.008272.195
Forest land 53.965 53.965
Commercial land 39.126 39.126
Garden land 0.0268.3756.9711627.741643.112
Sum/km223.213272.21262.34246.0971627.7482031.611
Table 4. 2015–2018 Land use transfer matrix.
Table 4. 2015–2018 Land use transfer matrix.
Land Transfer2018
GrasslandArable LandForest LandCommercial LandGarden LandSum /km2
2015Grassland2.974.372 15.85723.199
Arable land3.47458.425.34115.495189.129271.859
Forest land0.3066.1913.3971.09240.50861.493
Commercial land1.1714.6910.0943.91826.09145.964
Garden land9.149195.14834.31234.5551349.441622.603
Sum/km217.069278.82153.14455.0601621.0252025.119
Table 5. Proportion of land use type area in sub-basins.
Table 5. Proportion of land use type area in sub-basins.
Proportion of Land Use Type (%)Sub-Basin No. 13Sub-Basin No. 17
Arable LandGarden LandForest LandGrasslandCommercial LandArable LandGarden LandForest LandGrasslandCommercial Land
20100.293.690.010.050.000.674.030.430.230.00
20150.293.590.090.050.030.673.900.510.230.05
20180.513.330.160.000.050.893.940.200.100.19
2010—20150.00−0.110.080.000.030.00−0.130.080.000.05
2015—20180.22−0.250.07−0.050.020.220.03−0.31−0.130.14
Proportion of Land Use Type (%)Sub-Basin No. 21Sub-Basin No. 26
Arable LandGarden LandForest LandGrasslandCommercial LandArable LandGarden LandForest LandGrasslandCommercial Land
20100.521.330.150.050.000.361.070.000.000.08
20150.521.330.150.050.000.361.070.000.000.08
20180.411.660.100.000.010.350.990.000.000.16
2010—20150.000.000.000.000.000.000.000.000.000.00
2015—2018−0.110.33−0.06−0.050.01−0.01−0.080.000.000.08
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

He, J.; Wan, Y.-R.; Chen, H.-T.; Wang, S.-L. Effects of Land Use Change on Rainfall Erosion in Luojiang River Basin, China. Sustainability 2022, 14, 8441. https://doi.org/10.3390/su14148441

AMA Style

He J, Wan Y-R, Chen H-T, Wang S-L. Effects of Land Use Change on Rainfall Erosion in Luojiang River Basin, China. Sustainability. 2022; 14(14):8441. https://doi.org/10.3390/su14148441

Chicago/Turabian Style

He, Ji, Yu-Rong Wan, Hai-Tao Chen, and Song-Lin Wang. 2022. "Effects of Land Use Change on Rainfall Erosion in Luojiang River Basin, China" Sustainability 14, no. 14: 8441. https://doi.org/10.3390/su14148441

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