Simulation of the Water Dynamics and Root Water Uptake of Winter Wheat in Irrigation at Different Soil Depths

Xianghong Guo 1,* ID , Xihuan Sun 1,2, Juanjuan Ma 1, Tao Lei 1, Lijian Zheng 1 and Pu Wang 1 1 College of Water Resource Science and Engineering, Taiyuan University of Technology, Taiyuan 030024, China; sunxihuan@tyut.edu.cn (X.S.); mjjsxty@163.com (J.M.); lcsyt@126.com (T.L.); 13623661562@163.com (L.Z.); wangpusxty@163.com (P.W.) 2 Jinzhong University, Jinzhong 030600, China * Correspondence: xianghong7920@126.com; Tel.: +86-0351-6010102


Introduction
The shortage of water resources seriously restricts the sustainable development of agriculture in China.The contradiction between the supply and demand of agricultural water is particularly prominent in northern China.Winter wheat is one of the main grain crops in China.It is mainly planted in northern China.The planting area and yield of winter wheat account for more than half of the total planting area and yield of winter wheat in China [1,2].However, the annual rainfall in this area is low, and the rainfall distribution is uneven.Approximately 70% of the rainfall is concentrated in the June-September period, which does not coincide with the winter wheat growth period (i.e., October-May).The water consumption of 2/3-4/5 in the entire growth period of winter wheat needs to be satisfied through irrigation [3].Surface irrigation is still the main irrigation method for winter wheat in China because of its simple operation.However, this method requires a large amount of irrigation water because of a high soil water evaporation, which accounts for 25-30% of the total winter wheat water consumption [4].The water use efficiency of the irrigation method is low, and it also leads to soil hardening and structural damage [4,5].Consequently, finding reasonable irrigation methods for winter wheat production in northern China, maximising the use of deep soil water, reducing the evaporation and irrigation water consumption, and improving water use efficiency are urgent problems that need to be addressed.
Different irrigation methods result in different water distributions in the soil profile and different impacts on the root growth distribution [6][7][8].Wang et al. [9] showed that with the same irrigation lower limit and irrigation quota, the irrigation method had a significant effect on the spatial distribution of winter wheat roots.They also demonstrated that surface drip irrigation and basin irrigation significantly promoted the distribution of roots in 0-50 cm soil, whereas the root distribution was relatively uniform in 0-100 cm soil under subsurface drip irrigation [9].Lv et al. [7,10] found that the main distribution layer of winter wheat roots under drip irrigation also moved upwards significantly compared with that under surface irrigation.The root length density (RLD) in deep soil was relatively low and the root water uptake (RWU) was mainly concentrated in the upper soil [7,10].Zhang et al. [11] compared the distributions of cotton roots between border irrigation and surface drip irrigation and found that surface drip irrigation was beneficial to the root growth in the shallow root zone (0-30 cm), whereas border irrigation was beneficial to the root growth in the deeper root zone (below 30 cm).Studies have shown that the root distribution is closely related to the soil water content (SWC), and an upper drought is conducive to the development of the root system in the deep layer, as well as the absorption and utilisation of deep soil water [12][13][14][15][16].
Deep irrigation is an irrigation method which can increase the SWC in the deep layer of the root area by artificially controlling the irrigation soil depth.This method can promote root growth in deep soil, increase the root uptake of deep soil water, and improve the efficiency of soil water use [17,18].Wang et al. [19], Di et al. [20], and Huang et al. [21] conducted experiments on the growth and water regulation of winter wheat in irrigation at different soil depths by using a soil column method.The results showed that the surface SWC under deep irrigation is lower than that under surface irrigation.Furthermore, a greater irrigation depth results in a greater soil moisture content in the deep soil layer, a greater depth of rooting of winter wheat roots, a greater total root length, and a greater root weight in deep layers; furthermore, moderate deep irrigation can increase wheat yield and water use efficiency [19,20].Studies have shown that the distribution pattern of RWU is an important indicator for explaining the increase in the yield of winter wheat.However, the distribution patterns and dynamic changes of RWU of winter wheat under deep irrigation are still unclear.In recent years, an increasing number of scholars have adopted the soil water movement model under the condition of RWU to study the soil water dynamics and RWU characteristics of crops [22][23][24][25].Therefore, establishing the model for soil water movement and RWU in irrigation at different soil depths and quantifying the effects of irrigation soil depth on SWC and RWU of winter wheat are of great significance for revealing the mechanism of water saving and yield increase of winter wheat under deep irrigation.
The main objectives of this study were (1) to analyse the effects of irrigation at different soil depths on the root distribution of winter wheat, (2) to establish the model for the water dynamics and RWU of winter wheat during irrigation at different soil depths and (3) to reveal the RWU rate distribution of winter wheat in irrigation at different soil depths.

