Next Article in Journal
Using the OKR Method and Fuzzy Logic to Determine the Level of Sustainability in Restaurants
Previous Article in Journal
Feasibility of Efficient, Direct, Butanol Production from Food Waste without Nutrient Supplement by Clostridium saccharoperbutylacetonicum N1-4
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variation Characteristics of Two Erosion Forces and Their Potential Risk Assessment in the Pisha Sandstone Area

1
State Key Laboratory of Eco–Hydraulics in Northwest Arid Region, Xi’an University of Technology, Xi’an 710048, China
2
Key Laboratory of National Forestry Administration on Ecological Hydrology and Disaster Prevention in Arid Regions, National Forestry and Grassland Administration, Xi’an 710048, China
3
Shaanxi Institute of Soil and Water Conservation Survey and Planning, Xi’an 710004, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(7), 6064; https://doi.org/10.3390/su15076064
Submission received: 20 February 2023 / Revised: 24 March 2023 / Accepted: 28 March 2023 / Published: 31 March 2023
(This article belongs to the Section Soil Conservation and Sustainability)

Abstract

:
Precipitation and wind, as the main external erosion forces in wind–water erosion crisscross regions, have profound impacts on water and soil loss. Meanwhile, with the intensification of climate change and human activities, the variation characteristics and risks caused by erosion forces need to be reassessed. In this study, we explored the time-varying characteristics, differences in action period and spatial distribution, and temporal evolution of risk for the compound events of two erosion forces, including precipitation and wind, in the Pisha sandstone area, one of the most seriously eroding and difficult-to-control areas in the Loess Plateau. The results indicated that: (1) the stationarity of regional precipitation was not destroyed, but the mean change existed in the five subseries divided by the detected change points in wind; (2) wind acted earlier than precipitation and increased from southeast to northwest, while precipitation did the opposite; and (3) precipitation-led erosion has become the main erosion type in this area. The above results reveal the evolution and dominant types of regional external erosion forces in a changing environment and thus have implications for regional erosion studies and policy adjustments.

1. Introduction

Water and wind erosion are the two main types of erosion that cause the loss of land and water resources and a decline in land production capacity [1,2,3]. In different locations and at different times, the degree of soil and water loss caused by water and wind is also different. Taking the spatial distribution of erosion types in China as an example, water erosion is widespread throughout the country [4,5], while wind erosion occurs mainly in the arid and semi-arid lands of northwestern, northern, and northeastern China [6,7]. The erosion types overlap in some local areas to form erosion crisscross regions in which two or more of these erosion forces interact temporally and spatially. Such areas are prone to erosion, and are also the focus of soil erosion control [8]. Among them, arid and semi-arid regions are prone to wind–water complex erosion, hindering the normal development of local production, human life, and ecology, thus forming wind–water erosion crisscross regions, which are widely distributed and are hotspots for erosion research [9,10,11].
Studies of soil erosion in such areas have tended to focus on two aspects: either adopting one or more models to evaluate individual types of erosion, or conducting specific experiments in a laboratory or the field. Approaches based on the former use of well-developed methods that focus on evaluating the rate of one particular type of erosion using the revised universal soil loss equation (RUSLE), revised wind erosion equation (RWEQ), etc. For instance, Du et al. [12] assessed the local erosion situation during the 2000s using the RWEQ and RUSLE. Meanwhile, approaches based on the latter are always small-scale studies investigating microscopic erosional processes for a given hypothesis. Yang et al. [13] estimated soil loss by wind erosion using beryllium-7 measurements in a wind tunnel experiment. Ren et al. [14] evaluated the effects of sand cover on runoff soil loss by conducting a series of rainfall simulation experiments. Guo et al. [15] measured the compound erosion from 2012 to 2014 in an agricultural field at the northern end of the Loess Plateau (LP). These previous studies have promoted our re-understanding of the distribution and mechanism of erosion in erosion crisscross regions. However, the limited number of experimental scenarios and the time length of model input data made these studies more targeted. In fact, the related experiments on soil erosion often needed to answer a specific question, that is, whether the setting of external forces is reasonable. Usually, researchers use the mean or maximum value of external forces to carry out relevant experiments. That is to say, it is very important to zonal erosion studies whether the time series of external forces are stationarity. Meanwhile, analyses on the temporal and spatial features of external forces urgently need to transition from static to dynamic analyses under the intensification of climate change and human activities, especially in intersecting erosion regions with severe erosion conditions where the vegetation habitat is no longer sufficiently supported [16,17].
After a long period of systematic management, the overall situation of soil erosion in the Loess Plateau, belonging to the main sand source of the Yellow River, has improved [12,18,19]. In the future, the focus of soil erosion control will shift to the ecologically weak areas in the region. Among them, the Pisha sandstone area, belonging to a typical wind–water erosion crisscross region, is the most seriously eroding and difficult-to-control area in the LP and is affected by the low strength of the rock structure, the interlacing of wind and water erosion in time and space, and the fragile ecological environment [20]. In addition, precipitation and wind are the primary drivers of water and wind erosion, respectively. The variation characteristics between them can be used to answer whether the erosion risks caused by external erosion forces have changed.
Above all, this study investigated the change characteristics of the precipitation and average wind speed in the Pisha sandstone area from 1960 to 2019 and revealed the probability combination change in the two external forces in the time domain. The aims of this paper were to (1) explore whether the stationarity of the two main external forces in the region is broken; (2) reveal the action period and spatial distribution of the above forces based on long-term data; and (3) evaluate the variability and risk associated with the coupled system of multiple erosion forces.

