Evaluation of Soil Water Availability (swa) Based on Hydrological Modelling in Arid and Semi-arid Areas: a Case Study in Handan City, China

To evaluate the spatial and temporal soil water availability (SWA) and improve the soil water use efficiency in agriculture, this study established a new index system for SWA evaluation considering the soil storage capacity, the temporal and spatial matching degree between soil water supply and crop water demand, and the soil water transformation and utilization efficiency. To quantitatively calculate the indexes, a self-developed hydrological model (MODCYCLE) was introduced and applied in Handan City. The results show that MODCYCLE performed well in simulating the soil hydrological processes. A comparison of the observed soil moisture data and the simulated values revealed relative errors below 10% and correlation coefficients all above 0.7. According to the evaluation results, the SWA for different typical years was very unevenly distributed in Handan City. The average SWA in a normal year is 0.71, which is higher than that in dry years at 0.62 and in wet years at 0.65. Regarding different regions, the average SWA was generally much higher in the eastern plain area than that in the western mountainous area. Generally, there is a great potential for improving SWA and further the sustainable development of soil water resource, particularly for dryland farming in similar areas.


Introduction
Influenced by both climate change and high-intensity human activities, the volume of runoff water resources is currently decreasing sharply [1].Under this background, water use competition among different sectors (e.g., urban, industry, agriculture, and ecology) is increasing.The agriculture sector is known to be vulnerable compared with other sectors that have better economic benefits [2].Therefore, agricultural water resources are largely diverted to other sectors to pursue higher benefits.This practice is definitely a great threat to grain production and food security.To guarantee sustainable agricultural development and to promote food security, an increasing number of studies have focused on soil water resources, particularly in arid and semi-arid areas.In general, water can be directly absorbed by crops and plants after it is transformed into soil moisture.Therefore, soil water, as a type of water resource, plays an important role in agricultural production and food security.However, agricultural specialists and water resource managers have traditionally paid much attention to the "visible" water, e.g., surface water and ground water, which is also known as blue water.Owing to the difficulty of monitoring and regulating soil water, the utilization value of soil water resources is largely ignored, particularly in Chinese agriculture.Recently, soil water monitoring and simulation technologies have developed significantly, which has provided a good foundation for the efficient use and management of soil water resources in rain-fed agriculture, particularly in arid and semiarid areas.
Soil moisture, known as a direct index of soil water availability (SWA), is an important indicator for farmland hydrological forecasts and irrigation scheduling.Therefore, most previous studies have focused on the monitoring and forecasting of soil moisture by using different techniques.Robinson et al. [3] concluded that there are two main approaches for studying the spatial-temporal distribution of soil moisture.The first method refers mainly to measurement methods by direct contact with the soil profile, such as gravimetric methods [4], tensiometers [5], the neutron probe method [6], and time domain reflectometry (TDR) [7,8].Tanaka and Aase [9] used the neutron scattering technique to measure soil moisture at the field scale to determine the season in which the fallow method might influence soil water storage.Zegelin et al. [10] improved the field probes to detect the soil water content based on TDR, which has been widely used owing to its convenience and economic efficiency.Gao et al. [11] studied the spatial variability of the root zone soil moisture (RZSM) in gullies by using the observed soil moisture data collected by the TDR method from nearby uplands.The soil moisture data collected by the traditional direct methods are accurate but are time-consuming and cost intensive if the sampling area is large [12].
The second method, which refers mainly to remote sensing techniques [13], cannot directly measure the RZSM.However, it can sense the RZSM without direct contact with the soil based on multi-spectral information from the soil layer.Remote sensing methods include passive microwave radiometers, synthetic aperture radars, scatterometers, and thermal methods.These methods are either ground based or operated from airborne or spaceborne platforms [14].Recent progress in remote sensing technologies has allowed us to obtain frequent soil moisture distribution maps easily and quickly [15,16].Kong and Dorling [17] retrieved near-surface soil moisture from advanced synthetic aperture radar imagery, which showed satisfactory accuracy for the shallow soil layer profile.Ford et al. [18] used an exponential decay filter to estimate RZSM using the near-surface soil moisture derived from the Soil Moisture and Ocean Salinity (SMOS) satellite.Martínez-Fernández et al. [19] calculated the soil water deficit index (SWDI) based on SMOS satellite remote sensing data for agricultural drought monitoring.The advantage of remote sensing is that it gives uniform spatial and temporal coverage of the estimated soil moisture easily and quickly.However, the limitations include confinement to the top few centimetres of the soil layer and coarse spatial and temporal resolutions.However, when integrated with a land surface process model, remote sensing can estimate the soil moisture both at the surface and deep in the root zone with high spatial and temporal resolutions.In addition, data assimilation provides a means for combining the advantages of modelling with those of remote sensing data.A host of studies based on remote sensing measurement have used data assimilation to derive an indirect estimate of the RZSM [20].Montaldo and Albertson [21] developed a multi-scale assimilation approach for RZSM estimation, which is shown to provide marked improvements.Reichle and Koster [22] assimilated the Scanning Multi-channel Microwave Radiometer (SMMR)-derived real surface soil moisture data into the National Aeronautics and Space Administration (NASA) Catchment land surface model at a large continental scale over a long period.The soil moisture estimates in this research are significantly better than those obtained with the model or from the satellite alone.Moreover, using statistical and distributed hydrological models to simulate the spatial and temporal pattern of soil water has become an attractive method because it is more time saving and cost effective than traditional methods.Kim [23] developed a time series model for soil moisture simulation and the results showed that the simulated values fit well with the observed soil moisture series.Gao et al. [24] mapped the spatial and temporal distribution of soil water moisture based on a self-developed hydrological model in Northern China.The physical-based distributed hydrological model showed good reliability for soil moisture simulation.
The studies mentioned above focus mainly on soil water observation, soil moisture mapping, and agricultural drought monitoring.Few scientists have sought to improve soil water use efficiency based on soil water content observation and simulation for agricultural water management.In the context of climate change and strong anthropogenic influence on the hydrological cycle, the available volume of blue water resources is shrinking significantly.Therefore, improving the utilization efficiency of green water resources and making full use of soil water resources in agriculture are of great concern, particularly in arid and semi-arid areas.It is widely accepted that evaluation of soil water resources availability is the first and most important step for high-efficiency soil water resources development and sustainable utilization.Therefore, the objective of this study is to introduce a distributed hydrological model to depict the soil water transformation and flow routes at a regional farmland scale.The soil water availability assessment index system is established to jointly consider factors such as soil storage capacity, soil water recharge and discharge rates, and soil water uptake efficiency by crops and plants.The main objectives of this study include (1) development of a spatially distributed hydrological model to simulate the soil water transformation and flow process based on the complex tillage and water management regimes of farmland in the North China Plain; (2) establishment of the soil water availability assessment index framework considering such parameters as soil properties, soil water reservoir characteristics, and soil water recharge and discharge efficiency; and (3) selection of Handan City in the North China Plain in a semi-arid area as the study region and quantitative calculation of the farmland soil water availability.In addition, we reveal the spatial and temporal variability of soil water use potential and efficiency in accordance with the evaluation results, and we provide specific and differential suggestions and strategies for different sub-regions in the study area.This study can serve as a guide for improving soil water use efficiency and maintaining sustainable agricultural and ecological development in arid and semi-arid areas.

