Quantitative Estimation of Soil-Ground Water Storage Utilization during the Crop Growing Season in Arid Regions with Shallow Water Table Depth

: Water storage in unsaturated and saturated zones during the crop non-growing season is one of the important supplementary water resources to meet crop water requirements in arid areas with shallow water table depth. It is necessary to analyze utilization of the soil-ground water storage during the crop growing season and its attribution to irrigation during the non-growing season. To facilitate the analysis, a new method based on measurements of soil moisture content and water table depth is developed. The measurements used in this study include (1) 15-year data of soil moisture content within a depth of 1 m from the land surface and water table depth measured in Jiefangzha, including its four subareas and (2) 4-year data of the same kind in Yonglian, located in arid northern China. The soil-ground water storage utilization is calculated as the di ﬀ erence of water storage between the beginning and end of the crop growing season in the whole computational soil proﬁle. The results of average soil-ground water storage utilization in Jiefangzha and its four subareas and Yonglian are 121 mm, 126 mm, 113 mm, 124 mm, 185 mm and 117 mm, and the corresponding average utilization e ﬃ ciencies in the non-growing season are 32.2%, 32.5%, 31.5%, 31.6%, 57.3% and 47.6%, respectively. Further, the water table ﬂuctuation method was used to estimate the variation in water storage. The coe ﬃ cients of soil-ground water storage utilization, soil-ground water storage utilization below 1 m soil depth and ground water utilization are deﬁned, and their average values are 0.271, 0.111 and 0.026 in Jiefangzha, respectively. Then, the contribution of soil-ground water storage utilization to actual evapotranspiration is evaluated, which are over 23.5% in Jiefangzha and Yonglian. These results indicate that the soil-ground water storage plays an important role in the ecological environment in arid areas with shallow water table depth. contribution of soil ‐ ground water storage utilization below 1 m soil depth. RW g is the contribution of ground water utilization.