2. Materials and Methods

2.1. Study Area

The Pisha sandstone area (38.6°–40.3° N, 108.8°–111.7° E) is located in the northeastern part of the LP (Figure 1a). It is a typical wind–water erosion crisscross region. The zone covers an area of 1.67 × 104 km2, of which 1.5% is in Shanxi, 28.0% is in Shaanxi, and 70.5% is in Inner Mongolia [20]. The bedrock in this area is Pisha sandstone, a continental clastic rock characterized by a low degree of diagenesis, poor cementation between the sand grains, and low structural strength, with the result being that it is readily stripped by external forces. The terrain in the area has high relief with a large height difference (up to 845 m) and a fragmented topography (Figure 1b). Based on the characteristics of its topsoil and parent material, the area is divided into zone I (sand-covered Pisha sandstone), zone II (soil-covered Pisha sandstone), and zone III (bare Pisha sandstone) [21]. The climate is temperate continental, with the most intense rainfall from June to September. The mean annual precipitation is 394.9 mm and average wind speed is 2.49 m/s (1960–2019 data). Due to the limited habitat in this area, drought-tolerant shrubs (Hippophae rhamnoides Linn. and Caragana korshinskii Kom.) and grasses (Astragalus adsurgens Pall., Stipa capillata Linn., and Leymus chinensis (Trin.) Tzvelev) are the main pioneer plants. Pinus tabuliformis Carr. and Armeniaca sibirica (L.) Lam. are the main trees growing on the gentle slopes in the area.

2.2. Data Requirement and Preprocessing

The precipitation and average wind speed were considered to be the two main external forces causing erosion by water and wind factors. For the analysis, data for daily precipitation and average wind speed monitored by the 12 national surface meteorological stations (Figure 1b) were downloaded from the Resource and Environment Science and Data Center [22]. Then, daily precipitation and average wind speed for the period of 1960–2019 were converted into annual and monthly series for further analysis. In these processes, a Voronoi diagram was used to transform point meteorological data into areal data, from which the kriging interpolation results reflected the spatial distribution of the meteorological data in the study area. A digital elevation model (DEM) with 30 m resolution was derived from the ASTER GDEM V3 dataset provided by the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences [23]. A slope with 30 m resolution was calculated from the DEM. Based on a Landsat 30 m remote sensing image, the land use data (30 m resolution) of 1980 and 2018 in the region were interpreted, including six first-class categories of arable land, woodland, grassland, water body, construction land, and unutilized land.

2.3. Data Analysis

2.3.1. Identifying Change Point by the Pettitt Test

The Pettitt test [24,25] is a method of detecting the change point in an observation sequence based on nonparametric statistics. The statistic U t , N is equivalent to a Mann–Whitney U test for the two samples X1, …, Xt and Xt+1, …, XN, which are from the same sample; the function is:
U t , N = U t 1 , N + j = 1 N sgn x t x j , t 2 , N
where U t , N is the sample statistics for t; U t 1 , N is the sample statistics for t–1; and sgn x t x j  is the sign function, defined as:
sgn x t x j = 1 i f   x t x j > 0 0 i f   x t x j = 0 1 i f   x t x j < 0
where for the series, xt is the sample value for t and xj is the sample value for j.
When t is 1, the sample statistic is calculated from:
U 1 , N = j = 1 N sgn x 1 x j
For the null hypothesis of the Pettitt test, no change points exist in the series. When |Ut,N| takes a maximum value, the corresponding Xt is considered to be a possible change point. Ut,N > 0 indicates a shift down in the level from the beginning of the series. Similarly, Ut,N < 0 indicates a shift up in the level. The significance probabilities associated with |Ut,N| are given by:
p = 2 exp 6 U t , N 2 / N 2 + N 3
where p is the significance level; Ut,N is the sample statistics for t; and N is the sample number. Effective mutation points exist in the data for p ≤ 0.5 in this method.

2.3.2. Testing Trend Using Modified Mann–Kendall Test

The trends presented by variables over time are generally divided into three categories: downward trend, upward trend, and no trend. Such trends often need to be explained in conjunction with the significance level, and are further classed as significant and non-significant trends. In this study, the modified Mann–Kendall trend test was selected to test the series trend.
The Mann–Kendall test [26] is a widely used non-parametric test for detecting trends in time series, and is recommended by the World Meteorological Organization. When this method is used to calculate the trends in a time series, it appears that autocorrelated data cause the results to fluctuate. Here, the modified Mann–Kendall trend test, which is suitable for autocorrelated data [27], was adopted.
The test statistic S is calculated from:
S = k = 1 n 1 j = k + 1 n sgn x j x k
where xj and xi are the observed j and k, n > j > k, and:
sgn x j x k = 1 i f   x j x k > 0 0 i f   x j x k = 0 1 i f   x j x k < 0
In the modified Mann−Kendall trend test, the variance is modified as follows:
V ( S ) = V A R ( S ) n / n S = n ( n 1 ) ( 2 n + 5 ) / 18 n / n S
where n is the actual number of observations and n/nS* corrects for autocorrelated data. The function is given as:
n / n S = 1 + 2 / n ( n 1 ) ( n 2 ) i = 1 n 1 ( n i ) ( n i 1 ) ( n i 2 ) ρ S ( i )
where n is the actual number of observations and ρS(i) is the autocorrelation function of the ranks of the observations, which is calculated after subtracting a suitable nonparametric trend estimator. Only significant values of ρS(i) are used in Equation (8).
The statistic Z is calculated from:
Z = S 1 / V ( S ) i f S > 0 0 i f S = 0 S + 1 / V ( S ) i f S < 0
A positive (negative) value of Z indicates that the sequence to be tested has an upward (downward) trend. The statistic Z obeys the standard normal distribution. In the two-tailed test, the null hypothesis is rejected if the statistic |Z| > Z1–α/2. The value of Z1–α/2 is obtained from the standard normal distribution table. The significance level of the test is often 0.001, 0.01, 0.05, or 0.1 (Table 1).