Experimental Site
The experiment was conducted at the training base of Shanxi Conservancy Technical Institute from September 2016 to June 2017.The training base is located in Yanhu District, Yuncheng City, Shanxi Province (latitude 34 • 48 27 N, longitude 110 • 41 23 E).The site has a typical continental semi-arid climate under the influence of the warm-temperate monsoon.The average annual rainfall is 559.3 mm, which is mainly concentrated in the July-September period.The average annual temperature is 13.6 • C. The frost-free period is approximately 220 days.The average annual sunshine hours are 2247.4h.The soil in the test area is silty clay loam.The mechanical composition and physical properties of the soil are shown in Table 1.

Experimental Design
The purpose of the experiment was to find out the effect of deep irrigation on the growth of winter wheat.Three treatment experiments were performed in irrigation at different soil depths, namely, T1 for a soil depth of irrigation of 0% of the root distribution depth (surface irrigation), T2 for 40% of the root distribution depth, and T3 for 70% of the root distribution depth.The winter wheat cultivar was Liangxing99 (Shandong LiangXing seeds Co., Ltd, Dezhou City, Shandong Province, China), a mid-late-maturing variety [18] sown on 12 October 2016 and harvested on 26 May 2017.The experiments were conducted with 48 buried PVC (Polyvinyl Chloride) pipe-soil columns (outer diameter 20 cm, internal diameter 18.6 cm, length 3 m) [21], with each treatment consisting of 16 soil columns.The test was carried out as follows to ensure that the soil column represents the actual conditions in the field.Firstly, 48 holes with a diameter of 20 cm and a depth of 3 m were dug in the field.The excavated soil was stacked in layers.The SWC was determined.Secondly, the weight of the soil to be filled up to 5 cm in height for every soil column was calculated on the basis of the soil bulk density and moisture content of each layer, the excavated soil was filled into the soil column layer by layer according to the calculated weight, and then the 48 soil columns were placed in the design position.Finally, the wheat was sown into the soil column.The experimental field layout is shown in Figure 1.A movable shelter was installed on the experimental plot to avoid the influence of rainfall on the experiment.

Experimental Design
The purpose of the experiment was to find out the effect of deep irrigation on the growth of winter wheat.Three treatment experiments were performed in irrigation at different soil depths, namely, T1 for a soil depth of irrigation of 0% of the root distribution depth (surface irrigation), T2 for 40% of the root distribution depth, and T3 for 70% of the root distribution depth.The winter wheat cultivar was Liangxing99 (Shandong LiangXing seeds Co., Ltd, Dezhou City, Shandong Province, China), a mid-late-maturing variety [18] sown on 12 October 2016 and harvested on 26 May 2017.The experiments were conducted with 48 buried PVC (Polyvinyl Chloride) pipe-soil columns (outer diameter 20 cm, internal diameter 18.6 cm, length 3 m) [21], with each treatment consisting of 16 soil columns.The test was carried out as follows to ensure that the soil column represents the actual conditions in the field.Firstly, 48 holes with a diameter of 20 cm and a depth of 3 m were dug in the field.The excavated soil was stacked in layers.The SWC was determined.Secondly, the weight of the soil to be filled up to 5 cm in height for every soil column was calculated on the basis of the soil bulk density and moisture content of each layer, the excavated soil was filled into the soil column layer by layer according to the calculated weight, and then the 48 soil columns were placed in the design position.Finally, the wheat was sown into the soil column.The experimental field layout is shown in Figure 1.A movable shelter was installed on the experimental plot to avoid the influence of rainfall on the experiment.