Soil Water Simulation Tool
Before SWA evaluation is conducted, the scientist should understand the soil water movement process in the farmland that is used to quantify the water availability index.However, for a large-scale basin or region, it is impossible to explore the soil water movement process though field experiments, which are time consuming and cost intensive.With the development of computer numerical simulation technology, the soil transformation and flow process at a large scale has been efficiently simulated by hydrological models in many studies [25].In this study, the object-oriented modularised model for basin-scale water cycle simulation (MODCYCLE) is introduced to simulate the farmland hydrological processes and to map the temporal and spatial distribution of the farmland soil water moisture [24,26].MODCYCLE, self-developed by the China Institute of Water Resources and Hydropower Research (IWHR), is a spatially distributed model properly combining a hydro-physical model with a conceptual model.Moreover, MODCYCLE avoids the difficulty in parameterization and overload of computation.
MODCYCLE is similar to the commonly used soil and water assessment tool (SWAT) model [27], particularly in its hydrological modelling and calculation.Both are semi-distributed and conceptual-physical models in which the main hydrological processes such as overland runoff, groundwater runoff, evapotranspiration, soil water infiltration, and percolation are all depicted with generalised physical methods and simplified empirical formulas.Numerous previous studies have proved that SWAT is an effective tool for simulating soil moisture distribution at the regional scale, particularly in natural basins or areas of limited anthropogenic influence [28,29].However, it performs poorly when simulating the hydrological processes in areas influenced by strong human activities, especially for the complex cropping systems and irrigation regimes in China.Therefore, MODCYCLE is designed to improve the simulation accuracy under the conditions of strong anthropogenic influences and complicated Chinese tillage systems.According to the simulation principle, MODCYCLE uses the core of SWAT because the simulation principles of SWAT are proven to be effective and practical by a series of case studies and applications.Moreover, a special module, the agricultural management module, has been embedded in MODCYCLE to depict the complex hydrological processes of Chinese farmland with complicated tilling practices and strong anthropogenic influences.Based on the agricultural management module, an agricultural activity command database in accordance with the current farmland management regimes of China has been created, which controls MODCYCLE to precisely and elaborately simulate the farmland hydrological processes.In summary, compared with SWAT, MODCYCLE performs much better in depicting the farmland soil water movement processes in China.Therefore, it was used as an important tool for SWA evaluation in this study.
The structural framework of MODCYCLE is similar to that of SWAT, the most widely used hydrological model.The first step in MODCYCLE is to divide the study area into sub-basins.The sub-basins are further subdivided into different hydrological response units (HRUs).The runoff generated in each sub-basin is from two sources: the area within the sub-basin and the upstream contributing areas.Based on a digital elevation model (DEM) and GIS analysis, the topological relationships among the sub-basins are established, and the runoff flow directions in the sub-basin can be depicted.For SWA evaluation, soil water movement simulation is of extreme importance.Water molecules are absorbed in the soil mass forced by the soil matric potential in an unsaturated zone.Therefore, soil water moves much slower than surface water.However, compared with soil water movement in a horizontal direction, the speed in the vertical direction is much faster because of gravity.As a result, soil water movement in the horizontal direction in an unsaturated zone can be neglected, although that in the vertical direction must be considered.In MODCYCLE, the main hydrological process in the vertical direction is simulated through a one-dimensional (1D) Soil Column Method at the HRU scale (Figure 1).
proven to be effective and practical by a series of case studies and applications.Moreover, a special module, the agricultural management module, has been embedded in MODCYCLE to depict the complex hydrological processes of Chinese farmland with complicated tilling practices and strong anthropogenic influences.Based on the agricultural management module, an agricultural activity command database in accordance with the current farmland management regimes of China has been created, which controls MODCYCLE to precisely and elaborately simulate the farmland hydrological processes.In summary, compared with SWAT, MODCYCLE performs much better in depicting the farmland soil water movement processes in China.Therefore, it was used as an important tool for SWA evaluation in this study.
The structural framework of MODCYCLE is similar to that of SWAT, the most widely used hydrological model.The first step in MODCYCLE is to divide the study area into sub-basins.The sub-basins are further subdivided into different hydrological response units (HRUs).The runoff generated in each sub-basin is from two sources: the area within the sub-basin and the upstream contributing areas.Based on a digital elevation model (DEM) and GIS analysis, the topological relationships among the sub-basins are established, and the runoff flow directions in the sub-basin can be depicted.For SWA evaluation, soil water movement simulation is of extreme importance.Water molecules are absorbed in the soil mass forced by the soil matric potential in an unsaturated zone.Therefore, soil water moves much slower than surface water.However, compared with soil water movement in a horizontal direction, the speed in the vertical direction is much faster because of gravity.As a result, soil water movement in the horizontal direction in an unsaturated zone can be neglected, although that in the vertical direction must be considered.In MODCYCLE, the main hydrological process in the vertical direction is simulated through a one-dimensional (1D) Soil Column Method at the HRU scale (Figure 1).As shown in Figure 1, HRU is regarded as the basic simulation unit in the model.Vertical water movement consists of three processes: the surface water process, soil water process, and ground water process.Because the soil water process is the main concern in our study, the simulation method for soil water movement is introduced in detail below.In the 1D soil column, the soil layer is divided into two parts.The upper part is the schematised layer, and the lower part is the seepage transition layer, as shown in Figure 1.Ordinarily, the schematised layer is the crop root zone in farmland, where the soil water content fluctuates frequently at different timescales, particularly in As shown in Figure 1, HRU is regarded as the basic simulation unit in the model.Vertical water movement consists of three processes: the surface water process, soil water process, and ground water process.Because the soil water process is the main concern in our study, the simulation method for soil water movement is introduced in detail below.In the 1D soil column, the soil layer is divided into two parts.The upper part is the schematised layer, and the lower part is the seepage transition layer, as shown in Figure 1.Ordinarily, the schematised layer is the crop root zone in farmland, where the soil water content fluctuates frequently at different timescales, particularly in the growing season.To precisely simulate the soil water content variation process in the crop root zone, the schematised layer should be also subdivided into different sub-layers (e.g., Layer A, Layer B, Layer C, etc.).In MODCYCLE, the schematised layer can be subdivided into 10 multiple soil layers.The seepage transition layer is directly beneath the schematised layer, where the soil water content is relatively stable.However, the seepage transition layer is an important linkage zone for water exchange between the ground water and soil water.That is, the seepage transition layer is a buffer zone in which the surplus water can be stored during the wet season and is supplied into the soil root zone for crop growth during the dry season.Therefore, the seepage transition layer plays a significant role in farmland cultivation and management.In the 1D soil column method, the water balance equation can be described as where ∆SW is the change in the soil water content in the simulation period, R is the equivalent precipitation depth adjusted for canopy interception, IR is the irrigation depth in the simulation period, GEV is the evaporation depth from phreatic water to the crop root zone, SF is the surface runoff loss depth, ∆LF is the change in the surface water ponding depth, UEV is the evaporation depth from the surface water, PEV is the crop evapotranspiration depth, SEV is the evaporation depth from the soil layer, and SEP is the water deep infiltration depth from the soil layer to the groundwater aquifer; all are in millimetres.As described in Equation ( 1), the entire soil layer can be considered as a water tank.Most processes of water exchange in the soil layer occur through the top boundary (farmland surface) and the bottom boundary (ground water table).Nine variables need to be determined to solve the equation.R and IR are input parameters in the model, and the calculation methods for the rest of the parameters can be obtained from related studies [26,30].