2.3.3. Distribution Fitting of Samples Considering Non-Stationarity

The distribution fitting of samples considering non-stationarity follows these steps: (a) whether the non-stationarity exists in the study sequence is judged by combining the analysis methods of Section 2.3.1 and Section 2.3.2; (b) if the stationarity of the sequence is broken, the sequence is restored; otherwise, the third step is carried out; (c) parameters of alternative distributions (Table 2) are estimated by maximum likelihood; (d) the K–S test is used to test the hypothesis of sample distribution; (e) the optimal distribution of variables is screened based on Akaike information criterion (AIC) and minimum root mean square error (RSME).
The AIC is calculated from:
A I C = 2 k + n ln ( SSR / n )
where AIC is the Akaike information criterion; k is the number of parameters; n is the number of samples; SSR is the sum of squared residuals.
The RSME is given by:
R M S E = i = 1 n F e x i F t x i 2 n
where RMSE is the root mean square error; F e x i is the empirical CDF at xi; F t x i is the theoretical CDF at xi; n is the number of samples.

2.3.4. Coupled Scenarios of Multiple Erosion Forces

According to the analysis methods of Section 2.3.3, the optimal distribution of variables (annual precipitation and average wind) was screened. Then, the theoretical values corresponding to different probabilities (0.25 and 0.75) under the specified distribution could be calculated (Figure 2a) and used as the critical values of the nine coupled scenarios of multiple erosion forces (Figure 2b).

3. Results

3.1. Variation Analysis and Stationarity Correction of Annual External Erosion Forces

The annual variation in the precipitation and average wind in the Pisha sandstone area from 1960 to 2019 (Figure 3) shows that the annual precipitation increased at a rate of 0.43 mm/a, while the annual average wind speed decreased at a rate of 0.13 (m/s)/10a. The highest annual precipitation (656.8 mm) occurred in 1967 and the lowest (167.2 mm) in 1965. The highest annual average wind speed (3.18 m/s) occurred in 1960 and the lowest (2.05 m/s) in 1986. In addition, a nonparametric test (MMK) was used to test the trend of two annual series, which showed a non-significant upward trend (Zmmk = 0.72) in precipitation and a non-significant downward trend (Zmmk = −1.35) in the average wind speed. The change point test results of the two annual series showed that although there were no significant change points in annual precipitation, three change point levels were detected in the annual average wind speed: a first-order change point in 1981 (p < 0.001), second-order points in 1976 (p < 0.05) and 2004 (p < 0.01), and a third-order change point in 2014 (p < 0.05).
The annual average wind speed from 1960 to 2019 was divided into five subsequences according to the detected change points, and then the mean, Cv, and slope of the linear regression were used to compare the differences in the subsequences (Table 3). The results showed weakly variable, non-significant decreases (p > 0.05) in the subseries, so the mean change was considered to be the main internal difference between the subseries. In order to keep the sequence consistent, the annual average wind speed was restored according to the mean difference in the subsequences (Figure 4). Subsequently, trend and change point tests were carried out on the modified series, and the test results failed to pass the significance test of 0.05 reliability, indicating that the series modified by the reduction method met the requirements of stationarity.

3.2. Differences in the Action Period and Spatial Distribution of External Erosion Forces

The intra-annual distribution of precipitation in the Pisha sandstone area from 1960 to 2019 (Figure 5) shows that the regional precipitation is composed of a single-peaked distribution, with June to September as the main peak, accounting for 75.8% of the annual precipitation. According to the detected change points, the intra-annual distribution of the average wind speed in the Pisha sandstone area from 1960 to 2019 was divided into five subsequences. In comparison, we found that the order of the monthly average wind speed in each period was consistent with the annual average wind speed, with March to June as the main peak. The variation coefficients of the monthly precipitation and average wind speed were calculated, and the monthly precipitation showed strong variation with a coefficient of variation of 1.08, while the average wind speed showed moderate variation with a coefficient of variation range of 0.11–0.17.
The multi-year average annual precipitation and average wind speed values of 12 meteorological stations in and around the study area were calculated, and the spatial distribution maps of the monitoring values (Figure 6) were obtained by kriging interpolation. At the spatial scale, the multi-year average annual precipitation in the basin ranged from 356.3 to 436.2 mm, with an obvious distribution feature of decreasing from southeast to northwest. Meanwhile, the multi-year average annual average wind speed in the basin ranged from 1.59 to 3.06 m/s, with an obvious distribution feature of increasing from southeast to northwest. The spatial distribution of the multi-year average annual precipitation and average wind speed were divided into seven classes by the natural breaks, and the distribution of the data in each class was counted (Figure 6). The results show that the spatial distribution of the multi-year average annual precipitation was more uneven than that of the multi-year average annual average wind speed, and the variation coefficients of the multi-year average annual precipitation and average wind speed were 0.55 and 0.19, respectively. The main distribution interval of the multi-year average annual precipitation was 395.5–403.6 mm, and that of the multi-year average annual average wind speed was 2.10–2.27 m/s.

