Soil Water Behavior of Sandy Soils under Semiarid Conditions in the Shendong Mining Area (China)

: The soil water behavior of sandy soils was studied under semiarid conditions in the Shendong mining area (China). The soil water content ( θ ) was measured under different depths and topographies using an HH2 moisture meter. The inﬁltration process was studied using a Guelph soil permeameter. A set of hydrodynamic variables was calculated in the laboratory. The θ of the ﬁrst 20 cm was the lowest and increased with depth. The content of soil water increased from the top slope to the bottom slope. The inﬁltration experiments showed that the steady state inﬁltration rate was >40 mm h − 1 in most cases. Owing to the higher contents of sand and soil macropores at the top of the slope and the top 0–20 cm of surface soil, the initial inﬁltration rate and steady inﬁltration rate were higher. The average available water capacity was 18.28%, which was consistent with the predominance of a sandy textural fraction. The results of a soil water retention curve and a rainfall simulation experiment showed that there was a low soil water retention capacity throughout the whole proﬁle. This study contributes to the understanding of several aspects of the soil water behavior of sandy soils and provides key information for environmental management and land reclamation under semiarid conditions in the Shendong mining area.


Introduction
Soil water resources, an important part of the water cycle, are the basis for the survival of plants. In arid and semiarid regions, owing to the extremely fragile ecological system, soil moisture serves as the key ecological and environmental factor of plant growth [1][2][3][4][5]. It is an important indicator for monitoring land degradation [6,7] and an important input factor of terrestrial ecosystems [8]. It is also the key element of agricultural systems and hydrological and climatic elements [9].
In arid and semiarid environments, water stress affects many properties of the soil, including the decay of the soil's organic matter (OM), the characteristics of soil infiltration, the composition of soil granulometry, and the mineralization of soil nutrients [10,11]. OM provides necessary nutrients and improves soil fertility [12]. In addition, the OM retains water, stabilizes the soil structure, and renders the soil more resistant to degradation [13]. The infiltration process of soil water determines the redistribution of rainwater on the surface, which has an important influence on the content of water in the topsoil [14][15][16]. Soil moisture is affected by the soil's granulometric composition. Particularly in semiarid areas, where the evaporation is greater than the precipitation in the summer, the soil granulometric composition has a significant impact on the soil water content [17,18]. Consequently,

Soil Moisture Measurements
The surface soil moisture (θ) was measured at the six most representative points of the sandy soils of the study area. Two points were installed at the top (TS1 and TS2), middle (MS3 and MS4), and bottom (BS5 and BS6) of the slope. The points where the measurements were conducted in the study area included the top of slope, which was primarily covered by S. psammophila and Caragana, and the middle and bottom of the slope, with a predominant cover of sea buckthorn and Caragana. In all cases, the surface was undulating, and the slope gradients did not exceed 5%. An HH2 Moisture Meter (Theta-probe Type ML2x; Delta-T Devices, Ltd., Cambridge, UK) [31] was used, with three probes that were 6.8 cm long and three different analogue outputs. One of the probes was used to measure the volumetric water content; the others can be used to measure the bulk electrical conductivity and temperature, but they were not used for this paper. The probes were installed vertically at each soil moisture measurement point and measured the soil moisture at 0-20, 20-40, and 40-60 cm below the soil surface, respectively. Three duplicate readings were measured at each measurement point, which was sampled every hour [24]. The positions of the measurement sites were recorded by GPS.

Soil Moisture Measurements
The surface soil moisture (θ) was measured at the six most representative points of the sandy soils of the study area. Two points were installed at the top (TS1 and TS2), middle (MS3 and MS4), and bottom (BS5 and BS6) of the slope. The points where the measurements were conducted in the study area included the top of slope, which was primarily covered by S. psammophila and Caragana, and the middle and bottom of the slope, with a predominant cover of sea buckthorn and Caragana. In all cases, the surface was undulating, and the slope gradients did not exceed 5%. An HH2 Moisture Meter (Theta-probe Type ML2x; Delta-T Devices, Ltd., Cambridge, UK) [31] was used, with three probes that were 6.8 cm long and three different analogue outputs. One of the probes was used to measure the volumetric water content; the others can be used to measure the bulk electrical conductivity and temperature, but they were not used for this paper. The probes were installed vertically at each soil moisture measurement point and measured the soil moisture at 0-20, 20-40, and 40-60 cm below the soil surface, respectively. Three duplicate readings were measured at each measurement point, which was sampled every hour [24]. The positions of the measurement sites were recorded by GPS.
To study the soil water retention capacity, the after-rainfall soil moisture also needed to be measured at each control point. During the period of 4 June 2016 to 5 June 2016, there was 22.2 mm of rainfall in the study area based on the records of the meteorological station located 10 km away from the study site. When it had finished raining, a 1 m × 1 m × 1 m profile was excavated from the slope at TS1 and TS2, MS3 and MS4, and BS5 and BS6 with To study the soil water retention capacity, the after-rainfall soil moisture also needed to be measured at each control point. During the period of 4 June 2016 to 5 June 2016, there was 22.2 mm of rainfall in the study area based on the records of the meteorological station located 10 km away from the study site. When it had finished raining, a 1 m × 1 m × 1 m profile was excavated from the slope at TS1 and TS2, MS3 and MS4, and BS5 and BS6 with a shovel. An HH2 Moisture Meter (Theta-probe Type ML2x; Delta-T Devices, Ltd., Cambridge, UK) was used to measure three replicates of soil moisture along the soil profile at 0-20, 20-40, and 40-60 cm deep. The measurements were made at intervals of 0, 1, 10, and 48 h after the rain had fallen.