Definition of Soil Water Availability
Improving the efficiency of agricultural water resources is the central topic for water-saving agriculture technologies.In arid and semi-arid areas, any type of water resource should be transformed into soil water before it is effectively absorbed and utilised by crops and plants.Based on the analysis of soil water resource properties, this study argues that three factors will primarily influence soil water use efficiency at the farmland scale.The first is availability of the soil water quantity, which is mainly affected by soil particles and root zone depth.Not all water stored in the soil layer can be directly absorbed by crops; essentially, only the soil moisture above the wilting point and below the field capacity can be used by crops.Here, we introduce two first-class indices to represent the soil water quantity: soil water storage capacity [31] and the operational characteristics of the farmland soil reservoir [32], the calculation methods of which will be introduced later.The second factor is the matching degree of soil water resources between the supply and demand.In some cases, although there is sufficient soil water supply throughout the year, the crops may also suffer from drought in some months (e.g., April and May) because the temporal distribution of available soil water is uneven.
Here, we introduce one first-class index to represent the matching degree of the soil water resource, which includes two sub-indices: the temporal matching degree index and the spatial matching degree index.The third factor is the efficiency of soil water transformation and uptake by plants.Here, we introduce one first-class index to represent the soil water efficiency [33,34], which also includes two sub-indices: soil water conversion efficiency and crop absorption efficiency.The details of the aforementioned indices are given in Table 1.
Previous studies have largely examined how to increase the efficiency of soil water transformation and plant absorption.However, they neglect to discuss the importance of the availability and quantity of soil water to be used by plants and crops and the degree of temporal and spatial matching degree between the soil water supply and the demand.In this study, we comprehensively consider these three factors and suggest a new concept for SWA evaluation.An optimal objective function is presented to quantitatively characterise the soil water availability at the farmland scale, as shown in Equation ( 2): where U represents the comprehensive SWA, F (X, Y, Z) is the soil water availability evaluation function, A represents the factor for the degree of difficulty of the soil water to be used by plants and crops, R represents the factor for the degree of matching between the soil water supply and the water demand, and E represents the factor for the efficiency of soil water transformation and plant uptake.

Framework of the Soil Water Availability Index System
As shown in Equation ( 2), the optimal objective function of soil water availability is only a conceptual expression and it cannot be used to quantitatively evaluate the soil water availability directly.Therefore, an index system was introduced to quantitatively characterise the three factors.Based on the index system, we used the analytic hierarchy process (AHP) method to calculate the comprehensive SWA.
The SWA index system, presented in Table 1, includes four first-class indices and seven second-class indicators.The first and second first-class indices (C1 and C2) mainly represent the factor for the degree of difficulty for the soil water to be used by plants and crops, which are further subdivided into three second-class indicators.The third first-class index (C3) represents the matching degree between the soil water supply and the water demand, which is further subdivided into two second-class indicators.The fourth first-class index (C4) represents the factor for the efficiency of the soil water transformation and plant absorption, which is further subdivided into two second-class indicators.

Factor
First-Class Index Second-Class Index Description

Availability of soil water quantity
The soil water storage capacity (C1) The valid capacity of soil reservoir (C11) The valid water supply capacity (equivalent water depth) of the soil in the crop root zone, mm.

Operational characteristics of farmland soil reservoir (C2)
Average value of annual empty storage capacity of the soil reservoir (C21) To characterise the average soil water shortage severity in the soil root zone within a year, mm.

Average annual variation of empty storage capacity of the soil reservoir (C22)
To characterise the fluctuation of water shortage severity in the root zone within a year, dimensionless.
The spatial and temporal matching of soil water resource between supply and demand Temporal and spatial matching degree (C3)

Temporal matching degree index (C31)
The temporal matching degree between soil water supply and crop water requirement, dimensionless.

Spatial matching degree index (C32)
The spatial matching degree between soil water supply and crop water requirement, dimensionless.
The efficiency of soil water transformation and plant uptake Soil water transformation and crop absorption efficiency (C4)

Soil water conversion efficiency (C41)
To characterise the efficiency of water transforming from other types to soil water, dimensionless.

Crop absorption efficiency (C42)
To characterise the efficiency the soil water is absorbed by crops and plants, dimensionless.

Index for Soil Water Storage Capacity (C1)
Soil water storage is a key parameter for farmland hydrology in arid and semiarid areas.In view of the water holding ability of the soil, scientists always regard this capacity as the soil water reservoir capacity.According to previous studies [31], soil water reservoir capacities include mainly the total soil water reservoir capacity (TC), soil water storage capacity (SC), unavailable soil water reservoir capacity (UC), available soil water reservoir capacity (AC), and soil dead water storage capacity (DC).In this study, we present a second-class index (C11) to represent the first-class index soil water storage capacity (C1).Here, C11 represents the available soil water reservoir capacity, which is essentially determined by the soil properties in the root zone of the farmland.With a total of n soil layers in the root zone, the valid capacity of the soil reservoir can be represented as where W p represents the valid storage capacity (equivalent water depth) of the soil water in the crop root zone; FC (i) represents the field soil water holding capacity (equivalent water depth) of the soil layer i; θ wp (i) is the volume of the soil water (equivalent water depth) at the wilting point of the soil layer i, and D i represents the thickness of the soil layer i, all of which are in millimetres.In addition, n represents the number of soil layers in the root zone.

Index for the Operational Characteristics of the Farmland Soil Reservoir (C2)
To characterise the average soil water shortage severity in the soil root zone within one year, the two second-class indices (C21 and C22) are introduced to represent the index for the operational characteristics of the farmland soil reservoir.The index C21 is the average value of the annual empty storage capacity of the soil reservoir: where W e is the average value of the annual empty storage capacity (equivalent water depth) of the soil reservoir, in millimetres; FC (i) represents the field soil water holding capacity (equivalent water depth) of the soil layer i, in millimetres; θ wp (i) represents the actual soil water content (equivalent water depth) of the soil layer i; D i represents the thickness of the soil layer i, in millimetres; m represents the total number of days in the evaluation period; and n represents the number of layers of soil in the root zone.The index C22 represents the average annual variation of the empty storage capacity of the soil reservoir for the root zone: where C w is the average annual variation of the empty storage capacity of the soil reservoir, in millimetres.The other terms are the same as those described in Equation ( 4).