3.3. Temporal Evolution of the Risk of Compound Events Consisting of Multiple External Erosion Forces

In this paper, six sampling distribution functions—GAM, EXP, WBL, GEV, GP, and P–III—were used to fit the annual precipitation and average wind speed, and then the optimal distribution of the variables was optimized (Table 4). The results of the distribution optimization show that the annual precipitation followed the P–III distribution and the annual average wind speed followed the GEV distribution. Meanwhile, the PP plots of the above elements (Figure 7) show that the empirical and theoretical CDF of each element were well distributed around the 1:1 line. Then, theoretical values under 0.25, 0.5, and 0.75 probabilities were calculated based on the optimal distribution of the selected forces (Table 5).
Based on the theoretical values under the probabilities of 0.25 and 0.75, nine combination scenarios were constructed. The frequency of different combination scenarios (Table 6) was calculated, and the results show eight actual occurrence types of regional scenarios in the whole period, with five actual occurrence categories of regional scenarios in the 1960s, six in the 1970s, and a sharp reduction to three after the 1980s due to the stepwise variation of the actual monitored wind speed, which was mainly concentrated in the interval of [0, 0.25]. At the same time, when the annual precipitation or average wind speed fell into the interval of (0.75, 1], it was considered that the risk of soil and water loss caused by the above regional erosion force was greater than that of the general scenario, and this broad category included four types (P1W3, P2W3, P3W1, and P3W3) in the whole period. The frequency ranking of the above scenarios shows that P3W1 was the main scenario, with 13 occurrences, and that this type of event had the highest probability of occurring in the 2010s, at 40%. Therefore, regional precipitation-led erosion is still the main target for regional erosion control.

4. Discussion

4.1. Causes of Temporal and Spatial Differences between Two External Erosion Forces

It was found that the external forces causing erosion by wind in the western Pisha sandstone area was more intense than in the east, and that water erosion force was most intense in the east. The northwesterly monsoon is affected by the Siberian high-pressure system (the ‘Siberian High’) in the winter, while the southeasterly monsoon is affected by the Western Pacific Subtropical High in the summer, influencing the wind speed and precipitation distribution on a large scale [28,29,30]. The elevation of the land gradually increases from east to west in this area, which indirectly affects the local meteorological factors [31]. The specific heat capacities of bare and sand-covered land are smaller than for the soil cover, and vegetation on the former type of land are always less prolific than on the latter, which mainly affect the local wind speed [32].
In our research, non-stationarity analysis of the external forces found that only the stationarity of the wind was broken. The average wind speed monitored in 1982–2014 was distinctly smaller than in 2015–2019, which is consistent with a study conducted in northwestern China by Ge et al. [33]. While the precipitation in the Pisha sandstone area showed tiny fluctuations, similar related precipitation results were reported by Zhang et al. [19]. In summary, the wind was more sensitive than the precipitation in this area. The results of the present study are related to the implementation of soil and water conservation measures such as the Three-North Shelter Forest Program implemented in 1978 [34] and the changes in the temperature gradients between middle and high latitudes [35]. “Non-stationarity” is a typical statistical concept that cannot be ignored when carrying out studies involving long-term sequential data, which can be used to prove that under the impact of the above reasons, zonal wind has changed. Hence, when carrying out wind erosion-related experiments, it is necessary to adjust the experimental settings according to the current regional conditions.

4.2. Impacts and Adjustment of the Measures According the Results