Irrigation Design
Irrigation was conducted at five different stages during the experiment, as follows: overwintering stage (20 December), returning green stage (8 March), jointing stage (4 April), heading stage (27 April), and filling stage (10 May).Surface irrigation was conducted at the overwintering stage, whereas irrigation was conducted according to the design irrigation depth at the four succeeding stages.Therefore, to determine the soil depth of irrigation at each irrigation stage, three soil columns were taken out for each treatment before irrigation, the soil column was cut from the middle, and the root distribution depth of winter wheat was measured.The soil depth of irrigation of each treatment is obtained by the following: where H Di is the depth of irrigation (cm); α i is the irrigation depth coefficient, i = 1, 2, 3 (the irrigation depth coefficients of T1, T2, and T3 were 0, 0.4, and 0.7 respectively).H mi is the root distribution depth of winter wheat (cm).
In this study, the effects of irrigation quota and irrigation time were not considered, and both were in accordance with the irrigation time and irrigation quota of the local winter wheat planting habits, that is, the irrigation quota was 67.5 mm.During irrigation, the holes were symmetrically drilled at regular intervals on both sides of the PVC pipe wall, and the water supply bottles and the holes in the pipe wall were connected by drop pipes with a 3 mm inner diameter to supply water to the deep soil layers.The irrigation amount of the deep irrigation hole (below 30 cm) is calculated by Equation ( 2).The residual water, which is obtained by subtracting the deep irrigation volume from the total irrigation amount, is irrigated from the surface.The irrigation time, irrigation location and irrigation quota of each treatment are shown in Table 2.
where θ i and θ 0i denote the initial soil moisture content (cm 3 /cm 3 ) and irrigation upper limit (cm 3 /cm 3 , 85% of field water capacity), respectively, and H is the thickness of the planned wetting layer of the irrigation hole (cm).

Measurement Methods
SWC was determined by the time-domain reflectometry (TRIME-PICO IPH, IMKO, Ettlingen, Germany), at soil depths from 0 m to 3 m at intervals of 0.2 m every 5-7 days.Before each irrigation, the soil column was opened, the roots of winter wheat were removed from the soil column every 10 cm layer, the depth of root distribution was obtained, and the RLD was measured by WinRHIZO software (Regent Instruments Inc., Québec, QC, Canada) [20].Every 2 weeks, the length and width of all the leaves of winter wheat in the soil column were measured with a ruler, and the leaf area and leaf area index were calculated [21].At maturity, the winter wheat plants in all soil columns were harvested and the grains were dried to a constant weight at 75 • C. Subsequently, the yield of each treatment was obtained by weighing.The meteorological data, including rainfall, temperature, and humidity were collected every 10 min by the Wireless Automatic Weather Monitoring Station (ADCON Telemetry, Klosterneuburg, Austria) located in the field.

Governing Equation
The field experiment showed that the soil moisture movement can be simulated with a one-dimensional soil water movement equation in irrigation at different soil depths and the key task was to simulate irrigation at different depths.In the experiment, irrigation at different depths was achieved by quantitatively supplying water to the irrigation holes at different depths of the soil column through the drip tube.For the irrigation simulation, the source item is added to the soil water movement equation at different depths of irrigation.On the assumption that the soil is homogeneous and isotropic in each soil layer, the governing equation for the water flow is the following 1D Richards' equation: where θ is the volumetric water content (cm 3 /cm 3 ); h is the soil water pressure head (cm); t is the time (hour); z is the vertical space coordinate (cm); K (hour) is the unsaturated hydraulic conductivity (cm/h); and Q is the intensity of water supply at different depths (1/hour), which is meaningful only at the irrigation depth of the irrigation period but 0 at other time and depths; and S is the RWU rate (l/hour).
where h 50 represents the pressure head at which the water extraction rate is reduced by 50% under negligible osmotic stress, and p is an experimental constant with the value of approximately 3.
where LAI is leaf area index; k and A are constants with values of 0.3973 and 0.0136 [29], respectively; ET c is the potential crop evapotranspiration, which was estimated as follows: where ET 0 is the reference crop evapotranspiration calculated using the Penman-Monteith formula [30], and K c is the crop coefficient, which is determined according to a previous study [30].

Initial Condition
The initial condition is expressed as follows: where h st (z) is the suction head that corresponds to the initial water content (cm).

Boundary Conditions
The top boundary conditions are described as follows: where E s is the surface evaporation intensity (cm/hour).
where h 300 is the pressure head that corresponds to the measured water content at the bottom (cm).
In a previous study [31], E s is calculated using the following: where h c and h cc are the critical values of water potential at the soil surface (cm), and E p is the potential evaporation on the soil surface (cm/hour) calculated by the following:

Model Parameter Calibration
The model includes eight unknown parameters: θ r , θ s , α, n, and K s in the soil water movement parameter VG model; h c and h cc in the soil evaporation model and h 50 in the water stress function.The rationality of the model parameters directly affects the simulation accuracy of the model.In this study, the inverse method was used to solve the model parameters.The inverse method was used to determine the model parameters by solving the minimum value of the objective function that represents the difference between the experimental value and the model predicted value: where θ j is the measured SWC, θ(h j , X) is the calculated SWC of the model, M is the number of measured soil moisture samples, and X is the parameter vector (θ r , θ s , α, n, K s , h c , h cc , h 50 ) to be solved.
In this study, the measured data of the three treatments were divided into two parts.The data of T1 were used for the calibration of the model parameters, and those of the other treatments (i.e., T2 and T3) were used for the model verification.The measured SWC data of T1 was substituted into Equation (15), and Equation ( 15) was solved using the parallel multi-start optimisation algorithm in the MATLAB Optimization Toolbox.The model parameters obtained are shown in Table 3.

Model Evaluation
The simulation precision of the model was evaluated using the correlation coefficient (R), root-mean-square error (RMSE), mean absolute error (MAE), and mean relative error (MRE), which were calculated as follows: where θ s is the simulated water content (cm 3 /cm 3 ), θ R is the measured moisture content (cm 3 /cm 3 ), θ is the average water content (cm 3 /cm 3 ), and l is the number of measuring points.

Spatial and Temporal Distributions of Winter Wheat at Different Irrigation Depths
Figure 2 shows the vertical distributions of the RLD of winter wheat during irrigation at different depths in different growth periods.As shown in the figure, the RLD under different treatments initially increased with the growth period, reached the maximum in the filling period, and then gradually decayed.From the jointing stage to the filling stage, the mean RLDs of T1, T2, and T3 increased by 37%, 26%, and 42% respectively, whereas from the filling stage to the maturity stage, their mean RLDs decreased by 15%, 13%, and 19% respectively.Therefore, the RLDs of winter wheat at different growth stages significantly differed (p < 0.05).The variation in the root system of winter wheat at different growth stages is consistent with the results reported by Zhang et al. [2].The RLDs of the three treatments decreased with an increase in the depth, but the effects of different irrigation depths on the RLD differed.At the heading stage, the mean RLD decreased with an increase in the irrigation depth above 50 cm, and the mean RLDs of T2 and T3 decreased by 16% and 26% compared with that of T1.At the same stage, the mean RLD increased with an increase in the irrigation soil depth below 50 cm, and the mean RLDs of T2 and T3 increased by 31% and 125% compared with that of the T1 treatment.At the three other growth stages, the mean RLD decreased with an increase in the irrigation soil depth above 30 cm; compared with the RLDs of T1, the mean RLDs of T2 and T3 decreased by 6% and 36% at the heading stage, 18% and 29% at the filling stage, and 15% and 30% at the maturity stage.The mean RLD increased with an increase in the depth of irrigation below 30 cm depth; compared with the RLDs of T1, the mean RLDs of T2 and T3 increased by 41% and 132% at the heading stage, 26% and 94% at the filling stage, and 19% and 59% at the maturity stage.Therefore, the soil depth of irrigation has a significant effect on the RLD of winter wheat (p < 0.01), and deep irrigation can increase the growth of deep soil roots and reduce the growth of surface soil roots.The spatial and temporal distributions of the root system in the soil profile are affected not only by the crop species and variety but also by environmental factors, especially soil moisture [32].Studies have shown that when the SWC is lower in the upper soil layer and higher in the lower soil layer, the root system penetrates into the deeper layers of soil in accordance with the hydrotropism of roots [24].In this study, the irrigation water of T2 and that of T3 were directly The RLDs of the three treatments decreased with an increase in the depth, but the effects of different irrigation depths on the RLD differed.At the heading stage, the mean RLD decreased with an increase in the irrigation depth above 50 cm, and the mean RLDs of T2 and T3 decreased by 16% and 26% compared with that of T1.At the same stage, the mean RLD increased with an increase in the irrigation soil depth below 50 cm, and the mean RLDs of T2 and T3 increased by 31% and 125% compared with that of the T1 treatment.At the three other growth stages, the mean RLD decreased with an increase in the irrigation soil depth above 30 cm; compared with the RLDs of T1, the mean RLDs of T2 and T3 decreased by 6% and 36% at the heading stage, 18% and 29% at the filling stage, and 15% and 30% at the maturity stage.The mean RLD increased with an increase in the depth of irrigation below 30 cm depth; compared with the RLDs of T1, the mean RLDs of T2 and T3 increased by 41% and 132% at the heading stage, 26% and 94% at the filling stage, and 19% and 59% at the maturity stage.Therefore, the soil depth of irrigation has a significant effect on the RLD of winter wheat (p < 0.01), and deep irrigation can increase the growth of deep soil roots and reduce the growth of surface soil roots.The spatial and temporal distributions of the root system in the soil profile are affected not only by the crop species and variety but also by environmental factors, especially soil moisture [32].Studies have shown that when the SWC is lower in the upper soil layer and higher in the lower soil layer, the root system penetrates into the deeper layers of soil in accordance with the hydrotropism of roots [24].In this study, the irrigation water of T2 and that of T3 were directly injected into the soil from the deeper layer, and the SWCs of T2 and T3 were higher than that of T1, so the RLDs of T2 and T3 in the deeper soil were higher than that of T1.
In order to study the spatial and temporal distributions of the RLD of winter wheat quantitatively, the spatial-temporal distribution of the RLD of winter wheat was fitted using Equation (20).The results of the fitting are shown in Table 3.
where L(z, t) is the RLD (cm/cm 3 ); z is the vertical space coordinate (cm); t is the number of days after winter wheat sowing(day); and a, b, c, d, and f are the fitting parameters.Table 4 shows that the correlation coefficients of the Equation ( 20) fitting under different treatments were above 0.86, indicating that Equation ( 20) can be applied to the spatial and temporal distributions of the RLD of winter wheat in irrigation at different soil depths.The parameter a represents the maximum RLD.Table 4 shows that the value of a decreased with an increase in the soil depth of irrigation, that is, the maximum RLD decreased with an increase in the soil depth of irrigation.The parameter b represents the RLD declining rate with depth, the b value decreased with an increase in the soil depth of irrigation, indicating that the RLD of winter wheat decreased with the soil depth as the depth of irrigation increased.Deep irrigation can reduce the RLD declining rate of winter wheat; the greater the depth of irrigation, the lower the deceleration rate.The parameter c indicates the position of the maximum RLD in the vertical direction.The c values under different treatments ranged between 5 and 7, indicating that the maximum RLD in the vertical direction appeared on the surface.The parameter f represents the time at which the maximum RLD occurs in the time direction.The values of f under different treatments were between 192 and 205, showing that the maximum RLD of winter wheat at different irrigation depths appeared at the filling stage.