Index of Degree of Matching between the Soil Water Supply and Crop Water Demand (C3)
In this study, the two second-class indexes (C31 and C32) were originally introduced to characterise the matching degree between the soil water supply and crop water demand.The index C31 represents the temporal matching degree between the soil water supply and crop water demand: where m ti represents the temporal matching degree between the soil water supply and crop water demand in month i in the growing period, with values ranging from 0 to 1; m t is the average temporal matching degree between the soil water supply and crop water demand in the entire growing period; and n is the number of months in the entire growing period for a specific planting structure.DE i represents the total crop water demand in month i, and IR i represents the total artificial irrigation amount in month i, both in millimetres, because the concept of temporal matching degree between soil water supply and crop water demand does not consider the influence of artificial irrigation.Therefore, the irrigation amount should be excluded to calculate the index (C31).RC i is the total available soil water supply in month i, in millimetres.When RC i = ET i + RES i , where ET i is the total water evaporated from the crop and soil to the air in month i, RES i is the effective soil water remaining in the soil layer profile at the end of month i.
The index C32 represents the spatial matching degree between the soil water supply and crop water demand: where m sj represents the spatial matching degree between the soil water supply and crop water demand in the sub-region j in the growing period, with values ranging from 0 to 1; A j is the area of the sub-region j, in square kilometres; A is the total area of the entire region, in square kilometres; and n is the number of sub-regions in the entire study region.DE j represents the total crop water demand in the sub-region jin the growing period, in millimetres; and IR j represents the total artificial irrigation amount for crop growth in the whole growing period in sub-region j.As in the temporal matching degree, the artificial irrigation amount should be excluded to calculate the index C32.RC j is the total available soil water supply in sub-region j in the entire growing period, in millimetres, which is described as RC j = ET j + RES j , where ET j is the total water evaporated from the crop and soil to the air in sub-region j in the entire growing period, in millimetres; and RES j is the effective soil water remaining in the soil layer profile in sub-region j at the end of the growing period, in millimetres.

Index for Soil Water Transformation and Crop Absorption Efficiency (C4)
The soil water transformation efficiency refers to the ratio between the effective soil water recharged in the root zone and the total volume of the other types of water (e.g., precipitation, irrigation water, and phreatic water) that can be transformed into soil water: where E ct is the daily soil water transformation efficiency during day t, which ranges from 0 to 1; E c is the average soil water transformation efficiency in the entire growing period, which also ranges from 0 to 1; T is the total number of days in the entire growing period; θ i (t + 1) is the actual soil water content in soil layer i during day (t + 1); θ i (t) is the actual soil water content in soil layer i during day t; TH i is the thickness of the soil layer i; ET (t) is the crop transpiration during day t; ES (t) is the soil water evaporation during day t; R (t) is the water loss through deep percolation during day t; NP (t) is the net precipitation adjusted for canopy interception; I (t) is the irrigation amount during day t and ε (t) is the volume of the ground water returned to the root zone through phreatic water evaporation during day t.
Based on the concept of water use efficiency (WUE) [34], this study proposes the index of the crop absorption efficiency to characterise the soil water use efficiency by crops and plants: where E dt represents the soil water absorption efficiency by crops and plants in day t and E d is the average soil water absorption efficiency in the entire growth period.The other terms are the same as those described previously.

Analytic Hierarchy Process
The AHP method was first developed in the 1970s by Saaty from the University of Pittsburgh [35].Owing to its simplicity and great flexibility, the AHP method has been widely applied in engineering, management, manufacturing, marketing, and other applications [36].Moreover, the AHP method has proven to be an efficient tool to help decision-makers achieve structured determinations, and it is used in multi-criteria evaluation for decomposing a complex and unstructured problem into its component parts and synthesizing the subjective judgment from experts into a hierarchical structure to determine the weights of different factors [37].In this method, the structure usually has three or four levels.The highest level (I) reflects the overall goal of the problem.The second and third levels (II and III) include the factors and sub-factors, also known as criteria and sub-criteria, and the fourth level (IV) contains the decision options [38].In this study, only three levels are used.The highest level is the overall goal, which is the total comprehensive soil water availability (SWA) for different regions.The second level includes the first-class indices (C1, C2, C3, and C4) from the index system, and the third level is composed of the second-class indices (C11, C21, C22, C31, C32, C41, and C42).
After the factor levels are established, the multi-criteria SWA evaluation using the AHP method can be conducted following four main steps.
First, a pair-wise comparison matrix is developed.The relative weight of importance or contribution of the factors (or sub-factors) with respect to an element in the upper level is evaluated via pair-wise comparison based on Saaty's nine-point scale [39].Saaty's nine-point scale transforms the subjective judgment of the factors' importance into numerical values.Therefore, the matrix is generated by placing the factors (or sub-factors) horizontally and vertically.The values are defined by comparing the horizontal factor with the vertical factor.On the basis of this method, the pair-wise comparison matrix of the four factors in this study (C1, C2, C3, and C4) was determined, as shown in Table 2.
Secondly, to calculate the weights and contribution of the factors to the overall goal, we normalised the matrix and calculated the priority vectors.In our study, MATLAB software was used to calculate the priority vectors.According to the calculation, the priority vectors of the pair-wise comparison matrix separately were 4.17, 0.042, −0.106 + 0.842i, and −0.106 + 0.842i.
Thirdly, the consistency was checked.The reliability of the results needs to be tested by using the consistency ratio (CR), where CR = CI/RI.CI = (λ max − n) / (n − 1), where CI is the consistency index; λ max is the maximum eigenvalue of the matrix, and λ max = 4.17 in this study; n is the number of factors; and RI is the random consistency index, which is the index of a pair-wise comparison matrix generated randomly according to the number of factors [40].Normally, if CR is higher than 0.1, the consistency check is not acceptable.In this study, CR = (4.17− 4)/(4 − 1)/0.9 = 0.063, which is <0.1.Therefore, the consistency is acceptable.
Finally, the normalised priority vector corresponding to the maximum eigenvalue of the pair-wise comparison matrix is determined as the weights of the different factors relative to the evaluation goal.According to the calculation through MATLAB, the priority vector W = (0.072, 0.183, 0.240, 0.505).Thus, U = 0.072C 1 + 0.183C 2 + 0.240C 3 + 0.505C 4 , where, U is the comprehensive soil water availability; and C 1 , C 2 , C 3 , and C 4 are the first-class SWA indices as shown in Table 1.Moreover, according to the evaluation index system given in Table  .It should be noted that the values of the second indices must be transformed to be standard values according to the requirements of the AHP method.There are two types of second-class indices.One is the positive index, which means that when the second index is larger, the first-class index is larger.The other is the negative index, which means that when the second-class index is larger, the first-class index is lower.Among the second-class indices in this study, the indices C 21 and C 22 are the negative indices, which are transformed to be standard values by the following method: where XC represents the index standard value, which ranges from 0 to 1; C represents the actual calculated value of the second-class index; X max represents the theoretical maximum value of the second-class index; and X min represents the theoretical minimum value of the second-class index.
The remaining five indices (C 11 , C 31 , C 32 , C 41 , and C 42 ) are the positive indexes, which are transformed to be standard values by the following method: The terms are all described in Equation ( 14).Based on the above analysis, the comprehensive SWA is calculated by the following equation: U = 0.07XC 11 + 0.13XC 21 + 0.06XC 22 + 0.12(XC 31 + XC 32 ) + 0.1XC 41 + 0.4XC 42 (16) where the terms are the same as those described previously.