Through the analysis of regional land use change (Figure 8), this study found that with the acceleration of urbanization and the implementation of the “Grain for Green Program (GFGP)” in 1999, all kinds of land use in the region have undergone corresponding changes, manifested in the fluctuations in area and the transformation of use types. The proportion of construction land and woodland increased by 4.5% and 0.8%, respectively, while the proportion of other land types decreased by 0.2–1.8% and was less than half of the increase rate of construction land. At the same time, the main land types being turned into grassland were unutilized land and arable land. In the present study, P3W1, a scenario with an annual precipitation that fell into the interval of (0.75, 1] and an annual average wind speed that fell into the interval of [0, 0.25], was the main scenario of the combination scenarios containing soil and water loss risk caused by the related regional erosion forces. Combined with regional land use change, in years with an abundance of water, specific weak links such as regional urban soil erosion and the stability of grassland that belongs to the area changed from unutilized land or farmland still need to be further studied and corresponding measures should be formulated.
In the present study, the two external erosion forces showed obvious spatial differentiation, as shown in Figure 6, which implies that the spatial distribution of the related erosion will be different in this area. Soil erosion is related to external erosion forces, but the response of external forces differs due to underlying surface conditions such as topography, land use, and soil [36,37]. To facilitate further erosion studies in the area, we overlaid the results of the grading of the external forces with the slope of the area (Figure 9a). Then, this study found that the average slope difference in the seven-level classification was the largest, and the slope under this level of annual precipitation was 6.1° larger than that of the annual average wind speed. We argue that the changes in geomorphological features caused by water erosion may be greater than those caused by wind erosion. When further comparing the land use change in the seven-level classification of external forces (Figure 9b,c), we found that the ratio of easily eroded land types such as arable and unutilized land decreased and the ratio of high soil and water conservation functional land types such as forests increased in this class for both external forces. At the same time, the ratio of better soil and water conservation functional land types such as grassland showed inconsistent changes in the seven-level classification of external forces, and the reason for this phenomenon was that the ratio of grassland converted to urban area in the annual average wind speed class 7 was large, and the transfer rate was 45%. To sum up, the erosion in the region undoubtedly decreased, even in areas where the external erosion force was large. Thus, constantly reduced sediment delivery, which can reflect the improvement of regional erosion in local basins, has been reported, although our research indicated that regional precipitation-led erosion has been the main target for regional control [38,39]. By increasing the surface cover, the GFGP has increased soil infiltration and reduced soil erosion [40]. The construction of check dams, one of the main projects of the Ministry of Water Resources of the People’s Republic of China implemented in 2003, has the function of retaining runoff and sediment and controlling gully erosion [41]. Combinations of the above management measures and changes in regional external forces have contributed to enhancing the capabilities of water storage and soil conservation, which have improved the local ecological living environment. However, to prevent potential erosion of the region, the corresponding measure of years with an abundance of water still need to be evaluated in terms of their validity, since P3W1 had been at the highest probability of occurring in the 2010s, at 40%.

5. Conclusions

Clarifying the variability characteristics between the primary drivers of water and wind erosion in time and space can contribute to the improvement of the scientific setting of external erosion forces and the adjustment of relevant erosion control measures during the changing environment affected by intensified climate and human activities. In the present study, the variability characteristics in the time and space of long-term precipitation and wind speed series in a typical wind–water erosion crisscross region located in the northeastern part of the LP were investigated using statistical analyses. In addition, the changing erosion risk caused by the regional erosion forces were revealed using the frequency changes in each decade of different combination scenarios constructed by the theoretical values under 0.25 and 0.75 probabilities, both in precipitation and modified wind speed series.
The key conclusions are summarized as follows: Interannually, the stationarity of the regional average wind speed over the last 60 years was disrupted compared to that of the precipitation, as evidenced by the presence of four significant change points (p < 0.05), namely, 1976, 1987, 2004, and 2014 within the series. Meanwhile, the mean change in the subseries was the main variation type that destroys the stationarity of the annual average wind speed series. Intra-annually, the precipitation and wind speed were both mainly distributed from March to June. Spatially, the precipitation decreased from southeast to northwest, while the wind speed showed the opposite. P3W1 was the main scenario of the combination scenarios containing soil and water loss risk, and had the highest probability of occurring in the 2010s, at 40%. This is why, considering only the primary drivers, regional precipitation-led erosion has become the main risk erosion type.

Author Contributions