Introduction
Agricultural irrigation has consumed the largest amount of fresh water in most countries over recent decades due to the rapid increase in agricultural areas [1], which may aggravate water resources scarcity and lead to a challenge for socio-economic development [2], especially in arid regions [3]. Water consumed by plants is not only supplied by irrigation water, but also by water storage in the unsaturated and saturated zones in arid areas with shallow water table depth [4,5]. The irrigation or precipitation water during the non-growing season can be stored as soil water and ground water, which when the soil moisture content has slight changes at the beginning and the end of the crop growing season. Otherwise, it may underestimate the quantity of water storage used by crops. The second method is to divide the soil profile into two zones at a certain depth. In the zone above the depth, the difference of the soil water storage at the beginning and end of the crop growing season is calculated by measured soil moisture data, since it is easy to be obtained. In the zone below the depth, the ground water utilization is obtained by using lysimeter experiments, tracer experiments, or numerical modelling approaches when there are no measured water table depth data [18,[27][28][29][30], or by the water table fluctuation method where water table depth data are available [31,32]. However, the ground water utilization obtained by this experimental method can only represent a point or a small region [3], and the numerical modelling approaches have limitations because of the lack of data [33].
The crux of using the water table fluctuation method is to determine the specific yield, which is an important hydrogeological parameter connecting flow processes in the unsaturated and saturated zones [34]. This is defined as the volume of water that an aquifer releases from or takes into storage, per unit aquifer area, per unit change in water table depth [35]. Pumping tests are usually carried out to obtain the specific yield. However, the value obtained by pumping tests is usually lower than that obtained by water balance calculations [36], which is caused by the non-instantaneous drainage of pores in the unsaturated zone. Nwankwor et al. [37] found that the specific yield obtained by the type-curve method of Boulton was 0.07, which was significantly smaller than the value of 0.25 obtained by the volume-balance method for a shallow sandy aquifer. The specific yield estimated by Malama [38] was smaller than 0.10 using the pumping test data of the Neuman model [39] at the Boise Hydro-geophysical Research Site in Boise, Idaho (US), while the value was 0.23 when using an alternative linearization of water table kinematic condition theory with the same pumping test data at this site. It is usually arbitrary to determine the value of specific yield when calculating the storage variation of ground water. The values of 0.15, 0.07 and 0.046, were adopted by Gao et al. [40], Chen et al. [41] and Yang et al. [25], respectively, when using the water table fluctuation method for the same area, resulting in the large variation of calculated ground water utilization. The specific yield is critical to the calculation of soil-ground water storage utilization, while its determination is usually arbitrary, resulting in an inaccurate soil-ground water storage utilization. Therefore, it is important to determine the reasonable specific yield when using the water table fluctuation method.
To obtain the soil-ground water storage utilization during the crop growing season, 15 years' soil moisture content within a depth of 1 m, and the water table depth dataset for Jiefangzha and its four subareas and a four year dataset from Yonglian, both located in arid northern China, are used in this study. The impact factors of soil-ground water storage utilization are discussed for helping in water management strategy. Then, to help in using the water table fluctuation method accurately, coefficients of soil-ground water storage utilization, ground water utilization and soil-ground water storage utilization below 1 m soil depth are analyzed. Finally, the contributions of soil-ground water storage utilization, ground water utilization and soil-ground water storage utilization below 1 m soil depth to the actual evapotranspiration (ET) are evaluated.

Materials and Methods
In this section, the observed data of the study areas used for calculating the soil-ground water storage utilization is presented in Section 2.1. The conception and calculation methods for soil-ground water storage utilization are then described in Section 2.2. The moving average method used for reducing random error and eliminating the noise of the soil-ground water storage utilization is shown in Section 2.3. To obtain the actual evapotranspiration for an analysis of the contribution of soil-ground water storage utilization to actual evapotranspiration, the water balance method is described in Section 2.4.
Specifically, the soil-ground water storage utilization is calculated in Jiefangzha and Yonglian. The soil moisture content is extended to the whole computational soil profile by a linear interpolation method, and then the soil-ground water storage utilization is calculated directly by the difference of water storage in the whole computational soil profile between the beginning and end of the crop growing season. Then the soil-ground water storage utilization is divided into soil water storage utilization and ground water utilization. The soil-ground water storage utilization below 1 m soil depth is also estimated since the maximum sampling depth is usually 1 m in many agricultural areas. Next, coefficients of soil-ground water storage utilization, ground water utilization and soil-ground water storage utilization below 1 m soil depth are analyzed. Finally, the contributions of soil-ground water storage utilization, ground water utilization and soil-ground water storage utilization below 1 m soil depth to the actual evapotranspiration (ET) are calculated.

Study Area and Data Description
The Hetao Irrigation District is located in the west of the Inner Mongolia Autonomous Region, China, as shown in Figure 1. The average annual precipitation and evaporation (measured by 20 cm-diameter evaporation pans) are 152 mm and 2400 mm, respectively. The evaporation is obtained from the China Meteorological Data Service Center (http://data.cma.gov.cn). Agricultural irrigation heavily relies on surface water resources diverted from the Yellow River. The water table depth of the Hetao Irrigation District varies from 0 m to 3.5 m over space and time, with an average annual value of 2.03 m. The water table depth increases during the crop growing season in this area, because the water consumption is larger than the irrigation and precipitation. It is assumed that the study area does not have lateral hydraulic connection with surrounding regions, due to the gentle land surface and the small variation of water table depth in space [42], and that the hydrological processes (e.g., infiltration, evapotranspiration and deep percolation) mainly occur in the vertical direction. There is a freeze-thaw period from early December to late April in the Hetao Irrigation District.
Water 2020, 12, 3351 4 of 20 depth is also estimated since the maximum sampling depth is usually 1 m in many agricultural areas. Next, coefficients of soil-ground water storage utilization, ground water utilization and soil-ground water storage utilization below 1 m soil depth are analyzed. Finally, the contributions of soil-ground water storage utilization, ground water utilization and soil-ground water storage utilization below 1 m soil depth to the actual evapotranspiration (ET) are calculated.

Study Area and Data Description
The Hetao Irrigation District is located in the west of the Inner Mongolia Autonomous Region, China, as shown in Figure 1. The average annual precipitation and evaporation (measured by 20 cmdiameter evaporation pans) are 152 mm and 2400 mm, respectively. The evaporation is obtained from the China Meteorological Data Service Center (http://data.cma.gov.cn). Agricultural irrigation heavily relies on surface water resources diverted from the Yellow River. The water table depth of the Hetao Irrigation District varies from 0 m to 3.5 m over space and time, with an average annual value of 2.03 m. The water table depth increases during the crop growing season in this area, because the water consumption is larger than the irrigation and precipitation. It is assumed that the study area does not have lateral hydraulic connection with surrounding regions, due to the gentle land surface and the small variation of water table depth in space [42], and that the hydrological processes (e.g., infiltration, evapotranspiration and deep percolation) mainly occur in the vertical direction. There is a freeze-thaw period from early December to late April in the Hetao Irrigation District. There are two experimental areas in the district, the Jiefangzha sub-district (JFZSD) and the Yonglian irrigation area (YIA). The Jiefangzha sub-district can be divided into four subareas, Wula, Yangjia, Huangji and Qinghui, by the trunk diversion canals, as shown in Figure 1a. The area of Jiefangzha is 2287 km 2 . The major soil textures in Jiefangzha sub-district are silt loam and sandy loam [9], and the saturated moisture contents are close in space. Therefore, the soil moisture content is treated as the same. The saturated soil moisture content is 0.46 cm 3 /cm 3 in Jiefangzha as measured by the Water Resources Research Institute of Inner Mongolia. Wheat and maize are the main crops, and their growing seasons are from late March to late July and from the end of April to the end of September, respectively. The number of ground water and soil moisture content observation sites in each study area are shown in Table 1. There are 57 ground water observation wells and 22 soil moisture content observation sites distributed across Jiefangzha, as shown in Figure 1a. All soil moisture content observation sites are adjacent to the ground water observation wells. The major crops are sown after the freeze-thaw period after the end of April in Hetao Irrigation District. Only wheat is sown at the end of March, and consumes little water, which can be ignored. So the end of April is chosen as the beginning of crop growing season in this study. The soil moisture content and water table depth at the end of the crop growing season are observed at the end of September. Five soil samples are collected at one monitoring site at the depths of 0~0.1 m, 0.1~0.2 m, 0.2~0.4 m, 0.4~0.7 There are two experimental areas in the district, the Jiefangzha sub-district (JFZSD) and the Yonglian irrigation area (YIA). The Jiefangzha sub-district can be divided into four subareas, Wula, Yangjia, Huangji and Qinghui, by the trunk diversion canals, as shown in Figure 1a. The area of Jiefangzha is 2287 km 2 . The major soil textures in Jiefangzha sub-district are silt loam and sandy loam [9], and the saturated moisture contents are close in space. Therefore, the soil moisture content is treated as the same. The saturated soil moisture content is 0.46 cm 3 /cm 3 in Jiefangzha as measured by the Water Resources Research Institute of Inner Mongolia. Wheat and maize are the main crops, and their growing seasons are from late March to late July and from the end of April to the end of September, respectively. The number of ground water and soil moisture content observation sites in each study area are shown in Table 1. There are 57 ground water observation wells and 22 soil moisture content observation sites distributed across Jiefangzha, as shown in Figure 1a. All soil moisture content observation sites are adjacent to the ground water observation wells. The major crops are sown after the  Figure 2b,d. The water diversion during the crop growing season and non-growing season (which is also referred to as autumn irrigation) of Jiefangzha and its four subareas from 2003 to 2017 was provided by the Hetao Administration as shown in Figure 2e,g, respectively. Five weather stations are distributed in the Hetao Irrigation District from west to east. There are small differences of the meteorological data in the five weather stations when we analyze the data from 1961 to 2013. The previous study also showed that the representative distance of one weather station can be longer than 44 km in Hetao Irrigation District [43]. Therefore, one weather station can represent the meteorological condition of the whole area in Jiefangzha (2287 km 2 ). The precipitation data for Jiefangzha during the crop growing season and non-growing season comes from the Hangjinghouqi meteorological station (shown in Figure 1b), which can be found in Figure 2f,h, respectively.  Figure 2e,g, respectively. Five weather stations are distributed in the Hetao Irrigation District from west to east. There are small differences of the meteorological data in the five weather stations when we analyze the data from 1961 to 2013. The previous study also showed that the representative distance of one weather station can be longer than 44 km in Hetao Irrigation District [43]. Therefore, one weather station can represent the meteorological condition of the whole area in Jiefangzha (2287 km 2 ). The precipitation data for Jiefangzha during the crop growing season and non-growing season comes from the Hangjinghouqi meteorological station (shown in Figure 1b), which can be found in Figure 2f,h, respectively.   The location of Yonglian irrigation area can be found in Figure 1c. The area of Yonglian is 29.7 km 2 . The major crop is sunflower, the growing season of which is from the end of May to late September. There are 10 soil moisture content observation sites in Yonglian, which are adjacent to the 10 ground water observation wells, as shown in Figure 1c Figure 2e,g, respectively. The precipitation data during the crop growing season and non-growing season comes from the Yonglian experiment station and the Wuyuan meteorological station (shown in Figure 1b), which can be found in Figure 2f,h, respectively. Figure 3a shows the typical soil moisture profile at the beginning and end of crop growing season in the study area. The soil moisture profile at the beginning of the crop growing season is higher due to the water storage caused by the autumn irrigation. The water storage is consumed by transpiration and evaporation during the crop growing season, which leads to a lower soil moisture profile after harvest. The difference in water storage between the beginning and end of the crop growing season is regarded as the soil-ground water storage utilization during the crop growing season, which can be calculated by the following equation,

Calculation of Soil-Ground Water Storage Utilization
where ∆W is the soil-ground water storage utilization, mm; and W 1 and W 2 are the water storages at the beginning and end of the crop growing season, mm.  The per unit area water storage W of the vertical soil profile at a specific time is calculated as, where Zb is the maximum computational soil profile depth, mm; θ is the soil moisture content at the soil depth z, cm 3 cm −3 ; θi is the soil moisture content of layer i, cm 3 cm −3 ; di is the thickness of the layer i, mm; and m is the number of calculation layers. The Zb value is set as 3.5 m in this study since the maximum water table depth in Jiefangzha and Yonglian during the crop growing season is shallower than 3.5 m. The θi value is the average soil water content of the soil moisture content monitoring sites in the study area.
The vertical soil profile, as shown in Figure 3b, is divided into three zones from top to bottom for calculating the water storage W. In the first zone, the calculation uses measured soil moisture content data, the second zone uses soil water content data obtained by a linear interpolation method, and the third zone below the water table depth uses the saturated soil moisture content. The soilground water storage utilization (marked as ΔW) can be divided into two parts by the initial water table depth, which are the soil water storage utilization marked as ΔWs and the ground water utilization ΔWg as shown in Figure 3c. Besides, the soil-ground water storage utilization below 1 m soil depth (marked as ΔW1) is estimated as shown in Figure 3d, since the maximum sampling depth is usually 1 m in many agricultural areas and researchers usually consider ΔW1 as the ground water utilization [31]. To help in using the water table fluctuation method accurately, three coefficients are defined as follows for estimating the soil-ground water storage utilization, the ground water utilization and the soil-ground water storage utilization below 1 m soil depth: Figure 3. The schematic diagrams of (a) the soil moisture profiles at the beginning and end of the crop growing season, (b) the three zones of the vertical soil profile for calculating the water storage (marked as W), (c) the soil water storage utilization (marked as ∆W s ) and the ground water utilization (marked as ∆W g ) and (d) the soil-ground water storage utilization below 1 m soil depth (marked as ∆W 1 ) and above 1 m soil depth (marked as ∆W 01 ). Note: ∆W is the soil-ground water storage utilization, mm. W 2 is the water storages at end of the crop growing season. wtd is the water table depth, mm. wtd 1 and wtd 2 are the water table depths at the beginning and end of the crop growing season, mm. ∆h is the variation of water table depth during the crop growing season, mm. Z s is the maximum sampling depth, mm. Z b is the maximum computational depth, mm.
The per unit area water storage W of the vertical soil profile at a specific time is calculated as, where Z b is the maximum computational soil profile depth, mm; θ is the soil moisture content at the soil depth z, cm 3 cm −3 ; θ i is the soil moisture content of layer i, cm 3 cm −3 ; d i is the thickness of the layer i, mm; and m is the number of calculation layers. The Z b value is set as 3.5 m in this study since the maximum water table depth in Jiefangzha and Yonglian during the crop growing season is shallower than 3.5 m. The θ i value is the average soil water content of the soil moisture content monitoring sites in the study area.
The vertical soil profile, as shown in Figure 3b, is divided into three zones from top to bottom for calculating the water storage W. In the first zone, the calculation uses measured soil moisture content data, the second zone uses soil water content data obtained by a linear interpolation method, and the third zone below the water table depth uses the saturated soil moisture content. The soil-ground water storage utilization (marked as ∆W) can be divided into two parts by the initial water table depth, which are the soil water storage utilization marked as ∆W s and the ground water utilization ∆W g as shown in Figure 3c. Besides, the soil-ground water storage utilization below 1 m soil depth (marked as ∆W 1 ) is estimated as shown in Figure 3d, since the maximum sampling depth is usually 1 m in many agricultural areas and researchers usually consider ∆W 1 as the ground water utilization [31]. To help in using the water table fluctuation method accurately, three coefficients are defined as follows for estimating the soil-ground water storage utilization, the ground water utilization and the soil-ground water storage utilization below 1 m soil depth: where η is the coefficient of soil-ground water storage utilization, η g is the coefficient of ground water utilization, which is used to quantify the variation of water storage in the ground water fluctuation zone with the water table depth during the crop growing season, η 1 is the coefficient of soil-ground water storage utilization below 1 m soil depth, and ∆h is the increase of water table depth during the crop growing season, mm. Equations (3)-(5) are applied to arid and semi-arid regions where the water table depth at the end of the crop growing season is larger than that at the beginning of the crop growing season. Since the soil and ground water storage at the beginning of the crop growing season is mainly from the autumn irrigation in the study area, a coefficient referred to as the utilization efficiency of autumn irrigation (λ) is proposed to quantify the contribution of the autumn irrigation to the water consumption of crops in the next year, where Q a is the amount of autumn irrigation, mm.

Moving Average Method
The moving average method is used to reduce random errors and eliminate the noise of temporal soil-ground water storage utilization [44]. For time series data x, the moving average value is calculated as follows [45], where f j is the moving average soil-ground water storage utilization, mm; K is the length of memory window, which is expressed as K = 2k + 1 and k is an integer. In this study, K is adopted as 5 years. x i is the original time series of soil-ground water storage utilization, mm, and N is the length/number of the time series x.

Water Balance Analysis to Estimate the Actual Evapotranspiration
The water balance method is used to estimate the actual evapotranspiration during the crop growing season as follows [46], where ET is the actual evapotranspiration in the field during the crop growing season, mm; I s is the amount of field irrigation water during the crop growing season, mm; and P s is the precipitation during the crop growing season, mm. The contributions of soil-ground water storage utilization, soil-ground water storage utilization below 1 m soil depth and ground water utilization to the actual evapotranspiration (marked as RW, RW 1 and RW g ) can be calculated as, These contributions (RW, RW 1 and RW g ) can quantitatively distinguish the contribution of the soil-ground water storage utilization within different depths to the actual evapotranspiration. The soil-ground water storage utilization, soil-ground water storage utilization below 1 m soil depth and ground water utilization play different roles in actual evapotranspiration. So it is helpful to calculate Water 2020, 12, 3351 9 of 19 different contributions for determining the contribution of every example of water storage utilization to the actual evapotranspiration.

Statistical Indicators
Classical statistical indicators, including average, maximum, minimum, standard deviation and coefficient of variation, are used to characterize the features of soil-ground water storage utilization and other variables. The coefficient of variation (C v ) is used to explain the extent of the variability of variables, which can be calculated by [47], where S x is the standard deviation of the dataset, and X is the average value of the dataset.

Results of the Soil-Ground Water Storage Utilization
The soil-ground water storage utilizations at Jiefangzha and its four subareas (2003~2017) and Yonglian (2005, 2006, 2007 and 2009) during the crop growing season are shown in Figure 4a,b, and the statistical results are listed in Table 2. The average soil-ground water storage utilization of Jiefangzha in 2003~2017 is 121 mm, varying from 57 mm to 229 mm and the standard deviation is 49 mm. The average soil-ground water storage utilizations in Jiefangzha, Wula, Yangjia, Huangji, Qinghui and Yonglian are 121 mm, 126 mm, 113 mm, 124 mm, 185 mm, and 117 mm, respectively. The Qinghui subarea has the largest average soil-ground water storage utilization (185 mm), standard deviation (86 mm) and variation range (307 mm). One reason is that the Qinghui subarea obtained smaller water diversion and precipitation during the crop growing season, the average value of which is 357 mm from 2003 to 2017, smaller than those in Jiefangzha (381 mm) and other subareas. The average utilization efficiencies of autumn irrigation (λ) in Jiefangzha, Wula, Yangjia, Huangji, Qinghui and Yonglian are 32.2%, 32.5%, 31.5%, 31.6%, 57.3% and 47.6% (shown in Table 2), respectively.
Water 2020, 12, 3351 9 of 20 different contributions for determining the contribution of every example of water storage utilization to the actual evapotranspiration.

Statistical Indicators
Classical statistical indicators, including average, maximum, minimum, standard deviation and coefficient of variation, are used to characterize the features of soil-ground water storage utilization and other variables. The coefficient of variation (Cv) is used to explain the extent of the variability of variables, which can be calculated by [47], where Sx is the standard deviation of the dataset, and is the average value of the dataset.

Results of the Soil-Ground Water Storage Utilization
The soil-ground water storage utilizations at Jiefangzha and its four subareas (2003~2017) and Yonglian (2005, 2006, 2007 and 2009) during the crop growing season are shown in Figure 4a,b, and the statistical results are listed in Table 2. The average soil-ground water storage utilization of Jiefangzha in 2003~2017 is 121 mm, varying from 57 mm to 229 mm and the standard deviation is 49 mm. The average soil-ground water storage utilizations in Jiefangzha, Wula, Yangjia, Huangji, Qinghui and Yonglian are 121 mm, 126 mm, 113 mm, 124 mm, 185 mm, and 117 mm, respectively. The Qinghui subarea has the largest average soil-ground water storage utilization (185 mm), standard deviation (86 mm) and variation range (307 mm). One reason is that the Qinghui subarea obtained smaller water diversion and precipitation during the crop growing season, the average value of which is 357 mm from 2003 to 2017, smaller than those in Jiefangzha (381 mm) and other subareas. The average utilization efficiencies of autumn irrigation (λ) in Jiefangzha, Wula, Yangjia, Huangji, Qinghui and Yonglian are 32.2%, 32.5%, 31.5%, 31.6%, 57.3% and 47.6% (shown in Table 2), respectively.    To analyze the temporal variation of soil-ground water storage utilization, only the soil-ground water storage utilization in Jiefangzha and its four subareas from 2003 to 2017 are used since the data series of Yonglian is insufficient. The moving average soil-ground water storage utilization in Jiefangzha and its four subareas is shown in Figure 5a, in which an obvious increasing trend can be found. The soil-ground water storage utilization is increased by 7.00 mm, 7.33 mm, 12.01 mm, 4.85 mm, and 13.32 mm per year in Jiefangzha, Wula, Yangjia, Huangji and Qinghui, respectively. The data of Jiefangzha is taken as an example to analyze the reasons for the increasing trend. The water storage at the beginning of the crop growing season increases by 3.8 mm per year, contributing 54.6% to the increase of the soil-ground water storage utilization, and the water storage at the end of the crop growing season decreases by 3.2 mm per year, contributing 45.4% to the increase of the soil-ground water storage utilization, as shown in Figure 5b. To analyze the temporal variation of soil-ground water storage utilization, only the soil-ground water storage utilization in Jiefangzha and its four subareas from 2003 to 2017 are used since the data series of Yonglian is insufficient. The moving average soil-ground water storage utilization in Jiefangzha and its four subareas is shown in Figure 5a, in which an obvious increasing trend can be found. The soil-ground water storage utilization is increased by 7.00 mm, 7.33 mm, 12.01 mm, 4.85 mm, and 13.32 mm per year in Jiefangzha, Wula, Yangjia, Huangji and Qinghui, respectively. The data of Jiefangzha is taken as an example to analyze the reasons for the increasing trend. The water storage at the beginning of the crop growing season increases by 3.8 mm per year, contributing 54.6% to the increase of the soil-ground water storage utilization, and the water storage at the end of the crop growing season decreases by 3.2 mm per year, contributing 45.4% to the increase of the soilground water storage utilization, as shown in Figure 5b.

The Impact Factors on Soil-Ground Water Storage Utilization
The correlations between the soil-ground water storage utilization in Jiefangzha and the eight impact factors are listed in Table 3. It can be found that there are five impact factors that have significant correlation with the soil-ground water storage utilization at the significance level of 0.01 or 0.05. The factors are the autumn irrigation time (Tai), the area of the spring irrigation (Si), the area

The Impact Factors on Soil-Ground Water Storage Utilization
The correlations between the soil-ground water storage utilization in Jiefangzha and the eight impact factors are listed in Table 3. It can be found that there are five impact factors that have significant correlation with the soil-ground water storage utilization at the significance level of 0.01 or 0.05. The factors are the autumn irrigation time (T ai ), the area of the spring irrigation (S i ), the area of summer crops (S s ), the area of autumn crops (S a ) and the increase of water table depth during the crop growing season (∆h). The autumn irrigation time has a significant positive relationship with the soil-ground water storage utilization (R = 0.647, P = 0.009). The autumn irrigation occurs in 2003~2007 and 2013~2017 from 23 September~31 October and 3 October~9 November as shown in Figure 5c. The delayed autumn irrigation time leads to the higher soil moisture content at the beginning of the next crop growing season, so the water storage at the beginning of crop growing season would increase causing the increase of soil-ground water storage. The area of the spring irrigation has a positive relationship with the soil-ground water storage utilization (R = 0.641, P = 0.010). The spring irrigation area increases from 42 km 2 in 2003~2007 to 497 km 2 in 2013~2017 as shown in Figure 5d, and the larger area of the spring irrigation leads to the larger water storage at the beginning of crop growing season. The area of the summer crops has a significant negative correlation with the soil-ground water storage utilization (R = −0.572, P = 0.026), while there is a significant positive correlation between the area of the autumn crops and the soil-ground water storage utilization (R = 0.604, P = 0.017). As shown in Figure 5e,f, the area of summer crops is decreasing and the area of autumn crops is increasing in Jiefangzha, causing the increase of water consumption in August and September and the decrease of water storage at the end of the crop growing season. This indicates that the adjustment of planting structure can significantly impact the soil-ground water storage utilization. There is a significant positive relationship between the increase of water table depth during the crop growing period and soil-ground water storage utilization (R = 0.554, P = 0.032), which also can be found in Figure 6a. Table 3. Correlation between the soil-ground water storage utilization and the eight impact factors. Note: (1) T ai is the autumn irrigation time. S i is the area of the spring irrigation, mm. S s is the area of summer crops, km 2 . S a is the area of autumn crops, km 2 . ∆h is the increase of water table depth during the crop growing season, mm. I a + P a is the total water supply during the non-growing season, mm; I a is the field irrigation water during the non-growing season, mm; P a is the precipitation during the non-growing season, mm. I s + P s is the total water supply during the crop growing season, mm; I s is the field irrigation water during the crop growing season, mm; P s is the precipitation during the crop growing season, mm. h is the average water table depth during the crop growing season, mm. (2) * Correlation is significant at the 0.05 level (2-tailed), and ** Correlation is significant at the 0.01 level (2-tailed).

Impact Factors
The correlation of the other three impact factors (the total water supply during the non-growing season I a + P a , the total water supply during the crop growing season I s + P s and the average water table depth during the crop growing season h) with the soil-ground water storage utilization are not significant at the 0.05 level. It can be found that there is a positive relationship between the total water supply during the non-growing season and soil-ground water storage utilization (R = 0.479, P = 0.071). When the total water supply during the non-growing season increases, the water storage at the beginning of crop growing season will increase. However, as shown in Figure 6b, the total water supply during the non-growing season is very stable at 200 mm-210 mm from 2006 to 2014, which caused the poor correlation with the increase of soil-ground water storage utilization in these years. The total water supply during the crop growing season has a negative relationship with the soil-ground water storage utilization (R = −0.403, P = 0.136).The linear negative correlation is obvious when the total water supply during the crop growing season is larger than 400 mm, close to the crop water requirement of major crops in this area, which means that crops can obtain sufficient water by irrigation and precipitation causing less soil-ground water storage utilization, as shown in Figure 6c. In the works of Chen et al. [31] and Liu et al. [3], there is also a negative relationship between the ground water contribution to evapotranspiration and the total water supply during the crop growing season. The positive relationship between the average water table depth during the crop growing season h and the soil-ground water storage utilization is not significant, with R = 0.333 and P = 0.225. The relationship between the soil-ground water storage utilization and the average water table depth is quadratic as shown in Figure 6d. A critical water table depth (2080 mm) exists, before which the relationship between the soil-ground water storage utilization and the average water table depth is positive, and then it becomes negative. The critical water table depth was also found by Chen et al. [41]. The existence of the critical water table depth indicates that the soil-ground water storage utilization would increase first, and then decrease in the future with a decrease of water diversion from Yellow river, which is considered to increase the water table depth [48,49].
water storage at the end of the crop growing season. This indicates that the adjustment of planting structure can significantly impact the soil-ground water storage utilization. There is a significant positive relationship between the increase of water table depth during the crop growing period and soil-ground water storage utilization (R = 0.554, P = 0.032), which also can be found in Figure 6a. Table 3. Correlation between the soil-ground water storage utilization and the eight impact factors. Note: (1) Tai is the autumn irrigation time. Si is the area of the spring irrigation, mm. Ss is the area of summer crops, km 2 . Sa is the area of autumn crops, km 2 . Δh is the increase of water table depth during the crop growing season, mm. Ia + Pa is the total water supply during the non-growing season, mm; Ia is the field irrigation water during the non-growing season, mm; Pa is the precipitation during the non-growing season, mm. Is + Ps is the total water supply during the crop growing season, mm; Is is the field irrigation water during the crop growing season, mm; Ps is the precipitation during the crop growing season, mm. ℎ is the average water table depth during the crop growing season, mm. (2) * Correlation is significant at the 0.05 level (2-tailed), and ** Correlation is significant at the 0.01 level (2-tailed). In general, the soil-ground water storage utilization is significantly affected by the autumn irrigation time, the area of the spring irrigation, the area of summer crops, the area of autumn crops, and the variation of water table depth during the crop growing season. Changing the autumn irrigation scheme and planting structure would be effective ways to increase soil-ground water storage utilization. When the water supply during the crop growing season is insufficient, delaying the autumn irrigation time, increasing the amount of autumn irrigation and increasing the area of spring irrigation are useful strategies to help increase the soil-ground water storage before the crop growing period. Meanwhile, our suggestions are to adjust the planting structure, increase the area of autumn crops, and reduce the area of summer crops, which can help crops use more soil-ground water resources for water consumption during the crop growing season.

Results of Soil Water Storage Utilization, Ground Water Utilization, Soil-Ground Water Storage Utilization below and above 1 m Soil Depth and the Three Coefficients
In this section, the two components of soil-ground water storage utilization are estimated, including the soil water storage utilization ∆W s and the ground water utilization ∆W g , and the soil-ground water storage utilization below 1 m soil depth ∆W 1 and above 1 m soil depth ∆W 01 . The water table depth at the end of the crop growing season is larger than that at the beginning of the crop growing season in Jiefangzha and Yonglian, and the ∆h is larger than 0, so Equations (3)-(5) can be used in this research. Then, the three coefficients, the coefficient of soil-ground water storage utilization η, the coefficient of ground water utilization η g and the coefficient of soil-ground water storage utilization below 1 m soil depth η 1 , are calculated. The results of ∆W s and ∆W g and their contributions to the soil-ground water storage utilization ∆W in Jiefangzha and Yonglian are shown in Figure 7, the statistical results of which are listed in Table 4. The average soil water storage utilization ∆W s is 112 mm, which accounts for 92.6% of the soil-ground water storage utilization ∆W, while the ground water utilization ∆W g is 9 mm (7.4% of the ∆W) in Jiefangzha. In Yonglian, the average ∆W s is 98 mm, which accounts for 83.0% of the ∆W, and the average ∆W g is 19 mm (17.0 % of the ∆W). The contribution of soil water storage utilization during the crop growing season is larger than 80% in this study, which proves that the soil water storage utilization is the main part of the soil-ground water storage utilization. The results illustrate that the variation in water storage in the ground water fluctuation zone during the crop growing season is small. due to the capillary siphoning phenomenon, which leads to slight soil moisture content variation near the saturated zone. The ground water utilization ∆W g has a significant positive relationship (R = 0.961, P = 0.000) with the increase of water table depth during the crop growing season in Jiefangzha, as shown in Figure 8a. The coefficient of ground water utilization η g is shown in Figure 9, the statistical results of which are listed in Table 4. The coefficient of ground water utilization η g ranges from 0.008 to 0.04 with the standard deviation of 0.01, and the average value is 0.026 in Jiefangzha. In Yonglian, the average coefficient of ground water utilization is 0.028. In the works of Xu et al. [9] and Yue et al. [50], the coefficients of ground water utilization were 0.039 and 0.086 for estimating the variation of the ground water storage in Jiefangzha and Yichang, both of which are located in the Hetao Irrigation District. They chose these values according to the specific yield based on the hydrogeological studies. However, it can be found that the adopted values in the works of Xu et al. [9] and Yue et al. [50] would overestimate the variation of ground water storage. Therefore, in calculating the variation of ground water storage, the coefficient of ground water utilization should be carefully selected, which is smaller than the specific yield.
Water 2020, 12, 3351 13 of 20 in Figure 7, the statistical results of which are listed in Table 4. The average soil water storage utilization ΔWs is 112 mm, which accounts for 92.6% of the soil-ground water storage utilization ΔW, while the ground water utilization ΔWg is 9 mm (7.4% of the ΔW) in Jiefangzha. In Yonglian, the average ΔWs is 98 mm, which accounts for 83.0% of the ΔW, and the average ΔWg is 19 mm (17.0 % of the ΔW). The contribution of soil water storage utilization during the crop growing season is larger than 80% in this study, which proves that the soil water storage utilization is the main part of the soilground water storage utilization. The results illustrate that the variation in water storage in the ground water fluctuation zone during the crop growing season is small. due to the capillary siphoning phenomenon, which leads to slight soil moisture content variation near the saturated zone. The ground water utilization ΔWg has a significant positive relationship (R = 0.961, P = 0.000) with the increase of water table depth during the crop growing season in Jiefangzha, as shown in Figure  8a. The coefficient of ground water utilization ηg is shown in Figure 9, the statistical results of which are listed in Table 4. The coefficient of ground water utilization ηg ranges from 0.008 to 0.04 with the standard deviation of 0.01, and the average value is 0.026 in Jiefangzha. In Yonglian, the average coefficient of ground water utilization is 0.028. In the works of Xu et al. [9] and Yue et al. [50], the coefficients of ground water utilization were 0.039 and 0.086 for estimating the variation of the ground water storage in Jiefangzha and Yichang, both of which are located in the Hetao Irrigation District. They chose these values according to the specific yield based on the hydrogeological studies. However, it can be found that the adopted values in the works of Xu et al. [9] and Yue et al. [50] would overestimate the variation of ground water storage. Therefore, in calculating the variation of ground water storage, the coefficient of ground water utilization should be carefully selected, which is smaller than the specific yield.    Note: ∆W s is the soil water storage utilization, mm. ∆W g is the ground water utilization, mm. ∆W 01 is the soil-ground water storage utilization above 1 m soil depth, mm. ∆W 1 is the soil-ground water storage utilization below 1 m soil depth, mm. η g is the coefficient of ground water utilization. η 1 is the coefficient of soil-ground water storage utilization below 1 m soil depth. η is the coefficient of soil-ground water storage utilization. C v is the coefficient of variation. Note: ΔWs is the soil water storage utilization, mm. ΔWg is the ground water utilization, mm. ΔW01 is the soil-ground water storage utilization above 1 m soil depth, mm. ΔW1 is the soil-ground water storage utilization below 1 m soil depth, mm. ηg is the coefficient of ground water utilization. η1 is the coefficient of soil-ground water storage utilization below 1 m soil depth. η is the coefficient of soilground water storage utilization. Cv is the coefficient of variation.  The results of average soil-ground water storage utilization above and below 1 m soil depth (ΔW01 and ΔW1) and their contributions to soil-ground water storage utilization ΔW in Jiefangzha and Yonglian are shown in Figure 7, the statistical results of which are listed in Table 4. The average ΔW01 is 74 mm, accounting for 60.6% of the soil-ground water storage utilization ΔW, and the ΔW1 is 47 mm, accounting for 39.4% of ΔW in Jiefangzha. In Yonglian, the average ΔW01 and ΔW1 are 74 mm and 44 mm, accounting for 62.5% and 37.5% of ΔW. The contribution of soil-ground water storage utilization above 1 m soil depth is larger than 60%, which indicates that the utilization of soil-ground water storage is mainly from the soil depth above 1 m. There is a significant positive relationship (R = 0.842, P = 0.000) between the soil-ground water utilization below 1 m soil depth ΔW1 and the  Note: ΔWs is the soil water storage utilization, mm. ΔWg is the ground water utilization, mm. ΔW01 is the soil-ground water storage utilization above 1 m soil depth, mm. ΔW1 is the soil-ground water storage utilization below 1 m soil depth, mm. ηg is the coefficient of ground water utilization. η1 is the coefficient of soil-ground water storage utilization below 1 m soil depth. η is the coefficient of soilground water storage utilization. Cv is the coefficient of variation.  The results of average soil-ground water storage utilization above and below 1 m soil depth (ΔW01 and ΔW1) and their contributions to soil-ground water storage utilization ΔW in Jiefangzha and Yonglian are shown in Figure 7, the statistical results of which are listed in Table 4. The average ΔW01 is 74 mm, accounting for 60.6% of the soil-ground water storage utilization ΔW, and the ΔW1 is 47 mm, accounting for 39.4% of ΔW in Jiefangzha. In Yonglian, the average ΔW01 and ΔW1 are 74 mm and 44 mm, accounting for 62.5% and 37.5% of ΔW. The contribution of soil-ground water storage utilization above 1 m soil depth is larger than 60%, which indicates that the utilization of soil-ground water storage is mainly from the soil depth above 1 m. There is a significant positive relationship (R = 0.842, P = 0.000) between the soil-ground water utilization below 1 m soil depth ΔW1 and the The results of average soil-ground water storage utilization above and below 1 m soil depth (∆W 01 and ∆W 1 ) and their contributions to soil-ground water storage utilization ∆W in Jiefangzha and Yonglian are shown in Figure 7, the statistical results of which are listed in Table 4. The average ∆W 01 is 74 mm, accounting for 60.6% of the soil-ground water storage utilization ∆W, and the ∆W 1 is 47 mm, accounting for 39.4% of ∆W in Jiefangzha. In Yonglian, the average ∆W 01 and ∆W 1 are 74 mm and 44 mm, accounting for 62.5% and 37.5% of ∆W. The contribution of soil-ground water storage utilization above 1 m soil depth is larger than 60%, which indicates that the utilization of soil-ground water storage is mainly from the soil depth above 1 m. There is a significant positive relationship (R = 0.842, P = 0.000) between the soil-ground water utilization below 1 m soil depth ∆W 1 and the increase of the water table depth during the crop growing season in Jiefangzha, as shown in Figure 8b. The coefficient of soil-ground water storage utilization below 1 m soil depth η 1 is shown in Figure 9, the statistical results of which are displayed in Table 4. The coefficient of the soil-ground water storage utilization below 1 m soil depth η 1 ranges from 0.069 to 0.152 with the C v value of 0.243 and the average value is 0.111 in Jiefangzha. In Yonglian, the average η 1 is 0.065. The maximum sampling depth is usually 1 m in many agricultural areas. So this coefficient is an important parameter to estimate the variation of the water storage resource below 1 m soil depth. The value used for calculating the soil-ground water storage utilization below 1 m soil depth by Chen et al. [41] and Liu et al. [3] was 0.07 in Jiefangzha. It can be found that the calculated soil-ground water storage utilization below 1 m soil depth would induce a 37% error.
The coefficient of soil-ground water storage utilization η is shown in Figure 9, the statistical results of which are listed in Table 4. The coefficient of the soil-ground water storage utilization η ranges from 0.104 to 0.497 with a standard deviation of 0.114. The coefficient of the soil-ground water storage utilization η varies greatly over the years. The average coefficients of soil-ground water storage utilization are 0.271 and 0.172 in Jiefangzha and Yonglian, respectively. This value is consistent with the work of Steinwand et al. [51], where the coefficients were 0.18 for the sandy site and 0.28 for the clay loam site. This value is much larger than the specific yield used for calculating the variation of the soil-ground water storage during the crop growing season, which is obtained by pumping tests due to the non-instantaneous drainage in the vadose zone and aquifer compressibility [36,52,53]. Therefore, the specific yield from the pumping test cannot be used directly for calculation of the soil-ground water storage utilization, which could underestimate the results.
In general, the three coefficients are ranked as η > η 1 > η g , and there is a great difference among them. So one should be careful to choose suitable coefficients when calculating the variation of soil-ground water storage.  Table 5. The average calculated actual evapotranspiration during the crop growing season is 502 mm from 2003 to 2017 in Jiefangzha, ranging from 413 mm to 591 mm. In Yonglian, the average calculated actual evapotranspiration is 485 mm from 2005 to 2009, ranging from 431 mm to 561 mm. The contributions of soil-ground water storage utilization (RW), soil-ground water storage utilization below 1 m soil depth (RW 1 ) and ground water utilization (RW g ) to the actual evapotranspiration in Jiefangzha and Yonglian are shown in Figure 10b,d, the statistical results of which are listed in Table 5. In Jiefangzha, the average contribution of soil-ground water storage utilization to the actual evapotranspiration (RW) is 23.7%, which demonstrates that the soil-ground water storage utilization is an important source for crop water consumption. The average contribution of soil-ground water storage utilization below 1 m soil depth to the actual evapotranspiration (RW 1 ) is 9.2% and the average contribution of ground water utilization to the actual evapotranspiration (RW g ) is 1.8% in Jiefangzha. The RW in Jiefangzha has an increasing trend from 2013 to 2017, which increases from 19.1% in 2003~2007 to 29.3% in 2013~2017. The RW has a significant negative relationship with precipitation during the crop growing season (R = −0.621, P = 0.013). Due to the larger precipitation during the crop growing season in 2012 and 2014, the RW of 2012 and 2014 are smaller than that in 2011~2017. In Yonglian, the average RW, RW 1 and RW g are 24.2%, 9.0% and 4.0%, respectively. It can be found that the average contribution of soil-ground water storage utilization below 1 m soil depth to the actual evapotranspiration is over 9.0% and the contribution of soil-ground water storage utilization to the actual evapotranspiration is over 23.5%. The results of Gao et al. [40] are consistent with the results of this study; the average contribution of soil-ground water storage utilization below 0.5 m soil depth to the actual evapotranspiration were 25.6%, 12.0% and 15.2% in upper, middle and lower reaches of Hetao Irrigation District, respectively, which is within the range of RW 1 (9.0%) to RW (23.5%) in this study. These results indicate that storing water resources during the non-growing season by irrigation or precipitation is a useful way to alleviate the water shortage during the crop growing season in arid agricultural regions with shallow water table depth.   Note: ET is the actual evapotranspiration, mm. RW is the contribution of soil-ground water storage utilization. RW1 is the contribution of soil-ground water storage utilization below 1 m soil depth. RWg is the contribution of ground water utilization.

Conclusions
In this study, a long-term soil moisture content and water table depth dataset are used to estimate soil-ground water storage utilization during the crop growing season and its attribution to irrigation of the non-growing season in arid agricultural areas with shallow water table depth; coefficients of soil-ground water storage utilization, ground water utilization and soil-ground water storage utilization below 1 m soil depth are also analyzed. The contributions of soil-ground water storage utilization, soil-ground water storage utilization below 1 m soil depth and ground water utilization to actual evapotranspiration are further evaluated. The major conclusions are as follows: (1) The soil-ground water storage is critical to crop water consumption during the crop growing season, which accounts for over 23.5% of the actual evapotranspiration.
(2) More than 30% of the autumn irrigation water during the non-growing season can be used as soil-ground water utilization in arid areas with shallow water table depth.
(3) Changing the autumn irrigation scheme and the planting structure can significantly impact the soil-ground water storage utilization.  Note: ET is the actual evapotranspiration, mm. RW is the contribution of soil-ground water storage utilization. RW 1 is the contribution of soil-ground water storage utilization below 1 m soil depth. RW g is the contribution of ground water utilization.

Conclusions
In this study, a long-term soil moisture content and water table depth dataset are used to estimate soil-ground water storage utilization during the crop growing season and its attribution to irrigation of the non-growing season in arid agricultural areas with shallow water table depth; coefficients of soil-ground water storage utilization, ground water utilization and soil-ground water storage utilization below 1 m soil depth are also analyzed. The contributions of soil-ground water storage utilization, soil-ground water storage utilization below 1 m soil depth and ground water utilization to actual evapotranspiration are further evaluated. The major conclusions are as follows: (1) The soil-ground water storage is critical to crop water consumption during the crop growing season, which accounts for over 23.5% of the actual evapotranspiration.
(2) More than 30% of the autumn irrigation water during the non-growing season can be used as soil-ground water utilization in arid areas with shallow water table depth.
(3) Changing the autumn irrigation scheme and the planting structure can significantly impact the soil-ground water storage utilization.
(4) The utilization of soil-ground water storage is mainly from the soil depth above 1 m, which contributes more than 60% of the soil-ground water storage utilization.
(5) For using the water table fluctuation method to estimate the variation of the soil-ground water storage, one should be careful in choosing the suitable coefficient. The commonly used specific yield obtained by pumping tests would overestimate the ground water utilization and underestimate the soil-ground water utilization.

Conflicts of Interest:
The authors declare no conflict of interest.