Evaluation of the Soil Water Movement Model in Irrigation at Different Soil Depths
Figure 3 shows the linear relationship between the measured and predicted values of SWC.The slopes of the linear equations of T1, T2, and T3 were 0.92, 1.02, and 0.93 respectively, and the correlation coefficients (R) reached 0.93, 0.91, and 0.91, respectively.The simulated values were significantly correlated with the experimental values (p < 0.01), indicating a good consistency between the experimental and simulated SWCs.The MAEs, MREs, and RMSEs of the experimental and simulated SWCs are shown in Table 5.Table 5 shows that the MAEs, MREs, and RMSEs of the parameter calibration treatment (T1) and model validation treatments (T2 and T3) were relatively small.The MAEs of the three treatments were not higher than 0.022, their MREs were not higher than 7.95%, and their RMSEs were not higher than 0.028.Therefore, the model of soil water movement in irrigation at different soil depths exhibit a high simulation accuracy and can be used to simulate soil water movement.

Simulation Analysis of Soil Water Dynamics at Different Irrigation Depths
Soil water dynamics are affected by many factors, including soil texture, soil bulk density, rainfall, evaporation, irrigation method, and RWU.Amongst these factors, the irrigation method is one of the most important factors.Figure 4 shows comparisons of the calculated and measured soil water dynamics of winter wheat at different irrigation depths.The simulation period is from the first day after the implementation of different irrigation depth (8 March, 148 days after sowing) to the harvest (26 May, 226 days after sowing).As illustrated in Figure 4, the soil moisture dynamics demonstrated a fluctuation change law with time, and the fluctuation change decreased gradually with increasing depth.This outcome is because the soil moisture change in this study was mainly caused by irrigation, evaporation, and RWU.The SWC increased after irrigation and then gradually decreased over time under the effects of soil evaporation and RWU.Moreover, the greater the irrigation soil depth, the greater the degree of soil moisture increase.The main depth ranges of SWC increase in T1, T2, and T3 were 0-50 cm, 0-110 cm, and 0-190 cm respectively.However, the greater the irrigation soil depth, the lower the SWC in the 0-50 cm soil layers.The average SWCs in T2 and T3 decreased by 9% and 13% compared with that of T1 in 0-50 cm.The calculated values are consistent with the measured values and are in good agreement, as shown in Figure 4. Therefore, the mathematical model for the soil water movement of winter wheat under different irrigation soil depth conditions can simulate the dynamic changes in the soil water content of winter wheat.