Overview of the Study Area
Handan City is located in the southern part of the North China Plain at latitude 36.07 • N to 37.02 • N and longitude 113.47 • E to 115.47 • E, which covers an area of 12,047 km 2 (Figure 2).Handan City has a semi-arid continental monsoon climate with four distinct seasons and high spatially and temporally distributed precipitation.According to the statistics, the average annual precipitation is 548.9 mm, and the average annual temperature ranges from 12.5 • C to 14.One of the most important objectives of this research is to simulate the soil moisture in the farmland.In MODCYCLE, the Green-Ampt method is used to simulate the soil water infiltration and percolation process using land use data, plant variety, and soil type as key parameters.In this research, the land use data in the study area were retrieved by remote sensing.Five main land use patterns are included in the study area: farmland, forest land, grassland, urban areas, and rural residential areas.Farmland covers the largest area of Handan City, accounting for 86% of the total area.The crop variety data in this research were collected from the Rural Statistical Yearbook of Hebei Province combined with the field survey.For the field survey, we visited more than 100 farmers to determine the crop varieties and planting areas, particularly in Daming and Quzhou County.The survey data were used to verify the reliability of the Rural Statistical Yearbook of Hebei Province.The statistics indicate nine crop varieties in the study area, with winter wheat, multiple cropping of wheat and summer maize, cotton, vegetables, and apple orchard as the five dominant crops.The soil-type data were collected from the China Soil Scientific Database and were validated through high-density sampling experiments.Generally, 22 soil types are included in the study region.The distribution maps of the planting varieties with land use and soil type in the study area are separately shown in Figure 3a,b.
In this study, the period of the simulation was 1998-2014.A warm-up period from 1996 to 1997 was also used for this model.All of the data were collected within this period.The data from 1998 to 2000 were used for model calibration, and those from 2001 to 2003 were used for model validation.Model calibration is an important process for improving the simulation accuracy and precision.Ten parameters were selected as the most sensitive parameters.Because the MODCYCLE is a newly developed model, the auto calibration module is not yet embedded.Therefore, we adopted the manual parameter adjustment method for parameter calibration.After more than 20 attempts and simulations, the optimal set of parameters were determined and selected, as shown in Table 4.

Model Validation
In our study, we compared the simulated and soil observed water contents to validate the model.The observed soil water content data were obtained by sampling in the soil profile, where the average depth of the soil profile was usually 600 mm.Four sampling points were selected in one soil profile at depths of 0-15 mm, 15-30 mm, 30-45 mm, and 45-60 mm.After the soil samples were collected, the soil water content was measured using the oven drying method.Based on the above measurements, the average soil water content in the entire soil profile was calculated as the weighted average value of the soil water content of the four sampling points.We collected the in situ observed soil moisture data from 12 soil moisture monitoring stations evenly distributed in the study area.A list of these stations is given in Table 5.For convenience, the 12 monitoring stations were also classified into four groups.
For validation, the simulated soil water content results in the HRUs in which the soil moisture monitoring stations are located were selected for comparison with the observed soil water content from the 12 monitoring stations.The depth of the soil water content simulated in MODCYCLE was as much as 3000 mm.However, to keep the depth of the simulation results consistent with the observed results, we extracted only the simulated soil water content of the shallow soil layer (0-60 mm) to make comparisons with the observed soil moisture data.According to the simulation and observed results, comparison charts were drawn separately for the four groups of soil moisture monitoring stations, as shown in Figure 4.As shown in Figure 4, the points in which the x-and y-coordinates represent the simulated and observed values, respectively, were essentially distributed along the 1:1 lines.This demonstrates an obvious consistency between the simulated and observed results.To further illustrate and validate the simulation results, the relative error (ER) and correlation coefficient (R) were introduced to quantitatively depict the consistency between the simulated and the measured results.As shown in Table 6, the statistical results show good agreement between the simulated and observed soil water contents as well.
Water 2016, 8, 360 14 of 25 As shown in Figure 4, the points in which the x-and y-coordinates represent the simulated and observed values, respectively, were essentially distributed along the 1:1 lines.This demonstrates an obvious consistency between the simulated and observed results.To further illustrate and validate the simulation results, the relative error (ER) and correlation coefficient (R) were introduced to quantitatively depict the consistency between the simulated and the measured results.As shown in Table 6, the statistical results show good agreement between the simulated and observed soil water contents as well.As shown in the table, the relative errors (ER) were all below 10%, and the correlation coefficients (R) were all above 0.7 except for the data of Long-wangmiao and Linzhang stations.Because these two stations are located in the east and south plain area of Handan City, where the anthropogenic influence on farmland hydrological process is stronger than that in other areas, the simulation accuracy is poorer compared with other stations.This problem is also a global concern in hydrology.Generally, the simulation results were convincible and acceptable.Therefore, it is feasible to use MODCYCLE for soil water availability evaluation in the study area.As shown in the table, the relative errors (ER) were all below 10%, and the correlation coefficients (R) were all above 0.7 except for the data of Long-wangmiao and Linzhang stations.Because these two stations are located in the east and south plain area of Handan City, where the anthropogenic influence on farmland hydrological process is stronger than that in other areas, the simulation accuracy is poorer compared with other stations.This problem is also a global concern in hydrology.Generally, the simulation results were convincible and acceptable.Therefore, it is feasible to use MODCYCLE for soil water availability evaluation in the study area.

Temporal Perspective
As previously mentioned, the simulation period in this study was from 2001 to 2014.Based on the simulation results, the comprehensive soil water availability in the entire study area from 2001 to 2014 was calculated, as shown in Figure 5.   Figure 5 shows that the comprehensive SWA ranged from 59% to 76.8% in different years.The maximum value of 76.8% occurred in 2014 with yearly precipitation of 507 mm, and the minimum value of 59% occurred in 2002 with yearly precipitation of 346 mm.To analyse the relationship between soil water availability and the yearly precipitation of different hydrological level years, we ranked the SWA from 2001 to 2014 in ascending order.In addition, we determined the hydrological level of the different years according to the frequency analysis of precipitation.The resultant SWA ranking chart is shown in Figure 6.
Water 2016, 8, 360 16 of 25 between soil water availability and the yearly precipitation of different hydrological level years, we ranked the SWA from 2001 to 2014 in ascending order.In addition, we determined the hydrological level of the different years according to the frequency analysis of precipitation.The resultant SWA ranking chart is shown in Figure 6. Figure 6 shows the SWA from 2001 to 2014 sorted in ascending order.The high-soilwater-availability years were concentrated mainly in normal hydrological years.Five normal hydrological years were indicated in the simulation period, in which the SWA, shown by yellow bars in the figure, was much higher than that in dry and wet years, shown by brown and blue bars, respectively.Preliminary analysis indicated that in dry years, the water supply from precipitation was deficient, which cannot satisfy plant water requirements during the growing period.Therefore, the poor degree of matching between the soil water supply and plant water demand is the primary reason for the low SWA in dry years.In wet years, the precipitation was concentrated mostly in summer.Therefore, the available soil water distribution was extremely uneven, which in turn resulted in significant water loss in the rainy season, leaving plants vulnerable to drought in the sowing season.

