Wind Erosion Changes in a Semi-Arid Sandy Area , Inner Mongolia , China

Wind erosion is one of the major environmental problems in drylands. Identifying the dominant natural factors of wind erosion and using targeted treatment measures are the key steps in wind erosion control. Using Horqin Left Back Banner in China as a case study, we applied the revised wind erosion equation to simulate the spatial distribution of wind erosion in the semi-arid sandy area. Contribution assessment and constraint line analysis were used to investigate the contributions of driving forces to wind erosion changes. The results showed that the wind erosion in the whole area was reduced by 0.35 t/hm2·a from 2005 to 2016. The wind factor and vegetation coverage factor had dominant contributions to the wind erosion modulus and accounted for the erosion in 49.87% and 50.13% of the total area, respectively. In addition, the average wind speed exceeding the threshold and the number of occurrences exhibited significant correlations with the wind erosion severity. Meanwhile, the mitigation effects of vegetation coverage on wind erosion decreased with the increase in wind speed. The temporal mismatch between the wind speed and vegetation coverage was the main reason for the frequent severe wind erosion in spring. Reducing the spring wind speed through adding windbreaks would be an effective method for decreasing wind erosion in semi-arid areas.


Introduction
As a major environmental problem in drylands, wind erosion has drawn considerable attention from researchers [1][2][3].Wind erosion reduces soil fertility [4,5], increases the frequency of sand storms [6], spreads diseases [7,8], and poses serious threats to the health of people and sustainable development.Northern China is one of the most severe wind erosion areas in the world [9].To alleviate wind erosion and enhance the capacity of sand-stabilization, the Chinese government has introduced a series of projects and policies since 2000 [10], including the Three Norths Shelter Forest System Project (Phase IV) [11], the Grain for Green Project [12], and the Natural Forest Conservation Program (Phase II) [13].Although the wind erosion in northern China declined after the implementation of these policies [14], it has still not been fully controlled in some areas, such as the heavily sandy areas in Mu Us Sandy Lands and Horqin Sandy Lands in Inner Mongolia, China.
Wind erosion is generally affected by both natural factors and human activities [15][16][17][18].However, since human activities are bounded by natural conditions and affect wind erosion by changing natural factors [19], the natural factors, such as vegetation coverage, wind speed, temperature, and precipitation, should be the direct contributing drivers to wind erosion.The role of different natural conditions in wind erosion changes needs to be further explored [9,20,21].Previously, researchers have reported significant negative correlations between vegetation coverage and wind erosion [22], effects of temperature change on wind erosion in Inner Mongolia [23], and the ability of wind speed and precipitation to affect wind erosion [14].In addition, studies have also shown that the trend of wind erosion and its dominant factors varied greatly due to regional heterogeneity [23,24].For example, Meng et al. [20] found that the most important natural factors that can influence aeolian sand transportation in Inner Mongolia were wind speed, vegetation coverage, and soil moisture.A finding by Jiang et al. [14] suggested that the spatial distribution pattern of wind erosion was consistent with the distribution patterns of erodible soil, vegetation cover, and precipitation.However, the contribution of each natural factor to regional wind erosion still remains unclear.Therefore, quantitatively identifying the dominant natural factors of wind erosion and using targeted treatment measures are the key steps in wind erosion control.
Various process-based empirical or semi-empirical models were built to assess wind erosion.The most widely used models are the Wind Erosion Equation (WEQ) [25], the Revised Wind Erosion Equation (RWEQ) [26], and the Wind Erosion Prediction System [27].The RWEQ model is refined based on the WEQ model and is widely used in wind erosion simulation and driving forces determination in dryland systems.For instance, Du et al. [28] used field observation data of 2011 and 2012 to calibrate the RWEQ model and to calculate the wind erosion in the Ningxia-Inner Mongolia reach of Yellow River from 1986 to 2013.The calibration results showed satisfactory simulation accuracy.Jiang et al. [14] analyzed the results of the RWEQ by using redundancy analysis and constrained linear ordination, and concluded that wind erosion reduction was related to the increase in precipitation and decrease in the number of days of strong winds in Inner Mongolia.Li et al. [22] found that climate change accounted for almost half of the soil retention increase by using a multiple linear regression model.However, few studies have focused their analysis on the differences in the contribution between different natural factors to wind erosion, which might be the key to adopting targeted wind erosion control measures.
In this study, we used the RWEQ model to simulate wind erosion in the Horqin Left Wing Back Banner, China (HLBB) and used the partial differential method to analyze the contribution of different driving factors to wind erosion changes.Our analysis also focused on the effects of natural factors on the changes in wind erosion at a semi-monthly scale.The main objectives of this study are as follows: (1) to analyze the spatial and temporal changes of wind erosion in HLBB from 2005 to 2016; (2) to investigate the spatial distribution of meteorological and vegetation factors' contributions to the wind erosion dynamics in HLBB; (3) to identify the dominant driving factors that influenced wind erosion changes; (4) to reveal the possible mismatches and differences in the temporal distributions among various natural factors; and (5) to provide a scientific basis for wind erosion control policy making in sandy areas.