Simulation Analysis of Soil Water Dynamics at Different Irrigation Depths
Soil water dynamics are affected by many factors, including soil texture, soil bulk density, rainfall, evaporation, irrigation method, and RWU.Amongst these factors, the irrigation method is one of the most important factors.Figure 4 shows comparisons of the calculated and measured soil water dynamics of winter wheat at different irrigation depths.The simulation period is from the first day after the implementation of different irrigation depth (8 March, 148 days after sowing) to the harvest (26 May, 226 days after sowing).As illustrated in Figure 4, the soil moisture dynamics demonstrated a fluctuation change law with time, and the fluctuation change decreased gradually with increasing depth.This outcome is because the soil moisture change in this study was mainly caused by irrigation, evaporation, and RWU.The SWC increased after irrigation and then gradually decreased over time under the effects of soil evaporation and RWU.Moreover, the greater the irrigation soil depth, the greater the degree of soil moisture increase.The main depth ranges of SWC increase in T1, T2, and T3 were 0-50 cm, 0-110 cm, and 0-190 cm respectively.However, the greater the irrigation soil depth, the lower the SWC in the 0-50 cm soil layers.The average SWCs in T2 and T3 decreased by 9% and 13% compared with that of T1 in 0-50 cm.The calculated values are consistent with the measured values and are in good agreement, as shown in Figure 4. Therefore, the mathematical model for the soil water movement of winter wheat under different irrigation soil depth conditions can simulate the dynamic changes in the soil water content of winter wheat.

Simulation Analysis of RWU of Winter Wheat During Irrigation at Different Soil Depths
Figure 5 illustrates that the RWU of winter wheat decreased sharply with depth and mainly occurred in the upper soil layer.The main range and intensity of RWU of winter wheat initially increased, then subsequently decreased over time, and reached the maximum during heading and filling stages.This phenomenon is because at these stages, the growth of winter wheat is optimal and the transpiration intensity is the highest [2].Moreover, the greater the soil depth of irrigation, the deeper the RWU of winter wheat.The soil depth of irrigation increased from T1 to T3, and the RWU depth ranges also increased from 0-100 cm to 0-250 cm.In the upper soil layer, the SWC and RLD of T1 were greater than those of T2 and T3, so the RWU rate of T1 was greater than those of T2 and T3.Studies have that the distribution pattern of root water absorption in soil profiles depends on the distribution of RLD and the soil moisture content, that is, a greater RLD results in a higher SWC and higher root water absorption rate [24].In the middle and deep soil layers, the SWC and RLD of T3 were greater than those of T1 and T2, so the RWU rate of T3 was greater than those of T1 and T2.This result is similar to the results reported by He et al. [33], although the depth of RWU in their study is shallower.Figure 6a shows that the RWU in each layer decreased with an increase in the soil depth.In the 0-50 cm soil layer, the RWU values of T2 and T3 were lower than that of T1.By contrast, in the 50-300 cm layers, the RWU values of T2 and T3 were higher than that of T1, with that of T3 being much higher than those of T1 and T2 in the 150 cm-300 cm soil layers.