Spatial Perspective
The spatial distribution of the comprehensive SWA was mapped based on the simulation results.As the first step, the original values of the second-class indices were quantitatively calculated on the basis of Equation (3) to Equation (13).In this study, three different typical hydrological-level years were selected (e.g., a wet year, normal year, and dry year) through the frequency analysis of precipitation.The years 2006, 2010, and 2012 were determined to be typical dry, normal, and wet years.As the second step, the original values of the second-class indices were transformed to the standard values as shown in Tables 7-9 according to Equations ( 14) and (15).As the final step, the comprehensive SWA in different administrative regions was calculated on the basis of Equation ( 16). Figure 6 shows the SWA from 2001 to 2014 sorted in ascending order.The high-soil-water-availability years were concentrated mainly in normal hydrological years.Five normal hydrological years were indicated in the simulation period, in which the SWA, shown by yellow bars in the figure, was much higher than that in dry and wet years, shown by brown and blue bars, respectively.Preliminary analysis indicated that in dry years, the water supply from precipitation was deficient, which cannot satisfy plant water requirements during the growing period.Therefore, the poor degree of matching between the soil water supply and plant water demand is the primary reason for the low SWA in dry years.In wet years, the precipitation was concentrated mostly in summer.Therefore, the available soil water distribution was extremely uneven, which in turn resulted in significant water loss in the rainy season, leaving plants vulnerable to drought in the sowing season.

Spatial Perspective
The spatial distribution of the comprehensive SWA was mapped based on the simulation results.As the first step, the original values of the second-class indices were quantitatively calculated on the basis of Equation (3) to Equation (13).In this study, three different typical hydrological-level years were selected (e.g., a wet year, normal year, and dry year) through the frequency analysis of precipitation.The years 2006, 2010, and 2012 were determined to be typical dry, normal, and wet years.As the second step, the original values of the second-class indices were transformed to the standard values as shown in Tables 7-9 according to Equations ( 14) and (15).As the final step, the comprehensive SWA in different administrative regions was calculated on the basis of Equation ( 16).As shown in Tables 7-9, because the index C11 represents the valid capacity of the soil reservoir, which is largely determined by the soil properties in the root zone, the standard values of the index in the three different hydrological years are the same.Spatial distribution maps of SWA in normal, wet, and dry years, based on the standard values of the second-class indexes, are shown in Figures 7-9, respectively.
According to Tables 7-9 and Figures 7-9, the average SWA of Handan City is 0.712, 0.649, and 0.619 in the hydrological normal, wet, and dry years, respectively.In general, when the SWA is below 0.7, the farmland soil water use efficiency (SWUE) is too low to be acceptable; when the SWA is above 0.8, the SWUE can be regarded as satisfactory.With this guidance, the average SWA of Handan City in the wet and dry years is not acceptable.Although the SWA in the normal year is higher than in the wet and dry years, there is still great potential for improvement.From a spatial distribution perspective, the SWA in the west and north mountainous area is lower than that in the east and south plain areas.This may be attributed to the following reasons.(1) The total precipitation in the south and east areas is higher than that in the north and west mountainous areas, and the annual temporal distribution of precipitation in the south and east areas is relatively even, which will increase the degree of matching between the soil water supply and crop water demand and will therefore improve the SWA in the south and east areas; (2) The root zone depth in the west and north mountainous area is thin.Therefore, the factor of the soil water storage capacity is smaller than that of the east and south plain areas; (3) Because the quality of arable land in the eastern and southern plain areas is better than that in the mountainous areas, multiple agricultural management measures such as irrigation and cultivation are implemented in the plain areas, which will improve the use efficiency of the agricultural water resources in the farmland.As shown in Tables 7-9, because the index C11 represents the valid capacity of the soil reservoir, which is largely determined by the soil properties in the root zone, the standard values of the index in the three different hydrological years are the same.Spatial distribution maps of SWA in normal, wet, and dry years, based on the standard values of the second-class indexes, are shown in Figures 7-9, respectively.
According to Tables 7-9 and Figures 7-9, the average SWA of Handan City is 0.712, 0.649, and 0.619 in the hydrological normal, wet, and dry years, respectively.In general, when the SWA is below 0.7, the farmland soil water use efficiency (SWUE) is too low to be acceptable; when the SWA is above 0.8, the SWUE can be regarded as satisfactory.With this guidance, the average SWA of Handan City in the wet and dry years is not acceptable.Although the SWA in the normal year is higher than in the wet and dry years, there is still great potential for improvement.From a spatial distribution perspective, the SWA in the west and north mountainous area is lower than that in the east and south plain areas.This may be attributed to the following reasons.(1) The total precipitation in the south and east areas is higher than that in the north and west mountainous areas, and the annual temporal distribution of precipitation in the south and east areas is relatively even, which will increase the degree of matching between the soil water supply and crop water demand and will therefore improve the SWA in the south and east areas; (2) The root zone depth in the west and north mountainous area is thin.Therefore, the factor of the soil water storage capacity is smaller than that of the east and south plain areas; (3) Because the quality of arable land in the eastern and southern plain areas is better than that in the mountainous areas, multiple agricultural management measures such as irrigation and cultivation are implemented in the plain areas, which will improve the use efficiency of the agricultural water resources in the farmland.Overall, the maps show that the variation in the SWA in different regions of the study area is large and that there is still great potential for improvement in the SWA and SWUE.The SWA distribution maps provide a visual display of its spatial and temporal distribution, allowing easy detection of the problematic regions for targeting specific measures for soil water regulation.Overall, the maps show that the variation in the SWA in different regions of the study area is large and that there is still great potential for improvement in the SWA and SWUE.The SWA distribution maps provide a visual display of its spatial and temporal distribution, allowing easy detection of the problematic regions for targeting specific measures for soil water regulation.Overall, the maps show that the variation in the SWA in different regions of the study area is large and that there is still great potential for improvement in the SWA and SWUE.The SWA distribution maps provide a visual display of its spatial and temporal distribution, allowing easy detection of the problematic regions for targeting specific measures for soil water regulation.