Conceptualization, M.X. and Z.R.; data curation, X.M. and R.Y.; methodology, M.X.; software, M.X.; writing—original draft, M.X.; writing—review and editing, Z.R., B.M. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Key Research and Development Program of China (grant number 2022YFF1300803) and the National Natural Science Foundation of China (grant number 42077074).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to extend our thanks to the anonymous reviewers for their suggestions, which helped to improve this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Evans, R.; Boardman, J. The New Assessment of Soil Loss by Water Erosion in Europe. Panagos P. et al., 2015 Environmental Science & Policy 54, 438-447-A Response. Environ. Sci. Policy 2016, 58, 11–15. [Google Scholar]
  2. Mezösi, G.; Blanka, V.; Bata, T.; Kovacs, F.; Meyer, B. Estimation of Regional Differences in Wind Erosion Sensitivity in Hungary. Nat. Hazards Earth Syst. Sci. 2015, 15, 97–107. [Google Scholar] [CrossRef] [Green Version]
  3. Amundson, R.; Berhe, A.A.; Hopmans, J.; Olson, C. Soil and Human Security in the 21st Century. Science 2015, 348, 1261074. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Hou, X.; Shao, J.; Chen, X.; Li, J.; Lu, J. Changes in the Soil Erosion Status in the Middle and Lower Reaches of the Yangtze River Basin from 2001 to 2014 and the Impacts of Erosion on the Water Quality of Lakes and Reservoirs. Int. J. Remote Sens. 2020, 41, 3175–3196. [Google Scholar] [CrossRef]
  5. 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]
  6. Lin, J.; Guan, Q.; Pan, N.; Zhao, R.; Yang, L.; Xu, C. Spatiotemporal Variations and Driving Factors of the Potential Wind Erosion Rate in the Hexi Region, PR China. L. Degrad. Dev. 2021, 32, 139–157. [Google Scholar] [CrossRef]
  7. Shen, Y.; Zhang, C.; Wang, X.; Zou, X.; Kang, L. Statistical Characteristics of Wind Erosion Events in the Erosion Area of Northern China. Catena 2018, 167, 399–410. [Google Scholar] [CrossRef]
  8. Zhang, J.; Yang, M.; Sun, X.; Zhang, F. Estimation of Wind and Water Erosion Based on Slope Aspects in the Crisscross Region of the Chinese Loess Plateau. J. Soils Sediments 2018, 18, 1620–1631. [Google Scholar] [CrossRef]
  9. Ma, X.; Zhao, C.; Zhu, J. Aggravated Risk of Soil Erosion with Global Warming—A Global Meta-Analysis. Catena 2021, 200, 105129. [Google Scholar] [CrossRef]
  10. Jolivet, M.; Braucher, R.; Dovchintseren, D.; Hocquet, S.; Schmitt, J.M. Erosion around a Large-Scale Topographic High in a Semi-Arid Sedimentary Basin: Interactions between Fluvial Erosion, Aeolian Erosion and Aeolian Transport. Geomorphology 2021, 386, 107747. [Google Scholar] [CrossRef]
  11. Fenta, A.A.; Tsunekawa, A.; Haregeweyn, N.; Poesen, J.; Tsubo, M.; Borrelli, P.; Panagos, P.; Vanmaercke, M.; Broeckx, J.; Yasuda, H.; et al. Land Susceptibility to Water and Wind Erosion Risks in the East Africa Region. Sci. Total Environ. 2020, 703, 135016. [Google Scholar] [CrossRef] [PubMed]
  12. Du, H.; Dou, S.; Deng, X.; Xue, X.; Wang, T. Assessment of Wind and Water Erosion Risk in the Watershed of the Ningxia-Inner Mongolia Reach of the Yellow River, China. Ecol. Indic. 2016, 67, 117–131. [Google Scholar] [CrossRef]
  13. Yang, M.Y.; Walling, D.E.; Sun, X.J.; Zhang, F.B.; Zhang, B. A Wind Tunnel Experiment to Explore the Feasibility of Using Beryllium-7 Measurements to Estimate Soil Loss by Wind Erosion. Geochim. Cosmochim. Acta 2013, 114, 81–93. [Google Scholar] [CrossRef]
  14. Ren, Z.; Zhang, X.; Zhang, X.; Li, Z.; Li, P.; Zhou, Z. Sand Cover Enhances Rill Formation under Laboratory Rainfall Simulation. Catena 2021, 205, 105472. [Google Scholar] [CrossRef]
  15. Guo, Q.; Cheng, C.; Jiang, H.; Liu, B.; Wang, Y. Comparative Rates of Wind and Water Erosion on Typical Farmland at the Northern End of the Loess Plateau, China. Geoderma 2019, 352, 104–115. [Google Scholar] [CrossRef]
  16. Borrelli, P.; Robinson, D.A.; Panagos, P.; Lugato, E.; Yang, J.E.; Alewell, C.; Wuepper, D.; Montanarella, L.; Ballabio, C. Land Use and Climate Change Impacts on Global Soil Erosion by Water (2015–2070). Proc. Natl. Acad. Sci. USA 2020, 117, 21994–22001. [Google Scholar] [CrossRef]
  17. Qiu, B.; Chen, G.; Tang, Z.; Lu, D.; Wang, Z.; Chen, C. Assessing the Three-North Shelter Forest Program in China by a Novel Framework for Characterizing Vegetation Changes. ISPRS J. Photogramm. Remote Sens. 2017, 133, 75–88. [Google Scholar] [CrossRef]
  18. Xiao, L.; Li, G.; Zhao, R.; Zhang, L. Effects of Soil Conservation Measures on Wind Erosion Control in China: A Synthesis. Sci. Total Environ. 2021, 778, 146308. [Google Scholar] [CrossRef]
  19. Zhang, X.; Zhang, L.; Zhao, J.; Rustomji, P.; Hairsine, P. Responses of Streamflow to Changes in Climate and Land Use/Cover in the Loess Plateau, China. Water Resour. Res. 2008, 44, W00A07. [Google Scholar] [CrossRef]
  20. Zhang, P.; Yao, W.Y.; Liu, G.B.; Xiao, P.Q.; Sun, W.Y. Experimental Study of Sediment Transport Processes and Size Selectivity of Eroded Sediment on Steep Pisha Sandstone Slopes. Geomorphology 2020, 363, 107211. [Google Scholar] [CrossRef]
  21. Liang, Z.; Wu, Z.; Yao, W.; Noori, M.; Yang, C.; Xiao, P.; Leng, Y.; Deng, L. Pisha Sandstone: Causes, Processes and Erosion Options for Its Control and Prospects. Int. Soil Water Conserv. Res. 2019, 7, 1–8. [Google Scholar] [CrossRef]
  22. Resource and Environment Science and Data Center. Available online: http://www.resdc.cn (accessed on 13 March 2020).
  23. Geospatial Data Cloud Site, Computer Network Information Center, Chinese Academy of Sciences. Available online: http://www.gscloud.cn (accessed on 24 February 2022).
  24. Pettitt, A.N. A Non-Parametric Approach to the Change-Point Problem. Appl. Stat. 1979, 28, 126–135. [Google Scholar] [CrossRef]
  25. Pettitt, A.N. A Simple Cumulative Sum Type Statistic for the Change-Point Problem with Zero-One Observations. Biometrika 1980, 67, 79–84. [Google Scholar] [CrossRef]
  26. Mann, H.B. Non-Parametric Test Against Trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  27. Hamed, K.H.; Rao, A.R. A Modified Mann-Kendall Trend Test for Autocorrelated Data. J. Hydrol. 1998, 204, 182–196. [Google Scholar] [CrossRef]
  28. Ding, T.; Gao, H.; Li, X. Universal Pause of the Human-Perceived Winter Warming in the 21st Century over China. Environ. Res. Lett. 2021, 16, 064070. [Google Scholar] [CrossRef]
  29. Gu, X.; Zhang, Q.; Singh, V.P.; Shi, P. Nonstationarity in Timing of Extreme Precipitation across China and Impact of Tropical Cyclones. Glob. Planet. Change 2017, 149, 153–165. [Google Scholar] [CrossRef] [Green Version]
  30. Zuo, Z.; Zhang, R.; Huang, Y.; Xiao, D.; Guo, D. Extreme Cold and Warm Events over China in Wintertime. Int. J. Climatol. 2015, 35, 3568–3581. [Google Scholar] [CrossRef]
  31. Micu, D.M.; Dumitrescu, A.; Cheval, S.; Nita, I.A.; Birsan, M.V. Temperature Changes and Elevation-Warming Relationships in the Carpathian Mountains. Int. J. Climatol. 2021, 41, 2154–2172. [Google Scholar] [CrossRef]
  32. Duman, T.; Huang, C.W.; Litvak, M.E. Recent Land Cover Changes in the Southwestern US Lead to an Increase in Surface Temperature. Agric. For. Meteorol. 2021, 297, 108246. [Google Scholar] [CrossRef]
  33. Ge, J.; Feng, D.; You, Q.; Zhang, W.; Zhang, Y. Characteristics and Causes of Surface Wind Speed Variations in Northwest China from 1979 to 2019. Atmos. Res. 2021, 254, 105527. [Google Scholar] [CrossRef]
  34. Jiang, B.; Liang, S.; Yuan, W. Observational Evidence for Impacts of Vegetation Change on Local Surface Climate over Northern China Using the Granger Causality Test. J. Geophys. Res. Biogeosciences 2015, 120, 1–12. [Google Scholar] [CrossRef]
  35. Jung, C.; Schindler, D. Changing Wind Speed Distributions under Future Global Climate. Energy Convers. Manag. 2019, 198, 111841. [Google Scholar] [CrossRef]
  36. Brosens, L.; Robinet, J.; Pelckmans, I.; Ameijeiras-Mariño, Y.; Govers, G.; Opfergelt, S.; Minella, J.P.G.; Vanderborght, J. Have Land Use and Land Cover Change Affected Soil Thickness and Weathering Degree in a Subtropical Region in Southern Brazil? Insights from Applied Mid-Infrared Spectroscopy. Catena 2021, 207, 105698. [Google Scholar] [CrossRef]
  37. Ren, Z.; Pan, J.; Li, Z.; Xiao, P.; Shen, Z.; Jia, L.; Li, X. The Interaction of Aeolian Sand and Slope on Runoff and Soil Loss on a Loess Slope via Simulated Rainfall under Laboratory Conditions. Water 2023, 15, 888. [Google Scholar] [CrossRef]
  38. Zheng, H.; Miao, C.; Jiao, J.; Borthwick, A.G.L. Complex Relationships between Water Discharge and Sediment Concentration across the Loess Plateau, China. J. Hydrol. 2021, 596, 126078. [Google Scholar] [CrossRef]
  39. Xie, M.Y.; Ren, Z.P.; Li, Z.B.; Li, P.; Shi, P.; Zhang, X.M. Changes in Runoff and Sediment Load of the Huangfuchuan River Following a Water and Soil Conservation Project. J. Soil Water Conserv. 2020, 75, 590–600. [Google Scholar] [CrossRef]
  40. Yang, S.; Shi, P.; Li, P.; Li, Z.; Niu, H.; Zu, P.; Cui, L. Ecosystem Services Trade-Offs and Synergies Following Vegetation Restoration on the Loess Plateau of China. Sustainability 2023, 15, 229. [Google Scholar] [CrossRef]
  41. Yuan, S.; Li, Z.; Chen, L.; Li, P.; Zhang, Z.; Zhang, J.; Wang, A.; Yu, K. Effects of a Check Dam System on the Runoff Generation and Concentration Processes of a Catchment on the Loess Plateau. Int. Soil Water Conserv. Res. 2022, 10, 86–98. [Google Scholar] [CrossRef]