Simulation Analysis of RWU of Winter Wheat During Irrigation at Different Soil Depths
Figure 5 illustrates that the RWU of winter wheat decreased sharply with depth and mainly occurred in the upper soil layer.The main range and intensity of RWU of winter wheat initially increased, then subsequently decreased over time, and reached the maximum during heading and filling stages.This phenomenon is because at these stages, the growth of winter wheat is optimal and the transpiration intensity is the highest [2].Moreover, the greater the soil depth of irrigation, the deeper the RWU of winter wheat.The soil depth of irrigation increased from T1 to T3, and the RWU depth ranges also increased from 0-100 cm to 0-250 cm.In the upper soil layer, the SWC and RLD of T1 were greater than those of T2 and T3, so the RWU rate of T1 was greater than those of T2 and T3.Studies have shown that the distribution pattern of root water absorption in soil profiles depends on the distribution of RLD and the soil moisture content, that is, a greater RLD results in a higher SWC and higher root water absorption rate [24].In the middle and deep soil layers, the SWC and RLD of T3 were greater than those of T1 and T2, so the RWU rate of T3 was greater than those of T1 and T2.This result is similar to the results reported by He et al. [33], although the depth of RWU in their study is shallower.Figure 6a shows that the RWU in each layer decreased with an increase in the soil depth.In the 0-50 cm soil layer, the RWU values of T2 and T3 were lower than that of T1.By contrast, in the 50-300 cm layers, the RWU values of T2 and T3 were higher than that of T1, with that of T3 being much higher than those of T1 and T2 in the 150 cm-300 cm soil layers.Figure 6b shows that from the jointing stage to harvest, the total RWU (Transpiration, Tr) of T2 and T3 in the 0-300 cm layers increased by 0.5% and 3.8% compared with that of T1, respectively, and the total evaporation (E) of T2 and T3 decreased by 20.9% and 63.3% compared with that of T1, Figure 6b shows that from the jointing stage to harvest, the total RWU (Transpiration, Tr) of T2 and T3 in the 0-300 cm layers increased by 0.5% and 3.8% compared with that of T1, respectively, and the total evaporation (E) of T2 and T3 decreased by 20.9% and 63.3% compared with that of T1, respectively, and the evapotranspiration (ET) of T2 and T3 decreased by 1.6% and 1.7% compared with that of T1, respectively.A large number of studies have shown that under certain conditions crop yield is positively correlated with transpiration [34,35].Figure 7 shows that the winter wheat yields of T2 and T3 increased by 10% and 33% compared with that of T1, respectively.He et al. [33] studied the effect of under subsurface drip irrigation on winter wheat yield under a buried depths of 20 and 40 cm and found that the yield for the buried depth of 40 cm was higher than that for the buried depth of 20 cm.This result is consistent with the results of the current study.Thus, deep irrigation can improve soil water use efficiency and increase winter wheat yield.
Water 2018, 10, x FOR PEER REVIEW 14 of 17 respectively, and the evapotranspiration (ET) of T2 and T3 decreased by 1.6% and 1.7% compared with that of T1, respectively.A large number of studies have shown that under certain conditions crop yield is positively correlated with transpiration [34,35].Figure 7 shows that the winter wheat yields of T2 and T3 increased by 10% and 33% compared with that of T1, respectively.He et al. studied the effect of under subsurface drip irrigation on winter wheat yield under a buried depths of 20 and 40 cm and found that the yield for the buried depth of 40 cm was higher than that for the buried depth of 20 cm.This result is consistent with the results of the current study.Thus, deep irrigation can improve soil water use efficiency and increase winter wheat yield.It can be concluded that T3 treatment is an appropriate irrigation treatment for winter wheat in the study area.It can reduce the evapotranspiration and increase the yield of winter wheat.Therefore, for field winter wheat irrigation, seepage pipes can be buried at soil depths from 0 m to 1.8 m at intervals of 0.3 m.Then, when the winter wheat needs irrigation at different growth stages, the seepage pipe of corresponding depth will be opened for irrigation based on the depth distribution of the winter wheat root system.However, the specific techniques of field irrigation need to be further studied in detail.

Conclusions
The soil water movement model of winter wheat was established by considering the RWU and different depths of irrigation and using the source term of the soil water movement equation to respectively, and the evapotranspiration (ET) of T2 and T3 decreased by 1.6% and 1.7% compared with that of T1, respectively.A large number of studies have shown that under certain conditions crop yield is positively correlated with transpiration [34,35].Figure 7 shows that the winter wheat yields of T2 and T3 increased by 10% and 33% compared with that of T1, respectively.He et al. [33] studied the effect of under subsurface drip irrigation on winter wheat yield under a buried depths of 20 and 40 cm and found that the yield for the buried depth of 40 cm was higher than that for the buried depth of 20 cm.This result is consistent with the results of the current study.Thus, deep irrigation can improve soil water use efficiency and increase winter wheat yield.It can be concluded that T3 treatment is an appropriate irrigation treatment for winter wheat in the study area.It can reduce the evapotranspiration and increase the yield of winter wheat.Therefore, for field winter wheat irrigation, seepage pipes can be buried at soil depths from 0 m to 1.8 m at intervals of 0.3 m.Then, when the winter wheat needs irrigation at different growth stages, the seepage pipe of corresponding depth will be opened for irrigation based on the depth distribution of the winter wheat root system.However, the specific techniques of field irrigation need to be further studied in detail.