Soil Infiltration Experiments
The soil infiltration experiments were performed near six soil moisture control points to conduct a detailed study of the topsoil's water infiltration features and saturated hydraulic conductivity, using a Guelph soil permeameter (Zeal Quest Scientific Technology Co., Ltd., Shanghai, China) as described by Reynolds et al. [32], Li et al. [16], and Kuráž [33]. A detailed description of the Guelph soil permeameter has been provided by Kuráž [33] and Reynolds et al. [32]. At least three duplicate measurements of the topsoil were taken at each Water 2022, 14, 2159 4 of 15 measurement site. A constant depth of 20 cm was selected for the boreholes, as described by Kuráž [33] and Reynolds et al. [32]. A brush was used to clean the wall of borehole to eliminate the non-permeable membrane produced by the auger before measuring. The saturated hydraulic conductivity (Ksat) was calculated according to the following equation: where C is a dimensionless factor that is a function of the soil texture and the H/a ratio; Q (m 3 ·s −1 ) is the quasi-steady flow infiltrating into the vertical borehole; a (m) is the radius of the borehole; and H (m) is the height of the water level maintained above the borehole bottom [33].

Soil Sampling and Laboratory Analyses
To calculate the soil's physical and hydrodynamic properties in the laboratory, each plot was divided into two equal subplots to sample disturbed and undisturbed soil [34]. Three duplicate soil samples were taken from the center of each subplot at depths of 0-20 cm, 20-40, and 40-60 cm in June 2016.
The bulk density (BD) was determined using a 100 cm 3 (5.1 cm in diameter × 5 cm high) steel cylinder on the undisturbed soil samples, close to where the soil moisture had been measured [35][36][37][38]. Approximately 100 g of disturbed soil samples were collected for analyses of the soil's granulometric composition (GC) as previously described [24,[39][40][41] and the OM content using the Walkley-Black wet digestion method [24,42]. Based on a study by Zhang et al. [43], the following equation was used to fit the soil water retention curve: where θ (%) is the volumetric water content; ϕ (kPa) is the soil water suction, and a and b are fitting parameters. The soil's volumetric water content under different amounts of soil water suction (six different matric potentials ϕ: 100, 200, 500, 700, 1300 and 1500 kPa) was determined using an HH2 Moisture Meter (Delta-T, Inc., Makati, Philippines). The value of θ was measured at nine soil water potential values, which ranged from −33 to −1500 kPa [24]. The field capacity (FC) equaled the content of soil moisture at a soil water potential of −33 kPa, and the wilting point (WP) equaled the content of soil moisture at a soil water potential of −1500 kPa [10,24,44]. Thus, the available water capacity (AWC) of each collected soil sample was calculated using the following equation: AWC = FC − WP.

Statistical Analysis
The data were analyzed using Microsoft Excel (Redmond, WA, USA) and Origin 9.1 for Windows (OriginLab, Northampton, MA, USA). A Pearson's correlation coefficient analysis was used to determine the relationships between the soil moisture and the physical properties. The spatial distribution of the soil moisture was estimated using the geostatistical analytical software ArcGIS version 10.1 (ESRI, Redlands, CA, USA).