Figure 1. Geographical location of the Pisha sandstone area at different scales: (a) within the Loess Plateau; (b) inset showing the Pisha sandstone area (I = sand-covered; II = soil-covered; III = bare sandstone).
Figure 1. Geographical location of the Pisha sandstone area at different scales: (a) within the Loess Plateau; (b) inset showing the Pisha sandstone area (I = sand-covered; II = soil-covered; III = bare sandstone).
Sustainability 15 06064 g001
Figure 2. Scenario classification basis: (a) theoretical value selection; (b) coupled scenarios of multiple erosion forces (color-coded scenarios are at higher potential risk of soil and water loss).
Figure 2. Scenario classification basis: (a) theoretical value selection; (b) coupled scenarios of multiple erosion forces (color-coded scenarios are at higher potential risk of soil and water loss).
Sustainability 15 06064 g002
Figure 3. Trend regression test of annual precipitation and average wind speed.
Figure 3. Trend regression test of annual precipitation and average wind speed.
Sustainability 15 06064 g003
Figure 4. Restored annual average wind speed.
Figure 4. Restored annual average wind speed.
Sustainability 15 06064 g004
Figure 5. The intra-annual distribution of the precipitation and average wind speed.
Figure 5. The intra-annual distribution of the precipitation and average wind speed.
Sustainability 15 06064 g005
Figure 6. Spatial distribution of the annual precipitation and average wind speed.
Figure 6. Spatial distribution of the annual precipitation and average wind speed.
Sustainability 15 06064 g006
Figure 7. P–P plot of the optimal distribution function of annual elements.
Figure 7. P–P plot of the optimal distribution function of annual elements.
Sustainability 15 06064 g007
Figure 8. Regional land use change from 1980 to 2018.
Figure 8. Regional land use change from 1980 to 2018.
Sustainability 15 06064 g008
Figure 9. Discussion on risk considering the spatial differentiation of external erosion forces: (a) slope mean in external erosion force classes; (b) land use change in annual precipitation class 7; (c) land use change in annual average wind speed class 7.
Figure 9. Discussion on risk considering the spatial differentiation of external erosion forces: (a) slope mean in external erosion force classes; (b) land use change in annual precipitation class 7; (c) land use change in annual average wind speed class 7.
Sustainability 15 06064 g009
Table 1. Test value at the given significance level.
Table 1. Test value at the given significance level.
Significance LevelTest ValueDegree of Confidence
0.0013.29199.90%
0.012.57699.00%
0.051.9695.00%
0.11.64590.00%
Table 2. Probability functions adopted in this study.
Table 2. Probability functions adopted in this study.
DistributionProbability Density FunctionParameters
Gamma
(GAM)
f ( x α , β ) = 1 β α Γ α x α 1 e x β α, β
exponential
(EXP)
f x μ = 1 μ e x μ μ
Weibull
(WBL)
f x α , β = β α x α β 1 e x α β α, β
generalized extreme value
(GEV)
f x k , μ , σ = 1 σ exp 1 + k x μ σ 1 k 1 + k x μ σ 1 1 k k 0 f x 0 , μ , σ = 1 σ exp exp x μ σ x μ σ k = 0 k, μ, σ
generalized Pareto
(GP)
f x k , σ , θ = 1 σ 1 + k x θ σ 1 1 k k 0 f x 0 , σ , θ = 1 σ e x θ σ k = 0 k, σ, μ
Pearson type III
(P–III)
f x α , β , ξ = 1 β α Γ α x ξ α 1 e x ξ β α = 4 γ 2 β = 1 2 σ γ ξ = μ 2 σ γ μ, σ, γ
Table 3. Internal variability of annual average wind speed.
Table 3. Internal variability of annual average wind speed.
Statistical ParametersPeriod
P1P2P3P4P5
Mean /(m/s)2.962.582.192.302.58
Cv0.050.050.040.040.02
Slope of linear regression (m/s)/a−0.0024−0.0721−0.0041−0.0041−0.0115
Note: the significance of the regression equations was greater than 0.05.
Table 4. Optimization of the distribution of annual elements.
Table 4. Optimization of the distribution of annual elements.
Annual ElementsDistributionK–S TestRSMEAIC
hp_Value
PrecipitationGAM00.820.027−428.4
EXP10.000.227−175.8
WBL00.740.032−409.3
GEV00.930.026−434.5
GP10.000.163−213.4
P–III00.920.025−436.0
Average wind speedGAM00.840.032−408.3
EXP10.000.302−141.8
WBL00.750.033−405.8
GEV00.930.028−424.8
GP10.000.345−123.6
P–III00.920.028−421.5
Table 5. Theoretical values of annual elements.
Table 5. Theoretical values of annual elements.
Cumulative ProbabilityPrecipitation/mmAverage Wind Speed/(m/s)
0.25323.32.89
0.5389.02.96
0.75460.13.03
Table 6. Frequency of the different combination scenarios in each decade.
Table 6. Frequency of the different combination scenarios in each decade.
Sort1960s1970s1980s1990s2000s2010sSum
P1W11234111
P1W222
P1W3112
P2W112755525
P2W233
P2W333
P3W132121413
P3W2
P3W311
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

Xie, M.; Ren, Z.; Miao, B.; Li, Z.; Ma, X.; Yan, R. Variation Characteristics of Two Erosion Forces and Their Potential Risk Assessment in the Pisha Sandstone Area. Sustainability 2023, 15, 6064. https://doi.org/10.3390/su15076064

AMA Style

Xie M, Ren Z, Miao B, Li Z, Ma X, Yan R. Variation Characteristics of Two Erosion Forces and Their Potential Risk Assessment in the Pisha Sandstone Area. Sustainability. 2023; 15(7):6064. https://doi.org/10.3390/su15076064

Chicago/Turabian Style

Xie, Mengyao, Zongping Ren, Binxia Miao, Zhanbin Li, Xiaoni Ma, and Rui Yan. 2023. "Variation Characteristics of Two Erosion Forces and Their Potential Risk Assessment in the Pisha Sandstone Area" Sustainability 15, no. 7: 6064. https://doi.org/10.3390/su15076064

APA Style

Xie, M., Ren, Z., Miao, B., Li, Z., Ma, X., & Yan, R. (2023). Variation Characteristics of Two Erosion Forces and Their Potential Risk Assessment in the Pisha Sandstone Area. Sustainability, 15(7), 6064. https://doi.org/10.3390/su15076064

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