Conclusions
The soil water movement model of winter wheat was established by considering the RWU and different depths of irrigation and using the source term of the soil water movement equation to It can be concluded that T3 treatment is an appropriate irrigation treatment for winter wheat in the study area.It can reduce the evapotranspiration and increase the yield of winter wheat.Therefore, for field winter wheat irrigation, seepage pipes can be buried at soil depths from 0 m to 1.8 m at intervals of 0.3 m.Then, when the winter wheat needs irrigation at different growth stages, the seepage pipe of corresponding depth will be opened for irrigation based on the depth distribution of the winter wheat root system.However, the specific techniques of field irrigation need to be further studied in detail.

Conclusions
The soil water movement model of winter wheat was established by considering the RWU and different depths of irrigation and using the source term of the soil water movement equation to simulate irrigation of different depths.The model was calibrated and validated by comparing the simulated and experimental SWCs.The results showed that the soil water movement model of winter wheat at different irrigation depths has high simulation accuracy and can be used to simulate soil water movement.The experimental results of winter wheat growth at different irrigation depths showed that deep irrigation has a significant influence on root distribution and can reduce the RLD in the upper soil layers and increase the RLD in the deeper soil layers.The dynamic simulation of RWU and SWC showed that deep irrigation can increase the SWC of deep soil and promote the RWU of deep soil.As a result, deep irrigation can increase the transpiration of winter wheat and increase the yield of winter wheat.Under the experimental conditions, T3 treatment (namely, an irrigation soil depth that is 70% that of the root distribution depth, irrigation 5 times in the whole growth period, and each irrigation quota of 67.5 mm) is a suitable irrigation treatment for winter wheat.The results of this study can provide a reference for the reasonable irrigation of winter wheat in northern China.

Figure 1 .
Figure 1.The schematic diagram of the test layout.(a) Field layout; (b) Test equipment.TRIME refers to the time domain reflectometry with intelligent microElements.

Figure 1 .
Figure 1.The schematic diagram of the test layout.(a) Field layout; (b) Test equipment.TRIME refers to the time domain reflectometry with intelligent microElements.

Figure 4 .
Figure 4.The comparison of the calculated and measured soil water dynamics of winter wheat during irrigation at different soil depths.(a) T1 treatment; (b) T2 treatment; (c) T3 treatment.

Figure 5 .
Figure 5.The dynamics of the RWU rate of winter wheat at different irrigation depths.(a) T1 treatment; (b) T2 treatment; (c) T3 treatment.

Figure 5 .
Figure 5.The dynamics of the RWU rate of winter wheat at different irrigation depths.(a) T1 treatment; (b) T2 treatment; (c) T3 treatment.

Figure 6 .Figure 7 .
Figure 6.The comparison of the total RWU of winter wheat in different soil layers (a) and the water consumption composition of winter wheat (b) at different soil depths irrigation.Tr is the transpiration of winter wheat, E is the evaporation of winter wheat, and ET is the evapotranspiration of winter wheat.

Figure 6 .
Figure 6.The comparison of the total RWU of winter wheat in different soil layers (a) and the water consumption composition of winter wheat (b) at different soil depths irrigation.Tr is the transpiration of winter wheat, E is the evaporation of winter wheat, and ET is the evapotranspiration of winter wheat.

Figure 6 .Figure 7 .
Figure 6.The comparison of the total RWU of winter wheat in different soil layers (a) and the water consumption composition of winter wheat (b) at different soil depths irrigation.Tr is the transpiration of winter wheat, E is the evaporation of winter wheat, and ET is the evapotranspiration of winter wheat.

Figure 7 .
Figure 7.The comparison of winter wheat yield under different irrigation depth conditions.

Table 1 .
The soil mechanical composition and physical parameters.

Table 1 .
The soil in the test area is silty clay loam.The mechanical composition and physical properties of the soil are shown in Table1.The soil mechanical composition and physical parameters.

Table 2 .
The irrigation time, irrigation hole setting, irrigation scheme depth, and irrigation quota of winter wheat.

Table 3 .
The results of the model parameter solution.

Table 4 .
The fitting results of the RLD distribution function.

Table 5 .
The model accuracy evaluation indicators.

Table 5 .
The model accuracy evaluation indicators.