Analysis Main Factors Influencing the Soil Water Availability for Different Regions
As shown in Figures 7-9 1, the factor of C11 represents the valid physical capacity of the soil reservoir, C21 represents the average value of the annual empty storage capacity of the soil reservoir, C22 is the factor for the average annual variation of the empty storage capacity of the soil reservoir, C31 represents the factor for the temporal matching degree between the soil water supply and crop water demand, C32 represents the factor for the spatial matching degree between the soil water supply and crop water demand, C41 is the factor for the soil water conversion efficiency, and C42 is the factor for the crop use efficiency of the soil water.1, the factor of C11 represents the valid physical capacity of the soil reservoir, C21 represents the average value of the annual empty storage capacity of the soil reservoir, C22 is the factor for the average annual variation of the empty storage capacity of the soil reservoir, C31 represents the factor for the temporal matching degree between the soil water supply and crop water demand, C32 represents the factor for the spatial matching degree between the soil water supply and crop water demand, C41 is the factor for the soil water conversion efficiency, and C42 is the factor for the crop use efficiency of the soil water.1, the factor of C11 represents the valid physical capacity of the soil reservoir, C21 represents the average value of the annual empty storage capacity of the soil reservoir, C22 is the factor for the average annual variation of the empty storage capacity of the soil reservoir, C31 represents the factor for the temporal matching degree between the soil water supply and crop water demand, C32 represents the factor for the spatial matching degree between the soil water supply and crop water demand, C41 is the factor for the soil water conversion efficiency, and C42 is the factor for the crop use efficiency of the soil water.
According to the radar maps, the contribution rate of different factors to the SWA in different regions can be clearly illustrated.Among the seven main factors, C42 is the most important with a contribution rate above 0.3 except for Wu'an, Cixian, and urban regions.The factors C21, C31, and C41 are also important with contribution rates between 0.1 and 0.2.Comparatively, C11, C22, and C32 contribute less to the total SWA with contribution rates around or below 0.1.According to the above analysis, to improve the total SWA, we should focus on the four factors (C42, C21, C31, and C41) and take measures to target increases in the contribution rates of these four factors.To analyse the potential for improvement in these contribution rates, we defined an index K to indicate the difference between the current value of the factor and the desired value of the factor, as shown in Equation ( 17): where is the potential of the factor to be improved to a higher contribution rate to SWA, and is the maximum value of the contribution rate of the factor , which is determined on the basis of Equation ( 16).The coefficient of each factor in Equation ( 16) can be regarded as the maximum As shown in Table 1, the factor of C11 represents the valid physical capacity of the soil reservoir, C21 represents the average value of the annual empty storage capacity of the soil reservoir, C22 is the factor for the average annual variation of the empty storage capacity of the soil reservoir, C31 represents the factor for the temporal matching degree between the soil water supply and crop water demand, C32 represents the factor for the spatial matching degree between the soil water supply and crop water demand, C41 is the factor for the soil water conversion efficiency, and C42 is the factor for the crop use efficiency of the soil water.
According to the radar maps, the contribution rate of different factors to the SWA in different regions can be clearly illustrated.Among the seven main factors, C42 is the most important with a contribution rate above 0.3 except for Wu'an, Cixian, and urban regions.The factors C21, C31, and C41 are also important with contribution rates between 0.1 and 0.2.Comparatively, C11, C22, and C32 contribute less to the total SWA with contribution rates around or below 0.1.According to the above analysis, to improve the total SWA, we should focus on the four factors (C42, C21, C31, and C41) and take measures to target increases in the contribution rates of these four factors.To analyse the potential for improvement in these contribution rates, we defined an index K to indicate the difference between the current value of the factor and the desired value of the factor, as shown in Equation ( 17): where K n is the potential of the factor n to be improved to a higher contribution rate to SWA, and V d is the maximum value of the contribution rate of the factor n, which is determined on the basis of Equation ( 16).The each factor in Equation ( 16) can be regarded as the maximum value of the contribution rate of the corresponding factor; V c is the current value of the contribution rate of the factor n.
The potential of each factor in normal hydrological years for different regions was calculated on the basis of Equation ( 17), as shown in Table 10.According to the analysis in Table 10, it is concluded that the potential for improvement of the factors differs among regions.For example, in the Wu'an Region, the indices K C11 and K C42 were much higher those of other factors, which means that the factors C11 and C42 have a large potential to be increased for a higher contribution to SWA.However, in the Jize Region, the indices K C32 and K C42 were much higher than those of other factors, which means that the factors C32 and C42 have large potential to be increased for a higher contribution to SWA.The values of the index K changed dramatically among the regions.Therefore, to improve the SWA in the study area, different measures must be taken according to the different potential of the factors in different regions.

Strategies to Improve Soil Water Availability
The concept of SWA is a new comprehensive indicator to quantitatively characterise the SWUE of farmland at a regional scale.This indicator not only considers the traditional index of field water use efficiency but also considers the soil water storage capacity of the soil layer, the actual capacity of the soil reservoir, and the spatial and temporal matching degree between the soil water supply and crop water demand.Therefore, SWA is a better index for SWUE evaluation and soil water management.
According to the analysis in Table 10, to improve the SWA, it is necessary to take different measures in different regions because the potential of the factors to be improved for a higher contribution to SWA is different.
First, in the Wu'an, Shexian, Cixian, Handan, and Fengqu regions, the factor C11 has a large potential for improvement because the index K is above 40% in the five regions and is as high as 90% in the Shexian Region.Therefore, using agronomic measures such as deep ploughing and straw returning to increase the valid capacity of the soil layer in the five regions is a key step for improving the SWA.
Second, in the Jize, Quzhou, Handan, Feixiang, Guantao, Guangping, Cheng'an, Weixian, and Linzhang regions, the factor C32 has a large potential for improvement because the index K is above 40% in the nine regions.This indicates that the degree of spatial matching between the soil water supply and crop water very poor in the nine regions.Therefore, the planting structure in these regions should be adjusted to reduce the area of water-consuming crops such as winter wheat or to enlarge the area of water-efficient crops such as summer maize and potatoes.
Finally, as shown in Table 10, K C42 was mostly above 30% in all 17 regions, although the index was only 29.6% in the Linzhang Region.Obviously, the potential of improving the C42 factor is quite large in all of the regions in the entire study area.Therefore, increasing the soil water use efficiency for crop growth, such as mulching to reduce invalid evaporation and adopting regulated deficit irrigation, is a key approach for increasing the SWA.

Conclusions
Under the joint effects of climate change and human activities, the availability of blue water resources is severely reduced for agricultural production, particularly in arid and semi-arid areas.Therefore, soil water, known as green water resources, is becoming more and more important for sustainable development in agriculture and ecology.In this context, to improve the SWUE, this study put forward a new method to quantitatively evaluate the SWA based on a self-developed hydrological model at the regional scale and presented related strategies according to the main factors leading to low efficiency of soil water utilization.In summary, the main innovations in this research include two aspects.Firstly, this study originally proposed the concept of SWA and established a new evaluation index system to quantitatively calculate the SWA through the analytic hierarchy process method.Secondly, a self-developed distributed hydrological model (MODCYCLE) was introduced to simulate the soil moisture change process for SWA evaluation.However, modelling farmland soil moisture through a hydrological model is still challenging because of the high variability of underlying surface conditions and intensive human management measures, particularly in the Chinese farming system.Therefore, this study attempted the novel development of MODCYCLE considering the complicated tilling practices and management schemes in the North China Plain.Different from the SWAT model, MODCYCLE includes a special agricultural management module to depict the complex hydrological processes of Chinese farmland.Therefore, the simulation accuracy is greatly improved by using MODECYCLE to evaluate the SWA in the study area with complicated tilling practices and strong anthropogenic influences.
To verify the SWA evaluation method and the MODCYCLE model, Handan City in North China Plain was selected as a case study.The results show that the MODCYCLE model is an efficient tool for simulating the farmland hydrological process in the study area.According to the comparison between the observed soil moisture and simulated soil moisture results, the relative errors are all below 10%, and the correlation coefficients are almost above 0.7 except for Long-wangmiao and Linzhang stations.The SWA evaluation results indicate that the average SWA in the entire area ranges from 0.62 to 0.71 in different hydrological level years, which demonstrates that the potential is large for SWA improvement for a higher SWUE.However, the soil water availability changes greatly in different regions because the main factors influencing the SWA are different.Based on these findings, we suggest that managers can use the SWA index system to identify the main factors causing low soil water utilization efficiency in different regions and to select and determine the strategies for improving the SWA.Improving the reliability of the hydrological model is the foundation for SWA evaluation.However, owing to the complicated tillage practices and irrigation schemes in China, the farmland hydrological characteristics vary extensively among different regions, which poses a great challenge for model calibration and validation.Therefore, more attention should be paid to model improvement and validation when the model is adapted to other regions in future studies.