Study Areas
The Horqin Left Wing Back Banner (121.5 • E-123.7 • E, 42.7 • N-43.7 • N) is located at the southeast side of the Horqin Sandy Land (Figure 1) and covers an area of 11,476 km 2 .HLBB is a typical semi-arid sand area with dry and windy conditions in spring.The average annual precipitation is approximately 350-450 mm, and approximately 70% of the total annual precipitation is concentrated from July to September [29].The local terrain consists of flat to undulating sand and plains (Figure 1).The topography of HLBB is dominated by sand dunes but also contains alluvial valleys and flat sand.The main soil types include aeolian sandy soil, meadow soil, alkaline soil, and marsh soil, each accounting for approximately 68.9%, 23.1%, 7.2%, and 0.78%, respectively.The surface layer of sandy soil is structureless and is accompanied by coarse sand, which is extremely erosive [30].HLBB is also a typical agro-pastoral area in northern China.The vegetation community is dominated by secondary grassy plant species and sparse shrubs.Trees can only be found in artificially planted economic forests and windbreaks.Grazing and farming are the traditional land use types with distinct boundaries between cultivated land and pasture land (Figure 2).The desertification is more severe in pasture-dominated central and western regions.The areas with relatively low wind erosion are the eastern plains that are dominated by cultivated land and the Daqinggou nature reserve with high forest cover in the southwest region [31].
Sustainability 2018, 10, x FOR PEER REVIEW 3 of 20 boundaries between cultivated land and pasture land (Figure 2).The desertification is more severe in pasture-dominated central and western regions.The areas with relatively low wind erosion are the eastern plains that are dominated by cultivated land and the Daqinggou nature reserve with high forest cover in the southwest region [31].boundaries between cultivated land and pasture land (Figure 2).The desertification is more severe in pasture-dominated central and western regions.The areas with relatively low wind erosion are the eastern plains that are dominated by cultivated land and the Daqinggou nature reserve with high forest cover in the southwest region [31].Since the 1980s, the local government has implemented a series of policies and measures to control desertification, such as the "Two Million Mu" Sand Control Demonstration Project, the "Converting Farmland to Forest Project", and the "Three Norths Shelter Forest System Project (Phase V)".The construction of these projects is expected to control land desertification and slow down wind erosion by planting trees and grasses in areas under severe desertification in the northwest and north HLBB.The trend of desertification has started to reverse in recent years [32].However, some remote and broad sandy areas that are relatively difficult to treat are not controlled yet.Identifying the driving factors of wind erosion is the key to overcoming the obstacles of wind erosion control in HLBB.