Soil Moisture and Physical Properties
The spatial distribution (0-20, 20-40, and 40-60 cm) of the soil moisture is shown in Figure 2. The amount of moisture in the first 20 cm of topsoil was between 8.88% and 16.30%, and that of the 20-40 cm and 40-60 cm topsoil was from 9.10% to 29.0% and 10.0% to 36.10%, respectively. The soil moisture of the first 20 cm was obviously lower, since it increased with depth. The mean water content of the 0-20, 20-40, and 40-60 cm topsoil was 13.73%, 16.97%, and 18.26%, respectively. The soil moisture increased with a decrease in the slope height. The highest soil moisture was observed at the bottom of the slope (BS). The values of soil moisture were in the following order: BS > MS > TS. Figure 2. The amount of moisture in the first 20 cm of topsoil was between 8.88% and 16.30%, and that of the 20-40 cm and 40-60 cm topsoil was from 9.10% to 29.0% and 10.0% to 36.10%, respectively. The soil moisture of the first 20 cm was obviously lower, since it increased with depth. The mean water content of the 0-20, 20-40, and 40-60 cm topsoil was 13.73%, 16.97%, and 18.26%, respectively. The soil moisture increased with a decrease in the slope height. The highest soil moisture was observed at the bottom of the slope (BS). The values of soil moisture were in the following order: BS > MS > TS. The soil texture classification standards of the USDA indicated that the soil texture at the soil sampling points was loamy sand and sandy loam ( Figure 3). The mean fractions of sand, silt, and clay were 80.50%, 9.77%, and 9.73%, respectively ( Table 1). The clay fraction was highly variable and ranged from 3.73% to 13.35%. The silt fraction was also quite variable and ranged between 7.66% and 13.89%. The sand fraction was less variable (coefficient of variation, CV = 3.92%) compared with the higher spatial variability of the silt (CV = 14.22%) and clay (CV = 29.32%) fractions. An important decrease in the sand fraction was observed at the sampling sites from the top of the slope to the bottom of the slope, and the fraction of clay was higher at the bottom of the slope, accordingly. Table 1 also shows the marked regularity of the soil texture along the profile (0-20, 20-40, and 40-60 cm). The sand fraction decreased with depth, while the clay fraction increased. These phenomena are primarily attributed to the migration of fine earth. This result was similar to that of Ceballos et al. [10].
The OM content was between 0.29% and 0.60% with a mean value of 0.45% and a CV value of 22.85%. At each sampling site, the amount of OM in the first layer (0-20 cm) was obviously higher than those of the other two layers (20-40 and 40-60 cm). Liu et al. [44] and Ceballos et al. [10] found similar results in sandy soils. At the top of a slope, the soil's OM content was lower than that at the bottom of the slope. This could be explained by the multiple influences of the decomposition of animal and plant residues, root exudates, microbial synthesis, and surface runoff. As shown in Table 1, the values of the soil BD were very similar, with a CV value of 4.24% and an average measured value of 1.69 g cm -3 . The The soil texture classification standards of the USDA indicated that the soil texture at the soil sampling points was loamy sand and sandy loam ( Figure 3). The mean fractions of sand, silt, and clay were 80.50%, 9.77%, and 9.73%, respectively ( Table 1). The clay fraction was highly variable and ranged from 3.73% to 13.35%. The silt fraction was also quite variable and ranged between 7.66% and 13.89%. The sand fraction was less variable (coefficient of variation, CV = 3.92%) compared with the higher spatial variability of the silt (CV = 14.22%) and clay (CV = 29.32%) fractions. An important decrease in the sand fraction was observed at the sampling sites from the top of the slope to the bottom of the slope, and the fraction of clay was higher at the bottom of the slope, accordingly. Table 1 also shows the marked regularity of the soil texture along the profile (0-20, 20-40, and 40-60 cm). The sand fraction decreased with depth, while the clay fraction increased. These phenomena are primarily attributed to the migration of fine earth. This result was similar to that of Ceballos et al. [10].     The OM content was between 0.29% and 0.60% with a mean value of 0.45% and a CV value of 22.85%. At each sampling site, the amount of OM in the first layer (0-20 cm) was obviously higher than those of the other two layers (20-40 and 40-60 cm). Liu et al. [44] and Ceballos et al. [10] found similar results in sandy soils. At the top of a slope, the soil's OM content was lower than that at the bottom of the slope. This could be explained by the multiple influences of the decomposition of animal and plant residues, root exudates, microbial synthesis, and surface runoff. As shown in Table 1, the values of the soil BD were very similar, with a CV value of 4.24% and an average measured value of 1.69 g cm -3 . The values of the saturated soil hydraulic conductivity (Ksat) were between 5.01 and 54.0 cm h -1 with a mean value of 20.38 cm h -1 and high spatial variability (CV = 79.91%). The Ksat of the top of the slope was higher than those in the rest of the study area, with a maximum value of 54.0 cm h -1 . The lowest Ksat was found at the middle of the slope, with a value of 5.01 cm h -1 .
The soil moisture content is the key environmental factor for semiarid areas. It is generally influenced by the soil's physical properties, such as the soil texture, BD, the OM content, and the hydraulic conductivity. In addition, the soil moisture content is also related to topographic factors. The influence of these factors cannot be studied independently of each other. To better understand these interactions in a semiarid area of north China, the relationships between the soil moisture and soil texture, BD, OM content, hydraulic conductivity, and topographic factors were analyzed in this study. Figure 4 shows the relationships between soil moisture and soil texture, BD, OM content, and hydraulic conductivity. Figure 4a-c show the relationship between the soil moisture and soil particles in the study area. These analyses revealed that the soil moisture content and the sand and silt content were negatively correlated, with correlation coefficients of −0.3426 and −0.1945, respectively. The content of clay in the soil was positively correlated with the soil moisture content, with a correlation coefficient of 0.4735. This phenomenon is supported by previous research [11,17]. Zhang et al. [11] found that the soil moisture was primarily affected by the content of clay, and the relationship between contents of soil moisture and clay was highly significant. Li [17] showed that the soil water content and the soil particle composition are directly related. Smaller soil particles result in a greater surface area of the soil and a higher maximum moisture content. This further explains the relationship between the soil moisture content and the soil particles. The relationship between soil moisture and BD is shown in Figure 4d. A negative linear correlation was used to fit the relationship between soil moisture and BD, but the low R value indicated that such a relationship might not exist. The value equaled t −0.0988. The BD showed very similar values and a mean value of 1.69 g cm -3 . The BD is a less important variable for the soil moisture when compared with the soil texture, the OM content, and the hydraulic conductivity in this study. Figure 4e shows the relationship between soil moisture and the OM content. The fitted linear curve for the soil moisture and organic matter was y = 0.3921 + 0.0036x, with R = 0.1505. As previously described [38,45], the soil OM was sensitive to soil moisture. Research by Zhang and Shao [38] showed a significant effect of the soil OM content on soil moisture. The soil moisture content increased with an increase in OM content, and an appropriate increase in OM content could improve the soil water storage capacity. This study identified a positive correlation for the relationship between soil moisture and OM, but a low R value indicated that such a relationship might not exist. The low value of R was primarily attributed to the low OM of the surface soil samplings in this study area. The relationship between soil moisture and the saturated soil hydraulic conductivity is shown in Figure 4f. The fitted linear curve for soil moisture and saturated soil hydraulic conductivity was y = 53.9477 − 2.2049x, with R = −0.5803. This revealed that the content of soil water was negatively correlated with the soil's saturated hydraulic conductivity. A higher saturated soil hydraulic conductivity (Ksat) resulted in a weakening of the soil's ability to retain water. In contrast, the soil had a strong water retention capacity. This result is supported by a previous study [46]. In addition, the soil particle size had an important influence on the Ksat. Soil particles with larger diameters were not conducive to the formation of capillary pores, since the content of fine particles, such as clay, silt, and OM, was low, the ability to fill the space was weak, and the fractal dimension of the soil was small. Moreover, the sandy soil was better at infiltrating, which was not conducive to the retention of water.
The relationship between soil moisture and topography was also examined in this study. As shown in Figure 2, the soil moisture content of the TS was obviously lower than those of MS and BS. The soil moisture increased with a decrease in slope height. In addition, the soil moisture of BS1 was lower than that of BS2 ( Figure 2). There were essentially two reasons for this. On the one hand, the elevation of BS2 was lower compared with that of BS1 (Figure 1), which constituted the bottom of slope. This further proved that the soil moisture content increased with a decrease in slope. Alternatively, the content of sand was higher in BS1, and the content of clay was higher in BS2 (Table 1), which explained why the soil moisture was higher at BS2. This phenomenon is supported by a previous study [46], but contrasts with the findings of another study [24]. This difference is primarily owing to the higher content of sand in the middle and bottom of the slope and the greater lichen cover near the top. A study on the spatial variability of soil moisture and its influence factor on the Loess Plateau by Liu and Shao [47] found that the soil moisture content from the top to the bottom of a slope tended to increase. In addition, the soil moisture content was also influenced by the groundwater level, land use pattern, and surface vegetation [24,48,49]. The effects of different vegetation cover, land use patterns, and groundwater levels on the content of soil moisture in the study area merit further study.
influence factor on the Loess Plateau by Liu and Shao [47] found that the soil moisture content from the top to the bottom of a slope tended to increase. In addition, the soil moisture content was also influenced by the groundwater level, land use pattern, and surface vegetation [24,48,49]. The effects of different vegetation cover, land use patterns, and groundwater levels on the content of soil moisture in the study area merit further study.