Figure 1 .
Figure 1.Hydrological processes in the 1D soil column method.

Figure 1 .
Figure 1.Hydrological processes in the 1D soil column method.

25 Figure 2 .
Figure 2. Location map of the study area.

Figure 2 .
Figure 2. Location map of the study area.

Figure 4 .
Figure 4. Comparison of observed values and simulated values of the 12 monitoring stations.Figure 4. Comparison of observed values and simulated values of the 12 monitoring stations.

Figure 4 .
Figure 4. Comparison of observed values and simulated values of the 12 monitoring stations.Figure 4. Comparison of observed values and simulated values of the 12 monitoring stations.

3. 3 . 1 .
Temporal PerspectiveAs previously mentioned, the simulation period in this study was from 2001 to 2014.Based on the simulation results, the comprehensive soil water availability in the entire study area from 2001 to 2014 was calculated, as shown in Figure5.

Figure 5 .
Figure 5. Comprehensive soil water availability in the simulation period.

Figure 5
Figure5shows that the comprehensive SWA ranged from 59% to 76.8% in different years.The maximum value of 76.8% occurred in 2014 with yearly precipitation of 507 mm, and the minimum value of 59% occurred in 2002 with yearly precipitation of 346 mm.To analyse the relationship

Figure 5 .
Figure 5. Comprehensive soil water availability in the simulation period.

Figure 6 .
Figure 6.SWA ranking in different hydrological level years.

Figure 6 .
Figure 6.SWA ranking in different hydrological level years.

Figure 7 .
Figure 7. Spatial distribution of SWA in the normal hydrological year (2010).

Figure 7 .
Figure 7. Spatial distribution of SWA in the normal hydrological year (2010).

Figure 9 .
Figure 9. Spatial distribution of SWA in the hydrological dry year (2006).

Figure 8 .
Figure 8. Spatial distribution of SWA in the hydrological wet year (2012).

Figure 9 .
Figure 9. Spatial distribution of SWA in the hydrological dry year (2006).

Figure 9 .
Figure 9. Spatial distribution of SWA in the hydrological dry year (2006).

1 .
, the SWA varies greatly in different regions because the values of the influencing factors (C 11 , C 21 , C 22 , C 31 , C 32 , C 41 , and C 42 ) differ among regions.To analyze the main factors influencing the SWA for each region, we separately calculated the contribution ratio of each factor to the SWA in each region, and we plotted the statistical radar maps as shown in Figures 10 and 11.Analysis of the Main Factors Influencing the Soil Water Availability for Different Regions As shown in Figures 7-9, the SWA varies greatly in different regions because the values of the influencing factors (C , C , C , C , C , C , and C ) differ among regions.To analyze the main factors influencing the SWA for each region, we separately calculated the contribution ratio of each factor to the SWA in each region, and we plotted the statistical radar maps as shown in Figures 10 and 11 .

Figure 10 .
Figure 10.Contribution ratio of different factors to the comprehensive SWA in Wu'an, Jize, Qiuxian, Yongnian, Quzhou, Handan, Feixiang, Guantao, and Shexian regions in different hydrological-level years.As shown in Table1, the factor of C11 represents the valid physical capacity of the soil reservoir, C21 represents the average value of the annual empty storage capacity of the soil reservoir, C22 is the factor for the average annual variation of the empty storage capacity of the soil reservoir, C31 represents the factor for the temporal matching degree between the soil water supply and crop water demand, C32 represents the factor for the spatial matching degree between the soil water supply and crop water demand, C41 is the factor for the soil water conversion efficiency, and C42 is the factor for the crop use efficiency of the soil water.

Figure 10 .
Figure 10.Contribution ratio of different factors to the comprehensive SWA in Wu'an, Jize, Qiuxian, Yongnian, Quzhou, Handan, Feixiang, Guantao, and Shexian regions in different hydrological-level years.As shown in Table1, the factor of C11 represents the valid physical capacity of the soil reservoir, C21 represents the average value of the annual empty storage capacity of the soil reservoir, C22 is the factor for the average annual variation of the empty storage capacity of the soil reservoir, C31 represents the factor for the temporal matching degree between the soil water supply and crop water demand, C32 represents the factor for the spatial matching degree between the soil water supply and crop water demand, C41 is the factor for the soil water conversion efficiency, and C42 is the factor for the crop use efficiency of the soil water.

Figure 11 .
Figure 11.Contribution ratio of different factors to the comprehensive SWA in Guangping, Cheng'an, Weixian, Cixian, Linzhang, Daming, Fengqu, and urban regions in different hydrological-level years.As shown in Table1, the factor of C11 represents the valid physical capacity of the soil reservoir, C21 represents the average value of the annual empty storage capacity of the soil reservoir, C22 is the factor for the average annual variation of the empty storage capacity of the soil reservoir, C31 represents the factor for the temporal matching degree between the soil water supply and crop water demand, C32 represents the factor for the spatial matching degree between the soil water supply and crop water demand, C41 is the factor for the soil water conversion efficiency, and C42 is the factor for the crop use efficiency of the soil water.

Figure 11 .
Figure 11.Contribution ratio of different factors to the comprehensive SWA in Guangping, Cheng'an, Weixian, Cixian, Linzhang, Daming, Fengqu, and urban regions in different hydrological-level years.As shown in Table1, the factor of C11 represents the valid physical capacity of the soil reservoir, C21 represents the average value of the annual empty storage capacity of the soil reservoir, C22 is the factor for the average annual variation of the empty storage capacity of the soil reservoir, C31 represents the factor for the temporal matching degree between the soil water supply and crop water demand, C32 represents the factor for the spatial matching degree between the soil water supply and crop water demand, C41 is the factor for the soil water conversion efficiency, and C42 is the factor for the crop use efficiency of the soil water.

Table 1 .
Evaluation index system for farmland soil water availability.

Table 4 .
Key parameters and their optimal values used in MODCYCLE.

Table 5 .
Soil moisture monitoring stations in the study area.

Table 5 .
Soil moisture monitoring stations in the study area.

Table 6 .
Statistical results showing the consistency between observed and simulated values.

Table 6 .
Statistical results showing the consistency between observed and simulated values.

Table 7 .
Standard values of the second-class indices and the total SWA of different regions in the normal year (2010).

Table 8 .
Standard values of the second-class indices and the total SWA of different regions in the wet year (2012).

Table 9 .
Standard values of the second-class indices and the total SWA of different regions in the dry year (2006).

Table 10 .
Potentials of the factors to be improved for a higher contribution to SWA in normal hydrological years.