Data Acquisition and Processing
Meteorological data were obtained from the meteorological stations in HLBB and its surrounding areas (14 stations in total) from the China Meteorological Data Service Center (http://data.cma.cn/).The dataset includes temperature, precipitation, wind speed, wind direction, solar radiation, and sunlight hours data from 2005 to 2016.Given the flat terrain of HLBB, the spatial grid of meteorological data was directly interpolated by using the Kriging method in ArcGIS 10.2.The snow depth data were provided by the Environmental and Ecological Science Data Center for West China, National Natural Science Foundation of China (http://westdc.westgis.ac.cn).
Normalized Difference Vegetation Index (NDVI) data were derived from the MOD13Q1: MODIS/Terra Vegetation Indices 16-Day L3 Global SIN Grid V006 (https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/mod13q1_v006).The Digital Elevation Model (DEM) and soil structure dataset, such as sand, silt and clay, CaCO 3 , and organic carbon content data, were obtained from the Geospatial Data Cloud site, Computer Network Information Center, and Chinese Academy of Sciences (http://www.gscloud.cn).All grid cells were resampled into 100 m × 100 m resolution.

The RWEQ Model
The RWEQ model is an empirical-based model designed by Fryrear et al. [26] to estimate soil eroded and transported by wind between the soil surface and a height of 2 m at a field scale.The RWEQ model has been widely used for estimating soil loss at a regional scale [33][34][35][36].The equation combines the factors that can affect wind erosion processes, which include wind speed, snow depth, soil, and vegetation [37].The rate of soil loss at a specific location can be computed by the model as follows: where SL in kg/m 2 is the rate of soil loss caused by wind erosion, which is also referred to as the wind erosion modulus; x in m is the distance from the upwind edge of the field; Q max in kg/m is the maximum transport capacity; S in m is the critical field length that is defined as the distance at which 63% of the maximum transport capacity is reached; the Q max and S are calculated by the weather factor (WF) in kg/m; the dimensionless factors are soil erodibility factor (EF), soil crust factor (SCF), soil roughness factor (K ), and combined vegetation factor (C).
WF is one of the most significant factors affecting wind erosion, which is determined by wind speed, precipitation, temperature, radiation, and snow cover.The weather factor can be calculated as follows: where U 2 in m/s is 2 m wind speed; U t is the threshold wind speed assumed as a constant of 5.0 m/s [26]; N d is the number of days during the evaluation period; N is the observation frequency of wind speeds during a period; ρ is the air density in kg/m 3 ; g is the acceleration due to gravity in m/s 2 ; SW and SD are the dimensionless soil moisture and snow-cover factor, respectively.Wind speed is usually the most important driving factor affecting wind erosion [38].The actual measured wind speed is generally used in the RWEQ model [26].If the wind speed data are missing or fail to meet the needs of this model, a wind speed simulation add-in named Wind Generator is often used to simulate wind speeds [38].However, the Wind Generator tends to underestimate wind speed data, especially for the high wind section [39].In this study, the actual measured hourly wind speed from China Meteorological Center is used in the model calculation.
The threshold wind speed of sand detachment increases with the increase in soil moisture.In June 2016, the research team obtained 12 soil profiles in Wudan Tara Village, HLBB.The soil moisture content was measured by a TDR100 time domain reflectometer in these 12 soil profiles.The results showed that the average soil water content was between 0.8% and 1%.Accordingly, the threshold wind speed should be approximately 4.5-5.0m/s at this water content [40].Therefore, the U t is determined to be 5.0 m/s, which means that wind erosion is not estimated when the wind speed is below 5.0 m/s.
RWEQ computes the soil EF and SCF based on soil physical and chemical properties [41], which can be calculated as follows: where Sa, Si, and Cl are the proportion of sand (5.50-93.60%),silt (0.50-69.50%), and clay (5.00-39.30%),respectively; Sa/Cl is the ratio of sand to clay (1.20-53.00%);OM is the organic matter content (0.32-4.74%);CaCO 3 is calcium carbonate content (0-25.20%).The brackets show the ranges of data required by the RWEQ model.The K factor is related to chain random roughness (Crr) and soil ridge roughness (K r ) and can be calculated as follows: where Crr is considered 0 in this study, because the effect of random roughness on wind erosion is generally negligible [42], and the spatial data of random roughness in HLBB is difficult to obtain.The K r is generally calculated by the Smith-Carson Equation as follows: where L and ∆H are the relief parameter and the elevation differences within the distance L in cm, respectively.Vegetation cover can also effectively reduce wind erosion [38].In this study, the crop seedlings that provide partial canopy cover to protect the soil (SLR c ) are calculated as the effects of vegetation cover on wind erosion: where CC in % is the vegetation coverage of each pixel; NDVI is the normalized difference vegetation index, which is one of the main characteristic descriptors indicating land surface vegetation cover; NDVI soil is the value of bare soil pixels; NDVI veg is the value of vegetated pixels.Their selection mainly depends on 95% confidence of NDVI pixel statistic values.

Slope Analysis
The non-parametric Mann-Kendall test proposed by Mann [43] and improved by Kendall [44] is widely used when the probability distribution of data is undefined [45].Since the distribution of wind erosion in a time series also does not follow a particular probability distribution, the Mann-Kendall test is used for trend and catastrophe analyses of wind erosion in this study.However, the Mann-Kendall test focuses on the temporal analysis rather than spatial distribution.Both wind erosion and its changes exhibit significant spatial heterogeneity.Thus, the linear regression analysis is used in each grid to detect the spatial characteristics of wind erosion change trends.The regression coefficient represented the trend slope, which is calculated as follows: where θ is the trend slope; n is the total number of years; SL i is wind erosion modulus in the ith year with i set to 1-12; SL is the average wind erosion modulus in each grid during this period.A negative θ indicates the grid exhibiting a decrease in wind erosion over the period, whereas a positive θ indicates an increase trend.Since the absolute amount of the average wind erosion affects the trend, a relative slope (γ) is calculated in each pixel to eliminate the average wind erosion effects.

Sensitivity Analysis and Contribution Assessment
Partial differentials are used to calculate the sensitivity of each driving factor during the wind erosion process [46]: where S v m is the sensitivity coefficient of SL that related to the mth climatic variable v m ; S v m indicates the sensitivity of SL under the condition when all variable changes are the same.However, the change rate in each grid varies for each wind erosion variable, such as wind speed, precipitation, temperature, and vegetation coverage.To determine the relative contribution of each variable to wind erosion, we used a contribution assessment method proposed by Yin et al. [47], which studies the sensitivity of SL and the change in trends of these variables.The equations are as follows: where CB v m is the contribution rate of meteorological variable v m to SL; RC v m is the multi-year relative change of v m ; v m is the multi-year absolute average of v m ; β v m is the magnitude of the trend of meteorological variable v m , which is expressed by Sen's slope estimator.n 0 is the number of periods we calculated in the wind erosion modulus (23 periods per year for 12 years), which equals 276.

Constraint Line Analysis
In bivariate scattergrams, data points sometimes show clouds bounded by an informative edge, which implies that the independent variable may act as a limiting factor constraining the response of the dependent variable (i.e., constraint effect).In this case, constraint line analysis is suggested instead of the correlation and regression methods [48][49][50].The constraint line helps researchers extract effective information from scattered data and understand the limiting effect among variables.Therefore, we quantified the effects of vegetation cover on wind erosion using the constraint line method as suggested by Zhao et al. [51].

Temporal-Spatial Changes of Wind Erosion from 2005 to 2016
The RWEQ model simulation results showed that the average wind erosion rate in HLBB from 2005 to 2016 was 2.88 t/hm 2 •a.The spatial distribution of wind erosion displayed a decreasing pattern from west to east (Figure 3).The highest and lowest annual wind erosion values occurred in 2006 and 2014 and were 8.The slope analysis results also showed that the wind erosion exhibited a clearly downward trend in the whole region.The declining slope (θ) of the absolute amount of wind erosion modulus was about 0.35 t/hm 2 •a in HLBB from 2005 to 2016 (Figure 5).From the spatial perspective, the declining slope generally decreased from the west to the east.Its spatial pattern was similar to that of the wind erosion modulus.The relative slope (γ) of wind erosion was −22.24%/a, which excluded the effects of absolute wind erosion modulus on the slope calculation.The value of γ was low in areas with high vegetation coverage.Its spatial distribution exhibited a high correlation with the spatial distribution of vegetation coverage.Figure 6 shows the proportion changes of different wind erosion intensity areas.From 2005 to 2016, the area of tolerable erosion intensity increased, and the area of slight erosion intensity significantly decreased.Moreover, the proportion of other types of wind erosion was less than 1% of the total area except for 2005 and 2006.The wind erosion modulus was evidently reduced from 2005 to 2016.The slope analysis results also showed that the wind erosion exhibited a clearly downward trend in the whole region.The declining slope (θ) of the absolute amount of wind erosion modulus was about 0.35 t/hm 2 •a in HLBB from 2005 to 2016 (Figure 5).From the spatial perspective, the declining slope generally decreased from the west to the east.Its spatial pattern was similar to that of the wind erosion modulus.The relative slope (γ) of wind erosion was −22.24%/a, which excluded the effects of absolute wind erosion modulus on the slope calculation.The value of γ was low in areas with high vegetation coverage.Its spatial distribution exhibited a high correlation with the spatial distribution of vegetation coverage.Figure 6 shows the proportion changes of different wind erosion intensity areas.From 2005 to 2016, the area of tolerable erosion intensity increased, and the area of slight erosion intensity significantly decreased.Moreover, the proportion of other types of wind erosion was less than 1% of the total area except for 2005 and 2006.The wind erosion modulus was evidently reduced from 2005 to 2016.The slope analysis results also showed that the wind erosion exhibited a clearly downward trend in the whole region.The declining slope (θ) of the absolute amount of wind erosion modulus was about 0.35 t/hm 2 •a in HLBB from 2005 to 2016 (Figure 5).From the spatial perspective, the declining slope generally decreased from the west to the east.Its spatial pattern was similar to that of the wind erosion modulus.The relative slope (γ) of wind erosion was −22.24%/a, which excluded the effects of absolute wind erosion modulus on the slope calculation.The value of γ was low in areas with high vegetation coverage.Its spatial distribution exhibited a high correlation with the spatial distribution of vegetation coverage.Figure 6 shows the proportion changes of different wind erosion intensity areas.From 2005 to 2016, the area of tolerable erosion intensity increased, and the area of slight erosion intensity significantly decreased.Moreover, the proportion of other types of wind erosion was less than 1% of the total area except for 2005 and 2006.The wind erosion modulus was evidently reduced from 2005 to 2016.

Contribution of Driving Factors to the Changes of Wind Erosion
In the RWEQ model, K′, SCF, and EF are constants, whereas WF and C are variables.To reveal the correlation between the spatial distribution of wind erosion, vegetation, and meteorological factors, we represented C by the index of vegetation coverage factor (VCF).In addition, we also divided WF into three indexes, namely wind speed factors (WSF), soil moisture factors (SMF), and air density factors (ADF).
The results showed that the decrease in WSF and the increase in VCF played key roles in wind erosion reduction.The contributions from the other two factors were relatively small.Firstly, the average contribution of WSF in HLBB was 47.40%, which decreased from west to east.This spatial pattern was similar to the trend of wind erosion.Secondly, the average contribution of VCF was 45.93%.The spatial distribution exhibited a strong correlation with vegetation coverage.Thirdly, the average contribution of SMF was about 5.43%.The high contribution was mainly distributed in the eastern part of HLBB.Finally, the average contribution of ADF was 0.96%, which increased from west to east (Figure 7).
To further identify the dominant factors of wind erosion changes in different regions, we compared the contributions of each factor at a grid scale.The results illustrated that wind erosion is mainly dominated by WSF and VCF collectively.However, the dominant factors vary in different regions.We then subtracted the contributions from WSF and VCF at a grid scale and determined the dominant factors by the difference.In the southwestern region, which accounts for 49.87% of HLBB, the WSF was the dominant factor contributing to the decline of wind erosion.The decline of wind erosion in the northeast region, which accounts for 50.13% of the study area, was mainly caused by the VCF (Figure 8).

Contribution of Driving Factors to the Changes of Wind Erosion
In the RWEQ model, K , SCF, and EF are constants, whereas WF and C are variables.To reveal the correlation between the spatial distribution of wind erosion, vegetation, and meteorological factors, we represented C by the index of vegetation coverage factor (VCF).In addition, we also divided WF into three indexes, namely wind speed factors (WSF), soil moisture factors (SMF), and air density factors (ADF).
The results showed that the decrease in WSF and the increase in VCF played key roles in wind erosion reduction.The contributions from the other two factors were relatively small.Firstly, the average contribution of WSF in HLBB was 47.40%, which decreased from west to east.This spatial pattern was similar to the trend of wind erosion.Secondly, the average contribution of VCF was 45.93%.The spatial distribution exhibited a strong correlation with vegetation coverage.Thirdly, the average contribution of SMF was about 5.43%.The high contribution was mainly distributed in the eastern part of HLBB.Finally, the average contribution of ADF was 0.96%, which increased from west to east (Figure 7).
To further identify the dominant factors of wind erosion changes in different regions, we compared the contributions of each factor at a grid scale.The results illustrated that wind erosion is mainly dominated by WSF and VCF collectively.However, the dominant factors vary in different regions.We then subtracted the contributions from WSF and VCF at a grid scale and determined the dominant factors by the difference.In the southwestern region, which accounts for 49.87% of HLBB, the WSF was the dominant factor contributing to the decline of wind erosion.The decline of wind erosion in the northeast region, which accounts for 50.13% of the study area, was mainly caused by the VCF (Figure 8).

Correlation Between the Wind-Related Characteristics and Wind Erosion
Wind speed is widely regarded as the driving force of wind erosion [21].Our results of the contribution analysis also confirmed that the decline of WSF was one of the main contributing factors to the declining trend of wind erosion in HLBB.The change in WSF exhibited critical impacts on the wind erosion modulus.Thus, it is necessary to explore relationship between wind-related characteristics and wind erosion.We selected three indicators to analyze, namely the average wind speed (AWS), the number of days when the wind speed exceeded the threshold of the sand detachment wind speed (NWSE), and the average wind speed exceeding the threshold of sand detachment wind speed (AWSE).At the annual scale, all the indicators exhibited clear downward trends from 2005 to 2016 (Figure 9).The annual AWS decreased 9.89% from 2006 to 2015, and the annual NWSE decreased 45.20%.At the same time, the annual AWSE dropped by 46.35%, which was the largest decline among the three.
However, the annual indicators could not reflect the subtle temporal relationship between wind speed and wind erosion since most wind erosion occurred in spring in HLBB.Thus, the correlations between the wind erosion modulus and the semi-monthly AWS, AWSE, and NWSE were calculated through the Spearman's rank correlation coefficients (SRCC).The results showed that the SRCC between the wind erosion modulus and the semi-monthly AWSE of 0.926 was greater than the SRCC between the wind erosion modulus and semi-monthly AWS of 0.821.Furthermore, the SRCC between the semi-monthly NWSE and wind erosion modulus reached 0.920, which was similar to the SRCC of semi-monthly AWSE and higher than the SRCC of semi-monthly AWS.The Spearman rank correlation coefficients of wind erosion modulus and different wind-related indicators suggested that people should pay more attention to the intensity and frequency of windy weather rather than annual AWS in wind erosion simulation studies.Using AWSE and NWSE as the wind indicators might be able to explain better the effects of wind on the wind erosion modulus.

Correlation Between the Wind-Related Characteristics and Wind Erosion
Wind speed is widely regarded as the driving force of wind erosion [21].Our results of the contribution analysis also confirmed that the decline of WSF was one of the main contributing factors to the declining trend of wind erosion in HLBB.The change in WSF exhibited critical impacts on the wind erosion modulus.Thus, it is necessary to explore the relationship between wind-related characteristics and wind erosion.We selected three indicators to analyze, namely the average wind speed (AWS), the number of days when the wind speed exceeded the threshold of the sand detachment wind speed (NWSE), and the average wind speed exceeding the threshold of sand detachment wind speed (AWSE).At the annual scale, all the indicators exhibited clear downward trends from 2005 to 2016 (Figure 9).The annual AWS decreased 9.89% from 2006 to 2015, and the annual NWSE decreased 45.20%.At the same time, the annual AWSE dropped by 46.35%, which was the largest decline among the three.
However, the annual indicators could not reflect the subtle temporal relationship between wind speed and wind erosion since most wind erosion occurred in spring in HLBB.Thus, the correlations between the wind erosion modulus and the semi-monthly AWS, AWSE, and NWSE were calculated through the Spearman's rank correlation coefficients (SRCC).The results showed that the SRCC between the wind erosion modulus and the semi-monthly AWSE of 0.926 was greater than the SRCC between the wind erosion modulus and semi-monthly AWS of 0.821.Furthermore, the SRCC between the semi-monthly NWSE and wind erosion modulus reached 0.920, which was similar to the SRCC of semi-monthly AWSE and higher than the SRCC of semi-monthly AWS.The Spearman rank correlation coefficients of wind erosion modulus and different wind-related indicators suggested that people should pay more attention to the intensity and frequency of windy weather rather than annual AWS in wind erosion simulation studies.Using AWSE and NWSE as the wind indicators might be able to explain better the effects of wind on the wind erosion modulus.

Relationship Between Vegetation Coverage and Wind Erosion
Vegetation coverage is widely considered to be effective in mitigating wind erosion [52].High vegetation coverage can protect the soil from wind erosion.According to the contribution assessment results, the increment of vegetation coverage is another main contributor to wind erosion reduction.In the past 12 years, approximately 78.06% of HLBB experienced increased vegetation coverage, which mainly occurred in the western region (Figure 10).However, the increase in vegetation cannot reduce wind erosion indefinitely.When vegetation coverage exceeds 60%, any further increase in vegetation coverage will no longer result in wind erosion reduction [51].Moreover, since vegetation has different morphologies at different wind speeds [53], the effects of vegetation on wind erosion mitigation might vary accordingly.Therefore, the relationship between vegetation coverage and wind erosion modulus was analyzed at different wind speeds.

Relationship Between Vegetation Coverage and Wind Erosion
Vegetation coverage is widely considered to be effective in mitigating wind erosion [52].High vegetation coverage can protect the soil from wind erosion.According to the contribution assessment results, the increment of vegetation coverage is another main contributor to wind erosion reduction.In the past 12 years, approximately 78.06% of HLBB experienced increased vegetation coverage, which mainly occurred in the western region (Figure 10).However, the increase in vegetation cannot reduce wind erosion indefinitely.When vegetation coverage exceeds 60%, any further increase in vegetation coverage will no longer result in wind erosion reduction [51].Moreover, since vegetation has different morphologies at different wind speeds [53], the effects of vegetation on wind erosion mitigation might vary accordingly.Therefore, the relationship between vegetation coverage and wind erosion modulus was analyzed at different wind speeds.The constraint line method suggested that the increase in vegetation coverage contributed to the slowdown of wind erosion in HLBB.A maximum slowing effect of vegetation on wind erosion was also present.When the vegetation coverage reached approximately 40% and the semi-monthly AWS was greater than 3 m/s, further increase in vegetation would not lead to obvious declines in wind erosion.However, the constraint line changed at different wind speeds.When the semi-monthly AWS was between 2 and 3 m/s, the threshold of vegetation coverage was significantly reduced to 20%.When the semi-monthly AWS was less than 2 m/s, there was no significant threshold of vegetation coverage since the wind erosion was already low (Figure 11).The constraint line moved from top left to bottom right with decreasing semi-monthly AWS, and the threshold of vegetation coverage also decreased.
In summary, the increase in vegetation slows down wind erosion within a certain threshold.Meanwhile, the same vegetation coverage has different effects on wind erosion at different wind speeds.Approximately 25% of vegetation coverage should be enough to resist erosion caused by lowwind speed (less than 2 m/s).However, once the wind speed increases, this coverage of vegetation is The constraint line method suggested that the increase in vegetation coverage contributed to the slowdown of wind erosion in HLBB.A maximum slowing effect of vegetation on wind erosion was also present.When the vegetation coverage reached approximately 40% and the semi-monthly AWS was greater than 3 m/s, further increase in vegetation would not lead to obvious declines in wind erosion.However, the constraint line changed at different wind speeds.When the semi-monthly AWS was between 2 and 3 m/s, the threshold of vegetation coverage was significantly reduced to 20%.When the semi-monthly AWS was less than 2 m/s, there was no significant threshold of vegetation coverage since the wind erosion was already low (Figure 11).The constraint line moved from top left to bottom right with decreasing semi-monthly AWS, and the threshold of vegetation coverage also decreased.
In summary, the increase in vegetation slows down wind erosion within a certain threshold.Meanwhile, the same vegetation coverage has different effects on wind erosion at different wind speeds.Approximately 25% of vegetation coverage should be enough to resist erosion caused by low-wind speed (less than 2 m/s).However, once the wind speed increases, this coverage of vegetation is insufficient.Vegetation coverage should reach approximately 40% to achieve similar protective effectiveness at higher wind speeds.
insufficient.Vegetation coverage should reach approximately 40% to achieve similar protective effectiveness at higher wind speeds.

Influence of Natural Driving Factors on the Wind Erosion Modulus
Previous research has generally studied wind erosion and its correlations between different driving factors based on annual statistics [14,24].However, the results of this study suggest that wind erosion has a concentrated distribution during the year, which cannot be explained by interannual statistical results.The occurrence of wind erosion should be regarded as an event rather than a longterm state.The wind erosion modulus measures long-term regional wind erosion conditions.However, the wind erosion modulus of HLBB exhibited temporal heterogeneity within a year.The period of high wind erosion mainly concentrated in spring, which coincided with the period of the windy season.The proportion of the wind erosion modulus in spring accounted for more than 81% of the total annual erosion every year from 2005 to 2016 except for 2014.The highest value reached 98.33% in 2012.Therefore, it is necessary to explore the relationship between wind erosion events and its driving factors on a higher time-resolution to reveal the concealed relationships in the interannual analysis.In the following sections, the temporal distributions of wind erosion and driving factors within a year were analyzed on a semi-monthly scale.The temporal distributions were the same as those in the period during which we calculated the wind erosion modulus.

Influence of Natural Driving Factors on the Wind Erosion Modulus
Previous research has generally studied wind erosion and its correlations between different driving factors based on annual statistics [14,24].However, the results of this study suggest that wind erosion has a concentrated distribution during the year, which cannot be explained by interannual statistical results.The occurrence of wind erosion should be regarded as an event rather than a long-term state.The wind erosion modulus measures long-term regional wind erosion conditions.However, the wind erosion modulus of HLBB exhibited temporal heterogeneity within a year.The period of high wind erosion mainly concentrated in spring, which coincided with the period of the windy season.The proportion of the wind erosion modulus in spring accounted for more than 81% of the total annual erosion every year from 2005 to 2016 except for 2014.The highest value reached 98.33% in 2012.Therefore, it is necessary to explore the relationship between wind erosion events and its driving factors on a higher time-resolution to reveal the concealed relationships in the interannual analysis.In the following sections, the temporal distributions of wind erosion and driving factors within a year were analyzed on a semi-monthly scale.The temporal distributions were the same as those in the period during which we calculated the wind erosion modulus.

Mismatch Between Wind Speed and Vegetation Coverage
Wind speed and vegetation coverage are the most important factors affecting wind erosion.The reduction of wind speed and increment of vegetation coverage exhibited positive effects on annual wind erosion reduction [22].However, there is a temporal mismatch between wind speed and vegetation coverage within a year.The windy period in HLBB is from mid-February to mid-June.Wind speed just starts to increase from spring.Meanwhile, vegetation growth has not yet occurred (Figure 12).Wind speed in the windy season is high enough to detach the sand particles, but low vegetation coverage can only provide minimal protection for bare sand.In the windy periods, the semi-monthly NWSE accounts for 66.81% of the whole year.The semi-monthly AWSE is higher than 95.16% of the annual AWSE.In spring, when the low vegetation coverage in HLBB is not able to resist erosion, the frequently occurring high wind speed results in a significant increase in wind erosion.Under the dual influences from wind and vegetation, the wind erosion during this period is very severe.In mid-June, local plants begin their growth.However, the wind speed also begins to decline at this time, which will reduce the wind erosion.Therefore, even after the vegetation cover is fully established, the plants will not fully exert their slowing effects on wind erosion.The temporal mismatch between the occurrences of high wind speed and high vegetation coverage challenges regional wind erosion control.
Wind speed and vegetation coverage are the most important factors affecting wind erosion.The reduction of wind speed and increment of vegetation coverage exhibited positive effects on annual wind erosion reduction [22].However, there is a temporal mismatch between wind speed and vegetation coverage within a year.The windy period in HLBB is from mid-February to mid-June.Wind speed just starts to increase from spring.Meanwhile, vegetation growth has not yet occurred (Figure 12).Wind speed in the windy season is high enough to detach the sand particles, but low vegetation coverage can only provide minimal protection for bare sand.In the windy periods, the semi-monthly NWSE accounts for 66.81% of the whole year.The semi-monthly AWSE is higher than 95.16% of the annual AWSE.In spring, when the low vegetation coverage in HLBB is not able to resist erosion, the frequently occurring high wind speed results in a significant increase in wind erosion.Under the dual influences from wind and vegetation, the wind erosion during this period is very severe.In mid-June, local plants begin their growth.However, the wind speed also begins to decline at this time, which will reduce the wind erosion.Therefore, even after the vegetation cover is fully established, the plants will not fully exert their slowing effects on wind erosion.The temporal mismatch between the occurrences of high wind speed and high vegetation coverage challenges regional wind erosion control.Precipitation and temperature mainly affect wind erosion by influencing soil moisture and vegetation growth.Precipitation increases soil moisture, which makes sand detachment more difficult.In contrast, an increase in temperature increases evapotranspiration and reduces soil moisture.In HLBB, precipitation is mainly concentrated from the end of April to the end of August.During the rainy season, soil moisture is at a high level even when the air temperature is also elevated.After mid-September, the precipitation drops sharply.Consequently, the soil moisture is kept at a low level until the next rainy season.The soil then becomes dry and breaks easily by wind during the winter and spring, which exacerbates the wind erosion situation in spring.
More importantly, vegetation growth is often determined by the rise in temperature and precipitation with a certain hysteresis effect at the local scale [54].By the end of May, the vegetation coverage increases after the increase in precipitation and temperature at the end of April (Figure 12).As a result, vegetation coverage is not able to effectively slow down wind erosion before June.In contrast, an increase in temperature increases evapotranspiration and reduces soil moisture.In HLBB, precipitation is mainly concentrated from the end of April to the end of August.During the rainy season, soil moisture is at a high level even when the air temperature is also elevated.After mid-September, the precipitation drops sharply.Consequently, the soil moisture is kept at a low level until the next rainy season.The soil then becomes dry and breaks easily by wind during the winter and spring, which exacerbates the wind erosion situation in spring.

Targeted Wind Erosion Prevention Strategy
More importantly, vegetation growth is often determined by the rise in temperature and precipitation with a certain hysteresis effect at the local scale [54].By the end of May, the vegetation coverage increases after the increase in precipitation and temperature at the end of April (Figure 12).As a result, vegetation coverage is not able to effectively slow down wind erosion before June.

Targeted Wind Erosion Prevention Strategy
The results show that the increase in vegetation coverage is one of the most important driving forces of the wind erosion decline.Although it is difficult to greatly increase vegetation coverage in spring due to the low temperature and precipitation in this season, the vegetation coverage can survive in the spring of HLBB.The increase in vegetation coverage was mainly caused by large-scale artificial afforestation since there was no significant change in local temperature and precipitation during the period from 2005 to 2016.In order to alleviate land degradation and soil losses caused by wind erosion, the local government launched a series of afforestation projects, such as the "5820 Forestry Ecological Project" in 1999, which built concentrated and contiguous forest within five years, and eventually afforested 133,333 hectares of forest and 66,666 hectares of grassland.Recently, the "Double Ten Million Mu Project" (1 mu equals 0.0667 hectares), which plans to effectively manage 10 million mu of forest and 10 million mu of grassland in seven years, has been implemented.The purposes of this project are to reduce soil wind erosion and land degradation as well as improve the quality of the local environment by gradually restoring the sparse forest ecosystems.These afforestation projects play a positive role in the decline of wind erosion in HLBB.However, the wind erosion problem has not yet been fully managed.For example, vegetation restoration in areas under severe desertification is still difficult, which makes spring sandstorms still frequent.
From the above discussion, we can find the key factor in wind erosion control in HLBB, namely, the domination of the high wind erosion period in spring, which accounts for more than 80% of annual wind erosion.The results of this study provide two recommendations for wind erosion control.Firstly, vegetation coverage must reach the minimum thresholds as soon as possible to effectively reduce wind erosion.Secondly, at high wind speeds, vegetation coverage must be further increased to reach the corresponding threshold.Therefore, the ideal situation is the presence of enough vegetation coverage in the windy season.However, in reality, the mismatch between wind speed, precipitation, and temperature impeded the difficulties of wind erosion control in HLBB.If people can solve this mismatch, wind erosion in spring should be greatly alleviated.Given the insufficient precipitation and low temperature, increasing vegetation coverage via planting more trees and grasses might not be the most cost-effective method.Some management techniques, such as maintaining the residues on farmlands in winter, creating windbreaks, and having drought-tolerant bush covers, can be considered as erosion control methods during the high wind erosion season.These measures can alleviative wind erosion in spring through improving the surface roughness, reducing the wind speed on the surface, and capturing airborne dust, which hinders the process of sand transmission [55].However, HLBB is located in a semi-arid area with limited water, which makes the survival rate of existing windbreak plants low [56][57][58][59][60].Therefore, strategies such as the cultivation of optimal windbreak and shrub species, and careful selection of planting locations are needed to improve the mitigating effects of windbreaks.

Limitations
The RWEQ model is an empirical model based on various empirical parameters.These parameters are obtained from the statistical analyses of multiple data [37].Given that the RWEQ model was originally designed at a field scale, some problems might be raised during regional-scale wind erosion simulation, such as the adjustment of empirical parameters and the determination of the wind speed threshold of sand detachment.A uniform parameter setting might also cause simulation result deviations due to the spatial heterogeneity of the natural conditions.Therefore, applying the RWEQ model to simulate the wind erosion modulus in other regions might require adjustments of the model [61].
The RWEQ model has been widely used in arid regions of northern China, especially in Inner Mongolia.Some studies used the measured data to verify the model simulation accuracy.The results showed that the RWEQ model exhibited good adaptability in China [22,62].In the setting of wind speed parameters, we used the method of TDR (Time Domain Reflectometry) to determine the soil water content, which provided the basis to determining the local wind speed threshold.We also used the hourly wind speed with higher time resolution as the data source for the wind factor calculation, which is more consistent with the original RWEQ model.Compared with the results obtained from other studies [14,19,22], the wind erosion simulation results of this study are consistent with the modulus and trends of wind erosion reported there.However, there are also some limitations that need to be improved in future studies, for example, revised K and C parameters based on different land cover types.The Cs-137 tracing technique can be used to verify the wind erosion simulation results in the study area.

Conclusions
Based on the RWEQ model, we simulated the wind erosion modulus and mapped the temporal and spatial variations of wind erosion in HLBB from 2005 to 2016.Furthermore, the contribution analysis method was applied to analyze the dominant driving factors of wind erosion.The main conclusions are as follows: (1) The average modulus of wind erosion was 2.88 t/hm 2 •a in HLBB from 2005 to 2016, which generally increased from west to east.The wind erosion modulus significantly dropped by 0.35 t/hm 2 •a from 2005 to 2016.The spatial pattern of the relative decline rate was similar to the distribution of vegetation.The main causes for the wind erosion reduction in HLBB were the decrease in wind speed and the increase in vegetation coverage.The decline in wind erosion in the southeastern region, which accounts for 49.87% of the study area, was mainly attributed to wind speed reduction.The decline in wind erosion in the northeast region, which accounts for 50.13% of the study area, was mainly due to the increase in vegetation coverage.(2) Wind speed exhibited a significant correlation with the wind erosion modulus.However, the correlation between the annual AWSE and NWSE, and the wind erosion modulus was higher than that between the annual AWS and the wind erosion modulus.Vegetation coverage had different thresholds for effective wind erosion mitigation at different wind speeds.The threshold of vegetation coverage was approximately 40% when the semi-monthly AWS was greater than 3 m/s.The threshold became approximately 20% when the semi-monthly AWS was between 2 m/s and 3 m/s.(3) Wind erosion was mainly concentrated in mid-February to mid-June in HLBB because of the high wind speed, low precipitation, and low vegetation coverage.The temporal mismatch between wind speed and vegetation coverage was the main reason for the severe wind erosion in spring.
Because of this mismatch, more targeted strategies need to be applied to mitigate the high spring wind erosion.Given that improving the vegetation coverage in spring in HLBB is very difficult, reducing the wind speed in spring by adding drought-tolerant windbreaks might be a more cost-effective method for wind erosion reduction.

20 Figure 3 .
Figure 3.The spatial distribution of wind erosion per year in HLBB from 2005 to 2016.

Figure 3 .
Figure 3.The spatial distribution of wind erosion per year in HLBB from 2005 to 2016.

Figure 4 .
Figure 4.The annual wind erosion from 2005 to 2016 in HLBB.

Figure 5 .
Figure 5.The spatial distribution of wind erosion change trends in HLBB from 2005 to 2016.(a) the slope of wind erosion (θ), (b) the relative slope of wind erosion (γ).

Figure 4 .
Figure 4.The annual wind erosion from 2005 to 2016 in HLBB.

Figure 5 .
Figure 5.The spatial distribution of wind erosion change trends in HLBB from 2005 to 2016.(a) the slope of wind erosion (θ), (b) the relative slope of wind erosion (γ).

Figure 5 .
Figure 5.The spatial distribution of wind erosion change trends in HLBB from 2005 to 2016.(a) the slope of wind erosion (θ), (b) the relative slope of wind erosion (γ).

Figure 7 .
Figure 7.The spatial distribution of the contribution of different factors to average annual wind erosion changes from 2005 to 2016; (a) contribution of wind speed factor, (b) contribution of vegetation coverage factor, (c) contribution of soil moisture factor, and (d) contribution of air density factor.

Figure 8 .Figure 7 . 20 Figure 7 .
Figure 8.The distribution of driving factors to the average annual change of wind erosion in HLBB from 2005 to 2016.

Figure 8 .Figure 8 .
Figure 8.The distribution of driving factors to the average annual change of wind erosion in HLBB from 2005 to 2016.

Figure 9 .
Figure 9.The change in wind-related indicators from 2005 to 2016 in HLBB.Figure 9.The change in wind-related indicators from 2005 to 2016 in HLBB.

Figure 9 .
Figure 9.The change in wind-related indicators from 2005 to 2016 in HLBB.Figure 9.The change in wind-related indicators from 2005 to 2016 in HLBB.

Figure 10 .
Figure 10.The average trend of annual changes of vegetation coverage in HLBB from 2005 to 2016.

Figure 10 .
Figure 10.The average trend of annual changes of vegetation coverage in HLBB from 2005 to 2016.

Figure 11 .
Figure 11.The relationships between the wind erosion modulus and vegetation coverage under different semi-monthly average wind speed (AWS).Note: Constraint lines 1 and 2 are the relationship between the wind erosion modulus and the semi-monthly AWS above 3 m/s and between 2 and 3 m/s, respectively.

Figure 11 .
Figure 11.The relationships between the wind erosion modulus and vegetation coverage under different semi-monthly average wind speed (AWS).Note: Constraint lines 1 and 2 are the relationship between the wind erosion modulus and the semi-monthly AWS above 3 m/s and between 2 and 3 m/s, respectively.

Figure 12 .
Figure 12.The semi-monthly temporal characteristics of wind erosion and four natural factors.Note: all factors are normalized by the maximum and minimum values.

Figure 12 .
Figure 12.The semi-monthly temporal characteristics of wind erosion and four natural factors.Note: all factors are normalized by the maximum and minimum values.4.1.2.Relationships Between Precipitation, Temperature, and Wind Erosion Precipitation and temperature mainly affect wind erosion by influencing soil moisture and vegetation growth.Precipitation increases soil moisture, which makes sand detachment more difficult.In contrast, an increase in temperature increases evapotranspiration and reduces soil moisture.In HLBB, precipitation is mainly concentrated from the end of April to the end of August.During the rainy season, soil moisture is at a high level even when the air temperature is also elevated.After mid-September, the precipitation drops sharply.Consequently, the soil moisture is kept at a low level until the next rainy season.The soil then becomes dry and breaks easily by wind during the winter and spring, which exacerbates the wind erosion situation in spring.More importantly, vegetation growth is often determined by the rise in temperature and precipitation with a certain hysteresis effect at the local scale[54].By the end of May, the vegetation coverage increases after the increase in precipitation and temperature at the end of April (Figure12).As a result, vegetation coverage is not able to effectively slow down wind erosion before June.