Soil Infiltration Process
Soil infiltration is the process by which water on the soil surface penetrates the soil. The quantification of soil infiltrability is very important for determining the components of hydrological modeling, irrigation design, and many other processes [50]. The results of soil infiltration ( Figure 5) show that there were some differences during the process of soil infiltration at different layers. The slope of the curve indicates the infiltration rate of the soil. As a whole, the lowest steady infiltration rate was 25.8 mm h -1 (MS2, 0-20 cm). In most cases, the steady infiltration rates were >40 mm h -1 . These results are similar to those obtained by Ceballos [10] in a study on the soil water behavior of sandy soils under semiarid conditions in the Duero Basin. Ceballos [10] found significant relationships between the soil texture and steady infiltration rate. Typically, the rate of steady infiltration of the soil in which the sand content is >80% and the clay content is <10% is always >40 mm h -1 . In this study, the initial infiltration rate of TS1 for the first layer was much higher than

Soil Infiltration Process
Soil infiltration is the process by which water on the soil surface penetrates the soil. The quantification of soil infiltrability is very important for determining the components of hydrological modeling, irrigation design, and many other processes [50]. The results of soil infiltration ( Figure 5) show that there were some differences during the process of soil infiltration at different layers. The slope of the curve indicates the infiltration rate of the soil. As a whole, the lowest steady infiltration rate was 25.8 mm h -1 (MS2, 0-20 cm). In most cases, the steady infiltration rates were >40 mm h -1 . These results are similar to those obtained by Ceballos [10] in a study on the soil water behavior of sandy soils under semiarid conditions in the Duero Basin. Ceballos [10] found significant relationships between the soil texture and steady infiltration rate. Typically, the rate of steady infiltration of the soil in which the sand content is >80% and the clay content is <10% is always >40 mm h -1 . In this study, the initial infiltration rate of TS1 for the first layer was much higher than that of the other points measured. The steady infiltration rate of TS2 was the highest. The initial rate of infiltration of the 20-40 cm layer and the steady rate of infiltration of the top slope were higher than those of the middle and bottom slopes. The rate of initial infiltration and the steady infiltration rate of MS2 were the lowest in the 40-60 cm layer. that of the other points measured. The steady infiltration rate of TS2 was the highest. The initial rate of infiltration of the 20-40 cm layer and the steady rate of infiltration of the top slope were higher than those of the middle and bottom slopes. The rate of initial infiltration and the steady infiltration rate of MS2 were the lowest in the 40-60 cm layer. There are many factors that affect the soil infiltration characteristics, such as soil texture, particularly the content of sand; the gravel fraction and rock fragments in the soil; the soil use pattern; and the vegetation type [10,51,52]. The highest initial infiltration rate and steady infiltration rate were observed at the top of the slope ( Figure 5). Field observations found almost no gravel or anthropogenic interference in the top slope soil. In the analysis of soil infiltration, the highest content of sand was observed in the plots with the highest initial infiltration rates and steady infiltration rates (TS1, TS2, and BS1), where the mean content of sand in these plots was 80.85%, 83.43%, and 84.42%, respectively ( Table  1). The lowest initial infiltration rate and steady infiltration rate were found in plot MS2 ( Figure 5), where the average contents of sand and clay were 75.80% and 12.15%, respectively ( Table 1). The coarse pores of sand provide a natural condition for the production of preferential flow in the soil, which is conducive to the infiltration of water. Thus, higher initial infiltration rates and steady infiltration rates were related to a higher content of sand in this study. Similarly, a study by Vries and Chow [53] found that soil macropores were the primary reason for the rapid production of interflow. Under the same hydraulic conductivity, the water moved more quickly owing to the soil macropores, and the infiltration depth was deeper [53]. Meek [54] found that greater infiltration rates probably resulted from the flow of water through macropores in a study of infiltration rates as affected by an alfalfa (Medicago sativa) and no-till cotton (Gossypium hirsutum) cropping system. that of the other points measured. The steady infiltration rate of TS2 was the highest. The initial rate of infiltration of the 20-40 cm layer and the steady rate of infiltration of the top slope were higher than those of the middle and bottom slopes. The rate of initial infiltration and the steady infiltration rate of MS2 were the lowest in the 40-60 cm layer. There are many factors that affect the soil infiltration characteristics, such as soil texture, particularly the content of sand; the gravel fraction and rock fragments in the soil; the soil use pattern; and the vegetation type [10,51,52]. The highest initial infiltration rate and steady infiltration rate were observed at the top of the slope ( Figure 5). Field observations found almost no gravel or anthropogenic interference in the top slope soil. In the analysis of soil infiltration, the highest content of sand was observed in the plots with the highest initial infiltration rates and steady infiltration rates (TS1, TS2, and BS1), where the mean content of sand in these plots was 80.85%, 83.43%, and 84.42%, respectively ( Table  1). The lowest initial infiltration rate and steady infiltration rate were found in plot MS2 ( Figure 5), where the average contents of sand and clay were 75.80% and 12.15%, respectively ( Table 1). The coarse pores of sand provide a natural condition for the production of preferential flow in the soil, which is conducive to the infiltration of water. Thus, higher initial infiltration rates and steady infiltration rates were related to a higher content of sand in this study. Similarly, a study by Vries and Chow [53] found that soil macropores were the primary reason for the rapid production of interflow. Under the same hydraulic conductivity, the water moved more quickly owing to the soil macropores, and the infiltration depth was deeper [53]. Meek [54] found that greater infiltration rates probably resulted from the flow of water through macropores in a study of infiltration rates as affected by an alfalfa (Medicago sativa) and no-till cotton (Gossypium hirsutum) cropping system.

Available Soil Water
) TS1; ( that of the other points measured. The steady infiltration rate of TS2 was the highest. The initial rate of infiltration of the 20-40 cm layer and the steady rate of infiltration of the top slope were higher than those of the middle and bottom slopes. The rate of initial infiltration and the steady infiltration rate of MS2 were the lowest in the 40-60 cm layer. There are many factors that affect the soil infiltration characteristics, such as soil texture, particularly the content of sand; the gravel fraction and rock fragments in the soil; the soil use pattern; and the vegetation type [10,51,52]. The highest initial infiltration rate and steady infiltration rate were observed at the top of the slope ( Figure 5). Field observations found almost no gravel or anthropogenic interference in the top slope soil. In the analysis of soil infiltration, the highest content of sand was observed in the plots with the highest initial infiltration rates and steady infiltration rates (TS1, TS2, and BS1), where the mean content of sand in these plots was 80.85%, 83.43%, and 84.42%, respectively ( Table  1). The lowest initial infiltration rate and steady infiltration rate were found in plot MS2 ( Figure 5), where the average contents of sand and clay were 75.80% and 12.15%, respectively ( Table 1). The coarse pores of sand provide a natural condition for the production of preferential flow in the soil, which is conducive to the infiltration of water. Thus, higher initial infiltration rates and steady infiltration rates were related to a higher content of sand in this study. Similarly, a study by Vries and Chow [53] found that soil macropores were the primary reason for the rapid production of interflow. Under the same hydraulic conductivity, the water moved more quickly owing to the soil macropores, and the infiltration depth was deeper [53]. Meek [54] found that greater infiltration rates probably resulted from the flow of water through macropores in a study of infiltration rates as affected by an alfalfa (Medicago sativa) and no-till cotton (Gossypium hirsutum) cropping system.

Available Soil Water
) TS2; ( that of the other points measured. The steady infiltration rate of TS2 was the highest. The initial rate of infiltration of the 20-40 cm layer and the steady rate of infiltration of the top slope were higher than those of the middle and bottom slopes. The rate of initial infiltration and the steady infiltration rate of MS2 were the lowest in the 40-60 cm layer. There are many factors that affect the soil infiltration characteristics, such as soil texture, particularly the content of sand; the gravel fraction and rock fragments in the soil; the soil use pattern; and the vegetation type [10,51,52]. The highest initial infiltration rate and steady infiltration rate were observed at the top of the slope ( Figure 5). Field observations found almost no gravel or anthropogenic interference in the top slope soil. In the analysis of soil infiltration, the highest content of sand was observed in the plots with the highest initial infiltration rates and steady infiltration rates (TS1, TS2, and BS1), where the mean content of sand in these plots was 80.85%, 83.43%, and 84.42%, respectively ( Table  1). The lowest initial infiltration rate and steady infiltration rate were found in plot MS2 ( Figure 5), where the average contents of sand and clay were 75.80% and 12.15%, respectively ( Table 1). The coarse pores of sand provide a natural condition for the production of preferential flow in the soil, which is conducive to the infiltration of water. Thus, higher initial infiltration rates and steady infiltration rates were related to a higher content of sand in this study. Similarly, a study by Vries and Chow [53] found that soil macropores were the primary reason for the rapid production of interflow. Under the same hydraulic conductivity, the water moved more quickly owing to the soil macropores, and the infiltration depth was deeper [53]. Meek [54] found that greater infiltration rates probably resulted from the flow of water through macropores in a study of infiltration rates as affected by an alfalfa (Medicago sativa) and no-till cotton (Gossypium hirsutum) cropping system.

Available Soil Water
that of the other points measured. The steady infiltration rate of TS2 was the highest. The initial rate of infiltration of the 20-40 cm layer and the steady rate of infiltration of the top slope were higher than those of the middle and bottom slopes. The rate of initial infiltration and the steady infiltration rate of MS2 were the lowest in the 40-60 cm layer. There are many factors that affect the soil infiltration characteristics, such as soil texture, particularly the content of sand; the gravel fraction and rock fragments in the soil; the soil use pattern; and the vegetation type [10,51,52]. The highest initial infiltration rate and steady infiltration rate were observed at the top of the slope ( Figure 5). Field observations found almost no gravel or anthropogenic interference in the top slope soil. In the analysis of soil infiltration, the highest content of sand was observed in the plots with the highest initial infiltration rates and steady infiltration rates (TS1, TS2, and BS1), where the mean content of sand in these plots was 80.85%, 83.43%, and 84.42%, respectively ( Table  1). The lowest initial infiltration rate and steady infiltration rate were found in plot MS2 ( Figure 5), where the average contents of sand and clay were 75.80% and 12.15%, respectively ( Table 1). The coarse pores of sand provide a natural condition for the production of preferential flow in the soil, which is conducive to the infiltration of water. Thus, higher initial infiltration rates and steady infiltration rates were related to a higher content of sand in this study. Similarly, a study by Vries and Chow [53] found that soil macropores were the primary reason for the rapid production of interflow. Under the same hydraulic conductivity, the water moved more quickly owing to the soil macropores, and the infiltration depth was deeper [53]. Meek [54] found that greater infiltration rates probably resulted from the flow of water through macropores in a study of infiltration rates as affected by an alfalfa (Medicago sativa) and no-till cotton (Gossypium hirsutum) cropping system.

Available Soil Water
) MS2; ( that of the other points measured. The steady infiltration rate of TS2 was the highest. The initial rate of infiltration of the 20-40 cm layer and the steady rate of infiltration of the top slope were higher than those of the middle and bottom slopes. The rate of initial infiltration and the steady infiltration rate of MS2 were the lowest in the 40-60 cm layer. There are many factors that affect the soil infiltration characteristics, such as soil texture, particularly the content of sand; the gravel fraction and rock fragments in the soil; the soil use pattern; and the vegetation type [10,51,52]. The highest initial infiltration rate and steady infiltration rate were observed at the top of the slope ( Figure 5). Field observations found almost no gravel or anthropogenic interference in the top slope soil. In the analysis of soil infiltration, the highest content of sand was observed in the plots with the highest initial infiltration rates and steady infiltration rates (TS1, TS2, and BS1), where the mean content of sand in these plots was 80.85%, 83.43%, and 84.42%, respectively ( Table  1). The lowest initial infiltration rate and steady infiltration rate were found in plot MS2 ( Figure 5), where the average contents of sand and clay were 75.80% and 12.15%, respectively ( Table 1). The coarse pores of sand provide a natural condition for the production of preferential flow in the soil, which is conducive to the infiltration of water. Thus, higher initial infiltration rates and steady infiltration rates were related to a higher content of sand in this study. Similarly, a study by Vries and Chow [53] found that soil macropores were the primary reason for the rapid production of interflow. Under the same hydraulic conductivity, the water moved more quickly owing to the soil macropores, and the infiltration depth was deeper [53]. Meek [54] found that greater infiltration rates probably resulted from the flow of water through macropores in a study of infiltration rates as affected by an alfalfa (Medicago sativa) and no-till cotton (Gossypium hirsutum) cropping system.

Available Soil Water
) BS1; and ( that of the other points measured. The steady infiltration rate of TS2 was the highest. The initial rate of infiltration of the 20-40 cm layer and the steady rate of infiltration of the top slope were higher than those of the middle and bottom slopes. The rate of initial infiltration and the steady infiltration rate of MS2 were the lowest in the 40-60 cm layer. There are many factors that affect the soil infiltration characteristics, such as soil texture, particularly the content of sand; the gravel fraction and rock fragments in the soil; the soil use pattern; and the vegetation type [10,51,52]. The highest initial infiltration rate and steady infiltration rate were observed at the top of the slope ( Figure 5). Field observations found almost no gravel or anthropogenic interference in the top slope soil. In the analysis of soil infiltration, the highest content of sand was observed in the plots with the highest initial infiltration rates and steady infiltration rates (TS1, TS2, and BS1), where the mean content of sand in these plots was 80.85%, 83.43%, and 84.42%, respectively ( Table  1). The lowest initial infiltration rate and steady infiltration rate were found in plot MS2 ( Figure 5), where the average contents of sand and clay were 75.80% and 12.15%, respectively ( Table 1). The coarse pores of sand provide a natural condition for the production of preferential flow in the soil, which is conducive to the infiltration of water. Thus, higher initial infiltration rates and steady infiltration rates were related to a higher content of sand in this study. Similarly, a study by Vries and Chow [53] found that soil macropores were the primary reason for the rapid production of interflow. Under the same hydraulic conductivity, the water moved more quickly owing to the soil macropores, and the infiltration depth was deeper [53]. Meek [54] found that greater infiltration rates probably resulted from the flow of water through macropores in a study of infiltration rates as affected by an alfalfa (Medicago sativa) and no-till cotton (Gossypium hirsutum) cropping system.

Available Soil Water
There are many factors that affect the soil infiltration characteristics, such as soil texture, particularly the content of sand; the gravel fraction and rock fragments in the soil; the soil use pattern; and the vegetation type [10,51,52]. The highest initial infiltration rate and steady infiltration rate were observed at the top of the slope ( Figure 5). Field observations found almost no gravel or anthropogenic interference in the top slope soil. In the analysis of soil infiltration, the highest content of sand was observed in the plots with the highest initial infiltration rates and steady infiltration rates (TS1, TS2, and BS1), where the mean content of sand in these plots was 80.85%, 83.43%, and 84.42%, respectively ( Table 1). The lowest initial infiltration rate and steady infiltration rate were found in plot MS2 ( Figure 5), where the average contents of sand and clay were 75.80% and 12.15%, respectively ( Table 1). The coarse pores of sand provide a natural condition for the production of preferential flow in the soil, which is conducive to the infiltration of water. Thus, higher initial infiltration rates and steady infiltration rates were related to a higher content of sand in this study. Similarly, a study by Vries and Chow [53] found that soil macropores were the primary reason for the rapid production of interflow. Under the same hydraulic conductivity, the water moved more quickly owing to the soil macropores, and the infiltration depth was deeper [53]. Meek [54] found that greater infiltration rates probably resulted from the flow of water through macropores in a study of infiltration rates as affected by an alfalfa (Medicago sativa) and no-till cotton (Gossypium hirsutum) cropping system.

Available Soil Water
The soil water retention curve, which reflects the relationship between the soil moisture content and the energy and condition of the soil to retain water, is an important curve of soil water characteristics [55,56]. The fitted results are shown in Table 2 and Figure 6. The soil water retention curve reflects the calculation of the water content at FC, WP, and AWC for the sampling sites ( Table 2). The AWC increased from the TS to BS. The average FC was 25.81% with a variation coefficient of 10.50%. The WP was more variable, and varied by approximately 6.81%. The average AWC was 19.0%. These results are similar to those reported by Zhou [57], who showed that the WP of the northwest arid area of northwest China and the Tibetan Plateau was primarily between 5.5 and 8.3%, and the FC was < 32.0%. The AWC increased as the fraction of sand decreased and the fraction of clay increased. This phenomenon is supported by the research of Ceballos et al. [10] and Fares and Alva [58]. Ceballos et al. [10] obtained a highly significant correlation coefficient between the AWC and the clay fraction (R 2 = 0.77; p = 0.005). The relationship between AWC and the soil particle composition, particularly the clay fraction, explains why a higher AWC was found in the 40-60 cm layer below the soil surface and at the BS (owing to the migration of fine soil particles) ( Table 1). A highly significant correlation between the AWC and the OM content was obtained in some studies [10,59]. However, this relationship was not obvious in this study, owing to the low contents of soil OM at each sampling point.

Soil Water Retention Capacity
To estimate the soil water retention capacity of sandy soils under natural conditions, a soil water retention curve was obtained using Equation (2). The soil water retention curve reflects the soil water retention capacity. A higher curve indicates a stronger retention capacity, and a lower curve indicates a weaker retention capacity [43].
The soil water retention curve for the six sampling points is shown in Figure 6. The highest soil water retention capacity was in the 40-60 cm layer, where the decrease in the soil water content with the water potential was the slowest. The weakest soil water retention capacity was found in the 0-20 cm layer, where the decrease in the soil water content with the water potential was the fastest. Under the same suction conditions, the soil moisture content at the different layers of the soil varied. The 40-60 cm soil layer maintained the maximum amount of soil water. Correspondingly, the amount of soil water maintained by the 0-20 cm layer was the lowest. It is also apparent in Figure 6 that the soil water retention curves of MS1, MS2, and BS2 were higher than those of TS1, TS2, and BS1, which indicates that the water retention capacity of MS1, MS2, and BS2 was higher in the study area. This was primarily related to the lower soil hydraulic conductivity, and the lower contents of sand and higher contents of clay in MS1, MS2, and BS2 ( Table 1). The soil texture is another factor that affects the soil water retention capacity [27,[60][61][62][63]. A study by Franco-Vizcaíno [64] on the water regime in soils and plants along an aridity gradient in central Baja California showed that the soil water in sandy soils was depleted much more quickly than that in loamy soils.
The soil water retention curve for the six sampling points is shown in Figure 6. The highest soil water retention capacity was in the 40-60 cm layer, where the decrease in the soil water content with the water potential was the slowest. The weakest soil water retention capacity was found in the 0-20 cm layer, where the decrease in the soil water content with the water potential was the fastest. Under the same suction conditions, the soil moisture content at the different layers of the soil varied. The 40-60 cm soil layer maintained the maximum amount of soil water. Correspondingly, the amount of soil water maintained by the 0-20 cm layer was the lowest. It is also apparent in Figure 6 that the soil water retention curves of MS1, MS2, and BS2 were higher than those of TS1, TS2, and BS1, which indicates that the water retention capacity of MS1, MS2, and BS2 was higher in the study area. This was primarily related to the lower soil hydraulic conductivity, and the lower contents of sand and higher contents of clay in MS1, MS2, and BS2 ( Table 1). The soil texture is another factor that affects the soil water retention capacity [27,[60][61][62][63]. A study by Franco-Vizcaıńo [64] on the water regime in soils and plants along an aridity gradient in central Baja California showed that the soil water in sandy soils was depleted much more quickly than that in loamy soils. A rainfall experiment was used to study the retention characteristics of soil water. The soil moisture was measured during and after rainfall. The mean soil moisture content was 12.24% at the beginning of rainfall. After 22.2 mm of rainfall in approximately 1 h, the A rainfall experiment was used to study the retention characteristics of soil water. The soil moisture was measured during and after rainfall. The mean soil moisture content was 12.24% at the beginning of rainfall. After 22.2 mm of rainfall in approximately 1 h, the soil moisture content reached its maximum and decreased with time. The mean increase in the soil moisture content was 199.12% at 1 h, which decreased to 79.64% 48 h later (Table 3). At 48 h after the end of the rainfall, the soil water contents of TS1, TS2, MS1, MS2, and BS1 were less than the corresponding soil water contents at field capacity, indicating that the soil lost water at a faster rate than that of the water supplied in these sites. However, the soil water content of BS2 was higher than the soil water content at FC, revealing that the soil water retention capacity of BS2 was higher than that of the other sampling points. The soil water distribution pattern along the profile during the rainfall simulation experiment is shown in Figure 7. The water content of the whole profile ranged from 12.24% to 14.92% before the rainfall experiment, which was lower than the FC. One hour after the end of the experiment, the soil water content in the 0-20 and 20-40 cm layers below the soil surface exceeded the water content at FC, which increased to 36.61% and 29.70%, respectively. However, the 40-60 cm layer below the soil surface only showed small changes in the soil water content. At 48 h after the end of the rainfall, the soil water content of the 0-20 and 20-40 cm layers below the soil surface decreased to 21.99% and 24.62%, respectively. The relatively higher soil hydraulic conductivity of the 0-20 and 20-40 cm soil layers (Table 1) explained the faster flux of the wetting front at these layers. Similarly, the lower soil hydraulic conductivity of the 40-60 cm soil was the primary reason of the slower flux of the wetting front at this layer [10]. However, 48 h after the end of the rainfall, the soil water content of the whole profile was ≤FC, which reveals that the soil water retention capacity was low throughout the whole profile. During the experiment, the study area experienced a rainfall intensity of 22 mm in an hour. If the influence of temperature during this period is ignored, the rainfall would substantially increase the moisture content of the surface soil. In fact, owing to the high summer temperatures in this area and the sunny days after the rain, there is no doubt that a substantial amount of the surface water evaporated, which weakened the supplementary effect of rainfall on the soil water to some extent. Under the conditions of sufficient sunlight and strong solar radiation, with an increase in light intensity and sunshine hours, the evaporation of surface soil water will increase, which will further reduce the content of soil water. Therefore, during the experiment, particularly within 48 h after rainfall infiltration, the evaporation effect of soil water loss cannot be ignored. However, considering that the study area was located on a small hillside, the evaporation rates of all the measured points were considered almost the same, and the influence of soil water evaporation on the change in soil water content was almost the same. The soil water evaporation at the monitoring points was not measured during the experiment. The relationship between evaporation intensity and soil water loss will be studied in more detail in subsequent research.

Conclusions
This study examined the soil water behavior of sandy soils under semiarid conditions in the Shendong mining area (China). The values of soil moisture at different positions of the slope were in the following order: BS > MS > TS. The soil moisture of the first 20 cm of

Conclusions
This study examined the soil water behavior of sandy soils under semiarid conditions in the Shendong mining area (China). The values of soil moisture at different positions of the slope were in the following order: BS > MS > TS. The soil moisture of the first 20 cm of soil was the lowest, and the moisture levels increased with depth. These results can be explained by the distribution of the soil's clay fraction, owing to its positive influence on the soil moisture in the case studied. In most cases, the steady state infiltration rates were >40 mm h -1 . Owing to the higher sand content and soil macropores at the TS and 0-20 cm surface soil, the initial infiltration rate and steady infiltration rate were higher. The average AWC was 18.28%, which was consistent with the predominance of sandy textural fractions. The results of a soil water retention curve and a rainfall simulation experiment revealed a low soil water retention capacity throughout the whole profile. Several considerations on the soil water behavior of sandy soils have been reported in this study that will enable an understanding of the heterogeneity in soil water behavior in these arid ecosystems and provide data for optimal land management.